Hello again,
I have a few related questions, and help is, always, much appreciated!
I am using both pygmsh and meshpy to generate meshes and read them in, but
I'll focus on the former, pygmsh.
*Question 1*
For pygmsh, I can generate and read in (into libmesh) order-1 triangular
(and, I think, tetrahedral) meshes just fine, but when I generate order-2
elements, I receive an error that refers to an unrecognized type 'EDGE'
I've attached a small example of the output (points, cells and cell_data)
in pygmsh for a disk, which I save in exodus ('.e') format (using meshpy's
write function).
But in summary:
order-1 triangle elements:
'cells' is a dictionary with only 1 key, 'triangle'
order-2 triangle elements:
cells is a dictionary with keys:
cells['line3']
cells['triangle6']
cells['vertex']
Do I just need to change the key names to say, 'Edge3'?
*Question 2*
For the cases above, I noticed that I am no longer running libmesh in
parallel when I read in the mesh. I found, and tried to use
libMesh::Partitioner part();
// Partitioner()::partition_unpartitioned_elements(mesh);
// Partitioner.partition_unpartitioned_elements(mesh);
but I am not sure how exactly this works. Also, I am not sure if using
this the right approach?
*Question 3*
This is less important, as I may have sidestepped this for now, but I'm
sure it will return with a vengeance.
I looked through the forums to find out how VTKIO::write_nodal_data works,
and from there looked into the source for VTKIO::write_equation_system to
see if I could figure out how the former is used in the latter, but to no
avail. I'm writing using exo_io.write_timestep for now, but my screen
(rightly?) reprimands me for doing so.
I have a sense that this will involve a host of other functions
(allgather?) so I may send this as a separate email question.
------------------------------------------------------------------------------
Check out the vibrant tech community on one of the world's most
engaging tech sites, Slashdot.org! http://sdm.link/slashdot
_______________________________________________
Libmesh-users mailing list
Libmesh-users@lists.sourceforge.net
https://lists.sourceforge.net/lists/listinfo/libmesh-users
