1+ #include <assert.h>
12#include <stdio.h>
2- #include <stdlib.h>
3- #include <math.h>
4- #include <unistd.h>
53#include <sys/time.h>
6- #include <assert .h>
4+ #include <unistd .h>
75
86#define nx 300
97#define ny 300
108#define nz 300
119// #define T 20
1210// #define 9 9
1311
14- #define f60 0.2
12+ #define f60 0.2
1513#define f61 0.5
1614#define f62 0.3
1715
@@ -85,7 +83,7 @@ int main()
8583 init_array () ;
8684
8785#ifdef PERFCTR
88- PERF_INIT ;
86+ PERF_INIT ;
8987#endif
9088
9189 IF_TIME (t_start = rtclock ());
@@ -95,24 +93,26 @@ int main()
9593 for (j = 4 ; j <= ny + 9 - 2 ; j ++ )
9694 for (i = 4 ; i <= nx + 9 - 3 ; i ++ )
9795 for (k = 4 ; k <= nz + 9 - 3 ; k ++ )
98- ab [j ][i ][k ] = (f60 * (a [j - 1 ][i ][k ] + a [j ][i ][k ]) + f61
99- * (a [j - 2 ][i ][k ] + a [j + 1 ][i ][k ]) + f62 * (a [j - 3 ][i ][k ]
100- + a [j + 2 ][i ][k ])) * thirddtbydy * uyb [j ][i ][k ];
96+ ab [j ][i ][k ] = (f60 * (a [j - 1 ][i ][k ] + a [j ][i ][k ]) +
97+ f61 * (a [j - 2 ][i ][k ] + a [j + 1 ][i ][k ]) +
98+ f62 * (a [j - 3 ][i ][k ] + a [j + 2 ][i ][k ])) *
99+ thirddtbydy * uyb [j ][i ][k ];
101100
102101 for (j = 4 ; j <= ny + 9 - 3 ; j ++ )
103102 for (i = 4 ; i <= nx + 9 - 2 ; i ++ )
104103 for (k = 4 ; k <= nz + 9 - 3 ; k ++ )
105- al [j ][i ][k ] = (f60 * (a [j ][i - 1 ][k ] + a [j ][i ][k ]) + f61
106- * (a [j ][i - 2 ][k ] + a [j ][i + 1 ][k ]) + f62 * (a [j ][i - 3 ][k ] + a [j ][i + 2 ][k ]))
107- * thirddtbydx * uxl [j ][i ][k ];
104+ al [j ][i ][k ] = (f60 * (a [j ][i - 1 ][k ] + a [j ][i ][k ]) +
105+ f61 * (a [j ][i - 2 ][k ] + a [j ][i + 1 ][k ]) +
106+ f62 * (a [j ][i - 3 ][k ] + a [j ][i + 2 ][k ])) *
107+ thirddtbydx * uxl [j ][i ][k ];
108108
109109 for (j = 4 ; j <= ny + 9 - 3 ; j ++ )
110110 for (i = 4 ; i <= nx + 9 - 3 ; i ++ )
111111 for (k = 4 ; k <= nz + 9 - 2 ; k ++ )
112- af [j ][i ][k ] = (f60 * (a [j ][i ][k - 1 ] + a [j ][i ][k ]) + f61
113- * ( a [ j ][ i ][ k - 2 ] + a [ j ][ i ][ k + 1 ]) + f62 * (a [j ][i ][k - 3 ] + a [j ][i ][k + 2 ]))
114- * thirddtbydz * uzf [j ][i ][k ];
115-
112+ af [j ][i ][k ] = (f60 * (a [j ][i ][k - 1 ] + a [j ][i ][k ]) +
113+ f61 * (a [j ][i ][k - 2 ] + a [j ][i ][k + 1 ]) +
114+ f62 * ( a [j ][i ][k - 3 ] + a [ j ][ i ][ k + 2 ])) *
115+ thirddtbydz * uzf [ j ][ i ][ k ];
116116
117117 for (j = 4 ; j <= ny + 9 - 3 ; j ++ )
118118 for (i = 4 ; i <= nx + 9 - 3 ; i ++ )
@@ -125,7 +125,7 @@ int main()
125125 IF_TIME (fprintf (stderr , "%0.6lfs\n" , t_end - t_start ));
126126
127127#ifdef PERFCTR
128- PERF_EXIT ;
128+ PERF_EXIT ;
129129#endif
130130
131131#ifdef TEST
0 commit comments