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.