Hi,<br><br>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 <br><br>Is it necessary to make any changes in other PETSc files or build the PETSc library again? <br>
<br>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. <br><br>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. <br>
<br>Thanks,<br><br>Gaurish. <br>