<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>