<head><!-- BaNnErBlUrFlE-HeAdEr-start -->
<style>
#pfptBannerlhdr3m3 { all: revert !important; display: block !important;
visibility: visible !important; opacity: 1 !important;
background-color: #D0D8DC !important;
max-width: none !important; max-height: none !important }
.pfptPrimaryButtonlhdr3m3:hover, .pfptPrimaryButtonlhdr3m3:focus {
background-color: #b4c1c7 !important; }
.pfptPrimaryButtonlhdr3m3: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;">
On Mon, Jul 1, 2024 at 9: 48 AM Blauth, Sebastian <sebastian. blauth@ itwm. fraunhofer. de> wrote: Dear Matt, thanks a lot for your help. Unfortunately, for me these extra options do not have any effect, I still get the “u” and “p” fieldnames.
</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="pfptBannerlhdr3m3" 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="pfptBannerlhdr3m3" style="all: unset !important; float:left !important; display:block !important; margin: 0px 0px 1px 0px !important; max-width: 600px !important;">
<div id="pfptBannerlhdr3m3" 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="pfptBannerlhdr3m3" 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 dir="ltr">On Mon, Jul 1, 2024 at 9:48 AM Blauth, Sebastian <<a href="mailto:sebastian.blauth@itwm.fraunhofer.de">sebastian.blauth@itwm.fraunhofer.de</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 class="msg-552471448270057200"><div lang="DE" style="overflow-wrap: break-word;"><div class="m_-552471448270057200WordSection1"><p class="MsoNormal"><span>Dear Matt,<u></u><u></u></span></p><p class="MsoNormal"><span><u></u> <u></u></span></p><p class="MsoNormal"><span lang="EN-US">thanks a lot for your help. Unfortunately, for me these extra options do not have any effect, I still get the “u” and “p” fieldnames. Also, this would not help me to get rid of the “c” fieldname – on that level of the fieldsplit I am basically using your approach already, and still it does show up. The output of the -ksp_view is unchanged, so that I do not attach it here again. Maybe I misunderstood you?</span></p></div></div></div></blockquote><div><br></div><div>Oh, we make an exception for single fields, since we think you would want to use the name. I have to make an extra option to shut off naming.</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 class="msg-552471448270057200"><div lang="DE" style="overflow-wrap: break-word;"><div class="m_-552471448270057200WordSection1"><p class="MsoNormal"><span lang="EN-US">Thanks for the help and best regards,<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US">Sebastian<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US"><u></u> <u></u></span></p><p class="MsoNormal">--<u></u><u></u></p><p class="MsoNormal">Dr. Sebastian Blauth<u></u><u></u></p><p class="MsoNormal">Fraunhofer-Institut für<u></u><u></u></p><p class="MsoNormal">Techno- und Wirtschaftsmathematik ITWM<u></u><u></u></p><p class="MsoNormal">Abteilung Transportvorgänge<u></u><u></u></p><p class="MsoNormal">Fraunhofer-Platz 1, 67663 Kaiserslautern<u></u><u></u></p><p class="MsoNormal">Telefon: +49 631 31600-4968<u></u><u></u></p><p class="MsoNormal"><a href="mailto:sebastian.blauth@itwm.fraunhofer.de" target="_blank">sebastian.blauth@itwm.fraunhofer.de</a><u></u><u></u></p><p class="MsoNormal"><a href="https://urldefense.us/v3/__https://www.itwm.fraunhofer.de__;!!G_uCfscf7eWS!e9SeVBzwCGTmMl3gOc_WG4S_zL5JNMZYyiUGhfrhQAN_re34sVynQzfQyxsY8DUyFdPf1HHTM-Pas_5OTyRd$" target="_blank">https://www.itwm.fraunhofer.de</a><u></u><u></u></p><p class="MsoNormal"><span><u></u> <u></u></span></p><div style="border-top:none;border-right:none;border-bottom:none;border-left:1.5pt solid blue;padding:0cm 0cm 0cm 4pt"><div><div style="border-right:none;border-bottom:none;border-left:none;border-top:1pt solid rgb(225,225,225);padding:3pt 0cm 0cm"><p class="MsoNormal"><b><span lang="EN-US">From:</span></b><span lang="EN-US"> Matthew Knepley <<a href="mailto:knepley@gmail.com" target="_blank">knepley@gmail.com</a>> <br><b>Sent:</b> Monday, July 1, 2024 2:27 PM<br><b>To:</b> Blauth, Sebastian <<a href="mailto:sebastian.blauth@itwm.fraunhofer.de" target="_blank">sebastian.blauth@itwm.fraunhofer.de</a>><br><b>Cc:</b> <a href="mailto:petsc-users@mcs.anl.gov" target="_blank">petsc-users@mcs.anl.gov</a><br><b>Subject:</b> Re: [petsc-users] Question regarding naming of fieldsplit splits<u></u><u></u></span></p></div></div><p class="MsoNormal"><u></u> <u></u></p><div><div><p class="MsoNormal">On Fri, Jun 28, 2024 at 4:05 AM Blauth, Sebastian <<a href="mailto:sebastian.blauth@itwm.fraunhofer.de" target="_blank">sebastian.blauth@itwm.fraunhofer.de</a>> wrote:<u></u><u></u></p></div><div><blockquote style="border-top:none;border-right:none;border-bottom:none;border-left:1pt solid rgb(204,204,204);padding:0cm 0cm 0cm 6pt;margin-left:4.8pt;margin-right:0cm"><div><div><div><p class="MsoNormal">Hello everyone,<u></u><u></u></p><p class="MsoNormal"> <u></u><u></u></p><p class="MsoNormal"><span lang="EN-US">I have a question regarding the naming convention using PETSc’s PCFieldsplit. I have been following <a href="https://urldefense.us/v3/__https://lists.mcs.anl.gov/pipermail/petsc-users/2019-January/037262.html__;!!G_uCfscf7eWS!e9SeVBzwCGTmMl3gOc_WG4S_zL5JNMZYyiUGhfrhQAN_re34sVynQzfQyxsY8DUyFdPf1HHTM-Paswa3c8E2$" target="_blank">https://lists.mcs.anl.gov/pipermail/petsc-users/2019-January/037262.html</a> to create a DMShell with FEniCS in order to customize PCFieldsplit for my application. </span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US">I am using the following options, which work nicely for me:</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US"> </span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US">-ksp_type fgmres</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US">-pc_type fieldsplit</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US">-pc_fieldsplit_0_fields 0, 1</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US">-pc_fieldsplit_1_fields 2</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US">-pc_fieldsplit_type additive</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US">-fieldsplit_0_ksp_type fgmres</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US">-fieldsplit_0_pc_type fieldsplit</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US">-fieldsplit_0_pc_fieldsplit_type schur</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US">-fieldsplit_0_pc_fieldsplit_schur_fact_type full</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US">-fieldsplit_0_pc_fieldsplit_schur_precondition selfp</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US">-fieldsplit_0_fieldsplit_u_ksp_type preonly</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US">-fieldsplit_0_fieldsplit_u_pc_type lu</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US">-fieldsplit_0_fieldsplit_p_ksp_type cg</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US">-fieldsplit_0_fieldsplit_p_ksp_rtol 1e-14</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US">-fieldsplit_0_fieldsplit_p_ksp_atol 1e-30</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US">-fieldsplit_0_fieldsplit_p_pc_type icc</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US">-fieldsplit_0_ksp_rtol 1e-14</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US">-fieldsplit_0_ksp_atol 1e-30</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US">-fieldsplit_0_ksp_monitor_true_residual</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US">-fieldsplit_c_ksp_type preonly</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US">-fieldsplit_c_pc_type lu</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US">-ksp_view</span><u></u><u></u></p></div></div></div></blockquote><div><p class="MsoNormal"><u></u> <u></u></p></div><div><p class="MsoNormal">By default, we use the field names, but you can prevent this by specifying the fields by hand, so<u></u><u></u></p></div><div><p class="MsoNormal"><u></u> <u></u></p></div><div><p class="MsoNormal">-fieldsplit_0_pc_fieldsplit_0_fields 0<br>-fieldsplit_0_pc_fieldsplit_1_fields 1<u></u><u></u></p></div><div><p class="MsoNormal"><u></u> <u></u></p></div><div><p class="MsoNormal">should remove the 'u' and 'p' fieldnames. It is somewhat hacky, but I think easier to remember than<u></u><u></u></p></div><div><p class="MsoNormal">some extra option.<u></u><u></u></p></div><div><p class="MsoNormal"><u></u> <u></u></p></div><div><p class="MsoNormal"> Thanks,<u></u><u></u></p></div><div><p class="MsoNormal"><u></u> <u></u></p></div><div><p class="MsoNormal"> Matt<u></u><u></u></p></div><div><p class="MsoNormal"> <u></u><u></u></p></div><blockquote style="border-top:none;border-right:none;border-bottom:none;border-left:1pt solid rgb(204,204,204);padding:0cm 0cm 0cm 6pt;margin-left:4.8pt;margin-right:0cm"><div><div><div><p class="MsoNormal"><span lang="EN-US">Note that this is just an academic example (sorry for the low solver tolerances) to test the approach, consisting of a Stokes equation and some concentration equation (which is not even coupled to Stokes, just for testing).</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US">Completely analogous to <a href="https://urldefense.us/v3/__https://lists.mcs.anl.gov/pipermail/petsc-users/2019-January/037262.html__;!!G_uCfscf7eWS!e9SeVBzwCGTmMl3gOc_WG4S_zL5JNMZYyiUGhfrhQAN_re34sVynQzfQyxsY8DUyFdPf1HHTM-Paswa3c8E2$" target="_blank">https://lists.mcs.anl.gov/pipermail/petsc-users/2019-January/037262.html</a>, I translate my IS’s to a PETSc Section, which is then supplied to a DMShell and assigned to a KSP. I am not so familiar with the code or how / why this works, but it seems to do so perfectly. I name my sections with petsc4py using</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US"> </span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US">section.setFieldName(0, "u")</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US">section.setFieldName(1, "p")</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US">section.setFieldName(2, "c")</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US"> </span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US">However, this is also reflected in the way I can access the fieldsplit options from the command line. My question is: Is there any way of not using the FieldNames specified in python but use the index of the field as defined with “-pc_fieldsplit_0_fields 0, 1” and “-pc_fieldsplit_1_fields 2”, i.e., instead of the prefix “fieldsplit_0_fieldsplit_u” I want to write “fieldsplit_0_fieldsplit_0”, instead of “fieldsplit_0_fieldsplit_p” I want to use “fieldsplit_0_fieldsplit_1”, and instead of “fieldsplit_c” I want to use “fieldsplit_1”. Just changing the names of the fields to</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US"> </span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US">section.setFieldName(0, "0")</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US">section.setFieldName(1, "1")</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US">section.setFieldName(2, "2")</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US"> </span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US">does obviously not work as expected, as it works for velocity and pressure, but not for the concentration – the prefix there is then “fieldsplit_2” and not “fieldsplit_1”. In the docs, I have found <a href="https://urldefense.us/v3/__https://petsc.org/main/manualpages/PC/PCFieldSplitSetFields/__;!!G_uCfscf7eWS!e9SeVBzwCGTmMl3gOc_WG4S_zL5JNMZYyiUGhfrhQAN_re34sVynQzfQyxsY8DUyFdPf1HHTM-Pas1RiSjwn$" target="_blank">https://petsc.org/main/manualpages/PC/PCFieldSplitSetFields/</a> which seems to suggest that the fieldname can potentially be supplied, but I don’t see how to do so from the command line. Also, for the sake of completeness, I attach the output of the solve with “-ksp_view” below. </span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US"> </span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US">Thanks a lot in advance and best regards,</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US">Sebastian</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US"> </span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US"> </span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US">The output of ksp_view is the following:</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US">KSP Object: 1 MPI processes</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US"> type: fgmres</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US"> restart=30, using Classical (unmodified) Gram-Schmidt Orthogonalization with no iterative refinement</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US"> happy breakdown tolerance 1e-30</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US"> maximum iterations=10000, initial guess is zero</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US"> tolerances: relative=1e-05, absolute=1e-11, divergence=10000.</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US"> right preconditioning</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US"> using UNPRECONDITIONED norm type for convergence test</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US">PC Object: 1 MPI processes</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US"> type: fieldsplit</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US"> FieldSplit with ADDITIVE composition: total splits = 2</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US"> Solver info for each split is in the following KSP objects:</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US"> Split number 0 Defined by IS</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US"> KSP Object: (fieldsplit_0_) 1 MPI processes</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US"> type: fgmres</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US"> restart=30, using Classical (unmodified) Gram-Schmidt Orthogonalization with no iterative refinement</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US"> happy breakdown tolerance 1e-30</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US"> maximum iterations=10000, initial guess is zero</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US"> tolerances: relative=1e-14, absolute=1e-30, divergence=10000.</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US"> right preconditioning</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US"> using UNPRECONDITIONED norm type for convergence test</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US"> PC Object: (fieldsplit_0_) 1 MPI processes</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US"> type: fieldsplit</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US"> FieldSplit with Schur preconditioner, factorization FULL</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US"> Preconditioner for the Schur complement formed from Sp, an assembled approximation to S, which uses A00's diagonal's inverse</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US"> Split info:</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US"> Split number 0 Defined by IS</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US"> Split number 1 Defined by IS</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US"> KSP solver for A00 block</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US"> KSP Object: (fieldsplit_0_fieldsplit_u_) 1 MPI processes</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US"> type: preonly</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US"> maximum iterations=10000, initial guess is zero</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US"> tolerances: relative=1e-05, absolute=1e-50, divergence=10000.</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US"> left preconditioning</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US"> using NONE norm type for convergence test</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US"> PC Object: (fieldsplit_0_fieldsplit_u_) 1 MPI processes</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US"> type: lu</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US"> out-of-place factorization</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US"> tolerance for zero pivot 2.22045e-14</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US"> matrix ordering: nd</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US"> factor fill ratio given 5., needed 3.92639</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US"> Factored matrix follows:</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US"> Mat Object: 1 MPI processes</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US"> type: seqaij</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US"> rows=4290, cols=4290</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US"> package used to perform factorization: petsc</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US"> total: nonzeros=375944, allocated nonzeros=375944</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US"> using I-node routines: found 2548 nodes, limit used is 5</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US"> linear system matrix = precond matrix:</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US"> Mat Object: (fieldsplit_0_fieldsplit_u_) 1 MPI processes</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US"> type: seqaij</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US"> rows=4290, cols=4290</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US"> total: nonzeros=95748, allocated nonzeros=95748</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US"> total number of mallocs used during MatSetValues calls=0</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US"> using I-node routines: found 3287 nodes, limit used is 5</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US"> KSP solver for S = A11 - A10 inv(A00) A01 </span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US"> KSP Object: (fieldsplit_0_fieldsplit_p_) 1 MPI processes</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US"> type: cg</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US"> maximum iterations=10000, initial guess is zero</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US"> tolerances: relative=1e-14, absolute=1e-30, divergence=10000.</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US"> left preconditioning</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US"> using PRECONDITIONED norm type for convergence test</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US"> PC Object: (fieldsplit_0_fieldsplit_p_) 1 MPI processes</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US"> type: icc</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US"> out-of-place factorization</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US"> 0 levels of fill</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US"> tolerance for zero pivot 2.22045e-14</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US"> using Manteuffel shift [POSITIVE_DEFINITE]</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US"> matrix ordering: natural</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US"> factor fill ratio given 1., needed 1.</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US"> Factored matrix follows:</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US"> Mat Object: 1 MPI processes</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US"> type: seqsbaij</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US"> rows=561, cols=561</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US"> package used to perform factorization: petsc</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US"> total: nonzeros=5120, allocated nonzeros=5120</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US"> block size is 1</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US"> linear system matrix followed by preconditioner matrix:</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US"> Mat Object: (fieldsplit_0_fieldsplit_p_) 1 MPI processes</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US"> type: schurcomplement</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US"> rows=561, cols=561</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US"> Schur complement A11 - A10 inv(A00) A01</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US"> A11</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US"> Mat Object: (fieldsplit_0_fieldsplit_p_) 1 MPI processes</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US"> type: seqaij</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US"> rows=561, cols=561</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US"> total: nonzeros=3729, allocated nonzeros=3729</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US"> total number of mallocs used during MatSetValues calls=0</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US"> not using I-node routines</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US"> A10</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US"> Mat Object: 1 MPI processes</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US"> type: seqaij</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US"> rows=561, cols=4290</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US"> total: nonzeros=19938, allocated nonzeros=19938</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US"> total number of mallocs used during MatSetValues calls=0</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US"> not using I-node routines</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US"> KSP of A00</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US"> KSP Object: (fieldsplit_0_fieldsplit_u_) 1 MPI processes</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US"> type: preonly</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US"> maximum iterations=10000, initial guess is zero</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US"> tolerances: relative=1e-05, absolute=1e-50, divergence=10000.</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US"> left preconditioning</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US"> using NONE norm type for convergence test</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US"> PC Object: (fieldsplit_0_fieldsplit_u_) 1 MPI processes</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US"> type: lu</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US"> out-of-place factorization</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US"> tolerance for zero pivot 2.22045e-14</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US"> matrix ordering: nd</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US"> factor fill ratio given 5., needed 3.92639</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US"> Factored matrix follows:</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US"> Mat Object: 1 MPI processes</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US"> type: seqaij</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US"> rows=4290, cols=4290</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US"> package used to perform factorization: petsc</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US"> total: nonzeros=375944, allocated nonzeros=375944</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US"> using I-node routines: found 2548 nodes, limit used is 5</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US"> linear system matrix = precond matrix:</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US"> Mat Object: (fieldsplit_0_fieldsplit_u_) 1 MPI processes</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US"> type: seqaij</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US"> rows=4290, cols=4290</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US"> total: nonzeros=95748, allocated nonzeros=95748</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US"> total number of mallocs used during MatSetValues calls=0</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US"> using I-node routines: found 3287 nodes, limit used is 5</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US"> A01</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US"> Mat Object: 1 MPI processes</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US"> type: seqaij</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US"> rows=4290, cols=561</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US"> total: nonzeros=19938, allocated nonzeros=19938</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US"> total number of mallocs used during MatSetValues calls=0</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US"> using I-node routines: found 3287 nodes, limit used is 5</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US"> Mat Object: 1 MPI processes</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US"> type: seqaij</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US"> rows=561, cols=561</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US"> total: nonzeros=9679, allocated nonzeros=9679</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US"> total number of mallocs used during MatSetValues calls=0</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US"> not using I-node routines</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US"> linear system matrix = precond matrix:</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US"> Mat Object: (fieldsplit_0_) 1 MPI processes</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US"> type: seqaij</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US"> rows=4851, cols=4851</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US"> total: nonzeros=139353, allocated nonzeros=139353</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US"> total number of mallocs used during MatSetValues calls=0</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US"> using I-node routines: found 3830 nodes, limit used is 5</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US"> Split number 1 Defined by IS</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US"> KSP Object: (fieldsplit_c_) 1 MPI processes</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US"> type: preonly</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US"> maximum iterations=10000, initial guess is zero</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US"> tolerances: relative=1e-05, absolute=1e-50, divergence=10000.</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US"> left preconditioning</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US"> using NONE norm type for convergence test</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US"> PC Object: (fieldsplit_c_) 1 MPI processes</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US"> type: lu</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US"> out-of-place factorization</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US"> tolerance for zero pivot 2.22045e-14</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US"> matrix ordering: nd</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US"> factor fill ratio given 5., needed 4.24323</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US"> Factored matrix follows:</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US"> Mat Object: 1 MPI processes</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US"> type: seqaij</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US"> rows=561, cols=561</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US"> package used to perform factorization: petsc</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US"> total: nonzeros=15823, allocated nonzeros=15823</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US"> not using I-node routines</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US"> linear system matrix = precond matrix:</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US"> Mat Object: (fieldsplit_c_) 1 MPI processes</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US"> type: seqaij</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US"> rows=561, cols=561</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US"> total: nonzeros=3729, allocated nonzeros=3729</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US"> total number of mallocs used during MatSetValues calls=0</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US"> not using I-node routines</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US"> linear system matrix = precond matrix:</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US"> Mat Object: 1 MPI processes</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US"> type: seqaij</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US"> rows=5412, cols=5412</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US"> total: nonzeros=190416, allocated nonzeros=190416</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US"> total number of mallocs used during MatSetValues calls=0</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US"> using I-node routines: found 3833 nodes, limit used is 5</span><u></u><u></u></p><p class="MsoNormal"><span lang="EN-US"> </span><u></u><u></u></p><p class="MsoNormal">--<u></u><u></u></p><p class="MsoNormal">Dr. Sebastian Blauth<u></u><u></u></p><p class="MsoNormal">Fraunhofer-Institut für<u></u><u></u></p><p class="MsoNormal">Techno- und Wirtschaftsmathematik ITWM<u></u><u></u></p><p class="MsoNormal">Abteilung Transportvorgänge<u></u><u></u></p><p class="MsoNormal">Fraunhofer-Platz 1, 67663 Kaiserslautern<u></u><u></u></p><p class="MsoNormal">Telefon: +49 631 31600-4968<u></u><u></u></p><p class="MsoNormal"><a href="mailto:sebastian.blauth@itwm.fraunhofer.de" target="_blank">sebastian.blauth@itwm.fraunhofer.de</a><u></u><u></u></p><p class="MsoNormal"><a href="https://urldefense.us/v3/__https://www.itwm.fraunhofer.de__;!!G_uCfscf7eWS!e9SeVBzwCGTmMl3gOc_WG4S_zL5JNMZYyiUGhfrhQAN_re34sVynQzfQyxsY8DUyFdPf1HHTM-Pas_5OTyRd$" target="_blank">https://www.itwm.fraunhofer.de</a><u></u><u></u></p><p class="MsoNormal"> <u></u><u></u></p></div></div></div></blockquote></div><p class="MsoNormal"><br clear="all"><u></u><u></u></p><div><p class="MsoNormal"><u></u> <u></u></p></div><p class="MsoNormal"><span class="m_-552471448270057200gmailsignatureprefix">-- </span><u></u><u></u></p><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><p class="MsoNormal"><u></u> <u></u></p></div><div><p class="MsoNormal"><a href="https://urldefense.us/v3/__http://www.cse.buffalo.edu/*knepley/__;fg!!G_uCfscf7eWS!e9SeVBzwCGTmMl3gOc_WG4S_zL5JNMZYyiUGhfrhQAN_re34sVynQzfQyxsY8DUyFdPf1HHTM-PasyHZNNNw$" 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></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!e9SeVBzwCGTmMl3gOc_WG4S_zL5JNMZYyiUGhfrhQAN_re34sVynQzfQyxsY8DUyFdPf1HHTM-PasyHZNNNw$" target="_blank">https://www.cse.buffalo.edu/~knepley/</a><br></div></div></div></div></div></div></div></div>