<div dir="ltr"><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">Il giorno mar 1 feb 2022 alle ore 18:34 Jed Brown <<a href="mailto:jed@jedbrown.org">jed@jedbrown.org</a>> ha scritto:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">Patrick Sanan <<a href="mailto:patrick.sanan@gmail.com" target="_blank">patrick.sanan@gmail.com</a>> writes:<br>
<br>
> Am Di., 1. Feb. 2022 um 16:20 Uhr schrieb Jed Brown <<a href="mailto:jed@jedbrown.org" target="_blank">jed@jedbrown.org</a>>:<br>
><br>
>> Patrick Sanan <<a href="mailto:patrick.sanan@gmail.com" target="_blank">patrick.sanan@gmail.com</a>> writes:<br>
>><br>
>> > Sorry about the delay on this. I can reproduce.<br>
>> ><br>
>> > This regression appears to be a result of this optimization:<br>
>> > <a href="https://gitlab.com/petsc/petsc/-/merge_requests/4273" rel="noreferrer" target="_blank">https://gitlab.com/petsc/petsc/-/merge_requests/4273</a><br>
>><br>
>> Thanks for tracking this down. Is there a reason to prefer preallocating<br>
>> twice<br>
>><br>
>>    ierr =<br>
>> MatPreallocatorPreallocate(preallocator,PETSC_TRUE,A);CHKERRQ(ierr);<br>
>>    ierr =<br>
>> MatPreallocatorPreallocate(preallocator,PETSC_TRUE,A_duplicate);CHKERRQ(ierr);<br>
>><br>
>> versus using MatDuplicate() or MatConvert()?<br>
>><br></blockquote><div><br></div><div>Jed</div><div><br></div><div>this is not the point. Suppose you pass around only a preallocator, but do not pass around the matrices. Reusing the preallocator should be allowed. </div><div><br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
><br>
> Maybe if your preallocation is an overestimate for each of two different<br>
> post-assembly non-zero structures in A and A_duplicate?<br>
<br>
Even then, why not preallocate A and duplicate immediately, before compressing out zeros?<br>
</blockquote></div><br clear="all"><div><br></div>-- <br><div dir="ltr" class="gmail_signature">Stefano</div></div>