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