DoFHandler<dim> dof_handler (tria);
dof_handler.distribute_dofs (fe);
- QGauss3<dim-1> q_face;
+ QGauss<dim-1> q_face(3);
FEFaceValues<dim> fe_face_values (mapping, fe, q_face,
update_normal_vectors);
DoFHandler<dim> dof_handler (tria);
dof_handler.distribute_dofs (fe);
- QGauss3<dim-1> q_face;
+ QGauss<dim-1> q_face(3);
FEFaceValues<dim> fe_face_values (fe, q_face,
update_normal_vectors | update_JxW_values);
dof_handler.distribute_dofs (fe);
deallog.threshold_double(1.e-10);
- QGauss3<dim-1> q_face;
+ QGauss<dim-1> q_face(3);
FEFaceValues<dim> fe_face_values (mapping, fe, q_face,
update_normal_vectors | update_JxW_values);
DoFHandler<dim> dof_handler (tria);
dof_handler.distribute_dofs (fe);
- QGauss3<dim-1> q_face;
+ QGauss<dim-1> q_face(3);
FEFaceValues<dim> fe_face_values (mapping, fe, q_face,
update_normal_vectors | update_JxW_values);
DoFHandler<dim> dof_handler (tria);
dof_handler.distribute_dofs (fe);
- QGauss3<dim-1> q_face;
+ QGauss<dim-1> q_face(3);
FEFaceValues<dim> fe_face_values (mapping, fe, q_face,
update_normal_vectors | update_JxW_values);
DoFHandler<dim> dof_handler (tria);
dof_handler.distribute_dofs (fe);
- QGauss3<dim-1> q_face;
+ QGauss<dim-1> q_face(3);
FEFaceValues<dim> fe_face_values (fe, q_face,
update_q_points | update_JxW_values);
DoFHandler<dim> dof_handler (tria);
dof_handler.distribute_dofs (fe);
- QGauss3<dim-1> q_face;
+ QGauss<dim-1> q_face(3);
FEFaceValues<dim> fe_face_values (mapping, fe, q_face,
update_q_points | update_JxW_values);
DoFHandler<dim> dof_handler (tria);
dof_handler.distribute_dofs (fe);
- QGauss3<dim-1> q_face;
+ QGauss<dim-1> q_face(3);
FEFaceValues<dim> fe_face_values (mapping, fe, q_face,
update_q_points | update_JxW_values);
DoFHandler<dim> dof_handler (tria);
dof_handler.distribute_dofs (fe);
- QGauss3<dim-1> q_face;
+ QGauss<dim-1> q_face(3);
FEFaceValues<dim> fe_face_values (mapping, fe, q_face,
update_q_points | update_JxW_values);
DoFHandler<dim> dof_handler (tria);
dof_handler.distribute_dofs (fe);
- QGauss3<dim-1> q_face;
+ QGauss<dim-1> q_face(3);
FEFaceValues<dim> fe_face_values (fe, q_face,
update_normal_vectors |
DoFHandler<dim> dof_handler (tria);
dof_handler.distribute_dofs (fe);
- QGauss3<dim-1> q_face;
+ QGauss<dim-1> q_face(3);
FEFaceValues<dim> fe_face_values (mapping, fe, q_face,
update_normal_vectors |
DoFHandler<dim> dof_handler (tria);
dof_handler.distribute_dofs (fe);
- QGauss3<dim-1> q_face;
+ QGauss<dim-1> q_face(3);
FEFaceValues<dim> fe_face_values (mapping, fe, q_face,
update_normal_vectors |
DoFHandler<dim> dof_handler (tria);
dof_handler.distribute_dofs (fe);
- QGauss3<dim-1> q_face;
+ QGauss<dim-1> q_face(3);
FEFaceValues<dim> fe_face_values (mapping, fe, q_face,
update_normal_vectors |
// formula suited to the elements
// we have here
MappingQ<dim> mapping (3);
- QGauss6<dim> quadrature;
+ QGauss<dim> quadrature(6);
// create sparsity pattern. note
// that different components should
// formula suited to the elements
// we have here
MappingQ<dim> mapping (3);
- QGauss6<dim> quadrature;
+ QGauss<dim> quadrature(6);
// create sparsity pattern. note
// that different components should
// formula suited to the elements
// we have here
MappingQ<dim> mapping (3);
- QGauss6<dim> quadrature;
+ QGauss<dim> quadrature(6);
// create sparsity pattern. note
// that different components should
// formula suited to the elements
// we have here
MappingQ<dim> mapping (3);
- QGauss6<dim> quadrature;
+ QGauss<dim> quadrature(6);
// create sparsity pattern. note
// that different components should
// formula suited to the elements
// we have here
MappingQ<dim> mapping (3);
- QGauss6<dim> quadrature;
+ QGauss<dim> quadrature(6);
// create sparsity pattern. note
// that different components should
// formula suited to the elements
// we have here
MappingQ<dim> mapping (3);
- QGauss6<dim> quadrature;
+ QGauss<dim> quadrature(6);
// create sparsity pattern. note
// that different components should
// formula suited to the elements
// we have here
MappingQ<dim> mapping (3);
- QGauss6<dim> quadrature;
+ QGauss<dim> quadrature(6);
// create sparsity pattern. note
// that different components should
// formula suited to the elements
// we have here
MappingQ<dim> mapping (3);
- QGauss6<dim> quadrature;
+ QGauss<dim> quadrature(6);
// create sparsity pattern. note
// that different components should
// formula suited to the elements
// we have here
MappingQ<dim> mapping (3);
- QGauss6<dim> quadrature;
+ QGauss<dim> quadrature(6);
// create sparsity pattern. note
// that different components should
// formula suited to the elements
// we have here
MappingQ<dim> mapping (3);
- QGauss6<dim> quadrature;
+ QGauss<dim> quadrature(6);
// create sparsity pattern. note
// that different components should
// formula suited to the elements
// we have here
MappingQ<dim> mapping (3);
- QGauss6<dim> quadrature;
+ QGauss<dim> quadrature(6);
// create sparsity pattern. note
// that different components should
// formula suited to the elements
// we have here
MappingQ<dim> mapping (3);
- QGauss6<dim> quadrature;
+ QGauss<dim> quadrature(6);
// create sparsity pattern. note
// that different components should
// formula suited to the elements
// we have here
MappingQ<dim> mapping (3);
- QGauss6<dim> quadrature;
+ QGauss<dim> quadrature(6);
// create sparsity pattern. note
// that different components should
// formula suited to the elements
// we have here
MappingQ<dim> mapping (3);
- QGauss6<dim> quadrature;
+ QGauss<dim> quadrature(6);
// create sparsity pattern. note
// that different components should
// formula suited to the elements
// we have here
MappingQ<dim> mapping (3);
- QGauss6<dim> quadrature;
+ QGauss<dim> quadrature(6);
// create sparsity pattern. note
// that different components should
// formula suited to the elements
// we have here
MappingQ<dim> mapping (3);
- QGauss6<dim> quadrature;
+ QGauss<dim> quadrature(6);
// create sparsity pattern. note
// that different components should
// formula suited to the elements
// we have here
MappingQ<dim> mapping (3);
- QGauss6<dim> quadrature;
+ QGauss<dim> quadrature(6);
// create sparsity pattern. note
// that different blocks should
// formula suited to the elements
// we have here
MappingQ<dim> mapping (3);
- QGauss6<dim> quadrature;
+ QGauss<dim> quadrature(6);
// create sparsity pattern. note
// that different components should
hn_constraints);
hn_constraints.close ();
VectorTools::project (mapping, *dof_handler, hn_constraints,
- QGauss6<2> (), pol,
+ QGauss<2> (6), pol,
solution);
// Now evaluate error ...
{
// Rotate element
VectorTools::project (*dof_handler_def, hn_constraints_def,
- QGauss6<2> (), TestDef1<2>(2, rotat),
+ QGauss<2> (6), TestDef1<2>(2, rotat),
deformation);
sprintf (buf, "phi = %e\n", rotat);
deallog << buf;
for (double scale = -0.75; scale < 4.0; scale += 0.25)
{
VectorTools::project (*dof_handler_def, hn_constraints_def,
- QGauss6<2> (), TestDef2<2>(2, scale),
+ QGauss<2> (6), TestDef2<2>(2, scale),
deformation);
sprintf (buf, "Scale = %e\n", scale);
deallog << buf;
for (double scale = -1.0; scale < 1.0; scale += 0.25)
{
VectorTools::project (*dof_handler_def, hn_constraints_def,
- QGauss6<2> (), TestDef3<2>(2, scale),
+ QGauss<2> (6), TestDef3<2>(2, scale),
deformation);
sprintf (buf, "Scale = %e\n", scale);
deallog << buf;
hn_constraints);
hn_constraints.close ();
VectorTools::project (*dof_handler, hn_constraints,
- QGauss6<2> (), TestMap1<2>(2),
+ QGauss<2> (6), TestMap1<2>(2),
solution);
// Project reference solution onto RT FE field
hn_constraints_def.close ();
VectorTools::project (*dof_handler_def, hn_constraints_def,
- QGauss6<2> (), TestMap1<2>(2),
+ QGauss<2> (6), TestMap1<2>(2),
solution_q);
MappingQ1Eulerian<2> *mapping_euler
{
// Rotate element
VectorTools::project (*dof_handler_def, hn_constraints_def,
- QGauss6<2> (), TestDef1<2>(2, rotat),
+ QGauss<2> (6), TestDef1<2>(2, rotat),
deformation);
// Project 1 function to element
VectorTools::project (*mapping_euler, *dof_handler, hn_constraints,
- QGauss6<2> (), TestMap1<2>(2),
+ QGauss<2> (6), TestMap1<2>(2),
solution);
// Write output files
for (double scale = -0.75; scale < 4.0; scale += 0.25)
{
VectorTools::project (*dof_handler_def, hn_constraints_def,
- QGauss6<2> (), TestDef2<2>(2, scale),
+ QGauss<2> (6), TestDef2<2>(2, scale),
deformation);
// Project 1 function to element
VectorTools::project (*mapping_euler, *dof_handler, hn_constraints,
- QGauss6<2> (), TestMap1<2>(2),
+ QGauss<2> (6), TestMap1<2>(2),
solution);
char buf[1000];
for (double scale = -1.0; scale < 1.0; scale += 0.25)
{
VectorTools::project (*dof_handler_def, hn_constraints_def,
- QGauss6<2> (), TestDef3<2>(2, scale),
+ QGauss<2> (6), TestDef3<2>(2, scale),
deformation);
// Project 1 function to element
VectorTools::project (*mapping_euler, *dof_handler, hn_constraints,
- QGauss6<2> (), TestMap1<2>(2),
+ QGauss<2> (6), TestMap1<2>(2),
solution);
char buf[1000];
hn_constraints);
hn_constraints.close ();
VectorTools::project (dof_handler, hn_constraints,
- QGauss6<2> (), TestMap1<2>(2),
+ QGauss<2> (6), TestMap1<2>(2),
solution);
// Project reference solution onto RT FE field
hn_constraints_def.close ();
VectorTools::project (dof_handler_def, hn_constraints_def,
- QGauss6<2> (), TestMap1<2>(2),
+ QGauss<2> (6), TestMap1<2>(2),
solution_q);
MappingQ1Eulerian<2> mapping_euler (deformation, dof_handler_def);
{
// Rotate element
VectorTools::project (dof_handler_def, hn_constraints_def,
- QGauss6<2> (), TestDef1<2>(2, rotat),
+ QGauss<2> (6), TestDef1<2>(2, rotat),
deformation);
sprintf (buf,"phi = %e FE_RT Area %e FE_Q Area %e\n",
rotat,
for (double scale = -0.75; scale < 4.0; scale += 0.25)
{
VectorTools::project (dof_handler_def, hn_constraints_def,
- QGauss6<2> (), TestDef2<2>(2, scale),
+ QGauss<2> (6), TestDef2<2>(2, scale),
deformation);
sprintf (buf,"Scale = %e FE_RT Area %e FE_Q Area %e\n",
scale,
for (double scale = -1.0; scale < 1.0; scale += 0.25)
{
VectorTools::project (dof_handler_def, hn_constraints_def,
- QGauss6<2> (), TestDef3<2>(2, scale),
+ QGauss<2> (6), TestDef3<2>(2, scale),
deformation);
sprintf (buf,"Scale = %e FE_RT Area %e FE_Q Area %e\n",
scale,
hn_constraints.close ();
MappingQ1<2> map_default;
project (map_default, *dof_handler, hn_constraints,
- QGauss6<2> (), ConstantFunction<2>(1., 2),
+ QGauss<2> (6), ConstantFunction<2>(1., 2),
solution);
EvaluateDerivative (dof_handler, solution);
hn_constraints);
hn_constraints.close ();
VectorTools::project (mapping, *dof_handler, hn_constraints,
- QGauss6<2> (), pol,
+ QGauss<2> (6), pol,
solution);
// Now evaluate error ...
{
// Rotate element
VectorTools::project (*dof_handler_def, hn_constraints_def,
- QGauss6<2> (), TestDef1<2>(2, rotat),
+ QGauss<2> (6), TestDef1<2>(2, rotat),
deformation);
deallog << "phi = " << rotat << std::endl;
TestProjection (mapping_euler, dof_handler);
for (double scale = -0.75; scale < 4.0; scale += 0.25)
{
VectorTools::project (*dof_handler_def, hn_constraints_def,
- QGauss6<2> (), TestDef2<2>(2, scale),
+ QGauss<2> (6), TestDef2<2>(2, scale),
deformation);
deallog << "scale = " << scale << std::endl;
TestProjection (mapping_euler, dof_handler);
for (double scale = -1.0; scale < 1.0; scale += 0.25)
{
VectorTools::project (*dof_handler_def, hn_constraints_def,
- QGauss6<2> (), TestDef3<2>(2, scale),
+ QGauss<2> (6), TestDef3<2>(2, scale),
deformation);
deallog << "scale = " << scale << std::endl;
TestProjection (mapping_euler, dof_handler);
// formula suited to the elements
// we have here
MappingQ<dim> mapping (3);
- QGauss6<dim> quadrature;
+ QGauss<dim> quadrature(6);
// create sparsity pattern. note
// that different components should
// formula suited to the elements
// we have here
MappingQ<dim> mapping (3);
- QGauss6<dim> quadrature;
+ QGauss<dim> quadrature(6);
// create sparsity pattern. note
// that different components should
// formula suited to the elements
// we have here
MappingQ<dim> mapping (3);
- QGauss6<dim> quadrature;
+ QGauss<dim> quadrature(6);
// create sparsity pattern. note
// that different components should
// formula suited to the elements
// we have here
MappingQ<dim> mapping (3);
- QGauss6<dim> quadrature;
+ QGauss<dim> quadrature(6);
// create sparsity pattern. note
// that different components should
// formula suited to the elements
// we have here
MappingQ<dim> mapping (3);
- QGauss6<dim> quadrature;
+ QGauss<dim> quadrature(6);
// create sparsity pattern. note
// that different components should
// formula suited to the elements
// we have here
MappingQ<dim> mapping (3);
- QGauss6<dim> quadrature;
+ QGauss<dim> quadrature(6);
// create sparsity pattern. note
// that different components should
// formula suited to the elements
// we have here
MappingQ<dim> mapping (3);
- QGauss6<dim> quadrature;
+ QGauss<dim> quadrature(6);
// create sparsity pattern. note
// that different components should
// formula suited to the elements
// we have here
MappingQ<dim> mapping (3);
- QGauss6<dim> quadrature;
+ QGauss<dim> quadrature(6);
// create sparsity pattern. note
// that different components should
// formula suited to the elements
// we have here
MappingQ<dim> mapping (3);
- QGauss6<dim> quadrature;
+ QGauss<dim> quadrature(6);
// create sparsity pattern. note
// that different components should
// formula suited to the elements
// we have here
MappingQ<dim> mapping (3);
- QGauss6<dim> quadrature;
+ QGauss<dim> quadrature(6);
// create sparsity pattern. note
// that different components should
// formula suited to the elements
// we have here
MappingQ<dim> mapping (3);
- QGauss6<dim> quadrature;
+ QGauss<dim> quadrature(6);
// create sparsity pattern. note
// that different components should
// formula suited to the elements
// we have here
MappingQ<dim> mapping (3);
- QGauss6<dim> quadrature;
+ QGauss<dim> quadrature(6);
// create sparsity pattern. note
// that different components should
// formula suited to the elements
// we have here
MappingQ<dim> mapping (3);
- QGauss6<dim> quadrature;
+ QGauss<dim> quadrature(6);
// create sparsity pattern. note
// that different components should
// formula suited to the elements
// we have here
MappingQ<dim> mapping (3);
- QGauss6<dim> quadrature;
+ QGauss<dim> quadrature(6);
// create sparsity pattern. note
// that different components should
// formula suited to the elements
// we have here
MappingQ<dim> mapping (3);
- QGauss6<dim> quadrature;
+ QGauss<dim> quadrature(6);
// create sparsity pattern. note
// that different components should
// formula suited to the elements
// we have here
MappingQ<dim> mapping (3);
- QGauss6<dim> quadrature;
+ QGauss<dim> quadrature(6);
// create sparsity pattern. note
// that different components should
// formula suited to the elements
// we have here
MappingQ<dim> mapping (3);
- QGauss6<dim> quadrature;
+ QGauss<dim> quadrature(6);
// create sparsity pattern. note
// that different blocks should