<div class="gmail_quote">2012/3/16 Xiangze Zeng <span dir="ltr"><<a href="mailto:zengshixiangze@163.com">zengshixiangze@163.com</a>></span><br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div style="line-height:1.7;color:#000000;font-size:14px;font-family:arial">I replace MATAIJCUSP with MATMPIAIJCUSP. It works well. The attachments are the outputs of log_summary.<div>MATAIJCUSP is used for serial matrices?</div>
</div></blockquote></div><br><div>MATAIJCUSP automatically chooses MATSEQAIJCUSP if on serial communicator and MATMPIAIJCUSP if on a parallel communicator. Practically speaking, either will work, but when you run in serial, there is more overhead going through the MPIAIJCUSP container and fewer preconditioners will be available. So it's better to call bath MatSeqAIJSetPreallocation() and MatMPIAIJSetPreallocation().</div>