Re: [Libmesh-users] Mesh with subdomains problem

2016-04-28 Thread Pepijn Kessels
Ah perfect, thought it had to do with libMesh. Thanks! On Thu, Apr 28, 2016 at 6:45 PM, Roy Stogner wrote: > > On Thu, 28 Apr 2016, Pepijn Kessels wrote: > > One more question. If it is possible, how can I print out the nodal >> positions in double precision? Right now I print these out to a fil

Re: [Libmesh-users] Mesh with subdomains problem

2016-04-28 Thread Roy Stogner
On Thu, 28 Apr 2016, Pepijn Kessels wrote: > One more question. If it is possible, how can I print out the nodal > positions in double precision? Right now I print these out to a file but I > only get 6 significant digits. http://en.cppreference.com/w/cpp/io/manip/setprecision or http://en.cppre

Re: [Libmesh-users] Mesh with subdomains problem

2016-04-28 Thread Pepijn Kessels
Thanks, that did the trick! Though I used add_node instead of add_side, and I used the nodal position to determine if the node should be flagged as boundary node. The stitch_meshes function then worked perfectly fine. One more question. If it is possible, how can I print out the nodal positions in

Re: [Libmesh-users] Mesh with subdomains problem

2016-04-28 Thread David Knezevic
As John said, the problem is presumably that your mesh doesn't have any boundary IDs in it. The mesh stitching uses the boundary IDs in order to figure out what needs to be stitched. You can try to add that data in yourself, if you like. You can do this within libMesh by looping over element bound

Re: [Libmesh-users] Mesh with subdomains problem

2016-04-28 Thread Pepijn Kessels
Thanks for the suggestions. I did try to stitch the meshes together using misc example 10, but it failed to find any boundary nodes. For the 2D mesh I used boundary_id_type 0, 1, 2, 3 for bottom, right, top and left, respectively. I note when loading the mesh it says "Mesh Dimension: 3" though, whi

Re: [Libmesh-users] Mesh with subdomains problem

2016-04-27 Thread John Peterson
On Wed, Apr 27, 2016 at 3:46 PM, Pepijn Kessels wrote: > Hi there, > > I'm trying to import a structured mesh for the 2D backward facing step > problem given on the nasa site. I transformed the original plot3d format > into exodus format using paraview. I could import the mesh into libMesh > with

Re: [Libmesh-users] Mesh with subdomains problem

2016-04-27 Thread David Knezevic
Have a look at miscellaneous_ex10. That stitches separate meshes together to form one mesh. It sounds very similar to what you're asking about, so maybe you can use that approach. David On Wed, Apr 27, 2016 at 4:46 PM, Pepijn Kessels wrote: > Hi there, > > I'm trying to import a structured mes

[Libmesh-users] Mesh with subdomains problem

2016-04-27 Thread Pepijn Kessels
Hi there, I'm trying to import a structured mesh for the 2D backward facing step problem given on the nasa site. I transformed the original plot3d format into exodus format using paraview. I could import the mesh into libMesh without problems, but I ran into some issues when trying to run a simula