Hi Hong,<br><br>On Mon, Jul 23, 2012 at 11:53 AM, Hong Zhang <span dir="ltr"><<a href="mailto:hzhang@mcs.anl.gov" target="_blank">hzhang@mcs.anl.gov</a>></span> wrote:<br><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Jack:<div class="gmail_quote"><div class="im"><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><br>I'm getting started on creating documentation right now; it should end up looking very similar to Elemental's.<br>
</blockquote><div> </div></div><div>The detailed instruction will be very helpful, especially</div><div>the sparse matrix data structures. Notify us when the documentation is out.</div></div></blockquote><div><br>Will do. I just pushed out documentation for the build system:<br>
<a href="http://poulson.github.com/Clique/build.html">http://poulson.github.com/Clique/build.html</a><br><br>The rest should follow soon.<br> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div class="gmail_quote"><div>Is this a sparse Cholesky direct solver?</div></div></blockquote><div><br>More or less. It currently supports LDL^T and LDL^H factorizations without pivoting, and the latter certainly works for Hermitian positive-definite matrices. The current implementation is slightly more general than Cholesky and should achieve roughly the same performance.<br>
</div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div class="gmail_quote">
<div>Currently, only mumps supports parallel sparse Cholesky direct solver.</div><div>Yours will be a good addition to the petsc external solvers.</div></div></blockquote><div><br>Thanks! I'll do my best to help with the integration when the time comes.<br>
</div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div class="gmail_quote"><div class="im"><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<br>Please keep in mind that there is not _yet_ any support for pivoting, as I needed Clique's functionality for indefinite complex symmetric matrices which are nice enough to be factored accurately without pivoting. </blockquote>
<div> </div></div><div>Petsc users likely use it as a preconditioner. </div><div>In case of zero pivot, you may provide a routine/option for adding a small</div><div>shift? </div></div></blockquote><div><br>Yes, that is what I was thinking. For instance, for Helmholtz equations, one can get away with very small positive imaginary shifts (the physically corresponds to slightly damping waves).<br>
</div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div class="gmail_quote"><div class="im"><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<br>Clique's nested dissection is built on top of a nodal graph bisection routine which uses a custom interface to parmetis; I suspect that the graph partitioning is currently the least scalable part of the black-box interface due to KLFM refinement being so hard to parallelize.</blockquote>
<div><br></div></div><div>Xuan has almost finished the petsc-elemental interface</div><div>(see petsc-dev/src/mat/impls/elemental/).</div><div>We would appreciate if you can take a look at it and give us</div><div>your comments/suggests. </div>
<div class="im">
<div><br></div></div></div></blockquote><div><br>I just took a quick look and will make time to go over it in detail later today. I'll send out patch if I have any concrete suggestions.<br><br>Jack<br><br></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div class="gmail_quote"><div class="im"><div></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><span><font color="#888888"><br><br></font></span><div class="gmail_quote">
<div>On Mon, Jul 23, 2012 at 10:12 AM, Hong Zhang <span dir="ltr"><<a href="mailto:hzhang@mcs.anl.gov" target="_blank">hzhang@mcs.anl.gov</a>></span> wrote:<br></div><div><div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div>Xuan will take a look at this.</div><span><font color="#888888"><div>Hong </div></font></span><div><div><br><div class="gmail_quote">On Mon, Jul 23, 2012 at 9:47 AM, Matthew Knepley <span dir="ltr"><<a href="mailto:knepley@gmail.com" target="_blank">knepley@gmail.com</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div>I think we just need a small converter from AIJ:</div><div><br></div><a href="https://bitbucket.org/poulson/clique/src/dc417c7e9403/tests/DistSparseMatrix.cpp" target="_blank">https://bitbucket.org/poulson/clique/src/dc417c7e9403/tests/DistSparseMatrix.cpp</a><div>
<br></div><div> Matt<span><font color="#888888"><br clear="all"><div><br></div>-- <br>What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.<br>
-- Norbert Wiener<br>
</font></span></div>
</blockquote></div><br>
</div></div></blockquote></div></div></div><br>
</blockquote></div></div><br>
</blockquote></div><br>