I've got a very weird problem. I used the vdsqr function in a loop with k as index, the program was like this:
do k = 1, kk
call vdSQR(32, XI(K,:), RSQ(K, :))
kk is a very large number, xi and rsq are two matrix with kk rows and 32 columns.
The problem comes when k reaches a certain number, say 10189. At this iteration, after the call, the xi will change its value, and k will become 0. The other iterations when k is less than this number seem to be right.
Also, if I use l instead of k, say,
do l = 1, kk
the problem doesn't show up.
I'm just curious is there any reported problem like this, or is it just my programing mistakes? And any hint on what happend? Thank you very much!