Hello Daniel,
with your definition the programme uses the default instantiations for
VECTOR and DH, which are Vector<double> and DoFHandler<dim>. But, since
you want to use hp::DoFHandler, you have to define it as
SolutionTransfer<dim, Vector<double>, hp::DoFHandler<dim> >
solution_transfer (dof_handler).
Best Regards,
Markus
On 25.04.2012 14:44, Daniel Riedlbauer wrote:
Hallo Markus,
I defined it this way:
triangulation.prepare_coarsening_and_refinement();
//Transfer actual solution to the refined grid
SolutionTransfer<dim> solution_transfer(dof_handler);
solution_transfer.prepare_for_pure_refinement();
//Transfer old solution to the refined grid
SolutionTransfer<dim> old_solution_transfer(dof_handler);
old_solution_transfer.prepare_for_pure_refinement();
//Execute refinement and coarsening
triangulation.execute_coarsening_and_refinement ();
set_which_element_gets_which_fe_type ();
//Distribute dofs on the new mesh
dof_handler.distribute_dofs (fe_collection);
Thanks
Daniel
_______________________________________________
dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
