[petsc-users] Modification of a routine

[Gaurish Telang]

Hi,

I have tried to implement a recent least squares algorithm called LSMR, by
making modifications in the file  src/ksp/ksp/impls/lsqr/lsqr.c

Is it necessary to make any changes in other PETSc files or build the PETSc
library again?

I solved a simple least squares problem by supplying the flags -ksp_type
lsqr  -pc_type none and the problem seems to get solved correctly.

However, I had placed a couple of PetscPrintf statements inside the main
do-while loop of the algorithm in lsqr.c  but those statements are not
getting printed to standard output.

Thanks,

Gaurish.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20110404/75c15069/attachment.htm>