]> https://gitweb.dealii.org/ - dealii-svn.git/commitdiff
Fix one border case in initialization on non-constant Jacobians.
authorkronbichler <kronbichler@0785d39b-7218-0410-832d-ea1e28bc413d>
Mon, 10 Dec 2012 10:33:22 +0000 (10:33 +0000)
committerkronbichler <kronbichler@0785d39b-7218-0410-832d-ea1e28bc413d>
Mon, 10 Dec 2012 10:33:22 +0000 (10:33 +0000)
git-svn-id: https://svn.dealii.org/trunk@27789 0785d39b-7218-0410-832d-ea1e28bc413d

deal.II/include/deal.II/matrix_free/mapping_info.templates.h

index 7d315d7d7646cbb9fc76cb1a53c2759fcb1f93e3..2d9ab249413bdb3210ad6a88cb2e559450e6662e 100644 (file)
@@ -734,21 +734,18 @@ end_set:
                             break;
                           }
 
-                  // in case we have only one quadrature point,
-                  // we can have non-constant Jacobians, but we
-                  // cannot detect it by comparison from one
-                  // quadrature point to the next: in that case,
-                  // need to look at second derivatives and see
-                  // whether there are some non-zero entries
-                  // (this is necessary since we determine the
-                  // constness of the Jacobian for the first
-                  // quadrature formula and might not look at
-                  // them any more for the second, third
-                  // quadrature formula
-                  if (cell_cartesian == false && n_q_points == 1)
+                  // in case we have only one quadrature point, we can have
+                  // non-constant Jacobians, but we cannot detect it by
+                  // comparison from one quadrature point to the next: in that
+                  // case, need to look at second derivatives and see whether
+                  // there are some non-zero entries (this is necessary since
+                  // we determine the constness of the Jacobian for the first
+                  // quadrature formula and might not look at them any more
+                  // for the second, third quadrature formula). in any case,
+                  // the flag update_jacobian_grads will be set in that case
+                  if (cell_cartesian == false && n_q_points == 1 &&
+                      update_flags & update_jacobian_grads)
                     {
-                      Assert (update_flags & update_jacobian_grads,
-                              ExcInternalError());
                       const DerivativeForm<1,dim,dim> &jac = fe_val.jacobian(0);
                       const DerivativeForm<2,dim,dim> &jacobian_grad =
                         fe_val.jacobian_grad(0);

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.