Hi all,

I've post some days ago about efficient methods for solving poisson equations, and someone advice me to try the Stone SIP algorithm, as described in Peric book.

I've read about it, implement it in c++, but it seems to me that I've made a mistake somewhere and I can't find it. I've also compare my code with Peric's one, and I still didn't found the error.

I've doubts about my incomplete LU implementation, but there's no validation case to compare with...

So can someone help me?

i can only say this that one source of mistake during swiching between C++ and FOR is related to dimensions (formaly index is started from 1 in for and 0 in C++) it is especially important as Peric used vectorization (converting 2d array in 1D)

e.g. for 2d case

index (i,j)= i+(j-1)*imax in fortran index (i,j)= i+j*imax in C++

hope it help.

Are you in my class

NSPDE-1