From: Jiaqi Zhang Date: Sun, 20 Sep 2020 02:47:19 +0000 (-0400) Subject: put nu and sigma at the front X-Git-Tag: v9.3.0-rc1~687^2~12 X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=89e97bf782a379cae986d4514b6114fe1f4b0bf0;p=dealii.git put nu and sigma at the front --- diff --git a/examples/step-74/step-74.cc b/examples/step-74/step-74.cc index 92b06eb813..04080a4bc9 100644 --- a/examples/step-74/step-74.cc +++ b/examples/step-74/step-74.cc @@ -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