1*404b540aSrobert /* { dg-do run } */
2*404b540aSrobert
3*404b540aSrobert #include <omp.h>
4*404b540aSrobert extern void abort (void);
5*404b540aSrobert void
subdomain(float * x,int istart,int ipoints)6*404b540aSrobert subdomain (float *x, int istart, int ipoints)
7*404b540aSrobert {
8*404b540aSrobert int i;
9*404b540aSrobert for (i = 0; i < ipoints; i++)
10*404b540aSrobert x[istart + i] = 123.456;
11*404b540aSrobert }
12*404b540aSrobert
13*404b540aSrobert void
sub(float * x,int npoints)14*404b540aSrobert sub (float *x, int npoints)
15*404b540aSrobert {
16*404b540aSrobert int iam, nt, ipoints, istart;
17*404b540aSrobert #pragma omp parallel default(shared) private(iam,nt,ipoints,istart)
18*404b540aSrobert {
19*404b540aSrobert iam = omp_get_thread_num ();
20*404b540aSrobert nt = omp_get_num_threads ();
21*404b540aSrobert ipoints = npoints / nt; /* size of partition */
22*404b540aSrobert istart = iam * ipoints; /* starting array index */
23*404b540aSrobert if (iam == nt - 1) /* last thread may do more */
24*404b540aSrobert ipoints = npoints - istart;
25*404b540aSrobert subdomain (x, istart, ipoints);
26*404b540aSrobert }
27*404b540aSrobert }
28*404b540aSrobert int
main()29*404b540aSrobert main ()
30*404b540aSrobert {
31*404b540aSrobert int i;
32*404b540aSrobert float array[10000];
33*404b540aSrobert sub (array, 10000);
34*404b540aSrobert for (i = 0; i < 10000; i++)
35*404b540aSrobert if (array[i] < 123.45 || array[i] > 123.46)
36*404b540aSrobert abort ();
37*404b540aSrobert return 0;
38*404b540aSrobert }
39