cylinder example: gmsh2sod2d.py not working in converting msh to h5

[soaringxmc]

Hi, Bernat:

It seems gmsh2sod2d.py contained in sod2d_gitlab master branch commit 433356db does not work for the cylinder example.

gmsh version 4.12.1

python3 /home/mchao/code/sod2d_gitlab/utils/gmsh2sod2d/gmsh2sod2d.py cyl_p4_psenv -p 1

The error is
File "/home/mchao/code/sod2d_gitlab/utils/gmsh2sod2d/gmsh2sod2d.py", line 313, in
raiseError('Error! len(conec)!=npbou in boundary element!!')

Could you please resolve the issue? I think the cylinder case will be a very good simple example to test the coupling of the framework with SOD2D.

[b-fg]

Is your SOD2D up to date? The error suggests a problem on the setup of the boundary conditions in GMsh, but I assume you have used the example .geo file. @RicardMontala is the .geo up to date with the latest SOD2D version?

Also, maybe is better to directly open issues in SOD2D for this sort of trouble. This library (SmartSOD2D) is for the coupling, not the solver itself.

[RicardMontala]

Hello,

Just to be sure, could you try generating your mesh using GMSH version 4.10.3? Sometimes, I have experimented some issues regarding the gmsh version.

[soaringxmc]

The issue has been resolved. Thank you!