Hi Chris,
On 9/9/22 20:48, Christopher Walker wrote:
Hi Robert, Many thanks for your reply. I tried the commands shown at the top of the convert_mesh.py and it worked fine. e.g. python ./script/convert_mesh.py meshes/3d/cylinder.mesh new.vtk sfepy: reading mesh (meshes/3d/cylinder.mesh)... sfepy: number of vertices: 354 sfepy: number of cells: sfepy: 3_4: 1348 sfepy: ...done in 0.02 s sfepy: writing new.vtk... sfepy: ...done
However, when I tried with the screwdriver example, I get this: python ./script/convert_mesh.py -d 3 screwdriver_handle.msh screwdriver_handle.vtk sfepy: reading mesh (screwdriver_handle.msh)... Traceback (most recent call last): File "C:\Users\ChrisWalker\anaconda3\Lib\site-packages\sfepy\script\convert_mesh.py", line 237, in <module> main() File "C:\Users\ChrisWalker\anaconda3\Lib\site-packages\sfepy\script\convert_mesh.py", line 162, in main mesh = Mesh.from_file(filename_in) File "C:\Users\ChrisWalker\anaconda3\lib\site-packages\sfepy\discrete\fem\mesh.py", line 220, in from_file mesh = io.read(mesh, omit_facets=omit_facets) File "C:\Users\ChrisWalker\anaconda3\lib\site-packages\sfepy\discrete\fem\meshio.py", line 421, in read cgdata0[m.cell_sets[k][ic]] = cgrps[ik] IndexError: index 5 is out of bounds for axis 0 with size 2
I am in no rush, so maybe its best to wait until the end of September for the new updates, but I thought it best to mention it in case you think it is something you need to check before the next release. If you want me to send my version of screwdriver_handle.msh then I can do so.
Thanks for the report, I can reproduce the issue. You can try generate another format, e.g.:
gmsh -3 -format mesh -o screwdriver_handle.mesh screwdriver_handle.geo
Then the subsequent conversion (stripping lower dimensional entities) works for me.
Cheers, r.