<div dir="ltr">Thanks, Hong,<div><br></div><div>I will try the code with bs=1, and report back to you.</div><div><br></div><div>Fande,</div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Tue, Mar 31, 2020 at 9:51 PM Zhang, Hong <<a href="mailto:hzhang@mcs.anl.gov">hzhang@mcs.anl.gov</a>> wrote:<br></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 dir="ltr">
<div style="font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
Fande,</div>
<div style="font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
Checking aij.result:</div>
<div style="font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
<span>Mat Object: () 1 MPI processes<br>
</span>
<div>          type: seqaij<br>
</div>
<div>          rows=25816, cols=25816, bs=4<br>
</div>
<div>          total: nonzeros=1297664, allocated nonzeros=1297664<br>
</div>
<div>          total number of mallocs used during MatSetValues calls=0<br>
</div>
<span>            using I-node routines: found 6454 nodes, limit used is 5</span><br>
</div>
<div style="font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
<span><br>
</span></div>
<div style="font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
<span>i.e., it uses bs=4 with I-node. The implementation of MatSolve() is similar to baij with bs=4. What happens if you try aij with '-matload_block_size 1 -mat_no_inode true'?</span></div>
<div style="font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
<span>Hong</span></div>
<div id="gmail-m_-8852262757798538868appendonsend"></div>
<div style="font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
<br>
</div>
<hr style="display:inline-block;width:98%">
<div id="gmail-m_-8852262757798538868divRplyFwdMsg" dir="ltr"><font face="Calibri, sans-serif" color="#000000" style="font-size:11pt"><b>From:</b> petsc-users <<a href="mailto:petsc-users-bounces@mcs.anl.gov" target="_blank">petsc-users-bounces@mcs.anl.gov</a>> on behalf of Fande Kong <<a href="mailto:fdkong.jd@gmail.com" target="_blank">fdkong.jd@gmail.com</a>><br>
<b>Sent:</b> Monday, March 30, 2020 12:25 PM<br>
<b>To:</b> PETSc users list <<a href="mailto:petsc-users@mcs.anl.gov" target="_blank">petsc-users@mcs.anl.gov</a>><br>
<b>Subject:</b> [petsc-users] AIJ vs BAIJ when using ILU factorization</font>
<div> </div>
</div>
<div>
<div dir="ltr">Hi All,<br>
<div><br>
</div>
<div>There is a system of equations arising from the discretization of 3D incompressible Navier-Stoke equations using a finite element method. 4 unknowns are placed on each mesh point, and then there is a 4x4 saddle point block on each mesh vertex.  I was thinking
 to solve the linear equations using an incomplete LU factorization (that will be eventually used as a subdomain solver for ASM).  </div>
<div><br>
</div>
<div>Right now, I am trying to study the ILU performance using AIJ and BAIJ, respectively. From my understanding, BAIJ should give me better results since it inverses the 4x4 blocks exactly, while AIJ does not. However, I found that both BAIJ and AIJ gave me
 identical results in terms of the number of iterations.  Was that just coincident?  Or in theory, they are just identical.  I understand the runtimes may be different because BAIJ has a better data locality.</div>
<div><br>
</div>
<div><br>
</div>
<div>Please see the attached files for the results and solver configuration.</div>
<div><br>
</div>
<div><br>
</div>
<div>Thanks,</div>
<div><br>
</div>
<div>Fande,</div>
</div>
</div>
</div>

</blockquote></div>