<head><!-- BaNnErBlUrFlE-HeAdEr-start -->
<style>
#pfptBanner999h97f { all: revert !important; display: block !important;
visibility: visible !important; opacity: 1 !important;
background-color: #D0D8DC !important;
max-width: none !important; max-height: none !important }
.pfptPrimaryButton999h97f:hover, .pfptPrimaryButton999h97f:focus {
background-color: #b4c1c7 !important; }
.pfptPrimaryButton999h97f: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 Fri, Jun 28, 2024 at 4: 05 AM Blauth, Sebastian <sebastian. blauth@ itwm. fraunhofer. de> wrote: Hello everyone, I have a question regarding the naming convention using PETSc’s PCFieldsplit. I have been following https: //lists. mcs. anl. gov/pipermail/petsc-users/2019-January/037262. html
</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="pfptBanner999h97f" 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="pfptBanner999h97f" style="all: unset !important; float:left !important; display:block !important; margin: 0px 0px 1px 0px !important; max-width: 600px !important;">
<div id="pfptBanner999h97f" 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="pfptBanner999h97f" 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 Fri, Jun 28, 2024 at 4:05 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-360804226658516381"><div lang="DE" style="overflow-wrap: break-word;"><div class="m_6064567489365680837WordSection1"><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!bGTaf64ibyuvBn-Qy-UQpxjLdOqFq44f6kBHzEDsXKc0htzQNw1MabtoK463uwb95Pupw_BcLMNwOHdcKldy$" 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. <u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US">I am using the following options, which work nicely for me:<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US"><u></u> <u></u></span></p><p class="MsoNormal"><span lang="EN-US">-ksp_type fgmres<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US">-pc_type fieldsplit<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US">-pc_fieldsplit_0_fields 0, 1<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US">-pc_fieldsplit_1_fields 2<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US">-pc_fieldsplit_type additive<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US">-fieldsplit_0_ksp_type fgmres<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US">-fieldsplit_0_pc_type fieldsplit<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US">-fieldsplit_0_pc_fieldsplit_type schur<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US">-fieldsplit_0_pc_fieldsplit_schur_fact_type full<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US">-fieldsplit_0_pc_fieldsplit_schur_precondition selfp<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US">-fieldsplit_0_fieldsplit_u_ksp_type preonly<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US">-fieldsplit_0_fieldsplit_u_pc_type lu<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US">-fieldsplit_0_fieldsplit_p_ksp_type cg<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US">-fieldsplit_0_fieldsplit_p_ksp_rtol 1e-14<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US">-fieldsplit_0_fieldsplit_p_ksp_atol 1e-30<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US">-fieldsplit_0_fieldsplit_p_pc_type icc<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US">-fieldsplit_0_ksp_rtol 1e-14<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US">-fieldsplit_0_ksp_atol 1e-30<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US">-fieldsplit_0_ksp_monitor_true_residual<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US">-fieldsplit_c_ksp_type preonly<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US">-fieldsplit_c_pc_type lu<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US">-ksp_view</span></p></div></div></div></blockquote><div><br></div><div>By default, we use the field names, but you can prevent this by specifying the fields by hand, so</div><div><br></div><div>-fieldsplit_0_pc_fieldsplit_0_fields 0<br>-fieldsplit_0_pc_fieldsplit_1_fields 1<br></div><div><br></div><div>should remove the 'u' and 'p' fieldnames. It is somewhat hacky, but I think easier to remember than</div><div>some extra option.</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-360804226658516381"><div lang="DE" style="overflow-wrap: break-word;"><div class="m_6064567489365680837WordSection1"><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).<u></u><u></u></span></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!bGTaf64ibyuvBn-Qy-UQpxjLdOqFq44f6kBHzEDsXKc0htzQNw1MabtoK463uwb95Pupw_BcLMNwOHdcKldy$" 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<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US"><u></u> <u></u></span></p><p class="MsoNormal"><span lang="EN-US">section.setFieldName(0, "u")<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US">section.setFieldName(1, "p")<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US">section.setFieldName(2, "c")<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US"><u></u> <u></u></span></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<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US"><u></u> <u></u></span></p><p class="MsoNormal"><span lang="EN-US">section.setFieldName(0, "0")<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US">section.setFieldName(1, "1")<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US">section.setFieldName(2, "2")<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US"><u></u> <u></u></span></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!bGTaf64ibyuvBn-Qy-UQpxjLdOqFq44f6kBHzEDsXKc0htzQNw1MabtoK463uwb95Pupw_BcLMNwOD6iRa_k$" 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. <u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US"><u></u> <u></u></span></p><p class="MsoNormal"><span lang="EN-US">Thanks a lot in advance 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"><span lang="EN-US"><u></u> <u></u></span></p><p class="MsoNormal"><span lang="EN-US">The output of ksp_view is the following:<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US">KSP Object: 1 MPI processes<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US"> type: fgmres<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US"> restart=30, using Classical (unmodified) Gram-Schmidt Orthogonalization with no iterative refinement<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US"> happy breakdown tolerance 1e-30<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US"> maximum iterations=10000, initial guess is zero<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US"> tolerances: relative=1e-05, absolute=1e-11, divergence=10000.<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US"> right preconditioning<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US"> using UNPRECONDITIONED norm type for convergence test<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US">PC Object: 1 MPI processes<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US"> type: fieldsplit<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US"> FieldSplit with ADDITIVE composition: total splits = 2<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US"> Solver info for each split is in the following KSP objects:<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US"> Split number 0 Defined by IS<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US"> KSP Object: (fieldsplit_0_) 1 MPI processes<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US"> type: fgmres<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US"> restart=30, using Classical (unmodified) Gram-Schmidt Orthogonalization with no iterative refinement<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US"> happy breakdown tolerance 1e-30<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US"> maximum iterations=10000, initial guess is zero<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US"> tolerances: relative=1e-14, absolute=1e-30, divergence=10000.<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US"> right preconditioning<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US"> using UNPRECONDITIONED norm type for convergence test<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US"> PC Object: (fieldsplit_0_) 1 MPI processes<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US"> type: fieldsplit<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US"> FieldSplit with Schur preconditioner, factorization FULL<u></u><u></u></span></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<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US"> Split info:<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US"> Split number 0 Defined by IS<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US"> Split number 1 Defined by IS<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US"> KSP solver for A00 block<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US"> KSP Object: (fieldsplit_0_fieldsplit_u_) 1 MPI processes<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US"> type: preonly<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US"> maximum iterations=10000, initial guess is zero<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US"> tolerances: relative=1e-05, absolute=1e-50, divergence=10000.<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US"> left preconditioning<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US"> using NONE norm type for convergence test<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US"> PC Object: (fieldsplit_0_fieldsplit_u_) 1 MPI processes<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US"> type: lu<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US"> out-of-place factorization<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US"> tolerance for zero pivot 2.22045e-14<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US"> matrix ordering: nd<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US"> factor fill ratio given 5., needed 3.92639<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US"> Factored matrix follows:<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US"> Mat Object: 1 MPI processes<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US"> type: seqaij<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US"> rows=4290, cols=4290<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US"> package used to perform factorization: petsc<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US"> total: nonzeros=375944, allocated nonzeros=375944<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US"> using I-node routines: found 2548 nodes, limit used is 5<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US"> linear system matrix = precond matrix:<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US"> Mat Object: (fieldsplit_0_fieldsplit_u_) 1 MPI processes<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US"> type: seqaij<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US"> rows=4290, cols=4290<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US"> total: nonzeros=95748, allocated nonzeros=95748<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US"> total number of mallocs used during MatSetValues calls=0<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US"> using I-node routines: found 3287 nodes, limit used is 5<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US"> KSP solver for S = A11 - A10 inv(A00) A01 <u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US"> KSP Object: (fieldsplit_0_fieldsplit_p_) 1 MPI processes<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US"> type: cg<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US"> maximum iterations=10000, initial guess is zero<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US"> tolerances: relative=1e-14, absolute=1e-30, divergence=10000.<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US"> left preconditioning<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US"> using PRECONDITIONED norm type for convergence test<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US"> PC Object: (fieldsplit_0_fieldsplit_p_) 1 MPI processes<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US"> type: icc<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US"> out-of-place factorization<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US"> 0 levels of fill<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US"> tolerance for zero pivot 2.22045e-14<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US"> using Manteuffel shift [POSITIVE_DEFINITE]<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US"> matrix ordering: natural<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US"> factor fill ratio given 1., needed 1.<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US"> Factored matrix follows:<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US"> Mat Object: 1 MPI processes<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US"> type: seqsbaij<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US"> rows=561, cols=561<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US"> package used to perform factorization: petsc<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US"> total: nonzeros=5120, allocated nonzeros=5120<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US"> block size is 1<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US"> linear system matrix followed by preconditioner matrix:<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US"> Mat Object: (fieldsplit_0_fieldsplit_p_) 1 MPI processes<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US"> type: schurcomplement<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US"> rows=561, cols=561<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US"> Schur complement A11 - A10 inv(A00) A01<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US"> A11<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US"> Mat Object: (fieldsplit_0_fieldsplit_p_) 1 MPI processes<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US"> type: seqaij<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US"> rows=561, cols=561<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US"> total: nonzeros=3729, allocated nonzeros=3729<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US"> total number of mallocs used during MatSetValues calls=0<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US"> not using I-node routines<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US"> A10<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US"> Mat Object: 1 MPI processes<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US"> type: seqaij<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US"> rows=561, cols=4290<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US"> total: nonzeros=19938, allocated nonzeros=19938<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US"> total number of mallocs used during MatSetValues calls=0<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US"> not using I-node routines<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US"> KSP of A00<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US"> KSP Object: (fieldsplit_0_fieldsplit_u_) 1 MPI processes<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US"> type: preonly<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US"> maximum iterations=10000, initial guess is zero<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US"> tolerances: relative=1e-05, absolute=1e-50, divergence=10000.<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US"> left preconditioning<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US"> using NONE norm type for convergence test<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US"> PC Object: (fieldsplit_0_fieldsplit_u_) 1 MPI processes<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US"> type: lu<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US"> out-of-place factorization<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US"> tolerance for zero pivot 2.22045e-14<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US"> matrix ordering: nd<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US"> factor fill ratio given 5., needed 3.92639<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US"> Factored matrix follows:<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US"> Mat Object: 1 MPI processes<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US"> type: seqaij<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US"> rows=4290, cols=4290<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US"> package used to perform factorization: petsc<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US"> total: nonzeros=375944, allocated nonzeros=375944<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US"> using I-node routines: found 2548 nodes, limit used is 5<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US"> linear system matrix = precond matrix:<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US"> Mat Object: (fieldsplit_0_fieldsplit_u_) 1 MPI processes<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US"> type: seqaij<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US"> rows=4290, cols=4290<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US"> total: nonzeros=95748, allocated nonzeros=95748<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US"> total number of mallocs used during MatSetValues calls=0<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US"> using I-node routines: found 3287 nodes, limit used is 5<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US"> A01<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US"> Mat Object: 1 MPI processes<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US"> type: seqaij<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US"> rows=4290, cols=561<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US"> total: nonzeros=19938, allocated nonzeros=19938<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US"> total number of mallocs used during MatSetValues calls=0<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US"> using I-node routines: found 3287 nodes, limit used is 5<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US"> Mat Object: 1 MPI processes<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US"> type: seqaij<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US"> rows=561, cols=561<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US"> total: nonzeros=9679, allocated nonzeros=9679<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US"> total number of mallocs used during MatSetValues calls=0<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US"> not using I-node routines<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US"> linear system matrix = precond matrix:<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US"> Mat Object: (fieldsplit_0_) 1 MPI processes<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US"> type: seqaij<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US"> rows=4851, cols=4851<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US"> total: nonzeros=139353, allocated nonzeros=139353<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US"> total number of mallocs used during MatSetValues calls=0<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US"> using I-node routines: found 3830 nodes, limit used is 5<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US"> Split number 1 Defined by IS<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US"> KSP Object: (fieldsplit_c_) 1 MPI processes<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US"> type: preonly<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US"> maximum iterations=10000, initial guess is zero<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US"> tolerances: relative=1e-05, absolute=1e-50, divergence=10000.<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US"> left preconditioning<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US"> using NONE norm type for convergence test<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US"> PC Object: (fieldsplit_c_) 1 MPI processes<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US"> type: lu<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US"> out-of-place factorization<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US"> tolerance for zero pivot 2.22045e-14<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US"> matrix ordering: nd<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US"> factor fill ratio given 5., needed 4.24323<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US"> Factored matrix follows:<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US"> Mat Object: 1 MPI processes<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US"> type: seqaij<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US"> rows=561, cols=561<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US"> package used to perform factorization: petsc<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US"> total: nonzeros=15823, allocated nonzeros=15823<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US"> not using I-node routines<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US"> linear system matrix = precond matrix:<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US"> Mat Object: (fieldsplit_c_) 1 MPI processes<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US"> type: seqaij<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US"> rows=561, cols=561<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US"> total: nonzeros=3729, allocated nonzeros=3729<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US"> total number of mallocs used during MatSetValues calls=0<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US"> not using I-node routines<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US"> linear system matrix = precond matrix:<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US"> Mat Object: 1 MPI processes<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US"> type: seqaij<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US"> rows=5412, cols=5412<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US"> total: nonzeros=190416, allocated nonzeros=190416<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US"> total number of mallocs used during MatSetValues calls=0<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US"> using I-node routines: found 3833 nodes, limit used is 5<u></u><u></u></span></p><p class="MsoNormal"><span lang="EN-US"><u></u> <u></u></span></p><p class="MsoNormal"><span>--<u></u><u></u></span></p><p class="MsoNormal"><span>Dr. Sebastian Blauth<u></u><u></u></span></p><p class="MsoNormal"><span>Fraunhofer-Institut für<u></u><u></u></span></p><p class="MsoNormal"><span>Techno- und Wirtschaftsmathematik ITWM<u></u><u></u></span></p><p class="MsoNormal"><span>Abteilung Transportvorgänge<u></u><u></u></span></p><p class="MsoNormal"><span>Fraunhofer-Platz 1, 67663 Kaiserslautern<u></u><u></u></span></p><p class="MsoNormal"><span>Telefon: +49 631 31600-4968<u></u><u></u></span></p><p class="MsoNormal"><span><a href="mailto:sebastian.blauth@itwm.fraunhofer.de" target="_blank"><span style="color:blue">sebastian.blauth@itwm.fraunhofer.de</span></a><u></u><u></u></span></p><p class="MsoNormal"><span><a href="https://urldefense.us/v3/__https://www.itwm.fraunhofer.de__;!!G_uCfscf7eWS!bGTaf64ibyuvBn-Qy-UQpxjLdOqFq44f6kBHzEDsXKc0htzQNw1MabtoK463uwb95Pupw_BcLMNwOFn7-WhZ$" target="_blank"><span style="color:blue">https://www.itwm.fraunhofer.de</span></a><u></u><u></u></span></p><p class="MsoNormal"><u></u> <u></u></p></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!bGTaf64ibyuvBn-Qy-UQpxjLdOqFq44f6kBHzEDsXKc0htzQNw1MabtoK463uwb95Pupw_BcLMNwOBtG4j6G$" target="_blank">https://www.cse.buffalo.edu/~knepley/</a><br></div></div></div></div></div></div></div></div>