g_ipt[T][LX+1][x2][x3] = i_odd + (LY*LZ); g_ipt[T+1][LX][x2][x3] = i_odd + 2*(LY*LZ);
g_ipt[T][LX+1][x2][x3] = i_even2; g_ipt[T+1][LX][x2][x3] = i_even2 + (LY*LZ);;
void set_ipt() {  int x0,x1,x2,x3;#ifndef _NEW_GEOMETRY  int ix=0;  for (x0 = 0; x0 < T; x0++) {    for (x1 = 0; x1 < LX; x1++) {      for (x2 = 0; x2 < LY; x2++) {	for (x3 = 0; x3 < LZ; x3++) {	  g_ipt[x0][x1][x2][x3] = ix;	  ix++;	}      }    }  }#if ((defined PARALLELT) || (defined PARALLELXT))  /* the time boundary */  for (x0 = T; x0 < T+2; x0++) {    for (x1 = 0; x1 < LX; x1++) {      for (x2 = 0; x2 < LY; x2++) {	for (x3 = 0; x3 < LZ; x3++) {	  g_ipt[x0][x1][x2][x3] = ix;	  ix++;	}      }    }  }#endif#if (defined PARALLELXT)  /* the x boundary */  for (x1 = LX; x1 < LX+2; x1++) {    for (x0 = 0; x0 < T; x0++) {      for (x2 = 0; x2 < LY; x2++) {	for (x3 = 0; x3 < LZ; x3++) {	  g_ipt[x0][x1][x2][x3] = ix;	  ix++;	}      }    }  }  /* The edges */  for (x0 = T; x0 < T+2; x0++) {    for (x1 = LX; x1 < LX+2; x1++) {      for (x2 = 0; x2 < LY; x2++) {	for (x3 = 0; x3 < LZ; x3++) {	  g_ipt[x0][x1][x2][x3] = ix;	  ix++;	}      }    }  }#endif#else  int i_even = 0,i_odd = VOLUME/2;  for (x0 = 0; x0 < T; x0++) {    for (x1 = 0; x1 < LX; x1++) {      for (x2 = 0; x2 < LY; x2++) {	for (x3 = 0; x3 < LZ; x3++) {	  if((x0+x1+x2+x3+g_proc_coords[0]*T+g_proc_coords[1]*LX)%2==0) {	    g_ipt[x0][x1][x2][x3] = i_even;	    i_even++;	  } 	  else {	    g_ipt[x0][x1][x2][x3] = i_odd;	    i_odd++;	  }	}      }    }  }#if ((defined PARALLELT) || (defined PARALLELXT))  i_even = VOLUME; i_odd = VOLUME+(LX*LY*LZ)/2;  /* the time boundary */  for (x1 = 0; x1 < LX; x1++) {    for (x2 = 0; x2 < LY; x2++) {      for (x3 = 0; x3 < LZ; x3++) {	if((x0+x1+x2+x3+g_proc_coords[0]*T+g_proc_coords[1]*LX)%2==0) {	  g_ipt[T][x1][x2][x3] = i_even;	  g_ipt[T+1][x1][x2][x3] = i_even+(LX*LY*LZ);	  i_even++;	} 	else {	  g_ipt[T][x1][x2][x3] = i_odd;	  g_ipt[T+1][x1][x2][x3] = i_odd+(LX*LY*LZ);	  i_odd++;	}      }    }  }#endif#if (defined PARALLELXT)  i_even = VOLUME+2*LX*LY*LZ; i_odd = i_even + (T*LY*LZ)/2;  /* the x boundary */  for (x0 = 0; x0 < T; x0++) {    for (x2 = 0; x2 < LY; x2++) {      for (x3 = 0; x3 < LZ; x3++) {	if((x0+x1+x2+x3+g_proc_coords[0]*T+g_proc_coords[1]*LX)%2==0) {	  g_ipt[x0][LX][x2][x3] = i_even;	  g_ipt[x0][LX+1][x2][x3] = i_even + (T*LY*LZ);	  i_even++;	} 	else {	  g_ipt[x0][LX][x2][x3] = i_odd;	  g_ipt[x0][LX+1][x2][x3] = i_odd + (T*LY*LZ);	  i_odd++;	}      }    }  }  i_even = VOLUME+RAND; i_odd = i_even+(LY*LZ)/2;  /* The edges */  for (x2 = 0; x2 < LY; x2++) {    for (x3 = 0; x3 < LZ; x3++) {      if((x0+x1+x2+x3+g_proc_coords[0]*T+g_proc_coords[1]*LX)%2==0) {	g_ipt[T][LX][x2][x3] = i_even;	g_ipt[T][LX+1][x2][x3] = i_even + (LY*LZ);	g_ipt[T+1][LX][x2][x3] = i_even + 2*(LY*LZ);	g_ipt[T+1][LX+1][x2][x3] = i_even + 3*(LY*LZ);	i_even++;      }       else {	g_ipt[T][LX][x2][x3] = i_odd;	g_ipt[T][LX+1][x2][x3] = i_odd + (LY*LZ);	g_ipt[T+1][LX][x2][x3] = i_odd + 2*(LY*LZ);	g_ipt[T+1][LX+1][x2][x3] = i_odd + 3*(LY*LZ);	i_odd++;      }    }  }#endif#endif}