Subtract rather than add elastic force term in RHS assembly

geodynamics · Nov 16, 2017 · 4068d6a · 4068d6a
1 parent c076a8f
commit 4068d6a
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/source/simulator/assemblers/stokes.cc b/source/simulator/assemblers/stokes.cc
@@ -258,7 +258,7 @@ namespace aspect
 
               if (force != NULL)
                 data.local_rhs(i) += (force->rhs_u[q] * scratch.phi_u[i]
-                                      + double_contract(force->rhs_e[q],scratch.grad_phi_u[i])
+                                      - double_contract(force->rhs_e[q],Tensor<2,dim>(scratch.grads_phi_u[i]))
                                       + pressure_scaling * force->rhs_p[q] * scratch.phi_p[i])
                                      * JxW;