<html><head></head><body><div style="color:#000; background-color:#fff; font-family:Helvetica Neue, Helvetica, Arial, Lucida Grande, sans-serif;font-size:13px">Thanks.<div id="yui_3_16_0_ym19_1_1493214000562_20941"><span></span></div> <div class="qtdSeparateBR"><br><br></div><div class="yahoo_quoted" style="display: block;"> <div style="font-family: Helvetica Neue, Helvetica, Arial, Lucida Grande, sans-serif; font-size: 13px;"> <div style="font-family: HelveticaNeue, Helvetica Neue, Helvetica, Arial, Lucida Grande, sans-serif; font-size: 16px;"> <div dir="ltr"><font size="2" face="Arial"> On Wednesday, April 26, 2017 1:28 PM, Matthew Knepley <knepley@gmail.com> wrote:<br></font></div> <br><br> <div class="y_msg_container"><div id="yiv5668315034"><div><div dir="ltr"><div class="yiv5668315034gmail_extra"><div class="yiv5668315034gmail_quote">On Wed, Apr 26, 2017 at 12:10 PM, D D <span dir="ltr"><<a rel="nofollow" shape="rect" ymailto="mailto:driver.dan12@yahoo.com" target="_blank" href="mailto:driver.dan12@yahoo.com">driver.dan12@yahoo.com</a>></span> wrote:<br clear="none"><blockquote class="yiv5668315034gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;"><div><div style="color:#000;background-color:#fff;font-family:Helvetica Neue, Helvetica, Arial, Lucida Grande, sans-serif;font-size:13px;"><div id="yiv5668315034m_8172148363751200006yui_3_16_0_ym19_1_1493214000562_18520">Why does MatAlloc_AIJ_CSR make two preallocation calls?<br clear="none"></div><div id="yiv5668315034m_8172148363751200006yui_3_16_0_ym19_1_1493214000562_18735"><br clear="none"></div><div dir="ltr" id="yiv5668315034m_8172148363751200006yui_3_16_0_ym19_1_1493214000562_18719"> CHKERR( MatSeqAIJSetPreallocationCSR( A, i, j, v) )<br id="yiv5668315034m_8172148363751200006yui_3_16_0_ym19_1_1493214000562_18724" clear="none"> CHKERR( MatMPIAIJSetPreallocationCSR( A, i, j, v) )<br id="yiv5668315034m_8172148363751200006yui_3_16_0_ym19_1_1493214000562_18725" clear="none"></div><div dir="ltr" id="yiv5668315034m_8172148363751200006yui_3_16_0_ym19_1_1493214000562_18726"><br clear="none"></div><div dir="ltr" id="yiv5668315034m_8172148363751200006yui_3_16_0_ym19_1_1493214000562_18754">I am assuming memory if preallocated twice, once for sequential and another for MPI. So each matrix created with createAIJ will have a sequential and MPI structure.</div><div dir="ltr" id="yiv5668315034m_8172148363751200006yui_3_16_0_ym19_1_1493214000562_18769"><br clear="none"></div><div dir="ltr" id="yiv5668315034m_8172148363751200006yui_3_16_0_ym19_1_1493214000562_18783">Is this for convenience when switching from sequential to MPI?</div></div></div></blockquote><div><br clear="none"></div><div>Only one function will actually take affect, depending on the type of matrix A is. This is like Objective-C rather than C++.</div><div><br clear="none"></div><div> Thanks,</div><div><br clear="none"></div><div> Matt</div><div class="yiv5668315034yqt2715218642" id="yiv5668315034yqtfd83713"><div> </div></div><blockquote class="yiv5668315034gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;"><div><div style="color:#000;background-color:#fff;font-family:Helvetica Neue, Helvetica, Arial, Lucida Grande, sans-serif;font-size:13px;"><div class="yiv5668315034yqt2715218642" id="yiv5668315034yqtfd35514"><div dir="ltr" id="yiv5668315034m_8172148363751200006yui_3_16_0_ym19_1_1493214000562_18787">Thanks,</div></div><div dir="ltr"><div class="yiv5668315034yqt2715218642" id="yiv5668315034yqtfd27587">Dale</div><br clear="none"></div></div></div></blockquote></div><br clear="none"><br clear="all"><div><br clear="none"></div>-- <br clear="none"><div class="yiv5668315034gmail_signature">What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.<br clear="none">-- Norbert Wiener</div><div class="yiv5668315034yqt2715218642" id="yiv5668315034yqtfd79317">
</div></div></div></div></div><br><br></div> </div> </div> </div></div></body></html>