From 210f2eaa60dd4070acea9b16630988ac3d58c199 Mon Sep 17 00:00:00 2001 From: ddundo Date: Thu, 14 Nov 2024 11:00:04 +0000 Subject: [PATCH] #223: Don't update forward_old soln if steady --- goalie/go_mesh_seq.py | 5 +++-- 1 file changed, 3 insertions(+), 2 deletions(-) diff --git a/goalie/go_mesh_seq.py b/goalie/go_mesh_seq.py index 47a1066..d79706c 100644 --- a/goalie/go_mesh_seq.py +++ b/goalie/go_mesh_seq.py @@ -212,7 +212,7 @@ def indicate_errors( } for f, fs_e in enriched_spaces.items(): if self.field_types[f] == "steady": - u[f] = u_[f] = enriched_mesh_seq.fields[f] + u[f] = enriched_mesh_seq.fields[f] else: u[f], u_[f] = enriched_mesh_seq.fields[f] u_star[f] = Function(fs_e) @@ -233,7 +233,8 @@ def indicate_errors( for f in self.fields: # Transfer solutions associated with the current field f transfer(self.solutions[f][FWD][i][j], u[f]) - transfer(self.solutions[f][FWD_OLD][i][j], u_[f]) + if self.field_types[f] == "unsteady": + transfer(self.solutions[f][FWD_OLD][i][j], u_[f]) transfer(self.solutions[f][ADJ][i][j], u_star[f]) transfer(self.solutions[f][ADJ_NEXT][i][j], u_star_next[f])