<div class="gmail_quote">On Thu, Jun 16, 2011 at 17:57, Matthew Knepley <span dir="ltr">&lt;<a href="mailto:knepley@gmail.com">knepley@gmail.com</a>&gt;</span> wrote: <blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">
<div class="im"></div><div class="gmail_quote"><div class="im"><div><br></div></div><div>Doesn&#39;t he need A&#39; A (normal equations)</div></div></blockquote><div><br></div><div>Note that A has many more columns than rows, therefore A&#39;*A has huge null space.</div>
<div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;"><div class="gmail_quote"><div> for the solve, so a MatShell?</div></div></blockquote><div><br></div><div>I&#39;m sorry, yes. Or MatComposite:</div>
<div><br></div><div><a href="http://www.mcs.anl.gov/petsc/petsc-as/snapshots/petsc-dev/docs/manualpages/Mat/MatCreateComposite.html">http://www.mcs.anl.gov/petsc/petsc-as/snapshots/petsc-dev/docs/manualpages/Mat/MatCreateComposite.html</a></div>
<div><br></div><div>But there is also likely some regularization so MatShell makes sense.</div></div>