[petsc-users] PCDiagonalScale

[Klaij, Christiaan]

OK. So if I understand correctly (?) in Fortan all I need is:

call KSPGetPC
call PCDiagonalSet
call KSPSolve

I'm using a MatShell and PCShell but I guess that doesn't matter?


> 3) at line 253 I added :
>      call PCDiagonalScale(pc,PETSC_TRUE,ierr)
>

This is not the correct interface. You want PCDiagonalScaleSet(pc,X,ierr).
PCDiagonalScale() is a getter with petsc-3.1 and is not available from
Fortran (but you don't need it). The naming has been made consistent in
petsc-dev.


dr. ir. Christiaan Klaij
CFD Researcher
Research & Development
E mailto:C.Klaij at marin.nl
T +31 317 49 33 44

MARIN
2, Haagsteeg, P.O. Box 28, 6700 AA Wageningen, The Netherlands
T +31 317 49 39 11, F +31 317 49 32 45, I www.marin.nl