11#include <stdio.h>
2- #include <stdlib.h>
3- #include <math.h>
42#include <unistd.h>
53#include <sys/time.h>
64#include <assert.h>
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,21 +93,21 @@ 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 ]
96+ ab [j ][i ][k ] = (f60 * (a [j - 1 ][i ][k ] + a [j ][i ][k ]) + f61
97+ * (a [j - 2 ][i ][k ] + a [j + 1 ][i ][k ]) + f62 * (a [j - 3 ][i ][k ]
10098 + a [j + 2 ][i ][k ])) * thirddtbydy * uyb [j ][i ][k ];
10199
102100 for (j = 4 ; j <= ny + 9 - 3 ; j ++ )
103101 for (i = 4 ; i <= nx + 9 - 2 ; i ++ )
104102 for (k = 4 ; k <= nz + 9 - 3 ; k ++ )
105- al [j ][i ][k ] = (f60 * (a [j ][i - 1 ][k ] + a [j ][i ][k ]) + f61
103+ al [j ][i ][k ] = (f60 * (a [j ][i - 1 ][k ] + a [j ][i ][k ]) + f61
106104 * (a [j ][i - 2 ][k ] + a [j ][i + 1 ][k ]) + f62 * (a [j ][i - 3 ][k ] + a [j ][i + 2 ][k ]))
107105 * thirddtbydx * uxl [j ][i ][k ];
108106
109107 for (j = 4 ; j <= ny + 9 - 3 ; j ++ )
110108 for (i = 4 ; i <= nx + 9 - 3 ; i ++ )
111109 for (k = 4 ; k <= nz + 9 - 2 ; k ++ )
112- af [j ][i ][k ] = (f60 * (a [j ][i ][k - 1 ] + a [j ][i ][k ]) + f61
110+ af [j ][i ][k ] = (f60 * (a [j ][i ][k - 1 ] + a [j ][i ][k ]) + f61
113111 * (a [j ][i ][k - 2 ] + a [j ][i ][k + 1 ]) + f62 * (a [j ][i ][k - 3 ] + a [j ][i ][k + 2 ]))
114112 * thirddtbydz * uzf [j ][i ][k ];
115113
@@ -125,7 +123,7 @@ int main()
125123 IF_TIME (fprintf (stderr , "%0.6lfs\n" , t_end - t_start ));
126124
127125#ifdef PERFCTR
128- PERF_EXIT ;
126+ PERF_EXIT ;
129127#endif
130128
131129#ifdef TEST
0 commit comments