Skip to content

Commit 7701dc0

Browse files
committed
Add missing entries to FEAbstract::build()
Even if we can't use these dimensions later, FEMContext still wants to build them if it sees a variable of those types, just in case it hits an Elem of those dimensions.
1 parent 5a12fc5 commit 7701dc0

1 file changed

Lines changed: 24 additions & 0 deletions

File tree

src/fe/fe_abstract.C

Lines changed: 24 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -140,6 +140,18 @@ std::unique_ptr<FEAbstract> FEAbstract::build(const unsigned int dim,
140140
case SCALAR:
141141
return std::make_unique<FEScalar<0>>(fet);
142142

143+
case NEDELEC_ONE:
144+
return std::make_unique<FENedelecOne<0>>(fet);
145+
146+
case RAVIART_THOMAS:
147+
return std::make_unique<FERaviartThomas<0>>(fet);
148+
149+
case L2_RAVIART_THOMAS:
150+
return std::make_unique<FEL2RaviartThomas<0>>(fet);
151+
152+
case SUBDIVISION:
153+
return std::make_unique<FESubdivision>(fet);
154+
143155
default:
144156
libmesh_error_msg("ERROR: Bad FEType.family= " << Utility::enum_to_string(fet.family));
145157
}
@@ -205,6 +217,18 @@ std::unique_ptr<FEAbstract> FEAbstract::build(const unsigned int dim,
205217
case SCALAR:
206218
return std::make_unique<FEScalar<1>>(fet);
207219

220+
case NEDELEC_ONE:
221+
return std::make_unique<FENedelecOne<1>>(fet);
222+
223+
case RAVIART_THOMAS:
224+
return std::make_unique<FERaviartThomas<1>>(fet);
225+
226+
case L2_RAVIART_THOMAS:
227+
return std::make_unique<FEL2RaviartThomas<1>>(fet);
228+
229+
case SUBDIVISION:
230+
return std::make_unique<FESubdivision>(fet);
231+
208232
default:
209233
libmesh_error_msg("ERROR: Bad FEType.family= " << Utility::enum_to_string(fet.family));
210234
}

0 commit comments

Comments
 (0)