Make the following code work.

#pragma omp parallel \

private(myid,istart,iend,nthreads,nper)

nthreads = omp_get_num_threads();

nper = imax/nthreads;

myid = omp_get_thread_num();

istart = myid*nper;

iend = istart + nper 1;

do_work(istart,iend);

for(i=istart; i<=iend; i++){

c[i] = a[i]*b[i] + ...

}