Dear All,
I have a fourth order parabolic equation, which is split into two second
order equations. Hence I solve using components by
declaring FESystem<dim> fe(FE_Q<dim>(degree,2). I have in total three such
sub-problems, which are coupled to each other in a semi-implicit manner.
Therefore I have a std::vector of solution for the entire system.
std::vector<Vector<double>> solution.
In addition, I am employing mesh adaptivity. After estimating the error
using Kelly, I would like to perform a solution transfer from old to new
mesh. Do I need to create a std::vector of SolutionTransfer objects; one
for each solution?
The below code copied from step-26 seems to work. Is this the correct
approach?
std::vector<SolutionTransfer<dim>> solution_trans(3,dof_handler);
std::vector<Vector<double>> previous_solution(num_index);
for(unsigned int i=0; i < num_index; i++){previous_solution[i] =
solution[i];}
triangulation.prepare_coarsening_and_refinement();
for(unsigned int i=0; i < num_index;
i++){solution_trans[i].prepare_for_coarsening_and_refinement(previous_solution[i]);}
triangulation.execute_coarsening_and_refinement();
setup_system();
for(unsigned int i=0; i < num_index; i++){
solution_trans[i].interpolate(previous_solution[i], solution[i]);
constraints.distribute(solution[i]);}
I look forward to your response.
Best regards,
Karthi.
--
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 dealii+unsubscr...@googlegroups.com.
To view this discussion on the web visit
https://groups.google.com/d/msgid/dealii/c10bed1e-d654-49ab-ba86-da3a11eaf5f0n%40googlegroups.com.
