> > > Thank you. Is this class available in the release version 8.3.0? > I set separate dofhandler (history_dof_handler) for FE_DGQ. I just > confusing because the history_field vector has different size respect to > the main solution vector and it also has different DoFHandler. I attached > dof_handler defined for the main problem for FE_Q, to the > solution_transfer, so I am confusing how we can use this solution_transfer > which defined based on the dof_handler to transfer history_field?! >
On Tuesday, June 14, 2016 at 7:26:06 PM UTC+4:30, Denis Davydov wrote: > > Hi > > On 14 Jun 2016, at 16:51, dealii...@gmail.com <javascript:> wrote: > > 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? > > yes, you can. You would need to set up a separate DoFHandler and FE_DGQ > etc for this. > Alternatively, you can use this class to achieve the same: > https://www.dealii.org/developer/doxygen/deal.II/classparallel_1_1distributed_1_1ContinuousQuadratureDataTransfer.html > > > 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. > > > The filed you solve for in FEM and the auxiliary field for transferring > are not related. > > Regards, > Denis. > > > > > 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: bang...@math.tamu.edu >> 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 a topic in the > Google Groups "deal.II User Group" group. > To unsubscribe from this topic, visit > https://groups.google.com/d/topic/dealii/OSaxZJR19W8/unsubscribe. > To unsubscribe from this group and all its topics, send an email to > dealii+un...@googlegroups.com <javascript:>. > For more options, visit https://groups.google.com/d/optout. > > > -- 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. For more options, visit https://groups.google.com/d/optout.
