// this gives a Cartesian cell but in non-standard orientation (x-coordinate
// is gone through backwards)
- MappingQEulerian<dim> mapping(1,displacements, dof_h);
+ MappingQEulerian<dim> mapping(1, dof_h, displacements);
QGauss<dim> quad(1);
FEValues<dim> fe_val (mapping, fe, quad, update_JxW_values);
double integral = 0.;
grid_out.write_ucd (tria, logfile);
QTrapez<dim> quad;
- MappingQEulerian<dim,Vector<double>,spacedim> mapping(degree, shift, shift_dh);
+ MappingQEulerian<dim,Vector<double>,spacedim> mapping(degree, shift_dh, shift);
typename Triangulation<dim,spacedim>::active_cell_iterator cell=tria.begin_active(),
endc=tria.end() ;
{
QGauss<dim> quadrature_formula(degree+1);
- MappingQEulerian<dim> mapping(degree,displacements,dof_handler);
+ MappingQEulerian<dim> mapping(degree, dof_handler, displacements);
FEValues<dim> fe_values (mapping, fe, quadrature_formula,
update_JxW_values);
Displacement<dim>(),
displacements);
- MappingQEulerian<dim> euler(2, displacements, dof_handler);
+ MappingQEulerian<dim> euler(2, dof_handler, displacements);
// now the actual test
DataOut<dim> data_out;
data_out.attach_dof_handler(dof_handler);
Displacement<dim>(),
displacements);
- MappingQEulerian<dim> euler(2, displacements, dof_handler);
+ MappingQEulerian<dim> euler(2, dof_handler, displacements);
// now the actual test
DataOut<dim> data_out;
data_out.attach_dof_handler(dof_handler);
Displacement<dim>(),
displacements);
- MappingQEulerian<dim> euler(2, displacements, dof_handler);
+ MappingQEulerian<dim> euler(2, dof_handler, displacements);
// now the actual test
DataOut<dim> data_out;
data_out.attach_dof_handler(dof_handler);
grid_out.write_ucd (tria, logfile);
QTrapez<dim> quad;
- MappingQEulerian<dim,Vector<double>,spacedim> mapping(degree,shift, shift_dh);
+ MappingQEulerian<dim,Vector<double>,spacedim> mapping(degree, shift_dh, shift);
typename Triangulation<dim,spacedim>::active_cell_iterator cell=tria.begin_active(),
endc=tria.end() ;
shift(i) = 0.1;
QGauss<dim> quad(degree+1);
- MappingQEulerian<dim,Vector<double>,spacedim> mapping(degree,shift, shift_dh);
+ MappingQEulerian<dim,Vector<double>,spacedim> mapping(degree, shift_dh, shift);
Triangulation<dim,spacedim>::active_cell_iterator cell=tria.begin_active(),
endc=tria.end() ;
shift(0) = 0.121;
shift(1) = -0.005;
}
- MappingQEulerian<dim> mapping(2,shift,dofh_eulerian);
+ MappingQEulerian<dim> mapping(2, dofh_eulerian, shift);
{
const QGauss<1> quad (fe_degree+1);
GridGenerator::hyper_cube (triangulation);
triangulation.refine_global (2);
- SparsityPattern cell_connectivity;
+ DynamicSparsityPattern cell_connectivity;
GridTools::get_face_connectivity_of_cells (triangulation, cell_connectivity);
- partition (cell_connectivity, 5);
+ SparsityPattern sp_cell_connectivity;
+ sp_cell_connectivity.copy_from (cell_connectivity);
+ partition (sp_cell_connectivity, 5);
}
else
break;
}
- catch (const typename Functions::FEFieldFunction<dim,DoFHandler<dim>,TrilinosWrappers::MPI::Vector>::ExcPointNotAvailableHere &)
+ catch (const VectorTools::ExcPointNotAvailableHere &)
{
deallog << " ExcPointNotAvailableHere" << std::endl;
}
combined_cell->set_dof_values (dxy, displacements);
}
// and create the MappingQEulerian
- MappingQEulerian<2> euler(2, displacements, cdh);
+ MappingQEulerian<2> euler(2, cdh, displacements);
// now the actual test
DataOut<2> data_out;
data_out.attach_dof_handler(cdh);