fixing petsc output vector error in the MPhys wrapper (#300)

* fixing petsc output vector issue * removing the comments * minor fix for LBGS * minor fix for LBGS * minor fix for LBGS * removing unnecessary copy --------- Co-authored-by: Alasdair Gray <[email protected]>
mdolab · Aug 2, 2023 · d161960 · d161960
1 parent 4dda7b3
commit d161960
Showing 1 changed file with 3 additions and 2 deletions.
diff --git a/adflow/mphys/mphys_adflow.py b/adflow/mphys/mphys_adflow.py
@@ -608,8 +608,9 @@ def solve_linear(self, d_outputs, d_residuals, mode):
         if mode == "fwd":
             d_outputs["adflow_states"] = solver.solveDirectForRHS(d_residuals["adflow_states"])
         elif mode == "rev":
-            # d_residuals['adflow_states'] = solver.solveAdjointForRHS(d_outputs['adflow_states'])
-            solver.adflow.adjointapi.solveadjoint(d_outputs["adflow_states"], d_residuals["adflow_states"], True)
+            phi = d_residuals["adflow_states"].copy()
+            solver.adflow.adjointapi.solveadjoint(d_outputs["adflow_states"], phi, True)
+            d_residuals["adflow_states"] = phi
 
         return True, 0, 0