add extra information when overconstrained situation

This commit is contained in:
Jukka Aho
2017-01-09 09:32:02 +02:00
parent 9a382d9f7d
commit dcac3a9b53
+3
View File
@@ -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")