<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=utf-8">
</head>
<body style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;" class="">
Hi Dave,
<div class=""><br class="">
</div>
<div class="">Yes it is parallel so the preallocation calls are not lowered by the allocation.</div>
<div class=""><br class="">
</div>
<div class="">I am trying to use MatXAIJSetPreallocation, but not sure how, since the following link does not give an example:</div>
<div class=""><br class="">
</div>
<div class=""><a href="https://www.mcs.anl.gov/petsc/petsc-current/docs/manualpages/Mat/MatXAIJSetPreallocation.html" class="">https://www.mcs.anl.gov/petsc/petsc-current/docs/manualpages/Mat/MatXAIJSetPreallocation.html</a></div>
<div class=""><br class="">
</div>
<div class="">If I have the following matrix:</div>
<div class=""><br class="">
</div>
<div class="">0 1 2 0 </div>
<div class="">1 0 0 0 </div>
<div class="">2 0 1 3</div>
<div class="">0 0 3 2 </div>
<div class=""><br class="">
</div>
<div class="">How should I put in the parameters of MatXAIJSetPreallocation?</div>
<div class=""><br class="">
</div>
<div class="">Thanks!</div>
<div class=""><br class="">
</div>
<div class="">Cheers,</div>
<div class=""><br class="">
</div>
<div class="">Yang Bo</div>
<div class=""><br class="">
<div><br class="">
<blockquote type="cite" class="">
<div class="">On 21 May 2020, at 5:42 PM, Dave May <<a href="mailto:dave.mayhem23@gmail.com" class="">dave.mayhem23@gmail.com</a>> wrote:</div>
<br class="Apple-interchange-newline">
<div class=""><br class="Apple-interchange-newline">
<br style="caret-color: rgb(0, 0, 0); font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none;" class="">
<div class="gmail_quote" style="caret-color: rgb(0, 0, 0); font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none;">
<div dir="ltr" class="">On Thu 21. May 2020 at 10:49, Yang Bo (Asst Prof) <<a href="mailto:yang.bo@ntu.edu.sg" class="">yang.bo@ntu.edu.sg</a>> wrote:<br class="">
</div>
<blockquote class="gmail_quote" style="margin: 0px 0px 0px 0.8ex; border-left-width: 1px; border-left-style: solid; border-left-color: rgb(204, 204, 204); padding-left: 1ex;">
<div style="word-wrap: break-word; line-break: after-white-space;" class="">Hi Dave,
<div class=""><br class="">
</div>
<div class="">Thank you very much for your reply. That is indeed the problem. I have been working with matrices in Slepc but I don’t really understand it. I tried to preallocate but it still does not work.<span class="Apple-converted-space"> </span></div>
</div>
</blockquote>
<div dir="auto" class=""><br class="">
</div>
<div dir="auto" class="">Meaning the number of reported mallocs is still non-zero?</div>
<div dir="auto" class="">Is the number reported with you preallocation calls lower than what you originally saw?</div>
<div dir="auto" class=""><br class="">
</div>
<blockquote class="gmail_quote" style="margin: 0px 0px 0px 0.8ex; border-left-width: 1px; border-left-style: solid; border-left-color: rgb(204, 204, 204); padding-left: 1ex;">
<div style="word-wrap: break-word; line-break: after-white-space;" class="">
<div class="">If you look at my code below:</div>
<div class=""><br class="">
</div>
<div class="">
<div class="">ierr = MatCreate(PETSC_COMM_WORLD,&A);CHKERRQ(ierr);</div>
<div class="">ierr = MatSetSizes(A,PETSC_DECIDE,PETSC_DECIDE,h_dim,h_dim);<span style="white-space: pre-wrap;" class="">
</span>                     // h_dim is the dimension of the square matrix A</div>
<div class="">ierr = MatSetFromOptions(A);CHKERRQ(ierr);</div>
<div class="">ierr = MatSetUp(A);CHKERRQ(ierr);</div>
<div class="">ierr = MatGetOwnershipRange(A,&Istart,&Iend);CHKERRQ(ierr);</div>
<div class=""><br class="">
</div>
<div class="">MatSeqAIJSetPreallocation(A,0,nnz);                                      <span class="Apple-converted-space"> </span><span style="white-space: pre-wrap;" class=""></span>// I try to preallocate here, where nnz is the array containing the number
 of non-zero entries each row</div>
<div class=""><br class="">
</div>
<div class="">for (int i=0;i<row.size();i++) {</div>
<div class=""><span style="white-space: pre-wrap;" class=""></span>MatSetValue(A,row[i],column[i],h[i],INSERT_VALUES);</div>
<div class="">}</div>
</div>
<div class=""><br class="">
</div>
<div class="">I am not sure what other information I need to give for the pre-allocation…</div>
<div class=""></div>
</div>
</blockquote>
<div dir="auto" class=""><br class="">
</div>
<div dir="auto" class="">This looks fine. However MatSeqAIJSetPreallocation() has no effect if the Mat type is not SEQAIJ.</div>
<div dir="auto" class=""><br class="">
</div>
<div dir="auto" class="">Are you running in parallel? If yes then the Mat type will be MATMPIAIJ and you either have to call the MPI specific preallocator or use the generic one I pointed you too.</div>
<div dir="auto" class=""><br class="">
</div>
<div dir="auto" class="">Thanks</div>
<div dir="auto" class="">Dave</div>
<div dir="auto" class=""><br class="">
</div>
<div dir="auto" class=""><br class="">
</div>
<blockquote class="gmail_quote" style="margin: 0px 0px 0px 0.8ex; border-left-width: 1px; border-left-style: solid; border-left-color: rgb(204, 204, 204); padding-left: 1ex;">
<div style="word-wrap: break-word; line-break: after-white-space;" class="">
<div class=""><br class="">
</div>
<div class="">Cheers,</div>
<div class=""><br class="">
</div>
<div class="">Yang Bo</div>
<div class=""><br class="">
</div>
<div class=""><br class="">
<div class=""><br class="">
<blockquote type="cite" class="">
<div class="">On 21 May 2020, at 4:08 PM, Dave May <<a href="mailto:dave.mayhem23@gmail.com" target="_blank" class="">dave.mayhem23@gmail.com</a>> wrote:</div>
<br class="">
<div class=""><b style="font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-caps: normal; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; text-decoration: none;" class="">-info
 | grep malloc</b></div>
</blockquote>
</div>
<br class="">
</div>
</div>
<div style="word-wrap: break-word; line-break: after-white-space;" class="">
<hr class="">
<p style="font-size: 10pt; color: rgb(128, 128, 128); font-family: Arial;" class="">
CONFIDENTIALITY: This email is intended solely for the person(s) named and may be confidential and/or privileged. If you are not the intended recipient, please delete it, notify us and do not copy, use, or disclose its contents.<span class="Apple-converted-space"> </span><br class="">
Towards a sustainable earth: Print only when necessary. Thank you.<span class="Apple-converted-space"> </span></p>
</div>
</blockquote>
</div>
</div>
</blockquote>
</div>
<br class="">
</div>
</body>
</html>