<head><!-- BaNnErBlUrFlE-HeAdEr-start -->
<style>
#pfptBannerhiz612y { all: revert !important; display: block !important;
visibility: visible !important; opacity: 1 !important;
background-color: #D0D8DC !important;
max-width: none !important; max-height: none !important }
.pfptPrimaryButtonhiz612y:hover, .pfptPrimaryButtonhiz612y:focus {
background-color: #b4c1c7 !important; }
.pfptPrimaryButtonhiz612y:active {
background-color: #90a4ae !important; }
</style>
<!-- BaNnErBlUrFlE-HeAdEr-end -->
</head><!-- BaNnErBlUrFlE-BoDy-start -->
<!-- Preheader Text : BEGIN -->
<div style="display:none !important;display:none;visibility:hidden;mso-hide:all;font-size:1px;color:#ffffff;line-height:1px;height:0px;max-height:0px;opacity:0;overflow:hidden;">
Hi, Yongzhong, Since the two kernels of KSPGMRESOrthog are VecMDot and VecMAXPY, if we can speed up the two with OpenMP threads, then we can speed up KSPGMRESOrthog. We recently added an optimization to do VecMDot/MAXPY() in dense matrix-vector
</div>
<!-- Preheader Text : END -->
<!-- Email Banner : BEGIN -->
<div style="display:none !important;display:none;visibility:hidden;mso-hide:all;font-size:1px;color:#ffffff;line-height:1px;height:0px;max-height:0px;opacity:0;overflow:hidden;">ZjQcmQRYFpfptBannerStart</div>
<!--[if ((ie)|(mso))]>
<table border="0" cellspacing="0" cellpadding="0" width="100%" style="padding: 16px 0px 16px 0px; direction: ltr" ><tr><td>
<table border="0" cellspacing="0" cellpadding="0" style="padding: 0px 10px 5px 6px; width: 100%; border-radius:4px; border-top:4px solid #90a4ae;background-color:#D0D8DC;"><tr><td valign="top">
<table align="left" border="0" cellspacing="0" cellpadding="0" style="padding: 4px 8px 4px 8px">
<tr><td style="color:#000000; font-family: 'Arial', sans-serif; font-weight:bold; font-size:14px; direction: ltr">
This Message Is From an External Sender
</td></tr>
<tr><td style="color:#000000; font-weight:normal; font-family: 'Arial', sans-serif; font-size:12px; direction: ltr">
This message came from outside your organization.
</td></tr>
</table>
</td></tr></table>
</td></tr></table>
<![endif]-->
<![if !((ie)|(mso))]>
<div dir="ltr" id="pfptBannerhiz612y" style="all: revert !important; display:block !important; text-align: left !important; margin:16px 0px 16px 0px !important; padding:8px 16px 8px 16px !important; border-radius: 4px !important; min-width: 200px !important; background-color: #D0D8DC !important; background-color: #D0D8DC; border-top: 4px solid #90a4ae !important; border-top: 4px solid #90a4ae;">
<div id="pfptBannerhiz612y" style="all: unset !important; float:left !important; display:block !important; margin: 0px 0px 1px 0px !important; max-width: 600px !important;">
<div id="pfptBannerhiz612y" style="all: unset !important; display:block !important; visibility: visible !important; background-color: #D0D8DC !important; color:#000000 !important; color:#000000; font-family: 'Arial', sans-serif !important; font-family: 'Arial', sans-serif; font-weight:bold !important; font-weight:bold; font-size:14px !important; line-height:18px !important; line-height:18px">
This Message Is From an External Sender
</div>
<div id="pfptBannerhiz612y" style="all: unset !important; display:block !important; visibility: visible !important; background-color: #D0D8DC !important; color:#000000 !important; color:#000000; font-weight:normal; font-family: 'Arial', sans-serif !important; font-family: 'Arial', sans-serif; font-size:12px !important; line-height:18px !important; line-height:18px; margin-top:2px !important;">
This message came from outside your organization.
</div>
</div>
<div style="clear: both !important; display: block !important; visibility: hidden !important; line-height: 0 !important; font-size: 0.01px !important; height: 0px"> </div>
</div>
<![endif]>
<div style="display:none !important;display:none;visibility:hidden;mso-hide:all;font-size:1px;color:#ffffff;line-height:1px;height:0px;max-height:0px;opacity:0;overflow:hidden;">ZjQcmQRYFpfptBannerEnd</div>
<!-- Email Banner : END -->
<!-- BaNnErBlUrFlE-BoDy-end -->
<div dir="ltr"><div>Hi, Yongzhong,</div><div> Since the two kernels of KSPGMRESOrthog are VecMDot and VecMAXPY, if we can speed up the two with OpenMP threads, then we can speed up KSPGMRESOrthog. We recently added an optimization to do VecMDot/MAXPY() in dense matrix-vector multiplication (i.e., BLAS2 GEMV, with tall-and-skinny matrices ). So with MKL_VERBOSE=1, you should see something like "MKL_VERBOSE ZGEMV ..." in output. If not, could you try again with petsc/main?</div><div> petsc has a microbenchmark (vec/vec/tests/ex2k.c) to test them. I ran VecMDot with multithreaded oneMKL (via setting MKL_NUM_THREADS), it was strange to see no speedup. I then configured petsc with openblas, I did see better performance with more threads</div><div><br></div><div><font face="monospace">$ OMP_PROC_BIND=spread OMP_NUM_THREADS=1 ./ex2k -n 15 -m 4<br>Vector(N) VecMDot-3 VecMDot-8 VecMDot-30 (us)<br>--------------------------------------------------------------------------<br> 128 2.0 2.5 6.1 <br> 256 1.8 2.7 7.0 <br> 512 2.1 3.1 8.6 <br> 1024 2.7 4.0 12.3 <br> 2048 3.8 6.3 28.0 <br> 4096 6.1 10.6 42.4 <br> 8192 10.9 21.8 79.5 <br> 16384 21.2 39.4 149.6 <br> 32768 45.9 75.7 224.6 <br> 65536 142.2 215.8 732.1 <br> 131072 169.1 233.2 1729.4 <br> 262144 367.5 830.0 4159.2 <br> 524288 999.2 1718.1 8538.5 <br> 1048576 2113.5 4082.1 18274.8 <br> 2097152 5392.6 10273.4 43273.4 </font><br></div><div><br></div><div><br></div><div><font face="monospace">$ OMP_PROC_BIND=spread OMP_NUM_THREADS=8 ./ex2k -n 15 -m 4<br>Vector(N) VecMDot-3 VecMDot-8 VecMDot-30 (us)<br>--------------------------------------------------------------------------<br> 128 2.0 2.5 6.0 <br> 256 1.8 2.7 15.0 <br> 512 2.1 9.0 16.6 <br> 1024 2.6 8.7 16.1 <br> 2048 7.7 10.3 20.5 <br> 4096 9.9 11.4 25.9 <br> 8192 14.5 22.1 39.6 <br> 16384 25.1 27.8 67.8 <br> 32768 44.7 95.7 91.5 <br> 65536 82.1 156.8 165.1 <br> 131072 194.0 335.1 341.5 <br> 262144 388.5 380.8 612.9 <br> 524288 1046.7 967.1 1653.3 <br> 1048576 1997.4 2169.0 4034.4 <br> 2097152 5502.9 5787.3 12608.1</font></div><div> <br></div><div>The tall-and-skinny matrices in KSPGMRESOrthog vary in width. The average speedup depends on components. So I suggest you run ex2k to see in your environment whether oneMKL can speedup the kernels. <br></div><div><br></div><div><div dir="ltr" class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr">--Junchao Zhang</div></div></div><br></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Mon, Jun 24, 2024 at 11:35 AM Junchao Zhang <<a href="mailto:junchao.zhang@gmail.com">junchao.zhang@gmail.com</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr">Let me run some examples on our end to see whether the code calls expected functions. <div><br clear="all"><div><div dir="ltr" class="gmail_signature"><div dir="ltr">--Junchao Zhang</div></div></div><br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Mon, Jun 24, 2024 at 10:46 AM Matthew Knepley <<a href="mailto:knepley@gmail.com" target="_blank">knepley@gmail.com</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div>
<div style="font-size:1px;color:rgb(255,255,255);line-height:1px;height:0px;max-height:0px;opacity:0;overflow:hidden;display:none">
On Mon, Jun 24, 2024 at 11: 21 AM Yongzhong Li <yongzhong. li@ mail. utoronto. ca> wrote: Thank you Pierre for your information. Do we have a conclusion for my original question about the parallelization efficiency for different stages of
</div>
<div style="font-size:1px;color:rgb(255,255,255);line-height:1px;height:0px;max-height:0px;opacity:0;overflow:hidden;display:none">ZjQcmQRYFpfptBannerStart</div>
<u></u>
<div dir="ltr" id="m_-7305460844614688919m_2643581689528628203pfptBanner0gnzkk2" style="display:block;text-align:left;margin:16px 0px;padding:8px 16px;border-radius:4px;min-width:200px;background-color:rgb(208,216,220);border-top:4px solid rgb(144,164,174)">
<div id="m_-7305460844614688919m_2643581689528628203pfptBanner0gnzkk2" style="float:left;display:block;margin:0px 0px 1px;max-width:600px">
<div id="m_-7305460844614688919m_2643581689528628203pfptBanner0gnzkk2" style="display:block;background-color:rgb(208,216,220);color:rgb(0,0,0);font-family:Arial,sans-serif;font-weight:bold;font-size:14px;line-height:18px">
This Message Is From an External Sender
</div>
<div id="m_-7305460844614688919m_2643581689528628203pfptBanner0gnzkk2" style="font-weight:normal;display:block;background-color:rgb(208,216,220);color:rgb(0,0,0);font-family:Arial,sans-serif;font-size:12px;line-height:18px;margin-top:2px">
This message came from outside your organization.
</div>
</div>
<div style="height:0px;clear:both;display:block;line-height:0;font-size:0.01px"> </div>
</div>
<u></u>
<div style="font-size:1px;color:rgb(255,255,255);line-height:1px;height:0px;max-height:0px;opacity:0;overflow:hidden;display:none">ZjQcmQRYFpfptBannerEnd</div>
<div dir="ltr"><div dir="ltr">On Mon, Jun 24, 2024 at 11:21 AM Yongzhong Li <<a href="mailto:yongzhong.li@mail.utoronto.ca" target="_blank">yongzhong.li@mail.utoronto.ca</a>> wrote:<br></div><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div>
<div style="font-size:1px;color:rgb(255,255,255);line-height:1px;height:0px;max-height:0px;opacity:0;overflow:hidden;display:none">
Thank you Pierre for your information. Do we have a conclusion for my original question about the parallelization efficiency for different stages of KSP Solve? Do we need to do more testing to figure out the issues? Thank you, Yongzhong From:
</div>
<div style="font-size:1px;color:rgb(255,255,255);line-height:1px;height:0px;max-height:0px;opacity:0;overflow:hidden;display:none">ZjQcmQRYFpfptBannerStart</div>
<u></u>
<div dir="ltr" id="m_-7305460844614688919m_2643581689528628203m_6159968039748304261pfptBanner7d6k4n9" style="display:block;text-align:left;margin:16px 0px;padding:8px 16px;border-radius:4px;min-width:200px;background-color:rgb(208,216,220);border-top:4px solid rgb(144,164,174)">
<div id="m_-7305460844614688919m_2643581689528628203m_6159968039748304261pfptBanner7d6k4n9" style="float:left;display:block;margin:0px 0px 1px;max-width:600px">
<div id="m_-7305460844614688919m_2643581689528628203m_6159968039748304261pfptBanner7d6k4n9" style="display:block;background-color:rgb(208,216,220);color:rgb(0,0,0);font-family:Arial,sans-serif;font-weight:bold;font-size:14px;line-height:18px">
This Message Is From an External Sender
</div>
<div id="m_-7305460844614688919m_2643581689528628203m_6159968039748304261pfptBanner7d6k4n9" style="font-weight:normal;display:block;background-color:rgb(208,216,220);color:rgb(0,0,0);font-family:Arial,sans-serif;font-size:12px;line-height:18px;margin-top:2px">
This message came from outside your organization.
</div>
</div>
<div style="height:0px;clear:both;display:block;line-height:0;font-size:0.01px"> </div>
</div>
<u></u>
<div style="font-size:1px;color:rgb(255,255,255);line-height:1px;height:0px;max-height:0px;opacity:0;overflow:hidden;display:none">ZjQcmQRYFpfptBannerEnd</div>
<div lang="en-CN">
<div>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11pt">Thank you Pierre for your information. Do we have a conclusion for my original question about the parallelization efficiency for different stages of KSP Solve? Do we need to do more testing to
figure out the issues?<br></span></p></div></div></div></blockquote><div><br></div><div>We have an extended discussion of this here: <a href="https://urldefense.us/v3/__https://petsc.org/release/faq/*what-kind-of-parallel-computers-or-clusters-are-needed-to-use-petsc-or-why-do-i-get-little-speedup__;Iw!!G_uCfscf7eWS!aQJpmm5W6l6FUiumnIPmkouzwzNUfx-Dyq04i1O2KS_InQGk6qjI7wUir0Hx6QEUQE2AMiJDsez3x4zRO7V_$" target="_blank">https://petsc.org/release/faq/#what-kind-of-parallel-computers-or-clusters-are-needed-to-use-petsc-or-why-do-i-get-little-speedup</a></div><div><br></div><div>The kinds of operations you are talking about (SpMV, VecDot, VecAXPY, etc) are memory bandwidth limited. If there is no more bandwidth to be marshalled on your board, then adding more processes does nothing at all. This is why people were asking about how many "nodes" you are running on, because that is the unit of memory bandwidth, not "cores" which make little difference.</div><div><br></div><div> Thanks,</div><div><br></div><div> Matt</div><div> </div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div><div lang="en-CN"><div><p class="MsoNormal"><span lang="EN-US" style="font-size:11pt">
Thank you,<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11pt">Yongzhong<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:11pt"><u></u> <u></u></span></p>
<div id="m_-7305460844614688919m_2643581689528628203m_6159968039748304261mail-editor-reference-message-container">
<div>
<div style="border-right:none;border-bottom:none;border-left:none;border-top:1pt solid rgb(181,196,223);padding:3pt 0cm 0cm">
<p class="MsoNormal" style="margin-bottom:12pt"><b><span style="color:black">From:
</span></b><span style="color:black">Pierre Jolivet <<a href="mailto:pierre@joliv.et" target="_blank">pierre@joliv.et</a>><br>
<b>Date: </b>Sunday, June 23, 2024 at 12:41</span><span style="font-family:Arial,sans-serif;color:black"> </span><span style="color:black">AM<br>
<b>To: </b>Yongzhong Li <<a href="mailto:yongzhong.li@mail.utoronto.ca" target="_blank">yongzhong.li@mail.utoronto.ca</a>><br>
<b>Cc: </b><a href="mailto:petsc-users@mcs.anl.gov" target="_blank">petsc-users@mcs.anl.gov</a> <<a href="mailto:petsc-users@mcs.anl.gov" target="_blank">petsc-users@mcs.anl.gov</a>><br>
<b>Subject: </b>Re: [petsc-users] [petsc-maint] Assistance Needed with PETSc KSPSolve Performance Issue<u></u><u></u></span></p>
</div>
<p class="MsoNormal"><u></u> <u></u></p>
<div>
<p class="MsoNormal"><br>
<br>
</p>
<blockquote style="margin-top:5pt;margin-bottom:5pt">
<div>
<p class="MsoNormal">On 23 Jun 2024, at 4:07<span style="font-family:Arial,sans-serif"> </span>AM, Yongzhong Li <<a href="mailto:yongzhong.li@mail.utoronto.ca" target="_blank">yongzhong.li@mail.utoronto.ca</a>> wrote:</p>
</div>
<p class="MsoNormal"><u></u> <u></u></p>
<div>
<div id="m_-7305460844614688919m_2643581689528628203m_6159968039748304261pfptBannerfa2qtuo">
<div id="m_-7305460844614688919m_2643581689528628203m_6159968039748304261pfptBannerfa2qtuo">
<div id="m_-7305460844614688919m_2643581689528628203m_6159968039748304261pfptBannerfa2qtuo">
<p class="MsoNormal"><span style="font-family:Arial,sans-serif">This Message Is From an External Sender<u></u><u></u></span></p>
</div>
<div id="m_-7305460844614688919m_2643581689528628203m_6159968039748304261pfptBannerfa2qtuo">
<p class="MsoNormal"><span style="font-family:Arial,sans-serif">This message came from outside your organization.<u></u><u></u></span></p>
</div>
</div>
</div>
<div>
<div>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11pt">Yeah, I ran my program again using -mat_view::ascii_info and set MKL_VERBOSE to be 1, then I noticed the outputs suggested that the matrix to be seqaijmkl type (I’ve attached a few as below)<br>
<br>
--> Setting up matrix-vector products...</span></p>
</div>
<div>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11pt"> </span></p>
</div>
<div>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11pt">Mat Object: 1 MPI process</span></p>
</div>
<div>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11pt"> type: seqaijmkl</span></p>
</div>
<div>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11pt"> rows=16490, cols=35937</span></p>
</div>
<div>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11pt"> total: nonzeros=128496, allocated nonzeros=128496</span></p>
</div>
<div>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11pt"> total number of mallocs used during MatSetValues calls=0</span></p>
</div>
<div>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11pt"> not using I-node routines</span></p>
</div>
<div>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11pt">Mat Object: 1 MPI process</span></p>
</div>
<div>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11pt"> type: seqaijmkl</span></p>
</div>
<div>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11pt"> rows=16490, cols=35937</span></p>
</div>
<div>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11pt"> total: nonzeros=128496, allocated nonzeros=128496</span></p>
</div>
<div>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11pt"> total number of mallocs used during MatSetValues calls=0</span></p>
</div>
<div>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11pt"> not using I-node routines</span></p>
</div>
<div>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11pt"> </span></p>
</div>
<div>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11pt">--> Solving the system...</span></p>
</div>
<div>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11pt"> </span></p>
</div>
<div>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11pt">Excitation 1 of 1...</span></p>
</div>
<div>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11pt"> </span></p>
</div>
<div>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11pt">================================================</span></p>
</div>
<div>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11pt">Iterative solve completed in 7435 ms.</span></p>
</div>
<div>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11pt">CONVERGED: rtol.</span></p>
</div>
<div>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11pt">Iterations: 72</span></p>
</div>
<div>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11pt">Final relative residual norm: 9.22287e-07</span></p>
</div>
<div>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11pt">================================================</span></p>
</div>
<div>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11pt">[CPU TIME] System solution: 2.27160000e+02 s.</span></p>
</div>
<div>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11pt">[WALL TIME] System solution: 7.44387218e+00 s.<br>
<br>
However, it seems to me that there were still no MKL outputs even I set MKL_VERBOSE to be 1. Although, I think it should be many spmv operations when doing KSPSolve(). Do you see the possible reasons?</span></p>
</div>
</div>
</div>
</blockquote>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<p class="MsoNormal">SPMV are not reported with MKL_VERBOSE (last I checked), only dense BLAS is.</p>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<p class="MsoNormal">Thanks,</p>
</div>
<div>
<p class="MsoNormal">Pierre</p>
</div>
<p class="MsoNormal"><br>
<br>
</p>
<blockquote style="margin-top:5pt;margin-bottom:5pt">
<div>
<div>
<div>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11pt">Thanks,</span></p>
</div>
<div>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11pt">Yongzhong</span></p>
</div>
<div>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11pt"> </span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-size:11pt"> </span></p>
</div>
<div id="m_-7305460844614688919m_2643581689528628203m_6159968039748304261mail-editor-reference-message-container">
<div>
<div style="border-right:none currentcolor;border-bottom:none currentcolor;border-left:none currentcolor;border-top:1pt solid currentcolor;padding:3pt 0cm 0cm">
<p class="MsoNormal" style="margin-bottom:12pt"><b>From:<span> </span></b>Matthew Knepley <<a href="mailto:knepley@gmail.com" target="_blank">knepley@gmail.com</a>><br>
<b>Date:<span> </span></b>Saturday, June 22, 2024 at 5:56<span style="font-family:Arial,sans-serif"> </span>PM<br>
<b>To:<span> </span></b>Yongzhong Li <<a href="mailto:yongzhong.li@mail.utoronto.ca" target="_blank">yongzhong.li@mail.utoronto.ca</a>><br>
<b>Cc:<span> </span></b>Junchao Zhang <<a href="mailto:junchao.zhang@gmail.com" target="_blank">junchao.zhang@gmail.com</a>>, Pierre Jolivet <<a href="mailto:pierre@joliv.et" target="_blank">pierre@joliv.et</a>>,<span> </span><a href="mailto:petsc-users@mcs.anl.gov" target="_blank">petsc-users@mcs.anl.gov</a><span> </span><<a href="mailto:petsc-users@mcs.anl.gov" target="_blank">petsc-users@mcs.anl.gov</a>><br>
<b>Subject:<span> </span></b>Re: [petsc-users] [petsc-maint] Assistance Needed with PETSc KSPSolve Performance Issue<u></u><u></u></p>
</div>
<table border="0" cellspacing="0" cellpadding="0" align="left" width="100%" style="width:100%;display:table;float:none">
<tbody>
<tr>
<td style="background:rgb(166,166,166);padding:5.25pt 1.5pt"></td>
<td width="100%" style="width:100%;background:rgb(234,234,234);padding:5.25pt 3.75pt 5.25pt 11.25pt">
<div>
<div>
<p class="MsoNormal">
<span lang="ZH-CN" style="font-size:9pt;font-family:DengXian;color:rgb(33,33,33)">你通常不会收到来自</span><span><span style="font-size:9pt;font-family:"Segoe UI",sans-serif;color:rgb(33,33,33)"> </span></span><span style="color:black"><a href="mailto:knepley@gmail.com" target="_blank"><span style="font-size:9pt;font-family:"Segoe UI",sans-serif">knepley@gmail.com</span></a></span><span><span style="font-size:9pt;font-family:"Segoe UI",sans-serif;color:rgb(33,33,33)"> </span></span><span lang="ZH-CN" style="font-size:9pt;font-family:DengXian;color:rgb(33,33,33)">的电子邮件。</span><span style="color:black"><a href="https://urldefense.us/v3/__https://aka.ms/LearnAboutSenderIdentification__;!!G_uCfscf7eWS!fVvbGldqcUV5ju4jpu5oGmt-VjITi5JpCJzhHxpbgsERLVYZzglpxKOOyrBRGxjRxp7vWHwt3SnINFOQErR1Z8kcDcf3qwbYRxM$" target="_blank"><span lang="ZH-CN" style="font-size:9pt;font-family:DengXian">了解这一点为什么很重要</span></a></span></p>
</div>
</div>
</td>
<td width="75" style="width:56.25pt;background:rgb(234,234,234);padding:5.25pt 3.75pt">
</td>
</tr>
</tbody>
</table>
<div>
<div>
<div>
<div>
<p class="MsoNormal">On Sat, Jun 22, 2024 at 5:03<span style="font-family:Arial,sans-serif"> </span>PM Yongzhong Li <<a href="mailto:yongzhong.li@mail.utoronto.ca" target="_blank">yongzhong.li@mail.utoronto.ca</a>> wrote:<u></u><u></u></p>
</div>
</div>
<div>
<blockquote style="border-top:none currentcolor;border-right:none currentcolor;border-bottom:none currentcolor;border-left:1pt solid rgb(204,204,204);padding:0cm 0cm 0cm 6pt;margin:5pt 0cm 5pt 4.8pt">
<div>
<div>
<div>
<p class="MsoNormal"><span style="font-size:1pt;color:white">MKL_VERBOSE=1 ./ex1 matrix nonzeros = 100, allocated nonzeros = 100 MKL_VERBOSE Intel(R) MKL 2019.</span><span style="font-size:1pt;font-family:Arial,sans-serif;color:white"> </span><span style="font-size:1pt;color:white">0
Update 4 Product build 20190411 for Intel(R) 64 architecture Intel(R) Advanced Vector Extensions 512 (Intel(R) AVX-512) with support of Vector</span></p>
</div>
</div>
<div>
<div>
<p class="MsoNormal"><span style="font-size:1pt;color:white">ZjQcmQRYFpfptBannerStart</span></p>
</div>
</div>
<div style="border-right:none currentcolor;border-bottom:none currentcolor;border-left:none currentcolor;border-top:3pt solid currentcolor;padding:6pt 0cm 0cm;margin-top:12pt;margin-bottom:12pt;border-radius:4px;min-width:200px" id="m_-7305460844614688919m_2643581689528628203m_6159968039748304261m_1825240119000957322pfptBanner3e1hrxr">
<div style="margin-bottom:0.75pt;float:left" id="m_-7305460844614688919m_2643581689528628203m_6159968039748304261m_1825240119000957322pfptBanner3e1hrxr">
<div id="m_-7305460844614688919m_2643581689528628203m_6159968039748304261m_1825240119000957322pfptBanner3e1hrxr">
<div>
<p class="MsoNormal" style="line-height:13.5pt;background:rgb(208,216,220)"><b><span style="font-size:10.5pt;font-family:Arial,sans-serif;color:black">This Message Is From an External Sender</span></b></p>
</div>
</div>
<div style="margin-top:1.5pt" id="m_-7305460844614688919m_2643581689528628203m_6159968039748304261m_1825240119000957322pfptBanner3e1hrxr">
<div>
<p class="MsoNormal" style="line-height:13.5pt;background:rgb(208,216,220)"><span style="font-size:9pt;font-family:Arial,sans-serif;color:black">This message came from outside your organization.</span></p>
</div>
</div>
</div>
<div>
<div>
<p class="MsoNormal" style="background:rgb(208,216,220)"><span style="font-size:1pt;color:black"> </span></p>
</div>
</div>
</div>
<div>
<div>
<p class="MsoNormal"><span style="font-size:1pt;color:white">ZjQcmQRYFpfptBannerEnd</span></p>
</div>
</div>
<div>
<div>
<div>
<p class="MsoNormal"><span style="font-family:"Courier New"">MKL_VERBOSE=1 ./ex1</span></p>
</div>
<div>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11pt"><br>
</span><span style="font-size:11pt">matrix nonzeros = 100, allocated nonzeros = 100</span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-size:11pt">MKL_VERBOSE Intel(R) MKL 2019.0 Update 4 Product build 20190411 for Intel(R) 64 architecture Intel(R) Advanced Vector Extensions 512 (Intel(R) AVX-512) with support of Vector Neural Network Instructions enabled
processors, Lnx 2.50GHz lp64 gnu_thread</span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-size:11pt">MKL_VERBOSE ZGEMV(N,10,10,0x7ffd9d7078f0,0x187eb20,10,0x187f7c0,1,0x7ffd9d707900,0x187ff70,1) 167.34ms CNR:OFF Dyn:1 FastMM:1 TID:0 NThr:1</span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-size:11pt">MKL_VERBOSE ZSYTRF(L,10,0x1894b50,10,0x1893df0,0x7ffd9d7078c0,-1,0) 77.19ms CNR:OFF Dyn:1 FastMM:1 TID:0 NThr:1</span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-size:11pt">MKL_VERBOSE ZSYTRF(L,10,0x1894b50,10,0x1893df0,0x1894490,10,0) 83.97ms CNR:OFF Dyn:1 FastMM:1 TID:0 NThr:1</span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-size:11pt">MKL_VERBOSE ZSYTRS(L,10,1,0x1894b50,10,0x1893df0,0x1880720,10,0) 44.94ms CNR:OFF Dyn:1 FastMM:1 TID:0 NThr:1</span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-size:11pt">MKL_VERBOSE ZAXPY(10,0x7ffd9d7078f0,0x187f7c0,1,0x1880720,1) 20.72us CNR:OFF Dyn:1 FastMM:1 TID:0 NThr:1</span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-size:11pt">MKL_VERBOSE ZSYTRS(L,10,2,0x1894b50,10,0x1893df0,0x187d2a0,10,0) 4.22us CNR:OFF Dyn:1 FastMM:1 TID:0 NThr:1</span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-size:11pt">MKL_VERBOSE ZGEMM(N,N,10,2,10,0x7ffd9d707790,0x187eb20,10,0x187d2a0,10,0x7ffd9d7077a0,0x1896a70,10) 1.41ms CNR:OFF Dyn:1 FastMM:1 TID:0 NThr:1</span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-size:11pt">MKL_VERBOSE ZAXPY(20,0x7ffd9d7078a0,0x1896a70,1,0x187b650,1) 381ns CNR:OFF Dyn:1 FastMM:1 TID:0 NThr:1</span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-size:11pt">MKL_VERBOSE ZSYTRF(L,10,0x1894b50,10,0x1893df0,0x7ffd9d707840,-1,0) 742ns CNR:OFF Dyn:1 FastMM:1 TID:0 NThr:1</span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-size:11pt">MKL_VERBOSE ZSYTRF(L,10,0x1894b50,10,0x1893df0,0x18951a0,10,0) 4.20us CNR:OFF Dyn:1 FastMM:1 TID:0 NThr:1</span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-size:11pt">MKL_VERBOSE ZSYTRS(L,10,1,0x1894b50,10,0x1893df0,0x1880720,10,0) 2.94us CNR:OFF Dyn:1 FastMM:1 TID:0 NThr:1</span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-size:11pt">MKL_VERBOSE ZAXPY(10,0x7ffd9d7078f0,0x187f7c0,1,0x1880720,1) 292ns CNR:OFF Dyn:1 FastMM:1 TID:0 NThr:1</span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-size:11pt">MKL_VERBOSE ZGEMV(N,10,10,0x7ffd9d7078f0,0x187eb20,10,0x187f7c0,1,0x7ffd9d707900,0x187ff70,1) 1.17us CNR:OFF Dyn:1 FastMM:1 TID:0 NThr:1</span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-size:11pt">MKL_VERBOSE ZGETRF(10,10,0x1894b50,10,0x1893df0,0) 202.48ms CNR:OFF Dyn:1 FastMM:1 TID:0 NThr:1</span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-size:11pt">MKL_VERBOSE ZGETRS(N,10,1,0x1894b50,10,0x1893df0,0x1880720,10,0) 20.78ms CNR:OFF Dyn:1 FastMM:1 TID:0 NThr:1</span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-size:11pt">MKL_VERBOSE ZAXPY(10,0x7ffd9d7078f0,0x187f7c0,1,0x1880720,1) 954ns CNR:OFF Dyn:1 FastMM:1 TID:0 NThr:1</span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-size:11pt">MKL_VERBOSE ZGETRS(N,10,2,0x1894b50,10,0x1893df0,0x187d2a0,10,0) 30.74ms CNR:OFF Dyn:1 FastMM:1 TID:0 NThr:1</span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-size:11pt">MKL_VERBOSE ZGEMM(N,N,10,2,10,0x7ffd9d707790,0x187eb20,10,0x187d2a0,10,0x7ffd9d7077a0,0x18969c0,10) 3.95us CNR:OFF Dyn:1 FastMM:1 TID:0 NThr:1</span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-size:11pt">MKL_VERBOSE ZAXPY(20,0x7ffd9d7078a0,0x18969c0,1,0x187b650,1) 995ns CNR:OFF Dyn:1 FastMM:1 TID:0 NThr:1</span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-size:11pt">MKL_VERBOSE ZGETRF(10,10,0x1894b50,10,0x1893df0,0) 4.09us CNR:OFF Dyn:1 FastMM:1 TID:0 NThr:1</span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-size:11pt">MKL_VERBOSE ZGETRS(N,10,1,0x1894b50,10,0x1893df0,0x1880720,10,0) 3.92us CNR:OFF Dyn:1 FastMM:1 TID:0 NThr:1</span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-size:11pt">MKL_VERBOSE ZAXPY(10,0x7ffd9d7078f0,0x187f7c0,1,0x1880720,1) 274ns CNR:OFF Dyn:1 FastMM:1 TID:0 NThr:1</span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-size:11pt">MKL_VERBOSE ZGEMV(N,15,10,0x7ffd9d7078f0,0x187ec70,15,0x187fc30,1,0x7ffd9d707900,0x1880400,1) 1.59us CNR:OFF Dyn:1 FastMM:1 TID:0 NThr:1</span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-size:11pt">MKL_VERBOSE ZGEQRF(15,10,0x1894b40,15,0x1894550,0x7ffd9d707900,-1,0) 47.07us CNR:OFF Dyn:1 FastMM:1 TID:0 NThr:1</span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-size:11pt">MKL_VERBOSE ZGEQRF(15,10,0x1894b40,15,0x1894550,0x1895cb0,10,0) 26.62us CNR:OFF Dyn:1 FastMM:1 TID:0 NThr:1</span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-size:11pt">MKL_VERBOSE ZUNMQR(L,C,15,1,10,0x1894b40,15,0x1894550,0x1895b00,15,0x7ffd9d7078b0,-1,0) 35.32us CNR:OFF Dyn:1 FastMM:1 TID:0 NThr:1</span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-size:11pt">MKL_VERBOSE ZUNMQR(L,C,15,1,10,0x1894b40,15,0x1894550,0x1895b00,15,0x1895cb0,10,0) 42.33ms CNR:OFF Dyn:1 FastMM:1 TID:0 NThr:1</span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-size:11pt">MKL_VERBOSE ZTRTRS(U,N,N,10,1,0x1894b40,15,0x1895b00,15,0) 16.11us CNR:OFF Dyn:1 FastMM:1 TID:0 NThr:1</span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-size:11pt">MKL_VERBOSE ZAXPY(10,0x7ffd9d7078f0,0x187fc30,1,0x1880c70,1) 395ns CNR:OFF Dyn:1 FastMM:1 TID:0 NThr:1</span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-size:11pt">MKL_VERBOSE ZGEMM(N,N,15,2,10,0x7ffd9d707790,0x187ec70,15,0x187d310,10,0x7ffd9d7077a0,0x187b5b0,15) 3.22us CNR:OFF Dyn:1 FastMM:1 TID:0 NThr:1</span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-size:11pt">MKL_VERBOSE ZUNMQR(L,C,15,2,10,0x1894b40,15,0x1894550,0x1897760,15,0x7ffd9d7078c0,-1,0) 730ns CNR:OFF Dyn:1 FastMM:1 TID:0 NThr:1</span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-size:11pt">MKL_VERBOSE ZUNMQR(L,C,15,2,10,0x1894b40,15,0x1894550,0x1897760,15,0x1895cb0,10,0) 4.42us CNR:OFF Dyn:1 FastMM:1 TID:0 NThr:1</span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-size:11pt">MKL_VERBOSE ZTRTRS(U,N,N,10,2,0x1894b40,15,0x1897760,15,0) 5.96us CNR:OFF Dyn:1 FastMM:1 TID:0 NThr:1</span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-size:11pt">MKL_VERBOSE ZAXPY(20,0x7ffd9d7078a0,0x187d310,1,0x1897610,1) 222ns CNR:OFF Dyn:1 FastMM:1 TID:0 NThr:1</span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-size:11pt">MKL_VERBOSE ZGEQRF(15,10,0x1894b40,15,0x18954b0,0x7ffd9d707820,-1,0) 685ns CNR:OFF Dyn:1 FastMM:1 TID:0 NThr:1</span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-size:11pt">MKL_VERBOSE ZGEQRF(15,10,0x1894b40,15,0x18954b0,0x1895d60,10,0) 6.11us CNR:OFF Dyn:1 FastMM:1 TID:0 NThr:1</span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-size:11pt">MKL_VERBOSE ZUNMQR(L,C,15,1,10,0x1894b40,15,0x18954b0,0x1895bb0,15,0x7ffd9d7078b0,-1,0) 390ns CNR:OFF Dyn:1 FastMM:1 TID:0 NThr:1</span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-size:11pt">MKL_VERBOSE ZUNMQR(L,C,15,1,10,0x1894b40,15,0x18954b0,0x1895bb0,15,0x1895d60,10,0) 3.09us CNR:OFF Dyn:1 FastMM:1 TID:0 NThr:1</span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-size:11pt">MKL_VERBOSE ZTRTRS(U,N,N,10,1,0x1894b40,15,0x1895bb0,15,0) 1.05us CNR:OFF Dyn:1 FastMM:1 TID:0 NThr:1</span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-size:11pt">MKL_VERBOSE ZAXPY(10,0x7ffd9d7078f0,0x187fc30,1,0x1880c70,1) 257ns CNR:OFF Dyn:1 FastMM:1 TID:0 NThr:1</span><span lang="EN-US" style="font-size:11pt"><br>
<br>
Yes, for petsc example, there are MKL outputs, but for my own program. All I did is to change the matrix type from MATAIJ to MATAIJMKL to get optimized performance for spmv from MKL. Should I expect to see any MKL outputs in this case?</span></p>
</div>
</div>
</div>
</div>
</blockquote>
<div>
<div>
<p class="MsoNormal"> </p>
</div>
</div>
<div>
<div>
<p class="MsoNormal">Are you sure that the type changed? You can MatView() the matrix with format ascii_info to see.</p>
</div>
</div>
<div>
<div>
<p class="MsoNormal"> </p>
</div>
</div>
<div>
<div>
<p class="MsoNormal"> Thanks,</p>
</div>
</div>
<div>
<div>
<p class="MsoNormal"> </p>
</div>
</div>
<div>
<div>
<p class="MsoNormal"> Matt</p>
</div>
</div>
<div>
<div>
<p class="MsoNormal"> </p>
</div>
</div>
<div>
<div>
<p class="MsoNormal"> <u></u><u></u></p>
</div>
</div>
<blockquote style="border-top:none currentcolor;border-right:none currentcolor;border-bottom:none currentcolor;border-left:1pt solid rgb(204,204,204);padding:0cm 0cm 0cm 6pt;margin:5pt 0cm 5pt 4.8pt">
<div>
<div>
<div>
<div>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11pt">Thanks,</span></p>
</div>
<div>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11pt">Yongzhong</span></p>
</div>
<div>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11pt"> </span></p>
</div>
<div id="m_-7305460844614688919m_2643581689528628203m_6159968039748304261m_1825240119000957322mail-editor-reference-message-container">
<div>
<div style="border-right:none currentcolor;border-bottom:none currentcolor;border-left:none currentcolor;border-top:1pt solid currentcolor;padding:3pt 0cm 0cm">
<p class="MsoNormal" style="margin-bottom:12pt"><b>From:<span> </span></b>Junchao Zhang <<a href="mailto:junchao.zhang@gmail.com" target="_blank">junchao.zhang@gmail.com</a>><br>
<b>Date:<span> </span></b>Saturday, June 22, 2024 at 9:40<span style="font-family:Arial,sans-serif"> </span>AM<br>
<b>To:<span> </span></b>Yongzhong Li <<a href="mailto:yongzhong.li@mail.utoronto.ca" target="_blank">yongzhong.li@mail.utoronto.ca</a>><br>
<b>Cc:<span> </span></b>Pierre Jolivet <<a href="mailto:pierre@joliv.et" target="_blank">pierre@joliv.et</a>>,<span> </span><a href="mailto:petsc-users@mcs.anl.gov" target="_blank">petsc-users@mcs.anl.gov</a><span> </span><<a href="mailto:petsc-users@mcs.anl.gov" target="_blank">petsc-users@mcs.anl.gov</a>><br>
<b>Subject:<span> </span></b>Re: [petsc-users] [petsc-maint] Assistance Needed with PETSc KSPSolve Performance Issue</p>
</div>
<div>
<div>
<p class="MsoNormal">No, you don't. It is strange. Perhaps you can you run a petsc example first and see if MKL is really used<u></u><u></u></p>
</div>
<div>
<div>
<p class="MsoNormal"><span style="font-family:"Courier New"">$ cd src/mat/tests</span></p>
</div>
</div>
<div>
<div>
<p class="MsoNormal"><span style="font-family:"Courier New"">$ make ex1</span></p>
</div>
</div>
<div>
<div>
<div>
<p class="MsoNormal"><span style="font-family:"Courier New"">$ MKL_VERBOSE=1 ./ex1</span></p>
</div>
</div>
<div>
<div>
<p class="MsoNormal"><br clear="all">
</p>
</div>
<div>
<div>
<div>
<div>
<p class="MsoNormal">--Junchao Zhang</p>
</div>
</div>
</div>
</div>
<div>
<p class="MsoNormal"> </p>
</div>
</div>
</div>
</div>
<div>
<p class="MsoNormal"> </p>
</div>
<div>
<div>
<div>
<p class="MsoNormal">On Fri, Jun 21, 2024 at 4:03<span style="font-family:Arial,sans-serif"> </span>PM Yongzhong Li <<a href="mailto:yongzhong.li@mail.utoronto.ca" target="_blank">yongzhong.li@mail.utoronto.ca</a>> wrote:<u></u><u></u></p>
</div>
</div>
<blockquote style="border-top:none currentcolor;border-right:none currentcolor;border-bottom:none currentcolor;border-left:1pt solid rgb(204,204,204);padding:0cm 0cm 0cm 6pt;margin:5pt 0cm 5pt 4.8pt">
<div>
<div>
<div>
<div>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11pt">I am using<br>
<br>
export MKL_VERBOSE=1</span></p>
</div>
<div>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11pt">./xx<br>
<br>
in the bash file, do I have to use -</span><span><span lang="EN-US"> </span></span>ksp_converged_reason<span lang="EN-US">?<br>
<br>
Thanks,</span></p>
</div>
<div>
<p class="MsoNormal"><span lang="EN-US">Yongzhong</span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-size:11pt"> </span></p>
</div>
<div id="m_-7305460844614688919m_2643581689528628203m_6159968039748304261m_1825240119000957322m_-6075196701903926580mail-editor-reference-message-container">
<div>
<div style="border-right:none currentcolor;border-bottom:none currentcolor;border-left:none currentcolor;border-top:1pt solid currentcolor;padding:3pt 0cm 0cm">
<p class="MsoNormal" style="margin-bottom:12pt"><b>From:<span> </span></b>Pierre Jolivet <<a href="mailto:pierre@joliv.et" target="_blank">pierre@joliv.et</a>><br>
<b>Date:<span> </span></b>Friday, June 21, 2024 at 1:47<span style="font-family:Arial,sans-serif"> </span>PM<br>
<b>To:<span> </span></b>Yongzhong Li <<a href="mailto:yongzhong.li@mail.utoronto.ca" target="_blank">yongzhong.li@mail.utoronto.ca</a>><br>
<b>Cc:<span> </span></b>Junchao Zhang <<a href="mailto:junchao.zhang@gmail.com" target="_blank">junchao.zhang@gmail.com</a>>,<span> </span><a href="mailto:petsc-users@mcs.anl.gov" target="_blank">petsc-users@mcs.anl.gov</a><span> </span><<a href="mailto:petsc-users@mcs.anl.gov" target="_blank">petsc-users@mcs.anl.gov</a>><br>
<b>Subject:<span> </span></b>Re: [petsc-users] [petsc-maint] Assistance Needed with PETSc KSPSolve Performance Issue<u></u><u></u></p>
</div>
<table border="0" cellspacing="0" cellpadding="0" align="left" width="100%" style="width:100%;display:table;float:none">
<tbody>
<tr>
<td style="background:rgb(166,166,166);padding:5.25pt 1.5pt"></td>
<td width="100%" style="width:100%;background:rgb(234,234,234);padding:5.25pt 3.75pt 5.25pt 11.25pt">
<div>
<div>
<p class="MsoNormal">
<span lang="ZH-CN" style="font-size:9pt;font-family:DengXian;color:rgb(33,33,33)">你通常不会收到来自</span><span><span style="font-size:9pt;font-family:"Segoe UI",sans-serif;color:rgb(33,33,33)"> </span></span><span style="color:black"><a href="mailto:pierre@joliv.et" target="_blank"><span style="font-size:9pt;font-family:"Segoe UI",sans-serif">pierre@joliv.et</span></a></span><span><span style="font-size:9pt;font-family:"Segoe UI",sans-serif;color:rgb(33,33,33)"> </span></span><span lang="ZH-CN" style="font-size:9pt;font-family:DengXian;color:rgb(33,33,33)">的电子邮件。</span><span style="color:black"><a href="https://urldefense.us/v3/__https://aka.ms/LearnAboutSenderIdentification__;!!G_uCfscf7eWS!flsZMI97ne0yyxHhLda3hROB9qsgstuZS-jPinxGIzFCCSdn1ujdoMR8dyz-5_kVqqMM-12Lt0dTdjKrx3wXhHZmBhNydvFQeSY$" target="_blank"><span lang="ZH-CN" style="font-size:9pt;font-family:DengXian">了解这一点为什么很重要</span></a></span></p>
</div>
</div>
</td>
<td width="75" style="width:56.25pt;background:rgb(234,234,234);padding:5.25pt 3.75pt">
</td>
</tr>
</tbody>
</table>
<div>
<div>
<p class="MsoNormal">How do you set the variable?</p>
</div>
<div>
<div>
<p class="MsoNormal"> </p>
</div>
</div>
<div>
<div>
<div>
<p class="MsoNormal">$ MKL_VERBOSE=1 ./ex1 -ksp_converged_reason</p>
</div>
</div>
<div>
<div>
<p class="MsoNormal">MKL_VERBOSE oneMKL 2024.0 Update 1 Product build 20240215 for Intel(R) 64 architecture Intel(R) Advanced Vector Extensions 2 (Intel(R) AVX2) enabled processors, Lnx 2.80GHz lp64 intel_thread</p>
</div>
</div>
<div>
<div>
<p class="MsoNormal">MKL_VERBOSE DDOT(10,0x22127c0,1,0x22127c0,1) 2.02ms CNR:OFF Dyn:1 FastMM:1 TID:0 NThr:1</p>
</div>
</div>
<div>
<div>
<p class="MsoNormal">MKL_VERBOSE DSCAL(10,0x7ffc9fb4ff08,0x22127c0,1) 12.67us CNR:OFF Dyn:1 FastMM:1 TID:0 NThr:1</p>
</div>
</div>
<div>
<div>
<p class="MsoNormal">MKL_VERBOSE DDOT(10,0x22127c0,1,0x2212840,1) 1.52us CNR:OFF Dyn:1 FastMM:1 TID:0 NThr:1</p>
</div>
</div>
<div>
<div>
<p class="MsoNormal">MKL_VERBOSE DDOT(10,0x2212840,1,0x2212840,1) 167ns CNR:OFF Dyn:1 FastMM:1 TID:0 NThr:1</p>
</div>
</div>
<div>
<div>
<p class="MsoNormal">[...]</p>
</div>
</div>
<div>
<p class="MsoNormal" style="margin-bottom:12pt"> </p>
<blockquote style="margin-top:5pt;margin-bottom:5pt">
<div>
<div>
<p class="MsoNormal">On 21 Jun 2024, at 7:37<span style="font-family:Arial,sans-serif"> </span>PM, Yongzhong Li <<a href="mailto:yongzhong.li@mail.utoronto.ca" target="_blank">yongzhong.li@mail.utoronto.ca</a>> wrote:<u></u><u></u></p>
</div>
</div>
<div>
<p class="MsoNormal"> </p>
</div>
<div>
<div id="m_-7305460844614688919m_2643581689528628203m_6159968039748304261m_1825240119000957322m_-6075196701903926580pfptBannermu5fqka">
<div id="m_-7305460844614688919m_2643581689528628203m_6159968039748304261m_1825240119000957322m_-6075196701903926580pfptBannermu5fqka">
<div id="m_-7305460844614688919m_2643581689528628203m_6159968039748304261m_1825240119000957322m_-6075196701903926580pfptBannermu5fqka">
<div>
<p class="MsoNormal"><span style="font-family:Arial,sans-serif">This Message Is From an External Sender<span> </span></span></p>
</div>
</div>
<div id="m_-7305460844614688919m_2643581689528628203m_6159968039748304261m_1825240119000957322m_-6075196701903926580pfptBannermu5fqka">
<div>
<p class="MsoNormal"><span style="font-family:Arial,sans-serif">This message came from outside your organization.</span></p>
</div>
</div>
</div>
</div>
<div>
<div>
<div>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11pt">Hello all,<br>
<br>
I set MKL_VERBOSE = 1, but observed no print output specific to the use of MKL. Does PETSc enable this verbose output?<br>
<br>
Best,</span></p>
</div>
</div>
<div>
<p class="MsoNormal" style="margin-bottom:12pt"><span lang="EN-US" style="font-size:11pt">Yongzhong</span></p>
</div>
<div>
<div>
<p class="MsoNormal"><span style="font-size:11pt"> </span><u></u><u></u></p>
</div>
</div>
<div id="m_-7305460844614688919m_2643581689528628203m_6159968039748304261m_1825240119000957322m_-6075196701903926580mail-editor-reference-message-container">
<div>
<div>
<p class="MsoNormal" style="margin-bottom:12pt"><b>From:<span> </span></b>Pierre Jolivet <<a href="mailto:pierre@joliv.et" target="_blank">pierre@joliv.et</a>><br>
<b>Date:<span> </span></b>Friday, June 21, 2024 at 1:36<span style="font-family:Arial,sans-serif"> </span>AM<br>
<b>To:<span> </span></b>Junchao Zhang <<a href="mailto:junchao.zhang@gmail.com" target="_blank">junchao.zhang@gmail.com</a>><br>
<b>Cc:<span> </span></b>Yongzhong Li <<a href="mailto:yongzhong.li@mail.utoronto.ca" target="_blank">yongzhong.li@mail.utoronto.ca</a>>,<span> </span><a href="mailto:petsc-users@mcs.anl.gov" target="_blank">petsc-users@mcs.anl.gov</a><span> </span><<a href="mailto:petsc-users@mcs.anl.gov" target="_blank">petsc-users@mcs.anl.gov</a>><br>
<b>Subject:<span> </span></b>Re: [petsc-users] [petsc-maint] Assistance Needed with PETSc KSPSolve Performance Issue<u></u><u></u></p>
</div>
<table border="0" cellspacing="0" cellpadding="0" align="left" width="100%" style="width:100%;display:table;float:none">
<tbody>
<tr>
<td style="background:rgb(166,166,166);padding:5.25pt 1.5pt"></td>
<td width="100%" style="width:100%;background:rgb(234,234,234);padding:5.25pt 3.75pt 5.25pt 11.25pt">
<div>
<div>
<div>
<p class="MsoNormal">
<span lang="ZH-CN" style="font-size:9pt;font-family:DengXian;color:rgb(33,33,33)">你通常不会收到来自</span><span><span style="font-size:9pt;font-family:"Segoe UI",sans-serif;color:rgb(33,33,33)"> </span></span><span style="color:black"><a href="mailto:pierre@joliv.et" target="_blank"><span style="font-size:9pt;font-family:"Segoe UI",sans-serif">pierre@joliv.et</span></a></span><span><span style="font-size:9pt;font-family:"Segoe UI",sans-serif;color:rgb(33,33,33)"> </span></span><span lang="ZH-CN" style="font-size:9pt;font-family:DengXian;color:rgb(33,33,33)">的电子邮件。</span><span style="color:black"><a href="https://urldefense.us/v3/__https://aka.ms/LearnAboutSenderIdentification__;!!G_uCfscf7eWS!eXBeeIXo9Yqgp2nypqwKYimLnGBZXnF4dXxgLM1UoOIO6n8nt3XlfgjVWLPWJh4UOa5NNpx-nrJb_H828XRQKUREfR2m69oCbxI$" target="_blank"><span lang="ZH-CN" style="font-size:9pt;font-family:DengXian">了解这一点为什么很重要</span></a></span></p>
</div>
</div>
</div>
</td>
<td width="75" style="width:56.25pt;background:rgb(234,234,234);padding:5.25pt 3.75pt">
</td>
</tr>
</tbody>
</table>
<div>
<div>
<div>
<p class="MsoNormal"> </p>
</div>
</div>
<div>
<div>
<p class="MsoNormal" style="margin-bottom:12pt"> </p>
</div>
<blockquote style="margin-top:5pt;margin-bottom:5pt">
<div>
<div>
<div>
<p class="MsoNormal">On 21 Jun 2024, at 6:42<span style="font-family:Arial,sans-serif"> </span>AM, Junchao Zhang <<a href="mailto:junchao.zhang@gmail.com" target="_blank">junchao.zhang@gmail.com</a>> wrote:<u></u><u></u></p>
</div>
</div>
</div>
<div>
<div>
<p class="MsoNormal"> </p>
</div>
</div>
<div>
<div id="m_-7305460844614688919m_2643581689528628203m_6159968039748304261m_1825240119000957322m_-6075196701903926580pfptBanners71iuax">
<div id="m_-7305460844614688919m_2643581689528628203m_6159968039748304261m_1825240119000957322m_-6075196701903926580pfptBanners71iuax">
<div id="m_-7305460844614688919m_2643581689528628203m_6159968039748304261m_1825240119000957322m_-6075196701903926580pfptBanners71iuax">
<div>
<div>
<p class="MsoNormal"><span style="font-family:Arial,sans-serif">This Message Is From an External Sender</span></p>
</div>
</div>
</div>
<div id="m_-7305460844614688919m_2643581689528628203m_6159968039748304261m_1825240119000957322m_-6075196701903926580pfptBanners71iuax">
<div>
<div>
<p class="MsoNormal"><span style="font-family:Arial,sans-serif">This message came from outside your organization.</span></p>
</div>
</div>
</div>
</div>
</div>
<div>
<div>
<div>
<p class="MsoNormal"><span style="font-size:9pt;font-family:Helvetica">I remember there are some MKL env vars to print MKL routines called. </span></p>
</div>
</div>
</div>
</div>
</blockquote>
<div>
<div>
<div>
<p class="MsoNormal"> </p>
</div>
</div>
</div>
<div>
<div>
<div>
<p class="MsoNormal">The environment variable is MKL_VERBOSE</p>
</div>
</div>
</div>
<div>
<div>
<div>
<p class="MsoNormal"> </p>
</div>
</div>
</div>
<div>
<div>
<div>
<p class="MsoNormal">Thanks,</p>
</div>
</div>
</div>
<div>
<div>
<div>
<p class="MsoNormal">Pierre</p>
</div>
</div>
</div>
<div>
<p class="MsoNormal" style="margin-bottom:12pt"> </p>
</div>
<blockquote style="margin-top:5pt;margin-bottom:5pt">
<div>
<div>
<div>
<div>
<p class="MsoNormal"><span style="font-size:9pt;font-family:Helvetica">Maybe we can try it to see what MKL routines are really used and then we can understand why some petsc functions did not speed up </span><u></u><u></u></p>
</div>
</div>
<div>
<div>
<div>
<p class="MsoNormal"><span style="font-size:9pt;font-family:Helvetica"><br clear="all">
</span></p>
</div>
</div>
<div>
<div>
<div>
<div>
<div>
<p class="MsoNormal"><span style="font-size:9pt;font-family:Helvetica">--Junchao Zhang</span></p>
</div>
</div>
</div>
</div>
</div>
<div>
<div>
<p class="MsoNormal"><span style="font-size:9pt;font-family:Helvetica"> </span></p>
</div>
</div>
</div>
</div>
<div>
<div>
<p class="MsoNormal"> </p>
</div>
</div>
<div>
<div>
<div>
<div>
<p class="MsoNormal"><span style="font-size:9pt;font-family:Helvetica">On Thu, Jun 20, 2024 at 10:39 PM Yongzhong Li <</span><a href="mailto:yongzhong.li@mail.utoronto.ca" target="_blank"><span style="font-size:9pt;font-family:Helvetica">yongzhong.li@mail.utoronto.ca</span></a><span style="font-size:9pt;font-family:Helvetica">>
wrote:</span><u></u><u></u></p>
</div>
</div>
</div>
<blockquote style="border-top:none currentcolor;border-right:none currentcolor;border-bottom:none currentcolor;border-left:1pt solid rgb(204,204,204);padding:0cm 0cm 0cm 6pt;margin:5pt 0cm 5pt 4.8pt">
<div>
<div style="margin-top:12pt;margin-bottom:12pt;border-width:initial;border-style:initial;border-color:currentcolor;border-radius:4px;min-width:200px" id="m_-7305460844614688919m_2643581689528628203m_6159968039748304261m_1825240119000957322m_-6075196701903926580m_-4861803116733450987pfptBannerqdehbik">
<div style="margin-bottom:0.75pt;float:left" id="m_-7305460844614688919m_2643581689528628203m_6159968039748304261m_1825240119000957322m_-6075196701903926580m_-4861803116733450987pfptBannerqdehbik">
<div id="m_-7305460844614688919m_2643581689528628203m_6159968039748304261m_1825240119000957322m_-6075196701903926580m_-4861803116733450987pfptBannerqdehbik">
<div>
<div>
<p class="MsoNormal" style="line-height:13.5pt;background:rgb(208,216,220)"><b><span style="font-size:10.5pt;font-family:Arial,sans-serif;color:black">This Message Is From an External Sender</span></b></p>
</div>
</div>
</div>
<div style="margin-top:1.5pt" id="m_-7305460844614688919m_2643581689528628203m_6159968039748304261m_1825240119000957322m_-6075196701903926580m_-4861803116733450987pfptBannerqdehbik">
<div>
<div>
<p class="MsoNormal" style="line-height:13.5pt;background:rgb(208,216,220)"><span style="font-size:9pt;font-family:Arial,sans-serif;color:black">This message came from outside your organization.</span></p>
</div>
</div>
</div>
</div>
<div>
<div>
<div>
<p class="MsoNormal" style="background:rgb(208,216,220)"><span style="font-size:1pt;font-family:Helvetica;color:black"> </span></p>
</div>
</div>
</div>
</div>
<div>
<div>
<div>
<div>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11pt;font-family:Helvetica">Hi Barry, sorry for my last results. I didn’t fully understand the stage profiling and logging in PETSc, now I only record KSPSolve() stage of my program. Some sample codes
are as follow,<br>
<br>
<span> </span>// Static variable to keep track of the stage counter</span></p>
</div>
</div>
<div>
<div>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11pt;font-family:Helvetica"> <span> </span>static int stageCounter = 1;</span></p>
</div>
</div>
<div>
<div>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11pt;font-family:Helvetica"> </span></p>
</div>
</div>
<div>
<div>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11pt;font-family:Helvetica"> <span> </span>// Generate a unique stage name</span></p>
</div>
</div>
<div>
<div>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11pt;font-family:Helvetica"> <span> </span>std::ostringstream oss;</span></p>
</div>
</div>
<div>
<div>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11pt;font-family:Helvetica"> <span> </span>oss << "Stage " << stageCounter << " of Code";</span></p>
</div>
</div>
<div>
<div>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11pt;font-family:Helvetica"> <span> </span>std::string stageName = oss.str();</span></p>
</div>
</div>
<div>
<div>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11pt;font-family:Helvetica"> </span></p>
</div>
</div>
<div>
<div>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11pt;font-family:Helvetica"> <span> </span>// Register the stage</span></p>
</div>
</div>
<div>
<div>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11pt;font-family:Helvetica"> <span> </span>PetscLogStage stagenum;</span></p>
</div>
</div>
<div>
<div>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11pt;font-family:Helvetica"> </span></p>
</div>
</div>
<div>
<div>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11pt;font-family:Helvetica"> <span> </span>PetscLogStageRegister(stageName.c_str(), &stagenum);</span></p>
</div>
</div>
<div>
<div>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11pt;font-family:Helvetica"> <span> </span>PetscLogStagePush(stagenum);</span></p>
</div>
</div>
<div>
<div>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11pt;font-family:Helvetica"> </span></p>
</div>
</div>
<div>
<div>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11pt;font-family:Helvetica"> <span> </span><b>KSPSolve(*ksp_ptr, b, x);</b></span></p>
</div>
</div>
<div>
<div>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11pt;font-family:Helvetica"> </span></p>
</div>
</div>
<div>
<div>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11pt;font-family:Helvetica"> <span> </span>PetscLogStagePop();</span></p>
</div>
</div>
<div>
<div>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11pt;font-family:Helvetica"> <span> </span>stageCounter++;<br>
<br>
I have attached my new logging results, there are 1 main stage and 4 other stages where each one is KSPSolve() call.<br>
<br>
To provide some additional backgrounds, if you recall, I have been trying to get efficient iterative solution using multithreading. I found out by compiling PETSc with Intel MKL library instead of OpenBLAS, I am able to perform sparse matrix-vector multiplication
faster, I am using MATSEQAIJMKL. This makes the shell matrix vector product in each iteration scale well with the #of threads. However, I found out the total GMERS solve time (~KSPSolve() time) is not scaling well the #of threads.<br>
<br>
>From the logging results I learned that when performing KSPSolve(), there are some CPU overheads in PCApply() and KSPGMERSOrthog(). I ran my programs using different number of threads and plotted the time consumption for PCApply() and KSPGMERSOrthog() against
#of thread. I found out these two operations are not scaling with the threads at all! My results are attached as the pdf to give you a clear view.<br>
<br>
My questions is,<br>
<br>
>From my understanding, in PCApply, MatSolve() is involved, KSPGMERSOrthog() will have many vector operations, so why these two parts can’t scale well with the # of threads when the intel MKL library is linked?<br>
<br>
Thank you,<br>
Yongzhong</span></p>
</div>
</div>
<div>
<div>
<p class="MsoNormal"><span style="font-size:11pt;font-family:Helvetica"> </span></p>
</div>
</div>
<div id="m_-7305460844614688919m_2643581689528628203m_6159968039748304261m_1825240119000957322m_-6075196701903926580m_-4861803116733450987mail-editor-reference-message-container">
<div>
<div>
<p class="MsoNormal" style="margin-bottom:12pt"><b><span style="font-size:9pt;font-family:Helvetica">From:<span> </span></span></b><span style="font-size:9pt;font-family:Helvetica">Barry
Smith <</span><a href="mailto:bsmith@petsc.dev" target="_blank"><span style="font-size:9pt;font-family:Helvetica">bsmith@petsc.dev</span></a><span style="font-size:9pt;font-family:Helvetica">><br>
<b>Date:<span> </span></b>Friday, June 14, 2024 at 11:36</span><span style="font-size:9pt;font-family:Arial,sans-serif"> </span><span style="font-size:9pt;font-family:Helvetica">AM<br>
<b>To:<span> </span></b>Yongzhong Li <</span><a href="mailto:yongzhong.li@mail.utoronto.ca" target="_blank"><span style="font-size:9pt;font-family:Helvetica">yongzhong.li@mail.utoronto.ca</span></a><span style="font-size:9pt;font-family:Helvetica">><br>
<b>Cc:<span> </span></b></span><a href="mailto:petsc-users@mcs.anl.gov" target="_blank"><span style="font-size:9pt;font-family:Helvetica">petsc-users@mcs.anl.gov</span></a><span><span style="font-size:9pt;font-family:Helvetica"> </span></span><span style="font-size:9pt;font-family:Helvetica"><</span><a href="mailto:petsc-users@mcs.anl.gov" target="_blank"><span style="font-size:9pt;font-family:Helvetica">petsc-users@mcs.anl.gov</span></a><span style="font-size:9pt;font-family:Helvetica">>,<span> </span></span><a href="mailto:petsc-maint@mcs.anl.gov" target="_blank"><span style="font-size:9pt;font-family:Helvetica">petsc-maint@mcs.anl.gov</span></a><span><span style="font-size:9pt;font-family:Helvetica"> </span></span><span style="font-size:9pt;font-family:Helvetica"><</span><a href="mailto:petsc-maint@mcs.anl.gov" target="_blank"><span style="font-size:9pt;font-family:Helvetica">petsc-maint@mcs.anl.gov</span></a><span style="font-size:9pt;font-family:Helvetica">>,
Piero Triverio <</span><a href="mailto:piero.triverio@utoronto.ca" target="_blank"><span style="font-size:9pt;font-family:Helvetica">piero.triverio@utoronto.ca</span></a><span style="font-size:9pt;font-family:Helvetica">><br>
<b>Subject:<span> </span></b>Re: [petsc-maint] Assistance Needed with PETSc KSPSolve Performance Issue</span></p>
</div>
<div>
<div>
<div>
<p class="MsoNormal"><span style="font-size:9pt;font-family:Helvetica"> </span><u></u><u></u></p>
</div>
</div>
</div>
<div>
<div>
<p class="MsoNormal"><span style="font-size:9pt;font-family:Helvetica"> I am a bit confused. Without the initial guess computation, there are still a bunch of events I don't understand </span></p>
</div>
</div>
<div>
<div>
<div>
<p class="MsoNormal"><span style="font-size:9pt;font-family:Helvetica"> </span></p>
</div>
</div>
</div>
<div>
<div>
<div>
<div>
<p class="MsoNormal"><span style="font-size:9pt;font-family:Helvetica">MatTranspose 79 1.0 4.0598e+01 1.0 0.00e+00 0.0 0.0e+00 0.0e+00 0.0e+00 0 0 0 0 0 0 0 0 0 0 0</span></p>
</div>
</div>
</div>
<div>
<div>
<div>
<p class="MsoNormal"><span style="font-size:9pt;font-family:Helvetica">MatMatMultSym 110 1.0 1.7419e+02 1.0 0.00e+00 0.0 0.0e+00 0.0e+00 0.0e+00 1 0 0 0 0 1 0 0 0 0 0</span></p>
</div>
</div>
</div>
<div>
<div>
<div>
<p class="MsoNormal"><span style="font-size:9pt;font-family:Helvetica">MatMatMultNum 90 1.0 1.2640e+02 1.0 0.00e+00 0.0 0.0e+00 0.0e+00 0.0e+00 1 0 0 0 0 1 0 0 0 0 0</span></p>
</div>
</div>
</div>
<div>
<div>
<div>
<p class="MsoNormal"><span style="font-size:9pt;font-family:Helvetica">MatMatMatMultSym 20 1.0 1.3049e+02 1.0 0.00e+00 0.0 0.0e+00 0.0e+00 0.0e+00 1 0 0 0 0 1 0 0 0 0 0</span></p>
</div>
</div>
</div>
<div>
<div>
<div>
<p class="MsoNormal"><span style="font-size:9pt;font-family:Helvetica">MatRARtSym 25 1.0 1.2492e+02 1.0 0.00e+00 0.0 0.0e+00 0.0e+00 0.0e+00 1 0 0 0 0 1 0 0 0 0 0</span></p>
</div>
</div>
</div>
<div>
<div>
<div>
<p class="MsoNormal"><span style="font-size:9pt;font-family:Helvetica">MatMatTrnMultSym 25 1.0 8.8265e+01 1.0 0.00e+00 0.0 0.0e+00 0.0e+00 0.0e+00 0 0 0 0 0 0 0 0 0 0 0</span></p>
</div>
</div>
</div>
<div>
<div>
<div>
<p class="MsoNormal"><span style="font-size:9pt;font-family:Helvetica">MatMatTrnMultNum 25 1.0 2.4820e+02 1.0 6.83e+10 1.0 0.0e+00 0.0e+00 0.0e+00 1 0 0 0 0 1 0 0 0 0 275</span></p>
</div>
</div>
</div>
<div>
<div>
<div>
<p class="MsoNormal"><span style="font-size:9pt;font-family:Helvetica">MatTrnMatMultSym 10 1.0 7.2984e-01 1.0 0.00e+00 0.0 0.0e+00 0.0e+00 0.0e+00 0 0 0 0 0 0 0 0 0 0 0</span></p>
</div>
</div>
</div>
<div>
<div>
<div>
<p class="MsoNormal"><span style="font-size:9pt;font-family:Helvetica">MatTrnMatMultNum 10 1.0 9.3128e-01 1.0 0.00e+00 0.0 0.0e+00 0.0e+00 0.0e+00 0 0 0 0 0 0 0 0 0 0 0</span></p>
</div>
</div>
</div>
<div>
<div>
<div>
<p class="MsoNormal"><span style="font-size:9pt;font-family:Helvetica"> </span></p>
</div>
</div>
</div>
<div>
<div>
<div>
<p class="MsoNormal"><span style="font-size:9pt;font-family:Helvetica">in addition there are many more VecMAXPY then VecMDot (in GMRES they are each done the same number of times)</span></p>
</div>
</div>
</div>
<div>
<div>
<div>
<p class="MsoNormal"><span style="font-size:9pt;font-family:Helvetica"> </span></p>
</div>
</div>
</div>
<div>
<div>
<div>
<div>
<p class="MsoNormal"><span style="font-size:9pt;font-family:Helvetica">VecMDot 5588 1.0 1.7183e+03 1.0 2.06e+13 1.0 0.0e+00 0.0e+00 0.0e+00 8 10 0 0 0 8 10 0 0 0 12016</span></p>
</div>
</div>
</div>
</div>
<div>
<div>
<div>
<div>
<p class="MsoNormal"><span style="font-size:9pt;font-family:Helvetica">VecMAXPY 22412 1.0 8.4898e+03 1.0 4.17e+13 1.0 0.0e+00 0.0e+00 0.0e+00 39 20 0 0 0 39 20 0 0 0 4913</span></p>
</div>
</div>
</div>
</div>
<div>
<div>
<div>
<p class="MsoNormal"><span style="font-size:9pt;font-family:Helvetica"> </span></p>
</div>
</div>
</div>
<div>
<div>
<div>
<p class="MsoNormal"><span style="font-size:9pt;font-family:Helvetica">Finally there are a huge number of </span></p>
</div>
</div>
</div>
<div>
<div>
<div>
<p class="MsoNormal"><span style="font-size:9pt;font-family:Helvetica"> </span></p>
</div>
</div>
</div>
<div>
<div>
<div>
<div>
<p class="MsoNormal"><span style="font-size:9pt;font-family:Helvetica">MatMultAdd 258048 1.0 1.4178e+03 1.0 6.10e+13 1.0 0.0e+00 0.0e+00 0.0e+00 7 29 0 0 0 7 29 0 0 0 43025</span></p>
</div>
</div>
</div>
</div>
<div>
<div>
<div>
<p class="MsoNormal"><span style="font-size:9pt;font-family:Helvetica"> </span></p>
</div>
</div>
</div>
<div>
<div>
<div>
<p class="MsoNormal"><span style="font-size:9pt;font-family:Helvetica">Are you making calls to all these routines? Are you doing this inside your MatMult() or before you call KSPSolve?</span></p>
</div>
</div>
</div>
<div>
<div>
<div>
<p class="MsoNormal"><span style="font-size:9pt;font-family:Helvetica"> </span></p>
</div>
</div>
</div>
<div>
<div>
<div>
<p class="MsoNormal"><span style="font-size:9pt;font-family:Helvetica">The reason I wanted you to make a simpler run without the initial guess code is that your events are far more complicated than would be produced by GMRES alone so it is not possible to
understand the behavior you are seeing without fully understanding all the events happening in the code.</span></p>
</div>
</div>
</div>
<div>
<div>
<div>
<p class="MsoNormal"><span style="font-size:9pt;font-family:Helvetica"> </span></p>
</div>
</div>
</div>
<div>
<div>
<div>
<p class="MsoNormal"><span style="font-size:9pt;font-family:Helvetica"> <span> </span>Barry</span></p>
</div>
</div>
</div>
<div>
<div>
<div>
<p class="MsoNormal"><span style="font-size:9pt;font-family:Helvetica"> </span></p>
</div>
</div>
</div>
<div>
<p class="MsoNormal" style="margin-bottom:12pt"><span style="font-size:9pt;font-family:Helvetica"> </span></p>
<blockquote style="margin-top:5pt;margin-bottom:5pt">
<div>
<div>
<div>
<p class="MsoNormal"><span style="font-size:9pt;font-family:Helvetica">On Jun 14, 2024, at 1:19</span><span style="font-size:9pt;font-family:Arial,sans-serif"> </span><span style="font-size:9pt;font-family:Helvetica">AM, Yongzhong Li <</span><a href="mailto:yongzhong.li@mail.utoronto.ca" target="_blank"><span style="font-size:9pt;font-family:Helvetica">yongzhong.li@mail.utoronto.ca</span></a><span style="font-size:9pt;font-family:Helvetica">>
wrote:</span><u></u><u></u></p>
</div>
</div>
</div>
<div>
<div>
<p class="MsoNormal"><span style="font-size:9pt;font-family:Helvetica"> </span></p>
</div>
</div>
<div>
<div>
<div>
<div>
<div>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11pt;font-family:Helvetica">Thanks, I have attached the results without using any KSPGuess. At low frequency, the iteration steps are quite close to the one with KSPGuess, specifically<span> </span><br>
<br>
KSPGuess Object: 1 MPI process</span></p>
</div>
</div>
</div>
<div>
<div>
<div>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11pt;font-family:Helvetica"> <span> </span>type: fischer</span></p>
</div>
</div>
</div>
<div>
<div>
<div>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11pt;font-family:Helvetica"> <span> </span>Model 1, size 200<br>
<br>
However, I found at higher frequency, the # of iteration steps are significant higher than the one with KSPGuess, I have attahced both of the results for your reference.<br>
<br>
Moreover, could I ask why the one without the KSPGuess options can be used for a baseline comparsion? What are we comparing here? How does it relate to the performance issue/bottleneck I found? “</span><b><span style="font-size:9pt;font-family:Helvetica">I
have noticed that the time taken by </span></b><b><span style="font-size:9pt;font-family:"Courier New"">KSPSolve</span></b><b><span style="font-size:9pt;font-family:Helvetica"> is<span> </span></span></b><b><span lang="EN-US" style="font-size:9pt;font-family:Helvetica">almost
two times<span> </span></span></b><b><span style="font-size:9pt;font-family:Helvetica">greater than the CPU time for matrix-vector product multiplied by the number
of iteration</span></b><span lang="EN-US" style="font-size:11pt;font-family:Helvetica">”<span> </span><br>
<br>
Thank you!<br>
Yongzhong</span></p>
</div>
</div>
</div>
<div>
<div>
<div>
<div>
<p class="MsoNormal"><span style="font-size:11pt;font-family:Helvetica"> </span></p>
</div>
</div>
</div>
</div>
<div id="m_-7305460844614688919m_2643581689528628203m_6159968039748304261m_1825240119000957322m_-6075196701903926580m_-4861803116733450987mail-editor-reference-message-container">
<div>
<div>
<p class="MsoNormal" style="margin-bottom:12pt"><b><span style="font-size:9pt;font-family:Helvetica">From:<span> </span></span></b><span style="font-size:9pt;font-family:Helvetica">Barry
Smith <</span><a href="mailto:bsmith@petsc.dev" target="_blank"><span style="font-size:9pt;font-family:Helvetica">bsmith@petsc.dev</span></a><span style="font-size:9pt;font-family:Helvetica">><br>
<b>Date:<span> </span></b>Thursday, June 13, 2024 at 2:14</span><span style="font-size:9pt;font-family:Arial,sans-serif"> </span><span style="font-size:9pt;font-family:Helvetica">PM<br>
<b>To:<span> </span></b>Yongzhong Li <</span><a href="mailto:yongzhong.li@mail.utoronto.ca" target="_blank"><span style="font-size:9pt;font-family:Helvetica">yongzhong.li@mail.utoronto.ca</span></a><span style="font-size:9pt;font-family:Helvetica">><br>
<b>Cc:<span> </span></b></span><a href="mailto:petsc-users@mcs.anl.gov" target="_blank"><span style="font-size:9pt;font-family:Helvetica">petsc-users@mcs.anl.gov</span></a><span><span style="font-size:9pt;font-family:Helvetica"> </span></span><span style="font-size:9pt;font-family:Helvetica"><</span><a href="mailto:petsc-users@mcs.anl.gov" target="_blank"><span style="font-size:9pt;font-family:Helvetica">petsc-users@mcs.anl.gov</span></a><span style="font-size:9pt;font-family:Helvetica">>,<span> </span></span><a href="mailto:petsc-maint@mcs.anl.gov" target="_blank"><span style="font-size:9pt;font-family:Helvetica">petsc-maint@mcs.anl.gov</span></a><span><span style="font-size:9pt;font-family:Helvetica"> </span></span><span style="font-size:9pt;font-family:Helvetica"><</span><a href="mailto:petsc-maint@mcs.anl.gov" target="_blank"><span style="font-size:9pt;font-family:Helvetica">petsc-maint@mcs.anl.gov</span></a><span style="font-size:9pt;font-family:Helvetica">>,
Piero Triverio <</span><a href="mailto:piero.triverio@utoronto.ca" target="_blank"><span style="font-size:9pt;font-family:Helvetica">piero.triverio@utoronto.ca</span></a><span style="font-size:9pt;font-family:Helvetica">><br>
<b>Subject:<span> </span></b>Re: [petsc-maint] Assistance Needed with PETSc KSPSolve Performance Issue</span></p>
</div>
<div>
<div>
<div>
<div>
<p class="MsoNormal"><span style="font-size:9pt;font-family:Helvetica"> </span><u></u><u></u></p>
</div>
</div>
</div>
</div>
<div>
<div>
<div>
<p class="MsoNormal"><span style="font-size:9pt;font-family:Helvetica"> <span> </span>Can you please run the same thing without the KSPGuess option(s) for a baseline comparison?</span></p>
</div>
</div>
</div>
<div>
<div>
<div>
<div>
<div>
<p class="MsoNormal"><span style="font-size:9pt;font-family:Helvetica"> </span></p>
</div>
</div>
</div>
</div>
</div>
<div>
<div>
<div>
<div>
<p class="MsoNormal"><span style="font-size:9pt;font-family:Helvetica"> Thanks</span></p>
</div>
</div>
</div>
</div>
<div>
<div>
<div>
<div>
<div>
<p class="MsoNormal"><span style="font-size:9pt;font-family:Helvetica"> </span></p>
</div>
</div>
</div>
</div>
</div>
<div>
<div>
<div>
<div>
<p class="MsoNormal"><span style="font-size:9pt;font-family:Helvetica"> Barry</span></p>
</div>
</div>
</div>
</div>
<div>
<div>
<div>
<p class="MsoNormal" style="margin-bottom:12pt"><span style="font-size:9pt;font-family:Helvetica"> </span></p>
</div>
<blockquote style="margin-top:5pt;margin-bottom:5pt">
<div>
<div>
<div>
<div>
<p class="MsoNormal"><span style="font-size:9pt;font-family:Helvetica">On Jun 13, 2024, at 1:27</span><span style="font-size:9pt;font-family:Arial,sans-serif"> </span><span style="font-size:9pt;font-family:Helvetica">PM, Yongzhong Li <</span><a href="mailto:yongzhong.li@mail.utoronto.ca" target="_blank"><span style="font-size:9pt;font-family:Helvetica">yongzhong.li@mail.utoronto.ca</span></a><span style="font-size:9pt;font-family:Helvetica">>
wrote:</span><u></u><u></u></p>
</div>
</div>
</div>
</div>
<div>
<div>
<div>
<div>
<p class="MsoNormal"><span style="font-size:9pt;font-family:Helvetica"> </span></p>
</div>
</div>
</div>
</div>
<div>
<div id="m_-7305460844614688919m_2643581689528628203m_6159968039748304261m_1825240119000957322m_-6075196701903926580m_-4861803116733450987pfptBannerxs2204y">
<div id="m_-7305460844614688919m_2643581689528628203m_6159968039748304261m_1825240119000957322m_-6075196701903926580m_-4861803116733450987pfptBannerxs2204y">
<div id="m_-7305460844614688919m_2643581689528628203m_6159968039748304261m_1825240119000957322m_-6075196701903926580m_-4861803116733450987pfptBannerxs2204y">
<div>
<div>
<div>
<p class="MsoNormal"><span style="font-size:9pt;font-family:Arial,sans-serif">This Message Is From an External Sender</span></p>
</div>
</div>
</div>
</div>
<div id="m_-7305460844614688919m_2643581689528628203m_6159968039748304261m_1825240119000957322m_-6075196701903926580m_-4861803116733450987pfptBannerxs2204y">
<div>
<div>
<div>
<p class="MsoNormal"><span style="font-size:9pt;font-family:Arial,sans-serif">This message came from outside your organization.</span></p>
</div>
</div>
</div>
</div>
</div>
</div>
<div>
<div>
<div>
<div>
<div>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11pt;font-family:Helvetica">Hi Matt,<br>
<br>
I have rerun the program with the keys you provided. The system output when performing ksp solve and the final petsc log output were stored in a .txt file attached for your reference.<br>
<br>
Thanks!<br>
Yongzhong</span></p>
</div>
</div>
</div>
</div>
<div>
<div>
<div>
<div>
<div>
<p class="MsoNormal"><span style="font-size:11pt;font-family:Helvetica"> </span></p>
</div>
</div>
</div>
</div>
</div>
<div id="m_-7305460844614688919m_2643581689528628203m_6159968039748304261m_1825240119000957322m_-6075196701903926580m_-4861803116733450987mail-editor-reference-message-container">
<div>
<div>
<p class="MsoNormal" style="margin-bottom:12pt"><b><span style="font-size:9pt;font-family:Helvetica">From:<span> </span></span></b><span style="font-size:9pt;font-family:Helvetica">Matthew
Knepley <</span><a href="mailto:knepley@gmail.com" target="_blank"><span style="font-size:9pt;font-family:Helvetica">knepley@gmail.com</span></a><span style="font-size:9pt;font-family:Helvetica">><br>
<b>Date:<span> </span></b>Wednesday, June 12, 2024 at 6:46</span><span style="font-size:9pt;font-family:Arial,sans-serif"> </span><span style="font-size:9pt;font-family:Helvetica">PM<br>
<b>To:<span> </span></b>Yongzhong Li <</span><a href="mailto:yongzhong.li@mail.utoronto.ca" target="_blank"><span style="font-size:9pt;font-family:Helvetica">yongzhong.li@mail.utoronto.ca</span></a><span style="font-size:9pt;font-family:Helvetica">><br>
<b>Cc:<span> </span></b></span><a href="mailto:petsc-users@mcs.anl.gov" target="_blank"><span style="font-size:9pt;font-family:Helvetica">petsc-users@mcs.anl.gov</span></a><span><span style="font-size:9pt;font-family:Helvetica"> </span></span><span style="font-size:9pt;font-family:Helvetica"><</span><a href="mailto:petsc-users@mcs.anl.gov" target="_blank"><span style="font-size:9pt;font-family:Helvetica">petsc-users@mcs.anl.gov</span></a><span style="font-size:9pt;font-family:Helvetica">>,<span> </span></span><a href="mailto:petsc-maint@mcs.anl.gov" target="_blank"><span style="font-size:9pt;font-family:Helvetica">petsc-maint@mcs.anl.gov</span></a><span><span style="font-size:9pt;font-family:Helvetica"> </span></span><span style="font-size:9pt;font-family:Helvetica"><</span><a href="mailto:petsc-maint@mcs.anl.gov" target="_blank"><span style="font-size:9pt;font-family:Helvetica">petsc-maint@mcs.anl.gov</span></a><span style="font-size:9pt;font-family:Helvetica">>,
Piero Triverio <</span><a href="mailto:piero.triverio@utoronto.ca" target="_blank"><span style="font-size:9pt;font-family:Helvetica">piero.triverio@utoronto.ca</span></a><span style="font-size:9pt;font-family:Helvetica">><br>
<b>Subject:<span> </span></b>Re: [petsc-maint] Assistance Needed with PETSc KSPSolve Performance Issue</span><u></u><u></u></p>
</div>
<table border="0" cellspacing="0" cellpadding="0" align="left" width="100%" style="width:100%;display:table;float:none">
<tbody>
<tr>
<td style="background:rgb(166,166,166);padding:5.25pt 1.5pt"></td>
<td width="100%" style="width:100%;background:rgb(234,234,234);padding:5.25pt 3.75pt 5.25pt 11.25pt">
<div>
<div>
<div>
<div>
<div>
<p class="MsoNormal">
<span lang="ZH-CN" style="font-size:9pt;font-family:DengXian;color:rgb(33,33,33)">你通常不会收到来自</span><span><span style="font-size:9pt;font-family:"Segoe UI",sans-serif;color:rgb(33,33,33)"> </span></span><span style="color:black"><a href="mailto:knepley@gmail.com" target="_blank"><span style="font-size:9pt;font-family:"Segoe UI",sans-serif">knepley@gmail.com</span></a></span><span><span style="font-size:9pt;font-family:"Segoe UI",sans-serif;color:rgb(33,33,33)"> </span></span><span lang="ZH-CN" style="font-size:9pt;font-family:DengXian;color:rgb(33,33,33)">的电子邮件。</span><span style="color:black"><a href="https://urldefense.us/v3/__https://aka.ms/LearnAboutSenderIdentification__;!!G_uCfscf7eWS!djGfJnEhNJROfsMsBJy5u_KoRKbug55xZ64oHKUFnH2cWku_Th1hwt4TDdoMd8pWYVDzJeqJslMNZwpO3y0Et94d31qk-oCEwo4$" target="_blank"><span lang="ZH-CN" style="font-size:9pt;font-family:DengXian">了解这一点为什么很重要</span></a></span></p>
</div>
</div>
</div>
</div>
</div>
</td>
<td width="75" style="width:56.25pt;background:rgb(234,234,234);padding:5.25pt 3.75pt">
</td>
</tr>
</tbody>
</table>
<div>
<div>
<div>
<div>
<div>
<div>
<div>
<p class="MsoNormal"><span style="font-size:9pt;font-family:Helvetica">On Wed, Jun 12, 2024 at 6:36</span><span style="font-size:9pt;font-family:Arial,sans-serif"> </span><span style="font-size:9pt;font-family:Helvetica">PM Yongzhong Li <</span><a href="mailto:yongzhong.li@mail.utoronto.ca" target="_blank"><span style="font-size:9pt;font-family:Helvetica">yongzhong.li@mail.utoronto.ca</span></a><span style="font-size:9pt;font-family:Helvetica">>
wrote:</span><u></u><u></u></p>
</div>
</div>
</div>
</div>
</div>
<div>
<blockquote style="border-top:none currentcolor;border-right:none currentcolor;border-bottom:none currentcolor;border-left:1pt solid rgb(204,204,204);padding:0cm 0cm 0cm 6pt;margin:5pt 0cm 5pt 4.8pt">
<div>
<div>
<div>
<div>
<div>
<div>
<p class="MsoNormal"><span style="font-size:1pt;font-family:Helvetica;color:white">Dear PETSc’s developers, I hope this email finds you well. I am currently working on a project using PETSc and have encountered a performance issue with the KSPSolve function.
Specifically, I have noticed that the time taken by KSPSolve is<span> </span></span></p>
</div>
</div>
</div>
</div>
</div>
<div>
<div>
<div>
<div>
<div>
<p class="MsoNormal"><span style="font-size:1pt;font-family:Helvetica;color:white">ZjQcmQRYFpfptBannerStart</span></p>
</div>
</div>
</div>
</div>
</div>
<div style="margin-top:12pt;margin-bottom:12pt;border-width:initial;border-style:initial;border-color:currentcolor;border-radius:4px;min-width:200px" id="m_-7305460844614688919m_2643581689528628203m_6159968039748304261m_1825240119000957322m_-6075196701903926580m_-4861803116733450987m_5223539340070699115pfptBannerpeukvww">
<div style="margin-bottom:0.75pt;float:left" id="m_-7305460844614688919m_2643581689528628203m_6159968039748304261m_1825240119000957322m_-6075196701903926580m_-4861803116733450987m_5223539340070699115pfptBannerpeukvww">
<div id="m_-7305460844614688919m_2643581689528628203m_6159968039748304261m_1825240119000957322m_-6075196701903926580m_-4861803116733450987m_5223539340070699115pfptBannerpeukvww">
<div>
<div>
<div>
<div>
<p class="MsoNormal" style="line-height:13.5pt;background:rgb(208,216,220)"><b><span style="font-size:10.5pt;font-family:Arial,sans-serif;color:black">This Message Is From an External Sender</span></b></p>
</div>
</div>
</div>
</div>
</div>
<div style="margin-top:1.5pt" id="m_-7305460844614688919m_2643581689528628203m_6159968039748304261m_1825240119000957322m_-6075196701903926580m_-4861803116733450987m_5223539340070699115pfptBannerpeukvww">
<div>
<div>
<div>
<div>
<p class="MsoNormal" style="line-height:13.5pt;background:rgb(208,216,220)"><span style="font-size:9pt;font-family:Arial,sans-serif;color:black">This message came from outside your organization.</span></p>
</div>
</div>
</div>
</div>
</div>
</div>
<div>
<div>
<div>
<div>
<div>
<div>
<p class="MsoNormal" style="background:rgb(208,216,220)"><span style="font-size:1pt;font-family:Helvetica;color:black"> </span></p>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
<div>
<div>
<div>
<div>
<div>
<p class="MsoNormal"><span style="font-size:1pt;font-family:Helvetica;color:white">ZjQcmQRYFpfptBannerEnd</span></p>
</div>
</div>
</div>
</div>
</div>
<div>
<div>
<div>
<div>
<div>
<div>
<p class="MsoNormal"><span style="font-size:9pt;font-family:Helvetica">Dear<span> </span></span><span lang="EN-US" style="font-size:9pt;font-family:Helvetica">PETSc’s
developers</span><span style="font-size:9pt;font-family:Helvetica">,</span></p>
</div>
</div>
</div>
</div>
<div>
<div>
<div>
<div>
<p class="MsoNormal"><span style="font-size:9pt;font-family:Helvetica">I hope this email finds you well.</span></p>
</div>
</div>
</div>
</div>
<div>
<div>
<div>
<div>
<p class="MsoNormal"><span style="font-size:9pt;font-family:Helvetica">I am currently working on a project using PETSc and have encountered a performance issue with the </span><span style="font-size:9pt;font-family:"Courier New"">KSPSolve</span><span style="font-size:9pt;font-family:Helvetica"> function.
Specifically,<span> </span><b>I have noticed that the time taken by </b></span><b><span style="font-size:9pt;font-family:"Courier New"">KSPSolve</span></b><b><span style="font-size:9pt;font-family:Helvetica"> is<span> </span></span></b><b><span lang="EN-US" style="font-size:9pt;font-family:Helvetica">almost
two times<span> </span></span></b><b><span style="font-size:9pt;font-family:Helvetica">greater than the CPU time for matrix-vector product multiplied by the number
of iteration steps</span></b><span style="font-size:9pt;font-family:Helvetica">.<span> </span></span><span lang="EN-US" style="font-size:9pt;font-family:Helvetica">I
use C++ chrono to record CPU time.</span></p>
</div>
</div>
</div>
</div>
<div>
<div>
<div>
<div>
<p class="MsoNormal"><span style="font-size:9pt;font-family:Helvetica">For context, I am using a shell system matrix</span><span><span lang="EN-US" style="font-size:9pt;font-family:Helvetica"> </span></span><span lang="EN-US" style="font-size:9pt;font-family:Helvetica">A</span><span style="font-size:9pt;font-family:Helvetica">.
Despite my efforts to<span> </span></span><span lang="EN-US" style="font-size:9pt;font-family:Helvetica">parallelize</span><span><span style="font-size:9pt;font-family:Helvetica"> </span></span><span style="font-size:9pt;font-family:Helvetica">the
matrix-vector product</span><span><span lang="EN-US" style="font-size:9pt;font-family:Helvetica"> </span></span><span lang="EN-US" style="font-size:9pt;font-family:Helvetica">(Ax)</span><span style="font-size:9pt;font-family:Helvetica">,
the overall solve time remains higher</span><span><span lang="EN-US" style="font-size:9pt;font-family:Helvetica"> </span></span><span lang="EN-US" style="font-size:9pt;font-family:Helvetica">than
the matrix vector product per iteration indicates when multiple threads were used</span><span style="font-size:9pt;font-family:Helvetica">. Here are a few details of my setup:</span></p>
</div>
</div>
</div>
</div>
<ul style="margin-top:0cm" type="disc">
<li class="MsoNormal"><b><span style="font-size:9pt;font-family:Helvetica">Matrix Type</span></b><span style="font-size:9pt;font-family:Helvetica">: Shell system matrix</span><u></u><u></u></li><li class="MsoNormal"><b><span style="font-size:9pt;font-family:Helvetica">Preconditioner</span></b><span style="font-size:9pt;font-family:Helvetica">:<span> </span></span><span lang="EN-US" style="font-size:9pt;font-family:Helvetica">Shell
PC</span><u></u><u></u></li><li class="MsoNormal"><b><span style="font-size:9pt;font-family:Helvetica">Parallel Environment</span></b><span style="font-size:9pt;font-family:Helvetica">:<span> </span></span><span lang="EN-US" style="font-size:9pt;font-family:Helvetica">Using
Intel MKL as PETSc’s BLAS/LAPACK library, multithreading is enabled</span><u></u><u></u></li></ul>
<div>
<div>
<div>
<div>
<p class="MsoNormal"><span lang="EN-US" style="font-size:9pt;font-family:Helvetica">I have considered several potential reasons, such as preconditioner setup, additional solver operations, and the inherent overhead of using a shell system matrix.<span> </span><b>However,
since KSPSolve is a high-level API, I have been unable to pinpoint the exact cause of the increased solve time.</b></span><u></u><u></u></p>
</div>
</div>
</div>
</div>
<div>
<div>
<div>
<div>
<p class="MsoNormal"><span lang="EN-US" style="font-size:9pt;font-family:Helvetica">Have you observed the same issue?<span> </span></span><span style="font-size:9pt;font-family:Helvetica">Could
you</span><span><span lang="EN-US" style="font-size:9pt;font-family:Helvetica"> </span></span><span lang="EN-US" style="font-size:9pt;font-family:Helvetica">please<span> </span></span><span style="font-size:9pt;font-family:Helvetica">provide
some<span> </span></span><span lang="EN-US" style="font-size:9pt;font-family:Helvetica">experience<span> </span></span><span style="font-size:9pt;font-family:Helvetica">on
how to diagnose and address this performance discrepancy? Any insights or recommendations you could offer would be greatly appreciated.</span><u></u><u></u></p>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</blockquote>
<div>
<div>
<div>
<div>
<div>
<p class="MsoNormal"><span style="font-size:9pt;font-family:Helvetica"> </span><u></u><u></u></p>
</div>
</div>
</div>
</div>
</div>
<div>
<div>
<div>
<div>
<div>
<p class="MsoNormal"><span style="font-size:9pt;font-family:Helvetica">For any performance question like this, we need to see the output of your code run with</span><u></u><u></u></p>
</div>
</div>
</div>
</div>
</div>
<div>
<div>
<div>
<div>
<div>
<p class="MsoNormal"><span style="font-size:9pt;font-family:Helvetica"> </span><u></u><u></u></p>
</div>
</div>
</div>
</div>
</div>
<div>
<div>
<div>
<div>
<div>
<p class="MsoNormal"><span style="font-size:9pt;font-family:Helvetica"> <span> </span>-ksp_view -ksp_monitor_true_residual -ksp_converged_reason -log_view</span><u></u><u></u></p>
</div>
</div>
</div>
</div>
</div>
<div>
<div>
<div>
<div>
<div>
<p class="MsoNormal"><span style="font-size:9pt;font-family:Helvetica"> </span><u></u><u></u></p>
</div>
</div>
</div>
</div>
</div>
<div>
<div>
<div>
<div>
<div>
<p class="MsoNormal"><span style="font-size:9pt;font-family:Helvetica"> <span> </span>Thanks,</span><u></u><u></u></p>
</div>
</div>
</div>
</div>
</div>
<div>
<div>
<div>
<div>
<div>
<p class="MsoNormal"><span style="font-size:9pt;font-family:Helvetica"> </span><u></u><u></u></p>
</div>
</div>
</div>
</div>
</div>
<div>
<div>
<div>
<div>
<div>
<p class="MsoNormal"><span style="font-size:9pt;font-family:Helvetica"> Matt</span><u></u><u></u></p>
</div>
</div>
</div>
</div>
</div>
<div>
<div>
<div>
<div>
<div>
<p class="MsoNormal"><span style="font-size:9pt;font-family:Helvetica"> </span><u></u><u></u></p>
</div>
</div>
</div>
</div>
</div>
<blockquote style="border-top:none currentcolor;border-right:none currentcolor;border-bottom:none currentcolor;border-left:1pt solid rgb(204,204,204);padding:0cm 0cm 0cm 6pt;margin:5pt 0cm 5pt 4.8pt">
<div>
<div>
<div>
<div>
<div>
<div>
<div>
<p class="MsoNormal"><span style="font-size:9pt;font-family:Helvetica">Thank you for your time and assistance.</span><u></u><u></u></p>
</div>
</div>
</div>
</div>
<div>
<div>
<div>
<div>
<p class="MsoNormal"><span style="font-size:9pt;font-family:Helvetica">Best</span><span><span lang="EN-US" style="font-size:9pt;font-family:Helvetica"> </span></span><span lang="EN-US" style="font-size:9pt;font-family:Helvetica">regards,</span><u></u><u></u></p>
</div>
</div>
</div>
</div>
<div>
<div>
<div>
<div>
<p class="MsoNormal"><span lang="EN-US" style="font-size:9pt;font-family:Helvetica">Yongzhong</span><u></u><u></u></p>
</div>
</div>
</div>
</div>
<div>
<div>
<div>
<div>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11pt;font-family:Calibri,sans-serif">-----------------------------------------------------------</span><u></u><u></u></p>
</div>
</div>
</div>
</div>
<div>
<div>
<div>
<div>
<p class="MsoNormal"><b><span lang="EN-US" style="font-size:11pt;font-family:Calibri,sans-serif">Yongzhong Li</span></b><u></u><u></u></p>
</div>
</div>
</div>
</div>
<div>
<div>
<div>
<div>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11pt;font-family:Calibri,sans-serif">PhD student | Electromagnetics Group</span><u></u><u></u></p>
</div>
</div>
</div>
</div>
<div>
<div>
<div>
<div>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11pt;font-family:Calibri,sans-serif">Department of Electrical & Computer Engineering</span><u></u><u></u></p>
</div>
</div>
</div>
</div>
<div>
<div>
<div>
<div>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11pt;font-family:Calibri,sans-serif">University of Toronto</span><u></u><u></u></p>
</div>
</div>
</div>
</div>
<div>
<div>
<div>
<div>
<p class="MsoNormal"><a href="https://urldefense.us/v3/__http://www.modelics.org__;!!G_uCfscf7eWS!cuLttMJEcegaqu461Bt4QLsO4fASfLM5vjRbtyNhWJQiInbjgNwkGNdkFE1ebSbFjOUatYB0-jd2yQWMWzqkDFFjwMvNl3ZKAr8$" target="_blank"><span lang="EN-US" style="font-size:11pt;font-family:Calibri,sans-serif;color:rgb(5,99,193)">http://www.modelics.org</span></a><u></u><u></u></p>
</div>
</div>
</div>
</div>
<div>
<div>
<div>
<div>
<p class="MsoNormal"><span style="font-size:9pt;font-family:Helvetica"> </span><u></u><u></u></p>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</blockquote>
</div>
<div>
<div>
<div>
<div>
<p class="MsoNormal"><span style="font-size:9pt;font-family:Helvetica"><br clear="all">
</span><u></u><u></u></p>
</div>
</div>
</div>
</div>
<div>
<div>
<div>
<div>
<div>
<p class="MsoNormal"><span style="font-size:9pt;font-family:Helvetica"> </span><u></u><u></u></p>
</div>
</div>
</div>
</div>
</div>
<div>
<div>
<div>
<div>
<p class="MsoNormal"><span><span style="font-size:9pt;font-family:Helvetica">--</span></span><span><span style="font-size:9pt;font-family:Helvetica"> </span></span><u></u><u></u></p>
</div>
</div>
</div>
</div>
<div>
<div>
<div>
<div>
<div>
<div>
<div>
<div>
<div>
<div>
<div>
<p class="MsoNormal"><span style="font-size:9pt;font-family:Helvetica">What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.<br>
-- Norbert Wiener</span><u></u><u></u></p>
</div>
</div>
</div>
</div>
</div>
<div>
<div>
<div>
<div>
<div>
<p class="MsoNormal"><span style="font-size:9pt;font-family:Helvetica"> </span><u></u><u></u></p>
</div>
</div>
</div>
</div>
</div>
<div>
<div>
<div>
<div>
<div>
<p class="MsoNormal"><a href="https://urldefense.us/v3/__http://www.cse.buffalo.edu/*knepley/__;fg!!G_uCfscf7eWS!djGfJnEhNJROfsMsBJy5u_KoRKbug55xZ64oHKUFnH2cWku_Th1hwt4TDdoMd8pWYVDzJeqJslMNZwpO3y0Et94d31qkNOuenGA$" target="_blank"><span style="font-size:9pt;font-family:Helvetica">https://www.cse.buffalo.edu/~knepley/</span></a><u></u><u></u></p>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
<div>
<div>
<div>
<p class="MsoNormal"><span style="font-size:9pt;font-family:Helvetica"><ksp_petsc_log.txt></span><u></u><u></u></p>
</div>
</div>
</div>
</div>
</blockquote>
</div>
<div>
<div>
<div>
<p class="MsoNormal"><span style="font-size:9pt;font-family:Helvetica"> </span><u></u><u></u></p>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
<div>
<div>
<p class="MsoNormal"><span style="font-size:9pt;font-family:Helvetica"><ksp_petsc_log.txt><ksp_petsc_log_noguess.txt></span><u></u><u></u></p>
</div>
</div>
</div>
</blockquote>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</blockquote>
</div>
</div>
</blockquote>
</div>
</div>
</div>
</div>
</div>
</div>
</blockquote>
</div>
<div>
<p class="MsoNormal"> <u></u><u></u></p>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</blockquote>
</div>
</div>
</div>
</div>
</div>
</div>
</blockquote>
</div>
<div>
<p class="MsoNormal"><br clear="all">
<u></u><u></u></p>
</div>
<div>
<div>
<p class="MsoNormal"> <u></u><u></u></p>
</div>
</div>
<div>
<p class="MsoNormal"><span>--</span><span> </span><u></u><u></u></p>
</div>
<div>
<div>
<div>
<div>
<div>
<div>
<div>
<div>
<p class="MsoNormal">What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.<br>
-- Norbert Wiener<u></u><u></u></p>
</div>
</div>
<div>
<div>
<p class="MsoNormal"> <u></u><u></u></p>
</div>
</div>
<div>
<div>
<p class="MsoNormal"><a href="https://urldefense.us/v3/__http://www.cse.buffalo.edu/*knepley/__;fg!!G_uCfscf7eWS!fVvbGldqcUV5ju4jpu5oGmt-VjITi5JpCJzhHxpbgsERLVYZzglpxKOOyrBRGxjRxp7vWHwt3SnINFOQErR1Z8kcDcf3cNeD9Gw$" target="_blank">https://www.cse.buffalo.edu/~knepley/</a><u></u><u></u></p>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</blockquote>
</div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
</div>
</div>
</div>
</div></blockquote></div><br clear="all"><div><br></div><span class="gmail_signature_prefix">-- </span><br><div dir="ltr" class="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div>What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.<br>-- Norbert Wiener</div><div><br></div><div><a href="https://urldefense.us/v3/__http://www.cse.buffalo.edu/*knepley/__;fg!!G_uCfscf7eWS!aQJpmm5W6l6FUiumnIPmkouzwzNUfx-Dyq04i1O2KS_InQGk6qjI7wUir0Hx6QEUQE2AMiJDsez3x2Os2C2d$" target="_blank">https://www.cse.buffalo.edu/~knepley/</a><br></div></div></div></div></div></div></div></div>
</div></blockquote></div>
</blockquote></div>