Hi
I am solving a problem using Trilinos and parallel::distributed <https://dealii.org/8.4.1/doxygen/deal.II/namespaceparallel_1_1distributed_1_1GridRefinement.html#a7f0ebf71078b6da3e0cbb895fcf15b12> triangulation. It works well. I transfer solution data and old solution data during mesh refinement as follow, std::vector<const LA::MPI::BlockVector *> x(2); x[0] = &relevant_solution; x[1] = &solution_old; parallel::distributed::SolutionTransfer<dim, LA::MPI::BlockVector> solution_transfer(dof_handler); solution_transfer.prepare_for_coarsening_and_refinement(x); triangulation.execute_coarsening_and_refinement(); setup_system(); LA::MPI::BlockVector tmp_v(partition); LA::MPI::BlockVector tmp_vv(partition); std::vector<LA::MPI::BlockVector *> tmp(2); tmp[0] = &solution; tmp[1] = &tmp_2; solution_transfer.interpolate(tmp); solution_old = tmp_2; Now, I also want to transfer quadrature point's data to the new mesh based on the method explained at the end of step-18 in which we can use a discontinuous field that matches the values in the quadrature points. Now using this method, I can prepare a global field (history_field) and I want to transfer the history_field vector as usual using the SolutionTransfer <https://dealii.org/8.4.1/doxygen/deal.II/classSolutionTransfer.html> class. Can I attach the history_field with relevant_solution and solution_old during solution transfer method? Or I must use two solution transfers, the first one for relevant_solution and solution using dof_handler and the second one for history_field using history_dof_handler? I can't understand this step and the relation between FE_DGQ for history field and FE_Q for the main problem during refinement. On Monday, June 2, 2014 at 3:57:20 AM UTC+4:30, Wolfgang Bangerth wrote: > > On 05/26/2014 03:43 AM, Ali Karrech wrote: > > Hi Timo, > > Thank you very much. Your comments helped me. I though it's a good idea > to > > share the solution that I adopted with other users. Apart from > > solution_transfer, I used FESystem and FE_DGQ to create a Finite element > field > > which has the same number of degrees as my integration points. I used > this new > > field (which I called history) to store my variable and recuperate them > when I > > need it. It seems to work properly so far. > > This is what I would have done as well. There is a function > FETools::compute_projection_from_quadrature_points_matrix that can help > you > with this, in case you haven't already found it. > > Best > W. > > -- > ------------------------------------------------------------------------ > Wolfgang Bangerth email: [email protected] > <javascript:> > www: http://www.math.tamu.edu/~bangerth/ > > -- 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.
