From 552d137f9e09f2d86175972529eef4f94827b669 Mon Sep 17 00:00:00 2001 From: jvwilliams23 Date: Fri, 4 Oct 2024 11:53:42 +0100 Subject: [PATCH] fix case when order_h1 != order_hcurl --- src/sources/closed_coil.cpp | 3 ++- 1 file changed, 2 insertions(+), 1 deletion(-) diff --git a/src/sources/closed_coil.cpp b/src/sources/closed_coil.cpp index 0b17cd1e1..ce0380292 100644 --- a/src/sources/closed_coil.cpp +++ b/src/sources/closed_coil.cpp @@ -62,7 +62,6 @@ ClosedCoilSolver::Init(hephaestus::GridFunctions & gridfunctions, _mesh_parent = _h_curl_fe_space_parent->GetParMesh(); _order_hcurl = _h_curl_fe_space_parent->FEColl()->GetOrder(); - _order_h1 = _order_hcurl; // Optional FE Spaces and parameters if (!fespaces.Has(_h1_fespace_name)) @@ -70,6 +69,7 @@ ClosedCoilSolver::Init(hephaestus::GridFunctions & gridfunctions, logger.info("{} not found in fespaces when creating {}. Creating from mesh.", _h1_fespace_name, typeid(this).name()); + _order_h1 = _order_hcurl; // Need to free this memory after use. FEC not freed by ParFiniteElementSpace destructor! _h1_fe_space_parent_fec = @@ -81,6 +81,7 @@ ClosedCoilSolver::Init(hephaestus::GridFunctions & gridfunctions, else { _h1_fe_space_parent = fespaces.GetShared(_h1_fespace_name); + _order_h1 = _h1_fe_space_parent->FEColl()->GetOrder(); } if (!gridfunctions.Has(_source_electric_field_name))