Hi all, I want to be able to start a calculation using a grid that was generated and used in a previous calculation. The obvious choice seems to be to use GridOut.write_msh()/GridIn.read_msh(). However, while write_msh() has an option to save surface ids (write_faces), for some strange, strange reason, write_msh does not seem able to write cell material ids at all, which means I can't reuse the written mesh.
According to the GridOut documentation, the only format for which material ids are written is DX; but then the GridIn class doesn't *read* that. I can't even seem to find a way to dump the internal format to disk and read it back raw. My attempt at a solution has been to write a separate ".materials" file with all the cell->material_id()s, but writing and reading a GMSH-file seems to reorder the cells, and I can't seem to find unique cell ids with which to associate them. I can't really believe that it actually is this difficult to write the full mesh to disk and read it back later - I feel it must be hiding under my nose somewhere - but I've been staring and staring at the documentation and can't seem to find anything. Any help will be greatly welcomed! Thanks a lot in advance, -- Med venlig hilsen, James Avery <[email protected]> _______________________________________________ dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
