I am having issues when using shared memory arrays in fortran. If someone can look at the code and identify the errors in the code, that would be great.
I uploaded the files to mediafire.
In folder v1, I didn’t use shared memory and the code is working fine.
To get more performance, I changed the structure of code and used 2 dimensional blocks and with that I started using shared memory. But the code is giving garbage results. Logically the code looks sound, but I am making some mistake unknown to me. The modified code is in folder v2
The troublesome files are xi.90 and enrgy_eval.f90. And when I am mixing xi and enrgy_eval files from the two folders, I am getting memory errors.
Any help is greatly appreciated