]> https://gitweb.dealii.org/ - dealii-svn.git/commitdiff
Remove an assertion that I now works now, but that takes quite some
authorwolf <wolf@0785d39b-7218-0410-832d-ea1e28bc413d>
Fri, 17 Oct 2003 20:08:42 +0000 (20:08 +0000)
committerwolf <wolf@0785d39b-7218-0410-832d-ea1e28bc413d>
Fri, 17 Oct 2003 20:08:42 +0000 (20:08 +0000)
computing time to verify.

git-svn-id: https://svn.dealii.org/trunk@8104 0785d39b-7218-0410-832d-ea1e28bc413d

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

index ba804ac29d5d94c550d07abc344042b8cb1eefde..907dc6edb85afc0505a24b2357d9ee77ebf4c030 100644 (file)
@@ -399,9 +399,6 @@ estimate_some (const Mapping<dim>                  &mapping,
                                            ((!neumann_bc.empty() ||
                                              (coefficients != 0))  ?
                                             update_q_points : 0) |
-#ifdef DEBUG
-                                           update_q_points |
-#endif
                                            update_normal_vectors));
   FEFaceValues<dim> fe_face_values_neighbor (mapping,
                                             dof_handler.get_fe(),
@@ -410,12 +407,7 @@ estimate_some (const Mapping<dim>                  &mapping,
   FESubfaceValues<dim> fe_subface_values (mapping,
                                          dof_handler.get_fe(),
                                          quadrature,
-                                         UpdateFlags(
-                                            update_gradients
-#ifdef DEBUG
-                                            | update_q_points
-#endif
-                                            ));
+                                          update_gradients);
 
 
   active_cell_iterator cell = dof_handler.begin_active();
@@ -1036,25 +1028,6 @@ integrate_over_irregular_face (const DoFHandler<dim>               &dof_handler,
                                       // common @p{subface}.
       fe_face_values.reinit (neighbor_child, neighbor_neighbor);
 
-                                       // make sure that quadrature
-                                       // 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())
-                ||
-                cell->has_boundary_lines(),
-                ExcInternalError());
-
                                        // store the gradient of the
                                       // solution in psi
       for (unsigned int n=0; n<n_solution_vectors; ++n)

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.