[petsc-users] Fieldsplit with LSC for constrained elasticity/poroelasticity?
Tabrez Ali
stali at geology.wisc.edu
Tue Oct 28 05:39:42 CDT 2014
Mark
When I replace "-fieldsplit_0_ksp_type preonly" with
"-fieldsplit_0_ksp_type cg" then it becomes very slow (had to kill it).
With MUMPS , i.e., with '-fieldsplit_0_pc_type lu
-fieldsplit_0_pc_factor_mat_solver_package mumps -fieldsplit_0_ksp_type
preonly' it works fine but takes more time, and will be an issue for
larger problems. The output for this run is attached.
I will work on passing rigid body modes (as Matt mentioned) but short of
that what is the best set of options for solving the following problem
(i.e., linear elasticity with constraints):
|K cG'| | u | = |F|
|G 0 | |l/c| |d|
where c is a scaling factor (so that cG' terms are more or less of the
same order as K)? The constraints are used to impose slip between
surfaces and so on.
Tabrez
On 10/27/2014 01:17 PM, Mark Adams wrote:
> The null space for GAMG is not critical but useful for elasticity. If
> you in fact have an indefinite operator (eg, not "pinned) the you
> need to use an iterative coarse grid solver. You are
> using '-fieldsplit_0_pc_type gamg -fieldsplit_0_ksp_type preonly'.
> And you have a hard elasticity problem. You are going to want to
> start with a stronger solver. Use cg instead of preonly. As Matt said
> start with MUMPS, then go to CG/GAMG, then you can see how far you can
> cut the _0_ solver down.
> Mark
>
> On Thu, Oct 23, 2014 at 11:51 AM, Matthew Knepley <knepley at gmail.com
> <mailto:knepley at gmail.com>> wrote:
>
> On Thu, Oct 23, 2014 at 10:48 AM, Tabrez Ali
> <stali at geology.wisc.edu <mailto:stali at geology.wisc.edu>> wrote:
>
> Matt
>
> On 10/23/2014 09:54 AM, Matthew Knepley wrote:
>> On Thu, Oct 23, 2014 at 9:27 AM, Tabrez Ali
>> <stali at geology.wisc.edu <mailto:stali at geology.wisc.edu>> wrote:
>>
>> Matt
>>
>> Sorry about that (I always forget it). The output for the
>> smallest problem is now attached (see log.txt). I am also
>> attaching some results that compare results obtained
>> using FS/LSC and the direct solver (MUMPS), again for the
>> smallest problem. The difference, as you can see is
>> insignificant O(1E-6).
>>
>>
>> 1) How do you use MUMPS if you have a saddle point
> I simply used -pc_type lu -pc_factor_mat_solver_package mumps.
>
>>
>> 2) You can see from the output that something is seriously
>> wrong with the preconditioner. It looks like it has a null space.
>> Did you add the elastic null modes to GAMG? Without this,
>> it is not going to work. We have helper functions for this:
>>
>> http://www.mcs.anl.gov/petsc/petsc-current/docs/manualpages/DM/DMPlexCreateRigidBody.html
>>
>> you could just copy that code. And then use
>>
>> http://www.mcs.anl.gov/petsc/petsc-current/docs/manualpages/Mat/MatSetNearNullSpace.html
>>
>> I don't see it in the output, so I think this is your problem.
>>
>> In order to test, I would first use MUMPS as the A00 solver
>> and get the Schur stuff worked out. Then I would
>> replace MUMPS with GAMG and tune it until I get back my
>> original convergence.
> I will try this with MatNullSpaceCreateRigidBody. Btw does it
> matter if some nodes are pinned?
>
>
> No these are null modes of the operator, not of the particular
> problem.
>
> Matt
>
> Tabrez
>
>>
>> Thanks,
>>
>> Matt
>>
>> Also, I did pass 'upper' and 'full' to
>> '-pc_fieldsplit_schur_factorization_type' but the
>> iteration count doesn't improve (in fact, it increases
>> slightly). The attached log is with 'upper'.
>>
>> Regards,
>>
>> Tabrez
>>
>> On 10/23/2014 07:46 AM, Matthew Knepley wrote:
>>> On Thu, Oct 23, 2014 at 7:20 AM, Tabrez Ali
>>> <stali at geology.wisc.edu <mailto:stali at geology.wisc.edu>>
>>> wrote:
>>>
>>> Hello
>>>
>>> I am using the following options (below) for solving
>>> linear elasticity/poroelasticity problems involving
>>> slip between two surfaces involving non-trivial
>>> geometries, i.e., elements with high aspect ratios,
>>> large contrasts in material properties etc. The
>>> constraints are imposed using Lagrange Multipliers.
>>>
>>> A picture (shows displacement magnitude) is
>>> attached. The boundary nodes, i.e., the base and the
>>> four side are pinned.
>>>
>>> The following options appear to work well for the
>>> saddle point problem:
>>>
>>> -pc_type fieldsplit -pc_fieldsplit_type schur
>>> -pc_fieldsplit_detect_saddle_point
>>> -fieldsplit_0_pc_type gamg -fieldsplit_0_ksp_type
>>> preonly -fieldsplit_1_pc_type lsc
>>> -fieldsplit_1_ksp_type preonly
>>> -pc_fieldsplit_schur_fact_type lower -ksp_monitor
>>>
>>> However, the number of iterations keep on increasing
>>> with the problems size (see attached plot), e.g.,
>>>
>>> 120K Tets *507* Iterations (KSP Residual norm
>>> 8.827362494659e-05)in 17 secs on 3 cores
>>> 1 Million Tets *1374* Iterations (KSP Residual norm
>>> 7.164704416296e-05)in 117 secs on 20 cores
>>> 8 Million Tets *2495* Iterations (KSP Residual norm
>>> 9.101247550026e-05) in 225 secs on 160 cores
>>>
>>> So what other options should I try to improve solver
>>> performance? Any tips/insights would be appreciated
>>> as preconditioning is black magic to me.
>>>
>>>
>>> For reports, always run with
>>>
>>> -ksp_view -ksp_monitor_true_residual -ksp_converged_reason
>>>
>>> so that we can see exactly what you used.
>>>
>>> I believe the default is a diagonal factorization. Since
>>> your outer iterates are increasing, I would strengthen this
>>> to either upper or full
>>>
>>> -pc_fieldsplit_schur_factorization_type <upper, full>
>>>
>>> Thanks,
>>>
>>> Matt
>>>
>>> Thanks in advance.
>>>
>>> Tabrez
>>>
>>>
>>>
>>>
>>> --
>>> What most experimenters take for granted before they
>>> begin their experiments is infinitely more interesting
>>> than any results to which their experiments lead.
>>> -- Norbert Wiener
>>
>>
>>
>>
>> --
>> What most experimenters take for granted before they begin
>> their experiments is infinitely more interesting than any
>> results to which their experiments lead.
>> -- Norbert Wiener
>
>
>
>
> --
> What most experimenters take for granted before they begin their
> experiments is infinitely more interesting than any results to
> which their experiments lead.
> -- Norbert Wiener
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20141028/4c52a78e/attachment-0001.html>
-------------- next part --------------
salloc: Granted job allocation 78969
Reading input ...
Partitioning mesh ...
Reading mesh data ...
Forming [K] ...
Applying constraints ...
Forming RHS ...
Setting up solver ...
Solving ...
0 KSP preconditioned resid norm 9.320055451716e+05 true resid norm 1.755998647494e+02 ||r(i)||/||b|| 1.000000000000e+00
1 KSP preconditioned resid norm 9.127510702346e+05 true resid norm 1.362125072634e+14 ||r(i)||/||b|| 7.756982470220e+11
2 KSP preconditioned resid norm 8.626930780275e+05 true resid norm 7.140706098955e+14 ||r(i)||/||b|| 4.066464464051e+12
3 KSP preconditioned resid norm 6.917392651646e+05 true resid norm 2.788588220631e+15 ||r(i)||/||b|| 1.588035517346e+13
4 KSP preconditioned resid norm 5.948180320596e+05 true resid norm 3.828406259238e+15 ||r(i)||/||b|| 2.180187476056e+13
5 KSP preconditioned resid norm 5.365265868551e+05 true resid norm 4.509015181942e+15 ||r(i)||/||b|| 2.567778277265e+13
6 KSP preconditioned resid norm 5.020807291688e+05 true resid norm 4.938429989057e+15 ||r(i)||/||b|| 2.812319927527e+13
7 KSP preconditioned resid norm 4.691107247107e+05 true resid norm 5.504613096315e+15 ||r(i)||/||b|| 3.134747913486e+13
8 KSP preconditioned resid norm 4.500805942145e+05 true resid norm 5.950219821547e+15 ||r(i)||/||b|| 3.388510480939e+13
9 KSP preconditioned resid norm 4.250032501744e+05 true resid norm 6.670720166034e+15 ||r(i)||/||b|| 3.798818510228e+13
10 KSP preconditioned resid norm 4.051645463810e+05 true resid norm 7.414211051141e+15 ||r(i)||/||b|| 4.222219112596e+13
11 KSP preconditioned resid norm 3.880013874032e+05 true resid norm 8.151575238507e+15 ||r(i)||/||b|| 4.642130704452e+13
12 KSP preconditioned resid norm 3.741408602796e+05 true resid norm 8.948407390773e+15 ||r(i)||/||b|| 5.095907906048e+13
13 KSP preconditioned resid norm 3.532148969165e+05 true resid norm 1.028219119134e+16 ||r(i)||/||b|| 5.855466464062e+13
14 KSP preconditioned resid norm 3.306246816868e+05 true resid norm 1.207912805635e+16 ||r(i)||/||b|| 6.878779817734e+13
15 KSP preconditioned resid norm 3.077008490960e+05 true resid norm 1.412194821609e+16 ||r(i)||/||b|| 8.042117934570e+13
16 KSP preconditioned resid norm 2.709436838755e+05 true resid norm 1.756809508461e+16 ||r(i)||/||b|| 1.000461766282e+14
17 KSP preconditioned resid norm 2.491233804162e+05 true resid norm 1.976045004357e+16 ||r(i)||/||b|| 1.125311233683e+14
18 KSP preconditioned resid norm 2.169453185969e+05 true resid norm 2.315544602061e+16 ||r(i)||/||b|| 1.318648283338e+14
19 KSP preconditioned resid norm 1.853345270881e+05 true resid norm 2.647069144539e+16 ||r(i)||/||b|| 1.507443726291e+14
20 KSP preconditioned resid norm 1.684382143134e+05 true resid norm 2.847996525788e+16 ||r(i)||/||b|| 1.621867152263e+14
21 KSP preconditioned resid norm 1.487793180708e+05 true resid norm 3.106944054176e+16 ||r(i)||/||b|| 1.769331689754e+14
22 KSP preconditioned resid norm 1.324108267213e+05 true resid norm 3.364687707291e+16 ||r(i)||/||b|| 1.916110648544e+14
23 KSP preconditioned resid norm 1.154783368235e+05 true resid norm 3.664744963424e+16 ||r(i)||/||b|| 2.086986210755e+14
24 KSP preconditioned resid norm 9.690290840840e+04 true resid norm 4.059423491278e+16 ||r(i)||/||b|| 2.311746365563e+14
25 KSP preconditioned resid norm 8.541667474325e+04 true resid norm 4.382080823729e+16 ||r(i)||/||b|| 2.495492140602e+14
26 KSP preconditioned resid norm 6.738137388808e+04 true resid norm 4.935828200762e+16 ||r(i)||/||b|| 2.810838270181e+14
27 KSP preconditioned resid norm 5.559761934660e+04 true resid norm 5.372591255719e+16 ||r(i)||/||b|| 3.059564575057e+14
28 KSP preconditioned resid norm 4.606357787962e+04 true resid norm 5.783519870706e+16 ||r(i)||/||b|| 3.293578772945e+14
29 KSP preconditioned resid norm 3.846753020730e+04 true resid norm 6.227329505424e+16 ||r(i)||/||b|| 3.546317939545e+14
30 KSP preconditioned resid norm 3.331905232453e+04 true resid norm 6.606924388231e+16 ||r(i)||/||b|| 3.762488312653e+14
31 KSP preconditioned resid norm 3.290201988582e+04 true resid norm 6.618297473845e+16 ||r(i)||/||b|| 3.768965017878e+14
32 KSP preconditioned resid norm 3.162870028104e+04 true resid norm 6.654393352106e+16 ||r(i)||/||b|| 3.789520773039e+14
33 KSP preconditioned resid norm 2.991966449152e+04 true resid norm 6.709255863611e+16 ||r(i)||/||b|| 3.820763685204e+14
34 KSP preconditioned resid norm 2.922952746973e+04 true resid norm 6.743695964917e+16 ||r(i)||/||b|| 3.840376514265e+14
35 KSP preconditioned resid norm 2.798614888145e+04 true resid norm 6.838110987651e+16 ||r(i)||/||b|| 3.894143652907e+14
36 KSP preconditioned resid norm 2.701665220461e+04 true resid norm 6.929642686679e+16 ||r(i)||/||b|| 3.946268806396e+14
37 KSP preconditioned resid norm 2.601811428056e+04 true resid norm 7.044504608577e+16 ||r(i)||/||b|| 4.011679974032e+14
38 KSP preconditioned resid norm 2.531056406767e+04 true resid norm 7.153872989114e+16 ||r(i)||/||b|| 4.073962698846e+14
39 KSP preconditioned resid norm 2.332617529703e+04 true resid norm 7.497512777257e+16 ||r(i)||/||b|| 4.269657489747e+14
40 KSP preconditioned resid norm 2.254584604090e+04 true resid norm 7.640212636531e+16 ||r(i)||/||b|| 4.350921709100e+14
41 KSP preconditioned resid norm 2.082879229685e+04 true resid norm 7.986664049248e+16 ||r(i)||/||b|| 4.548217654180e+14
42 KSP preconditioned resid norm 1.914379086156e+04 true resid norm 8.367676210994e+16 ||r(i)||/||b|| 4.765195134368e+14
43 KSP preconditioned resid norm 1.815822267369e+04 true resid norm 8.615765209492e+16 ||r(i)||/||b|| 4.906475994038e+14
44 KSP preconditioned resid norm 1.677428596939e+04 true resid norm 8.993611308328e+16 ||r(i)||/||b|| 5.121650475736e+14
45 KSP preconditioned resid norm 1.534247135375e+04 true resid norm 9.414497408723e+16 ||r(i)||/||b|| 5.361335227769e+14
46 KSP preconditioned resid norm 1.412608783430e+04 true resid norm 9.852630241654e+16 ||r(i)||/||b|| 5.610841589039e+14
47 KSP preconditioned resid norm 1.177668978657e+04 true resid norm 1.077880272442e+17 ||r(i)||/||b|| 6.138275072022e+14
48 KSP preconditioned resid norm 1.065678918714e+04 true resid norm 1.124882415134e+17 ||r(i)||/||b|| 6.405941238848e+14
49 KSP preconditioned resid norm 9.962410636884e+03 true resid norm 1.158017305255e+17 ||r(i)||/||b|| 6.594636658224e+14
50 KSP preconditioned resid norm 9.371192462863e+03 true resid norm 1.195573007708e+17 ||r(i)||/||b|| 6.808507565847e+14
51 KSP preconditioned resid norm 8.608163058470e+03 true resid norm 1.259685427089e+17 ||r(i)||/||b|| 7.173612741030e+14
52 KSP preconditioned resid norm 7.720924959940e+03 true resid norm 1.350861356613e+17 ||r(i)||/||b|| 7.692838252129e+14
53 KSP preconditioned resid norm 6.628413884374e+03 true resid norm 1.478135305653e+17 ||r(i)||/||b|| 8.417633508790e+14
54 KSP preconditioned resid norm 5.583298606838e+03 true resid norm 1.607393601806e+17 ||r(i)||/||b|| 9.153729156340e+14
55 KSP preconditioned resid norm 4.854896331327e+03 true resid norm 1.721936891600e+17 ||r(i)||/||b|| 9.806026297673e+14
56 KSP preconditioned resid norm 4.262399411003e+03 true resid norm 1.837125351729e+17 ||r(i)||/||b|| 1.046199753258e+15
57 KSP preconditioned resid norm 3.584024669247e+03 true resid norm 2.026508285264e+17 ||r(i)||/||b|| 1.154048887313e+15
58 KSP preconditioned resid norm 3.302611896991e+03 true resid norm 2.155979885039e+17 ||r(i)||/||b|| 1.227779923473e+15
59 KSP preconditioned resid norm 3.106452236021e+03 true resid norm 2.294151324082e+17 ||r(i)||/||b|| 1.306465313830e+15
60 KSP preconditioned resid norm 2.894632302133e+03 true resid norm 2.500169543959e+17 ||r(i)||/||b|| 1.423787852871e+15
61 KSP preconditioned resid norm 2.892427608148e+03 true resid norm 2.503572781245e+17 ||r(i)||/||b|| 1.425725916599e+15
62 KSP preconditioned resid norm 2.881855812307e+03 true resid norm 2.520690203632e+17 ||r(i)||/||b|| 1.435473886742e+15
63 KSP preconditioned resid norm 2.872986380486e+03 true resid norm 2.536302493244e+17 ||r(i)||/||b|| 1.444364719110e+15
64 KSP preconditioned resid norm 2.865635949494e+03 true resid norm 2.551732099856e+17 ||r(i)||/||b|| 1.453151517797e+15
65 KSP preconditioned resid norm 2.860836435783e+03 true resid norm 2.565769302199e+17 ||r(i)||/||b|| 1.461145374947e+15
66 KSP preconditioned resid norm 2.860026342234e+03 true resid norm 2.574294644148e+17 ||r(i)||/||b|| 1.466000357017e+15
67 KSP preconditioned resid norm 2.860026336219e+03 true resid norm 2.574330960365e+17 ||r(i)||/||b|| 1.466021038250e+15
68 KSP preconditioned resid norm 2.859446352247e+03 true resid norm 2.564102967021e+17 ||r(i)||/||b|| 1.460196436188e+15
69 KSP preconditioned resid norm 2.855291192171e+03 true resid norm 2.530743143761e+17 ||r(i)||/||b|| 1.441198800109e+15
70 KSP preconditioned resid norm 2.843683581105e+03 true resid norm 2.467887076760e+17 ||r(i)||/||b|| 1.405403745773e+15
71 KSP preconditioned resid norm 2.825924372436e+03 true resid norm 2.382078643091e+17 ||r(i)||/||b|| 1.356537857527e+15
72 KSP preconditioned resid norm 2.803340022921e+03 true resid norm 2.288941643098e+17 ||r(i)||/||b|| 1.303498522829e+15
73 KSP preconditioned resid norm 2.768188650590e+03 true resid norm 2.170599057122e+17 ||r(i)||/||b|| 1.236105198726e+15
74 KSP preconditioned resid norm 2.728040723580e+03 true resid norm 2.050447447412e+17 ||r(i)||/||b|| 1.167681678080e+15
75 KSP preconditioned resid norm 2.663723517309e+03 true resid norm 1.886813389061e+17 ||r(i)||/||b|| 1.074495923874e+15
76 KSP preconditioned resid norm 2.580490699705e+03 true resid norm 1.697855415273e+17 ||r(i)||/||b|| 9.668887944168e+14
77 KSP preconditioned resid norm 2.507785655255e+03 true resid norm 1.549544648278e+17 ||r(i)||/||b|| 8.824292948570e+14
78 KSP preconditioned resid norm 2.381407885460e+03 true resid norm 1.327524813313e+17 ||r(i)||/||b|| 7.559942117309e+14
79 KSP preconditioned resid norm 2.275211715962e+03 true resid norm 1.165811444580e+17 ||r(i)||/||b|| 6.639022451665e+14
80 KSP preconditioned resid norm 2.161471490944e+03 true resid norm 1.015723094108e+17 ||r(i)||/||b|| 5.784304535530e+14
81 KSP preconditioned resid norm 1.982790705425e+03 true resid norm 8.141089236173e+16 ||r(i)||/||b|| 4.636159172328e+14
82 KSP preconditioned resid norm 1.740676429181e+03 true resid norm 5.840122410223e+16 ||r(i)||/||b|| 3.325812590208e+14
83 KSP preconditioned resid norm 1.587150102469e+03 true resid norm 4.629492748351e+16 ||r(i)||/||b|| 2.636387422598e+14
84 KSP preconditioned resid norm 1.322911469852e+03 true resid norm 2.928137569559e+16 ||r(i)||/||b|| 1.667505595029e+14
85 KSP preconditioned resid norm 1.181708294070e+03 true resid norm 2.197349343473e+16 ||r(i)||/||b|| 1.251338858722e+14
86 KSP preconditioned resid norm 9.902319656225e+02 true resid norm 1.423216418217e+16 ||r(i)||/||b|| 8.104883339450e+13
87 KSP preconditioned resid norm 8.359003992039e+02 true resid norm 9.303072101775e+15 ||r(i)||/||b|| 5.297881131659e+13
88 KSP preconditioned resid norm 6.981353830695e+02 true resid norm 5.935654404148e+15 ||r(i)||/||b|| 3.380215817717e+13
89 KSP preconditioned resid norm 5.800715479524e+02 true resid norm 3.708718763345e+15 ||r(i)||/||b|| 2.112028257333e+13
90 KSP preconditioned resid norm 4.883815028666e+02 true resid norm 2.388206870627e+15 ||r(i)||/||b|| 1.360027739222e+13
91 KSP preconditioned resid norm 4.789803991475e+02 true resid norm 2.397467761477e+15 ||r(i)||/||b|| 1.365301599121e+13
92 KSP preconditioned resid norm 4.571508470005e+02 true resid norm 2.421108357131e+15 ||r(i)||/||b|| 1.378764363279e+13
93 KSP preconditioned resid norm 4.504410932746e+02 true resid norm 2.430412703020e+15 ||r(i)||/||b|| 1.384062969803e+13
94 KSP preconditioned resid norm 4.405021383131e+02 true resid norm 2.456970033422e+15 ||r(i)||/||b|| 1.399186745917e+13
95 KSP preconditioned resid norm 4.293421518077e+02 true resid norm 2.491683359417e+15 ||r(i)||/||b|| 1.418955170025e+13
96 KSP preconditioned resid norm 4.173639525791e+02 true resid norm 2.540094021569e+15 ||r(i)||/||b|| 1.446523905468e+13
97 KSP preconditioned resid norm 4.071406525465e+02 true resid norm 2.585692506997e+15 ||r(i)||/||b|| 1.472491172295e+13
98 KSP preconditioned resid norm 3.923770601661e+02 true resid norm 2.664481005387e+15 ||r(i)||/||b|| 1.517359372224e+13
99 KSP preconditioned resid norm 3.696685994238e+02 true resid norm 2.790356355448e+15 ||r(i)||/||b|| 1.589042428609e+13
100 KSP preconditioned resid norm 3.522302545921e+02 true resid norm 2.891433818038e+15 ||r(i)||/||b|| 1.646603670319e+13
101 KSP preconditioned resid norm 3.289451418467e+02 true resid norm 3.043740223578e+15 ||r(i)||/||b|| 1.733338592215e+13
102 KSP preconditioned resid norm 3.041026598143e+02 true resid norm 3.223224072832e+15 ||r(i)||/||b|| 1.835550430197e+13
103 KSP preconditioned resid norm 2.721797727812e+02 true resid norm 3.434699332293e+15 ||r(i)||/||b|| 1.955980625153e+13
104 KSP preconditioned resid norm 2.495065013936e+02 true resid norm 3.595571585792e+15 ||r(i)||/||b|| 2.047593596341e+13
105 KSP preconditioned resid norm 2.226326085921e+02 true resid norm 3.801867914908e+15 ||r(i)||/||b|| 2.165074512064e+13
106 KSP preconditioned resid norm 2.073865356308e+02 true resid norm 3.926005989758e+15 ||r(i)||/||b|| 2.235768231007e+13
107 KSP preconditioned resid norm 1.889509629828e+02 true resid norm 4.115952048686e+15 ||r(i)||/||b|| 2.343938051752e+13
108 KSP preconditioned resid norm 1.749949837765e+02 true resid norm 4.278800226377e+15 ||r(i)||/||b|| 2.436676265374e+13
109 KSP preconditioned resid norm 1.618193234597e+02 true resid norm 4.465831577902e+15 ||r(i)||/||b|| 2.543186228688e+13
110 KSP preconditioned resid norm 1.465181121950e+02 true resid norm 4.714656055446e+15 ||r(i)||/||b|| 2.684885926407e+13
111 KSP preconditioned resid norm 1.318107867392e+02 true resid norm 5.037611920944e+15 ||r(i)||/||b|| 2.868801709007e+13
112 KSP preconditioned resid norm 1.242084520074e+02 true resid norm 5.234477391902e+15 ||r(i)||/||b|| 2.980911972440e+13
113 KSP preconditioned resid norm 1.155780632405e+02 true resid norm 5.582223270103e+15 ||r(i)||/||b|| 3.178945085220e+13
114 KSP preconditioned resid norm 1.021440284034e+02 true resid norm 6.169174994158e+15 ||r(i)||/||b|| 3.513200310810e+13
115 KSP preconditioned resid norm 9.495516699860e+01 true resid norm 6.557142438207e+15 ||r(i)||/||b|| 3.734138660964e+13
116 KSP preconditioned resid norm 8.974660320529e+01 true resid norm 6.966234434027e+15 ||r(i)||/||b|| 3.967106947360e+13
117 KSP preconditioned resid norm 8.668620937542e+01 true resid norm 7.370496578775e+15 ||r(i)||/||b|| 4.197324746972e+13
118 KSP preconditioned resid norm 8.544499762555e+01 true resid norm 7.721074005539e+15 ||r(i)||/||b|| 4.396970360176e+13
119 KSP preconditioned resid norm 8.540357047452e+01 true resid norm 7.783986871592e+15 ||r(i)||/||b|| 4.432797760238e+13
120 KSP preconditioned resid norm 8.509182838139e+01 true resid norm 7.567150037172e+15 ||r(i)||/||b|| 4.309314274229e+13
121 KSP preconditioned resid norm 8.501310700513e+01 true resid norm 7.594940200282e+15 ||r(i)||/||b|| 4.325140119625e+13
122 KSP preconditioned resid norm 8.493474829240e+01 true resid norm 7.626607956863e+15 ||r(i)||/||b|| 4.343174163458e+13
123 KSP preconditioned resid norm 8.488333741555e+01 true resid norm 7.656182294989e+15 ||r(i)||/||b|| 4.360016054635e+13
124 KSP preconditioned resid norm 8.487226574331e+01 true resid norm 7.669732385858e+15 ||r(i)||/||b|| 4.367732513237e+13
125 KSP preconditioned resid norm 8.487223372854e+01 true resid norm 7.668341900577e+15 ||r(i)||/||b|| 4.366940664517e+13
126 KSP preconditioned resid norm 8.485866512717e+01 true resid norm 7.639207309985e+15 ||r(i)||/||b|| 4.350349199236e+13
127 KSP preconditioned resid norm 8.480919141670e+01 true resid norm 7.584336229378e+15 ||r(i)||/||b|| 4.319101407169e+13
128 KSP preconditioned resid norm 8.461824281885e+01 true resid norm 7.465571488217e+15 ||r(i)||/||b|| 4.251467675600e+13
129 KSP preconditioned resid norm 8.425430504400e+01 true resid norm 7.263958891880e+15 ||r(i)||/||b|| 4.136654035724e+13
130 KSP preconditioned resid norm 8.380956322835e+01 true resid norm 7.064381851222e+15 ||r(i)||/||b|| 4.022999597012e+13
131 KSP preconditioned resid norm 8.314252321333e+01 true resid norm 6.808447921810e+15 ||r(i)||/||b|| 3.877251233381e+13
132 KSP preconditioned resid norm 8.226839402109e+01 true resid norm 6.516419146749e+15 ||r(i)||/||b|| 3.710947702636e+13
133 KSP preconditioned resid norm 8.078979668705e+01 true resid norm 6.095156115716e+15 ||r(i)||/||b|| 3.471048297455e+13
134 KSP preconditioned resid norm 7.933264464976e+01 true resid norm 5.736990471656e+15 ||r(i)||/||b|| 3.267081372667e+13
135 KSP preconditioned resid norm 7.629871013849e+01 true resid norm 5.098475273910e+15 ||r(i)||/||b|| 2.903461959488e+13
136 KSP preconditioned resid norm 7.384550751976e+01 true resid norm 4.620263860149e+15 ||r(i)||/||b|| 2.631131787455e+13
137 KSP preconditioned resid norm 7.019186335982e+01 true resid norm 3.995768284315e+15 ||r(i)||/||b|| 2.275496219783e+13
138 KSP preconditioned resid norm 6.549021374298e+01 true resid norm 3.271622487708e+15 ||r(i)||/||b|| 1.863112191104e+13
139 KSP preconditioned resid norm 6.216358544968e+01 true resid norm 2.825089859247e+15 ||r(i)||/||b|| 1.608822343502e+13
140 KSP preconditioned resid norm 5.819590402733e+01 true resid norm 2.365235725586e+15 ||r(i)||/||b|| 1.346946211469e+13
141 KSP preconditioned resid norm 5.322146614464e+01 true resid norm 1.862541582703e+15 ||r(i)||/||b|| 1.060673700040e+13
142 KSP preconditioned resid norm 4.838864093542e+01 true resid norm 1.460907819000e+15 ||r(i)||/||b|| 8.319527017206e+12
143 KSP preconditioned resid norm 4.288902076387e+01 true resid norm 1.075226017243e+15 ||r(i)||/||b|| 6.123159711867e+12
144 KSP preconditioned resid norm 3.774714312484e+01 true resid norm 7.813012695376e+14 ||r(i)||/||b|| 4.449327285375e+12
145 KSP preconditioned resid norm 3.210844216241e+01 true resid norm 5.182913258583e+14 ||r(i)||/||b|| 2.951547409208e+12
146 KSP preconditioned resid norm 2.633903472489e+01 true resid norm 3.154832514028e+14 ||r(i)||/||b|| 1.796603043249e+12
147 KSP preconditioned resid norm 2.195326834454e+01 true resid norm 1.974701889296e+14 ||r(i)||/||b|| 1.124546361191e+12
148 KSP preconditioned resid norm 1.932058797915e+01 true resid norm 1.401755026989e+14 ||r(i)||/||b|| 7.982665755408e+11
149 KSP preconditioned resid norm 1.720408988995e+01 true resid norm 1.028932403836e+14 ||r(i)||/||b|| 5.859528452969e+11
150 KSP preconditioned resid norm 1.559612800537e+01 true resid norm 8.129959387556e+13 ||r(i)||/||b|| 4.629820984862e+11
151 KSP preconditioned resid norm 1.528043313778e+01 true resid norm 8.168488419825e+13 ||r(i)||/||b|| 4.651762364102e+11
152 KSP preconditioned resid norm 1.473378998995e+01 true resid norm 8.241472575976e+13 ||r(i)||/||b|| 4.693325127407e+11
153 KSP preconditioned resid norm 1.446554996245e+01 true resid norm 8.283307200666e+13 ||r(i)||/||b|| 4.717148963917e+11
154 KSP preconditioned resid norm 1.433069476310e+01 true resid norm 8.320633972939e+13 ||r(i)||/||b|| 4.738405684318e+11
155 KSP preconditioned resid norm 1.398451011406e+01 true resid norm 8.447592357675e+13 ||r(i)||/||b|| 4.810705503521e+11
156 KSP preconditioned resid norm 1.338238273977e+01 true resid norm 8.694936201334e+13 ||r(i)||/||b|| 4.951562015007e+11
157 KSP preconditioned resid norm 1.285596644400e+01 true resid norm 8.942599864914e+13 ||r(i)||/||b|| 5.092600656428e+11
158 KSP preconditioned resid norm 1.259640614474e+01 true resid norm 9.096283789856e+13 ||r(i)||/||b|| 5.180120043281e+11
159 KSP preconditioned resid norm 1.150402605336e+01 true resid norm 9.770715990572e+13 ||r(i)||/||b|| 5.564193346343e+11
160 KSP preconditioned resid norm 1.091253930358e+01 true resid norm 1.013341154326e+14 ||r(i)||/||b|| 5.770739947734e+11
161 KSP preconditioned resid norm 1.005190049488e+01 true resid norm 1.071270024087e+14 ||r(i)||/||b|| 6.100631259688e+11
162 KSP preconditioned resid norm 9.249803666300e+00 true resid norm 1.126347826847e+14 ||r(i)||/||b|| 6.414286414481e+11
163 KSP preconditioned resid norm 8.498257679037e+00 true resid norm 1.184209546140e+14 ||r(i)||/||b|| 6.743795320290e+11
164 KSP preconditioned resid norm 7.683747445800e+00 true resid norm 1.250523348873e+14 ||r(i)||/||b|| 7.121436856789e+11
165 KSP preconditioned resid norm 7.203608831235e+00 true resid norm 1.292064073553e+14 ||r(i)||/||b|| 7.358001530340e+11
166 KSP preconditioned resid norm 6.554687342538e+00 true resid norm 1.355010457000e+14 ||r(i)||/||b|| 7.716466404650e+11
167 KSP preconditioned resid norm 6.040626373592e+00 true resid norm 1.410685680534e+14 ||r(i)||/||b|| 8.033523730485e+11
168 KSP preconditioned resid norm 5.602778891604e+00 true resid norm 1.472337491989e+14 ||r(i)||/||b|| 8.384616321263e+11
169 KSP preconditioned resid norm 5.137313672639e+00 true resid norm 1.546288196744e+14 ||r(i)||/||b|| 8.805748221682e+11
170 KSP preconditioned resid norm 4.772073867658e+00 true resid norm 1.617510890095e+14 ||r(i)||/||b|| 9.211344737672e+11
171 KSP preconditioned resid norm 4.429551648207e+00 true resid norm 1.699873935479e+14 ||r(i)||/||b|| 9.680382942806e+11
172 KSP preconditioned resid norm 4.058041912274e+00 true resid norm 1.807915228232e+14 ||r(i)||/||b|| 1.029565273761e+12
173 KSP preconditioned resid norm 3.659353082764e+00 true resid norm 1.979585093233e+14 ||r(i)||/||b|| 1.127327231179e+12
174 KSP preconditioned resid norm 3.345837180280e+00 true resid norm 2.142616541693e+14 ||r(i)||/||b|| 1.220169813201e+12
175 KSP preconditioned resid norm 3.144285418302e+00 true resid norm 2.281976140099e+14 ||r(i)||/||b|| 1.299531832417e+12
176 KSP preconditioned resid norm 3.022396101784e+00 true resid norm 2.394374033937e+14 ||r(i)||/||b|| 1.363539793925e+12
177 KSP preconditioned resid norm 2.937050870466e+00 true resid norm 2.511983521409e+14 ||r(i)||/||b|| 1.430515635644e+12
178 KSP preconditioned resid norm 2.894843034244e+00 true resid norm 2.630993320056e+14 ||r(i)||/||b|| 1.498288921698e+12
179 KSP preconditioned resid norm 2.891858751192e+00 true resid norm 2.664327984692e+14 ||r(i)||/||b|| 1.517272230531e+12
180 KSP preconditioned resid norm 2.881378316904e+00 true resid norm 2.599860822183e+14 ||r(i)||/||b|| 1.480559695131e+12
181 KSP preconditioned resid norm 2.879541972597e+00 true resid norm 2.605905686926e+14 ||r(i)||/||b|| 1.484002103672e+12
182 KSP preconditioned resid norm 2.877331259337e+00 true resid norm 2.614134463899e+14 ||r(i)||/||b|| 1.488688198951e+12
183 KSP preconditioned resid norm 2.874456654944e+00 true resid norm 2.627919787038e+14 ||r(i)||/||b|| 1.496538616809e+12
184 KSP preconditioned resid norm 2.874349383523e+00 true resid norm 2.629889220630e+14 ||r(i)||/||b|| 1.497660162998e+12
185 KSP preconditioned resid norm 2.874315970654e+00 true resid norm 2.632240084014e+14 ||r(i)||/||b|| 1.498998924499e+12
186 KSP preconditioned resid norm 2.873426488098e+00 true resid norm 2.616629867267e+14 ||r(i)||/||b|| 1.490109272579e+12
187 KSP preconditioned resid norm 2.871828875691e+00 true resid norm 2.600477959076e+14 ||r(i)||/||b|| 1.480911140101e+12
188 KSP preconditioned resid norm 2.866311172083e+00 true resid norm 2.569199084288e+14 ||r(i)||/||b|| 1.463098555318e+12
189 KSP preconditioned resid norm 2.846490699873e+00 true resid norm 2.470737050314e+14 ||r(i)||/||b|| 1.407026738796e+12
190 KSP preconditioned resid norm 2.819730703425e+00 true resid norm 2.362845722118e+14 ||r(i)||/||b|| 1.345585160610e+12
191 KSP preconditioned resid norm 2.792935326546e+00 true resid norm 2.273870208563e+14 ||r(i)||/||b|| 1.294915694729e+12
192 KSP preconditioned resid norm 2.756924664914e+00 true resid norm 2.167287516776e+14 ||r(i)||/||b|| 1.234219354251e+12
193 KSP preconditioned resid norm 2.712358962874e+00 true resid norm 2.044300382742e+14 ||r(i)||/||b|| 1.164181069080e+12
194 KSP preconditioned resid norm 2.634447232139e+00 true resid norm 1.861632836558e+14 ||r(i)||/||b|| 1.060156190448e+12
195 KSP preconditioned resid norm 2.554392023213e+00 true resid norm 1.695435830513e+14 ||r(i)||/||b|| 9.655108977062e+11
196 KSP preconditioned resid norm 2.423823542815e+00 true resid norm 1.453834575243e+14 ||r(i)||/||b|| 8.279246554760e+11
197 KSP preconditioned resid norm 2.301144929663e+00 true resid norm 1.254689122515e+14 ||r(i)||/||b|| 7.145159959579e+11
198 KSP preconditioned resid norm 2.127715354714e+00 true resid norm 1.001847989583e+14 ||r(i)||/||b|| 5.705289073048e+11
199 KSP preconditioned resid norm 2.009569485205e+00 true resid norm 8.505075510895e+13 ||r(i)||/||b|| 4.843440809612e+11
200 KSP preconditioned resid norm 1.877928091553e+00 true resid norm 7.067964496444e+13 ||r(i)||/||b|| 4.025039829348e+11
201 KSP preconditioned resid norm 1.742305371081e+00 true resid norm 5.778287055757e+13 ||r(i)||/||b|| 3.290598807695e+11
202 KSP preconditioned resid norm 1.552693582614e+00 true resid norm 4.295277252766e+13 ||r(i)||/||b|| 2.446059545033e+11
203 KSP preconditioned resid norm 1.405130331702e+00 true resid norm 3.341307717615e+13 ||r(i)||/||b|| 1.902796293370e+11
204 KSP preconditioned resid norm 1.207664877323e+00 true resid norm 2.297639302916e+13 ||r(i)||/||b|| 1.308451635880e+11
205 KSP preconditioned resid norm 1.053317227920e+00 true resid norm 1.622811796763e+13 ||r(i)||/||b|| 9.241532156523e+10
206 KSP preconditioned resid norm 8.888835507681e-01 true resid norm 1.056119524799e+13 ||r(i)||/||b|| 6.014352723484e+10
207 KSP preconditioned resid norm 7.509558493683e-01 true resid norm 6.806255926965e+12 ||r(i)||/||b|| 3.876002943783e+10
208 KSP preconditioned resid norm 6.625677908839e-01 true resid norm 4.847382413253e+12 ||r(i)||/||b|| 2.760470470847e+10
209 KSP preconditioned resid norm 5.988641055583e-01 true resid norm 3.692435549416e+12 ||r(i)||/||b|| 2.102755349320e+10
210 KSP preconditioned resid norm 5.549019306325e-01 true resid norm 3.047376495457e+12 ||r(i)||/||b|| 1.735409363673e+10
211 KSP preconditioned resid norm 5.455063254750e-01 true resid norm 3.068911085737e+12 ||r(i)||/||b|| 1.747672807218e+10
212 KSP preconditioned resid norm 5.256033146031e-01 true resid norm 3.117144722512e+12 ||r(i)||/||b|| 1.775140730866e+10
213 KSP preconditioned resid norm 5.188506724092e-01 true resid norm 3.135841625267e+12 ||r(i)||/||b|| 1.785788177993e+10
214 KSP preconditioned resid norm 5.064655860471e-01 true resid norm 3.177028328480e+12 ||r(i)||/||b|| 1.809243038435e+10
215 KSP preconditioned resid norm 4.942961912759e-01 true resid norm 3.227522585036e+12 ||r(i)||/||b|| 1.837998331971e+10
216 KSP preconditioned resid norm 4.748831446073e-01 true resid norm 3.319195929690e+12 ||r(i)||/||b|| 1.890204149318e+10
217 KSP preconditioned resid norm 4.515763724804e-01 true resid norm 3.439076398391e+12 ||r(i)||/||b|| 1.958473261525e+10
218 KSP preconditioned resid norm 4.373634838140e-01 true resid norm 3.528928306121e+12 ||r(i)||/||b|| 2.009641813311e+10
219 KSP preconditioned resid norm 4.009407832729e-01 true resid norm 3.761064602302e+12 ||r(i)||/||b|| 2.141837983572e+10
220 KSP preconditioned resid norm 3.740641614089e-01 true resid norm 3.936779171425e+12 ||r(i)||/||b|| 2.241903305019e+10
221 KSP preconditioned resid norm 3.517620162262e-01 true resid norm 4.097027645608e+12 ||r(i)||/||b|| 2.333161048533e+10
222 KSP preconditioned resid norm 3.229170697014e-01 true resid norm 4.313508830841e+12 ||r(i)||/||b|| 2.456442000680e+10
223 KSP preconditioned resid norm 2.994650896035e-01 true resid norm 4.502568103319e+12 ||r(i)||/||b|| 2.564106817363e+10
224 KSP preconditioned resid norm 2.773846058866e-01 true resid norm 4.710002710255e+12 ||r(i)||/||b|| 2.682235955578e+10
225 KSP preconditioned resid norm 2.501291804794e-01 true resid norm 4.969587839164e+12 ||r(i)||/||b|| 2.830063591596e+10
226 KSP preconditioned resid norm 2.304709168381e-01 true resid norm 5.185431689243e+12 ||r(i)||/||b|| 2.952981596337e+10
227 KSP preconditioned resid norm 2.100777905891e-01 true resid norm 5.450168294124e+12 ||r(i)||/||b|| 3.103742876968e+10
228 KSP preconditioned resid norm 1.902333022433e-01 true resid norm 5.764157521664e+12 ||r(i)||/||b|| 3.282552369782e+10
229 KSP preconditioned resid norm 1.696089510455e-01 true resid norm 6.179725445319e+12 ||r(i)||/||b|| 3.519208545029e+10
230 KSP preconditioned resid norm 1.581595533833e-01 true resid norm 6.486809236950e+12 ||r(i)||/||b|| 3.694085554227e+10
231 KSP preconditioned resid norm 1.508618057318e-01 true resid norm 6.725919801784e+12 ||r(i)||/||b|| 3.830253406734e+10
232 KSP preconditioned resid norm 1.423058782095e-01 true resid norm 7.059550776628e+12 ||r(i)||/||b|| 4.020248413462e+10
233 KSP preconditioned resid norm 1.315340285601e-01 true resid norm 7.576206972091e+12 ||r(i)||/||b|| 4.314471986014e+10
234 KSP preconditioned resid norm 1.237625992155e-01 true resid norm 8.062848589615e+12 ||r(i)||/||b|| 4.591602961154e+10
235 KSP preconditioned resid norm 1.167453140487e-01 true resid norm 8.626376541470e+12 ||r(i)||/||b|| 4.912518898453e+10
236 KSP preconditioned resid norm 1.118147547475e-01 true resid norm 9.193248341542e+12 ||r(i)||/||b|| 5.235339078799e+10
237 KSP preconditioned resid norm 1.093201267993e-01 true resid norm 9.616798592267e+12 ||r(i)||/||b|| 5.476541001892e+10
238 KSP preconditioned resid norm 1.083570175037e-01 true resid norm 9.931570457757e+12 ||r(i)||/||b|| 5.655796188643e+10
239 KSP preconditioned resid norm 1.082842938953e-01 true resid norm 1.003792921081e+13 ||r(i)||/||b|| 5.716365001267e+10
240 KSP preconditioned resid norm 1.078314038130e-01 true resid norm 9.849622726373e+12 ||r(i)||/||b|| 5.609128879701e+10
241 KSP preconditioned resid norm 1.077292227053e-01 true resid norm 9.870891496109e+12 ||r(i)||/||b|| 5.621240944688e+10
242 KSP preconditioned resid norm 1.075373948198e-01 true resid norm 9.914536589062e+12 ||r(i)||/||b|| 5.646095800365e+10
243 KSP preconditioned resid norm 1.074752343474e-01 true resid norm 9.944186669264e+12 ||r(i)||/||b|| 5.662980824875e+10
244 KSP preconditioned resid norm 1.074752341803e-01 true resid norm 9.944105914292e+12 ||r(i)||/||b|| 5.662934836815e+10
245 KSP preconditioned resid norm 1.074690883936e-01 true resid norm 9.924451054116e+12 ||r(i)||/||b|| 5.651741855428e+10
246 KSP preconditioned resid norm 1.073632398521e-01 true resid norm 9.811970192974e+12 ||r(i)||/||b|| 5.587686645988e+10
247 KSP preconditioned resid norm 1.072870146724e-01 true resid norm 9.751045941996e+12 ||r(i)||/||b|| 5.552991715520e+10
248 KSP preconditioned resid norm 1.067010059812e-01 true resid norm 9.474108561586e+12 ||r(i)||/||b|| 5.395282379692e+10
249 KSP preconditioned resid norm 1.055793587548e-01 true resid norm 9.010840121971e+12 ||r(i)||/||b|| 5.131461880584e+10
250 KSP preconditioned resid norm 1.047262495938e-01 true resid norm 8.722549728520e+12 ||r(i)||/||b|| 4.967287270391e+10
251 KSP preconditioned resid norm 1.029466337065e-01 true resid norm 8.186021332697e+12 ||r(i)||/||b|| 4.661746946320e+10
252 KSP preconditioned resid norm 1.004304639605e-01 true resid norm 7.495626733836e+12 ||r(i)||/||b|| 4.268583432301e+10
253 KSP preconditioned resid norm 9.795540774722e-02 true resid norm 6.920483082409e+12 ||r(i)||/||b|| 3.941052626826e+10
254 KSP preconditioned resid norm 9.555358753906e-02 true resid norm 6.409473015860e+12 ||r(i)||/||b|| 3.650044392123e+10
255 KSP preconditioned resid norm 9.212232118013e-02 true resid norm 5.754525445666e+12 ||r(i)||/||b|| 3.277067128656e+10
256 KSP preconditioned resid norm 8.758697690555e-02 true resid norm 4.960269002563e+12 ||r(i)||/||b|| 2.824756732952e+10
257 KSP preconditioned resid norm 8.129606182327e-02 true resid norm 4.021896818151e+12 ||r(i)||/||b|| 2.290375806320e+10
258 KSP preconditioned resid norm 7.531177419725e-02 true resid norm 3.231819088719e+12 ||r(i)||/||b|| 1.840445089939e+10
259 KSP preconditioned resid norm 7.147261886049e-02 true resid norm 2.780265759853e+12 ||r(i)||/||b|| 1.583296071339e+10
260 KSP preconditioned resid norm 6.629473978281e-02 true resid norm 2.268089888953e+12 ||r(i)||/||b|| 1.291623938430e+10
261 KSP preconditioned resid norm 6.208831287626e-02 true resid norm 1.905334327725e+12 ||r(i)||/||b|| 1.085043163584e+10
262 KSP preconditioned resid norm 5.627448281017e-02 true resid norm 1.483140584032e+12 ||r(i)||/||b|| 8.446137393950e+09
263 KSP preconditioned resid norm 4.989994432902e-02 true resid norm 1.096525121051e+12 ||r(i)||/||b|| 6.244453107160e+09
264 KSP preconditioned resid norm 4.497048672753e-02 true resid norm 8.467526728715e+11 ||r(i)||/||b|| 4.822057659783e+09
265 KSP preconditioned resid norm 3.747768169066e-02 true resid norm 5.348775283731e+11 ||r(i)||/||b|| 3.046001938193e+09
266 KSP preconditioned resid norm 3.272507533606e-02 true resid norm 3.761068823582e+11 ||r(i)||/||b|| 2.141840387491e+09
267 KSP preconditioned resid norm 2.833978687737e-02 true resid norm 2.584061281222e+11 ||r(i)||/||b|| 1.471562227517e+09
268 KSP preconditioned resid norm 2.427354160616e-02 true resid norm 1.695159727985e+11 ||r(i)||/||b|| 9.653536637995e+08
269 KSP preconditioned resid norm 2.223267559447e-02 true resid norm 1.346280124342e+11 ||r(i)||/||b|| 7.666749209995e+08
270 KSP preconditioned resid norm 2.051292146863e-02 true resid norm 1.110744657811e+11 ||r(i)||/||b|| 6.325430030351e+08
271 KSP preconditioned resid norm 2.015852923458e-02 true resid norm 1.116963414691e+11 ||r(i)||/||b|| 6.360844390658e+08
272 KSP preconditioned resid norm 1.928845567194e-02 true resid norm 1.133138718159e+11 ||r(i)||/||b|| 6.452958946045e+08
273 KSP preconditioned resid norm 1.901639059774e-02 true resid norm 1.139195524924e+11 ||r(i)||/||b|| 6.487451038471e+08
274 KSP preconditioned resid norm 1.867369139213e-02 true resid norm 1.148722990709e+11 ||r(i)||/||b|| 6.541707719128e+08
275 KSP preconditioned resid norm 1.820272269340e-02 true resid norm 1.164715379777e+11 ||r(i)||/||b|| 6.632780619961e+08
276 KSP preconditioned resid norm 1.749869482046e-02 true resid norm 1.194547750039e+11 ||r(i)||/||b|| 6.802668964145e+08
277 KSP preconditioned resid norm 1.687368081121e-02 true resid norm 1.224936003716e+11 ||r(i)||/||b|| 6.975722933867e+08
278 KSP preconditioned resid norm 1.639575587312e-02 true resid norm 1.252604862821e+11 ||r(i)||/||b|| 7.133290589993e+08
279 KSP preconditioned resid norm 1.496512503896e-02 true resid norm 1.337998974184e+11 ||r(i)||/||b|| 7.619590004206e+08
280 KSP preconditioned resid norm 1.418029241433e-02 true resid norm 1.386300487682e+11 ||r(i)||/||b|| 7.894655782685e+08
281 KSP preconditioned resid norm 1.331740374144e-02 true resid norm 1.444340453171e+11 ||r(i)||/||b|| 8.225179758723e+08
282 KSP preconditioned resid norm 1.222653019600e-02 true resid norm 1.521351041691e+11 ||r(i)||/||b|| 8.663736978743e+08
283 KSP preconditioned resid norm 1.131424510923e-02 true resid norm 1.595359004069e+11 ||r(i)||/||b|| 9.085194947878e+08
284 KSP preconditioned resid norm 1.037042507094e-02 true resid norm 1.680264407990e+11 ||r(i)||/||b|| 9.568711287950e+08
285 KSP preconditioned resid norm 9.688182631346e-03 true resid norm 1.743226122924e+11 ||r(i)||/||b|| 9.927263471485e+08
286 KSP preconditioned resid norm 9.052103564731e-03 true resid norm 1.811486798072e+11 ||r(i)||/||b|| 1.031599198927e+09
287 KSP preconditioned resid norm 8.110652749532e-03 true resid norm 1.926557561601e+11 ||r(i)||/||b|| 1.097129296968e+09
288 KSP preconditioned resid norm 7.508113362095e-03 true resid norm 2.010760992496e+11 ||r(i)||/||b|| 1.145081173818e+09
289 KSP preconditioned resid norm 6.917230248586e-03 true resid norm 2.109956414742e+11 ||r(i)||/||b|| 1.201570637741e+09
290 KSP preconditioned resid norm 6.270792144381e-03 true resid norm 2.235532529095e+11 ||r(i)||/||b|| 1.273083286417e+09
291 KSP preconditioned resid norm 5.723434422299e-03 true resid norm 2.366863243489e+11 ||r(i)||/||b|| 1.347873044702e+09
292 KSP preconditioned resid norm 5.384057623534e-03 true resid norm 2.476827107725e+11 ||r(i)||/||b|| 1.410494883501e+09
293 KSP preconditioned resid norm 5.001836333933e-03 true resid norm 2.646740880851e+11 ||r(i)||/||b|| 1.507256787828e+09
294 KSP preconditioned resid norm 4.705081336940e-03 true resid norm 2.809588439377e+11 ||r(i)||/||b|| 1.599994648849e+09
295 KSP preconditioned resid norm 4.424589025495e-03 true resid norm 3.023646578354e+11 ||r(i)||/||b|| 1.721895733046e+09
296 KSP preconditioned resid norm 4.204809966953e-03 true resid norm 3.239120797248e+11 ||r(i)||/||b|| 1.844603241506e+09
297 KSP preconditioned resid norm 4.020983239105e-03 true resid norm 3.491084935202e+11 ||r(i)||/||b|| 1.988090902111e+09
298 KSP preconditioned resid norm 3.962918432092e-03 true resid norm 3.653308451948e+11 ||r(i)||/||b|| 2.080473385992e+09
299 KSP preconditioned resid norm 3.957383148812e-03 true resid norm 3.700697096946e+11 ||r(i)||/||b|| 2.107460106662e+09
300 KSP preconditioned resid norm 3.939097482021e-03 true resid norm 3.622457576700e+11 ||r(i)||/||b|| 2.062904536897e+09
301 KSP preconditioned resid norm 3.936139794992e-03 true resid norm 3.628684964620e+11 ||r(i)||/||b|| 2.066450888102e+09
302 KSP preconditioned resid norm 3.931455705928e-03 true resid norm 3.639581395965e+11 ||r(i)||/||b|| 2.072656149912e+09
303 KSP preconditioned resid norm 3.926424749526e-03 true resid norm 3.655814698003e+11 ||r(i)||/||b|| 2.081900634275e+09
304 KSP preconditioned resid norm 3.926009575006e-03 true resid norm 3.661253634305e+11 ||r(i)||/||b|| 2.084997980796e+09
305 KSP preconditioned resid norm 3.925339982792e-03 true resid norm 3.649034775025e+11 ||r(i)||/||b|| 2.078039627327e+09
306 KSP preconditioned resid norm 3.923356505149e-03 true resid norm 3.626519814581e+11 ||r(i)||/||b|| 2.065217885991e+09
307 KSP preconditioned resid norm 3.917014224631e-03 true resid norm 3.580396053341e+11 ||r(i)||/||b|| 2.038951486923e+09
308 KSP preconditioned resid norm 3.903303707944e-03 true resid norm 3.514819854106e+11 ||r(i)||/||b|| 2.001607381146e+09
309 KSP preconditioned resid norm 3.868422381327e-03 true resid norm 3.361181915994e+11 ||r(i)||/||b|| 1.914114182714e+09
310 KSP preconditioned resid norm 3.830516524577e-03 true resid norm 3.221653520608e+11 ||r(i)||/||b|| 1.834656037581e+09
311 KSP preconditioned resid norm 3.796415788584e-03 true resid norm 3.115912412675e+11 ||r(i)||/||b|| 1.774438959348e+09
312 KSP preconditioned resid norm 3.690310261976e-03 true resid norm 2.821305118145e+11 ||r(i)||/||b|| 1.606667022308e+09
313 KSP preconditioned resid norm 3.593395172499e-03 true resid norm 2.581263239915e+11 ||r(i)||/||b|| 1.469968808688e+09
314 KSP preconditioned resid norm 3.480205774823e-03 true resid norm 2.333504532643e+11 ||r(i)||/||b|| 1.328876042116e+09
315 KSP preconditioned resid norm 3.353602368302e-03 true resid norm 2.081533340394e+11 ||r(i)||/||b|| 1.185384364256e+09
316 KSP preconditioned resid norm 3.194409146950e-03 true resid norm 1.800154987321e+11 ||r(i)||/||b|| 1.025145998768e+09
317 KSP preconditioned resid norm 3.071672481669e-03 true resid norm 1.608433892393e+11 ||r(i)||/||b|| 9.159653366981e+08
318 KSP preconditioned resid norm 2.839844293931e-03 true resid norm 1.286439364064e+11 ||r(i)||/||b|| 7.325970130446e+08
319 KSP preconditioned resid norm 2.631447463216e-03 true resid norm 1.035755023455e+11 ||r(i)||/||b|| 5.898381669786e+08
320 KSP preconditioned resid norm 2.421709878855e-03 true resid norm 8.242231041182e+10 ||r(i)||/||b|| 4.693757055534e+08
321 KSP preconditioned resid norm 2.249677749119e-03 true resid norm 6.796005580822e+10 ||r(i)||/||b|| 3.870165612326e+08
322 KSP preconditioned resid norm 2.051679744073e-03 true resid norm 5.362063892849e+10 ||r(i)||/||b|| 3.053569489078e+08
323 KSP preconditioned resid norm 1.852618445705e-03 true resid norm 4.140290200052e+10 ||r(i)||/||b|| 2.357798057510e+08
324 KSP preconditioned resid norm 1.586867536205e-03 true resid norm 2.821223751515e+10 ||r(i)||/||b|| 1.606620685922e+08
325 KSP preconditioned resid norm 1.397039923183e-03 true resid norm 2.041649071932e+10 ||r(i)||/||b|| 1.162671209825e+08
326 KSP preconditioned resid norm 1.221057962383e-03 true resid norm 1.446679933564e+10 ||r(i)||/||b|| 8.238502550265e+07
327 KSP preconditioned resid norm 1.033496887800e-03 true resid norm 9.360606892139e+09 ||r(i)||/||b|| 5.330645843890e+07
328 KSP preconditioned resid norm 9.116911789257e-04 true resid norm 6.661951473115e+09 ||r(i)||/||b|| 3.793824945494e+07
329 KSP preconditioned resid norm 8.282610842511e-04 true resid norm 5.144839330536e+09 ||r(i)||/||b|| 2.929865201137e+07
330 KSP preconditioned resid norm 7.530731044695e-04 true resid norm 4.075530635865e+09 ||r(i)||/||b|| 2.320919006220e+07
331 KSP preconditioned resid norm 7.379285385788e-04 true resid norm 4.106395550093e+09 ||r(i)||/||b|| 2.338495850184e+07
332 KSP preconditioned resid norm 7.208624023923e-04 true resid norm 4.143156222668e+09 ||r(i)||/||b|| 2.359430190098e+07
333 KSP preconditioned resid norm 7.065775764655e-04 true resid norm 4.177514984231e+09 ||r(i)||/||b|| 2.378996698086e+07
334 KSP preconditioned resid norm 6.922452814868e-04 true resid norm 4.218601053057e+09 ||r(i)||/||b|| 2.402394249607e+07
335 KSP preconditioned resid norm 6.744135487420e-04 true resid norm 4.282753533716e+09 ||r(i)||/||b|| 2.438927581083e+07
336 KSP preconditioned resid norm 6.469130429753e-04 true resid norm 4.403672957096e+09 ||r(i)||/||b|| 2.507788353585e+07
337 KSP preconditioned resid norm 6.112369725507e-04 true resid norm 4.573122433971e+09 ||r(i)||/||b|| 2.604285852098e+07
338 KSP preconditioned resid norm 5.788565956748e-04 true resid norm 4.751840775440e+09 ||r(i)||/||b|| 2.706061751370e+07
339 KSP preconditioned resid norm 5.410631845567e-04 true resid norm 4.973946362559e+09 ||r(i)||/||b|| 2.832545668334e+07
340 KSP preconditioned resid norm 5.009844336682e-04 true resid norm 5.235146248376e+09 ||r(i)||/||b|| 2.981292870497e+07
341 KSP preconditioned resid norm 4.659768311587e-04 true resid norm 5.473524476167e+09 ||r(i)||/||b|| 3.117043674253e+07
342 KSP preconditioned resid norm 4.188668678783e-04 true resid norm 5.805084466010e+09 ||r(i)||/||b|| 3.305859303647e+07
343 KSP preconditioned resid norm 3.879486344930e-04 true resid norm 6.063617165803e+09 ||r(i)||/||b|| 3.453087605994e+07
344 KSP preconditioned resid norm 3.714346767924e-04 true resid norm 6.230541109727e+09 ||r(i)||/||b|| 3.548146872789e+07
345 KSP preconditioned resid norm 3.513107945823e-04 true resid norm 6.444239345039e+09 ||r(i)||/||b|| 3.669843000299e+07
346 KSP preconditioned resid norm 3.188508242839e-04 true resid norm 6.808853174711e+09 ||r(i)||/||b|| 3.877482015393e+07
347 KSP preconditioned resid norm 2.976989911735e-04 true resid norm 7.088594142323e+09 ||r(i)||/||b|| 4.036787928305e+07
348 KSP preconditioned resid norm 2.760902621649e-04 true resid norm 7.404251197607e+09 ||r(i)||/||b|| 4.216547209859e+07
349 KSP preconditioned resid norm 2.491238069081e-04 true resid norm 7.856460563260e+09 ||r(i)||/||b|| 4.474069826006e+07
350 KSP preconditioned resid norm 2.261931392463e-04 true resid norm 8.300304774785e+09 ||r(i)||/||b|| 4.726828683286e+07
351 KSP preconditioned resid norm 2.108262997601e-04 true resid norm 8.660994767295e+09 ||r(i)||/||b|| 4.932233165246e+07
352 KSP preconditioned resid norm 1.952804314878e-04 true resid norm 9.171687969602e+09 ||r(i)||/||b|| 5.223060953203e+07
353 KSP preconditioned resid norm 1.857561963852e-04 true resid norm 9.612624044306e+09 ||r(i)||/||b|| 5.474163694844e+07
354 KSP preconditioned resid norm 1.755682138867e-04 true resid norm 1.022824097483e+10 ||r(i)||/||b|| 5.824743082478e+07
355 KSP preconditioned resid norm 1.629922761938e-04 true resid norm 1.128526888154e+10 ||r(i)||/||b|| 6.426695656995e+07
356 KSP preconditioned resid norm 1.523745456214e-04 true resid norm 1.237228264869e+10 ||r(i)||/||b|| 7.045724474987e+07
357 KSP preconditioned resid norm 1.479156994814e-04 true resid norm 1.303748936001e+10 ||r(i)||/||b|| 7.424544078444e+07
358 KSP preconditioned resid norm 1.460139994385e-04 true resid norm 1.352484512512e+10 ||r(i)||/||b|| 7.702081743870e+07
359 KSP preconditioned resid norm 1.458891367415e-04 true resid norm 1.366183537709e+10 ||r(i)||/||b|| 7.780094475922e+07
360 KSP preconditioned resid norm 1.455077198104e-04 true resid norm 1.345578836961e+10 ||r(i)||/||b|| 7.662755543012e+07
361 KSP preconditioned resid norm 1.453704662869e-04 true resid norm 1.350023349252e+10 ||r(i)||/||b|| 7.688065996966e+07
362 KSP preconditioned resid norm 1.452981914522e-04 true resid norm 1.353022367766e+10 ||r(i)||/||b|| 7.705144703254e+07
363 KSP preconditioned resid norm 1.452129596669e-04 true resid norm 1.358122759269e+10 ||r(i)||/||b|| 7.734190235326e+07
364 KSP preconditioned resid norm 1.452121241088e-04 true resid norm 1.358710317940e+10 ||r(i)||/||b|| 7.737536244000e+07
365 KSP preconditioned resid norm 1.452046967752e-04 true resid norm 1.356520699403e+10 ||r(i)||/||b|| 7.725066880542e+07
366 KSP preconditioned resid norm 1.450384393428e-04 true resid norm 1.342016311361e+10 ||r(i)||/||b|| 7.642467796180e+07
367 KSP preconditioned resid norm 1.448709937803e-04 true resid norm 1.330896122194e+10 ||r(i)||/||b|| 7.579140929825e+07
368 KSP preconditioned resid norm 1.440619441880e-04 true resid norm 1.296138556368e+10 ||r(i)||/||b|| 7.381204753305e+07
369 KSP preconditioned resid norm 1.424189519513e-04 true resid norm 1.234231668356e+10 ||r(i)||/||b|| 7.028659561425e+07
370 KSP preconditioned resid norm 1.405178593540e-04 true resid norm 1.171445058083e+10 ||r(i)||/||b|| 6.671104557822e+07
371 KSP preconditioned resid norm 1.384125397526e-04 true resid norm 1.109105276720e+10 ||r(i)||/||b|| 6.316094139954e+07
372 KSP preconditioned resid norm 1.349489907590e-04 true resid norm 1.017955796073e+10 ||r(i)||/||b|| 5.797019248993e+07
373 KSP preconditioned resid norm 1.299224645965e-04 true resid norm 9.012691544632e+09 ||r(i)||/||b|| 5.132516222319e+07
374 KSP preconditioned resid norm 1.244108688000e-04 true resid norm 7.863710058584e+09 ||r(i)||/||b|| 4.478198243380e+07
375 KSP preconditioned resid norm 1.196624061573e-04 true resid norm 6.954908378743e+09 ||r(i)||/||b|| 3.960657024804e+07
376 KSP preconditioned resid norm 1.150069930686e-04 true resid norm 6.146973261882e+09 ||r(i)||/||b|| 3.500556945562e+07
377 KSP preconditioned resid norm 1.078887656664e-04 true resid norm 5.102072199579e+09 ||r(i)||/||b|| 2.905510324203e+07
378 KSP preconditioned resid norm 1.020396071938e-04 true resid norm 4.354821250020e+09 ||r(i)||/||b|| 2.479968453413e+07
379 KSP preconditioned resid norm 9.618818827534e-05 true resid norm 3.697423694068e+09 ||r(i)||/||b|| 2.105595980580e+07
380 KSP preconditioned resid norm 9.035897971170e-05 true resid norm 3.119110900940e+09 ||r(i)||/||b|| 1.776260423316e+07
KSP Object: 4 MPI processes
type: gmres
GMRES: restart=30, using Classical (unmodified) Gram-Schmidt Orthogonalization with no iterative refinement
GMRES: happy breakdown tolerance 1e-30
maximum iterations=10000, initial guess is zero
tolerances: relative=1e-10, absolute=1e-50, divergence=10000
left preconditioning
using PRECONDITIONED norm type for convergence test
PC Object: 4 MPI processes
type: fieldsplit
FieldSplit with Schur preconditioner, blocksize = 1, factorization LOWER
Preconditioner for the Schur complement formed from S itself
Split info:
Split number 0 Defined by IS
Split number 1 Defined by IS
KSP solver for A00 block
KSP Object: (fieldsplit_0_) 4 MPI processes
type: preonly
maximum iterations=10000, initial guess is zero
tolerances: relative=1e-05, absolute=1e-50, divergence=10000
left preconditioning
using NONE norm type for convergence test
PC Object: (fieldsplit_0_) 4 MPI processes
type: lu
LU: out-of-place factorization
tolerance for zero pivot 2.22045e-14
matrix ordering: natural
factor fill ratio given 0, needed 0
Factored matrix follows:
Mat Object: 4 MPI processes
type: mpiaij
rows=71946, cols=71946
package used to perform factorization: mumps
total: nonzeros=6.47715e+07, allocated nonzeros=6.47715e+07
total number of mallocs used during MatSetValues calls =0
MUMPS run parameters:
SYM (matrix type): 0
PAR (host participation): 1
ICNTL(1) (output for error): 6
ICNTL(2) (output of diagnostic msg): 0
ICNTL(3) (output for global info): 0
ICNTL(4) (level of printing): 0
ICNTL(5) (input mat struct): 0
ICNTL(6) (matrix prescaling): 7
ICNTL(7) (sequentia matrix ordering):7
ICNTL(8) (scalling strategy): 77
ICNTL(10) (max num of refinements): 0
ICNTL(11) (error analysis): 0
ICNTL(12) (efficiency control): 1
ICNTL(13) (efficiency control): 0
ICNTL(14) (percentage of estimated workspace increase): 20
ICNTL(18) (input mat struct): 3
ICNTL(19) (Shur complement info): 0
ICNTL(20) (rhs sparse pattern): 0
ICNTL(21) (somumpstion struct): 1
ICNTL(22) (in-core/out-of-core facility): 0
ICNTL(23) (max size of memory can be allocated locally):0
ICNTL(24) (detection of null pivot rows): 0
ICNTL(25) (computation of a null space basis): 0
ICNTL(26) (Schur options for rhs or solution): 0
ICNTL(27) (experimental parameter): -8
ICNTL(28) (use parallel or sequential ordering): 1
ICNTL(29) (parallel ordering): 0
ICNTL(30) (user-specified set of entries in inv(A)): 0
ICNTL(31) (factors is discarded in the solve phase): 0
ICNTL(33) (compute determinant): 0
CNTL(1) (relative pivoting threshold): 0.01
CNTL(2) (stopping criterion of refinement): 1.49012e-08
CNTL(3) (absomumpste pivoting threshold): 0
CNTL(4) (vamumpse of static pivoting): -1
CNTL(5) (fixation for null pivots): 0
RINFO(1) (local estimated flops for the elimination after analysis):
[0] 1.6957e+10
[1] 1.67964e+10
[2] 1.70959e+10
[3] 1.7073e+10
RINFO(2) (local estimated flops for the assembly after factorization):
[0] 3.99956e+07
[1] 3.71513e+07
[2] 4.43658e+07
[3] 2.95262e+07
RINFO(3) (local estimated flops for the elimination after factorization):
[0] 1.732e+10
[1] 1.7234e+10
[2] 1.70959e+10
[3] 1.62723e+10
INFO(15) (estimated size of (in MB) MUMPS internal data for running numerical factorization):
[0] 243
[1] 237
[2] 256
[3] 212
INFO(16) (size of (in MB) MUMPS internal data used during numerical factorization):
[0] 243
[1] 237
[2] 256
[3] 212
INFO(23) (num of pivots eliminated on this processor after factorization):
[0] 23748
[1] 15630
[2] 20355
[3] 12213
RINFOG(1) (global estimated flops for the elimination after analysis): 6.79222e+10
RINFOG(2) (global estimated flops for the assembly after factorization): 1.51039e+08
RINFOG(3) (global estimated flops for the elimination after factorization): 6.79222e+10
(RINFOG(12) RINFOG(13))*2^INFOG(34) (determinant): (0,0)*(2^0)
INFOG(3) (estimated real workspace for factors on all processors after analysis): 64771542
INFOG(4) (estimated integer workspace for factors on all processors after analysis): 1209207
INFOG(5) (estimated maximum front size in the complete tree): 2562
INFOG(6) (number of nodes in the complete tree): 3832
INFOG(7) (ordering option effectively use after analysis): 5
INFOG(8) (structural symmetry in percent of the permuted matrix after analysis): 100
INFOG(9) (total real/complex workspace to store the matrix factors after factorization): 64771542
INFOG(10) (total integer space store the matrix factors after factorization): 1208659
INFOG(11) (order of largest frontal matrix after factorization): 2562
INFOG(12) (number of off-diagonal pivots): 0
INFOG(13) (number of delayed pivots after factorization): 0
INFOG(14) (number of memory compress after factorization): 0
INFOG(15) (number of steps of iterative refinement after solution): 0
INFOG(16) (estimated size (in MB) of all MUMPS internal data for factorization after analysis: value on the most memory consuming processor): 256
INFOG(17) (estimated size of all MUMPS internal data for factorization after analysis: sum over all processors): 948
INFOG(18) (size of all MUMPS internal data allocated during factorization: value on the most memory consuming processor): 256
INFOG(19) (size of all MUMPS internal data allocated during factorization: sum over all processors): 948
INFOG(20) (estimated number of entries in the factors): 64771542
INFOG(21) (size in MB of memory effectively used during factorization - value on the most memory consuming processor): 212
INFOG(22) (size in MB of memory effectively used during factorization - sum over all processors): 789
INFOG(23) (after analysis: value of ICNTL(6) effectively used): 0
INFOG(24) (after analysis: value of ICNTL(12) effectively used): 1
INFOG(25) (after factorization: number of pivots modified by static pivoting): 0
INFOG(28) (after factorization: number of null pivots encountered): 0
INFOG(29) (after factorization: effective number of entries in the factors (sum over all processors)): 58207698
INFOG(30, 31) (after solution: size in Mbytes of memory used during solution phase): 229, 844
INFOG(32) (after analysis: type of analysis done): 1
INFOG(33) (value used for ICNTL(8)): 7
INFOG(34) (exponent of the determinant if determinant is requested): 0
linear system matrix = precond matrix:
Mat Object: (fieldsplit_0_) 4 MPI processes
type: mpiaij
rows=71946, cols=71946
total: nonzeros=2.96023e+06, allocated nonzeros=2.96023e+06
total number of mallocs used during MatSetValues calls =0
using I-node (on process 0) routines: found 6259 nodes, limit used is 5
KSP solver for S = A11 - A10 inv(A00) A01
KSP Object: (fieldsplit_1_) 4 MPI processes
type: preonly
maximum iterations=10000, initial guess is zero
tolerances: relative=1e-05, absolute=1e-50, divergence=10000
left preconditioning
using NONE norm type for convergence test
PC Object: (fieldsplit_1_) 4 MPI processes
type: lsc
KSP Object: (fieldsplit_1_lsc_) 4 MPI processes
type: preonly
maximum iterations=10000, initial guess is zero
tolerances: relative=1e-05, absolute=1e-50, divergence=10000
left preconditioning
using NONE norm type for convergence test
PC Object: (fieldsplit_1_lsc_) 4 MPI processes
type: bjacobi
block Jacobi: number of blocks = 4
Local solve is same for all blocks, in the following KSP and PC objects:
KSP Object: (fieldsplit_1_lsc_sub_) 1 MPI processes
type: preonly
maximum iterations=10000, initial guess is zero
tolerances: relative=1e-05, absolute=1e-50, divergence=10000
left preconditioning
using NONE norm type for convergence test
PC Object: (fieldsplit_1_lsc_sub_) 1 MPI processes
type: ilu
ILU: out-of-place factorization
0 levels of fill
tolerance for zero pivot 2.22045e-14
using diagonal shift on blocks to prevent zero pivot [INBLOCKS]
matrix ordering: natural
factor fill ratio given 1, needed 1
Factored matrix follows:
Mat Object: 1 MPI processes
type: seqaij
rows=0, cols=0
package used to perform factorization: petsc
total: nonzeros=0, allocated nonzeros=0
total number of mallocs used during MatSetValues calls =0
not using I-node routines
linear system matrix = precond matrix:
Mat Object: 1 MPI processes
type: seqaij
rows=0, cols=0
total: nonzeros=0, allocated nonzeros=0
total number of mallocs used during MatSetValues calls =0
not using I-node routines
linear system matrix = precond matrix:
Mat Object: 4 MPI processes
type: mpiaij
rows=3183, cols=3183
total: nonzeros=9549, allocated nonzeros=9549
total number of mallocs used during MatSetValues calls =0
not using I-node (on process 0) routines
linear system matrix = precond matrix:
Mat Object: (fieldsplit_1_) 4 MPI processes
type: schurcomplement
rows=3183, cols=3183
Schur complement A11 - A10 inv(A00) A01
A11
Mat Object: (fieldsplit_1_) 4 MPI processes
type: mpiaij
rows=3183, cols=3183
total: nonzeros=3183, allocated nonzeros=3183
total number of mallocs used during MatSetValues calls =0
not using I-node (on process 0) routines
A10
Mat Object: 4 MPI processes
type: mpiaij
rows=3183, cols=71946
total: nonzeros=19098, allocated nonzeros=19098
total number of mallocs used during MatSetValues calls =0
not using I-node (on process 0) routines
KSP of A00
KSP Object: (fieldsplit_0_) 4 MPI processes
type: preonly
maximum iterations=10000, initial guess is zero
tolerances: relative=1e-05, absolute=1e-50, divergence=10000
left preconditioning
using NONE norm type for convergence test
PC Object: (fieldsplit_0_) 4 MPI processes
type: lu
LU: out-of-place factorization
tolerance for zero pivot 2.22045e-14
matrix ordering: natural
factor fill ratio given 0, needed 0
Factored matrix follows:
Mat Object: 4 MPI processes
type: mpiaij
rows=71946, cols=71946
package used to perform factorization: mumps
total: nonzeros=6.47715e+07, allocated nonzeros=6.47715e+07
total number of mallocs used during MatSetValues calls =0
MUMPS run parameters:
SYM (matrix type): 0
PAR (host participation): 1
ICNTL(1) (output for error): 6
ICNTL(2) (output of diagnostic msg): 0
ICNTL(3) (output for global info): 0
ICNTL(4) (level of printing): 0
ICNTL(5) (input mat struct): 0
ICNTL(6) (matrix prescaling): 7
ICNTL(7) (sequentia matrix ordering):7
ICNTL(8) (scalling strategy): 77
ICNTL(10) (max num of refinements): 0
ICNTL(11) (error analysis): 0
ICNTL(12) (efficiency control): 1
ICNTL(13) (efficiency control): 0
ICNTL(14) (percentage of estimated workspace increase): 20
ICNTL(18) (input mat struct): 3
ICNTL(19) (Shur complement info): 0
ICNTL(20) (rhs sparse pattern): 0
ICNTL(21) (somumpstion struct): 1
ICNTL(22) (in-core/out-of-core facility): 0
ICNTL(23) (max size of memory can be allocated locally):0
ICNTL(24) (detection of null pivot rows): 0
ICNTL(25) (computation of a null space basis): 0
ICNTL(26) (Schur options for rhs or solution): 0
ICNTL(27) (experimental parameter): -8
ICNTL(28) (use parallel or sequential ordering): 1
ICNTL(29) (parallel ordering): 0
ICNTL(30) (user-specified set of entries in inv(A)): 0
ICNTL(31) (factors is discarded in the solve phase): 0
ICNTL(33) (compute determinant): 0
CNTL(1) (relative pivoting threshold): 0.01
CNTL(2) (stopping criterion of refinement): 1.49012e-08
CNTL(3) (absomumpste pivoting threshold): 0
CNTL(4) (vamumpse of static pivoting): -1
CNTL(5) (fixation for null pivots): 0
RINFO(1) (local estimated flops for the elimination after analysis):
[0] 1.6957e+10
[1] 1.67964e+10
[2] 1.70959e+10
[3] 1.7073e+10
RINFO(2) (local estimated flops for the assembly after factorization):
[0] 3.99956e+07
[1] 3.71513e+07
[2] 4.43658e+07
[3] 2.95262e+07
RINFO(3) (local estimated flops for the elimination after factorization):
[0] 1.732e+10
[1] 1.7234e+10
[2] 1.70959e+10
[3] 1.62723e+10
INFO(15) (estimated size of (in MB) MUMPS internal data for running numerical factorization):
[0] 243
[1] 237
[2] 256
[3] 212
INFO(16) (size of (in MB) MUMPS internal data used during numerical factorization):
[0] 243
[1] 237
[2] 256
[3] 212
INFO(23) (num of pivots eliminated on this processor after factorization):
[0] 23748
[1] 15630
[2] 20355
[3] 12213
RINFOG(1) (global estimated flops for the elimination after analysis): 6.79222e+10
RINFOG(2) (global estimated flops for the assembly after factorization): 1.51039e+08
RINFOG(3) (global estimated flops for the elimination after factorization): 6.79222e+10
(RINFOG(12) RINFOG(13))*2^INFOG(34) (determinant): (0,0)*(2^0)
INFOG(3) (estimated real workspace for factors on all processors after analysis): 64771542
INFOG(4) (estimated integer workspace for factors on all processors after analysis): 1209207
INFOG(5) (estimated maximum front size in the complete tree): 2562
INFOG(6) (number of nodes in the complete tree): 3832
INFOG(7) (ordering option effectively use after analysis): 5
INFOG(8) (structural symmetry in percent of the permuted matrix after analysis): 100
INFOG(9) (total real/complex workspace to store the matrix factors after factorization): 64771542
INFOG(10) (total integer space store the matrix factors after factorization): 1208659
INFOG(11) (order of largest frontal matrix after factorization): 2562
INFOG(12) (number of off-diagonal pivots): 0
INFOG(13) (number of delayed pivots after factorization): 0
INFOG(14) (number of memory compress after factorization): 0
INFOG(15) (number of steps of iterative refinement after solution): 0
INFOG(16) (estimated size (in MB) of all MUMPS internal data for factorization after analysis: value on the most memory consuming processor): 256
INFOG(17) (estimated size of all MUMPS internal data for factorization after analysis: sum over all processors): 948
INFOG(18) (size of all MUMPS internal data allocated during factorization: value on the most memory consuming processor): 256
INFOG(19) (size of all MUMPS internal data allocated during factorization: sum over all processors): 948
INFOG(20) (estimated number of entries in the factors): 64771542
INFOG(21) (size in MB of memory effectively used during factorization - value on the most memory consuming processor): 212
INFOG(22) (size in MB of memory effectively used during factorization - sum over all processors): 789
INFOG(23) (after analysis: value of ICNTL(6) effectively used): 0
INFOG(24) (after analysis: value of ICNTL(12) effectively used): 1
INFOG(25) (after factorization: number of pivots modified by static pivoting): 0
INFOG(28) (after factorization: number of null pivots encountered): 0
INFOG(29) (after factorization: effective number of entries in the factors (sum over all processors)): 58207698
INFOG(30, 31) (after solution: size in Mbytes of memory used during solution phase): 229, 844
INFOG(32) (after analysis: type of analysis done): 1
INFOG(33) (value used for ICNTL(8)): 7
INFOG(34) (exponent of the determinant if determinant is requested): 0
linear system matrix = precond matrix:
Mat Object: (fieldsplit_0_) 4 MPI processes
type: mpiaij
rows=71946, cols=71946
total: nonzeros=2.96023e+06, allocated nonzeros=2.96023e+06
total number of mallocs used during MatSetValues calls =0
using I-node (on process 0) routines: found 6259 nodes, limit used is 5
A01
Mat Object: 4 MPI processes
type: mpiaij
rows=71946, cols=3183
total: nonzeros=19098, allocated nonzeros=19098
total number of mallocs used during MatSetValues calls =0
using I-node (on process 0) routines: found 3757 nodes, limit used is 5
linear system matrix = precond matrix:
Mat Object: 4 MPI processes
type: mpiaij
rows=75129, cols=75129
total: nonzeros=3.0016e+06, allocated nonzeros=1.26217e+07
total number of mallocs used during MatSetValues calls =0
using I-node (on process 0) routines: found 6259 nodes, limit used is 5
Recovering stress ...
Cleaning up ...
Finished
0.00user 0.02system 0:26.40elapsed 0%CPU (0avgtext+0avgdata 17056maxresident)k
1168inputs+0outputs (6major+1797minor)pagefaults 0swaps
salloc: Relinquishing job allocation 78969
salloc: Job allocation 78969 has been revoked.
More information about the petsc-users
mailing list