]> https://gitweb.dealii.org/ - dealii.git/commitdiff
Move the multiplication with the coefficient to the right place.
authorWolfgang Bangerth <bangerth@math.tamu.edu>
Wed, 5 Aug 1998 15:24:54 +0000 (15:24 +0000)
committerWolfgang Bangerth <bangerth@math.tamu.edu>
Wed, 5 Aug 1998 15:24:54 +0000 (15:24 +0000)
git-svn-id: https://svn.dealii.org/trunk@477 0785d39b-7218-0410-832d-ea1e28bc413d

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

index a03e61414747bd1e391128084b0c628069c932f0..647482a5017f00418e1f37a08b70753cb3ef425e 100644 (file)
@@ -319,6 +319,17 @@ integrate_over_regular_face (const active_cell_iterator &cell,
   for (unsigned int point=0; point<n_q_points; ++point)
     phi[point] = psi[point]*normal_vectors[point];
 
+                                  // if a coefficient was given: use that
+                                  // to scale the jump in the gradient
+  if (coefficient != 0)
+    {
+      vector<double>      coefficient_values (n_q_points);
+      coefficient->value_list (fe_face_values_cell.get_quadrature_points(),
+                              coefficient_values);
+      for (unsigned int point=0; point<n_q_points; ++point)
+       phi[point] *= coefficient_values[point];
+    };
+      
   
   if (face->at_boundary() == true)
                                     // neumann boundary face. compute
@@ -341,17 +352,6 @@ integrate_over_regular_face (const active_cell_iterator &cell,
        phi[point] -= g[point];
     };
   
-                                  // if a coefficient was given: use that
-                                  // to scale the jump in the gradient
-  if (coefficient != 0)
-    {
-      vector<double>      coefficient_values (n_q_points);
-      coefficient->value_list (fe_face_values_cell.get_quadrature_points(),
-                              coefficient_values);
-      for (unsigned int point=0; point<n_q_points; ++point)
-       phi[point] *= coefficient_values[point];
-    };
-      
   
                                   // now phi contains the following:
                                   // - for an internal face, phi=[a du/dn]

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.