Hi,<br><br>I am planning to implement the LSMR algorithm in PETSc which does least squares and is supposed to have more favourable mathematical properties than LSQR which has already been implemented in PETSc.  <br><br>But I am not really sure how to go about this,  since the implementations of the standard KSP methods themselves look quite complicated. The manual does not seeem to say much <br>
about how to go about  adding routines to the PETSc library.<br><br>Could you give me guidelines I should follow while implementing a KSP method? <br><br>Also is it necessary to build the PETSc library again after implementing this routine. If so is it necessary to make any changes to makefiles ?<br>
<br>Regards<br>