Command line option --with-indices not working

[phuschke]

Problem

Neither

gmshreader --with-indices gmshMWE.msh

nor

gmshreader --local-ordering --with-indices gmshMWE.msh

gives the expected result.

MWE

gmshMWE.msh

$MeshFormat
2.2 0 8
$EndMeshFormat
$Nodes
5
1 0 0 0
2 1 0 0
3 1 1 0
4 0 1 0
5 0.5 0.5 0
$EndNodes
$Elements
4
1 2 5 7 6 2 2 -1 1 2 5
2 2 5 7 6 2 1 -2 1 5 4
3 2 5 7 6 2 2 -1 2 3 5
4 2 5 7 6 2 1 -2 3 4 5
$EndElements

[dbeurle]

Okay thanks for the issue. I get the following:

{
   "Elements" : [
      {
         "Name" : "",
         "NodalConnectivity" : [
            [ 1, 5, 4 ],
            [ 3, 4, 5 ]
         ],
         "Type" : 2
      }
   ],
   "Interface" : [
      {
         "GlobalStartId" : [ 0 ],
         "Master" : [ 1 ],
         "NodeIds" : [
            [ 1, 3, 5 ]
         ],
         "Slave" : [ 2 ],
         "Value" : [ 1 ]
      }
   ],
   "LocalToGlobalMap" : [ 1, 3, 4, 5 ],
   "Nodes" : [
      {
         "Coordinates" : [
            [ 0, 0, 0 ],
            [ 1, 1, 0 ],
            [ 0, 1, 0 ],
            [ 0.5, 0.5, 0 ]
         ]
      }
   ],
   "NumInterfaceNodes" : [ 3 ]
}

for the first, and for the second partition

{
   "Elements" : [
      {
         "Name" : "",
         "NodalConnectivity" : [
            [ 1, 2, 5 ],
            [ 2, 3, 5 ]
         ],
         "Type" : 2
      }
   ],
   "Interface" : [
      {
         "GlobalStartId" : [ 0 ],
         "Master" : [ 1 ],
         "NodeIds" : [
            [ 1, 3, 5 ]
         ],
         "Slave" : [ 2 ],
         "Value" : [ -1 ]
      }
   ],
   "LocalToGlobalMap" : [ 1, 2, 3, 5 ],
   "Nodes" : [
      {
         "Coordinates" : [
            [ 0, 0, 0 ],
            [ 1, 0, 0 ],
            [ 1, 1, 0 ],
            [ 0.5, 0.5, 0 ]
         ]
      }
   ],
   "NumInterfaceNodes" : [ 3 ]
}

With default arguments to the program. When I use --with-indices (with a patch) I get

{
   "Elements" : [
      {
         "Indices" : [ 2, 4 ],
         "Name" : "",
         "NodalConnectivity" : [
            [ 1, 5, 4 ],
            [ 3, 4, 5 ]
         ],
         "Type" : 2
      }
   ],
   "Interface" : [
      {
         "GlobalStartId" : [ 0 ],
         "Master" : [ 1 ],
         "NodeIds" : [
            [ 1, 3, 5 ]
         ],
         "Slave" : [ 2 ],
         "Value" : [ 1 ]
      }
   ],
   "LocalToGlobalMap" : [ 1, 3, 4, 5 ],
   "Nodes" : [
      {
         "Coordinates" : [
            [ 0, 0, 0 ],
            [ 1, 1, 0 ],
            [ 0, 1, 0 ],
            [ 0.5, 0.5, 0 ]
         ],
         "Indices" : [ 1, 3, 4, 5 ]
      }
   ],
   "NumInterfaceNodes" : [ 3 ]
}

and

{
   "Elements" : [
      {
         "Indices" : [ 1, 3 ],
         "Name" : "",
         "NodalConnectivity" : [
            [ 1, 2, 5 ],
            [ 2, 3, 5 ]
         ],
         "Type" : 2
      }
   ],
   "Interface" : [
      {
         "GlobalStartId" : [ 0 ],
         "Master" : [ 1 ],
         "NodeIds" : [
            [ 1, 3, 5 ]
         ],
         "Slave" : [ 2 ],
         "Value" : [ -1 ]
      }
   ],
   "LocalToGlobalMap" : [ 1, 2, 3, 5 ],
   "Nodes" : [
      {
         "Coordinates" : [
            [ 0, 0, 0 ],
            [ 1, 0, 0 ],
            [ 1, 1, 0 ],
            [ 0.5, 0.5, 0 ]
         ],
         "Indices" : [ 1, 2, 3, 5 ]
      }
   ],
   "NumInterfaceNodes" : [ 3 ]
}

Is this now the expected behaviour?

[phuschke]

Perfect 👍