Merge pull request #6078 from quangx/bfbt_scaling

Change Weighted BFBT preconditioner pressure scaling factor
geodynamics · Oct 15, 2024 · 727d714 · 727d714
2 parents 03b739e + f4595aa
commit 727d714
Showing 1 changed file with 0 additions and 1 deletion.
diff --git a/source/simulator/assemblers/stokes.cc b/source/simulator/assemblers/stokes.cc
@@ -165,7 +165,6 @@ namespace aspect
                       if (scratch.dof_component_indices[i] == pressure_component_index && scratch.dof_component_indices[j] == pressure_component_index)
                         data.local_matrix(i, j) += (
                                                      1.0/sqrt_eta * pressure_scaling
-                                                     * pressure_scaling
                                                      * (scratch.grad_phi_p[i]
                                                         * scratch.grad_phi_p[j] + 1e-6*scratch.phi_p[i]*scratch.phi_p[j] ))
                                                    * JxW;