]> https://gitweb.dealii.org/ - dealii.git/commitdiff
put nu and sigma at the front
authorJiaqi Zhang <zjiaqi@vt.edu>
Sun, 20 Sep 2020 02:47:19 +0000 (22:47 -0400)
committerJiaqi Zhang <zjiaqi@vt.edu>
Mon, 30 Nov 2020 17:06:40 +0000 (12:06 -0500)
examples/step-74/step-74.cc

index 92b06eb813b14b81a0a6a20d9fa73129c2672b86..04080a4bc95c91900ff9808ae6e043665a82e0e0 100644 (file)
@@ -351,8 +351,8 @@ namespace Step74
           {
             for (unsigned int j = 0; j < fe_v.dofs_per_cell; ++j)
               copy_data.cell_matrix(i, j) +=
-                fe_v.shape_grad(i, point) * // grad v_h
                 diffusion_coefficient *     // nu
+                fe_v.shape_grad(i, point) * // grad v_h
                 fe_v.shape_grad(j, point) * // grad u_h
                 JxW[point];                 // dx
 
@@ -389,33 +389,33 @@ namespace Step74
           for (unsigned int i = 0; i < dofs_per_cell; ++i)
             for (unsigned int j = 0; j < dofs_per_cell; ++j)
               copy_data.cell_matrix(i, j) +=
-                (-fe_fv.shape_value(i, point) *  // - v_h
-                   diffusion_coefficient *       // nu
+                (-diffusion_coefficient *        // - nu
+                   fe_fv.shape_value(i, point) * // v_h
                    (fe_fv.shape_grad(j, point) * // (grad u_h .
                     normals[point])              //  n)
 
-                 - (fe_fv.shape_grad(i, point) * // - (grad v_h .
-                    normals[point]) *            //    n)
-                     diffusion_coefficient *     // nu
-                     fe_fv.shape_value(j, point) // u_h
+                 - diffusion_coefficient *         // - nu
+                     (fe_fv.shape_grad(i, point) * // (grad v_h .
+                      normals[point]) *            //  n)
+                     fe_fv.shape_value(j, point)   // u_h
 
-                 +                                   // +
-                 fe_fv.shape_value(i, point) *       // v_h
-                   diffusion_coefficient * penalty * // nu sigma
-                   fe_fv.shape_value(j, point)       // u_h
+                 + diffusion_coefficient * penalty * // + nu sigma
+                     fe_fv.shape_value(i, point) *   // v_h
+                     fe_fv.shape_value(j, point)     // u_h
 
                  ) *
                 JxW[point]; // dx
 
           for (unsigned int i = 0; i < dofs_per_cell; ++i)
             copy_data.cell_rhs(i) +=
-              (-(fe_fv.shape_grad(i, point) * // - (grad v_h .
-                 normals[point]) *            //    n)
-                 diffusion_coefficient *      // nu
-                 g[point]                     // g
+              (-diffusion_coefficient *        // - nu
+                 (fe_fv.shape_grad(i, point) * // (grad v_h .
+                  normals[point]) *            //  n)
+                 g[point]                      // g
+
 
-               + fe_fv.shape_value(i, point) *                // + v_h
-                   diffusion_coefficient * penalty * g[point] // nu sigma g
+               + diffusion_coefficient * penalty *        // + nu sigma
+                   fe_fv.shape_value(i, point) * g[point] // v_h g
 
                ) *
               JxW[point]; // dx
@@ -460,22 +460,19 @@ namespace Step74
           for (unsigned int i = 0; i < n_dofs_face; ++i)
             for (unsigned int j = 0; j < n_dofs_face; ++j)
               copy_data_face.cell_matrix(i, j) +=
-                (-fe_iv.jump(i, point) *               // - [v_h]
-                   diffusion_coefficient *             // nu
+                (-diffusion_coefficient *              // - nu
+                   fe_iv.jump(i, point) *              // [v_h]
                    (fe_iv.average_gradient(j, point) * // ({grad u_h} .
                     normals[point])                    //  n)
 
+                 - diffusion_coefficient *               // - nu
+                     (fe_iv.average_gradient(i, point) * // (grad v_h .
+                      normals[point]) *                  //  n)
+                     fe_iv.jump(j, point)                // [u_h]
 
-                 - (fe_iv.average_gradient(i, point) * // -(grad v_h .
-                    normals[point]) *                  //   n)
-                     diffusion_coefficient *           // nu
-                     fe_iv.jump(j, point)              // [u_h]
-
-
-                 + fe_iv.jump(i, point) *              // + [v_h]
-                     diffusion_coefficient * penalty * // nu sigma
-                     fe_iv.jump(j, point)              // [u_h]
-
+                 + diffusion_coefficient * penalty * // + nu sigma
+                     fe_iv.jump(i, point) *          // [v_h]
+                     fe_iv.jump(j, point)            // [u_h]
 
                  ) *
                 JxW[point]; // dx

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.