Dear Shawn, The GridIn::read_abaqus() function is currently not equipped to deal with the codimension 1 case. I might be able to find some time to look at this in the next few days though. In addition to the mesh file that you’ve already given us, may you please construct and give us a fully structured grid in codimension 1 (so something like a rectangle warped in the third dimension) that also has some boundary colouring? That would help a great deal with implementing this feature, and would also ensure that the assignment of boundary ID’s is possible.
Best, Jean-Paul > On 19 Sep 2018, at 01:38, Yuxiang Wang <[email protected]> wrote: > > Prof. Bangerth, > > Thank you for the prompt response! Yes please find attached the minimal > program. The main body is pretty short (see below), modified from step-49. > Please let me know if you cannot reproduce it! > > ``` > int main () > { > Triangulation<2, 3> triangulation; > > GridIn<2, 3> gridin; > gridin.attach_triangulation(triangulation); > std::ifstream f("mesh.inp"); > gridin.read_abaqus(f); > > GridOut gridout; > std::ofstream fout("mesh.vtu"); > gridout.write_vtu(triangulation, fout); > } > > ``` > > > Shawn > > On Tue, Sep 18, 2018 at 4:28 PM Wolfgang Bangerth <[email protected] > <mailto:[email protected]>> wrote: > On 09/18/2018 03:18 PM, Yuxiang Wang wrote: > > > > Thank you for the pointer! I changed my code as you suggested and it > > finished > > running. > > Great! > > > > However, when I opened the output file (mesh.vtu), the 3D mesh was > > projected > > to a 2D plane of z=0, and lost its 3D features. Was that expected behavior? > > No. Are you saying that if you read a mesh, then immediately output it again, > that the z-coordinate is lost? If that's correct, then that's a bug and it > would be nice to have a small program and input file that demonstrates this. > > Best > W. > > -- > ------------------------------------------------------------------------ > Wolfgang Bangerth email: [email protected] > <mailto:[email protected]> > www: http://www.math.colostate.edu/~bangerth/ > <http://www.math.colostate.edu/~bangerth/> > > -- > The deal.II project is located at http://www.dealii.org/ > <http://www.dealii.org/> > For mailing list/forum options, see > https://groups.google.com/d/forum/dealii?hl=en > <https://groups.google.com/d/forum/dealii?hl=en> > --- > You received this message because you are subscribed to a topic in the Google > Groups "deal.II User Group" group. > To unsubscribe from this topic, visit > https://groups.google.com/d/topic/dealii/Erk7pxwTEj4/unsubscribe > <https://groups.google.com/d/topic/dealii/Erk7pxwTEj4/unsubscribe>. > To unsubscribe from this group and all its topics, send an email to > [email protected] > <mailto:dealii%[email protected]>. > For more options, visit https://groups.google.com/d/optout > <https://groups.google.com/d/optout>. > > > -- > Yuxiang "Shawn" Wang, PhD > [email protected] <mailto:[email protected]> > +1 (434) 284-0836 > > -- > The deal.II project is located at http://www.dealii.org/ > <http://www.dealii.org/> > For mailing list/forum options, see > https://groups.google.com/d/forum/dealii?hl=en > <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] > <mailto:[email protected]>. > For more options, visit https://groups.google.com/d/optout > <https://groups.google.com/d/optout>. > <step-49.cc><mesh.inp><CMakeLists.txt><mesh.vtu> -- 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]. For more options, visit https://groups.google.com/d/optout.
