Hello,
It seems there is a small bug in the routine for reading in a GMSH mesh. If the
mesh has been saved in Version 2.0 format, and the mesh has Physical
Points/Physical Lines/Physical Surfaces, etc defined then the following error
is generated:
terminate called after throwing an instance of
'dealii::GridIn::ExcInvalidGMSHInput'
what(): --------------------------------------------------------
An error occurred in line of file in function
void dealii::GridIn::read_msh(std::istream&) [with int dim = 2, int
spacedim = 2]
The violated condition was:
line == "$Nodes"
The name and call sequence of the exception was:
ExcInvalidGMSHInput(line)
Additional Information:
The string is not recognized at the present position of a Gmsh Mesh file.
--------------------------------------------------------
..the error appears duw to the following block appearsing before the $Nodes
flag in the .msh file:
$PhysicalNames
2
1 1 "1"
2 2 "This surface"
$EndPhysicalNames
....in the grid_in.read_msh function, the "$Nodes" flag is read just after the
"$EndMeshFormat" flag, causing the error.
of course this can be worked arround by:
(1) deleting the $PhysicalNames block
(2) saving in version-1 format :)
Hope it helps somebody....
Abhishek_______________________________________________
dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
