<!-- 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;">
OFFICIAL USE ONLY / À USAGE EXCLUSIF I’d like to unsubscribe to this mailing list. Thanks. ~ Sharon From: petsc-users <petsc-users-bounces@ mcs. anl. gov> On Behalf Of Barry Smith Sent: Thursday, March 28, 2024 3: 26 PM To: Zou, Ling <lzou@ anl. gov>
</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="pfptBanner87md2uj" 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="pfptBanner87md2uj" style="all: unset !important; float:left !important; display:block !important; margin: 0px 0px 1px 0px !important; max-width: 600px !important;">
<div id="pfptBanner87md2uj" 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="pfptBanner87md2uj" 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 -->
<html xmlns:v="urn:schemas-microsoft-com:vml" xmlns:o="urn:schemas-microsoft-com:office:office" xmlns:w="urn:schemas-microsoft-com:office:word" xmlns:m="http://schemas.microsoft.com/office/2004/12/omml" xmlns="http://www.w3.org/TR/REC-html40">
<head><!-- BaNnErBlUrFlE-HeAdEr-start -->
<style>
#pfptBanner87md2uj { all: revert !important; display: block !important;
visibility: visible !important; opacity: 1 !important;
background-color: #D0D8DC !important;
max-width: none !important; max-height: none !important }
.pfptPrimaryButton87md2uj:hover, .pfptPrimaryButton87md2uj:focus {
background-color: #b4c1c7 !important; }
.pfptPrimaryButton87md2uj:active {
background-color: #90a4ae !important; }
</style>
<!-- BaNnErBlUrFlE-HeAdEr-end -->
<meta http-equiv="Content-Type" content="text/html; charset=utf-8">
<meta name="Generator" content="Microsoft Word 15 (filtered medium)">
<!--[if !mso]><style>v\:* {behavior:url(#default#VML);}
o\:* {behavior:url(#default#VML);}
w\:* {behavior:url(#default#VML);}
.shape {behavior:url(#default#VML);}
</style><![endif]--><style><!--
/* Font Definitions */
@font-face
{font-family:Helvetica;
panose-1:2 11 6 4 2 2 2 2 2 4;}
@font-face
{font-family:"Cambria Math";
panose-1:2 4 5 3 5 4 6 3 2 4;}
@font-face
{font-family:Calibri;
panose-1:2 15 5 2 2 2 4 3 2 4;}
@font-face
{font-family:Aptos;}
@font-face
{font-family:"Andale Mono";}
@font-face
{font-family:Garamond;
panose-1:2 2 4 4 3 3 1 1 8 3;}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
{margin:0cm;
font-size:10.0pt;
font-family:"Calibri",sans-serif;}
a:link, span.MsoHyperlink
{mso-style-priority:99;
color:blue;
text-decoration:underline;}
span.apple-converted-space
{mso-style-name:apple-converted-space;}
span.s1
{mso-style-name:s1;}
span.EmailStyle20
{mso-style-type:personal-reply;
font-family:"Calibri",sans-serif;
color:blue;
font-weight:normal;
font-style:normal;
text-decoration:none none;}
p.msipheader5f328e30, li.msipheader5f328e30, div.msipheader5f328e30
{mso-style-name:msipheader5f328e30;
mso-margin-top-alt:auto;
margin-right:0cm;
mso-margin-bottom-alt:auto;
margin-left:0cm;
font-size:11.0pt;
font-family:"Calibri",sans-serif;}
.MsoChpDefault
{mso-style-type:export-only;
font-size:10.0pt;
mso-ligatures:none;}
@page WordSection1
{size:612.0pt 792.0pt;
margin:72.0pt 72.0pt 72.0pt 72.0pt;}
div.WordSection1
{page:WordSection1;}
--></style><!--[if gte mso 9]><xml>
<o:shapedefaults v:ext="edit" spidmax="1026" />
</xml><![endif]--><!--[if gte mso 9]><xml>
<o:shapelayout v:ext="edit">
<o:idmap v:ext="edit" data="1" />
</o:shapelayout></xml><![endif]-->
</head>
<body lang="EN-CA" link="blue" vlink="purple" style="word-wrap:break-word">
<div class="WordSection1">
<p class="msipheader5f328e30" align="right" style="margin:0cm;text-align:right"><span style="font-size:10.0pt;color:black">OFFICIAL USE ONLY / À USAGE EXCLUSIF</span><o:p></o:p></p>
<p class="MsoNormal"><span style="font-size:11.0pt"><o:p> </o:p></span></p>
<p class="MsoNormal"><span style="font-size:12.0pt;color:blue;mso-fareast-language:EN-US">I’d like to unsubscribe to this mailing list.<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:12.0pt;color:blue;mso-fareast-language:EN-US"><o:p> </o:p></span></p>
<p class="MsoNormal"><span style="font-size:12.0pt;color:blue;mso-fareast-language:EN-US">Thanks.<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:12.0pt;color:blue;mso-fareast-language:EN-US"><o:p> </o:p></span></p>
<p class="MsoNormal"><i><span style="font-size:12.0pt;font-family:"Garamond",serif;color:#660033">~ Sharon<o:p></o:p></span></i></p>
<p class="MsoNormal"><span style="font-size:12.0pt;color:blue;mso-fareast-language:EN-US"><o:p> </o:p></span></p>
<div style="border:none;border-top:solid #E1E1E1 1.0pt;padding:3.0pt 0cm 0cm 0cm">
<p class="MsoNormal"><b><span lang="EN-US" style="font-size:11.0pt">From:</span></b><span lang="EN-US" style="font-size:11.0pt"> petsc-users <petsc-users-bounces@mcs.anl.gov>
<b>On Behalf Of </b>Barry Smith<br>
<b>Sent:</b> Thursday, March 28, 2024 3:26 PM<br>
<b>To:</b> Zou, Ling <lzou@anl.gov><br>
<b>Cc:</b> petsc-users@mcs.anl.gov<br>
<b>Subject:</b> [External] Re: [petsc-users] Does ILU(15) still make sense or should just use LU?<o:p></o:p></span></p>
</div>
<p class="MsoNormal"><o:p> </o:p></p>
<div>
<p class="MsoNormal" style="mso-line-height-alt:.75pt"><span style="font-size:1.0pt;color:white">You may benefit from a literature search on your model AND preconditioners to see what others have used. But I would try PETSc/MUMPS on the biggest size you want
and see how it goes (better it runs for a little longer and you don't waste months
<o:p></o:p></span></p>
</div>
<div>
<p class="MsoNormal" style="mso-line-height-alt:.75pt"><span style="font-size:1.0pt;color:white">ZjQcmQRYFpfptBannerStart<o:p></o:p></span></p>
</div>
<table class="MsoNormalTable" border="0" cellspacing="0" cellpadding="0" width="100%" style="width:100.0%;border-radius:4px">
<tbody>
<tr>
<td style="padding:12.0pt 0cm 12.0pt 0cm">
<table class="MsoNormalTable" border="1" cellspacing="0" cellpadding="0" width="100%" style="width:100.0%;background:#D0D8DC;border:none;border-top:solid #90A4AE 3.0pt">
<tbody>
<tr>
<td valign="top" style="border:none;padding:0cm 7.5pt 3.75pt 4.5pt">
<table class="MsoNormalTable" border="0" cellspacing="0" cellpadding="0" align="left">
<tbody>
<tr>
<td style="padding:3.0pt 6.0pt 3.0pt 6.0pt">
<p class="MsoNormal"><b><span style="font-size:10.5pt;font-family:"Arial",sans-serif;color:black">This Message Is From an External Sender
<o:p></o:p></span></b></p>
</td>
</tr>
<tr>
<td style="padding:3.0pt 6.0pt 3.0pt 6.0pt">
<p class="MsoNormal"><span style="font-size:9.0pt;font-family:"Arial",sans-serif;color:black">This message came from outside your organization. Do not open attachments or click links from unknown senders or unexpected email
<o:p></o:p></span></p>
</td>
</tr>
</tbody>
</table>
<table class="MsoNormalTable" border="0" cellspacing="0" cellpadding="0" align="right">
<tbody>
<tr>
<td style="padding:3.0pt 0cm 3.0pt 0cm">
<p class="MsoNormal"> <a href="https://us-phishalarm-ewt.proofpoint.com/EWT/v1/I4HAk28DiqjYeRyS71U!1R4lyukC3Ed5L2o7iduV5poB6rBS6qdB-U1FW_TqIgjA_xzIUeaCsMG8LclAzVjZorV7mIx3TwCgK1k7a-ev1Mjdfhq6a4ReIg_cpA$" target="_blank"><strong><span style="font-size:10.5pt;font-family:"Arial",sans-serif;color:black;border:solid #666666 1.0pt;padding:6.0pt;font-weight:normal;text-decoration:none"> Report Suspicious </span></strong></a>
<span style="font-size:11.0pt"><o:p></o:p></span></p>
</td>
</tr>
</tbody>
</table>
</td>
</tr>
</tbody>
</table>
</td>
</tr>
</tbody>
</table>
<div>
<p class="MsoNormal" style="mso-line-height-alt:.75pt"><span style="font-size:1.0pt;color:white">ZjQcmQRYFpfptBannerEnd<o:p></o:p></span></p>
</div>
<div>
<p class="MsoNormal" style="mso-line-height-alt:.75pt"><span style="font-size:1.0pt;color:white">You may benefit from a literature search on your model AND preconditioners to see what others have used. But I would try PETSc/MUMPS on the biggest size you want
and see how it goes (better it runs for a little longer and you don't waste months
<o:p></o:p></span></p>
</div>
<div>
<p class="MsoNormal" style="mso-line-height-alt:.75pt"><span style="font-size:1.0pt;color:white">ZjQcmQRYFpfptBannerStart<o:p></o:p></span></p>
</div>
<table class="MsoNormalTable" border="0" cellspacing="0" cellpadding="0" width="100%" style="width:100.0%;border-radius:4px">
<tbody>
<tr>
<td style="padding:12.0pt 0cm 12.0pt 0cm">
<table class="MsoNormalTable" border="1" cellspacing="0" cellpadding="0" width="100%" style="width:100.0%;background:#D0D8DC;border:none;border-top:solid #90A4AE 3.0pt">
<tbody>
<tr>
<td valign="top" style="border:none;padding:0cm 7.5pt 3.75pt 4.5pt">
<table class="MsoNormalTable" border="0" cellspacing="0" cellpadding="0" align="left">
<tbody>
<tr>
<td style="padding:3.0pt 6.0pt 3.0pt 6.0pt">
<p class="MsoNormal"><b><span style="font-size:10.5pt;font-family:"Arial",sans-serif;color:black">This Message Is From an External Sender
<o:p></o:p></span></b></p>
</td>
</tr>
<tr>
<td style="padding:3.0pt 6.0pt 3.0pt 6.0pt">
<p class="MsoNormal"><span style="font-size:9.0pt;font-family:"Arial",sans-serif;color:black">This message came from outside your organization.
<o:p></o:p></span></p>
</td>
</tr>
</tbody>
</table>
</td>
</tr>
</tbody>
</table>
</td>
</tr>
</tbody>
</table>
<div>
<p class="MsoNormal" style="mso-line-height-alt:.75pt"><span style="font-size:1.0pt;color:white">ZjQcmQRYFpfptBannerEnd<o:p></o:p></span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-size:11.0pt"><o:p> </o:p></span></p>
</div>
<p class="MsoNormal"><span style="font-size:11.0pt"> You may benefit from a literature search on your model AND preconditioners to see what others have used. But I would try PETSc/MUMPS on the biggest size you want and see how it goes (better it runs for
a little longer and you don't waste months trying to find a good preconditioner).<o:p></o:p></span></p>
<div>
<p class="MsoNormal"><span style="font-size:11.0pt"><o:p> </o:p></span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-size:11.0pt"><o:p> </o:p></span></p>
<div>
<p class="MsoNormal"><span style="font-size:11.0pt"><o:p> </o:p></span></p>
<div>
<p class="MsoNormal" style="margin-bottom:12.0pt"><span style="font-size:11.0pt"><o:p> </o:p></span></p>
<blockquote style="margin-top:5.0pt;margin-bottom:5.0pt">
<div>
<p class="MsoNormal"><span style="font-size:11.0pt">On Mar 28, 2024, at 2:20 PM, Zou, Ling <<a href="mailto:lzou@anl.gov">lzou@anl.gov</a>> wrote:<o:p></o:p></span></p>
</div>
<p class="MsoNormal"><span style="font-size:11.0pt"><o:p> </o:p></span></p>
<div>
<div>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Aptos",sans-serif">Thank you, Barry.</span><span style="font-size:12.0pt;font-family:"Aptos",sans-serif"><o:p></o:p></span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Aptos",sans-serif">Yes, I have tried different preconditioners, but in a naïve way, i.e., looping through possible options using `-pc_type <option>` command line.</span><span style="font-size:12.0pt;font-family:"Aptos",sans-serif"><o:p></o:p></span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Aptos",sans-serif">But no, not in a meaningful way because the lack of understanding of the connection between physics (the problem we are dealing with) to math (the correct combination of those
preconditioners).</span><span style="font-size:12.0pt;font-family:"Aptos",sans-serif"><o:p></o:p></span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Aptos",sans-serif"> </span><span style="font-size:12.0pt;font-family:"Aptos",sans-serif"><o:p></o:p></span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Aptos",sans-serif">-Ling</span><span style="font-size:12.0pt;font-family:"Aptos",sans-serif"><o:p></o:p></span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Aptos",sans-serif"> </span><span style="font-size:12.0pt;font-family:"Aptos",sans-serif"><o:p></o:p></span></p>
</div>
<div id="mail-editor-reference-message-container">
<div>
<div style="border:none;border-top:solid windowtext 1.0pt;padding:3.0pt 0cm 0cm 0cm;border-color:currentcolor currentcolor;border-image: none">
<p class="MsoNormal" style="margin-bottom:12.0pt"><b><span style="font-size:12.0pt;font-family:"Aptos",sans-serif">From:<span class="apple-converted-space"> </span></span></b><span style="font-size:12.0pt;font-family:"Aptos",sans-serif">Barry Smith <<a href="mailto:bsmith@petsc.dev">bsmith@petsc.dev</a>><br>
<b>Date:<span class="apple-converted-space"> </span></b>Thursday, March 28, 2024 at 1:09 PM<br>
<b>To:<span class="apple-converted-space"> </span></b>Zou, Ling <<a href="mailto:lzou@anl.gov">lzou@anl.gov</a>><br>
<b>Cc:<span class="apple-converted-space"> </span></b>Zhang, Hong <<a href="mailto:hzhang@mcs.anl.gov">hzhang@mcs.anl.gov</a>>,<span class="apple-converted-space"> </span><a href="mailto:petsc-users@mcs.anl.gov">petsc-users@mcs.anl.gov</a><span class="apple-converted-space"> </span><<a href="mailto:petsc-users@mcs.anl.gov">petsc-users@mcs.anl.gov</a>><br>
<b>Subject:<span class="apple-converted-space"> </span></b>Re: [petsc-users] Does ILU(15) still make sense or should just use LU?<o:p></o:p></span></p>
</div>
<div>
<div>
<p class="MsoNormal"><span style="font-size:1.0pt;font-family:"Aptos",sans-serif;color:white">1 million is possible for direct solvers using PETSc with the MUMPS direct solver when you cannot get a preconditioner to work well for your problems. ILU are not
very robust preconditioners and I would not rely on them. Have you investigated<span class="apple-converted-space"> </span></span><span style="font-size:12.0pt;font-family:"Aptos",sans-serif"><o:p></o:p></span></p>
</div>
</div>
<div>
<div>
<p class="MsoNormal"><span style="font-size:1.0pt;font-family:"Aptos",sans-serif;color:white">ZjQcmQRYFpfptBannerStart</span><span style="font-size:12.0pt;font-family:"Aptos",sans-serif"><o:p></o:p></span></p>
</div>
</div>
<div id="pfptBannerc5qmvvr">
<div id="pfptBannerc5qmvvr">
<div>
<p class="MsoNormal" style="line-height:13.5pt;background:#D0D8DC"><b><span style="font-size:12.0pt;font-family:"Arial",sans-serif;color:black">This Message Is From an External Sender</span></b><span style="font-size:12.0pt;font-family:"Aptos",sans-serif"><o:p></o:p></span></p>
</div>
</div>
<div id="pfptBannerc5qmvvr">
<div>
<p class="MsoNormal" style="line-height:13.5pt;background:#D0D8DC"><span style="font-size:12.0pt;font-family:"Arial",sans-serif;color:black">This message came from outside your organization.</span><span style="font-size:12.0pt;font-family:"Aptos",sans-serif"><o:p></o:p></span></p>
</div>
</div>
</div>
<div>
<div>
<p class="MsoNormal" style="background:#D0D8DC"><span style="font-size:12.0pt;font-family:"Aptos",sans-serif;color:black"> </span><span style="font-size:12.0pt;font-family:"Aptos",sans-serif"><o:p></o:p></span></p>
</div>
</div>
<div>
<div>
<p class="MsoNormal"><span style="font-size:1.0pt;font-family:"Aptos",sans-serif;color:white">ZjQcmQRYFpfptBannerEnd</span><span style="font-size:12.0pt;font-family:"Aptos",sans-serif"><o:p></o:p></span></p>
</div>
</div>
<div>
<div>
<p class="MsoNormal"><span style="font-size:12.0pt;font-family:"Aptos",sans-serif"> <o:p></o:p></span></p>
</div>
</div>
<div>
<p class="MsoNormal"><span style="font-size:12.0pt;font-family:"Aptos",sans-serif"> 1 million is possible for direct solvers using PETSc with the MUMPS direct solver when you cannot get a preconditioner to work well for your problems.<o:p></o:p></span></p>
</div>
<div>
<div>
<p class="MsoNormal"><span style="font-size:12.0pt;font-family:"Aptos",sans-serif"> <o:p></o:p></span></p>
</div>
</div>
<div>
<div>
<p class="MsoNormal"><span style="font-size:12.0pt;font-family:"Aptos",sans-serif"> ILU are not very robust preconditioners and I would not rely on them. Have you investigated other preconditioners in PETSc, PCGAMG, PCASM, PCFIELDSPLIT or some combination
of these preconditioners work for many problems, though certainly not all.<o:p></o:p></span></p>
</div>
</div>
<div>
<div>
<p class="MsoNormal"><span style="font-size:12.0pt;font-family:"Aptos",sans-serif"> <o:p></o:p></span></p>
</div>
<div>
<div>
<p class="MsoNormal" style="margin-bottom:12.0pt"><span style="font-size:12.0pt;font-family:"Aptos",sans-serif"><br>
<br>
<o:p></o:p></span></p>
</div>
<blockquote style="margin-top:5.0pt;margin-bottom:5.0pt">
<div>
<div>
<p class="MsoNormal"><span style="font-size:12.0pt;font-family:"Aptos",sans-serif">On Mar 28, 2024, at 1:14</span><span style="font-size:12.0pt;font-family:"Arial",sans-serif"> </span><span style="font-size:12.0pt;font-family:"Aptos",sans-serif">PM, Zou, Ling
<<a href="mailto:lzou@anl.gov">lzou@anl.gov</a>> wrote:<o:p></o:p></span></p>
</div>
</div>
<div>
<p class="MsoNormal"><span style="font-size:12.0pt;font-family:"Aptos",sans-serif"> <o:p></o:p></span></p>
</div>
<div>
<div>
<div>
<div>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Aptos",sans-serif">Thank you, Barry.</span><span style="font-size:12.0pt;font-family:"Aptos",sans-serif"><o:p></o:p></span></p>
</div>
</div>
<div>
<div>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Aptos",sans-serif">Yeah, this is unfortunate given that the problem we are handling is quite heterogeneous (in both mesh and physics).</span><span style="font-size:12.0pt;font-family:"Aptos",sans-serif"><o:p></o:p></span></p>
</div>
</div>
<div>
<div>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Aptos",sans-serif">I expect that our problem sizes will be mostly smaller than 1 million DOF, should LU still be a practical solution? Can it scale well if we choose to run the problem in a parallel
way?</span><span style="font-size:12.0pt;font-family:"Aptos",sans-serif"><o:p></o:p></span></p>
</div>
</div>
<div>
<div>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Aptos",sans-serif"> </span><span style="font-size:12.0pt;font-family:"Aptos",sans-serif"><o:p></o:p></span></p>
</div>
</div>
<div>
<div>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Aptos",sans-serif">PS1:<span class="apple-converted-space"> </span></span><span style="font-size:12.0pt;font-family:"Aptos",sans-serif">-ksp_norm_type unpreconditioned did not work as the true
residual did not go down, even with 300 linear iterations.<o:p></o:p></span></p>
</div>
</div>
<div>
<div>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Aptos",sans-serif">PS2: what do you think if it will be beneficial to have more detailed discussions (e.g., a presentation?) on the problem we are solving to seek more advice?</span><span style="font-size:12.0pt;font-family:"Aptos",sans-serif"><o:p></o:p></span></p>
</div>
</div>
<div>
<div>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Aptos",sans-serif"> </span><span style="font-size:12.0pt;font-family:"Aptos",sans-serif"><o:p></o:p></span></p>
</div>
</div>
<div>
<div>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Aptos",sans-serif">-Ling</span><span style="font-size:12.0pt;font-family:"Aptos",sans-serif"><o:p></o:p></span></p>
</div>
</div>
<div>
<div>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Aptos",sans-serif"> </span><span style="font-size:12.0pt;font-family:"Aptos",sans-serif"><o:p></o:p></span></p>
</div>
</div>
<div id="mail-editor-reference-message-container">
<div>
<div style="border:none;border-top:solid windowtext 1.0pt;padding:3.0pt 0cm 0cm 0cm;border-image: none;border-color:currentcolor">
<p class="MsoNormal" style="margin-bottom:12.0pt"><b><span style="font-size:12.0pt;font-family:"Aptos",sans-serif">From:<span class="apple-converted-space"> </span></span></b><span style="font-size:12.0pt;font-family:"Aptos",sans-serif">Barry Smith <<a href="mailto:bsmith@petsc.dev">bsmith@petsc.dev</a>><br>
<b>Date:<span class="apple-converted-space"> </span></b>Thursday, March 28, 2024 at 11:14 AM<br>
<b>To:<span class="apple-converted-space"> </span></b>Zou, Ling <<a href="mailto:lzou@anl.gov">lzou@anl.gov</a>><br>
<b>Cc:<span class="apple-converted-space"> </span></b>Zhang, Hong <<a href="mailto:hzhang@mcs.anl.gov">hzhang@mcs.anl.gov</a>>,<span class="apple-converted-space"> </span><a href="mailto:petsc-users@mcs.anl.gov">petsc-users@mcs.anl.gov</a><span class="apple-converted-space"> </span><<a href="mailto:petsc-users@mcs.anl.gov">petsc-users@mcs.anl.gov</a>><br>
<b>Subject:<span class="apple-converted-space"> </span></b>Re: [petsc-users] Does ILU(15) still make sense or should just use LU?<o:p></o:p></span></p>
</div>
<div>
<div>
<div>
<p class="MsoNormal"><span style="font-size:1.0pt;font-family:"Aptos",sans-serif;color:white">This is a bad situation, the solver is not really converging. This can happen with ILU() sometimes, it so badly scales things that the preconditioned residual decreases
a lot but the true residual is not really getting smaller. Since your matrices<span class="apple-converted-space"> </span></span><span style="font-size:12.0pt;font-family:"Aptos",sans-serif"><o:p></o:p></span></p>
</div>
</div>
</div>
<div>
<div>
<div>
<p class="MsoNormal"><span style="font-size:1.0pt;font-family:"Aptos",sans-serif;color:white">ZjQcmQRYFpfptBannerStart</span><span style="font-size:12.0pt;font-family:"Aptos",sans-serif"><o:p></o:p></span></p>
</div>
</div>
</div>
<div id="pfptBannerwwham91">
<div id="pfptBannerwwham91">
<div>
<div>
<p class="MsoNormal" style="line-height:13.5pt;background:#D0D8DC"><b><span style="font-size:12.0pt;font-family:"Arial",sans-serif;color:black">This Message Is From an External Sender</span></b><span style="font-size:12.0pt;font-family:"Aptos",sans-serif"><o:p></o:p></span></p>
</div>
</div>
</div>
<div id="pfptBannerwwham91">
<div>
<div>
<p class="MsoNormal" style="line-height:13.5pt;background:#D0D8DC"><span style="font-size:12.0pt;font-family:"Arial",sans-serif;color:black">This message came from outside your organization.</span><span style="font-size:12.0pt;font-family:"Aptos",sans-serif"><o:p></o:p></span></p>
</div>
</div>
</div>
</div>
<div>
<div>
<div>
<p class="MsoNormal" style="background:#D0D8DC"><span style="font-size:12.0pt;font-family:"Aptos",sans-serif;color:black"> </span><span style="font-size:12.0pt;font-family:"Aptos",sans-serif"><o:p></o:p></span></p>
</div>
</div>
</div>
<div>
<div>
<div>
<p class="MsoNormal"><span style="font-size:1.0pt;font-family:"Aptos",sans-serif;color:white">ZjQcmQRYFpfptBannerEnd</span><span style="font-size:12.0pt;font-family:"Aptos",sans-serif"><o:p></o:p></span></p>
</div>
</div>
</div>
<div>
<div>
<div>
<p class="MsoNormal"><span style="font-size:12.0pt;font-family:"Aptos",sans-serif"> <o:p></o:p></span></p>
</div>
</div>
</div>
<div>
<div>
<p class="MsoNormal"><span style="font-size:12.0pt;font-family:"Aptos",sans-serif"> This is a bad situation, the solver is not really converging. This can happen with ILU() sometimes, it so badly scales things that the preconditioned residual decreases a
lot but the true residual is not really getting smaller. Since your matrices are small best to stick to LU.<o:p></o:p></span></p>
</div>
</div>
<div>
<div>
<div>
<p class="MsoNormal"><span style="font-size:12.0pt;font-family:"Aptos",sans-serif"> <o:p></o:p></span></p>
</div>
</div>
</div>
<div>
<div>
<div>
<p class="MsoNormal"><span style="font-size:12.0pt;font-family:"Aptos",sans-serif"> You can use -ksp_norm_type unpreconditioned to force the convergence test to use the true residual for a convergence test and the solver will discover that it is not converging.<o:p></o:p></span></p>
</div>
</div>
</div>
<div>
<div>
<div>
<p class="MsoNormal"><span style="font-size:12.0pt;font-family:"Aptos",sans-serif"> <o:p></o:p></span></p>
</div>
</div>
</div>
<div>
<div>
<div>
<p class="MsoNormal"><span style="font-size:12.0pt;font-family:"Aptos",sans-serif"> Barry<o:p></o:p></span></p>
</div>
</div>
</div>
<div>
<div>
<div>
<p class="MsoNormal"><span style="font-size:12.0pt;font-family:"Aptos",sans-serif"> <o:p></o:p></span></p>
</div>
</div>
<div>
<div>
<p class="MsoNormal" style="margin-bottom:12.0pt"><span style="font-size:12.0pt;font-family:"Aptos",sans-serif"> <o:p></o:p></span></p>
</div>
<blockquote style="margin-top:5.0pt;margin-bottom:5.0pt">
<div>
<div>
<div>
<p class="MsoNormal"><span style="font-size:12.0pt;font-family:"Aptos",sans-serif">On Mar 28, 2024, at 11:43</span><span style="font-size:12.0pt;font-family:"Arial",sans-serif"> </span><span style="font-size:12.0pt;font-family:"Aptos",sans-serif">AM, Zou, Ling
via petsc-users <<a href="mailto:petsc-users@mcs.anl.gov">petsc-users@mcs.anl.gov</a>> wrote:<o:p></o:p></span></p>
</div>
</div>
</div>
<div>
<div>
<p class="MsoNormal"><span style="font-size:12.0pt;font-family:"Aptos",sans-serif"> <o:p></o:p></span></p>
</div>
</div>
<div>
<div>
<div>
<div>
<div>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Aptos",sans-serif">Hong, thanks! That makes perfect sense.</span><span style="font-size:12.0pt;font-family:"Aptos",sans-serif"><o:p></o:p></span></p>
</div>
</div>
</div>
<div>
<div>
<div>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Aptos",sans-serif">A follow up question about ILU.</span><span style="font-size:12.0pt;font-family:"Aptos",sans-serif"><o:p></o:p></span></p>
</div>
</div>
</div>
<div>
<div>
<div>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Aptos",sans-serif"> </span><span style="font-size:12.0pt;font-family:"Aptos",sans-serif"><o:p></o:p></span></p>
</div>
</div>
</div>
<div>
<div>
<div>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Aptos",sans-serif">The following is the performance of ILU(5). Note that each KPS solving reports converged but as the output shows, the preconditioned residual does while true residual does not.
Is there any way this performance could be improved?</span><span style="font-size:12.0pt;font-family:"Aptos",sans-serif"><o:p></o:p></span></p>
</div>
</div>
</div>
<div>
<div>
<div>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Aptos",sans-serif">Background: the preconditioning matrix is finite difference generated, and should be exact.</span><span style="font-size:12.0pt;font-family:"Aptos",sans-serif"><o:p></o:p></span></p>
</div>
</div>
</div>
<div>
<div>
<div>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Aptos",sans-serif"> </span><span style="font-size:12.0pt;font-family:"Aptos",sans-serif"><o:p></o:p></span></p>
</div>
</div>
</div>
<div>
<div>
<div>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Aptos",sans-serif">-Ling</span><span style="font-size:12.0pt;font-family:"Aptos",sans-serif"><o:p></o:p></span></p>
</div>
</div>
</div>
<div>
<div>
<div>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Aptos",sans-serif"> </span><span style="font-size:12.0pt;font-family:"Aptos",sans-serif"><o:p></o:p></span></p>
</div>
</div>
</div>
<div>
<div>
<div>
<p class="MsoNormal" style="background:black"><span class="s1"><span style="font-size:9.0pt;font-family:"Andale Mono";color:white">Time Step 21, time = -491.75, dt = 1</span></span><span style="font-size:12.0pt;font-family:"Aptos",sans-serif"><o:p></o:p></span></p>
</div>
</div>
</div>
<div>
<div>
<div>
<p class="MsoNormal" style="background:black"><span class="apple-converted-space"><span style="font-size:9.0pt;font-family:"Andale Mono";color:white"> </span></span><span class="s1"><span style="font-size:9.0pt;font-family:"Andale Mono";color:white">NL Step
=</span></span><span class="apple-converted-space"><span style="font-size:9.0pt;font-family:"Andale Mono";color:white"> </span></span><span class="s1"><span style="font-size:9.0pt;font-family:"Andale Mono";color:white">0, fnorm =</span></span><span class="apple-converted-space"><span style="font-size:9.0pt;font-family:"Andale Mono";color:white"> </span></span><span class="s1"><span style="font-size:9.0pt;font-family:"Andale Mono";color:white">6.98749E+01</span></span><span style="font-size:12.0pt;font-family:"Aptos",sans-serif"><o:p></o:p></span></p>
</div>
</div>
</div>
<div>
<div>
<div>
<p class="MsoNormal" style="background:black"><span class="apple-converted-space"><span style="font-size:9.0pt;font-family:"Andale Mono";color:white"> </span></span><span class="s1"><span style="font-size:9.0pt;font-family:"Andale Mono";color:white">0 KSP
preconditioned resid norm 1.684131526824e+04 true resid norm 6.987489798042e+01 ||r(i)||/||b|| 1.000000000000e+00</span></span><span style="font-size:12.0pt;font-family:"Aptos",sans-serif"><o:p></o:p></span></p>
</div>
</div>
</div>
<div>
<div>
<div>
<p class="MsoNormal" style="background:black"><span class="apple-converted-space"><span style="font-size:9.0pt;font-family:"Andale Mono";color:white"> </span></span><span class="s1"><span style="font-size:9.0pt;font-family:"Andale Mono";color:white">1 KSP
preconditioned resid norm 5.970568556551e+02 true resid norm 6.459553545222e+01 ||r(i)||/||b|| 9.244455064582e-01</span></span><span style="font-size:12.0pt;font-family:"Aptos",sans-serif"><o:p></o:p></span></p>
</div>
</div>
</div>
<div>
<div>
<div>
<p class="MsoNormal" style="background:black"><span class="apple-converted-space"><span style="font-size:9.0pt;font-family:"Andale Mono";color:white"> </span></span><span class="s1"><span style="font-size:9.0pt;font-family:"Andale Mono";color:white">2 KSP
preconditioned resid norm 3.349113985192e+02 true resid norm 7.250836872274e+01 ||r(i)||/||b|| 1.037688366186e+00</span></span><span style="font-size:12.0pt;font-family:"Aptos",sans-serif"><o:p></o:p></span></p>
</div>
</div>
</div>
<div>
<div>
<div>
<p class="MsoNormal" style="background:black"><span class="apple-converted-space"><span style="font-size:9.0pt;font-family:"Andale Mono";color:white"> </span></span><span class="s1"><span style="font-size:9.0pt;font-family:"Andale Mono";color:white">3 KSP
preconditioned resid norm 3.290585904777e+01 true resid norm 1.186282435163e+02 ||r(i)||/||b|| 1.697723316169e+00</span></span><span style="font-size:12.0pt;font-family:"Aptos",sans-serif"><o:p></o:p></span></p>
</div>
</div>
</div>
<div>
<div>
<div>
<p class="MsoNormal" style="background:black"><span class="apple-converted-space"><span style="font-size:9.0pt;font-family:"Andale Mono";color:white"> </span></span><span class="s1"><span style="font-size:9.0pt;font-family:"Andale Mono";color:white">4 KSP
preconditioned resid norm 8.530606201233e+00 true resid norm 4.088729421459e+01 ||r(i)||/||b|| 5.851499665310e-01</span></span><span style="font-size:12.0pt;font-family:"Aptos",sans-serif"><o:p></o:p></span></p>
</div>
</div>
</div>
<div>
<div>
<div>
<p class="MsoNormal" style="background:black"><span class="apple-converted-space"><span style="font-size:9.0pt;font-family:"Andale Mono";color:white"> </span></span><span class="s1"><span style="font-size:9.0pt;font-family:"Andale Mono";color:white">Linear
solve converged due to CONVERGED_RTOL iterations 4</span></span><span style="font-size:12.0pt;font-family:"Aptos",sans-serif"><o:p></o:p></span></p>
</div>
</div>
</div>
<div>
<div>
<div>
<p class="MsoNormal" style="background:black"><span class="apple-converted-space"><span style="font-size:9.0pt;font-family:"Andale Mono";color:white"> </span></span><span class="s1"><span style="font-size:9.0pt;font-family:"Andale Mono";color:white">NL Step
=</span></span><span class="apple-converted-space"><span style="font-size:9.0pt;font-family:"Andale Mono";color:white"> </span></span><span class="s1"><span style="font-size:9.0pt;font-family:"Andale Mono";color:white">1, fnorm =</span></span><span class="apple-converted-space"><span style="font-size:9.0pt;font-family:"Andale Mono";color:white"> </span></span><span class="s1"><span style="font-size:9.0pt;font-family:"Andale Mono";color:white">4.08788E+01</span></span><span style="font-size:12.0pt;font-family:"Aptos",sans-serif"><o:p></o:p></span></p>
</div>
</div>
</div>
<div>
<div>
<div>
<p class="MsoNormal" style="background:black"><span class="apple-converted-space"><span style="font-size:9.0pt;font-family:"Andale Mono";color:white"> </span></span><span class="s1"><span style="font-size:9.0pt;font-family:"Andale Mono";color:white">0 KSP
preconditioned resid norm 1.851047973094e+03 true resid norm 4.087882723223e+01 ||r(i)||/||b|| 1.000000000000e+00</span></span><span style="font-size:12.0pt;font-family:"Aptos",sans-serif"><o:p></o:p></span></p>
</div>
</div>
</div>
<div>
<div>
<div>
<p class="MsoNormal" style="background:black"><span class="apple-converted-space"><span style="font-size:9.0pt;font-family:"Andale Mono";color:white"> </span></span><span class="s1"><span style="font-size:9.0pt;font-family:"Andale Mono";color:white">1 KSP
preconditioned resid norm 3.696809614513e+01 true resid norm 2.720016413105e+01 ||r(i)||/||b|| 6.653851387793e-01</span></span><span style="font-size:12.0pt;font-family:"Aptos",sans-serif"><o:p></o:p></span></p>
</div>
</div>
</div>
<div>
<div>
<div>
<p class="MsoNormal" style="background:black"><span class="apple-converted-space"><span style="font-size:9.0pt;font-family:"Andale Mono";color:white"> </span></span><span class="s1"><span style="font-size:9.0pt;font-family:"Andale Mono";color:white">2 KSP
preconditioned resid norm 5.751891392534e+00 true resid norm 3.326338240872e+01 ||r(i)||/||b|| 8.137068663873e-01</span></span><span style="font-size:12.0pt;font-family:"Aptos",sans-serif"><o:p></o:p></span></p>
</div>
</div>
</div>
<div>
<div>
<div>
<p class="MsoNormal" style="background:black"><span class="apple-converted-space"><span style="font-size:9.0pt;font-family:"Andale Mono";color:white"> </span></span><span class="s1"><span style="font-size:9.0pt;font-family:"Andale Mono";color:white">3 KSP
preconditioned resid norm 8.540729397958e-01 true resid norm 8.672410748720e+00 ||r(i)||/||b|| 2.121492062249e-01</span></span><span style="font-size:12.0pt;font-family:"Aptos",sans-serif"><o:p></o:p></span></p>
</div>
</div>
</div>
<div>
<div>
<div>
<p class="MsoNormal" style="background:black"><span class="apple-converted-space"><span style="font-size:9.0pt;font-family:"Andale Mono";color:white"> </span></span><span class="s1"><span style="font-size:9.0pt;font-family:"Andale Mono";color:white">Linear
solve converged due to CONVERGED_RTOL iterations 3</span></span><span style="font-size:12.0pt;font-family:"Aptos",sans-serif"><o:p></o:p></span></p>
</div>
</div>
</div>
<div>
<div>
<div>
<p class="MsoNormal" style="background:black"><span class="apple-converted-space"><span style="font-size:9.0pt;font-family:"Andale Mono";color:white"> </span></span><span class="s1"><span style="font-size:9.0pt;font-family:"Andale Mono";color:white">NL Step
=</span></span><span class="apple-converted-space"><span style="font-size:9.0pt;font-family:"Andale Mono";color:white"> </span></span><span class="s1"><span style="font-size:9.0pt;font-family:"Andale Mono";color:white">2, fnorm =</span></span><span class="apple-converted-space"><span style="font-size:9.0pt;font-family:"Andale Mono";color:white"> </span></span><span class="s1"><span style="font-size:9.0pt;font-family:"Andale Mono";color:white">8.67124E+00</span></span><span style="font-size:12.0pt;font-family:"Aptos",sans-serif"><o:p></o:p></span></p>
</div>
</div>
</div>
<div>
<div>
<div>
<p class="MsoNormal" style="background:black"><span class="apple-converted-space"><span style="font-size:9.0pt;font-family:"Andale Mono";color:white"> </span></span><span class="s1"><span style="font-size:9.0pt;font-family:"Andale Mono";color:white">0 KSP
preconditioned resid norm 5.511333966852e+00 true resid norm 8.671237519593e+00 ||r(i)||/||b|| 1.000000000000e+00</span></span><span style="font-size:12.0pt;font-family:"Aptos",sans-serif"><o:p></o:p></span></p>
</div>
</div>
</div>
<div>
<div>
<div>
<p class="MsoNormal" style="background:black"><span class="apple-converted-space"><span style="font-size:9.0pt;font-family:"Andale Mono";color:white"> </span></span><span class="s1"><span style="font-size:9.0pt;font-family:"Andale Mono";color:white">1 KSP
preconditioned resid norm 1.174962622023e+00 true resid norm 8.731034658309e+00 ||r(i)||/||b|| 1.006896032842e+00</span></span><span style="font-size:12.0pt;font-family:"Aptos",sans-serif"><o:p></o:p></span></p>
</div>
</div>
</div>
<div>
<div>
<div>
<p class="MsoNormal" style="background:black"><span class="apple-converted-space"><span style="font-size:9.0pt;font-family:"Andale Mono";color:white"> </span></span><span class="s1"><span style="font-size:9.0pt;font-family:"Andale Mono";color:white">2 KSP
preconditioned resid norm 1.104604471016e+00 true resid norm 1.018397505468e+01 ||r(i)||/||b|| 1.174454630227e+00</span></span><span style="font-size:12.0pt;font-family:"Aptos",sans-serif"><o:p></o:p></span></p>
</div>
</div>
</div>
<div>
<div>
<div>
<p class="MsoNormal" style="background:black"><span class="apple-converted-space"><span style="font-size:9.0pt;font-family:"Andale Mono";color:white"> </span></span><span class="s1"><span style="font-size:9.0pt;font-family:"Andale Mono";color:white">3 KSP
preconditioned resid norm 4.257063674222e-01 true resid norm 4.023093124996e+00 ||r(i)||/||b|| 4.639583584126e-01</span></span><span style="font-size:12.0pt;font-family:"Aptos",sans-serif"><o:p></o:p></span></p>
</div>
</div>
</div>
<div>
<div>
<div>
<p class="MsoNormal" style="background:black"><span class="apple-converted-space"><span style="font-size:9.0pt;font-family:"Andale Mono";color:white"> </span></span><span class="s1"><span style="font-size:9.0pt;font-family:"Andale Mono";color:white">4 KSP
preconditioned resid norm 1.023038868263e-01 true resid norm 2.365298462869e+00 ||r(i)||/||b|| 2.727751901068e-01</span></span><span style="font-size:12.0pt;font-family:"Aptos",sans-serif"><o:p></o:p></span></p>
</div>
</div>
</div>
<div>
<div>
<div>
<p class="MsoNormal" style="background:black"><span class="apple-converted-space"><span style="font-size:9.0pt;font-family:"Andale Mono";color:white"> </span></span><span class="s1"><span style="font-size:9.0pt;font-family:"Andale Mono";color:white">5 KSP
preconditioned resid norm 4.073772638935e-02 true resid norm 2.302623112025e+00 ||r(i)||/||b|| 2.655472309255e-01</span></span><span style="font-size:12.0pt;font-family:"Aptos",sans-serif"><o:p></o:p></span></p>
</div>
</div>
</div>
<div>
<div>
<div>
<p class="MsoNormal" style="background:black"><span class="apple-converted-space"><span style="font-size:9.0pt;font-family:"Andale Mono";color:white"> </span></span><span class="s1"><span style="font-size:9.0pt;font-family:"Andale Mono";color:white">6 KSP
preconditioned resid norm 1.510323179379e-02 true resid norm 2.300216593521e+00 ||r(i)||/||b|| 2.652697020839e-01</span></span><span style="font-size:12.0pt;font-family:"Aptos",sans-serif"><o:p></o:p></span></p>
</div>
</div>
</div>
<div>
<div>
<div>
<p class="MsoNormal" style="background:black"><span class="apple-converted-space"><span style="font-size:9.0pt;font-family:"Andale Mono";color:white"> </span></span><span class="s1"><span style="font-size:9.0pt;font-family:"Andale Mono";color:white">7 KSP
preconditioned resid norm 1.337324816903e-02 true resid norm 2.300057733345e+00 ||r(i)||/||b|| 2.652513817259e-01</span></span><span style="font-size:12.0pt;font-family:"Aptos",sans-serif"><o:p></o:p></span></p>
</div>
</div>
</div>
<div>
<div>
<div>
<p class="MsoNormal" style="background:black"><span class="apple-converted-space"><span style="font-size:9.0pt;font-family:"Andale Mono";color:white"> </span></span><span class="s1"><span style="font-size:9.0pt;font-family:"Andale Mono";color:white">8 KSP
preconditioned resid norm 1.247384902656e-02 true resid norm 2.300456226062e+00 ||r(i)||/||b|| 2.652973374174e-01</span></span><span style="font-size:12.0pt;font-family:"Aptos",sans-serif"><o:p></o:p></span></p>
</div>
</div>
</div>
<div>
<div>
<div>
<p class="MsoNormal" style="background:black"><span class="apple-converted-space"><span style="font-size:9.0pt;font-family:"Andale Mono";color:white"> </span></span><span class="s1"><span style="font-size:9.0pt;font-family:"Andale Mono";color:white">9 KSP
preconditioned resid norm 1.247038855375e-02 true resid norm 2.300532560993e+00 ||r(i)||/||b|| 2.653061406512e-01</span></span><span style="font-size:12.0pt;font-family:"Aptos",sans-serif"><o:p></o:p></span></p>
</div>
</div>
</div>
<div>
<div>
<div>
<p class="MsoNormal" style="background:black"><span class="apple-converted-space"><span style="font-size:9.0pt;font-family:"Andale Mono";color:white"> </span></span><span class="s1"><span style="font-size:9.0pt;font-family:"Andale Mono";color:white">10 KSP
preconditioned resid norm 1.244611343317e-02 true resid norm 2.299441241514e+00 ||r(i)||/||b|| 2.651802855496e-01</span></span><span style="font-size:12.0pt;font-family:"Aptos",sans-serif"><o:p></o:p></span></p>
</div>
</div>
</div>
<div>
<div>
<div>
<p class="MsoNormal" style="background:black"><span class="apple-converted-space"><span style="font-size:9.0pt;font-family:"Andale Mono";color:white"> </span></span><span class="s1"><span style="font-size:9.0pt;font-family:"Andale Mono";color:white">11 KSP
preconditioned resid norm 1.227243209527e-02 true resid norm 2.273668115236e+00 ||r(i)||/||b|| 2.622080308720e-01</span></span><span style="font-size:12.0pt;font-family:"Aptos",sans-serif"><o:p></o:p></span></p>
</div>
</div>
</div>
<div>
<div>
<div>
<p class="MsoNormal" style="background:black"><span class="apple-converted-space"><span style="font-size:9.0pt;font-family:"Andale Mono";color:white"> </span></span><span class="s1"><span style="font-size:9.0pt;font-family:"Andale Mono";color:white">12 KSP
preconditioned resid norm 1.172621459354e-02 true resid norm 2.113927895437e+00 ||r(i)||/||b|| 2.437861828442e-01</span></span><span style="font-size:12.0pt;font-family:"Aptos",sans-serif"><o:p></o:p></span></p>
</div>
</div>
</div>
<div>
<div>
<div>
<p class="MsoNormal" style="background:black"><span class="apple-converted-space"><span style="font-size:9.0pt;font-family:"Andale Mono";color:white"> </span></span><span class="s1"><span style="font-size:9.0pt;font-family:"Andale Mono";color:white">13 KSP
preconditioned resid norm 2.880752338189e-03 true resid norm 1.076190247720e-01 ||r(i)||/||b|| 1.241103412620e-02</span></span><span style="font-size:12.0pt;font-family:"Aptos",sans-serif"><o:p></o:p></span></p>
</div>
</div>
</div>
<div>
<div>
<div>
<p class="MsoNormal" style="background:black"><span class="apple-converted-space"><span style="font-size:9.0pt;font-family:"Andale Mono";color:white"> </span></span><span class="s1"><span style="font-size:9.0pt;font-family:"Andale Mono";color:white">Linear
solve converged due to CONVERGED_RTOL iterations 13</span></span><span style="font-size:12.0pt;font-family:"Aptos",sans-serif"><o:p></o:p></span></p>
</div>
</div>
</div>
<div>
<div>
<div>
<p class="MsoNormal" style="background:black"><span class="apple-converted-space"><span style="font-size:9.0pt;font-family:"Andale Mono";color:white"> </span></span><span class="s1"><span style="font-size:9.0pt;font-family:"Andale Mono";color:white">NL Step
=</span></span><span class="apple-converted-space"><span style="font-size:9.0pt;font-family:"Andale Mono";color:white"> </span></span><span class="s1"><span style="font-size:9.0pt;font-family:"Andale Mono";color:white">3, fnorm =</span></span><span class="apple-converted-space"><span style="font-size:9.0pt;font-family:"Andale Mono";color:white"> </span></span><span class="s1"><span style="font-size:9.0pt;font-family:"Andale Mono";color:white">1.59729E-01</span></span><span style="font-size:12.0pt;font-family:"Aptos",sans-serif"><o:p></o:p></span></p>
</div>
</div>
</div>
<div>
<div>
<div>
<p class="MsoNormal" style="background:black"><span class="apple-converted-space"><span style="font-size:9.0pt;font-family:"Andale Mono";color:white"> </span></span><span class="s1"><span style="font-size:9.0pt;font-family:"Andale Mono";color:white">0 KSP
preconditioned resid norm 1.676948440854e+03 true resid norm 1.597288981238e-01 ||r(i)||/||b|| 1.000000000000e+00</span></span><span style="font-size:12.0pt;font-family:"Aptos",sans-serif"><o:p></o:p></span></p>
</div>
</div>
</div>
<div>
<div>
<div>
<p class="MsoNormal" style="background:black"><span class="apple-converted-space"><span style="font-size:9.0pt;font-family:"Andale Mono";color:white"> </span></span><span class="s1"><span style="font-size:9.0pt;font-family:"Andale Mono";color:white">1 KSP
preconditioned resid norm 2.266131510513e+00 true resid norm 1.819663943811e+00 ||r(i)||/||b|| 1.139220244542e+01</span></span><span style="font-size:12.0pt;font-family:"Aptos",sans-serif"><o:p></o:p></span></p>
</div>
</div>
</div>
<div>
<div>
<div>
<p class="MsoNormal" style="background:black"><span class="apple-converted-space"><span style="font-size:9.0pt;font-family:"Andale Mono";color:white"> </span></span><span class="s1"><span style="font-size:9.0pt;font-family:"Andale Mono";color:white">2 KSP
preconditioned resid norm 2.239911493901e+00 true resid norm 1.923976907755e+00 ||r(i)||/||b|| 1.204526501062e+01</span></span><span style="font-size:12.0pt;font-family:"Aptos",sans-serif"><o:p></o:p></span></p>
</div>
</div>
</div>
<div>
<div>
<div>
<p class="MsoNormal" style="background:black"><span class="apple-converted-space"><span style="font-size:9.0pt;font-family:"Andale Mono";color:white"> </span></span><span class="s1"><span style="font-size:9.0pt;font-family:"Andale Mono";color:white">3 KSP
preconditioned resid norm 1.446859034276e-01 true resid norm 8.692945031946e-01 ||r(i)||/||b|| 5.442312026225e+00</span></span><span style="font-size:12.0pt;font-family:"Aptos",sans-serif"><o:p></o:p></span></p>
</div>
</div>
</div>
<div>
<div>
<div>
<p class="MsoNormal" style="background:black"><span class="apple-converted-space"><span style="font-size:9.0pt;font-family:"Andale Mono";color:white"> </span></span><span class="s1"><span style="font-size:9.0pt;font-family:"Andale Mono";color:white">Linear
solve converged due to CONVERGED_RTOL iterations 3</span></span><span style="font-size:12.0pt;font-family:"Aptos",sans-serif"><o:p></o:p></span></p>
</div>
</div>
</div>
<div>
<div>
<div>
<p class="MsoNormal" style="background:black"><span class="apple-converted-space"><span style="font-size:9.0pt;font-family:"Andale Mono";color:white"> </span></span><span class="s1"><span style="font-size:9.0pt;font-family:"Andale Mono";color:white">NL Step
=</span></span><span class="apple-converted-space"><span style="font-size:9.0pt;font-family:"Andale Mono";color:white"> </span></span><span class="s1"><span style="font-size:9.0pt;font-family:"Andale Mono";color:white">4, fnorm =</span></span><span class="apple-converted-space"><span style="font-size:9.0pt;font-family:"Andale Mono";color:white"> </span></span><span class="s1"><span style="font-size:9.0pt;font-family:"Andale Mono";color:white">1.59564E-01</span></span><span style="font-size:12.0pt;font-family:"Aptos",sans-serif"><o:p></o:p></span></p>
</div>
</div>
</div>
<div>
<div>
<div>
<p class="MsoNormal" style="background:black"><span class="apple-converted-space"><span style="font-size:9.0pt;font-family:"Andale Mono";color:white"> </span></span><span class="s1"><span style="font-size:9.0pt;font-family:"Andale Mono";color:white">0 KSP
preconditioned resid norm 1.509663716414e+03 true resid norm 1.595641817504e-01 ||r(i)||/||b|| 1.000000000000e+00</span></span><span style="font-size:12.0pt;font-family:"Aptos",sans-serif"><o:p></o:p></span></p>
</div>
</div>
</div>
<div>
<div>
<div>
<p class="MsoNormal" style="background:black"><span class="apple-converted-space"><span style="font-size:9.0pt;font-family:"Andale Mono";color:white"> </span></span><span class="s1"><span style="font-size:9.0pt;font-family:"Andale Mono";color:white">1 KSP
preconditioned resid norm 1.995956587709e+00 true resid norm 1.712323298361e+00 ||r(i)||/||b|| 1.073125108390e+01</span></span><span style="font-size:12.0pt;font-family:"Aptos",sans-serif"><o:p></o:p></span></p>
</div>
</div>
</div>
<div>
<div>
<div>
<p class="MsoNormal" style="background:black"><span class="apple-converted-space"><span style="font-size:9.0pt;font-family:"Andale Mono";color:white"> </span></span><span class="s1"><span style="font-size:9.0pt;font-family:"Andale Mono";color:white">2 KSP
preconditioned resid norm 1.994336275847e+00 true resid norm 1.741263472491e+00 ||r(i)||/||b|| 1.091262119975e+01</span></span><span style="font-size:12.0pt;font-family:"Aptos",sans-serif"><o:p></o:p></span></p>
</div>
</div>
</div>
<div>
<div>
<div>
<p class="MsoNormal" style="background:black"><span class="apple-converted-space"><span style="font-size:9.0pt;font-family:"Andale Mono";color:white"> </span></span><span class="s1"><span style="font-size:9.0pt;font-family:"Andale Mono";color:white">3 KSP
preconditioned resid norm 1.268035008497e-01 true resid norm 8.197057317360e-01 ||r(i)||/||b|| 5.137153731769e+00</span></span><span style="font-size:12.0pt;font-family:"Aptos",sans-serif"><o:p></o:p></span></p>
</div>
</div>
</div>
<div>
<div>
<div>
<p class="MsoNormal" style="background:black"><span class="apple-converted-space"><span style="font-size:9.0pt;font-family:"Andale Mono";color:white"> </span></span><span class="s1"><span style="font-size:9.0pt;font-family:"Andale Mono";color:white">Linear
solve converged due to CONVERGED_RTOL iterations 3</span></span><span style="font-size:12.0pt;font-family:"Aptos",sans-serif"><o:p></o:p></span></p>
</div>
</div>
</div>
<div>
<div>
<div>
<p class="MsoNormal" style="background:black"><span class="s1"><span style="font-size:9.0pt;font-family:"Andale Mono";color:white">Nonlinear solve did not converge due to DIVERGED_LINE_SEARCH iterations 4</span></span><span style="font-size:12.0pt;font-family:"Aptos",sans-serif"><o:p></o:p></span></p>
</div>
</div>
</div>
<div>
<div>
<div>
<p class="MsoNormal" style="background:black"><span class="apple-converted-space"><span style="font-size:9.0pt;font-family:"Andale Mono";color:#B42419"> </span></span><span class="s1"><span style="font-size:9.0pt;font-family:"Andale Mono";color:#B42419">Solve
Did NOT Converge!</span></span><span style="font-size:12.0pt;font-family:"Aptos",sans-serif"><o:p></o:p></span></p>
</div>
</div>
</div>
<div>
<div>
<div>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Aptos",sans-serif"> </span><span style="font-size:12.0pt;font-family:"Aptos",sans-serif"><o:p></o:p></span></p>
</div>
</div>
</div>
<div>
<div>
<div>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Aptos",sans-serif"> </span><span style="font-size:12.0pt;font-family:"Aptos",sans-serif"><o:p></o:p></span></p>
</div>
</div>
</div>
<div>
<div>
<div>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Aptos",sans-serif"> </span><span style="font-size:12.0pt;font-family:"Aptos",sans-serif"><o:p></o:p></span></p>
</div>
</div>
</div>
<div id="mail-editor-reference-message-container">
<div>
<div style="border:none;border-top:solid windowtext 1.0pt;padding:3.0pt 0cm 0cm 0cm;border-image: none;border-color:currentcolor">
<p class="MsoNormal" style="margin-bottom:12.0pt"><b><span style="font-size:12.0pt;font-family:"Aptos",sans-serif">From:<span class="apple-converted-space"> </span></span></b><span style="font-size:12.0pt;font-family:"Aptos",sans-serif">Zhang, Hong <<a href="mailto:hzhang@mcs.anl.gov">hzhang@mcs.anl.gov</a>><br>
<b>Date:<span class="apple-converted-space"> </span></b>Wednesday, March 27, 2024 at 4:59 PM<br>
<b>To:<span class="apple-converted-space"> </span></b><a href="mailto:petsc-users@mcs.anl.gov">petsc-users@mcs.anl.gov</a><span class="apple-converted-space"> </span><<a href="mailto:petsc-users@mcs.anl.gov">petsc-users@mcs.anl.gov</a>>, Zou, Ling <<a href="mailto:lzou@anl.gov">lzou@anl.gov</a>><br>
<b>Subject:<span class="apple-converted-space"> </span></b>Re: Does ILU(15) still make sense or should just use LU?<o:p></o:p></span></p>
</div>
<div>
<div>
<div>
<div>
<p class="MsoNormal"><span style="font-size:12.0pt;font-family:"Aptos",sans-serif">Ling,<o:p></o:p></span></p>
</div>
</div>
</div>
</div>
<div>
<div>
<div>
<div>
<p class="MsoNormal"><span style="font-size:12.0pt;font-family:"Aptos",sans-serif">ILU(level) is used for saving storage space with more computations. Normally, we use level=1 or 2. It does not make sense to use level 15. If you have sufficient space, LU would
be the best.<o:p></o:p></span></p>
</div>
</div>
</div>
</div>
<div>
<div>
<div>
<div>
<p class="MsoNormal"><span style="font-size:12.0pt;font-family:"Aptos",sans-serif">Hong<o:p></o:p></span></p>
</div>
</div>
</div>
</div>
<div>
<div>
<div>
<div>
<p class="MsoNormal"><span style="font-size:12.0pt;font-family:"Aptos",sans-serif"> <o:p></o:p></span></p>
</div>
</div>
</div>
</div>
<div class="MsoNormal" align="center" style="text-align:center"><span style="font-size:12.0pt;font-family:"Aptos",sans-serif">
<hr size="1" width="100%" align="center">
</span></div>
<div id="divRplyFwdMsg">
<div>
<div>
<div>
<p class="MsoNormal"><b><span style="font-size:11.0pt">From:</span></b><span style="font-size:11.0pt"> petsc-users <<a href="mailto:petsc-users-bounces@mcs.anl.gov">petsc-users-bounces@mcs.anl.gov</a>> on behalf of Zou, Ling via petsc-users <<a href="mailto:petsc-users@mcs.anl.gov">petsc-users@mcs.anl.gov</a>><br>
<b>Sent:</b> Wednesday, March 27, 2024 4:24 PM<br>
<b>To:</b> <a href="mailto:petsc-users@mcs.anl.gov">petsc-users@mcs.anl.gov</a><span class="apple-converted-space"> </span><<a href="mailto:petsc-users@mcs.anl.gov">petsc-users@mcs.anl.gov</a>><br>
<b>Subject:</b> [petsc-users] Does ILU(15) still make sense or should just use LU?</span><span style="font-size:12.0pt;font-family:"Aptos",sans-serif"><o:p></o:p></span></p>
</div>
</div>
</div>
<div>
<div>
<div>
<div>
<p class="MsoNormal"><span style="font-size:12.0pt;font-family:"Aptos",sans-serif"> <o:p></o:p></span></p>
</div>
</div>
</div>
</div>
</div>
<p><span style="font-family:"Helvetica",sans-serif">Hi, I’d like to avoid using LU, but in some cases to use ILU and still converge, I have to go to ILU(15), i.e., `-pc_factor_levels 15`. Does it still make sense, or should I give it up and switch to LU?</span><span style="font-size:13.5pt;font-family:"Helvetica",sans-serif"><o:p></o:p></span></p>
<p><span style="font-family:"Helvetica",sans-serif"> </span><span style="font-size:13.5pt;font-family:"Helvetica",sans-serif"><o:p></o:p></span></p>
<p><span style="font-family:"Helvetica",sans-serif">For this particular case, ~2k DoF, and both ILU(15) and LU perform similarly in terms of wall time.</span><span style="font-size:13.5pt;font-family:"Helvetica",sans-serif"><o:p></o:p></span></p>
<p><span style="font-family:"Helvetica",sans-serif"> </span><span style="font-size:13.5pt;font-family:"Helvetica",sans-serif"><o:p></o:p></span></p>
<p><span style="font-family:"Helvetica",sans-serif">-Ling</span><span style="font-size:13.5pt;font-family:"Helvetica",sans-serif"><o:p></o:p></span></p>
</div>
</div>
</div>
</div>
</blockquote>
</div>
</div>
</div>
</div>
</div>
</div>
</blockquote>
</div>
</div>
</div>
</div>
</div>
</blockquote>
</div>
<p class="MsoNormal"><span style="font-size:11.0pt"><o:p> </o:p></span></p>
</div>
</div>
</div>
This email is confidential between the sender and the intended recipient specified in the message, and may contain proprietary or otherwise protected information of the intended parties. If you received this email in error, please notify the sender immediately
and confirm you have deleted this message. Any unauthorized review, disclosure, retransmission, dissemination or other use of or reliance on this information may be unlawful and is strictly prohibited.
</body>
</html>