On 11/27/22 19:14, 'yy.wayne' via deal.II User Group wrote:

I get an error from compute_component_wise with level=0. I reduced the problem to a minimal working example, where I write a mesh and read to create a dense coarse level. The error shows even there's only one triangulation level. I assume the write & read mesh is causing this bug, since other part is similiar to dealii/tests/multigrid/renumbering_06.cc.

Writing the mesh just outputs the active cells, and loses the hierarchy that resulted from refinement. When you read it back in, you get a mesh with just one level.

You need a way to write and read the mesh that preserves the refinement levels. You can either to this via serialization, or using the PersistentTriangulation class.

Best
 W.
--
------------------------------------------------------------------------
Wolfgang Bangerth          email:                 bange...@colostate.edu
                           www: http://www.math.colostate.edu/~bangerth/


--
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/3d04f74a-0df7-3433-bd0d-db3d5ef7d00b%40colostate.edu.

Reply via email to