Hi,
I like to link all the dofs of all the nodes along an edge of the body together. So I find the first dof (dof i) located on the edge and then I try to link all the dofs on that edge (except i) to dof i. It does not work if I have sufficiently large number of processes where that edge is chopped by different MPI processes. I know why it does not work, but I don't know any methods that can provide my desired dof constraints. Anyone has attempted to solve something similar? This is the error that I get with 8 MPI processes. -------------------------------------------------------- An error occurred in line <1607> of file </opt/local/dealii-9/deal.II-v9.0.1/include/deal.II/lac/constraint_matrix.h> in function void dealii::ConstraintMatrix::add_entry(dealii::ConstraintMatrix::size_type, dealii::ConstraintMatrix::size_type, double) The violated condition was: !local_lines.size() || local_lines.is_element(column) Additional information: The index set given to this constraint matrix indicates constraints using degree of freedom 580 should not be stored by this object, but a constraint for degree of freedom 736 uses it. -------------------------------------------------------- Regards, Reza -- The deal.II project is located at http://www.dealii.org/ For mailing list/forum options, see https://groups.google.com/d/forum/dealii?hl=en --- You received this message because you are subscribed to the Google Groups "deal.II User Group" group. To unsubscribe from this group and stop receiving emails from it, send an email to dealii+unsubscr...@googlegroups.com. To view this discussion on the web visit https://groups.google.com/d/msgid/dealii/31e7aedd-b087-41d4-a412-0e211bcb27bd%40googlegroups.com. For more options, visit https://groups.google.com/d/optout.