]> https://gitweb.dealii.org/ - dealii.git/commitdiff
Use MappingQ in aniso test.
authorRalf Hartmann <Ralf.Hartmann@dlr.de>
Fri, 12 Jun 2009 05:25:21 +0000 (05:25 +0000)
committerRalf Hartmann <Ralf.Hartmann@dlr.de>
Fri, 12 Jun 2009 05:25:21 +0000 (05:25 +0000)
git-svn-id: https://svn.dealii.org/trunk@18926 0785d39b-7218-0410-832d-ea1e28bc413d

tests/aniso/mesh_3d_21.cc

index 1725dbea60b1cc3020c15c55a27af417d0ce06fe..601a257c17655277aadb7e8aeffcea13acdc313f 100644 (file)
@@ -30,6 +30,7 @@
 #include <dofs/dof_accessor.h>
 #include <fe/fe_dgq.h>
 #include <fe/fe_values.h>
+#include <fe/mapping_q.h>
 
 #include <fstream>
 
 // these objects which is considerable
 FE_DGQ<3> fe(1);
 QGauss<2> quadrature(3);
-FEFaceValues<3> fe_face_values1 (fe, quadrature,
+MappingQ<3> mapping(2);
+FEFaceValues<3> fe_face_values1 (mapping, fe, quadrature,
                                 update_q_points | update_JxW_values |
                                 update_normal_vectors);
-FESubfaceValues<3> fe_face_values2 (fe, quadrature,
+FESubfaceValues<3> fe_face_values2 (mapping, fe, quadrature,
                                    update_q_points | update_JxW_values |
                                    update_normal_vectors);
 
@@ -77,7 +79,7 @@ void check_this (Triangulation<3> &tria)
                         ExcInternalError());
 
                 Assert (std::fabs(fe_face_values1.JxW(q)-
-                                  fe_face_values2.JxW(q)) < 1e-15,
+                                  fe_face_values2.JxW(q)) < 1e-14,
                         ExcInternalError());
                 Assert ((fe_face_values1.normal_vector(q) +
                          fe_face_values2.normal_vector(q)).square()

In the beginning the Universe was created. This has made a lot of people very angry and has been widely regarded as a bad move.

Douglas Adams


Typeset in Trocchi and Trocchi Bold Sans Serif.