#include $$stdio$.$h$>$

#define N $1000$

"define THREAD$\_$NO $10$

int main$()\{$

double $a[N],b[N],c[N],d[N]$;

int $i$;

$/*$ some instructions that initialize arrays a$,$ b$,$ c$,$ d

take place here. $*/$

#pragma omp parallel for num$\_$threads$($THREAD$\_$NO$)$

for $b[2^{**}i]+=a[2^{**}i]+d[2^{**}i]$;$c[i]+=b[2**i+1]-a[2**i](i=0$;$i$

$b[2^{**}i]+=a[2^{**}i]+d[2^{**}i]$;

$c[i]+=b[2**i+1]-a[2**i]$;

$\}$