Hi Daniel, In the source of ConstraintMatrix::shift, there is an assertion that the local_lines must be empty. Please see ConstraintMatrix::shift <http://www.dealii.org/developer/doxygen/deal.II/constraint__matrix_8cc_source.html#l00592> .
Best, Ce Qin 2016-12-30 23:34 GMT+08:00 Daniel Arndt <d.ar...@math.uni-goettingen.de>: > Ce Qin, > > Am Freitag, 30. Dezember 2016 02:46:05 UTC+1 schrieb Ce Qin: >> >> Dear all, >> >> I want to build a 2x2 block matrix, while all blocks are built from a >> smaller DoFHandler. I need to generate 2 ConstraintMatrix and shift one >> of them, then merge them together. In the source code, Timo Heister said >> that the merge function doesn’t work with IndexSet yet. Can anyone >> provide some hints on how to implement it? >> > You can merge two ConstraintMatrix objects using ConstraintMatrix::shift > [1] and ConstraintMatrix::merge [2]. Also, you can merge two IndexSet > objects using > IndexSet::add_indices[3]. > What exactly is the problem you are experiencing? Which source code do you > refer to? > > Best, > Daniel > > [1] https://www.dealii.org/8.4.0/doxygen/deal.II/ > classConstraintMatrix.html#ad971fb29f0a15bf8d30704d2a5ab445e > [2] https://www.dealii.org/8.4.0/doxygen/deal.II/ > classConstraintMatrix.html#a08b5378ad8b23e5fb86c71a2c66bff13 > [3] https://www.dealii.org/8.4.1/doxygen/deal.II/classIndexSet.html# > a4d7d0507018483fb625ecf1358187e26 > >> -- > 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. > -- 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.
