GridGenerator::hyper_cube(tria0, -.4, .3);
GridGenerator::hyper_cube(tria1, -1, 1);
- tria0.refine_global(3);
+ tria0.refine_global((dim0 == 3) ? 2 : 3);
tria1.refine_global(2);
FE_Q<dim0, spacedim> fe0(1);
<< "Dofs 1 : " << dh1.n_dofs() << std::endl
<< "Constrained dofs: " << constraints0.n_constraints() << std::endl;
- QGauss<dim0> quad0(3); // Quadrature for coupling
- QGauss<dim1> quad1(3); // Quadrature for coupling
+ QGauss<dim0> quad0(2); // Quadrature for coupling
+ QGauss<dim1> quad1(2); // Quadrature for coupling
const double epsilon = 2 * std::max(GridTools::maximal_cell_diameter(tria0),
GridTools::maximal_cell_diameter(tria1));
DEAL::Dofs 1 : 5
DEAL::Constrained dofs: 0
DEAL::Epsilon: 1.00000
-DEAL::Error on constants: 1.72490
+DEAL::Error on constants: 1.72494
DEAL::dim0: 2, dim1: 1, spacedim: 2
DEAL::FE0 : FE_Q<2>(1)
DEAL::FE1 : FE_Q<1,2>(1)
DEAL::Dofs 1 : 5
DEAL::Constrained dofs: 0
DEAL::Epsilon: 1.00000
-DEAL::Error on constants: 1.87222
+DEAL::Error on constants: 1.87214
DEAL::dim0: 2, dim1: 2, spacedim: 2
DEAL::FE0 : FE_Q<2>(1)
DEAL::FE1 : FE_Q<2>(1)
DEAL::Dofs 1 : 25
DEAL::Constrained dofs: 0
DEAL::Epsilon: 1.41421
-DEAL::Error on constants: 4.61731
+DEAL::Error on constants: 4.61805
DEAL::dim0: 3, dim1: 2, spacedim: 3
DEAL::FE0 : FE_Q<3>(1)
DEAL::FE1 : FE_Q<2,3>(1)
-DEAL::Dofs 0 : 729
+DEAL::Dofs 0 : 125
DEAL::Dofs 1 : 25
DEAL::Constrained dofs: 0
DEAL::Epsilon: 1.41421
-DEAL::Error on constants: 4.79570
+DEAL::Error on constants: 4.79607
DEAL::dim0: 3, dim1: 3, spacedim: 3
DEAL::FE0 : FE_Q<3>(1)
DEAL::FE1 : FE_Q<3>(1)
-DEAL::Dofs 0 : 729
+DEAL::Dofs 0 : 125
DEAL::Dofs 1 : 125
DEAL::Constrained dofs: 0
DEAL::Epsilon: 1.73205
-DEAL::Error on constants: 10.9402
\ No newline at end of file
+DEAL::Error on constants: 10.9407