[petsc-users] Inverting a Large Matrix with Fortran Code

Hello, I have a problem of the kind: G(i) = (invM.P).G(i+1) +invM.(d -R.b), where M, P and R are square tri-diagonal matrices, invM is the inverse of M, G(i) and G(i+1) are the linearised solution vectors (given G(i+1), G(i) is determined using above), and d and b are vectors. The problem

Re: [petsc-users] Inverting a Large Matrix with Fortran Code

Koki Imada koki.im...@york.ac.uk writes: Hello, I have a problem of the kind: G(i) = (invM.P).G(i+1) +invM.(d -R.b), where M, P and R are square tri-diagonal matrices, invM is the inverse of M, G(i) and G(i+1) are the linearised solution vectors (given G(i+1), G(i) is determined