10 double t_start
, t_end
;
23 IF_TIME(t_start
= rtclock());
25 /* pluto start (T,N) */
26 for (t
=0; t
<=T
-1; t
++) {
27 for (i
=1; i
<=N
-2; i
++) {
28 for (j
=1; j
<=N
-2; j
++) {
29 a
[i
][j
] = (a
[i
-1][j
-1] + a
[i
-1][j
] + a
[i
-1][j
+1]
30 + a
[i
][j
-1] + a
[i
][j
] + a
[i
][j
+1]
31 + a
[i
+1][j
-1] + a
[i
+1][j
] + a
[i
+1][j
+1])/9.0;
37 IF_TIME(t_end
= rtclock());
38 IF_TIME(fprintf(stderr
, "%0.6lfs\n", t_end
- t_start
));