Methods to Load Balance a GCR Pressure Solver Using a Stencil Framework on Multi- and Many-Core Architectures

Algorithm 2

Geometry geom = Geometry::init<

 DomainSize<NP,MP,LP>,

 HaloSize<HLEFT,HRIGHT,HBOT,

  HTOP,HGND,HSKY>,

 ProcessorGridSize<NPX,NPY,NPZ>>( );

Communicator comm =

 Communicator::init(geom);