On Fri, 8 Nov 2013, at 10:09, Svetlana Tkachenko wrote: > ex2f.F involves a band matrix with 4 on main diagonal and -1s on the > sub-diagonals. I have A_central matrix which stores the -4s, and A_left, > A_right, Atop, A_bottom matrices which store the -1s (they actually are more > complicated numbers than that). Would it be recommended to simply wrap the > "do 10, II=..." block like the one in ex2f.F into $OMP DO to speed up filling > of the matrix as each of the values set can be set independently of others? > (I tried doing so but got a segmentation fault and I'm working on a small > example to reproduce the problem.) > > Svetlana
I figured out this is called 'stencils' correctly and corrected the subject. Trying to use example ex22f.F but it's quite involved. Svetlana
