Ce Qin, I don't immediately see a reason why shifting local_lines in ConstraintMatrix::shift should cause problems as long as the IndexSet you initialize the ConstraintMatrix objects with accounts for all DoFs. What happens if you just remove that assertion?
Best, Daniel -- 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. For more options, visit https://groups.google.com/d/optout.
