Hello,
I changed backend of linear algebra from TrilinosWrappers
to PetscWrappers.
I need to interpolate(project) initial conditions on solution vector,
I have define InitialValuePhase<dim> () function, which
sets the condition.
This worked for trilinos
InitialValuePhase<dim> ()solution_phase_n
solution_phase_n.reinit (locally_owned_dofs_phase,
locally_relevant_dofs_phase,
MPI_COMM_WORLD);
TrilinosWrappers::Vector<double> vec_old_solution (dof_handler_phase.n_dofs
());
VectorTools::interpolate (dof_handler_phase, InitialValuePhase<dim> (),
vec_old_solution );
solution_phase_n =vec_old_solution;
However I needed to change the TrilinosWrappers for
PetscWrappers.
Now the interpolation does not work. I get errors like
"Ghosted elements present"
Could You give me an advice, how to interpolate initial
conditions on PetscWrappers::MPI::Vector?
Thank You
Marek
--
The deal.II project is located at http://www.dealii.org/
For mailing list/forum options, see
https://groups.google.com/d/forum/dealii?hl=en
---
You received this message because you are subscribed to the Google Groups
"deal.II User Group" group.
To unsubscribe from this group and stop receiving emails from it, send an email
to [email protected].
For more options, visit https://groups.google.com/d/optout.
