<div dir="ltr"><br><br><div class="gmail_quote"><div dir="ltr">On Thu, Sep 20, 2018 at 12:09 PM Matthew Knepley <<a href="mailto:knepley@gmail.com">knepley@gmail.com</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div class="gmail_quote"><div dir="ltr">On Thu, Sep 20, 2018 at 1:26 PM Fande Kong <<a href="mailto:fdkong.jd@gmail.com" target="_blank">fdkong.jd@gmail.com</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Hi Developers,<div><br></div><div>MATBAIJ actually assumes that the point-block is dense. It is fine if the block size is small for example less than 50, but we may suffer from a memory issue if the block is large such as 1000. The block is not necessary dense. In my application, neutron transport equations, the block is actually sparse (coming from direction and energy grids).  Now I am using MATAIJ, but I want to take into consideration there are many variables on each spatial grid when doing some ILU(0).  ILU(0) is used as a subdomain solver for ASM and BJACOBI. </div><div><br></div><div>I would like to support point-block ILU(0) in MATAIJ (not MATBAIJ) because the block is sparse. Any comments? It is worthwhile to do that? I think we might reduce the symbolic factor time when considering the block.</div></div></blockquote><div><br></div><div>What optimization do you hope to take advantage of? </div></div></div></blockquote><div><br></div><div>The inverse of sparse block can be approximated with the inverse of the block diagonal. </div><div><br></div><div>Does it help reduce the compute time in symbolic factorization? </div><div><br></div><div>Fande,</div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div class="gmail_quote"><div>The main optimization in sparse LU is finding dense blocks,</div><div>and it should do a good job of that. I do not see what defining huge, sparse blocks does for you. </div></div></div></blockquote><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div class="gmail_quote"><div><br></div><div>  Thanks,</div><div><br></div><div>     Matt</div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div>Thanks,</div><div><br></div><div>Fande,</div></div>
</blockquote></div><br clear="all"><div><br></div>-- <br><div dir="ltr" class="m_-6985450323502940962gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div>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</div><div><br></div><div><a href="http://www.cse.buffalo.edu/~knepley/" target="_blank">https://www.cse.buffalo.edu/~knepley/</a><br></div></div></div></div></div></div></div></div>
</blockquote></div></div>