]> https://gitweb.dealii.org/ - dealii-svn.git/commitdiff
Fix the assertion.
authorwolf <wolf@0785d39b-7218-0410-832d-ea1e28bc413d>
Fri, 17 Oct 2003 15:35:48 +0000 (15:35 +0000)
committerwolf <wolf@0785d39b-7218-0410-832d-ea1e28bc413d>
Fri, 17 Oct 2003 15:35:48 +0000 (15:35 +0000)
git-svn-id: https://svn.dealii.org/trunk@8099 0785d39b-7218-0410-832d-ea1e28bc413d

deal.II/deal.II/source/numerics/error_estimator.cc

index 0e0d61db9a4f7a3f9efd4a3f5a3a28c03b14a69f..ba804ac29d5d94c550d07abc344042b8cb1eefde 100644 (file)
@@ -1037,13 +1037,22 @@ integrate_over_irregular_face (const DoFHandler<dim>               &dof_handler,
       fe_face_values.reinit (neighbor_child, neighbor_neighbor);
 
                                        // make sure that quadrature
-                                       // points match
+                                       // points match. note that this
+                                       // won't hold in 3d if one of
+                                       // the lines of the cell is at
+                                       // the boundary and if this is
+                                       // a curved boundary, since
+                                       // then the four subfaces don't
+                                       // exactly make up the mother
+                                       // cell. so exclude this case
       for (unsigned int q=0; q<n_q_points; ++q)
         Assert ((fe_face_values.quadrature_point(q) -
                  fe_subface_values.quadrature_point(q)).square()
                 <
                 1.e-15 * (fe_face_values.quadrature_point(q).square() + 
-                          fe_subface_values.quadrature_point(q).square()),
+                          fe_subface_values.quadrature_point(q).square())
+                ||
+                cell->has_boundary_lines(),
                 ExcInternalError());
 
                                        // store the gradient of the

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.