4 extern void abort (void);
6 subdomain (float *x
, int istart
, int ipoints
)
9 for (i
= 0; i
< ipoints
; i
++)
10 x
[istart
+ i
] = 123.456;
14 sub (float *x
, int npoints
)
16 int iam
, nt
, ipoints
, istart
;
17 #pragma omp parallel default(shared) private(iam,nt,ipoints,istart)
19 iam
= omp_get_thread_num ();
20 nt
= omp_get_num_threads ();
21 ipoints
= npoints
/ nt
; /* size of partition */
22 istart
= iam
* ipoints
; /* starting array index */
23 if (iam
== nt
- 1) /* last thread may do more */
24 ipoints
= npoints
- istart
;
25 subdomain (x
, istart
, ipoints
);
34 for (i
= 0; i
< 10000; i
++)
35 if (array
[i
] < 123.45 || array
[i
] > 123.46)