/* { dg-do compile } */
/* { dg-require-effective-target size32plus } */
/* { dg-additional-options "-march=core-avx2" { target { x86_64-*-* i?86-*-* } } } */
/* { dg-require-effective-target size32plus } */

typedef enum {
 C = 0,               N, S, E, W, T, B,               NE, NW, SE, SW,               NT, NB, ST, SB,               ET, EB, WT, WB,               FLAGS, N_CELL_ENTRIES} CELL_ENTRIES;
typedef double LBM_Grid[(130)*100*100*N_CELL_ENTRIES];
void foo( LBM_Grid srcGrid )
{
  double ux , uy , uz , rho ,         ux1, uy1, uz1, rho1,         ux2, uy2, uz2, rho2,         u2, px, py;
  int i;
  for( i = 0;
       i < (N_CELL_ENTRIES*( 100*100));
       i += N_CELL_ENTRIES )
    {
      rho1 = + ((srcGrid)[((C)+N_CELL_ENTRIES*( 100*100))+(i)])
	  + ((srcGrid)[((N)+N_CELL_ENTRIES*( 100*100))+(i)])
	  + ((srcGrid)[((S)+N_CELL_ENTRIES*( 100*100))+(i)])
	  + ((srcGrid)[((E)+N_CELL_ENTRIES*( 100*100))+(i)])
	  + ((srcGrid)[((W)+N_CELL_ENTRIES*( 100*100))+(i)])
	  + ((srcGrid)[((T)+N_CELL_ENTRIES*( 100*100))+(i)])
	  + ((srcGrid)[((B)+N_CELL_ENTRIES*( 100*100))+(i)])
	  + ((srcGrid)[((NE)+N_CELL_ENTRIES*( 100*100))+(i)]) 
	  + ((srcGrid)[((NW)+N_CELL_ENTRIES*( 100*100))+(i)])
	  + ((srcGrid)[((SE)+N_CELL_ENTRIES*( 100*100))+(i)])
	  + ((srcGrid)[((SW)+N_CELL_ENTRIES*( 100*100))+(i)])
	  + ((srcGrid)[((NT)+N_CELL_ENTRIES*( 100*100))+(i)])
	  + ((srcGrid)[((NB)+N_CELL_ENTRIES*( 100*100))+(i)])
	  + ((srcGrid)[((ST)+N_CELL_ENTRIES*( 100*100))+(i)])
	  + ((srcGrid)[((SB)+N_CELL_ENTRIES*( 100*100))+(i)])
	  + ((srcGrid)[((ET)+N_CELL_ENTRIES*( 100*100))+(i)])
	  + ((srcGrid)[((EB)+N_CELL_ENTRIES*( 100*100))+(i)])
	  + ((srcGrid)[((WT)+N_CELL_ENTRIES*( 100*100))+(i)])
	  + ((srcGrid)[((WB)+N_CELL_ENTRIES*( 100*100))+(i)]);
      rho = 2.0*rho1 - rho2;
      px = (((i / N_CELL_ENTRIES) % 100) / (0.5*(100-1))) - 1.0;
      uz = 0.01 * (1.0-px*px) * (1.0-py*py);
      u2 = 1.5 * (ux*ux + uy*uy + uz*uz);
      (((srcGrid)[((C))+(i)])) = (1.0/ 3.0)*rho*(1.0 - u2);
      (((srcGrid)[((N))+(i)])) = (1.0/18.0)*rho*(1.0 + uy*(4.5*uy + 3.0) - u2);
    }
}