Hi Mat,
aj(im1)=0.d0
bj(im1)=0.d0
cj(im1)=0.d0
dj(im1)=0.d0
aj1=0.d0
bj1=0.d0
cj1=0.d0
dj1=0.d0
you are right, those matrices I put them there. I thought it’s a good idea to initialize. you could remove them. if you remove them, still the code does not work and the akn,akd will equal to zero with acc directives. Please advice.
thanks,
Dolf