From 3166048dc1dfe56e4fc2d8a3ea34293a5703c861 Mon Sep 17 00:00:00 2001 From: Edward Palmer Date: Fri, 24 May 2024 16:06:56 +0000 Subject: [PATCH] Rebuild a0 solver in ScalarPotentialSource when Apply called. Required for updates to work correctly. Refer to TestTEAM4HFormMeshUpdates test. --- src/sources/scalar_potential_source.cpp | 5 +---- 1 file changed, 1 insertion(+), 4 deletions(-) diff --git a/src/sources/scalar_potential_source.cpp b/src/sources/scalar_potential_source.cpp index c3aabaae0..7bb61eb9b 100644 --- a/src/sources/scalar_potential_source.cpp +++ b/src/sources/scalar_potential_source.cpp @@ -117,11 +117,8 @@ ScalarPotentialSource::Apply(mfem::ParLinearForm * lf) _a0->FormLinearSystem( poisson_ess_tdof_list, phi_gf, *_b0, *_diffusion_mat, *_p_tdofs, *_b0_tdofs); - if (_a0_solver == nullptr) - { - _a0_solver = std::make_unique(_solver_options, *_diffusion_mat); - } // Solve + _a0_solver = std::make_unique(_solver_options, *_diffusion_mat); _a0_solver->Mult(*_b0_tdofs, *_p_tdofs); // "undo" the static condensation saving result in grid function dP