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 <[email protected]> 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.
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!
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!
