<div dir="ltr">Hi,<div><br></div><div>I am trying to determine where there is malloc when I am doing a MatSetValuesBlocked.  The following is the error:</div><div><br></div><div><div>[1]PETSC ERROR: --------------------- Error Message ------------------------------------</div>
<div>[1]PETSC ERROR: Argument out of range!</div><div>[1]PETSC ERROR: New nonzero at (69216,95036) caused a malloc!</div><div><br></div><div style>However, when I look at the info during the assembly, I have </div><div style>
<br></div><div style><div>[0] MatAssemblyBegin_MPIBAIJ(): Stash has 0 entries,uses 0 mallocs.</div><div>[0] MatAssemblyBegin_MPIBAIJ(): Block-Stash has 0 entries, uses 0 mallocs.</div><div>[1] MatAssemblyEnd_SeqBAIJ(): Matrix size: 381300 X 381300, block size 4; storage space: 14143776 unneeded, 10259424 used</div>
<div>[1] MatAssemblyEnd_SeqBAIJ(): Number of mallocs during MatSetValues is 0</div><div>[1] MatAssemblyEnd_SeqBAIJ(): Most nonzeros blocks in any row is 16</div><div>[0] MatAssemblyEnd_SeqBAIJ(): Matrix size: 381300 X 381300, block size 4; storage space: 13696176 unneeded, 10707024 used</div>
<div>[0] MatAssemblyEnd_SeqBAIJ(): Number of mallocs during MatSetValues is 0</div><div>[0] MatAssemblyEnd_SeqBAIJ(): Most nonzeros blocks in any row is 11</div><div><br></div><div style>which indicates there is mallocs.</div>
<div style><br></div><div style>The following is how I set up the matrix:</div><div style><br></div><div style><div>call MatCreate(PETSC_COMM_WORLD,petsc_obj%jacmat, pierr) </div><div>call MatSetSizes(petsc_obj%jacmat,PETSC_DECIDE,PETSC_DECIDE,nel*neq,nel*neq,pierr)</div>
<div>call MatSetFromOptions(petsc_obj%jacmat,pierr)<br></div><div>call MatSetBlockSize(petsc_obj%jacmat,neq,pierr)<br></div><div>call MatSetType(petsc_obj%jacmat,MATBAIJ,pierr)</div><div>call MatSeqBAIJSetPreallocation(petsc_obj%jacmat,neq,(mncon+1)*neq*neq,PETSC_NULL_INTEGER,pierr)</div>
<div>call MatMPIBAIJSetPreallocation(petsc_obj%jacmat,   </div><div>neq,neq*neq,PETSC_NULL_INTEGER,mncon*neq*neq,PETSC_NULL_INTEGER,pierr)</div><div>    call MatSetOption(petsc_obj%jacmat,MAT_ROW_ORIENTED,PETSC_FALSE,pierr) !because of how aval is stored.</div>
</div></div><div><br></div>-- <br><span style="color:rgb(34,34,34);font-family:arial,sans-serif;font-size:13px;background-color:rgb(255,255,255)">George Pau</span><br style="color:rgb(34,34,34);font-family:arial,sans-serif;font-size:13px;background-color:rgb(255,255,255)">
<span style="color:rgb(34,34,34);font-family:arial,sans-serif;font-size:13px;background-color:rgb(255,255,255)">Earth Sciences Division</span><br style="color:rgb(34,34,34);font-family:arial,sans-serif;font-size:13px;background-color:rgb(255,255,255)">
<span style="color:rgb(34,34,34);font-family:arial,sans-serif;font-size:13px;background-color:rgb(255,255,255)">Lawrence Berkeley National Laboratory</span><br style="color:rgb(34,34,34);font-family:arial,sans-serif;font-size:13px;background-color:rgb(255,255,255)">
<span style="color:rgb(34,34,34);font-family:arial,sans-serif;font-size:13px;background-color:rgb(255,255,255)">One Cyclotron, MS 74-120</span><br style="color:rgb(34,34,34);font-family:arial,sans-serif;font-size:13px;background-color:rgb(255,255,255)">
<span style="color:rgb(34,34,34);font-family:arial,sans-serif;font-size:13px;background-color:rgb(255,255,255)">Berkeley, CA 94720</span><br style="color:rgb(34,34,34);font-family:arial,sans-serif;font-size:13px;background-color:rgb(255,255,255)">
<br style="color:rgb(34,34,34);font-family:arial,sans-serif;font-size:13px;background-color:rgb(255,255,255)">(510) 486-7196<br style="color:rgb(34,34,34);font-family:arial,sans-serif;font-size:13px;background-color:rgb(255,255,255)">
<a href="mailto:gpau@lbl.gov" style="color:rgb(17,85,204);font-family:arial,sans-serif;font-size:13px;background-color:rgb(255,255,255)" target="_blank">gpau@lbl.gov</a><br style="color:rgb(34,34,34);font-family:arial,sans-serif;font-size:13px;background-color:rgb(255,255,255)">
<a href="http://esd.lbl.gov/about/staff/georgepau/" style="color:rgb(17,85,204);font-family:arial,sans-serif;font-size:13px;background-color:rgb(255,255,255)" target="_blank">http://esd.lbl.gov/about/staff/georgepau/</a><br>

</div></div>