Richard, [...] > The very last part (reinit) is working only iff every processor owns parts > of fluid AND solid domain (which is so far not ensured). > This behaviour can be seen running exactly the same problem with mpirun -n > 2 or -n 3, which is attached. > > THE QUESTION IS NOW: > How may one reinit the blocks in a BlockSparseMatrix, if not every > subdomain has locally_owned_dofs in all blocks? > Or what is the point I am missing here? > > I added a chunk of code at the end, that shows, that the > block(block_row,block_col).reinit(...) fails, if the locally_owned indexset > is empty. >
This assertion was removed in https://github.com/dealii/dealii/pull/6106. Either use a recent release or simply remove it in your copy of the deal.II source files. 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. To view this discussion on the web visit https://groups.google.com/d/msgid/dealii/CAOYDWbJhT0hB%3DT0r8qBXZ5QAp4%3DqpYugwRZzwaAhTZF7F3S1pg%40mail.gmail.com.
