Fortran TDMA algorithm on wiki -
i visited wiki site:
it says a,b,c sub-diagonal, diagonal, , super-diagonal. if 'n' size of b, aren't sizes of 'a' , 'c' n-1? clear inspection both c(1) , c(n) accessed, contradicts sizes i've described.
what's going on here? algorithm wrong? matrix sizes i've assumed wrong? input matrix algorithm?
any appreciated!
cp(n) calculated avoid awkward if statement - never used in determining x. c(n) not used.
a(1) unused, expected.
Comments
Post a Comment