<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><meta http-equiv=Content-Type content="text/html; charset=iso-8859-1"><meta name=Generator content="Microsoft Word 15 (filtered medium)"><style><!--
/* Font Definitions */
@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;}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
{margin:0cm;
font-size:11.0pt;
font-family:"Calibri",sans-serif;
mso-ligatures:standardcontextual;
mso-fareast-language:EN-US;}
a:link, span.MsoHyperlink
{mso-style-priority:99;
color:#0563C1;
text-decoration:underline;}
.MsoChpDefault
{mso-style-type:export-only;
font-family:"Calibri",sans-serif;
mso-fareast-language:EN-US;}
@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=DE link="#0563C1" vlink="#954F72" style='word-wrap:break-word'><div class=WordSection1><p class=MsoNormal>Hello everyone,<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></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://lists.mcs.anl.gov/pipermail/petsc-users/2019-January/037262.html">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. <o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US>I am using the following options, which work nicely for me:<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US><o:p> </o:p></span></p><p class=MsoNormal><span lang=EN-US>-ksp_type fgmres<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US>-pc_type fieldsplit<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US>-pc_fieldsplit_0_fields 0, 1<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US>-pc_fieldsplit_1_fields 2<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US>-pc_fieldsplit_type additive<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US>-fieldsplit_0_ksp_type fgmres<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US>-fieldsplit_0_pc_type fieldsplit<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US>-fieldsplit_0_pc_fieldsplit_type schur<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US>-fieldsplit_0_pc_fieldsplit_schur_fact_type full<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US>-fieldsplit_0_pc_fieldsplit_schur_precondition selfp<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US>-fieldsplit_0_fieldsplit_u_ksp_type preonly<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US>-fieldsplit_0_fieldsplit_u_pc_type lu<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US>-fieldsplit_0_fieldsplit_p_ksp_type cg<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US>-fieldsplit_0_fieldsplit_p_ksp_rtol 1e-14<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US>-fieldsplit_0_fieldsplit_p_ksp_atol 1e-30<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US>-fieldsplit_0_fieldsplit_p_pc_type icc<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US>-fieldsplit_0_ksp_rtol 1e-14<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US>-fieldsplit_0_ksp_atol 1e-30<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US>-fieldsplit_0_ksp_monitor_true_residual<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US>-fieldsplit_c_ksp_type preonly<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US>-fieldsplit_c_pc_type lu<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US>-ksp_view<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US><o:p> </o:p></span></p><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).<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US>Completely analogous to <a href="https://lists.mcs.anl.gov/pipermail/petsc-users/2019-January/037262.html">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<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US><o:p> </o:p></span></p><p class=MsoNormal><span lang=EN-US>section.setFieldName(0, "u")<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US>section.setFieldName(1, "p")<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US>section.setFieldName(2, "c")<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US><o:p> </o:p></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<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US><o:p> </o:p></span></p><p class=MsoNormal><span lang=EN-US>section.setFieldName(0, "0")<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US>section.setFieldName(1, "1")<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US>section.setFieldName(2, "2")<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US><o:p> </o:p></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://petsc.org/main/manualpages/PC/PCFieldSplitSetFields/">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. <o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US><o:p> </o:p></span></p><p class=MsoNormal><span lang=EN-US>Thanks a lot in advance and best regards,<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US>Sebastian<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US><o:p> </o:p></span></p><p class=MsoNormal><span lang=EN-US><o:p> </o:p></span></p><p class=MsoNormal><span lang=EN-US>The output of ksp_view is the following:<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US>KSP Object: 1 MPI processes<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US> type: fgmres<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US> restart=30, using Classical (unmodified) Gram-Schmidt Orthogonalization with no iterative refinement<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US> happy breakdown tolerance 1e-30<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US> maximum iterations=10000, initial guess is zero<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US> tolerances: relative=1e-05, absolute=1e-11, divergence=10000.<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US> right preconditioning<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US> using UNPRECONDITIONED norm type for convergence test<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US>PC Object: 1 MPI processes<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US> type: fieldsplit<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US> FieldSplit with ADDITIVE composition: total splits = 2<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US> Solver info for each split is in the following KSP objects:<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US> Split number 0 Defined by IS<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US> KSP Object: (fieldsplit_0_) 1 MPI processes<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US> type: fgmres<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US> restart=30, using Classical (unmodified) Gram-Schmidt Orthogonalization with no iterative refinement<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US> happy breakdown tolerance 1e-30<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US> maximum iterations=10000, initial guess is zero<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US> tolerances: relative=1e-14, absolute=1e-30, divergence=10000.<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US> right preconditioning<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US> using UNPRECONDITIONED norm type for convergence test<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US> PC Object: (fieldsplit_0_) 1 MPI processes<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US> type: fieldsplit<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US> FieldSplit with Schur preconditioner, factorization FULL<o:p></o:p></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<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US> Split info:<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US> Split number 0 Defined by IS<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US> Split number 1 Defined by IS<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US> KSP solver for A00 block<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US> KSP Object: (fieldsplit_0_fieldsplit_u_) 1 MPI processes<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US> type: preonly<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US> maximum iterations=10000, initial guess is zero<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US> tolerances: relative=1e-05, absolute=1e-50, divergence=10000.<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US> left preconditioning<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US> using NONE norm type for convergence test<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US> PC Object: (fieldsplit_0_fieldsplit_u_) 1 MPI processes<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US> type: lu<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US> out-of-place factorization<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US> tolerance for zero pivot 2.22045e-14<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US> matrix ordering: nd<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US> factor fill ratio given 5., needed 3.92639<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US> Factored matrix follows:<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US> Mat Object: 1 MPI processes<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US> type: seqaij<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US> rows=4290, cols=4290<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US> package used to perform factorization: petsc<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US> total: nonzeros=375944, allocated nonzeros=375944<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US> using I-node routines: found 2548 nodes, limit used is 5<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US> linear system matrix = precond matrix:<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US> Mat Object: (fieldsplit_0_fieldsplit_u_) 1 MPI processes<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US> type: seqaij<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US> rows=4290, cols=4290<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US> total: nonzeros=95748, allocated nonzeros=95748<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US> total number of mallocs used during MatSetValues calls=0<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US> using I-node routines: found 3287 nodes, limit used is 5<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US> KSP solver for S = A11 - A10 inv(A00) A01 <o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US> KSP Object: (fieldsplit_0_fieldsplit_p_) 1 MPI processes<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US> type: cg<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US> maximum iterations=10000, initial guess is zero<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US> tolerances: relative=1e-14, absolute=1e-30, divergence=10000.<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US> left preconditioning<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US> using PRECONDITIONED norm type for convergence test<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US> PC Object: (fieldsplit_0_fieldsplit_p_) 1 MPI processes<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US> type: icc<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US> out-of-place factorization<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US> 0 levels of fill<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US> tolerance for zero pivot 2.22045e-14<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US> using Manteuffel shift [POSITIVE_DEFINITE]<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US> matrix ordering: natural<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US> factor fill ratio given 1., needed 1.<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US> Factored matrix follows:<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US> Mat Object: 1 MPI processes<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US> type: seqsbaij<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US> rows=561, cols=561<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US> package used to perform factorization: petsc<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US> total: nonzeros=5120, allocated nonzeros=5120<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US> block size is 1<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US> linear system matrix followed by preconditioner matrix:<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US> Mat Object: (fieldsplit_0_fieldsplit_p_) 1 MPI processes<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US> type: schurcomplement<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US> rows=561, cols=561<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US> Schur complement A11 - A10 inv(A00) A01<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US> A11<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US> Mat Object: (fieldsplit_0_fieldsplit_p_) 1 MPI processes<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US> type: seqaij<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US> rows=561, cols=561<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US> total: nonzeros=3729, allocated nonzeros=3729<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US> total number of mallocs used during MatSetValues calls=0<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US> not using I-node routines<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US> A10<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US> Mat Object: 1 MPI processes<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US> type: seqaij<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US> rows=561, cols=4290<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US> total: nonzeros=19938, allocated nonzeros=19938<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US> total number of mallocs used during MatSetValues calls=0<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US> not using I-node routines<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US> KSP of A00<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US> KSP Object: (fieldsplit_0_fieldsplit_u_) 1 MPI processes<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US> type: preonly<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US> maximum iterations=10000, initial guess is zero<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US> tolerances: relative=1e-05, absolute=1e-50, divergence=10000.<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US> left preconditioning<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US> using NONE norm type for convergence test<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US> PC Object: (fieldsplit_0_fieldsplit_u_) 1 MPI processes<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US> type: lu<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US> out-of-place factorization<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US> tolerance for zero pivot 2.22045e-14<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US> matrix ordering: nd<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US> factor fill ratio given 5., needed 3.92639<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US> Factored matrix follows:<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US> Mat Object: 1 MPI processes<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US> type: seqaij<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US> rows=4290, cols=4290<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US> package used to perform factorization: petsc<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US> total: nonzeros=375944, allocated nonzeros=375944<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US> using I-node routines: found 2548 nodes, limit used is 5<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US> linear system matrix = precond matrix:<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US> Mat Object: (fieldsplit_0_fieldsplit_u_) 1 MPI processes<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US> type: seqaij<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US> rows=4290, cols=4290<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US> total: nonzeros=95748, allocated nonzeros=95748<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US> total number of mallocs used during MatSetValues calls=0<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US> using I-node routines: found 3287 nodes, limit used is 5<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US> A01<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US> Mat Object: 1 MPI processes<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US> type: seqaij<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US> rows=4290, cols=561<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US> total: nonzeros=19938, allocated nonzeros=19938<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US> total number of mallocs used during MatSetValues calls=0<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US> using I-node routines: found 3287 nodes, limit used is 5<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US> Mat Object: 1 MPI processes<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US> type: seqaij<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US> rows=561, cols=561<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US> total: nonzeros=9679, allocated nonzeros=9679<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US> total number of mallocs used during MatSetValues calls=0<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US> not using I-node routines<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US> linear system matrix = precond matrix:<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US> Mat Object: (fieldsplit_0_) 1 MPI processes<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US> type: seqaij<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US> rows=4851, cols=4851<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US> total: nonzeros=139353, allocated nonzeros=139353<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US> total number of mallocs used during MatSetValues calls=0<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US> using I-node routines: found 3830 nodes, limit used is 5<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US> Split number 1 Defined by IS<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US> KSP Object: (fieldsplit_c_) 1 MPI processes<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US> type: preonly<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US> maximum iterations=10000, initial guess is zero<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US> tolerances: relative=1e-05, absolute=1e-50, divergence=10000.<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US> left preconditioning<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US> using NONE norm type for convergence test<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US> PC Object: (fieldsplit_c_) 1 MPI processes<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US> type: lu<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US> out-of-place factorization<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US> tolerance for zero pivot 2.22045e-14<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US> matrix ordering: nd<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US> factor fill ratio given 5., needed 4.24323<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US> Factored matrix follows:<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US> Mat Object: 1 MPI processes<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US> type: seqaij<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US> rows=561, cols=561<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US> package used to perform factorization: petsc<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US> total: nonzeros=15823, allocated nonzeros=15823<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US> not using I-node routines<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US> linear system matrix = precond matrix:<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US> Mat Object: (fieldsplit_c_) 1 MPI processes<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US> type: seqaij<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US> rows=561, cols=561<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US> total: nonzeros=3729, allocated nonzeros=3729<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US> total number of mallocs used during MatSetValues calls=0<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US> not using I-node routines<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US> linear system matrix = precond matrix:<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US> Mat Object: 1 MPI processes<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US> type: seqaij<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US> rows=5412, cols=5412<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US> total: nonzeros=190416, allocated nonzeros=190416<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US> total number of mallocs used during MatSetValues calls=0<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US> using I-node routines: found 3833 nodes, limit used is 5<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US><o:p> </o:p></span></p><p class=MsoNormal><span style='mso-ligatures:none;mso-fareast-language:DE'>--<o:p></o:p></span></p><p class=MsoNormal><span style='mso-ligatures:none;mso-fareast-language:DE'>Dr. Sebastian Blauth<o:p></o:p></span></p><p class=MsoNormal><span style='mso-ligatures:none;mso-fareast-language:DE'>Fraunhofer-Institut für<o:p></o:p></span></p><p class=MsoNormal><span style='mso-ligatures:none;mso-fareast-language:DE'>Techno- und Wirtschaftsmathematik ITWM<o:p></o:p></span></p><p class=MsoNormal><span style='mso-ligatures:none;mso-fareast-language:DE'>Abteilung Transportvorgänge<o:p></o:p></span></p><p class=MsoNormal><span style='mso-ligatures:none;mso-fareast-language:DE'>Fraunhofer-Platz 1, 67663 Kaiserslautern<o:p></o:p></span></p><p class=MsoNormal><span style='mso-ligatures:none;mso-fareast-language:DE'>Telefon: +49 631 31600-4968<o:p></o:p></span></p><p class=MsoNormal><span style='mso-ligatures:none;mso-fareast-language:DE'><a href="mailto:sebastian.blauth@itwm.fraunhofer.de"><span style='color:blue'>sebastian.blauth@itwm.fraunhofer.de</span></a><o:p></o:p></span></p><p class=MsoNormal><span style='mso-ligatures:none;mso-fareast-language:DE'><a href="https://www.itwm.fraunhofer.de"><span style='color:blue'>https://www.itwm.fraunhofer.de</span></a><o:p></o:p></span></p><p class=MsoNormal><o:p> </o:p></p></div></body></html>