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

Algorithm 6

Communicator comm =

 Communicator::init(geom);

comm->update(p_in,sizeLeft,shiftLeft,

 sizeRight,shiftRight,updateLeft);

// or

comm->update_beg(p_in,sizeLeft,shiftLeft,

 sizeRight,shiftRight,updateLeft);

// kernel execution

comm->update_end(p_in,sizeLeft,shiftLeft,

 sizeRight,shiftRight,updateLeft);