rcce/bt/header.h

49 #define JMAX MAX_CELL_DIM
68 #define us(i,j,k,c) us[(i+1)+(IMAX+2)*((j+1)+(JMAX+2)*((k+1)+(KMAX+2)*(c-1)))]
69 #define vs(i,j,k,c) vs[(i+1)+(IMAX+2)*((j+1)+(JMAX+2)*((k+1)+(KMAX+2)*(c-1)))]
70 #define ws(i,j,k,c) ws[(i+1)+(IMAX+2)*((j+1)+(JMAX+2)*((k+1)+(KMAX+2)*(c-1)))]
71 #define qs(i,j,k,c) qs[(i+1)+(IMAX+2)*((j+1)+(JMAX+2)*((k+1)+(KMAX+2)*(c-1)))]
72 #define rho_i(i,j,k,c) rho_i[(i+1)+(IMAX+2)*((j+1)+(JMAX+2)*((k+1)+(KMAX+2)*(c-1)))]
73 #define square(i,j,k,c) square[(i+1)+(IMAX+2)*((j+1)+(JMAX+2)*((k+1)+(KMAX+2)*(c-1)))]
74 #define forcing(m,i,j,k,c) forcing[(m-1)+5*(i+IMAX*(j+JMAX*(k+KMAX*(c-1))))]
75 #define u(m,i,j,k,c) u[(m-1)+5*((i+2)+(IMAX+4)*((j+2)+(JMAX+4)*((k+2)+(KMAX+4)*(c-1))))]
76 #define rhs(m,i,j,k,c) rhs[(m-1)+5*((i+1)+(IMAX+1)*((j+1)+(JMAX+1)*((k+1)+(KMAX+1)*(c-1))))]
77 #define lhsc(m,n,i,j,k,c) lhsc[(m-1)+5*((n-1)+5*((i+1)+(IMAX+1)*((j+1)+(JMAX+1)*((k+1)+(KMAX+1)*(c-1)))))]
78 #define backsub_info(m,i,j,c) backsub_info[(m-1)+5*((i)+(IMAX+1)*((j)+(JMAX+1)*(c-1)))]
119          us      [(IMAX+2)*(JMAX+2)*(KMAX+2)*MAXCELLS],
120          vs      [(IMAX+2)*(JMAX+2)*(KMAX+2)*MAXCELLS],
121          ws      [(IMAX+2)*(JMAX+2)*(KMAX+2)*MAXCELLS],
122          qs      [(IMAX+2)*(JMAX+2)*(KMAX+2)*MAXCELLS],
123          rho_i   [(IMAX+2)*(JMAX+2)*(KMAX+2)*MAXCELLS],
124          square  [(IMAX+2)*(JMAX+2)*(KMAX+2)*MAXCELLS],
125          forcing [5*IMAX*JMAX*KMAX*MAXCELLS],
126          u       [5*(IMAX+4)*(JMAX+4)*(KMAX+4)*MAXCELLS],
127          rhs     [5*(IMAX+1)*(JMAX+1)*(KMAX+1)*MAXCELLS],
128          lhsc    [5*5*(IMAX+1)*(JMAX+1)*(KMAX+1)*MAXCELLS],
178          us      [(IMAX+2)*(JMAX+2)*(KMAX+2)*MAXCELLS],
179          vs      [(IMAX+2)*(JMAX+2)*(KMAX+2)*MAXCELLS],
180          ws      [(IMAX+2)*(JMAX+2)*(KMAX+2)*MAXCELLS],
181          qs      [(IMAX+2)*(JMAX+2)*(KMAX+2)*MAXCELLS],
182          rho_i   [(IMAX+2)*(JMAX+2)*(KMAX+2)*MAXCELLS],
183          square  [(IMAX+2)*(JMAX+2)*(KMAX+2)*MAXCELLS],
184          forcing [5*IMAX*JMAX*KMAX*MAXCELLS],
185          u       [5*(IMAX+4)*(JMAX+4)*(KMAX+4)*MAXCELLS],
186          rhs     [5*(IMAX+1)*(JMAX+1)*(KMAX+1)*MAXCELLS],
187          lhsc    [5*5*(IMAX+1)*(JMAX+1)*(KMAX+1)*MAXCELLS],