Fix sign of Newton iterate - consistent with builtin solver. (#13)

a-latyshev · Oct 24, 2024 · 89dc592 · 89dc592
1 parent 41a9baa
commit 89dc592
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/doc/demo/demo_plasticity_von_mises.py b/doc/demo/demo_plasticity_von_mises.py
@@ -478,7 +478,7 @@ def sigma_external(derivatives):
         external_operator_problem.solve(du)
         du.x.scatter_forward()
 
-        Du.x.petsc_vec.axpy(-1.0, du.x.petsc_vec)
+        Du.x.petsc_vec.axpy(1.0, du.x.petsc_vec)
         Du.x.scatter_forward()
 
         evaluated_operands = evaluate_operands(F_external_operators)