Pastebin

  if (InsideMesh) {
    for (k = istart[2]; k <= iend[2]; k++) {
      k1 = (GridRank > 2) ? k % Mesh->Dimension[2] : k;
      for (j = istart[1]; j <= iend[1]; j++) {
	j1 = (GridRank > 1) ? j % Mesh->Dimension[1] : j;
	index = (k1*Mesh->Dimension[1] + j1)*Mesh->Dimension[0];
	for (i = istart[0]; i <= iend[0]; i++) {
	  i1 = i % Mesh->Dimension[0];
	  new_link = new ChainingMeshLink;
	  new_link->NextLink = Mesh->HeadOfChain[index+i1];
	  new_link->GridData = this;
	  Mesh->HeadOfChain[index+i1] = new_link;
	}
      }
    }
  }