Hi all,
I need to create a restart/resume checkpoint utility in my deal.ii based
code.
Looking at some previous posts, this possible by using:
parallel::distributed::SolutionTransfer
prepare_serialization/deserialize
triangulation.save/.load
This is working for me, except that the solution on the boundary nodes of
each partition is always zero while
for the interior nodes is OK.
I am using a parallel distributed triangulation and PETSc vector and here
is the extracted piece of the code:
/*
Save the solution
where sol_u is a PETScWrappers::MPI::Vector
*/
parallel::distributed::SolutionTransfer<dim,
PETScWrappers::MPI::Vector>system_trans (dof_handler);
system_trans.prepare_serialization (sol_u);
triangulation.save ( filename.c_str() );
//Reload the solution
parallel::distributed::SolutionTransfer<dim,
PETScWrappers::MPI::Vector>system_trans (dof_handler);
system_trans.deserialize(sol_u);
I know that there must be some issues with distributed vector/dofs but I
cannot figure it out.
--
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.
