diff --git a/src/solvers.jl b/src/solvers.jl index d64bdff..ddfceb9 100644 --- a/src/solvers.jl +++ b/src/solvers.jl @@ -199,6 +199,9 @@ function get_boundary_assembly(solver::Solver) new_constraints = get_nonzero_rows(C2_) overconstrained_dofs = intersect(already_constrained, new_constraints) if length(overconstrained_dofs) != 0 + warn("overconstrained dofs $overconstrained_dofs") + warn("already constrained = $already_constrained") + warn("new constraints = $new_constraints") overconstrained_dofs = sort(overconstrained_dofs) overconstrained_nodes = find_nodes_by_dofs(problem, overconstrained_dofs) warn("in overconstrained nodes $overconstrained_nodes")