> This results in failed assertion, since: rows[0].size() = 2145, > rows[1].size() = 2145, but bdsp.block(0,0).n_rows() = 1815, > bdsp.block(1,1).n_rows() = 330 etc...
The size() of each IndexSet needs to correspond to the size of each block of the sparsity pattern and this is what the assertion checks. I think that this check is correct. If you want your first block to be 1815 DoFs large, then rows[0].size() needs to be 1815. How are you creating these IndexSets? I find the function IndexSet::split_by_block() convenient to use. -- Timo Heister http://www.math.clemson.edu/~heister/ -- 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]. To view this discussion on the web visit https://groups.google.com/d/msgid/dealii/CAMRj59GiiUSyV2h7y1Y1nUTWm%3DWLvK9PTCQbXj5sAa57E%2BuiWg%40mail.gmail.com.
