Skip to content

Commit d9c91f7

Browse files
committed
Error in EigenSystem::solve if we have constraints
At least with SLEPc's eigen solver and the assembly routines in MOOSE, the result of a solve does not honor constraints such as hanging nodes. We do have `CondensedEigenSystem` which is a class intended for use with constraints. Unrelated to this PR that class will need to be updated in order to support shell matrices if MOOSE is to be able to use it
1 parent 7aa5322 commit d9c91f7

1 file changed

Lines changed: 4 additions & 0 deletions

File tree

src/systems/eigen_system.C

Lines changed: 4 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -208,6 +208,10 @@ void EigenSystem::reinit ()
208208

209209
void EigenSystem::solve ()
210210
{
211+
if (get_dof_map().n_constrained_dofs())
212+
libmesh_error_msg("EigenSystem does not support constrained degrees of freedom. If you wish to "
213+
"perform a solve on a system with constraints, then please use the "
214+
"CondensedEigenSystem class instead");
211215

212216
// A reference to the EquationSystems
213217
EquationSystems & es = this->get_equation_systems();

0 commit comments

Comments
 (0)