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

Algorithm 3

struct LaplcStencil

  DEFINE_DO(const T _ _restrict_ _ in_p,

  T _ _restrict_ _ out_p,const

  T &cCoeff)

 IN3D(out_p, 0, 0, 0) = cCoeff   (

  IN3D(in_p, -1, 0, 0, CACHED) +

  IN3D(in_p, 1, 0, 0, CACHED) +

  IN3D(in_p, 0, -1, 0, CACHED) +

  IN3D(in_p, 0, 1, 0, CACHED) +

  IN3D(in_p, 0, 0, -1, CACHED) +

  IN3D(in_p, 0, 0, 1, CACHED));

;