fix Grid rank checking

Author: jxy

src/grid/GridDefs.nim

@@ -102,6 +102,11 @@ proc newGridLatticeGaugeField*(x: ptr GridCartesian):
 template gauge*(x: GridCartesian): untyped =
   newGridLatticeGaugeField(unsafeaddr x)
 
+proc ThisRank*(x: GridBase): int =
+  var r = 0
+  {.emit:[r," = ",x,"->ThisRank();"].}
+  r
+
 proc lSites*(x: ptr GridBase): cint {.importcpp,gh.}
 proc RankIndexToGlobalCoor*(g: GridBase, rank,oixd,iidx: SomeInteger; coord: Coordinate) {.importcpp.}
 

src/grid/GridUtils.nim

@@ -8,13 +8,13 @@ template getGridPtr*(g: Field): auto =
   let simd_layout = GridDefaultSimd(lo.nDim, Nsimd(GridVComplex))
   let mpi_layout = newCoordinate(lo.rankGeom)
   let grid = newGridCartesian(latt_size,simd_layout,mpi_layout)
-  var coor = newSeq[cint](lo.nDim)
-  var gcoor = newCoordinate(coor)
-  grid.RankIndexToGlobalCoor(lo.myrank, 0, 0, gcoor)
-  coor := gcoor
-  let ri = lo.rankIndex(coor)
-  #echoAll lo.myrank, ": ", coor
-  if ri.rank != lo.myrank:
+  #var coor = newSeq[cint](lo.nDim)
+  #var gcoor = newCoordinate(coor)
+  #grid.RankIndexToGlobalCoor(lo.myrank, 0, 0, gcoor)
+  #coor := gcoor
+  #let ri = lo.rankIndex(coor)
+  #echoAll lo.myrank, ": ThisRank=",grid.ThisRank, ": ", coor
+  if grid.ThisRank != lo.myrank:
     qexFatal "Grid rank and QEX rank disagree"
   unsafeAddr grid