/home/eoseret/qaas_runs_CPU_9468/171-148-3214/intel/CoMD/build/CoMD/CoMD/src-openmp/linkCells.c: 352 - 378
--------------------------------------------------------------------------------

352:    int ix = (int)(floor((rr[0] - localMin[0])*boxes->invBoxSize[0]));
353:    int iy = (int)(floor((rr[1] - localMin[1])*boxes->invBoxSize[1]));
354:    int iz = (int)(floor((rr[2] - localMin[2])*boxes->invBoxSize[2]));
355: 
356: 
357:    // For each axis, if we are inside the local domain, make sure we get
358:    // a local link cell.  Otherwise, make sure we get a halo link cell.
359:    if (rr[0] < localMax[0]) 
360:    {
361:       if (ix == gridSize[0]) ix = gridSize[0] - 1;
362:    }
363:    else
364:       ix = gridSize[0]; // assign to halo cell
365:    if (rr[1] < localMax[1])
366:    {
367:       if (iy == gridSize[1]) iy = gridSize[1] - 1;
368:    }
369:    else
370:       iy = gridSize[1];
371:    if (rr[2] < localMax[2])
372:    {
373:       if (iz == gridSize[2]) iz = gridSize[2] - 1;
374:    }
375:    else
376:       iz = gridSize[2];
377:    
378:    return getBoxFromTuple(boxes, ix, iy, iz);
