template<int dim>
void test (const Triangulation<dim> &tr,
- const FiniteElement<dim> &fe)
+ const FiniteElement<dim> &fe,
+ const unsigned int degree)
{
deallog << "FE=" << fe.get_name()
<< std::endl;
// equals the vector function (0,x^2)
ConstraintMatrix cm;
cm.close ();
- VectorTools::project (dof, cm, QGauss<2>(4), F(), fe_function);
+ VectorTools::project (dof, cm, QGauss<2>(2+degree), F(), fe_function);
const QGauss<dim> quadrature(2);
FEValues<dim> fe_values (fe, quadrature,
Triangulation<dim> tr;
GridGenerator::hyper_cube(tr);
- FE_Nedelec<dim> fe(2);
- test(tr, fe);
+ for (unsigned int degree=2; degree<5; ++degree)
+ {
+ FE_Nedelec<dim> fe(degree);
+ test(tr, fe, degree);
+ }
}
DEAL:: grads[q]= 0.00 0.00 0.423 -2.00e-16
DEAL:: curls[q]= 1.58
DEAL:: grads[q]= 0.00 0.00 1.58 -3.38e-16
+DEAL::FE=FE_Nedelec<2>(3)
+DEAL:: curls[q]= 0.423
+DEAL:: grads[q]= 0.00 0.00 0.423 5.69e-17
+DEAL:: curls[q]= 1.58
+DEAL:: grads[q]= 0.00 0.00 1.58 1.29e-15
+DEAL:: curls[q]= 0.423
+DEAL:: grads[q]= 0.00 0.00 0.423 2.40e-15
+DEAL:: curls[q]= 1.58
+DEAL:: grads[q]= 0.00 0.00 1.58 6.52e-16
+DEAL::FE=FE_Nedelec<2>(4)
+DEAL:: curls[q]= 0.423
+DEAL:: grads[q]= 0.00 0.00 0.423 -1.06e-15
+DEAL:: curls[q]= 1.58
+DEAL:: grads[q]= 0.00 0.00 1.58 2.32e-15
+DEAL:: curls[q]= 0.423
+DEAL:: grads[q]= 0.00 0.00 0.423 1.64e-15
+DEAL:: curls[q]= 1.58
+DEAL:: grads[q]= 0.00 0.00 1.58 -2.01e-15