Dear Daniel, Removing the assertion works, but it will cause a problem. Since the ConstraintMatrix internally uses local_lines to determine the position of a ConstraintLine, we can not insert new line or query an existing line anymore.
If we shift the local_lines as well, another problem arises. We intend to merge the ConstraintMatrix after shifting, but the merge function requires the two ConstraintMatrices have same local_lines. And ConstraintMatrix does not provide a way to modify the local_lines without affecting the internal data structures. I haven't found a proper way to do it. Currently, I rebuild the whole ConstraintMatrix to implement the shift-and-merge procedure. It is not very efficient, but it just works. Best, Ce Qin 2016-12-31 18:53 GMT+08:00 Daniel Arndt <[email protected]>: > 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 [email protected]. > For more options, visit https://groups.google.com/d/optout. > -- 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 [email protected]. For more options, visit https://groups.google.com/d/optout.
