BTW, Here I used GmshIO to read the .msh file of version 2.2. By default gmsh count $PhysicalNames from 1. This is also inherited by GmshIO::read_mesh(..).
I guess inside ParsedFunction it is assumed the user given boundary_id_type should start from zero, right? If so, there may be incompatibility problem like what I showed in last email. I tried force the boundary_id 1 to zero, then the assetation failure never occur, meanwhile the code give the wrong final result since with improper boundary conditions. Zhenyu ------------------------------------------------------------------------------ Dive into the World of Parallel Programming The Go Parallel Website, sponsored by Intel and developed in partnership with Slashdot Media, is your hub for all things parallel software development, from weekly thought leadership blogs to news, videos, case studies, tutorials and more. Take a look and join the conversation now. http://goparallel.sourceforge.net/ _______________________________________________ Libmesh-users mailing list [email protected] https://lists.sourceforge.net/lists/listinfo/libmesh-users
