[petsc-users] Debugging failed solve (what's an acceptable upper bound to the condition number?)
Alex Lindsay
adlinds3 at ncsu.edu
Mon Nov 23 14:44:33 CST 2015
I've found that with a "full" variable set, I can get convergence.
However, if I remove two of my variables (the other variables have no
dependence on the variables that I remove; the coupling is one-way),
then I no longer get convergence. I've attached logs of one time-step
for both the converged and non-converged cases.
On 11/23/2015 01:29 PM, Alex Lindsay wrote:
> On 11/20/2015 02:33 PM, Jed Brown wrote:
>> Alex Lindsay <adlinds3 at ncsu.edu> writes:
>>> I'm almost ashamed to share my condition number because I'm sure it
>>> must
>>> be absurdly high. Without applying -ksp_diagonal_scale and
>>> -ksp_diagonal_scale_fix, the condition number is around 1e25. When I do
>>> apply those two parameters, the condition number is reduced to 1e17.
>>> Even after scaling all my variable residuals so that they were all on
>>> the order of unity (a suggestion on the Moose list), I still have a
>>> condition number of 1e12.
>> Double precision provides 16 digits of accuracy in the best case. When
>> you finite difference, the accuracy is reduced to 8 digits if the
>> differencing parameter is chosen optimally. With the condition numbers
>> you're reporting, your matrix is singular up to available precision.
>>
>>> I have no experience with condition numbers, but knowing that perfect
>>> condition number is unity, 1e12 seems unacceptable. What's an
>>> acceptable upper limit on the condition number? Is it problem
>>> dependent? Having already tried scaling the individual variable
>>> residuals, I'm not exactly sure what my next method would be for
>>> trying to reduce the condition number.
>> Singular operators are often caused by incorrect boundary conditions.
>> You should try a small and simple version of your problem and find out
>> why it's producing a singular (or so close to singular we can't tell)
>> operator.
> Could large variable values also create singular operators? I'm
> essentially solving an advection-diffusion-reaction problem for
> several species where the advection is driven by an electric field.
> The species concentrations are in a logarithmic form such that the
> true concentration is given by exp(u). With my current units (# of
> particles / m^3) exp(u) is anywhere from 1e13 to 1e20, and thus the
> initial residuals are probably on the same order of magnitude. After
> I've assembled the total residual for each variable and before the
> residual is passed to the solver, I apply scaling to the residuals
> such that the sum of the variable residuals is around 1e3. But perhaps
> I lose some accuracy during the residual assembly process?
>
> I'm equating "incorrect" boundary conditions to "unphysical" or
> "unrealistic" boundary conditions. Hopefully that's fair.
-------------- next part --------------
Time Step 1, time = 6.4e-11
dt = 6.4e-11
|residual|_2 of individual variables:
potential: 1250
potentialliq: 0
em: 5.73385e-12
emliq: 3.08187e-10
Arp: 5.73385e-12
0 Nonlinear |R| = 1.250000e+03
0 KSP unpreconditioned resid norm 1.250000000000e+03 true resid norm 1.250000000000e+03 ||r(i)||/||b|| 1.000000000000e+00
1 KSP unpreconditioned resid norm 1.250000000000e+03 true resid norm 1.250000000000e+03 ||r(i)||/||b|| 1.000000000000e+00
2 KSP unpreconditioned resid norm 2.182679427760e+02 true resid norm 7.819529716916e+08 ||r(i)||/||b|| 6.255623773533e+05
3 KSP unpreconditioned resid norm 1.011652745364e+02 true resid norm 4.461678857470e+01 ||r(i)||/||b|| 3.569343085976e-02
4 KSP unpreconditioned resid norm 8.125623676015e+00 true resid norm 8.053940519499e+08 ||r(i)||/||b|| 6.443152415599e+05
5 KSP unpreconditioned resid norm 5.805247155944e+00 true resid norm 8.054105876447e+08 ||r(i)||/||b|| 6.443284701157e+05
6 KSP unpreconditioned resid norm 4.756488143441e+00 true resid norm 8.054162537433e+08 ||r(i)||/||b|| 6.443330029946e+05
7 KSP unpreconditioned resid norm 4.126450902175e+00 true resid norm 8.054191165808e+08 ||r(i)||/||b|| 6.443352932646e+05
8 KSP unpreconditioned resid norm 3.694696196953e+00 true resid norm 8.054208439360e+08 ||r(i)||/||b|| 6.443366751488e+05
9 KSP unpreconditioned resid norm 3.375152117403e+00 true resid norm 8.054219995564e+08 ||r(i)||/||b|| 6.443375996451e+05
10 KSP unpreconditioned resid norm 3.126354693526e+00 true resid norm 8.054228269923e+08 ||r(i)||/||b|| 6.443382615939e+05
Linear solve converged due to CONVERGED_ITS iterations 10
Line search: gnorm after quadratic fit 2.835243959849e+07
Line search: Cubic step no good, shrinking lambda, current gnorm 2.692665331241e+06 lambda=1.0000000000000002e-02
Line search: Cubic step no good, shrinking lambda, current gnorm 2.678836859543e+05 lambda=1.0000000000000002e-03
Line search: Cubic step no good, shrinking lambda, current gnorm 2.680339963946e+04 lambda=1.0000000000000003e-04
Line search: Cubic step no good, shrinking lambda, current gnorm 2.954721111888e+03 lambda=1.0000000000000004e-05
Line search: Cubic step no good, shrinking lambda, current gnorm 1.278349531522e+03 lambda=1.0000000000000004e-06
Line search: Cubic step no good, shrinking lambda, current gnorm 1.250287112903e+03 lambda=1.0000000000000005e-07
Line search: Cubic step no good, shrinking lambda, current gnorm 1.250002830553e+03 lambda=1.0000000000000005e-08
Line search: Cubic step no good, shrinking lambda, current gnorm 1.250000030832e+03 lambda=1.0000000000000005e-09
Line search: Cubic step no good, shrinking lambda, current gnorm 1.250000000426e+03 lambda=1.0000000000000006e-10
Line search: Cubic step no good, shrinking lambda, current gnorm 1.250000000481e+03 lambda=1.0922440478434965e-11
Line search: Cubically determined step, current gnorm 1.249999999999e+03 lambda=1.0922440478434966e-12
|residual|_2 of individual variables:
potential: 1250
potentialliq: 8.09094e-26
em: 1.07274e-11
emliq: 3.08187e-10
Arp: 5.73508e-12
1 Nonlinear |R| = 1.250000e+03
0 KSP unpreconditioned resid norm 1.249999999999e+03 true resid norm 1.249999999999e+03 ||r(i)||/||b|| 1.000000000000e+00
1 KSP unpreconditioned resid norm 1.249999999999e+03 true resid norm 1.249999999999e+03 ||r(i)||/||b|| 1.000000000000e+00
2 KSP unpreconditioned resid norm 2.183465611429e+02 true resid norm 8.427616262825e+08 ||r(i)||/||b|| 6.742093010267e+05
3 KSP unpreconditioned resid norm 1.599839731484e+02 true resid norm 8.129044290551e+08 ||r(i)||/||b|| 6.503235432448e+05
4 KSP unpreconditioned resid norm 1.320218637568e+02 true resid norm 8.182745495747e+08 ||r(i)||/||b|| 6.546196396604e+05
5 KSP unpreconditioned resid norm 1.149512781910e+02 true resid norm 1.687718402197e+09 ||r(i)||/||b|| 1.350174721759e+06
6 KSP unpreconditioned resid norm 1.031509670396e+02 true resid norm 1.645498956462e+09 ||r(i)||/||b|| 1.316399165171e+06
7 KSP unpreconditioned resid norm 9.436932085660e+01 true resid norm 1.647788157739e+09 ||r(i)||/||b|| 1.318230526193e+06
8 KSP unpreconditioned resid norm 8.750587817698e+01 true resid norm 1.649436300901e+09 ||r(i)||/||b|| 1.319549040722e+06
9 KSP unpreconditioned resid norm 8.195066818457e+01 true resid norm 1.696531801097e+09 ||r(i)||/||b|| 1.357225440879e+06
10 KSP unpreconditioned resid norm 7.733475365117e+01 true resid norm 1.651650881484e+09 ||r(i)||/||b|| 1.321320705189e+06
Linear solve converged due to CONVERGED_ITS iterations 10
Line search: gnorm after quadratic fit 1.594871627961e+15
Line search: Cubic step no good, shrinking lambda, current gnorm 1.130742402617e+08 lambda=5.0000000000000003e-02
Line search: Cubic step no good, shrinking lambda, current gnorm 3.597173061575e+07 lambda=2.5000000000000001e-02
Line search: Cubic step no good, shrinking lambda, current gnorm 7.589836896768e+06 lambda=5.3308924183071470e-03
Line search: Cubic step no good, shrinking lambda, current gnorm 7.570310699594e+05 lambda=5.3308924183071468e-04
Line search: Cubic step no good, shrinking lambda, current gnorm 7.569384032198e+04 lambda=5.3308924183071473e-05
Line search: Cubic step no good, shrinking lambda, current gnorm 7.670683394580e+03 lambda=5.3308924183071480e-06
Line search: Cubic step no good, shrinking lambda, current gnorm 1.461254084896e+03 lambda=5.3308924183071482e-07
Line search: Cubic step no good, shrinking lambda, current gnorm 1.252290220827e+03 lambda=5.3308924183071487e-08
Line search: Cubic step no good, shrinking lambda, current gnorm 1.250022917003e+03 lambda=5.3308924183071492e-09
Line search: Cubic step no good, shrinking lambda, current gnorm 1.250000221997e+03 lambda=5.3308924183071499e-10
Line search: Cubic step no good, shrinking lambda, current gnorm 1.250000001910e+03 lambda=5.3308924183071504e-11
Line search: Cubically determined step, current gnorm 1.249999999992e+03 lambda=5.3308924183071507e-12
|residual|_2 of individual variables:
potential: 1250
potentialliq: 7.42554e-26
em: 1.22453e-09
emliq: 3.08187e-10
Arp: 1.38635e-11
2 Nonlinear |R| = 1.250000e+03
0 KSP unpreconditioned resid norm 1.249999999992e+03 true resid norm 1.249999999992e+03 ||r(i)||/||b|| 1.000000000000e+00
1 KSP unpreconditioned resid norm 1.249999999992e+03 true resid norm 1.249999999992e+03 ||r(i)||/||b|| 1.000000000000e+00
2 KSP unpreconditioned resid norm 2.182596482393e+02 true resid norm 8.253967021927e+08 ||r(i)||/||b|| 6.603173617584e+05
3 KSP unpreconditioned resid norm 1.947039771924e+01 true resid norm 2.075931294022e+01 ||r(i)||/||b|| 1.660745035228e-02
4 KSP unpreconditioned resid norm 1.436461464059e+01 true resid norm 2.083785869768e+01 ||r(i)||/||b|| 1.667028695825e-02
5 KSP unpreconditioned resid norm 1.190526501926e+01 true resid norm 2.086697733900e+01 ||r(i)||/||b|| 1.669358187131e-02
6 KSP unpreconditioned resid norm 1.038939480233e+01 true resid norm 2.088215650170e+01 ||r(i)||/||b|| 1.670572520146e-02
7 KSP unpreconditioned resid norm 9.335802160607e+00 true resid norm 2.089153252677e+01 ||r(i)||/||b|| 1.671322602152e-02
8 KSP unpreconditioned resid norm 8.549016490263e+00 true resid norm 2.089795904301e+01 ||r(i)||/||b|| 1.671836723451e-02
9 KSP unpreconditioned resid norm 7.932601421080e+00 true resid norm 2.090238701435e+01 ||r(i)||/||b|| 1.672190961159e-02
10 KSP unpreconditioned resid norm 7.432799238568e+00 true resid norm 2.090615112323e+01 ||r(i)||/||b|| 1.672492089869e-02
Linear solve converged due to CONVERGED_ITS iterations 10
Line search: gnorm after quadratic fit 1.142198665507e+07
Line search: Cubic step no good, shrinking lambda, current gnorm 5.550401765706e+06 lambda=5.0000000000000003e-02
Line search: Cubic step no good, shrinking lambda, current gnorm 5.409748027962e+05 lambda=5.0000000000000010e-03
Line search: Cubic step no good, shrinking lambda, current gnorm 5.397311265735e+04 lambda=5.0000000000000012e-04
Line search: Cubic step no good, shrinking lambda, current gnorm 5.537394894252e+03 lambda=5.0000000000000016e-05
Line search: Cubic step no good, shrinking lambda, current gnorm 1.361427787561e+03 lambda=5.0000000000000021e-06
Line search: Cubic step no good, shrinking lambda, current gnorm 1.251161609975e+03 lambda=5.0000000000000019e-07
Line search: Cubic step no good, shrinking lambda, current gnorm 1.250011611567e+03 lambda=5.0000000000000024e-08
Line search: Cubic step no good, shrinking lambda, current gnorm 1.250000107464e+03 lambda=5.0000000000000026e-09
Line search: Cubic step no good, shrinking lambda, current gnorm 1.250000001345e+03 lambda=5.0000000000000024e-10
Line search: Cubically determined step, current gnorm 1.249999999897e+03 lambda=7.6233253179964001e-11
|residual|_2 of individual variables:
potential: 1250
potentialliq: 3.15588e-24
em: 2.84656e-09
emliq: 3.08187e-10
Arp: 2.99446e-11
3 Nonlinear |R| = 1.250000e+03
0 KSP unpreconditioned resid norm 1.249999999897e+03 true resid norm 1.249999999897e+03 ||r(i)||/||b|| 1.000000000000e+00
1 KSP unpreconditioned resid norm 1.249999999897e+03 true resid norm 1.249999999897e+03 ||r(i)||/||b|| 1.000000000000e+00
2 KSP unpreconditioned resid norm 2.181586411076e+02 true resid norm 7.819766499349e+08 ||r(i)||/||b|| 6.255813199996e+05
3 KSP unpreconditioned resid norm 1.854352309003e+01 true resid norm 2.003012339454e+01 ||r(i)||/||b|| 1.602409871696e-02
4 KSP unpreconditioned resid norm 1.362119502191e+01 true resid norm 8.053751742035e+08 ||r(i)||/||b|| 6.443001394160e+05
5 KSP unpreconditioned resid norm 1.127240232708e+01 true resid norm 8.054050843536e+08 ||r(i)||/||b|| 6.443240675361e+05
6 KSP unpreconditioned resid norm 9.828809807849e+00 true resid norm 8.054206656511e+08 ||r(i)||/||b|| 6.443365325741e+05
7 KSP unpreconditioned resid norm 8.827500453712e+00 true resid norm 8.054302209816e+08 ||r(i)||/||b|| 6.443441768386e+05
8 KSP unpreconditioned resid norm 8.080717820462e+00 true resid norm 8.054366796615e+08 ||r(i)||/||b|| 6.443493437824e+05
9 KSP unpreconditioned resid norm 7.496176409465e+00 true resid norm 8.054413371079e+08 ||r(i)||/||b|| 6.443530697396e+05
10 KSP unpreconditioned resid norm 7.022528708078e+00 true resid norm 8.054448546154e+08 ||r(i)||/||b|| 6.443558837456e+05
Linear solve converged due to CONVERGED_ITS iterations 10
Line search: gnorm after quadratic fit 2.511396331623e+07
Line search: Cubic step no good, shrinking lambda, current gnorm 1.220228888591e+07 lambda=5.0000000000000003e-02
Line search: Cubic step no good, shrinking lambda, current gnorm 1.189164603546e+06 lambda=5.0000000000000010e-03
Line search: Cubic step no good, shrinking lambda, current gnorm 1.186167166565e+05 lambda=5.0000000000000012e-04
Line search: Cubic step no good, shrinking lambda, current gnorm 1.192365831756e+04 lambda=5.0000000000000016e-05
Line search: Cubic step no good, shrinking lambda, current gnorm 1.722937964326e+03 lambda=5.0000000000000021e-06
Line search: Cubic step no good, shrinking lambda, current gnorm 1.255613303497e+03 lambda=5.0000000000000019e-07
Line search: Cubic step no good, shrinking lambda, current gnorm 1.250056305277e+03 lambda=5.0000000000000024e-08
Line search: Cubic step no good, shrinking lambda, current gnorm 1.250000578304e+03 lambda=5.0000000000000026e-09
Line search: Cubic step no good, shrinking lambda, current gnorm 1.250000003888e+03 lambda=5.0000000000000024e-10
Line search: Cubically determined step, current gnorm 1.249999999834e+03 lambda=5.0000000000000028e-11
|residual|_2 of individual variables:
potential: 1250
potentialliq: 4.64775e-24
em: 3.8499e-09
emliq: 3.08187e-10
Arp: 4.01495e-11
4 Nonlinear |R| = 1.250000e+03
0 KSP unpreconditioned resid norm 1.249999999834e+03 true resid norm 1.249999999834e+03 ||r(i)||/||b|| 1.000000000000e+00
1 KSP unpreconditioned resid norm 1.249999999834e+03 true resid norm 1.249999999834e+03 ||r(i)||/||b|| 1.000000000000e+00
2 KSP unpreconditioned resid norm 2.180911377684e+02 true resid norm 7.819912747111e+08 ||r(i)||/||b|| 6.255930198518e+05
3 KSP unpreconditioned resid norm 1.858951272351e+01 true resid norm 2.008399501169e+01 ||r(i)||/||b|| 1.606719601148e-02
4 KSP unpreconditioned resid norm 1.366214280005e+01 true resid norm 8.054474955190e+08 ||r(i)||/||b|| 6.443579965007e+05
5 KSP unpreconditioned resid norm 1.130864341976e+01 true resid norm 2.012535124271e+01 ||r(i)||/||b|| 1.610028099630e-02
6 KSP unpreconditioned resid norm 9.862166661146e+00 true resid norm 2.013215117729e+01 ||r(i)||/||b|| 1.610572094397e-02
7 KSP unpreconditioned resid norm 8.858424099310e+00 true resid norm 8.055390402102e+08 ||r(i)||/||b|| 6.444312322537e+05
8 KSP unpreconditioned resid norm 8.109622559839e+00 true resid norm 2.013927312439e+01 ||r(i)||/||b|| 1.611141850165e-02
9 KSP unpreconditioned resid norm 7.523390433943e+00 true resid norm 2.014118037135e+01 ||r(i)||/||b|| 1.611294429922e-02
10 KSP unpreconditioned resid norm 7.048305660327e+00 true resid norm 2.014299220373e+01 ||r(i)||/||b|| 1.611439376512e-02
Linear solve converged due to CONVERGED_ITS iterations 10
Line search: gnorm after quadratic fit 2.391135555969e+07
Line search: Cubic step no good, shrinking lambda, current gnorm 1.161799021499e+07 lambda=5.0000000000000003e-02
Line search: Cubic step no good, shrinking lambda, current gnorm 1.132224103053e+06 lambda=5.0000000000000010e-03
Line search: Cubic step no good, shrinking lambda, current gnorm 1.129376510639e+05 lambda=5.0000000000000012e-04
Line search: Cubic step no good, shrinking lambda, current gnorm 1.135915546728e+04 lambda=5.0000000000000016e-05
Line search: Cubic step no good, shrinking lambda, current gnorm 1.684379468954e+03 lambda=5.0000000000000021e-06
Line search: Cubic step no good, shrinking lambda, current gnorm 1.255085688155e+03 lambda=5.0000000000000019e-07
Line search: Cubic step no good, shrinking lambda, current gnorm 1.250051037016e+03 lambda=5.0000000000000024e-08
Line search: Cubic step no good, shrinking lambda, current gnorm 1.250000478960e+03 lambda=5.0000000000000026e-09
Line search: Cubic step no good, shrinking lambda, current gnorm 1.250000003825e+03 lambda=5.0000000000000024e-10
Line search: Cubically determined step, current gnorm 1.249999999772e+03 lambda=5.0000000000000028e-11
|residual|_2 of individual variables:
potential: 1250
potentialliq: 5.64851e-24
em: 4.85699e-09
emliq: 3.08187e-10
Arp: 5.04572e-11
5 Nonlinear |R| = 1.250000e+03
0 KSP unpreconditioned resid norm 1.249999999772e+03 true resid norm 1.249999999772e+03 ||r(i)||/||b|| 1.000000000000e+00
1 KSP unpreconditioned resid norm 1.249999999772e+03 true resid norm 1.249999999772e+03 ||r(i)||/||b|| 1.000000000000e+00
2 KSP unpreconditioned resid norm 2.180245879330e+02 true resid norm 7.820056805196e+08 ||r(i)||/||b|| 6.256045445299e+05
3 KSP unpreconditioned resid norm 1.946205886553e+01 true resid norm 2.109700215024e+01 ||r(i)||/||b|| 1.687760172328e-02
4 KSP unpreconditioned resid norm 1.435814972977e+01 true resid norm 2.117498811000e+01 ||r(i)||/||b|| 1.693999049110e-02
5 KSP unpreconditioned resid norm 1.189985941072e+01 true resid norm 2.120350622398e+01 ||r(i)||/||b|| 1.696280498229e-02
6 KSP unpreconditioned resid norm 1.038474534940e+01 true resid norm 2.121871510654e+01 ||r(i)||/||b|| 1.697497208834e-02
7 KSP unpreconditioned resid norm 9.331644867005e+00 true resid norm 2.122782552136e+01 ||r(i)||/||b|| 1.698226042019e-02
8 KSP unpreconditioned resid norm 8.545199194359e+00 true resid norm 2.123399509830e+01 ||r(i)||/||b|| 1.698719608174e-02
9 KSP unpreconditioned resid norm 7.929057651795e+00 true resid norm 2.123866776923e+01 ||r(i)||/||b|| 1.699093421849e-02
10 KSP unpreconditioned resid norm 7.429473836034e+00 true resid norm 2.124192017166e+01 ||r(i)||/||b|| 1.699353614043e-02
Linear solve converged due to CONVERGED_ITS iterations 10
Line search: gnorm after quadratic fit 1.137552646194e+07
Line search: Cubic step no good, shrinking lambda, current gnorm 5.527796109975e+06 lambda=5.0000000000000003e-02
Line search: Cubic step no good, shrinking lambda, current gnorm 5.387690215403e+05 lambda=5.0000000000000010e-03
Line search: Cubic step no good, shrinking lambda, current gnorm 5.375314646450e+04 lambda=5.0000000000000012e-04
Line search: Cubic step no good, shrinking lambda, current gnorm 5.515964934581e+03 lambda=5.0000000000000016e-05
Line search: Cubic step no good, shrinking lambda, current gnorm 1.360556335854e+03 lambda=5.0000000000000021e-06
Line search: Cubic step no good, shrinking lambda, current gnorm 1.251153963289e+03 lambda=5.0000000000000019e-07
Line search: Cubic step no good, shrinking lambda, current gnorm 1.250011459902e+03 lambda=5.0000000000000024e-08
Line search: Cubic step no good, shrinking lambda, current gnorm 1.250000107244e+03 lambda=5.0000000000000026e-09
Line search: Cubic step no good, shrinking lambda, current gnorm 1.250000001236e+03 lambda=5.0000000000000024e-10
Line search: Cubic step no good, shrinking lambda, current gnorm 1.250000000233e+03 lambda=7.1940637925727652e-11
Line search: Cubically determined step, current gnorm 1.249999999763e+03 lambda=7.1940637925727654e-12
|residual|_2 of individual variables:
potential: 1250
potentialliq: 6.10417e-24
em: 5.0098e-09
emliq: 3.08187e-10
Arp: 5.20276e-11
6 Nonlinear |R| = 1.250000e+03
0 KSP unpreconditioned resid norm 1.249999999763e+03 true resid norm 1.249999999763e+03 ||r(i)||/||b|| 1.000000000000e+00
1 KSP unpreconditioned resid norm 1.249999999763e+03 true resid norm 1.249999999763e+03 ||r(i)||/||b|| 1.000000000000e+00
2 KSP unpreconditioned resid norm 2.180044693796e+02 true resid norm 7.802198484180e+08 ||r(i)||/||b|| 6.241758788529e+05
3 KSP unpreconditioned resid norm 1.732135088416e+02 true resid norm 1.563215239340e+09 ||r(i)||/||b|| 1.250572191709e+06
4 KSP unpreconditioned resid norm 1.061485630366e+02 true resid norm 1.577335315211e+09 ||r(i)||/||b|| 1.261868252409e+06
5 KSP unpreconditioned resid norm 8.357822500760e+01 true resid norm 1.580790144294e+09 ||r(i)||/||b|| 1.264632115675e+06
6 KSP unpreconditioned resid norm 7.114455942454e+01 true resid norm 1.582376132218e+09 ||r(i)||/||b|| 1.265900906015e+06
7 KSP unpreconditioned resid norm 6.299675123208e+01 true resid norm 1.583286480061e+09 ||r(i)||/||b|| 1.266629184289e+06
8 KSP unpreconditioned resid norm 5.712987155691e+01 true resid norm 1.583877062081e+09 ||r(i)||/||b|| 1.267101649905e+06
9 KSP unpreconditioned resid norm 5.264617509968e+01 true resid norm 2.420444849004e+09 ||r(i)||/||b|| 1.936355879571e+06
10 KSP unpreconditioned resid norm 4.907562644314e+01 true resid norm 2.420913038054e+09 ||r(i)||/||b|| 1.936730430811e+06
Linear solve converged due to CONVERGED_ITS iterations 10
Line search: gnorm after quadratic fit 2.034499848293e+08
Line search: Cubic step no good, shrinking lambda, current gnorm 9.871552786457e+07 lambda=5.0000000000000003e-02
Line search: Cubic step no good, shrinking lambda, current gnorm 9.608320969172e+06 lambda=5.0000000000000010e-03
Line search: Cubic step no good, shrinking lambda, current gnorm 9.582394863587e+05 lambda=5.0000000000000012e-04
Line search: Cubic step no good, shrinking lambda, current gnorm 9.580610519546e+04 lambda=5.0000000000000016e-05
Line search: Cubic step no good, shrinking lambda, current gnorm 9.660746121724e+03 lambda=5.0000000000000021e-06
Line search: Cubic step no good, shrinking lambda, current gnorm 1.574849255364e+03 lambda=5.0000000000000019e-07
Line search: Cubic step no good, shrinking lambda, current gnorm 1.253664335218e+03 lambda=5.0000000000000024e-08
Line search: Cubic step no good, shrinking lambda, current gnorm 1.250036823987e+03 lambda=5.0000000000000026e-09
Line search: Cubic step no good, shrinking lambda, current gnorm 1.250000366999e+03 lambda=5.0000000000000024e-10
Line search: Cubic step no good, shrinking lambda, current gnorm 1.250000001678e+03 lambda=5.0000000000000028e-11
Line search: Cubically determined step, current gnorm 1.249999999756e+03 lambda=5.0000000000000029e-12
|residual|_2 of individual variables:
potential: 1250
potentialliq: 6.45789e-24
em: 4.39941e-09
emliq: 3.08187e-10
Arp: 4.57523e-11
7 Nonlinear |R| = 1.250000e+03
0 KSP unpreconditioned resid norm 1.249999999756e+03 true resid norm 1.249999999756e+03 ||r(i)||/||b|| 1.000000000000e+00
1 KSP unpreconditioned resid norm 1.249999999756e+03 true resid norm 1.249999999756e+03 ||r(i)||/||b|| 1.000000000000e+00
2 KSP unpreconditioned resid norm 2.180081873135e+02 true resid norm 7.820092290130e+08 ||r(i)||/||b|| 6.256073833323e+05
3 KSP unpreconditioned resid norm 1.946129832795e+01 true resid norm 2.109616756847e+01 ||r(i)||/||b|| 1.687693405807e-02
4 KSP unpreconditioned resid norm 1.435756568789e+01 true resid norm 2.117410765832e+01 ||r(i)||/||b|| 1.693928612996e-02
5 KSP unpreconditioned resid norm 1.189936827192e+01 true resid norm 2.120240607210e+01 ||r(i)||/||b|| 1.696192486099e-02
6 KSP unpreconditioned resid norm 1.038420125742e+01 true resid norm 2.121758216460e+01 ||r(i)||/||b|| 1.697406573498e-02
7 KSP unpreconditioned resid norm 9.331121524454e+00 true resid norm 2.122668148932e+01 ||r(i)||/||b|| 1.698134519476e-02
8 KSP unpreconditioned resid norm 8.544713991399e+00 true resid norm 2.123285231780e+01 ||r(i)||/||b|| 1.698628185755e-02
9 KSP unpreconditioned resid norm 7.928603436376e+00 true resid norm 2.123726676664e+01 ||r(i)||/||b|| 1.698981341662e-02
10 KSP unpreconditioned resid norm 7.429045402840e+00 true resid norm 2.124077912149e+01 ||r(i)||/||b|| 1.699262330050e-02
Linear solve converged due to CONVERGED_ITS iterations 10
Line search: gnorm after quadratic fit 1.149326227191e+07
Line search: Cubic step no good, shrinking lambda, current gnorm 5.585045811793e+06 lambda=5.0000000000000003e-02
Line search: Cubic step no good, shrinking lambda, current gnorm 5.443521396286e+05 lambda=5.0000000000000010e-03
Line search: Cubic step no good, shrinking lambda, current gnorm 5.430987061393e+04 lambda=5.0000000000000012e-04
Line search: Cubic step no good, shrinking lambda, current gnorm 5.570191486810e+03 lambda=5.0000000000000016e-05
Line search: Cubic step no good, shrinking lambda, current gnorm 1.362768603478e+03 lambda=5.0000000000000021e-06
Line search: Cubic step no good, shrinking lambda, current gnorm 1.251176978273e+03 lambda=5.0000000000000019e-07
Line search: Cubic step no good, shrinking lambda, current gnorm 1.250011763766e+03 lambda=5.0000000000000024e-08
Line search: Cubic step no good, shrinking lambda, current gnorm 1.250000107229e+03 lambda=5.0000000000000026e-09
Line search: Cubic step no good, shrinking lambda, current gnorm 1.250000001109e+03 lambda=5.0000000000000024e-10
Line search: Cubically determined step, current gnorm 1.249999999661e+03 lambda=7.6233447816453448e-11
|residual|_2 of individual variables:
potential: 1250
potentialliq: 9.08761e-24
em: 6.01857e-09
emliq: 3.08187e-10
Arp: 6.23999e-11
8 Nonlinear |R| = 1.250000e+03
0 KSP unpreconditioned resid norm 1.249999999661e+03 true resid norm 1.249999999661e+03 ||r(i)||/||b|| 1.000000000000e+00
1 KSP unpreconditioned resid norm 1.249999999661e+03 true resid norm 1.249999999661e+03 ||r(i)||/||b|| 1.000000000000e+00
2 KSP unpreconditioned resid norm 2.179064067278e+02 true resid norm 7.820312571339e+08 ||r(i)||/||b|| 6.256250058767e+05
3 KSP unpreconditioned resid norm 1.945676272501e+01 true resid norm 2.109203237667e+01 ||r(i)||/||b|| 1.687362590591e-02
4 KSP unpreconditioned resid norm 1.435481284529e+01 true resid norm 2.117033993758e+01 ||r(i)||/||b|| 1.693627195466e-02
5 KSP unpreconditioned resid norm 1.189760945248e+01 true resid norm 8.068396959495e+08 ||r(i)||/||b|| 6.454717569346e+05
6 KSP unpreconditioned resid norm 1.038290484991e+01 true resid norm 2.121519397161e+01 ||r(i)||/||b|| 1.697215518189e-02
7 KSP unpreconditioned resid norm 9.330088374361e+00 true resid norm 2.122447838136e+01 ||r(i)||/||b|| 1.697958270969e-02
8 KSP unpreconditioned resid norm 8.543849617441e+00 true resid norm 2.123076797725e+01 ||r(i)||/||b|| 1.698461438641e-02
9 KSP unpreconditioned resid norm 7.927856131237e+00 true resid norm 2.123553406701e+01 ||r(i)||/||b|| 1.698842725821e-02
10 KSP unpreconditioned resid norm 7.428383953238e+00 true resid norm 2.123911387058e+01 ||r(i)||/||b|| 1.699129110107e-02
Linear solve converged due to CONVERGED_ITS iterations 10
Line search: gnorm after quadratic fit 1.044209236144e+07
Line search: Cubic step no good, shrinking lambda, current gnorm 5.074235494424e+06 lambda=5.0000000000000003e-02
Line search: Cubic step no good, shrinking lambda, current gnorm 4.945654291498e+05 lambda=5.0000000000000010e-03
Line search: Cubic step no good, shrinking lambda, current gnorm 4.934540687239e+04 lambda=5.0000000000000012e-04
Line search: Cubic step no good, shrinking lambda, current gnorm 5.087618901961e+03 lambda=5.0000000000000016e-05
Line search: Cubic step no good, shrinking lambda, current gnorm 1.343763724232e+03 lambda=5.0000000000000021e-06
Line search: Cubic step no good, shrinking lambda, current gnorm 1.250972379117e+03 lambda=5.0000000000000019e-07
Line search: Cubic step no good, shrinking lambda, current gnorm 1.250009557978e+03 lambda=5.0000000000000024e-08
Line search: Cubic step no good, shrinking lambda, current gnorm 1.250000092632e+03 lambda=5.0000000000000026e-09
Line search: Cubic step no good, shrinking lambda, current gnorm 1.250000001014e+03 lambda=5.0000000000000024e-10
Line search: Cubically determined step, current gnorm 1.249999999566e+03 lambda=7.5768409431551880e-11
|residual|_2 of individual variables:
potential: 1250
potentialliq: 1.95779e-23
em: 7.62773e-09
emliq: 3.08187e-10
Arp: 7.89936e-11
9 Nonlinear |R| = 1.250000e+03
0 KSP unpreconditioned resid norm 1.249999999566e+03 true resid norm 1.249999999566e+03 ||r(i)||/||b|| 1.000000000000e+00
1 KSP unpreconditioned resid norm 1.249999999566e+03 true resid norm 1.249999999566e+03 ||r(i)||/||b|| 1.000000000000e+00
2 KSP unpreconditioned resid norm 2.177996374504e+02 true resid norm 2.177992926509e+02 ||r(i)||/||b|| 1.742394341812e-01
3 KSP unpreconditioned resid norm 7.296868244457e+01 true resid norm 8.007747954267e+08 ||r(i)||/||b|| 6.406198365636e+05
4 KSP unpreconditioned resid norm 3.789417960370e+01 true resid norm 4.444792493019e+01 ||r(i)||/||b|| 3.555833995648e-02
5 KSP unpreconditioned resid norm 3.401642135341e+01 true resid norm 5.080285249605e+01 ||r(i)||/||b|| 4.064228201094e-02
6 KSP unpreconditioned resid norm 2.788071299035e+01 true resid norm 5.009980145722e+01 ||r(i)||/||b|| 4.007984117968e-02
7 KSP unpreconditioned resid norm 2.415622001384e+01 true resid norm 4.966990407429e+01 ||r(i)||/||b|| 3.973592327321e-02
8 KSP unpreconditioned resid norm 2.124505872401e+01 true resid norm 8.043149971256e+08 ||r(i)||/||b|| 6.434519979236e+05
9 KSP unpreconditioned resid norm 1.936316058299e+01 true resid norm 8.043550613379e+08 ||r(i)||/||b|| 6.434840492935e+05
10 KSP unpreconditioned resid norm 1.790786512505e+01 true resid norm 8.043834705860e+08 ||r(i)||/||b|| 6.435067766920e+05
Linear solve converged due to CONVERGED_ITS iterations 10
Line search: gnorm after quadratic fit 4.565500932636e+07
Line search: Cubic step no good, shrinking lambda, current gnorm 2.218074635601e+07 lambda=5.0000000000000003e-02
Line search: Cubic step no good, shrinking lambda, current gnorm 2.161434937923e+06 lambda=5.0000000000000010e-03
Line search: Cubic step no good, shrinking lambda, current gnorm 2.155886905194e+05 lambda=5.0000000000000012e-04
Line search: Cubic step no good, shrinking lambda, current gnorm 2.158915141138e+04 lambda=5.0000000000000016e-05
Line search: Cubic step no good, shrinking lambda, current gnorm 2.491477605019e+03 lambda=5.0000000000000021e-06
Line search: Cubic step no good, shrinking lambda, current gnorm 1.268443119324e+03 lambda=5.0000000000000019e-07
Line search: Cubic step no good, shrinking lambda, current gnorm 1.250186068288e+03 lambda=5.0000000000000024e-08
Line search: Cubic step no good, shrinking lambda, current gnorm 1.250001812872e+03 lambda=5.0000000000000026e-09
Line search: Cubic step no good, shrinking lambda, current gnorm 1.250000017071e+03 lambda=5.0000000000000024e-10
Line search: Cubically determined step, current gnorm 1.249999999504e+03 lambda=5.0000000000000028e-11
|residual|_2 of individual variables:
potential: 1250
potentialliq: 1.38731e-23
em: 5.22255e-09
emliq: 3.08187e-10
Arp: 5.42863e-11
10 Nonlinear |R| = 1.250000e+03
0 KSP unpreconditioned resid norm 1.249999999504e+03 true resid norm 1.249999999504e+03 ||r(i)||/||b|| 1.000000000000e+00
1 KSP unpreconditioned resid norm 1.249999999504e+03 true resid norm 1.249999999504e+03 ||r(i)||/||b|| 1.000000000000e+00
2 KSP unpreconditioned resid norm 2.177384421333e+02 true resid norm 7.820675831572e+08 ||r(i)||/||b|| 6.256540667740e+05
3 KSP unpreconditioned resid norm 1.081777973894e+02 true resid norm 4.201341309844e+01 ||r(i)||/||b|| 3.361073049209e-02
4 KSP unpreconditioned resid norm 2.640821579390e+01 true resid norm 2.805535044414e+01 ||r(i)||/||b|| 2.244428036422e-02
5 KSP unpreconditioned resid norm 1.905945882684e+01 true resid norm 8.069791145039e+08 ||r(i)||/||b|| 6.455832918593e+05
6 KSP unpreconditioned resid norm 1.567136385119e+01 true resid norm 8.070387542176e+08 ||r(i)||/||b|| 6.456310036303e+05
7 KSP unpreconditioned resid norm 1.361991788497e+01 true resid norm 8.070692097320e+08 ||r(i)||/||b|| 6.456553680418e+05
8 KSP unpreconditioned resid norm 1.220806886567e+01 true resid norm 8.070876916759e+08 ||r(i)||/||b|| 6.456701535969e+05
9 KSP unpreconditioned resid norm 1.116032600035e+01 true resid norm 8.071001011474e+08 ||r(i)||/||b|| 6.456800811741e+05
10 KSP unpreconditioned resid norm 1.034303066560e+01 true resid norm 8.071090086418e+08 ||r(i)||/||b|| 6.456872071697e+05
Linear solve converged due to CONVERGED_ITS iterations 10
Line search: gnorm after quadratic fit 1.089937861665e+07
Line search: Cubic step no good, shrinking lambda, current gnorm 1.035457409406e+06 lambda=1.0000000000000002e-02
Line search: Cubic step no good, shrinking lambda, current gnorm 1.030236386438e+05 lambda=1.0000000000000002e-03
Line search: Cubic step no good, shrinking lambda, current gnorm 1.037190239506e+04 lambda=1.0000000000000003e-04
Line search: Cubic step no good, shrinking lambda, current gnorm 1.619425850339e+03 lambda=1.0000000000000004e-05
Line search: Cubic step no good, shrinking lambda, current gnorm 1.254229242392e+03 lambda=1.0000000000000004e-06
Line search: Cubic step no good, shrinking lambda, current gnorm 1.250042364531e+03 lambda=1.0000000000000005e-07
Line search: Cubic step no good, shrinking lambda, current gnorm 1.250000412395e+03 lambda=1.0000000000000005e-08
Line search: Cubic step no good, shrinking lambda, current gnorm 1.250000002870e+03 lambda=1.0000000000000005e-09
Line search: Cubically determined step, current gnorm 1.249999999336e+03 lambda=1.3446088162829657e-10
|residual|_2 of individual variables:
potential: 1250
potentialliq: 1.43966e-23
em: 1.50979e-09
emliq: 3.08188e-10
Arp: 1.75738e-11
11 Nonlinear |R| = 1.250000e+03
0 KSP unpreconditioned resid norm 1.249999999336e+03 true resid norm 1.249999999336e+03 ||r(i)||/||b|| 1.000000000000e+00
1 KSP unpreconditioned resid norm 1.249999999336e+03 true resid norm 1.249999999336e+03 ||r(i)||/||b|| 1.000000000000e+00
2 KSP unpreconditioned resid norm 2.175597245573e+02 true resid norm 7.821062023649e+08 ||r(i)||/||b|| 6.256849622243e+05
3 KSP unpreconditioned resid norm 1.021043830525e+02 true resid norm 4.377529678057e+01 ||r(i)||/||b|| 3.502023744306e-02
4 KSP unpreconditioned resid norm 5.286981533874e+00 true resid norm 5.395069109689e+00 ||r(i)||/||b|| 4.316055290044e-03
5 KSP unpreconditioned resid norm 3.764995339701e+00 true resid norm 5.379330980748e+00 ||r(i)||/||b|| 4.303464786885e-03
6 KSP unpreconditioned resid norm 3.081422795569e+00 true resid norm 5.373914648069e+00 ||r(i)||/||b|| 4.299131720739e-03
7 KSP unpreconditioned resid norm 2.671775781134e+00 true resid norm 5.370824562979e+00 ||r(i)||/||b|| 4.296659652666e-03
8 KSP unpreconditioned resid norm 2.391423255607e+00 true resid norm 5.369627011698e+00 ||r(i)||/||b|| 4.295701611640e-03
9 KSP unpreconditioned resid norm 2.184105961898e+00 true resid norm 5.368237080665e+00 ||r(i)||/||b|| 4.294589666813e-03
10 KSP unpreconditioned resid norm 2.022781236997e+00 true resid norm 5.367390500679e+00 ||r(i)||/||b|| 4.293912402824e-03
Linear solve converged due to CONVERGED_ITS iterations 10
Line search: gnorm after quadratic fit 2.861209708657e+07
Line search: Cubic step no good, shrinking lambda, current gnorm 2.717339432175e+06 lambda=1.0000000000000002e-02
Line search: Cubic step no good, shrinking lambda, current gnorm 2.703384896028e+05 lambda=1.0000000000000002e-03
Line search: Cubic step no good, shrinking lambda, current gnorm 2.704852379134e+04 lambda=1.0000000000000003e-04
Line search: Cubic step no good, shrinking lambda, current gnorm 2.976958286838e+03 lambda=1.0000000000000004e-05
Line search: Cubic step no good, shrinking lambda, current gnorm 1.278865140534e+03 lambda=1.0000000000000004e-06
Line search: Cubic step no good, shrinking lambda, current gnorm 1.250291693482e+03 lambda=1.0000000000000005e-07
Line search: Cubic step no good, shrinking lambda, current gnorm 1.250002905364e+03 lambda=1.0000000000000005e-08
Line search: Cubic step no good, shrinking lambda, current gnorm 1.250000030612e+03 lambda=1.0000000000000005e-09
Line search: Cubic step no good, shrinking lambda, current gnorm 1.249999999762e+03 lambda=1.0000000000000006e-10
Line search: Cubically determined step, current gnorm 1.249999999322e+03 lambda=1.0930423985251739e-11
|residual|_2 of individual variables:
potential: 1250
potentialliq: 1.57218e-23
em: 1.45111e-09
emliq: 3.08188e-10
Arp: 1.70541e-11
12 Nonlinear |R| = 1.250000e+03
0 KSP unpreconditioned resid norm 1.249999999322e+03 true resid norm 1.249999999322e+03 ||r(i)||/||b|| 1.000000000000e+00
1 KSP unpreconditioned resid norm 1.249999999322e+03 true resid norm 1.249999999322e+03 ||r(i)||/||b|| 1.000000000000e+00
2 KSP unpreconditioned resid norm 2.175450305482e+02 true resid norm 7.821093769274e+08 ||r(i)||/||b|| 6.256875018812e+05
3 KSP unpreconditioned resid norm 1.856598278624e+01 true resid norm 8.053047975903e+08 ||r(i)||/||b|| 6.442438384215e+05
4 KSP unpreconditioned resid norm 1.364130029651e+01 true resid norm 8.053862298593e+08 ||r(i)||/||b|| 6.443089842368e+05
5 KSP unpreconditioned resid norm 1.128861070753e+01 true resid norm 2.009050004965e+01 ||r(i)||/||b|| 1.607240004843e-02
6 KSP unpreconditioned resid norm 9.844049423257e+00 true resid norm 2.009560936027e+01 ||r(i)||/||b|| 1.607648749694e-02
7 KSP unpreconditioned resid norm 8.841796624314e+00 true resid norm 8.054761688339e+08 ||r(i)||/||b|| 6.443809354165e+05
8 KSP unpreconditioned resid norm 8.094181518398e+00 true resid norm 8.054867356846e+08 ||r(i)||/||b|| 6.443893888970e+05
9 KSP unpreconditioned resid norm 7.508918360359e+00 true resid norm 2.010262861965e+01 ||r(i)||/||b|| 1.608210290444e-02
10 KSP unpreconditioned resid norm 7.034644000207e+00 true resid norm 2.010415273871e+01 ||r(i)||/||b|| 1.608332219969e-02
Linear solve converged due to CONVERGED_ITS iterations 10
Line search: gnorm after quadratic fit 2.456914540663e+07
Line search: Cubic step no good, shrinking lambda, current gnorm 1.193761683055e+07 lambda=5.0000000000000003e-02
Line search: Cubic step no good, shrinking lambda, current gnorm 1.163374916793e+06 lambda=5.0000000000000010e-03
Line search: Cubic step no good, shrinking lambda, current gnorm 1.160445421994e+05 lambda=5.0000000000000012e-04
Line search: Cubic step no good, shrinking lambda, current gnorm 1.166795782185e+04 lambda=5.0000000000000016e-05
Line search: Cubic step no good, shrinking lambda, current gnorm 1.705330567194e+03 lambda=5.0000000000000021e-06
Line search: Cubic step no good, shrinking lambda, current gnorm 1.255370779993e+03 lambda=5.0000000000000019e-07
Line search: Cubic step no good, shrinking lambda, current gnorm 1.250053955588e+03 lambda=5.0000000000000024e-08
Line search: Cubic step no good, shrinking lambda, current gnorm 1.250000542389e+03 lambda=5.0000000000000026e-09
Line search: Cubic step no good, shrinking lambda, current gnorm 1.250000003313e+03 lambda=5.0000000000000024e-10
Line search: Cubically determined step, current gnorm 1.249999999260e+03 lambda=5.0000000000000028e-11
|residual|_2 of individual variables:
potential: 1250
potentialliq: 1.7191e-23
em: 2.45625e-09
emliq: 3.08188e-10
Arp: 2.66652e-11
13 Nonlinear |R| = 1.250000e+03
0 KSP unpreconditioned resid norm 1.249999999260e+03 true resid norm 1.249999999260e+03 ||r(i)||/||b|| 1.000000000000e+00
1 KSP unpreconditioned resid norm 1.249999999260e+03 true resid norm 1.249999999260e+03 ||r(i)||/||b|| 1.000000000000e+00
2 KSP unpreconditioned resid norm 2.174776207655e+02 true resid norm 7.821239394476e+08 ||r(i)||/||b|| 6.256991519286e+05
3 KSP unpreconditioned resid norm 1.943727378418e+01 true resid norm 8.066314023718e+08 ||r(i)||/||b|| 6.453051222796e+05
4 KSP unpreconditioned resid norm 1.433914491640e+01 true resid norm 8.067188574536e+08 ||r(i)||/||b|| 6.453750863451e+05
5 KSP unpreconditioned resid norm 1.188388041985e+01 true resid norm 8.067515271785e+08 ||r(i)||/||b|| 6.454012221250e+05
6 KSP unpreconditioned resid norm 1.037058464541e+01 true resid norm 8.067686157693e+08 ||r(i)||/||b|| 6.454148929976e+05
7 KSP unpreconditioned resid norm 9.318830051429e+00 true resid norm 8.067791235915e+08 ||r(i)||/||b|| 6.454232992554e+05
8 KSP unpreconditioned resid norm 8.533423866217e+00 true resid norm 8.067862382342e+08 ||r(i)||/||b|| 6.454289909695e+05
9 KSP unpreconditioned resid norm 7.918104223305e+00 true resid norm 8.067913747994e+08 ||r(i)||/||b|| 6.454331002218e+05
10 KSP unpreconditioned resid norm 7.419191319918e+00 true resid norm 8.067952575497e+08 ||r(i)||/||b|| 6.454362064219e+05
Linear solve converged due to CONVERGED_ITS iterations 10
Line search: gnorm after quadratic fit 1.160995970104e+07
Line search: Cubic step no good, shrinking lambda, current gnorm 5.641751396461e+06 lambda=5.0000000000000003e-02
Line search: Cubic step no good, shrinking lambda, current gnorm 5.498787894103e+05 lambda=5.0000000000000010e-03
Line search: Cubic step no good, shrinking lambda, current gnorm 5.486101720170e+04 lambda=5.0000000000000012e-04
Line search: Cubic step no good, shrinking lambda, current gnorm 5.623924275130e+03 lambda=5.0000000000000016e-05
Line search: Cubic step no good, shrinking lambda, current gnorm 1.364972618139e+03 lambda=5.0000000000000021e-06
Line search: Cubic step no good, shrinking lambda, current gnorm 1.251201759974e+03 lambda=5.0000000000000019e-07
Line search: Cubic step no good, shrinking lambda, current gnorm 1.250012071103e+03 lambda=5.0000000000000024e-08
Line search: Cubic step no good, shrinking lambda, current gnorm 1.250000122223e+03 lambda=5.0000000000000026e-09
Line search: Cubic step no good, shrinking lambda, current gnorm 1.250000000724e+03 lambda=5.0000000000000024e-10
Line search: Cubically determined step, current gnorm 1.249999999169e+03 lambda=7.2394763111969298e-11
|residual|_2 of individual variables:
potential: 1250
potentialliq: 2.18503e-23
em: 3.99192e-09
emliq: 3.08188e-10
Arp: 4.21813e-11
14 Nonlinear |R| = 1.250000e+03
0 KSP unpreconditioned resid norm 1.249999999169e+03 true resid norm 1.249999999169e+03 ||r(i)||/||b|| 1.000000000000e+00
1 KSP unpreconditioned resid norm 1.249999999169e+03 true resid norm 1.249999999169e+03 ||r(i)||/||b|| 1.000000000000e+00
2 KSP unpreconditioned resid norm 2.173795908399e+02 true resid norm 7.821379143780e+08 ||r(i)||/||b|| 6.257103319182e+05
3 KSP unpreconditioned resid norm 1.943273247043e+01 true resid norm 2.106828788472e+01 ||r(i)||/||b|| 1.685463031898e-02
4 KSP unpreconditioned resid norm 1.433566777424e+01 true resid norm 2.114336582255e+01 ||r(i)||/||b|| 1.691469266928e-02
5 KSP unpreconditioned resid norm 1.188095918186e+01 true resid norm 2.117103567694e+01 ||r(i)||/||b|| 1.693682855281e-02
6 KSP unpreconditioned resid norm 1.036801734617e+01 true resid norm 2.118538176387e+01 ||r(i)||/||b|| 1.694830542236e-02
7 KSP unpreconditioned resid norm 9.316513162844e+00 true resid norm 2.119418433075e+01 ||r(i)||/||b|| 1.695534747587e-02
8 KSP unpreconditioned resid norm 8.531296064394e+00 true resid norm 2.120033082506e+01 ||r(i)||/||b|| 1.696026467132e-02
9 KSP unpreconditioned resid norm 7.916125715076e+00 true resid norm 2.120472884008e+01 ||r(i)||/||b|| 1.696378308334e-02
10 KSP unpreconditioned resid norm 7.417334541413e+00 true resid norm 2.120771141093e+01 ||r(i)||/||b|| 1.696616914002e-02
Linear solve converged due to CONVERGED_ITS iterations 10
Line search: gnorm after quadratic fit 1.163854034154e+07
Line search: Cubic step no good, shrinking lambda, current gnorm 5.655577109047e+06 lambda=5.0000000000000003e-02
Line search: Cubic step no good, shrinking lambda, current gnorm 5.512208074702e+05 lambda=5.0000000000000010e-03
Line search: Cubic step no good, shrinking lambda, current gnorm 5.499476969807e+04 lambda=5.0000000000000012e-04
Line search: Cubic step no good, shrinking lambda, current gnorm 5.636988389459e+03 lambda=5.0000000000000016e-05
Line search: Cubic step no good, shrinking lambda, current gnorm 1.365501184305e+03 lambda=5.0000000000000021e-06
Line search: Cubic step no good, shrinking lambda, current gnorm 1.251207931333e+03 lambda=5.0000000000000019e-07
Line search: Cubic step no good, shrinking lambda, current gnorm 1.250011916602e+03 lambda=5.0000000000000024e-08
Line search: Cubic step no good, shrinking lambda, current gnorm 1.250000106642e+03 lambda=5.0000000000000026e-09
Line search: Cubically determined step, current gnorm 1.249999999095e+03 lambda=5.0000000000000024e-10
|residual|_2 of individual variables:
potential: 1250
potentialliq: 8.77002e-23
em: 1.45953e-08
emliq: 0.0371094
Arp: 1.51539e-10
15 Nonlinear |R| = 1.250000e+03
Nonlinear solve did not converge due to DIVERGED_MAX_IT iterations 15
Solve Did NOT Converge!
-------------- next part --------------
Time Step 1, time = 6.4e-11
dt = 6.4e-11
|residual|_2 of individual variables:
potential: 1250
potentialliq: 0
em: 5.73385e-12
emliq: 3.08187e-10
Arp: 5.73385e-12
OHm: 5.14298e-23
H3Op: 5.14298e-23
0 Nonlinear |R| = 1.250000e+03
0 KSP unpreconditioned resid norm 1.250000000000e+03 true resid norm 1.250000000000e+03 ||r(i)||/||b|| 1.000000000000e+00
1 KSP unpreconditioned resid norm 1.247314056942e+03 true resid norm 1.247314056942e+03 ||r(i)||/||b|| 9.978512455538e-01
2 KSP unpreconditioned resid norm 4.860927105197e-02 true resid norm 4.776005315904e-02 ||r(i)||/||b|| 3.820804252723e-05
3 KSP unpreconditioned resid norm 4.787844580540e-02 true resid norm 4.752835483387e-02 ||r(i)||/||b|| 3.802268386709e-05
4 KSP unpreconditioned resid norm 4.437366678888e-02 true resid norm 2.756372625290e-01 ||r(i)||/||b|| 2.205098100232e-04
5 KSP unpreconditioned resid norm 1.696908986557e-05 true resid norm 5.105761919450e+00 ||r(i)||/||b|| 4.084609535560e-03
Linear solve converged due to CONVERGED_RTOL iterations 5
Line search: gnorm after quadratic fit 1.041879543556e+03
Line search: Quadratically determined step, lambda=1.6649745633364457e-01
|residual|_2 of individual variables:
potential: 1041.88
potentialliq: 5.8088e-17
em: 0.072466
emliq: 0.000102405
Arp: 0.000565487
OHm: 0.56393
H3Op: 1.01493
1 Nonlinear |R| = 1.041880e+03
0 KSP unpreconditioned resid norm 1.041879543556e+03 true resid norm 1.041879543556e+03 ||r(i)||/||b|| 1.000000000000e+00
1 KSP unpreconditioned resid norm 1.039102369962e+03 true resid norm 1.039102369962e+03 ||r(i)||/||b|| 9.973344580846e-01
2 KSP unpreconditioned resid norm 5.950166555328e-01 true resid norm 5.675973956612e-01 ||r(i)||/||b|| 5.447821671629e-04
3 KSP unpreconditioned resid norm 3.889464332642e-02 true resid norm 7.643703143268e-01 ||r(i)||/||b|| 7.336455726137e-04
4 KSP unpreconditioned resid norm 3.068800383788e-02 true resid norm 7.487202250394e-01 ||r(i)||/||b|| 7.186245566198e-04
5 KSP unpreconditioned resid norm 3.016106994371e-02 true resid norm 6.513904475182e-01 ||r(i)||/||b|| 6.252070611684e-04
6 KSP unpreconditioned resid norm 1.208157846839e-02 true resid norm 9.354083556895e+00 ||r(i)||/||b|| 8.978085436793e-03
7 KSP unpreconditioned resid norm 1.419355549112e-04 true resid norm 1.067868324404e+01 ||r(i)||/||b|| 1.024944131986e-02
Linear solve converged due to CONVERGED_RTOL iterations 7
Line search: Using full step: fnorm 1.041879543556e+03 gnorm 1.603618773063e+01
|residual|_2 of individual variables:
potential: 0.00484951
potentialliq: 1.5564e-05
em: 0.0536872
emliq: 0.125383
Arp: 0.00243548
OHm: 5.68961
H3Op: 14.9923
2 Nonlinear |R| = 1.603619e+01
0 KSP unpreconditioned resid norm 1.603618773063e+01 true resid norm 1.603618773063e+01 ||r(i)||/||b|| 1.000000000000e+00
1 KSP unpreconditioned resid norm 2.105432058619e+00 true resid norm 2.105432058619e+00 ||r(i)||/||b|| 1.312925549380e-01
2 KSP unpreconditioned resid norm 1.098264610573e-04 true resid norm 9.432163694076e-04 ||r(i)||/||b|| 5.881799248371e-05
Linear solve converged due to CONVERGED_RTOL iterations 2
Line search: gnorm after quadratic fit 1.443471881195e+01
Line search: Quadratically determined step, lambda=1.0000000000000001e-01
|residual|_2 of individual variables:
potential: 0.00436455
potentialliq: 1.40075e-05
em: 0.0483446
emliq: 0.112845
Arp: 0.00219194
OHm: 5.12895
H3Op: 13.4922
3 Nonlinear |R| = 1.443472e+01
0 KSP unpreconditioned resid norm 1.443471881195e+01 true resid norm 1.443471881195e+01 ||r(i)||/||b|| 1.000000000000e+00
1 KSP unpreconditioned resid norm 1.884409964965e+00 true resid norm 1.884409964965e+00 ||r(i)||/||b|| 1.305470504493e-01
2 KSP unpreconditioned resid norm 2.375784959854e-04 true resid norm 3.914568823513e-04 ||r(i)||/||b|| 2.711912074291e-05
3 KSP unpreconditioned resid norm 5.044976302331e-05 true resid norm 2.951075274163e-04 ||r(i)||/||b|| 2.044428653310e-05
Linear solve converged due to CONVERGED_RTOL iterations 3
Line search: Using full step: fnorm 1.443471881195e+01 gnorm 1.229473661441e+00
|residual|_2 of individual variables:
potential: 2.08202e-07
potentialliq: 6.92938e-10
em: 0.00232657
emliq: 2.5874e-05
Arp: 4.9962e-07
OHm: 0.90422
H3Op: 0.833058
4 Nonlinear |R| = 1.229474e+00
0 KSP unpreconditioned resid norm 1.229473661441e+00 true resid norm 1.229473661441e+00 ||r(i)||/||b|| 1.000000000000e+00
1 KSP unpreconditioned resid norm 2.888383837807e-03 true resid norm 2.888383837807e-03 ||r(i)||/||b|| 2.349284843095e-03
2 KSP unpreconditioned resid norm 2.027495745286e-03 true resid norm 2.038450555261e-03 ||r(i)||/||b|| 1.657986355618e-03
3 KSP unpreconditioned resid norm 1.984035502182e-04 true resid norm 2.605402233762e-04 ||r(i)||/||b|| 2.119120006775e-04
4 KSP unpreconditioned resid norm 3.857309277933e-05 true resid norm 1.325600705450e-02 ||r(i)||/||b|| 1.078185525257e-02
5 KSP unpreconditioned resid norm 5.087499897324e-07 true resid norm 1.454419650915e-02 ||r(i)||/||b|| 1.182961210580e-02
Linear solve converged due to CONVERGED_RTOL iterations 5
Line search: Using full step: fnorm 1.229473661441e+00 gnorm 2.408927331884e-01
|residual|_2 of individual variables:
potential: 1.00981e-07
potentialliq: 2.64367e-08
em: 1.33448e-05
emliq: 5.26969e-05
Arp: 9.8925e-07
OHm: 0.210391
H3Op: 0.117323
5 Nonlinear |R| = 2.408927e-01
0 KSP unpreconditioned resid norm 2.408927331884e-01 true resid norm 2.408927331884e-01 ||r(i)||/||b|| 1.000000000000e+00
1 KSP unpreconditioned resid norm 6.838178821298e-04 true resid norm 6.838178821298e-04 ||r(i)||/||b|| 2.838682068483e-03
2 KSP unpreconditioned resid norm 3.862120943656e-06 true resid norm 6.724384237429e-05 ||r(i)||/||b|| 2.791443373333e-04
3 KSP unpreconditioned resid norm 1.405230248854e-07 true resid norm 7.134981350742e-05 ||r(i)||/||b|| 2.961891484357e-04
Linear solve converged due to CONVERGED_RTOL iterations 3
Line search: Using full step: fnorm 2.408927331884e-01 gnorm 3.538570473573e-02
|residual|_2 of individual variables:
potential: 1.65546e-08
potentialliq: 8.63114e-11
em: 6.27176e-08
emliq: 3.26416e-07
Arp: 6.13426e-09
OHm: 0.0333149
H3Op: 0.0119275
6 Nonlinear |R| = 3.538570e-02
0 KSP unpreconditioned resid norm 3.538570473573e-02 true resid norm 3.538570473573e-02 ||r(i)||/||b|| 1.000000000000e+00
1 KSP unpreconditioned resid norm 3.691161652629e-05 true resid norm 3.691161652629e-05 ||r(i)||/||b|| 1.043122266519e-03
2 KSP unpreconditioned resid norm 5.495132451034e-06 true resid norm 2.671458554487e-05 ||r(i)||/||b|| 7.549541755458e-04
3 KSP unpreconditioned resid norm 7.269366746948e-08 true resid norm 2.624083112879e-05 ||r(i)||/||b|| 7.415658759593e-04
Linear solve converged due to CONVERGED_RTOL iterations 3
Line search: Using full step: fnorm 3.538570473573e-02 gnorm 1.647449494568e-03
|residual|_2 of individual variables:
potential: 1.43289e-08
potentialliq: 2.5783e-11
em: 2.34555e-08
emliq: 8.30641e-08
Arp: 1.43537e-09
OHm: 0.00163337
H3Op: 0.000214957
7 Nonlinear |R| = 1.647449e-03
0 KSP unpreconditioned resid norm 1.647449494568e-03 true resid norm 1.647449494568e-03 ||r(i)||/||b|| 1.000000000000e+00
1 KSP unpreconditioned resid norm 1.677765659074e-06 true resid norm 1.677765659074e-06 ||r(i)||/||b|| 1.018401877937e-03
2 KSP unpreconditioned resid norm 7.379250907830e-07 true resid norm 1.186104631721e-06 ||r(i)||/||b|| 7.199641844145e-04
3 KSP unpreconditioned resid norm 1.512636252550e-07 true resid norm 1.191859851853e-06 ||r(i)||/||b|| 7.234575965958e-04
4 KSP unpreconditioned resid norm 1.499662963136e-07 true resid norm 1.507725064021e-06 ||r(i)||/||b|| 9.151874269849e-04
5 KSP unpreconditioned resid norm 8.585859763891e-08 true resid norm 5.572122911654e-05 ||r(i)||/||b|| 3.382272373161e-02
6 KSP unpreconditioned resid norm 9.699214415589e-11 true resid norm 5.632188557407e-05 ||r(i)||/||b|| 3.418732152930e-02
Linear solve converged due to CONVERGED_RTOL iterations 6
Line search: Using full step: fnorm 1.647449494568e-03 gnorm 5.673172147266e-05
|residual|_2 of individual variables:
potential: 5.69344e-09
potentialliq: 1.05299e-10
em: 1.18164e-07
emliq: 3.30266e-07
Arp: 6.25748e-09
OHm: 2.53197e-05
H3Op: 5.07669e-05
8 Nonlinear |R| = 5.673172e-05
0 KSP unpreconditioned resid norm 5.673172147266e-05 true resid norm 5.673172147266e-05 ||r(i)||/||b|| 1.000000000000e+00
1 KSP unpreconditioned resid norm 2.283987852633e-05 true resid norm 2.283987852633e-05 ||r(i)||/||b|| 4.025944909382e-01
2 KSP unpreconditioned resid norm 1.620979330079e-09 true resid norm 9.175446330767e-09 ||r(i)||/||b|| 1.617339663347e-04
3 KSP unpreconditioned resid norm 3.189928680388e-11 true resid norm 8.666557521570e-09 ||r(i)||/||b|| 1.527638734838e-04
Linear solve converged due to CONVERGED_RTOL iterations 3
Line search: Using full step: fnorm 5.673172147266e-05 gnorm 5.106381275877e-09
|residual|_2 of individual variables:
potential: 6.59384e-12
potentialliq: 8.77051e-15
em: 9.64886e-12
emliq: 3.45014e-11
Arp: 6.84834e-13
OHm: 2.14381e-09
H3Op: 4.63442e-09
9 Nonlinear |R| = 5.106381e-09
Nonlinear solve converged due to CONVERGED_FNORM_RELATIVE iterations 9
Solve Converged!
More information about the petsc-users
mailing list