Description

According to the documentation

http://www.magpar.net/magpar/doc/html/preprocessing.html#gmsh

the "Physical Volumes" in Gmsh .geo files should be numbered consecutively starting from 0.

However, this is an error in the documentation. Their numbering of the "Physical Volumes" has to start with 1.

Bug reported by Rajmani Gaur - thank you!

Steps to reproduce

1. See attached file cylinder.geo

2. Mesh it with "gmsh -3 cylinder.geo", which gives cylinder.msh
3. "gmshtoucd cylinder.msh" gives cylinder.inp
4. Run simulation with cylinder.inp.

Example

magpar output/error message:

[...]
Opening file cylinder.inp
Number of grid vertices: 2937
Number of grid elements: 13696
Number of vertex data:   3
Number of element data:  0
Number of element data2: 0
Allocating array for vertices and reading their coordinates
Allocating array for element connectivity and reading it
magpar.exe: io/readinp.c:220: int ReadINP(GridData*, char*, _p_Vec*, int):
Assertion `a1>0' failed.
Aborted

Details

Workaround/Fix

Make sure "Physical Volumes" in .geo files are numbered consecutively starting from 1.

Here is a fixed version of the Gmsh file: cylinder_fixed.geo

Plan

• Priority: medium

• Assigned to: WernerScholz

• Status: fixed: documentation fixed in development version leading to 0.9

Attachments

There are 2 attachment(s) stored for this page.

CategoryMagparBugFixed