Hi Prof. Bangerth, Sorry that I forget to mention *it only breaks when runnning with MPI, run in serialization is good.* PersistentTriangulation Class may not address this error. The intention of write & read a mesh here is to create a not-to-coarse coarse grid (so it's refine several times before output). I read the output grid and expect all cells are on level = 0. Preserving the multigrid structure is not preferred in this case.
Do using a straight GridOut destroys the mesh structure even I expect to eliminate multigrid information? Best, Wayne 在2022年11月28日星期一 UTC+8 11:48:55<yy.wayne> 写道： > Thanks for your quick reply. > > I write mesh this way in order to create a refined coarse mesh. I tried > generate a distributed mesh in serial but it gives other bugs when further > refine it, therefore decide to write a mesh generated serialized and read > it back. > I'll learn the PersistentTriangulation class then, thank you. > > Best, > Wayne > > 在2022年11月28日星期一 UTC+8 11:26:04<Wolfgang Bangerth> 写道： > >> 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: bang...@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/9bb446b8-f36c-41d5-8d4d-ac1c313b2daan%40googlegroups.com.
