<div class="gmail_quote">On Thu, Jun 16, 2011 at 17:34, Matthew Knepley <span dir="ltr">&lt;<a href="mailto:knepley@gmail.com">knepley@gmail.com</a>&gt;</span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">
<div>Yes, you can do MF applications of A&#39; A. You could just make a MatShell that called MatMult and MatMultTranspose.</div></blockquote><div><br></div><div>MatShell is not necessary, just use MatCreateTranspose().</div>
<div><br></div><div><a href="http://www.mcs.anl.gov/petsc/petsc-as/snapshots/petsc-dev/docs/manualpages/Mat/MatCreateTranspose.html">http://www.mcs.anl.gov/petsc/petsc-as/snapshots/petsc-dev/docs/manualpages/Mat/MatCreateTranspose.html</a></div>
<div><br></div><div>But I think Alexander is clear on that part.</div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;"><div></div><div>For modifying the matrix, MatGetRow() IF you are using an AIJ matrix, return a pointer directly to the values.</div>
</blockquote><div><br></div><div>This is only for SeqAIJ.</div><div><br></div><div><br></div><div>Alexander, could you clarify where in the Gauss-Newton procedure you outlined is the dense matrix on which you need to do the nonlinear row operation?</div>
<div><br></div><div>Also, the matrix A&#39;*A has a huge null space. What sort of regularization are you using.</div></div>