This page was generated from unit-3.8-nonlmin/nonlmin.ipynb.
3.8 Nonlinear minimization problems¶
We consider problems of the form
\[\text{find } u \in V \text{ s.t. } E(u) \leq E(v) \quad \forall~ v \in V.\]
[1]:
from netgen import gui
from ngsolve import *
Scalar minimization problems¶
As a first example we take \(V = H^1_0\) and
\[E(u) = \int_{\Omega} \vert \nabla u \vert^2 + u^4 - fu ~dx.\]
The minimization is equivalent to solving the nonlinear PDE:
\[- \Delta u + 4 u^3 = f \text{ in } \Omega\]
We solve the PDE with a Newton iteration.
[2]:
from netgen.geom2d import unit_square
mesh = Mesh (unit_square.GenerateMesh(maxh=0.2))
V = H1(mesh, order=4, dirichlet=[1,2,3,4])
u = V.TrialFunction()
To solve the problem we use the Variation integrator. Based on the symbolic description of the energy functional, it is able to
evaluate the energy functional (
Energy)
\[E(u) \qquad (E:V \to \mathbb{R})\]
compute the Gateau derivative for a given \(u\) (
Apply):
\[A(u)(v) = E'(u)(v) \qquad (A(u): V \to \mathbb{R})\]
compute the second derivative (
AssembleLinearization)
\[(\delta A)(w)(u,v) \qquad (\delta A(w): V\times V \to \mathbb{R})\]
[3]:
a = BilinearForm (V, symmetric=True)
a += Variation ( (grad(u)*grad(u) + u**4-u) * dx)
Equivalent to:
a += (2 * grad(u) * grad(v) + 4*u*u*u*v - 1 * v)*dx
(which has the same form as the problems in the nonlinear example)
We recall the Newton iteration (cf. unit-3.7 ) we make the loop:
Given an initial guess \(u^0\)
loop over \(i=0,..\) until convergence:
Compute linearization: $A u^i + \delta `A(u^i) :nbsphinx-math:Delta `u^{i} = 0 $:
\(f^i = A u^i\)
\(B^i = \delta A(u^i)\)
Solve \(B^i \Delta u^i = -f^i\)
Update \(u^{i+1} = u^i + \Delta u^{i}\)
Evaluate stopping criteria
Evaluate \(E(u^{i+1})\)
As a stopping criteria we take \(\langle A u^i,\Delta u^i \rangle = \langle A u^i, A u^i \rangle_{(B^i)^{-1}}< \varepsilon\).
[4]:
def SolveNonlinearMinProblem(a,gfu,tol=1e-13,maxits=25):
res = gfu.vec.CreateVector()
du = gfu.vec.CreateVector()
for it in range(maxits):
print ("Newton iteration {:3}".format(it),end="")
print ("energy = {:16}".format(a.Energy(gfu.vec)),end="")
#solve linearized problem:
a.Apply (gfu.vec, res)
a.AssembleLinearization (gfu.vec)
inv = a.mat.Inverse(V.FreeDofs())
du.data = inv * res
#update iteration
gfu.vec.data -= du
#stopping criteria
stopcritval = sqrt(abs(InnerProduct(du,res)))
print ("<A u",it,", A u",it,">_{-1}^0.5 = ", stopcritval)
if stopcritval < tol:
break
Redraw(blocking=True)
[5]:
gfu = GridFunction (V)
gfu.vec[:] = 0
Draw(gfu,mesh,"u")
SolveNonlinearMinProblem(a,gfu)
print ("energy = ", a.Energy(gfu.vec))
Newton iteration 0energy = 0.0<A u 0 , A u 0 >_{-1}^0.5 = 0.1325595865311363
Newton iteration 1energy = -0.008785678644368577<A u 1 , A u 1 >_{-1}^0.5 = 1.110759852324267e-05
Newton iteration 2energy = -0.008785678706057953<A u 2 , A u 2 >_{-1}^0.5 = 2.807474658179167e-13
Newton iteration 3energy = -0.008785678706057953<A u 3 , A u 3 >_{-1}^0.5 = 3.112658076268718e-17
energy = -0.008785678706057952
Again, a Newton for minimization is shipped with NGSolve:
[6]:
from ngsolve.solvers import *
gfu.vec[:] = 0
NewtonMinimization(a,gfu)
Redraw()
Newton iteration 0
Energy: 0.0
err = 0.13255958653113625
Newton iteration 1
Energy: -0.008785678644368577
err = 1.1107598523186047e-05
Newton iteration 2
Energy: -0.008785678706057953
err = 2.807427909405764e-13
Nonlinear elasticity¶
We consider a beam which is fixed on one side and is subject to gravity only. We assume a Neo-Hookean hyperelastic material. The model is a nonlinear minimization problem with
\[E(v) := \int_{\Omega} \frac{\mu}{2} ( \operatorname{tr}(F^T F-I)+\frac{2 \mu}{\lambda} \operatorname{det}(F^T F)^{\frac{\lambda}{2\mu}} - 1) - \gamma ~ (f,v) ~~ dx\]
where \(\mu\) and \(\lambda\) are the Lamé parameters and \(F = I + D v\) where \(v: \Omega \to \mathbb{R}^2\) is the sought for displacement.
[7]:
import netgen.geom2d as geom2d
from ngsolve import *
geo = geom2d.SplineGeometry()
pnums = [ geo.AddPoint (x,y,maxh=0.01) for x,y in [(0,0), (1,0), (1,0.1), (0,0.1)] ]
for p1,p2,bc in [(0,1,"bot"), (1,2,"right"), (2,3,"top"), (3,0,"left")]:
geo.Append(["line", pnums[p1], pnums[p2]], bc=bc)
mesh = Mesh(geo.GenerateMesh(maxh=0.05))
[8]:
# E module and poisson number:
E, nu = 210, 0.2
# Lamé constants:
mu = E / 2 / (1+nu)
lam = E * nu / ((1+nu)*(1-2*nu))
V = H1(mesh, order=2, dirichlet="left", dim=mesh.dim)
u = V.TrialFunction()
#gravity:
force = CoefficientFunction( (0,-1) )
[9]:
def Pow(a, b):
return exp (log(a)*b)
def NeoHook (C):
return 0.5 * mu * (Trace(C-I) + 2*mu/lam * Pow(Det(C), -lam/2/mu) - 1)
I = Id(mesh.dim)
F = I + Grad(u)
C = F.trans * F.trans
factor = Parameter(1.0)
a = BilinearForm(V, symmetric=True)
a += Variation( NeoHook (C).Compile() * dx
-factor * (InnerProduct(force,u) ).Compile() * dx)
We want to solve the minimization problem for \(\gamma = 5\). Due to the high nonlinearity in the problem, the Newton iteration will not convergence with any initial guess. We approach the case \(\gamma = 5\) by solving problems with \(\gamma = i/10\) for \(i=1,..,50\) and taking the solution of the previous problem as an initial guess.
[10]:
gfu = GridFunction(V)
gfu.vec[:] = 0
Draw (gfu, mesh, "u")
SetVisualization (deformation=True)
res = gfu.vec.CreateVector()
du = gfu.vec.CreateVector()
for loadstep in range(50):
print ("loadstep", loadstep)
factor.Set ((loadstep+1)/10)
SolveNonlinearMinProblem(a,gfu)
Redraw()
loadstep 0
Newton iteration 0energy = 8.75<A u 0 , A u 0 >_{-1}^0.5 = 0.0011216716617545427
Newton iteration 1energy = 8.76961764499874<A u 1 , A u 1 >_{-1}^0.5 = 0.15680152051564286
Newton iteration 2energy = 8.76451981188013<A u 2 , A u 2 >_{-1}^0.5 = 0.14374491176191792
Newton iteration 3energy = 9.24297982875597<A u 3 , A u 3 >_{-1}^0.5 = 2.9073222521249513
Newton iteration 4energy = 4348.831600937326<A u 4 , A u 4 >_{-1}^0.5 = 93.39683046695056
Newton iteration 5energy = 28.320949192630028<A u 5 , A u 5 >_{-1}^0.5 = 4.847718610473548
Newton iteration 6energy = 15.060324244381892<A u 6 , A u 6 >_{-1}^0.5 = 3.9705920866694853
Newton iteration 7energy = -3.131979688116065<A u 7 , A u 7 >_{-1}^0.5 = 2.8581744301821517
Newton iteration 8energy = -26.686639262379913<A u 8 , A u 8 >_{-1}^0.5 = 7.212177209123498
Newton iteration 9energy = 61.315691428655555<A u 9 , A u 9 >_{-1}^0.5 = 9.56656547659594
Newton iteration 10energy = -89.7422394042369<A u 10 , A u 10 >_{-1}^0.5 = 14.731944195414332
Newton iteration 11energy = 5.635757625102997<A u 11 , A u 11 >_{-1}^0.5 = 8.57621922723777
Newton iteration 12energy = 1133.0<A u 12 , A u 12 >_{-1}^0.5 = 47.98674504144227
Newton iteration 13energy = -57.0<A u 13 , A u 13 >_{-1}^0.5 = 9.28044245744348
Newton iteration 14energy = 164.0<A u 14 , A u 14 >_{-1}^0.5 = 5.368684424616927
Newton iteration 15energy = -242.0<A u 15 , A u 15 >_{-1}^0.5 = 16.199237072492675
Newton iteration 16energy = 150.0<A u 16 , A u 16 >_{-1}^0.5 = 3.2336012757537365
Newton iteration 17energy = 370.0<A u 17 , A u 17 >_{-1}^0.5 = 11.00470846656969
Newton iteration 18energy = -1364.0<A u 18 , A u 18 >_{-1}^0.5 = 40.26432623263448
Newton iteration 19energy = -504.0<A u 19 , A u 19 >_{-1}^0.5 = 21.097533659806224
Newton iteration 20energy = nan<A u 20 , A u 20 >_{-1}^0.5 = nan
Newton iteration 21energy = nan<A u 21 , A u 21 >_{-1}^0.5 = nan
Newton iteration 22energy = nan<A u 22 , A u 22 >_{-1}^0.5 = nan
Newton iteration 23energy = nan<A u 23 , A u 23 >_{-1}^0.5 = nan
Newton iteration 24energy = nan<A u 24 , A u 24 >_{-1}^0.5 = nan
loadstep 1
Newton iteration 0energy = nan<A u 0 , A u 0 >_{-1}^0.5 = nan
Newton iteration 1energy = nan<A u 1 , A u 1 >_{-1}^0.5 = nan
Newton iteration 2energy = nan<A u 2 , A u 2 >_{-1}^0.5 = nan
Newton iteration 3energy = nan<A u 3 , A u 3 >_{-1}^0.5 = nan
Newton iteration 4energy = nan<A u 4 , A u 4 >_{-1}^0.5 = nan
Newton iteration 5energy = nan<A u 5 , A u 5 >_{-1}^0.5 = nan
Newton iteration 6energy = nan<A u 6 , A u 6 >_{-1}^0.5 = nan
Newton iteration 7energy = nan<A u 7 , A u 7 >_{-1}^0.5 = nan
Newton iteration 8energy = nan<A u 8 , A u 8 >_{-1}^0.5 = nan
Newton iteration 9energy = nan<A u 9 , A u 9 >_{-1}^0.5 = nan
Newton iteration 10energy = nan<A u 10 , A u 10 >_{-1}^0.5 = nan
Newton iteration 11energy = nan<A u 11 , A u 11 >_{-1}^0.5 = nan
Newton iteration 12energy = nan<A u 12 , A u 12 >_{-1}^0.5 = nan
Newton iteration 13energy = nan<A u 13 , A u 13 >_{-1}^0.5 = nan
Newton iteration 14energy = nan<A u 14 , A u 14 >_{-1}^0.5 = nan
Newton iteration 15energy = nan<A u 15 , A u 15 >_{-1}^0.5 = nan
Newton iteration 16energy = nan<A u 16 , A u 16 >_{-1}^0.5 = nan
Newton iteration 17energy = nan<A u 17 , A u 17 >_{-1}^0.5 = nan
Newton iteration 18energy = nan<A u 18 , A u 18 >_{-1}^0.5 = nan
Newton iteration 19energy = nan<A u 19 , A u 19 >_{-1}^0.5 = nan
Newton iteration 20energy = nan<A u 20 , A u 20 >_{-1}^0.5 = nan
Newton iteration 21energy = nan<A u 21 , A u 21 >_{-1}^0.5 = nan
Newton iteration 22energy = nan<A u 22 , A u 22 >_{-1}^0.5 = nan
Newton iteration 23energy = nan<A u 23 , A u 23 >_{-1}^0.5 = nan
Newton iteration 24energy = nan<A u 24 , A u 24 >_{-1}^0.5 = nan
loadstep 2
Newton iteration 0energy = nan<A u 0 , A u 0 >_{-1}^0.5 = nan
Newton iteration 1energy = nan<A u 1 , A u 1 >_{-1}^0.5 = nan
Newton iteration 2energy = nan<A u 2 , A u 2 >_{-1}^0.5 = nan
Newton iteration 3energy = nan<A u 3 , A u 3 >_{-1}^0.5 = nan
Newton iteration 4energy = nan<A u 4 , A u 4 >_{-1}^0.5 = nan
Newton iteration 5energy = nan<A u 5 , A u 5 >_{-1}^0.5 = nan
Newton iteration 6energy = nan<A u 6 , A u 6 >_{-1}^0.5 = nan
Newton iteration 7energy = nan<A u 7 , A u 7 >_{-1}^0.5 = nan
Newton iteration 8energy = nan<A u 8 , A u 8 >_{-1}^0.5 = nan
Newton iteration 9energy = nan<A u 9 , A u 9 >_{-1}^0.5 = nan
Newton iteration 10energy = nan<A u 10 , A u 10 >_{-1}^0.5 = nan
Newton iteration 11energy = nan<A u 11 , A u 11 >_{-1}^0.5 = nan
Newton iteration 12energy = nan<A u 12 , A u 12 >_{-1}^0.5 = nan
Newton iteration 13energy = nan<A u 13 , A u 13 >_{-1}^0.5 = nan
Newton iteration 14energy = nan<A u 14 , A u 14 >_{-1}^0.5 = nan
Newton iteration 15energy = nan<A u 15 , A u 15 >_{-1}^0.5 = nan
Newton iteration 16energy = nan<A u 16 , A u 16 >_{-1}^0.5 = nan
Newton iteration 17energy = nan<A u 17 , A u 17 >_{-1}^0.5 = nan
Newton iteration 18energy = nan<A u 18 , A u 18 >_{-1}^0.5 = nan
Newton iteration 19energy = nan<A u 19 , A u 19 >_{-1}^0.5 = nan
Newton iteration 20energy = nan<A u 20 , A u 20 >_{-1}^0.5 = nan
Newton iteration 21energy = nan<A u 21 , A u 21 >_{-1}^0.5 = nan
Newton iteration 22energy = nan<A u 22 , A u 22 >_{-1}^0.5 = nan
Newton iteration 23energy = nan<A u 23 , A u 23 >_{-1}^0.5 = nan
Newton iteration 24energy = nan<A u 24 , A u 24 >_{-1}^0.5 = nan
loadstep 3
Newton iteration 0energy = nan<A u 0 , A u 0 >_{-1}^0.5 = nan
Newton iteration 1energy = nan<A u 1 , A u 1 >_{-1}^0.5 = nan
Newton iteration 2energy = nan<A u 2 , A u 2 >_{-1}^0.5 = nan
Newton iteration 3energy = nan<A u 3 , A u 3 >_{-1}^0.5 = nan
Newton iteration 4energy = nan<A u 4 , A u 4 >_{-1}^0.5 = nan
Newton iteration 5energy = nan<A u 5 , A u 5 >_{-1}^0.5 = nan
Newton iteration 6energy = nan<A u 6 , A u 6 >_{-1}^0.5 = nan
Newton iteration 7energy = nan<A u 7 , A u 7 >_{-1}^0.5 = nan
Newton iteration 8energy = nan<A u 8 , A u 8 >_{-1}^0.5 = nan
Newton iteration 9energy = nan<A u 9 , A u 9 >_{-1}^0.5 = nan
Newton iteration 10energy = nan<A u 10 , A u 10 >_{-1}^0.5 = nan
Newton iteration 11energy = nan<A u 11 , A u 11 >_{-1}^0.5 = nan
Newton iteration 12energy = nan<A u 12 , A u 12 >_{-1}^0.5 = nan
Newton iteration 13energy = nan<A u 13 , A u 13 >_{-1}^0.5 = nan
Newton iteration 14energy = nan<A u 14 , A u 14 >_{-1}^0.5 = nan
Newton iteration 15energy = nan<A u 15 , A u 15 >_{-1}^0.5 = nan
Newton iteration 16energy = nan<A u 16 , A u 16 >_{-1}^0.5 = nan
Newton iteration 17energy = nan<A u 17 , A u 17 >_{-1}^0.5 = nan
Newton iteration 18energy = nan<A u 18 , A u 18 >_{-1}^0.5 = nan
Newton iteration 19energy = nan<A u 19 , A u 19 >_{-1}^0.5 = nan
Newton iteration 20energy = nan<A u 20 , A u 20 >_{-1}^0.5 = nan
Newton iteration 21energy = nan<A u 21 , A u 21 >_{-1}^0.5 = nan
Newton iteration 22energy = nan<A u 22 , A u 22 >_{-1}^0.5 = nan
Newton iteration 23energy = nan<A u 23 , A u 23 >_{-1}^0.5 = nan
Newton iteration 24energy = nan<A u 24 , A u 24 >_{-1}^0.5 = nan
loadstep 4
Newton iteration 0energy = nan<A u 0 , A u 0 >_{-1}^0.5 = nan
Newton iteration 1energy = nan<A u 1 , A u 1 >_{-1}^0.5 = nan
Newton iteration 2energy = nan<A u 2 , A u 2 >_{-1}^0.5 = nan
Newton iteration 3energy = nan<A u 3 , A u 3 >_{-1}^0.5 = nan
Newton iteration 4energy = nan<A u 4 , A u 4 >_{-1}^0.5 = nan
Newton iteration 5energy = nan<A u 5 , A u 5 >_{-1}^0.5 = nan
Newton iteration 6energy = nan<A u 6 , A u 6 >_{-1}^0.5 = nan
Newton iteration 7energy = nan<A u 7 , A u 7 >_{-1}^0.5 = nan
Newton iteration 8energy = nan<A u 8 , A u 8 >_{-1}^0.5 = nan
Newton iteration 9energy = nan<A u 9 , A u 9 >_{-1}^0.5 = nan
Newton iteration 10energy = nan<A u 10 , A u 10 >_{-1}^0.5 = nan
Newton iteration 11energy = nan<A u 11 , A u 11 >_{-1}^0.5 = nan
Newton iteration 12energy = nan<A u 12 , A u 12 >_{-1}^0.5 = nan
Newton iteration 13energy = nan<A u 13 , A u 13 >_{-1}^0.5 = nan
Newton iteration 14energy = nan<A u 14 , A u 14 >_{-1}^0.5 = nan
Newton iteration 15energy = nan<A u 15 , A u 15 >_{-1}^0.5 = nan
Newton iteration 16energy = nan<A u 16 , A u 16 >_{-1}^0.5 = nan
Newton iteration 17energy = nan<A u 17 , A u 17 >_{-1}^0.5 = nan
Newton iteration 18energy = nan<A u 18 , A u 18 >_{-1}^0.5 = nan
Newton iteration 19energy = nan<A u 19 , A u 19 >_{-1}^0.5 = nan
Newton iteration 20energy = nan<A u 20 , A u 20 >_{-1}^0.5 = nan
Newton iteration 21energy = nan<A u 21 , A u 21 >_{-1}^0.5 = nan
Newton iteration 22energy = nan<A u 22 , A u 22 >_{-1}^0.5 = nan
Newton iteration 23energy = nan<A u 23 , A u 23 >_{-1}^0.5 = nan
Newton iteration 24energy = nan<A u 24 , A u 24 >_{-1}^0.5 = nan
loadstep 5
Newton iteration 0energy = nan<A u 0 , A u 0 >_{-1}^0.5 = nan
Newton iteration 1energy = nan<A u 1 , A u 1 >_{-1}^0.5 = nan
Newton iteration 2energy = nan<A u 2 , A u 2 >_{-1}^0.5 = nan
Newton iteration 3energy = nan<A u 3 , A u 3 >_{-1}^0.5 = nan
Newton iteration 4energy = nan<A u 4 , A u 4 >_{-1}^0.5 = nan
Newton iteration 5energy = nan<A u 5 , A u 5 >_{-1}^0.5 = nan
Newton iteration 6energy = nan<A u 6 , A u 6 >_{-1}^0.5 = nan
Newton iteration 7energy = nan<A u 7 , A u 7 >_{-1}^0.5 = nan
Newton iteration 8energy = nan<A u 8 , A u 8 >_{-1}^0.5 = nan
Newton iteration 9energy = nan<A u 9 , A u 9 >_{-1}^0.5 = nan
Newton iteration 10energy = nan<A u 10 , A u 10 >_{-1}^0.5 = nan
Newton iteration 11energy = nan<A u 11 , A u 11 >_{-1}^0.5 = nan
Newton iteration 12energy = nan<A u 12 , A u 12 >_{-1}^0.5 = nan
Newton iteration 13energy = nan<A u 13 , A u 13 >_{-1}^0.5 = nan
Newton iteration 14energy = nan<A u 14 , A u 14 >_{-1}^0.5 = nan
Newton iteration 15energy = nan<A u 15 , A u 15 >_{-1}^0.5 = nan
Newton iteration 16energy = nan<A u 16 , A u 16 >_{-1}^0.5 = nan
Newton iteration 17energy = nan<A u 17 , A u 17 >_{-1}^0.5 = nan
Newton iteration 18energy = nan<A u 18 , A u 18 >_{-1}^0.5 = nan
Newton iteration 19energy = nan<A u 19 , A u 19 >_{-1}^0.5 = nan
Newton iteration 20energy = nan<A u 20 , A u 20 >_{-1}^0.5 = nan
Newton iteration 21energy = nan<A u 21 , A u 21 >_{-1}^0.5 = nan
Newton iteration 22energy = nan<A u 22 , A u 22 >_{-1}^0.5 = nan
Newton iteration 23energy = nan<A u 23 , A u 23 >_{-1}^0.5 = nan
Newton iteration 24energy = nan<A u 24 , A u 24 >_{-1}^0.5 = nan
loadstep 6
Newton iteration 0energy = nan<A u 0 , A u 0 >_{-1}^0.5 = nan
Newton iteration 1energy = nan<A u 1 , A u 1 >_{-1}^0.5 = nan
Newton iteration 2energy = nan<A u 2 , A u 2 >_{-1}^0.5 = nan
Newton iteration 3energy = nan<A u 3 , A u 3 >_{-1}^0.5 = nan
Newton iteration 4energy = nan<A u 4 , A u 4 >_{-1}^0.5 = nan
Newton iteration 5energy = nan<A u 5 , A u 5 >_{-1}^0.5 = nan
Newton iteration 6energy = nan<A u 6 , A u 6 >_{-1}^0.5 = nan
Newton iteration 7energy = nan<A u 7 , A u 7 >_{-1}^0.5 = nan
Newton iteration 8energy = nan<A u 8 , A u 8 >_{-1}^0.5 = nan
Newton iteration 9energy = nan<A u 9 , A u 9 >_{-1}^0.5 = nan
Newton iteration 10energy = nan<A u 10 , A u 10 >_{-1}^0.5 = nan
Newton iteration 11energy = nan<A u 11 , A u 11 >_{-1}^0.5 = nan
Newton iteration 12energy = nan<A u 12 , A u 12 >_{-1}^0.5 = nan
Newton iteration 13energy = nan<A u 13 , A u 13 >_{-1}^0.5 = nan
Newton iteration 14energy = nan<A u 14 , A u 14 >_{-1}^0.5 = nan
Newton iteration 15energy = nan<A u 15 , A u 15 >_{-1}^0.5 = nan
Newton iteration 16energy = nan<A u 16 , A u 16 >_{-1}^0.5 = nan
Newton iteration 17energy = nan<A u 17 , A u 17 >_{-1}^0.5 = nan
Newton iteration 18energy = nan<A u 18 , A u 18 >_{-1}^0.5 = nan
Newton iteration 19energy = nan<A u 19 , A u 19 >_{-1}^0.5 = nan
Newton iteration 20energy = nan<A u 20 , A u 20 >_{-1}^0.5 = nan
Newton iteration 21energy = nan<A u 21 , A u 21 >_{-1}^0.5 = nan
Newton iteration 22energy = nan<A u 22 , A u 22 >_{-1}^0.5 = nan
Newton iteration 23energy = nan<A u 23 , A u 23 >_{-1}^0.5 = nan
Newton iteration 24energy = nan<A u 24 , A u 24 >_{-1}^0.5 = nan
loadstep 7
Newton iteration 0energy = nan<A u 0 , A u 0 >_{-1}^0.5 = nan
Newton iteration 1energy = nan<A u 1 , A u 1 >_{-1}^0.5 = nan
Newton iteration 2energy = nan<A u 2 , A u 2 >_{-1}^0.5 = nan
Newton iteration 3energy = nan<A u 3 , A u 3 >_{-1}^0.5 = nan
Newton iteration 4energy = nan<A u 4 , A u 4 >_{-1}^0.5 = nan
Newton iteration 5energy = nan<A u 5 , A u 5 >_{-1}^0.5 = nan
Newton iteration 6energy = nan<A u 6 , A u 6 >_{-1}^0.5 = nan
Newton iteration 7energy = nan<A u 7 , A u 7 >_{-1}^0.5 = nan
Newton iteration 8energy = nan<A u 8 , A u 8 >_{-1}^0.5 = nan
Newton iteration 9energy = nan<A u 9 , A u 9 >_{-1}^0.5 = nan
Newton iteration 10energy = nan<A u 10 , A u 10 >_{-1}^0.5 = nan
Newton iteration 11energy = nan<A u 11 , A u 11 >_{-1}^0.5 = nan
Newton iteration 12energy = nan<A u 12 , A u 12 >_{-1}^0.5 = nan
Newton iteration 13energy = nan<A u 13 , A u 13 >_{-1}^0.5 = nan
Newton iteration 14energy = nan<A u 14 , A u 14 >_{-1}^0.5 = nan
Newton iteration 15energy = nan<A u 15 , A u 15 >_{-1}^0.5 = nan
Newton iteration 16energy = nan<A u 16 , A u 16 >_{-1}^0.5 = nan
Newton iteration 17energy = nan<A u 17 , A u 17 >_{-1}^0.5 = nan
Newton iteration 18energy = nan<A u 18 , A u 18 >_{-1}^0.5 = nan
Newton iteration 19energy = nan<A u 19 , A u 19 >_{-1}^0.5 = nan
Newton iteration 20energy = nan<A u 20 , A u 20 >_{-1}^0.5 = nan
Newton iteration 21energy = nan<A u 21 , A u 21 >_{-1}^0.5 = nan
Newton iteration 22energy = nan<A u 22 , A u 22 >_{-1}^0.5 = nan
Newton iteration 23energy = nan<A u 23 , A u 23 >_{-1}^0.5 = nan
Newton iteration 24energy = nan<A u 24 , A u 24 >_{-1}^0.5 = nan
loadstep 8
Newton iteration 0energy = nan<A u 0 , A u 0 >_{-1}^0.5 = nan
Newton iteration 1energy = nan<A u 1 , A u 1 >_{-1}^0.5 = nan
Newton iteration 2energy = nan<A u 2 , A u 2 >_{-1}^0.5 = nan
Newton iteration 3energy = nan<A u 3 , A u 3 >_{-1}^0.5 = nan
Newton iteration 4energy = nan<A u 4 , A u 4 >_{-1}^0.5 = nan
Newton iteration 5energy = nan<A u 5 , A u 5 >_{-1}^0.5 = nan
Newton iteration 6energy = nan<A u 6 , A u 6 >_{-1}^0.5 = nan
Newton iteration 7energy = nan<A u 7 , A u 7 >_{-1}^0.5 = nan
Newton iteration 8energy = nan<A u 8 , A u 8 >_{-1}^0.5 = nan
Newton iteration 9energy = nan<A u 9 , A u 9 >_{-1}^0.5 = nan
Newton iteration 10energy = nan<A u 10 , A u 10 >_{-1}^0.5 = nan
Newton iteration 11energy = nan<A u 11 , A u 11 >_{-1}^0.5 = nan
Newton iteration 12energy = nan<A u 12 , A u 12 >_{-1}^0.5 = nan
Newton iteration 13energy = nan<A u 13 , A u 13 >_{-1}^0.5 = nan
Newton iteration 14energy = nan<A u 14 , A u 14 >_{-1}^0.5 = nan
Newton iteration 15energy = nan<A u 15 , A u 15 >_{-1}^0.5 = nan
Newton iteration 16energy = nan<A u 16 , A u 16 >_{-1}^0.5 = nan
Newton iteration 17energy = nan<A u 17 , A u 17 >_{-1}^0.5 = nan
Newton iteration 18energy = nan<A u 18 , A u 18 >_{-1}^0.5 = nan
Newton iteration 19energy = nan<A u 19 , A u 19 >_{-1}^0.5 = nan
Newton iteration 20energy = nan<A u 20 , A u 20 >_{-1}^0.5 = nan
Newton iteration 21energy = nan<A u 21 , A u 21 >_{-1}^0.5 = nan
Newton iteration 22energy = nan<A u 22 , A u 22 >_{-1}^0.5 = nan
Newton iteration 23energy = nan<A u 23 , A u 23 >_{-1}^0.5 = nan
Newton iteration 24energy = nan<A u 24 , A u 24 >_{-1}^0.5 = nan
loadstep 9
Newton iteration 0energy = nan<A u 0 , A u 0 >_{-1}^0.5 = nan
Newton iteration 1energy = nan<A u 1 , A u 1 >_{-1}^0.5 = nan
Newton iteration 2energy = nan<A u 2 , A u 2 >_{-1}^0.5 = nan
Newton iteration 3energy = nan<A u 3 , A u 3 >_{-1}^0.5 = nan
Newton iteration 4energy = nan<A u 4 , A u 4 >_{-1}^0.5 = nan
Newton iteration 5energy = nan<A u 5 , A u 5 >_{-1}^0.5 = nan
Newton iteration 6energy = nan<A u 6 , A u 6 >_{-1}^0.5 = nan
Newton iteration 7energy = nan<A u 7 , A u 7 >_{-1}^0.5 = nan
Newton iteration 8energy = nan<A u 8 , A u 8 >_{-1}^0.5 = nan
Newton iteration 9energy = nan<A u 9 , A u 9 >_{-1}^0.5 = nan
Newton iteration 10energy = nan<A u 10 , A u 10 >_{-1}^0.5 = nan
Newton iteration 11energy = nan<A u 11 , A u 11 >_{-1}^0.5 = nan
Newton iteration 12energy = nan<A u 12 , A u 12 >_{-1}^0.5 = nan
Newton iteration 13energy = nan<A u 13 , A u 13 >_{-1}^0.5 = nan
Newton iteration 14energy = nan<A u 14 , A u 14 >_{-1}^0.5 = nan
Newton iteration 15energy = nan<A u 15 , A u 15 >_{-1}^0.5 = nan
Newton iteration 16energy = nan<A u 16 , A u 16 >_{-1}^0.5 = nan
Newton iteration 17energy = nan<A u 17 , A u 17 >_{-1}^0.5 = nan
Newton iteration 18energy = nan<A u 18 , A u 18 >_{-1}^0.5 = nan
Newton iteration 19energy = nan<A u 19 , A u 19 >_{-1}^0.5 = nan
Newton iteration 20energy = nan<A u 20 , A u 20 >_{-1}^0.5 = nan
Newton iteration 21energy = nan<A u 21 , A u 21 >_{-1}^0.5 = nan
Newton iteration 22energy = nan<A u 22 , A u 22 >_{-1}^0.5 = nan
Newton iteration 23energy = nan<A u 23 , A u 23 >_{-1}^0.5 = nan
Newton iteration 24energy = nan<A u 24 , A u 24 >_{-1}^0.5 = nan
loadstep 10
Newton iteration 0energy = nan<A u 0 , A u 0 >_{-1}^0.5 = nan
Newton iteration 1energy = nan<A u 1 , A u 1 >_{-1}^0.5 = nan
Newton iteration 2energy = nan<A u 2 , A u 2 >_{-1}^0.5 = nan
Newton iteration 3energy = nan<A u 3 , A u 3 >_{-1}^0.5 = nan
Newton iteration 4energy = nan<A u 4 , A u 4 >_{-1}^0.5 = nan
Newton iteration 5energy = nan<A u 5 , A u 5 >_{-1}^0.5 = nan
Newton iteration 6energy = nan<A u 6 , A u 6 >_{-1}^0.5 = nan
Newton iteration 7energy = nan<A u 7 , A u 7 >_{-1}^0.5 = nan
Newton iteration 8energy = nan<A u 8 , A u 8 >_{-1}^0.5 = nan
Newton iteration 9energy = nan<A u 9 , A u 9 >_{-1}^0.5 = nan
Newton iteration 10energy = nan<A u 10 , A u 10 >_{-1}^0.5 = nan
Newton iteration 11energy = nan<A u 11 , A u 11 >_{-1}^0.5 = nan
Newton iteration 12energy = nan<A u 12 , A u 12 >_{-1}^0.5 = nan
Newton iteration 13energy = nan<A u 13 , A u 13 >_{-1}^0.5 = nan
Newton iteration 14energy = nan<A u 14 , A u 14 >_{-1}^0.5 = nan
Newton iteration 15energy = nan<A u 15 , A u 15 >_{-1}^0.5 = nan
Newton iteration 16energy = nan<A u 16 , A u 16 >_{-1}^0.5 = nan
Newton iteration 17energy = nan<A u 17 , A u 17 >_{-1}^0.5 = nan
Newton iteration 18energy = nan<A u 18 , A u 18 >_{-1}^0.5 = nan
Newton iteration 19energy = nan<A u 19 , A u 19 >_{-1}^0.5 = nan
Newton iteration 20energy = nan<A u 20 , A u 20 >_{-1}^0.5 = nan
Newton iteration 21energy = nan<A u 21 , A u 21 >_{-1}^0.5 = nan
Newton iteration 22energy = nan<A u 22 , A u 22 >_{-1}^0.5 = nan
Newton iteration 23energy = nan<A u 23 , A u 23 >_{-1}^0.5 = nan
Newton iteration 24energy = nan<A u 24 , A u 24 >_{-1}^0.5 = nan
loadstep 11
Newton iteration 0energy = nan<A u 0 , A u 0 >_{-1}^0.5 = nan
Newton iteration 1energy = nan<A u 1 , A u 1 >_{-1}^0.5 = nan
Newton iteration 2energy = nan<A u 2 , A u 2 >_{-1}^0.5 = nan
Newton iteration 3energy = nan<A u 3 , A u 3 >_{-1}^0.5 = nan
Newton iteration 4energy = nan<A u 4 , A u 4 >_{-1}^0.5 = nan
Newton iteration 5energy = nan<A u 5 , A u 5 >_{-1}^0.5 = nan
Newton iteration 6energy = nan<A u 6 , A u 6 >_{-1}^0.5 = nan
Newton iteration 7energy = nan<A u 7 , A u 7 >_{-1}^0.5 = nan
Newton iteration 8energy = nan<A u 8 , A u 8 >_{-1}^0.5 = nan
Newton iteration 9energy = nan<A u 9 , A u 9 >_{-1}^0.5 = nan
Newton iteration 10energy = nan<A u 10 , A u 10 >_{-1}^0.5 = nan
Newton iteration 11energy = nan<A u 11 , A u 11 >_{-1}^0.5 = nan
Newton iteration 12energy = nan<A u 12 , A u 12 >_{-1}^0.5 = nan
Newton iteration 13energy = nan<A u 13 , A u 13 >_{-1}^0.5 = nan
Newton iteration 14energy = nan<A u 14 , A u 14 >_{-1}^0.5 = nan
Newton iteration 15energy = nan<A u 15 , A u 15 >_{-1}^0.5 = nan
Newton iteration 16energy = nan<A u 16 , A u 16 >_{-1}^0.5 = nan
Newton iteration 17energy = nan<A u 17 , A u 17 >_{-1}^0.5 = nan
Newton iteration 18energy = nan<A u 18 , A u 18 >_{-1}^0.5 = nan
Newton iteration 19energy = nan<A u 19 , A u 19 >_{-1}^0.5 = nan
Newton iteration 20energy = nan<A u 20 , A u 20 >_{-1}^0.5 = nan
Newton iteration 21energy = nan<A u 21 , A u 21 >_{-1}^0.5 = nan
Newton iteration 22energy = nan<A u 22 , A u 22 >_{-1}^0.5 = nan
Newton iteration 23energy = nan<A u 23 , A u 23 >_{-1}^0.5 = nan
Newton iteration 24energy = nan<A u 24 , A u 24 >_{-1}^0.5 = nan
loadstep 12
Newton iteration 0energy = nan<A u 0 , A u 0 >_{-1}^0.5 = nan
Newton iteration 1energy = nan<A u 1 , A u 1 >_{-1}^0.5 = nan
Newton iteration 2energy = nan<A u 2 , A u 2 >_{-1}^0.5 = nan
Newton iteration 3energy = nan<A u 3 , A u 3 >_{-1}^0.5 = nan
Newton iteration 4energy = nan<A u 4 , A u 4 >_{-1}^0.5 = nan
Newton iteration 5energy = nan<A u 5 , A u 5 >_{-1}^0.5 = nan
Newton iteration 6energy = nan<A u 6 , A u 6 >_{-1}^0.5 = nan
Newton iteration 7energy = nan<A u 7 , A u 7 >_{-1}^0.5 = nan
Newton iteration 8energy = nan<A u 8 , A u 8 >_{-1}^0.5 = nan
Newton iteration 9energy = nan<A u 9 , A u 9 >_{-1}^0.5 = nan
Newton iteration 10energy = nan<A u 10 , A u 10 >_{-1}^0.5 = nan
Newton iteration 11energy = nan<A u 11 , A u 11 >_{-1}^0.5 = nan
Newton iteration 12energy = nan<A u 12 , A u 12 >_{-1}^0.5 = nan
Newton iteration 13energy = nan<A u 13 , A u 13 >_{-1}^0.5 = nan
Newton iteration 14energy = nan<A u 14 , A u 14 >_{-1}^0.5 = nan
Newton iteration 15energy = nan<A u 15 , A u 15 >_{-1}^0.5 = nan
Newton iteration 16energy = nan<A u 16 , A u 16 >_{-1}^0.5 = nan
Newton iteration 17energy = nan<A u 17 , A u 17 >_{-1}^0.5 = nan
Newton iteration 18energy = nan<A u 18 , A u 18 >_{-1}^0.5 = nan
Newton iteration 19energy = nan<A u 19 , A u 19 >_{-1}^0.5 = nan
Newton iteration 20energy = nan<A u 20 , A u 20 >_{-1}^0.5 = nan
Newton iteration 21energy = nan<A u 21 , A u 21 >_{-1}^0.5 = nan
Newton iteration 22energy = nan<A u 22 , A u 22 >_{-1}^0.5 = nan
Newton iteration 23energy = nan<A u 23 , A u 23 >_{-1}^0.5 = nan
Newton iteration 24energy = nan<A u 24 , A u 24 >_{-1}^0.5 = nan
loadstep 13
Newton iteration 0energy = nan<A u 0 , A u 0 >_{-1}^0.5 = nan
Newton iteration 1energy = nan<A u 1 , A u 1 >_{-1}^0.5 = nan
Newton iteration 2energy = nan<A u 2 , A u 2 >_{-1}^0.5 = nan
Newton iteration 3energy = nan<A u 3 , A u 3 >_{-1}^0.5 = nan
Newton iteration 4energy = nan<A u 4 , A u 4 >_{-1}^0.5 = nan
Newton iteration 5energy = nan<A u 5 , A u 5 >_{-1}^0.5 = nan
Newton iteration 6energy = nan<A u 6 , A u 6 >_{-1}^0.5 = nan
Newton iteration 7energy = nan<A u 7 , A u 7 >_{-1}^0.5 = nan
Newton iteration 8energy = nan<A u 8 , A u 8 >_{-1}^0.5 = nan
Newton iteration 9energy = nan<A u 9 , A u 9 >_{-1}^0.5 = nan
Newton iteration 10energy = nan<A u 10 , A u 10 >_{-1}^0.5 = nan
Newton iteration 11energy = nan<A u 11 , A u 11 >_{-1}^0.5 = nan
Newton iteration 12energy = nan<A u 12 , A u 12 >_{-1}^0.5 = nan
Newton iteration 13energy = nan<A u 13 , A u 13 >_{-1}^0.5 = nan
Newton iteration 14energy = nan<A u 14 , A u 14 >_{-1}^0.5 = nan
Newton iteration 15energy = nan<A u 15 , A u 15 >_{-1}^0.5 = nan
Newton iteration 16energy = nan<A u 16 , A u 16 >_{-1}^0.5 = nan
Newton iteration 17energy = nan<A u 17 , A u 17 >_{-1}^0.5 = nan
Newton iteration 18energy = nan<A u 18 , A u 18 >_{-1}^0.5 = nan
Newton iteration 19energy = nan<A u 19 , A u 19 >_{-1}^0.5 = nan
Newton iteration 20energy = nan<A u 20 , A u 20 >_{-1}^0.5 = nan
Newton iteration 21energy = nan<A u 21 , A u 21 >_{-1}^0.5 = nan
Newton iteration 22energy = nan<A u 22 , A u 22 >_{-1}^0.5 = nan
Newton iteration 23energy = nan<A u 23 , A u 23 >_{-1}^0.5 = nan
Newton iteration 24energy = nan<A u 24 , A u 24 >_{-1}^0.5 = nan
loadstep 14
Newton iteration 0energy = nan<A u 0 , A u 0 >_{-1}^0.5 = nan
Newton iteration 1energy = nan<A u 1 , A u 1 >_{-1}^0.5 = nan
Newton iteration 2energy = nan<A u 2 , A u 2 >_{-1}^0.5 = nan
Newton iteration 3energy = nan<A u 3 , A u 3 >_{-1}^0.5 = nan
Newton iteration 4energy = nan<A u 4 , A u 4 >_{-1}^0.5 = nan
Newton iteration 5energy = nan<A u 5 , A u 5 >_{-1}^0.5 = nan
Newton iteration 6energy = nan<A u 6 , A u 6 >_{-1}^0.5 = nan
Newton iteration 7energy = nan<A u 7 , A u 7 >_{-1}^0.5 = nan
Newton iteration 8energy = nan<A u 8 , A u 8 >_{-1}^0.5 = nan
Newton iteration 9energy = nan<A u 9 , A u 9 >_{-1}^0.5 = nan
Newton iteration 10energy = nan<A u 10 , A u 10 >_{-1}^0.5 = nan
Newton iteration 11energy = nan<A u 11 , A u 11 >_{-1}^0.5 = nan
Newton iteration 12energy = nan<A u 12 , A u 12 >_{-1}^0.5 = nan
Newton iteration 13energy = nan<A u 13 , A u 13 >_{-1}^0.5 = nan
Newton iteration 14energy = nan<A u 14 , A u 14 >_{-1}^0.5 = nan
Newton iteration 15energy = nan<A u 15 , A u 15 >_{-1}^0.5 = nan
Newton iteration 16energy = nan<A u 16 , A u 16 >_{-1}^0.5 = nan
Newton iteration 17energy = nan<A u 17 , A u 17 >_{-1}^0.5 = nan
Newton iteration 18energy = nan<A u 18 , A u 18 >_{-1}^0.5 = nan
Newton iteration 19energy = nan<A u 19 , A u 19 >_{-1}^0.5 = nan
Newton iteration 20energy = nan<A u 20 , A u 20 >_{-1}^0.5 = nan
Newton iteration 21energy = nan<A u 21 , A u 21 >_{-1}^0.5 = nan
Newton iteration 22energy = nan<A u 22 , A u 22 >_{-1}^0.5 = nan
Newton iteration 23energy = nan<A u 23 , A u 23 >_{-1}^0.5 = nan
Newton iteration 24energy = nan<A u 24 , A u 24 >_{-1}^0.5 = nan
loadstep 15
Newton iteration 0energy = nan<A u 0 , A u 0 >_{-1}^0.5 = nan
Newton iteration 1energy = nan<A u 1 , A u 1 >_{-1}^0.5 = nan
Newton iteration 2energy = nan<A u 2 , A u 2 >_{-1}^0.5 = nan
Newton iteration 3energy = nan<A u 3 , A u 3 >_{-1}^0.5 = nan
Newton iteration 4energy = nan<A u 4 , A u 4 >_{-1}^0.5 = nan
Newton iteration 5energy = nan<A u 5 , A u 5 >_{-1}^0.5 = nan
Newton iteration 6energy = nan<A u 6 , A u 6 >_{-1}^0.5 = nan
Newton iteration 7energy = nan<A u 7 , A u 7 >_{-1}^0.5 = nan
Newton iteration 8energy = nan<A u 8 , A u 8 >_{-1}^0.5 = nan
Newton iteration 9energy = nan<A u 9 , A u 9 >_{-1}^0.5 = nan
Newton iteration 10energy = nan<A u 10 , A u 10 >_{-1}^0.5 = nan
Newton iteration 11energy = nan<A u 11 , A u 11 >_{-1}^0.5 = nan
Newton iteration 12energy = nan<A u 12 , A u 12 >_{-1}^0.5 = nan
Newton iteration 13energy = nan<A u 13 , A u 13 >_{-1}^0.5 = nan
Newton iteration 14energy = nan<A u 14 , A u 14 >_{-1}^0.5 = nan
Newton iteration 15energy = nan<A u 15 , A u 15 >_{-1}^0.5 = nan
Newton iteration 16energy = nan<A u 16 , A u 16 >_{-1}^0.5 = nan
Newton iteration 17energy = nan<A u 17 , A u 17 >_{-1}^0.5 = nan
Newton iteration 18energy = nan<A u 18 , A u 18 >_{-1}^0.5 = nan
Newton iteration 19energy = nan<A u 19 , A u 19 >_{-1}^0.5 = nan
Newton iteration 20energy = nan<A u 20 , A u 20 >_{-1}^0.5 = nan
Newton iteration 21energy = nan<A u 21 , A u 21 >_{-1}^0.5 = nan
Newton iteration 22energy = nan<A u 22 , A u 22 >_{-1}^0.5 = nan
Newton iteration 23energy = nan<A u 23 , A u 23 >_{-1}^0.5 = nan
Newton iteration 24energy = nan<A u 24 , A u 24 >_{-1}^0.5 = nan
loadstep 16
Newton iteration 0energy = nan<A u 0 , A u 0 >_{-1}^0.5 = nan
Newton iteration 1energy = nan<A u 1 , A u 1 >_{-1}^0.5 = nan
Newton iteration 2energy = nan<A u 2 , A u 2 >_{-1}^0.5 = nan
Newton iteration 3energy = nan<A u 3 , A u 3 >_{-1}^0.5 = nan
Newton iteration 4energy = nan<A u 4 , A u 4 >_{-1}^0.5 = nan
Newton iteration 5energy = nan<A u 5 , A u 5 >_{-1}^0.5 = nan
Newton iteration 6energy = nan<A u 6 , A u 6 >_{-1}^0.5 = nan
Newton iteration 7energy = nan<A u 7 , A u 7 >_{-1}^0.5 = nan
Newton iteration 8energy = nan<A u 8 , A u 8 >_{-1}^0.5 = nan
Newton iteration 9energy = nan<A u 9 , A u 9 >_{-1}^0.5 = nan
Newton iteration 10energy = nan<A u 10 , A u 10 >_{-1}^0.5 = nan
Newton iteration 11energy = nan<A u 11 , A u 11 >_{-1}^0.5 = nan
Newton iteration 12energy = nan<A u 12 , A u 12 >_{-1}^0.5 = nan
Newton iteration 13energy = nan<A u 13 , A u 13 >_{-1}^0.5 = nan
Newton iteration 14energy = nan<A u 14 , A u 14 >_{-1}^0.5 = nan
Newton iteration 15energy = nan<A u 15 , A u 15 >_{-1}^0.5 = nan
Newton iteration 16energy = nan<A u 16 , A u 16 >_{-1}^0.5 = nan
Newton iteration 17energy = nan<A u 17 , A u 17 >_{-1}^0.5 = nan
Newton iteration 18energy = nan<A u 18 , A u 18 >_{-1}^0.5 = nan
Newton iteration 19energy = nan<A u 19 , A u 19 >_{-1}^0.5 = nan
Newton iteration 20energy = nan<A u 20 , A u 20 >_{-1}^0.5 = nan
Newton iteration 21energy = nan<A u 21 , A u 21 >_{-1}^0.5 = nan
Newton iteration 22energy = nan<A u 22 , A u 22 >_{-1}^0.5 = nan
Newton iteration 23energy = nan<A u 23 , A u 23 >_{-1}^0.5 = nan
Newton iteration 24energy = nan<A u 24 , A u 24 >_{-1}^0.5 = nan
loadstep 17
Newton iteration 0energy = nan<A u 0 , A u 0 >_{-1}^0.5 = nan
Newton iteration 1energy = nan<A u 1 , A u 1 >_{-1}^0.5 = nan
Newton iteration 2energy = nan<A u 2 , A u 2 >_{-1}^0.5 = nan
Newton iteration 3energy = nan<A u 3 , A u 3 >_{-1}^0.5 = nan
Newton iteration 4energy = nan<A u 4 , A u 4 >_{-1}^0.5 = nan
Newton iteration 5energy = nan<A u 5 , A u 5 >_{-1}^0.5 = nan
Newton iteration 6energy = nan<A u 6 , A u 6 >_{-1}^0.5 = nan
Newton iteration 7energy = nan<A u 7 , A u 7 >_{-1}^0.5 = nan
Newton iteration 8energy = nan<A u 8 , A u 8 >_{-1}^0.5 = nan
Newton iteration 9energy = nan<A u 9 , A u 9 >_{-1}^0.5 = nan
Newton iteration 10energy = nan<A u 10 , A u 10 >_{-1}^0.5 = nan
Newton iteration 11energy = nan<A u 11 , A u 11 >_{-1}^0.5 = nan
Newton iteration 12energy = nan<A u 12 , A u 12 >_{-1}^0.5 = nan
Newton iteration 13energy = nan<A u 13 , A u 13 >_{-1}^0.5 = nan
Newton iteration 14energy = nan<A u 14 , A u 14 >_{-1}^0.5 = nan
Newton iteration 15energy = nan<A u 15 , A u 15 >_{-1}^0.5 = nan
Newton iteration 16energy = nan<A u 16 , A u 16 >_{-1}^0.5 = nan
Newton iteration 17energy = nan<A u 17 , A u 17 >_{-1}^0.5 = nan
Newton iteration 18energy = nan<A u 18 , A u 18 >_{-1}^0.5 = nan
Newton iteration 19energy = nan<A u 19 , A u 19 >_{-1}^0.5 = nan
Newton iteration 20energy = nan<A u 20 , A u 20 >_{-1}^0.5 = nan
Newton iteration 21energy = nan<A u 21 , A u 21 >_{-1}^0.5 = nan
Newton iteration 22energy = nan<A u 22 , A u 22 >_{-1}^0.5 = nan
Newton iteration 23energy = nan<A u 23 , A u 23 >_{-1}^0.5 = nan
Newton iteration 24energy = nan<A u 24 , A u 24 >_{-1}^0.5 = nan
loadstep 18
Newton iteration 0energy = nan<A u 0 , A u 0 >_{-1}^0.5 = nan
Newton iteration 1energy = nan<A u 1 , A u 1 >_{-1}^0.5 = nan
Newton iteration 2energy = nan<A u 2 , A u 2 >_{-1}^0.5 = nan
Newton iteration 3energy = nan<A u 3 , A u 3 >_{-1}^0.5 = nan
Newton iteration 4energy = nan<A u 4 , A u 4 >_{-1}^0.5 = nan
Newton iteration 5energy = nan<A u 5 , A u 5 >_{-1}^0.5 = nan
Newton iteration 6energy = nan<A u 6 , A u 6 >_{-1}^0.5 = nan
Newton iteration 7energy = nan<A u 7 , A u 7 >_{-1}^0.5 = nan
Newton iteration 8energy = nan<A u 8 , A u 8 >_{-1}^0.5 = nan
Newton iteration 9energy = nan<A u 9 , A u 9 >_{-1}^0.5 = nan
Newton iteration 10energy = nan<A u 10 , A u 10 >_{-1}^0.5 = nan
Newton iteration 11energy = nan<A u 11 , A u 11 >_{-1}^0.5 = nan
Newton iteration 12energy = nan<A u 12 , A u 12 >_{-1}^0.5 = nan
Newton iteration 13energy = nan<A u 13 , A u 13 >_{-1}^0.5 = nan
Newton iteration 14energy = nan<A u 14 , A u 14 >_{-1}^0.5 = nan
Newton iteration 15energy = nan<A u 15 , A u 15 >_{-1}^0.5 = nan
Newton iteration 16energy = nan<A u 16 , A u 16 >_{-1}^0.5 = nan
Newton iteration 17energy = nan<A u 17 , A u 17 >_{-1}^0.5 = nan
Newton iteration 18energy = nan<A u 18 , A u 18 >_{-1}^0.5 = nan
Newton iteration 19energy = nan<A u 19 , A u 19 >_{-1}^0.5 = nan
Newton iteration 20energy = nan<A u 20 , A u 20 >_{-1}^0.5 = nan
Newton iteration 21energy = nan<A u 21 , A u 21 >_{-1}^0.5 = nan
Newton iteration 22energy = nan<A u 22 , A u 22 >_{-1}^0.5 = nan
Newton iteration 23energy = nan<A u 23 , A u 23 >_{-1}^0.5 = nan
Newton iteration 24energy = nan<A u 24 , A u 24 >_{-1}^0.5 = nan
loadstep 19
Newton iteration 0energy = nan<A u 0 , A u 0 >_{-1}^0.5 = nan
Newton iteration 1energy = nan<A u 1 , A u 1 >_{-1}^0.5 = nan
Newton iteration 2energy = nan<A u 2 , A u 2 >_{-1}^0.5 = nan
Newton iteration 3energy = nan<A u 3 , A u 3 >_{-1}^0.5 = nan
Newton iteration 4energy = nan<A u 4 , A u 4 >_{-1}^0.5 = nan
Newton iteration 5energy = nan<A u 5 , A u 5 >_{-1}^0.5 = nan
Newton iteration 6energy = nan<A u 6 , A u 6 >_{-1}^0.5 = nan
Newton iteration 7energy = nan<A u 7 , A u 7 >_{-1}^0.5 = nan
Newton iteration 8energy = nan<A u 8 , A u 8 >_{-1}^0.5 = nan
Newton iteration 9energy = nan<A u 9 , A u 9 >_{-1}^0.5 = nan
Newton iteration 10energy = nan<A u 10 , A u 10 >_{-1}^0.5 = nan
Newton iteration 11energy = nan<A u 11 , A u 11 >_{-1}^0.5 = nan
Newton iteration 12energy = nan<A u 12 , A u 12 >_{-1}^0.5 = nan
Newton iteration 13energy = nan<A u 13 , A u 13 >_{-1}^0.5 = nan
Newton iteration 14energy = nan<A u 14 , A u 14 >_{-1}^0.5 = nan
Newton iteration 15energy = nan<A u 15 , A u 15 >_{-1}^0.5 = nan
Newton iteration 16energy = nan<A u 16 , A u 16 >_{-1}^0.5 = nan
Newton iteration 17energy = nan<A u 17 , A u 17 >_{-1}^0.5 = nan
Newton iteration 18energy = nan<A u 18 , A u 18 >_{-1}^0.5 = nan
Newton iteration 19energy = nan<A u 19 , A u 19 >_{-1}^0.5 = nan
Newton iteration 20energy = nan<A u 20 , A u 20 >_{-1}^0.5 = nan
Newton iteration 21energy = nan<A u 21 , A u 21 >_{-1}^0.5 = nan
Newton iteration 22energy = nan<A u 22 , A u 22 >_{-1}^0.5 = nan
Newton iteration 23energy = nan<A u 23 , A u 23 >_{-1}^0.5 = nan
Newton iteration 24energy = nan<A u 24 , A u 24 >_{-1}^0.5 = nan
loadstep 20
Newton iteration 0energy = nan<A u 0 , A u 0 >_{-1}^0.5 = nan
Newton iteration 1energy = nan<A u 1 , A u 1 >_{-1}^0.5 = nan
Newton iteration 2energy = nan<A u 2 , A u 2 >_{-1}^0.5 = nan
Newton iteration 3energy = nan<A u 3 , A u 3 >_{-1}^0.5 = nan
Newton iteration 4energy = nan<A u 4 , A u 4 >_{-1}^0.5 = nan
Newton iteration 5energy = nan<A u 5 , A u 5 >_{-1}^0.5 = nan
Newton iteration 6energy = nan<A u 6 , A u 6 >_{-1}^0.5 = nan
Newton iteration 7energy = nan<A u 7 , A u 7 >_{-1}^0.5 = nan
Newton iteration 8energy = nan<A u 8 , A u 8 >_{-1}^0.5 = nan
Newton iteration 9energy = nan<A u 9 , A u 9 >_{-1}^0.5 = nan
Newton iteration 10energy = nan<A u 10 , A u 10 >_{-1}^0.5 = nan
Newton iteration 11energy = nan<A u 11 , A u 11 >_{-1}^0.5 = nan
Newton iteration 12energy = nan<A u 12 , A u 12 >_{-1}^0.5 = nan
Newton iteration 13energy = nan<A u 13 , A u 13 >_{-1}^0.5 = nan
Newton iteration 14energy = nan<A u 14 , A u 14 >_{-1}^0.5 = nan
Newton iteration 15energy = nan<A u 15 , A u 15 >_{-1}^0.5 = nan
Newton iteration 16energy = nan<A u 16 , A u 16 >_{-1}^0.5 = nan
Newton iteration 17energy = nan<A u 17 , A u 17 >_{-1}^0.5 = nan
Newton iteration 18energy = nan<A u 18 , A u 18 >_{-1}^0.5 = nan
Newton iteration 19energy = nan<A u 19 , A u 19 >_{-1}^0.5 = nan
Newton iteration 20energy = nan<A u 20 , A u 20 >_{-1}^0.5 = nan
Newton iteration 21energy = nan<A u 21 , A u 21 >_{-1}^0.5 = nan
Newton iteration 22energy = nan<A u 22 , A u 22 >_{-1}^0.5 = nan
Newton iteration 23energy = nan<A u 23 , A u 23 >_{-1}^0.5 = nan
Newton iteration 24energy = nan<A u 24 , A u 24 >_{-1}^0.5 = nan
loadstep 21
Newton iteration 0energy = nan<A u 0 , A u 0 >_{-1}^0.5 = nan
Newton iteration 1energy = nan<A u 1 , A u 1 >_{-1}^0.5 = nan
Newton iteration 2energy = nan<A u 2 , A u 2 >_{-1}^0.5 = nan
Newton iteration 3energy = nan<A u 3 , A u 3 >_{-1}^0.5 = nan
Newton iteration 4energy = nan<A u 4 , A u 4 >_{-1}^0.5 = nan
Newton iteration 5energy = nan<A u 5 , A u 5 >_{-1}^0.5 = nan
Newton iteration 6energy = nan<A u 6 , A u 6 >_{-1}^0.5 = nan
Newton iteration 7energy = nan<A u 7 , A u 7 >_{-1}^0.5 = nan
Newton iteration 8energy = nan<A u 8 , A u 8 >_{-1}^0.5 = nan
Newton iteration 9energy = nan<A u 9 , A u 9 >_{-1}^0.5 = nan
Newton iteration 10energy = nan<A u 10 , A u 10 >_{-1}^0.5 = nan
Newton iteration 11energy = nan<A u 11 , A u 11 >_{-1}^0.5 = nan
Newton iteration 12energy = nan<A u 12 , A u 12 >_{-1}^0.5 = nan
Newton iteration 13energy = nan<A u 13 , A u 13 >_{-1}^0.5 = nan
Newton iteration 14energy = nan<A u 14 , A u 14 >_{-1}^0.5 = nan
Newton iteration 15energy = nan<A u 15 , A u 15 >_{-1}^0.5 = nan
Newton iteration 16energy = nan<A u 16 , A u 16 >_{-1}^0.5 = nan
Newton iteration 17energy = nan<A u 17 , A u 17 >_{-1}^0.5 = nan
Newton iteration 18energy = nan<A u 18 , A u 18 >_{-1}^0.5 = nan
Newton iteration 19energy = nan<A u 19 , A u 19 >_{-1}^0.5 = nan
Newton iteration 20energy = nan<A u 20 , A u 20 >_{-1}^0.5 = nan
Newton iteration 21energy = nan<A u 21 , A u 21 >_{-1}^0.5 = nan
Newton iteration 22energy = nan<A u 22 , A u 22 >_{-1}^0.5 = nan
Newton iteration 23energy = nan<A u 23 , A u 23 >_{-1}^0.5 = nan
Newton iteration 24energy = nan<A u 24 , A u 24 >_{-1}^0.5 = nan
loadstep 22
Newton iteration 0energy = nan<A u 0 , A u 0 >_{-1}^0.5 = nan
Newton iteration 1energy = nan<A u 1 , A u 1 >_{-1}^0.5 = nan
Newton iteration 2energy = nan<A u 2 , A u 2 >_{-1}^0.5 = nan
Newton iteration 3energy = nan<A u 3 , A u 3 >_{-1}^0.5 = nan
Newton iteration 4energy = nan<A u 4 , A u 4 >_{-1}^0.5 = nan
Newton iteration 5energy = nan<A u 5 , A u 5 >_{-1}^0.5 = nan
Newton iteration 6energy = nan<A u 6 , A u 6 >_{-1}^0.5 = nan
Newton iteration 7energy = nan<A u 7 , A u 7 >_{-1}^0.5 = nan
Newton iteration 8energy = nan<A u 8 , A u 8 >_{-1}^0.5 = nan
Newton iteration 9energy = nan<A u 9 , A u 9 >_{-1}^0.5 = nan
Newton iteration 10energy = nan<A u 10 , A u 10 >_{-1}^0.5 = nan
Newton iteration 11energy = nan<A u 11 , A u 11 >_{-1}^0.5 = nan
Newton iteration 12energy = nan<A u 12 , A u 12 >_{-1}^0.5 = nan
Newton iteration 13energy = nan<A u 13 , A u 13 >_{-1}^0.5 = nan
Newton iteration 14energy = nan<A u 14 , A u 14 >_{-1}^0.5 = nan
Newton iteration 15energy = nan<A u 15 , A u 15 >_{-1}^0.5 = nan
Newton iteration 16energy = nan<A u 16 , A u 16 >_{-1}^0.5 = nan
Newton iteration 17energy = nan<A u 17 , A u 17 >_{-1}^0.5 = nan
Newton iteration 18energy = nan<A u 18 , A u 18 >_{-1}^0.5 = nan
Newton iteration 19energy = nan<A u 19 , A u 19 >_{-1}^0.5 = nan
Newton iteration 20energy = nan<A u 20 , A u 20 >_{-1}^0.5 = nan
Newton iteration 21energy = nan<A u 21 , A u 21 >_{-1}^0.5 = nan
Newton iteration 22energy = nan<A u 22 , A u 22 >_{-1}^0.5 = nan
Newton iteration 23energy = nan<A u 23 , A u 23 >_{-1}^0.5 = nan
Newton iteration 24energy = nan<A u 24 , A u 24 >_{-1}^0.5 = nan
loadstep 23
Newton iteration 0energy = nan<A u 0 , A u 0 >_{-1}^0.5 = nan
Newton iteration 1energy = nan<A u 1 , A u 1 >_{-1}^0.5 = nan
Newton iteration 2energy = nan<A u 2 , A u 2 >_{-1}^0.5 = nan
Newton iteration 3energy = nan<A u 3 , A u 3 >_{-1}^0.5 = nan
Newton iteration 4energy = nan<A u 4 , A u 4 >_{-1}^0.5 = nan
Newton iteration 5energy = nan<A u 5 , A u 5 >_{-1}^0.5 = nan
Newton iteration 6energy = nan<A u 6 , A u 6 >_{-1}^0.5 = nan
Newton iteration 7energy = nan<A u 7 , A u 7 >_{-1}^0.5 = nan
Newton iteration 8energy = nan<A u 8 , A u 8 >_{-1}^0.5 = nan
Newton iteration 9energy = nan<A u 9 , A u 9 >_{-1}^0.5 = nan
Newton iteration 10energy = nan<A u 10 , A u 10 >_{-1}^0.5 = nan
Newton iteration 11energy = nan<A u 11 , A u 11 >_{-1}^0.5 = nan
Newton iteration 12energy = nan<A u 12 , A u 12 >_{-1}^0.5 = nan
Newton iteration 13energy = nan<A u 13 , A u 13 >_{-1}^0.5 = nan
Newton iteration 14energy = nan<A u 14 , A u 14 >_{-1}^0.5 = nan
Newton iteration 15energy = nan<A u 15 , A u 15 >_{-1}^0.5 = nan
Newton iteration 16energy = nan<A u 16 , A u 16 >_{-1}^0.5 = nan
Newton iteration 17energy = nan<A u 17 , A u 17 >_{-1}^0.5 = nan
Newton iteration 18energy = nan<A u 18 , A u 18 >_{-1}^0.5 = nan
Newton iteration 19energy = nan<A u 19 , A u 19 >_{-1}^0.5 = nan
Newton iteration 20energy = nan<A u 20 , A u 20 >_{-1}^0.5 = nan
Newton iteration 21energy = nan<A u 21 , A u 21 >_{-1}^0.5 = nan
Newton iteration 22energy = nan<A u 22 , A u 22 >_{-1}^0.5 = nan
Newton iteration 23energy = nan<A u 23 , A u 23 >_{-1}^0.5 = nan
Newton iteration 24energy = nan<A u 24 , A u 24 >_{-1}^0.5 = nan
loadstep 24
Newton iteration 0energy = nan<A u 0 , A u 0 >_{-1}^0.5 = nan
Newton iteration 1energy = nan<A u 1 , A u 1 >_{-1}^0.5 = nan
Newton iteration 2energy = nan<A u 2 , A u 2 >_{-1}^0.5 = nan
Newton iteration 3energy = nan<A u 3 , A u 3 >_{-1}^0.5 = nan
Newton iteration 4energy = nan<A u 4 , A u 4 >_{-1}^0.5 = nan
Newton iteration 5energy = nan<A u 5 , A u 5 >_{-1}^0.5 = nan
Newton iteration 6energy = nan<A u 6 , A u 6 >_{-1}^0.5 = nan
Newton iteration 7energy = nan<A u 7 , A u 7 >_{-1}^0.5 = nan
Newton iteration 8energy = nan<A u 8 , A u 8 >_{-1}^0.5 = nan
Newton iteration 9energy = nan<A u 9 , A u 9 >_{-1}^0.5 = nan
Newton iteration 10energy = nan<A u 10 , A u 10 >_{-1}^0.5 = nan
Newton iteration 11energy = nan<A u 11 , A u 11 >_{-1}^0.5 = nan
Newton iteration 12energy = nan<A u 12 , A u 12 >_{-1}^0.5 = nan
Newton iteration 13energy = nan<A u 13 , A u 13 >_{-1}^0.5 = nan
Newton iteration 14energy = nan<A u 14 , A u 14 >_{-1}^0.5 = nan
Newton iteration 15energy = nan<A u 15 , A u 15 >_{-1}^0.5 = nan
Newton iteration 16energy = nan<A u 16 , A u 16 >_{-1}^0.5 = nan
Newton iteration 17energy = nan<A u 17 , A u 17 >_{-1}^0.5 = nan
Newton iteration 18energy = nan<A u 18 , A u 18 >_{-1}^0.5 = nan
Newton iteration 19energy = nan<A u 19 , A u 19 >_{-1}^0.5 = nan
Newton iteration 20energy = nan<A u 20 , A u 20 >_{-1}^0.5 = nan
Newton iteration 21energy = nan<A u 21 , A u 21 >_{-1}^0.5 = nan
Newton iteration 22energy = nan<A u 22 , A u 22 >_{-1}^0.5 = nan
Newton iteration 23energy = nan<A u 23 , A u 23 >_{-1}^0.5 = nan
Newton iteration 24energy = nan<A u 24 , A u 24 >_{-1}^0.5 = nan
loadstep 25
Newton iteration 0energy = nan<A u 0 , A u 0 >_{-1}^0.5 = nan
Newton iteration 1energy = nan<A u 1 , A u 1 >_{-1}^0.5 = nan
Newton iteration 2energy = nan<A u 2 , A u 2 >_{-1}^0.5 = nan
Newton iteration 3energy = nan<A u 3 , A u 3 >_{-1}^0.5 = nan
Newton iteration 4energy = nan<A u 4 , A u 4 >_{-1}^0.5 = nan
Newton iteration 5energy = nan<A u 5 , A u 5 >_{-1}^0.5 = nan
Newton iteration 6energy = nan<A u 6 , A u 6 >_{-1}^0.5 = nan
Newton iteration 7energy = nan<A u 7 , A u 7 >_{-1}^0.5 = nan
Newton iteration 8energy = nan<A u 8 , A u 8 >_{-1}^0.5 = nan
Newton iteration 9energy = nan<A u 9 , A u 9 >_{-1}^0.5 = nan
Newton iteration 10energy = nan<A u 10 , A u 10 >_{-1}^0.5 = nan
Newton iteration 11energy = nan<A u 11 , A u 11 >_{-1}^0.5 = nan
Newton iteration 12energy = nan<A u 12 , A u 12 >_{-1}^0.5 = nan
Newton iteration 13energy = nan<A u 13 , A u 13 >_{-1}^0.5 = nan
Newton iteration 14energy = nan<A u 14 , A u 14 >_{-1}^0.5 = nan
Newton iteration 15energy = nan<A u 15 , A u 15 >_{-1}^0.5 = nan
Newton iteration 16energy = nan<A u 16 , A u 16 >_{-1}^0.5 = nan
Newton iteration 17energy = nan<A u 17 , A u 17 >_{-1}^0.5 = nan
Newton iteration 18energy = nan<A u 18 , A u 18 >_{-1}^0.5 = nan
Newton iteration 19energy = nan<A u 19 , A u 19 >_{-1}^0.5 = nan
Newton iteration 20energy = nan<A u 20 , A u 20 >_{-1}^0.5 = nan
Newton iteration 21energy = nan<A u 21 , A u 21 >_{-1}^0.5 = nan
Newton iteration 22energy = nan<A u 22 , A u 22 >_{-1}^0.5 = nan
Newton iteration 23energy = nan<A u 23 , A u 23 >_{-1}^0.5 = nan
Newton iteration 24energy = nan<A u 24 , A u 24 >_{-1}^0.5 = nan
loadstep 26
Newton iteration 0energy = nan<A u 0 , A u 0 >_{-1}^0.5 = nan
Newton iteration 1energy = nan<A u 1 , A u 1 >_{-1}^0.5 = nan
Newton iteration 2energy = nan<A u 2 , A u 2 >_{-1}^0.5 = nan
Newton iteration 3energy = nan<A u 3 , A u 3 >_{-1}^0.5 = nan
Newton iteration 4energy = nan<A u 4 , A u 4 >_{-1}^0.5 = nan
Newton iteration 5energy = nan<A u 5 , A u 5 >_{-1}^0.5 = nan
Newton iteration 6energy = nan<A u 6 , A u 6 >_{-1}^0.5 = nan
Newton iteration 7energy = nan<A u 7 , A u 7 >_{-1}^0.5 = nan
Newton iteration 8energy = nan<A u 8 , A u 8 >_{-1}^0.5 = nan
Newton iteration 9energy = nan<A u 9 , A u 9 >_{-1}^0.5 = nan
Newton iteration 10energy = nan<A u 10 , A u 10 >_{-1}^0.5 = nan
Newton iteration 11energy = nan<A u 11 , A u 11 >_{-1}^0.5 = nan
Newton iteration 12energy = nan<A u 12 , A u 12 >_{-1}^0.5 = nan
Newton iteration 13energy = nan<A u 13 , A u 13 >_{-1}^0.5 = nan
Newton iteration 14energy = nan<A u 14 , A u 14 >_{-1}^0.5 = nan
Newton iteration 15energy = nan<A u 15 , A u 15 >_{-1}^0.5 = nan
Newton iteration 16energy = nan<A u 16 , A u 16 >_{-1}^0.5 = nan
Newton iteration 17energy = nan<A u 17 , A u 17 >_{-1}^0.5 = nan
Newton iteration 18energy = nan<A u 18 , A u 18 >_{-1}^0.5 = nan
Newton iteration 19energy = nan<A u 19 , A u 19 >_{-1}^0.5 = nan
Newton iteration 20energy = nan<A u 20 , A u 20 >_{-1}^0.5 = nan
Newton iteration 21energy = nan<A u 21 , A u 21 >_{-1}^0.5 = nan
Newton iteration 22energy = nan<A u 22 , A u 22 >_{-1}^0.5 = nan
Newton iteration 23energy = nan<A u 23 , A u 23 >_{-1}^0.5 = nan
Newton iteration 24energy = nan<A u 24 , A u 24 >_{-1}^0.5 = nan
loadstep 27
Newton iteration 0energy = nan<A u 0 , A u 0 >_{-1}^0.5 = nan
Newton iteration 1energy = nan<A u 1 , A u 1 >_{-1}^0.5 = nan
Newton iteration 2energy = nan<A u 2 , A u 2 >_{-1}^0.5 = nan
Newton iteration 3energy = nan<A u 3 , A u 3 >_{-1}^0.5 = nan
Newton iteration 4energy = nan<A u 4 , A u 4 >_{-1}^0.5 = nan
Newton iteration 5energy = nan<A u 5 , A u 5 >_{-1}^0.5 = nan
Newton iteration 6energy = nan<A u 6 , A u 6 >_{-1}^0.5 = nan
Newton iteration 7energy = nan<A u 7 , A u 7 >_{-1}^0.5 = nan
Newton iteration 8energy = nan<A u 8 , A u 8 >_{-1}^0.5 = nan
Newton iteration 9energy = nan<A u 9 , A u 9 >_{-1}^0.5 = nan
Newton iteration 10energy = nan<A u 10 , A u 10 >_{-1}^0.5 = nan
Newton iteration 11energy = nan<A u 11 , A u 11 >_{-1}^0.5 = nan
Newton iteration 12energy = nan<A u 12 , A u 12 >_{-1}^0.5 = nan
Newton iteration 13energy = nan<A u 13 , A u 13 >_{-1}^0.5 = nan
Newton iteration 14energy = nan<A u 14 , A u 14 >_{-1}^0.5 = nan
Newton iteration 15energy = nan<A u 15 , A u 15 >_{-1}^0.5 = nan
Newton iteration 16energy = nan<A u 16 , A u 16 >_{-1}^0.5 = nan
Newton iteration 17energy = nan<A u 17 , A u 17 >_{-1}^0.5 = nan
Newton iteration 18energy = nan<A u 18 , A u 18 >_{-1}^0.5 = nan
Newton iteration 19energy = nan<A u 19 , A u 19 >_{-1}^0.5 = nan
Newton iteration 20energy = nan<A u 20 , A u 20 >_{-1}^0.5 = nan
Newton iteration 21energy = nan<A u 21 , A u 21 >_{-1}^0.5 = nan
Newton iteration 22energy = nan<A u 22 , A u 22 >_{-1}^0.5 = nan
Newton iteration 23energy = nan<A u 23 , A u 23 >_{-1}^0.5 = nan
Newton iteration 24energy = nan<A u 24 , A u 24 >_{-1}^0.5 = nan
loadstep 28
Newton iteration 0energy = nan<A u 0 , A u 0 >_{-1}^0.5 = nan
Newton iteration 1energy = nan<A u 1 , A u 1 >_{-1}^0.5 = nan
Newton iteration 2energy = nan<A u 2 , A u 2 >_{-1}^0.5 = nan
Newton iteration 3energy = nan<A u 3 , A u 3 >_{-1}^0.5 = nan
Newton iteration 4energy = nan<A u 4 , A u 4 >_{-1}^0.5 = nan
Newton iteration 5energy = nan<A u 5 , A u 5 >_{-1}^0.5 = nan
Newton iteration 6energy = nan<A u 6 , A u 6 >_{-1}^0.5 = nan
Newton iteration 7energy = nan<A u 7 , A u 7 >_{-1}^0.5 = nan
Newton iteration 8energy = nan<A u 8 , A u 8 >_{-1}^0.5 = nan
Newton iteration 9energy = nan<A u 9 , A u 9 >_{-1}^0.5 = nan
Newton iteration 10energy = nan<A u 10 , A u 10 >_{-1}^0.5 = nan
Newton iteration 11energy = nan<A u 11 , A u 11 >_{-1}^0.5 = nan
Newton iteration 12energy = nan<A u 12 , A u 12 >_{-1}^0.5 = nan
Newton iteration 13energy = nan<A u 13 , A u 13 >_{-1}^0.5 = nan
Newton iteration 14energy = nan<A u 14 , A u 14 >_{-1}^0.5 = nan
Newton iteration 15energy = nan<A u 15 , A u 15 >_{-1}^0.5 = nan
Newton iteration 16energy = nan<A u 16 , A u 16 >_{-1}^0.5 = nan
Newton iteration 17energy = nan<A u 17 , A u 17 >_{-1}^0.5 = nan
Newton iteration 18energy = nan<A u 18 , A u 18 >_{-1}^0.5 = nan
Newton iteration 19energy = nan<A u 19 , A u 19 >_{-1}^0.5 = nan
Newton iteration 20energy = nan<A u 20 , A u 20 >_{-1}^0.5 = nan
Newton iteration 21energy = nan<A u 21 , A u 21 >_{-1}^0.5 = nan
Newton iteration 22energy = nan<A u 22 , A u 22 >_{-1}^0.5 = nan
Newton iteration 23energy = nan<A u 23 , A u 23 >_{-1}^0.5 = nan
Newton iteration 24energy = nan<A u 24 , A u 24 >_{-1}^0.5 = nan
loadstep 29
Newton iteration 0energy = nan<A u 0 , A u 0 >_{-1}^0.5 = nan
Newton iteration 1energy = nan<A u 1 , A u 1 >_{-1}^0.5 = nan
Newton iteration 2energy = nan<A u 2 , A u 2 >_{-1}^0.5 = nan
Newton iteration 3energy = nan<A u 3 , A u 3 >_{-1}^0.5 = nan
Newton iteration 4energy = nan<A u 4 , A u 4 >_{-1}^0.5 = nan
Newton iteration 5energy = nan<A u 5 , A u 5 >_{-1}^0.5 = nan
Newton iteration 6energy = nan<A u 6 , A u 6 >_{-1}^0.5 = nan
Newton iteration 7energy = nan<A u 7 , A u 7 >_{-1}^0.5 = nan
Newton iteration 8energy = nan<A u 8 , A u 8 >_{-1}^0.5 = nan
Newton iteration 9energy = nan<A u 9 , A u 9 >_{-1}^0.5 = nan
Newton iteration 10energy = nan<A u 10 , A u 10 >_{-1}^0.5 = nan
Newton iteration 11energy = nan<A u 11 , A u 11 >_{-1}^0.5 = nan
Newton iteration 12energy = nan<A u 12 , A u 12 >_{-1}^0.5 = nan
Newton iteration 13energy = nan<A u 13 , A u 13 >_{-1}^0.5 = nan
Newton iteration 14energy = nan<A u 14 , A u 14 >_{-1}^0.5 = nan
Newton iteration 15energy = nan<A u 15 , A u 15 >_{-1}^0.5 = nan
Newton iteration 16energy = nan<A u 16 , A u 16 >_{-1}^0.5 = nan
Newton iteration 17energy = nan<A u 17 , A u 17 >_{-1}^0.5 = nan
Newton iteration 18energy = nan<A u 18 , A u 18 >_{-1}^0.5 = nan
Newton iteration 19energy = nan<A u 19 , A u 19 >_{-1}^0.5 = nan
Newton iteration 20energy = nan<A u 20 , A u 20 >_{-1}^0.5 = nan
Newton iteration 21energy = nan<A u 21 , A u 21 >_{-1}^0.5 = nan
Newton iteration 22energy = nan<A u 22 , A u 22 >_{-1}^0.5 = nan
Newton iteration 23energy = nan<A u 23 , A u 23 >_{-1}^0.5 = nan
Newton iteration 24energy = nan<A u 24 , A u 24 >_{-1}^0.5 = nan
loadstep 30
Newton iteration 0energy = nan<A u 0 , A u 0 >_{-1}^0.5 = nan
Newton iteration 1energy = nan<A u 1 , A u 1 >_{-1}^0.5 = nan
Newton iteration 2energy = nan<A u 2 , A u 2 >_{-1}^0.5 = nan
Newton iteration 3energy = nan<A u 3 , A u 3 >_{-1}^0.5 = nan
Newton iteration 4energy = nan<A u 4 , A u 4 >_{-1}^0.5 = nan
Newton iteration 5energy = nan<A u 5 , A u 5 >_{-1}^0.5 = nan
Newton iteration 6energy = nan<A u 6 , A u 6 >_{-1}^0.5 = nan
Newton iteration 7energy = nan<A u 7 , A u 7 >_{-1}^0.5 = nan
Newton iteration 8energy = nan<A u 8 , A u 8 >_{-1}^0.5 = nan
Newton iteration 9energy = nan<A u 9 , A u 9 >_{-1}^0.5 = nan
Newton iteration 10energy = nan<A u 10 , A u 10 >_{-1}^0.5 = nan
Newton iteration 11energy = nan<A u 11 , A u 11 >_{-1}^0.5 = nan
Newton iteration 12energy = nan<A u 12 , A u 12 >_{-1}^0.5 = nan
Newton iteration 13energy = nan<A u 13 , A u 13 >_{-1}^0.5 = nan
Newton iteration 14energy = nan<A u 14 , A u 14 >_{-1}^0.5 = nan
Newton iteration 15energy = nan<A u 15 , A u 15 >_{-1}^0.5 = nan
Newton iteration 16energy = nan<A u 16 , A u 16 >_{-1}^0.5 = nan
Newton iteration 17energy = nan<A u 17 , A u 17 >_{-1}^0.5 = nan
Newton iteration 18energy = nan<A u 18 , A u 18 >_{-1}^0.5 = nan
Newton iteration 19energy = nan<A u 19 , A u 19 >_{-1}^0.5 = nan
Newton iteration 20energy = nan<A u 20 , A u 20 >_{-1}^0.5 = nan
Newton iteration 21energy = nan<A u 21 , A u 21 >_{-1}^0.5 = nan
Newton iteration 22energy = nan<A u 22 , A u 22 >_{-1}^0.5 = nan
Newton iteration 23energy = nan<A u 23 , A u 23 >_{-1}^0.5 = nan
Newton iteration 24energy = nan<A u 24 , A u 24 >_{-1}^0.5 = nan
loadstep 31
Newton iteration 0energy = nan<A u 0 , A u 0 >_{-1}^0.5 = nan
Newton iteration 1energy = nan<A u 1 , A u 1 >_{-1}^0.5 = nan
Newton iteration 2energy = nan<A u 2 , A u 2 >_{-1}^0.5 = nan
Newton iteration 3energy = nan<A u 3 , A u 3 >_{-1}^0.5 = nan
Newton iteration 4energy = nan<A u 4 , A u 4 >_{-1}^0.5 = nan
Newton iteration 5energy = nan<A u 5 , A u 5 >_{-1}^0.5 = nan
Newton iteration 6energy = nan<A u 6 , A u 6 >_{-1}^0.5 = nan
Newton iteration 7energy = nan<A u 7 , A u 7 >_{-1}^0.5 = nan
Newton iteration 8energy = nan<A u 8 , A u 8 >_{-1}^0.5 = nan
Newton iteration 9energy = nan<A u 9 , A u 9 >_{-1}^0.5 = nan
Newton iteration 10energy = nan<A u 10 , A u 10 >_{-1}^0.5 = nan
Newton iteration 11energy = nan<A u 11 , A u 11 >_{-1}^0.5 = nan
Newton iteration 12energy = nan<A u 12 , A u 12 >_{-1}^0.5 = nan
Newton iteration 13energy = nan<A u 13 , A u 13 >_{-1}^0.5 = nan
Newton iteration 14energy = nan<A u 14 , A u 14 >_{-1}^0.5 = nan
Newton iteration 15energy = nan<A u 15 , A u 15 >_{-1}^0.5 = nan
Newton iteration 16energy = nan<A u 16 , A u 16 >_{-1}^0.5 = nan
Newton iteration 17energy = nan<A u 17 , A u 17 >_{-1}^0.5 = nan
Newton iteration 18energy = nan<A u 18 , A u 18 >_{-1}^0.5 = nan
Newton iteration 19energy = nan<A u 19 , A u 19 >_{-1}^0.5 = nan
Newton iteration 20energy = nan<A u 20 , A u 20 >_{-1}^0.5 = nan
Newton iteration 21energy = nan<A u 21 , A u 21 >_{-1}^0.5 = nan
Newton iteration 22energy = nan<A u 22 , A u 22 >_{-1}^0.5 = nan
Newton iteration 23energy = nan<A u 23 , A u 23 >_{-1}^0.5 = nan
Newton iteration 24energy = nan<A u 24 , A u 24 >_{-1}^0.5 = nan
loadstep 32
Newton iteration 0energy = nan<A u 0 , A u 0 >_{-1}^0.5 = nan
Newton iteration 1energy = nan<A u 1 , A u 1 >_{-1}^0.5 = nan
Newton iteration 2energy = nan<A u 2 , A u 2 >_{-1}^0.5 = nan
Newton iteration 3energy = nan<A u 3 , A u 3 >_{-1}^0.5 = nan
Newton iteration 4energy = nan<A u 4 , A u 4 >_{-1}^0.5 = nan
Newton iteration 5energy = nan<A u 5 , A u 5 >_{-1}^0.5 = nan
Newton iteration 6energy = nan<A u 6 , A u 6 >_{-1}^0.5 = nan
Newton iteration 7energy = nan<A u 7 , A u 7 >_{-1}^0.5 = nan
Newton iteration 8energy = nan<A u 8 , A u 8 >_{-1}^0.5 = nan
Newton iteration 9energy = nan<A u 9 , A u 9 >_{-1}^0.5 = nan
Newton iteration 10energy = nan<A u 10 , A u 10 >_{-1}^0.5 = nan
Newton iteration 11energy = nan<A u 11 , A u 11 >_{-1}^0.5 = nan
Newton iteration 12energy = nan<A u 12 , A u 12 >_{-1}^0.5 = nan
Newton iteration 13energy = nan<A u 13 , A u 13 >_{-1}^0.5 = nan
Newton iteration 14energy = nan<A u 14 , A u 14 >_{-1}^0.5 = nan
Newton iteration 15energy = nan<A u 15 , A u 15 >_{-1}^0.5 = nan
Newton iteration 16energy = nan<A u 16 , A u 16 >_{-1}^0.5 = nan
Newton iteration 17energy = nan<A u 17 , A u 17 >_{-1}^0.5 = nan
Newton iteration 18energy = nan<A u 18 , A u 18 >_{-1}^0.5 = nan
Newton iteration 19energy = nan<A u 19 , A u 19 >_{-1}^0.5 = nan
Newton iteration 20energy = nan<A u 20 , A u 20 >_{-1}^0.5 = nan
Newton iteration 21energy = nan<A u 21 , A u 21 >_{-1}^0.5 = nan
Newton iteration 22energy = nan<A u 22 , A u 22 >_{-1}^0.5 = nan
Newton iteration 23energy = nan<A u 23 , A u 23 >_{-1}^0.5 = nan
Newton iteration 24energy = nan<A u 24 , A u 24 >_{-1}^0.5 = nan
loadstep 33
Newton iteration 0energy = nan<A u 0 , A u 0 >_{-1}^0.5 = nan
Newton iteration 1energy = nan<A u 1 , A u 1 >_{-1}^0.5 = nan
Newton iteration 2energy = nan<A u 2 , A u 2 >_{-1}^0.5 = nan
Newton iteration 3energy = nan<A u 3 , A u 3 >_{-1}^0.5 = nan
Newton iteration 4energy = nan<A u 4 , A u 4 >_{-1}^0.5 = nan
Newton iteration 5energy = nan<A u 5 , A u 5 >_{-1}^0.5 = nan
Newton iteration 6energy = nan<A u 6 , A u 6 >_{-1}^0.5 = nan
Newton iteration 7energy = nan<A u 7 , A u 7 >_{-1}^0.5 = nan
Newton iteration 8energy = nan<A u 8 , A u 8 >_{-1}^0.5 = nan
Newton iteration 9energy = nan<A u 9 , A u 9 >_{-1}^0.5 = nan
Newton iteration 10energy = nan<A u 10 , A u 10 >_{-1}^0.5 = nan
Newton iteration 11energy = nan<A u 11 , A u 11 >_{-1}^0.5 = nan
Newton iteration 12energy = nan<A u 12 , A u 12 >_{-1}^0.5 = nan
Newton iteration 13energy = nan<A u 13 , A u 13 >_{-1}^0.5 = nan
Newton iteration 14energy = nan<A u 14 , A u 14 >_{-1}^0.5 = nan
Newton iteration 15energy = nan<A u 15 , A u 15 >_{-1}^0.5 = nan
Newton iteration 16energy = nan<A u 16 , A u 16 >_{-1}^0.5 = nan
Newton iteration 17energy = nan<A u 17 , A u 17 >_{-1}^0.5 = nan
Newton iteration 18energy = nan<A u 18 , A u 18 >_{-1}^0.5 = nan
Newton iteration 19energy = nan<A u 19 , A u 19 >_{-1}^0.5 = nan
Newton iteration 20energy = nan<A u 20 , A u 20 >_{-1}^0.5 = nan
Newton iteration 21energy = nan<A u 21 , A u 21 >_{-1}^0.5 = nan
Newton iteration 22energy = nan<A u 22 , A u 22 >_{-1}^0.5 = nan
Newton iteration 23energy = nan<A u 23 , A u 23 >_{-1}^0.5 = nan
Newton iteration 24energy = nan<A u 24 , A u 24 >_{-1}^0.5 = nan
loadstep 34
Newton iteration 0energy = nan<A u 0 , A u 0 >_{-1}^0.5 = nan
Newton iteration 1energy = nan<A u 1 , A u 1 >_{-1}^0.5 = nan
Newton iteration 2energy = nan<A u 2 , A u 2 >_{-1}^0.5 = nan
Newton iteration 3energy = nan<A u 3 , A u 3 >_{-1}^0.5 = nan
Newton iteration 4energy = nan<A u 4 , A u 4 >_{-1}^0.5 = nan
Newton iteration 5energy = nan<A u 5 , A u 5 >_{-1}^0.5 = nan
Newton iteration 6energy = nan<A u 6 , A u 6 >_{-1}^0.5 = nan
Newton iteration 7energy = nan<A u 7 , A u 7 >_{-1}^0.5 = nan
Newton iteration 8energy = nan<A u 8 , A u 8 >_{-1}^0.5 = nan
Newton iteration 9energy = nan<A u 9 , A u 9 >_{-1}^0.5 = nan
Newton iteration 10energy = nan<A u 10 , A u 10 >_{-1}^0.5 = nan
Newton iteration 11energy = nan<A u 11 , A u 11 >_{-1}^0.5 = nan
Newton iteration 12energy = nan<A u 12 , A u 12 >_{-1}^0.5 = nan
Newton iteration 13energy = nan<A u 13 , A u 13 >_{-1}^0.5 = nan
Newton iteration 14energy = nan<A u 14 , A u 14 >_{-1}^0.5 = nan
Newton iteration 15energy = nan<A u 15 , A u 15 >_{-1}^0.5 = nan
Newton iteration 16energy = nan<A u 16 , A u 16 >_{-1}^0.5 = nan
Newton iteration 17energy = nan<A u 17 , A u 17 >_{-1}^0.5 = nan
Newton iteration 18energy = nan<A u 18 , A u 18 >_{-1}^0.5 = nan
Newton iteration 19energy = nan<A u 19 , A u 19 >_{-1}^0.5 = nan
Newton iteration 20energy = nan<A u 20 , A u 20 >_{-1}^0.5 = nan
Newton iteration 21energy = nan<A u 21 , A u 21 >_{-1}^0.5 = nan
Newton iteration 22energy = nan<A u 22 , A u 22 >_{-1}^0.5 = nan
Newton iteration 23energy = nan<A u 23 , A u 23 >_{-1}^0.5 = nan
Newton iteration 24energy = nan<A u 24 , A u 24 >_{-1}^0.5 = nan
loadstep 35
Newton iteration 0energy = nan<A u 0 , A u 0 >_{-1}^0.5 = nan
Newton iteration 1energy = nan<A u 1 , A u 1 >_{-1}^0.5 = nan
Newton iteration 2energy = nan<A u 2 , A u 2 >_{-1}^0.5 = nan
Newton iteration 3energy = nan<A u 3 , A u 3 >_{-1}^0.5 = nan
Newton iteration 4energy = nan<A u 4 , A u 4 >_{-1}^0.5 = nan
Newton iteration 5energy = nan<A u 5 , A u 5 >_{-1}^0.5 = nan
Newton iteration 6energy = nan<A u 6 , A u 6 >_{-1}^0.5 = nan
Newton iteration 7energy = nan<A u 7 , A u 7 >_{-1}^0.5 = nan
Newton iteration 8energy = nan<A u 8 , A u 8 >_{-1}^0.5 = nan
Newton iteration 9energy = nan<A u 9 , A u 9 >_{-1}^0.5 = nan
Newton iteration 10energy = nan<A u 10 , A u 10 >_{-1}^0.5 = nan
Newton iteration 11energy = nan<A u 11 , A u 11 >_{-1}^0.5 = nan
Newton iteration 12energy = nan<A u 12 , A u 12 >_{-1}^0.5 = nan
Newton iteration 13energy = nan<A u 13 , A u 13 >_{-1}^0.5 = nan
Newton iteration 14energy = nan<A u 14 , A u 14 >_{-1}^0.5 = nan
Newton iteration 15energy = nan<A u 15 , A u 15 >_{-1}^0.5 = nan
Newton iteration 16energy = nan<A u 16 , A u 16 >_{-1}^0.5 = nan
Newton iteration 17energy = nan<A u 17 , A u 17 >_{-1}^0.5 = nan
Newton iteration 18energy = nan<A u 18 , A u 18 >_{-1}^0.5 = nan
Newton iteration 19energy = nan<A u 19 , A u 19 >_{-1}^0.5 = nan
Newton iteration 20energy = nan<A u 20 , A u 20 >_{-1}^0.5 = nan
Newton iteration 21energy = nan<A u 21 , A u 21 >_{-1}^0.5 = nan
Newton iteration 22energy = nan<A u 22 , A u 22 >_{-1}^0.5 = nan
Newton iteration 23energy = nan<A u 23 , A u 23 >_{-1}^0.5 = nan
Newton iteration 24energy = nan<A u 24 , A u 24 >_{-1}^0.5 = nan
loadstep 36
Newton iteration 0energy = nan<A u 0 , A u 0 >_{-1}^0.5 = nan
Newton iteration 1energy = nan<A u 1 , A u 1 >_{-1}^0.5 = nan
Newton iteration 2energy = nan<A u 2 , A u 2 >_{-1}^0.5 = nan
Newton iteration 3energy = nan<A u 3 , A u 3 >_{-1}^0.5 = nan
Newton iteration 4energy = nan<A u 4 , A u 4 >_{-1}^0.5 = nan
Newton iteration 5energy = nan<A u 5 , A u 5 >_{-1}^0.5 = nan
Newton iteration 6energy = nan<A u 6 , A u 6 >_{-1}^0.5 = nan
Newton iteration 7energy = nan<A u 7 , A u 7 >_{-1}^0.5 = nan
Newton iteration 8energy = nan<A u 8 , A u 8 >_{-1}^0.5 = nan
Newton iteration 9energy = nan<A u 9 , A u 9 >_{-1}^0.5 = nan
Newton iteration 10energy = nan<A u 10 , A u 10 >_{-1}^0.5 = nan
Newton iteration 11energy = nan<A u 11 , A u 11 >_{-1}^0.5 = nan
Newton iteration 12energy = nan<A u 12 , A u 12 >_{-1}^0.5 = nan
Newton iteration 13energy = nan<A u 13 , A u 13 >_{-1}^0.5 = nan
Newton iteration 14energy = nan<A u 14 , A u 14 >_{-1}^0.5 = nan
Newton iteration 15energy = nan<A u 15 , A u 15 >_{-1}^0.5 = nan
Newton iteration 16energy = nan<A u 16 , A u 16 >_{-1}^0.5 = nan
Newton iteration 17energy = nan<A u 17 , A u 17 >_{-1}^0.5 = nan
Newton iteration 18energy = nan<A u 18 , A u 18 >_{-1}^0.5 = nan
Newton iteration 19energy = nan<A u 19 , A u 19 >_{-1}^0.5 = nan
Newton iteration 20energy = nan<A u 20 , A u 20 >_{-1}^0.5 = nan
Newton iteration 21energy = nan<A u 21 , A u 21 >_{-1}^0.5 = nan
Newton iteration 22energy = nan<A u 22 , A u 22 >_{-1}^0.5 = nan
Newton iteration 23energy = nan<A u 23 , A u 23 >_{-1}^0.5 = nan
Newton iteration 24energy = nan<A u 24 , A u 24 >_{-1}^0.5 = nan
loadstep 37
Newton iteration 0energy = nan<A u 0 , A u 0 >_{-1}^0.5 = nan
Newton iteration 1energy = nan<A u 1 , A u 1 >_{-1}^0.5 = nan
Newton iteration 2energy = nan<A u 2 , A u 2 >_{-1}^0.5 = nan
Newton iteration 3energy = nan<A u 3 , A u 3 >_{-1}^0.5 = nan
Newton iteration 4energy = nan<A u 4 , A u 4 >_{-1}^0.5 = nan
Newton iteration 5energy = nan<A u 5 , A u 5 >_{-1}^0.5 = nan
Newton iteration 6energy = nan<A u 6 , A u 6 >_{-1}^0.5 = nan
Newton iteration 7energy = nan<A u 7 , A u 7 >_{-1}^0.5 = nan
Newton iteration 8energy = nan<A u 8 , A u 8 >_{-1}^0.5 = nan
Newton iteration 9energy = nan<A u 9 , A u 9 >_{-1}^0.5 = nan
Newton iteration 10energy = nan<A u 10 , A u 10 >_{-1}^0.5 = nan
Newton iteration 11energy = nan<A u 11 , A u 11 >_{-1}^0.5 = nan
Newton iteration 12energy = nan<A u 12 , A u 12 >_{-1}^0.5 = nan
Newton iteration 13energy = nan<A u 13 , A u 13 >_{-1}^0.5 = nan
Newton iteration 14energy = nan<A u 14 , A u 14 >_{-1}^0.5 = nan
Newton iteration 15energy = nan<A u 15 , A u 15 >_{-1}^0.5 = nan
Newton iteration 16energy = nan<A u 16 , A u 16 >_{-1}^0.5 = nan
Newton iteration 17energy = nan<A u 17 , A u 17 >_{-1}^0.5 = nan
Newton iteration 18energy = nan<A u 18 , A u 18 >_{-1}^0.5 = nan
Newton iteration 19energy = nan<A u 19 , A u 19 >_{-1}^0.5 = nan
Newton iteration 20energy = nan<A u 20 , A u 20 >_{-1}^0.5 = nan
Newton iteration 21energy = nan<A u 21 , A u 21 >_{-1}^0.5 = nan
Newton iteration 22energy = nan<A u 22 , A u 22 >_{-1}^0.5 = nan
Newton iteration 23energy = nan<A u 23 , A u 23 >_{-1}^0.5 = nan
Newton iteration 24energy = nan<A u 24 , A u 24 >_{-1}^0.5 = nan
loadstep 38
Newton iteration 0energy = nan<A u 0 , A u 0 >_{-1}^0.5 = nan
Newton iteration 1energy = nan<A u 1 , A u 1 >_{-1}^0.5 = nan
Newton iteration 2energy = nan<A u 2 , A u 2 >_{-1}^0.5 = nan
Newton iteration 3energy = nan<A u 3 , A u 3 >_{-1}^0.5 = nan
Newton iteration 4energy = nan<A u 4 , A u 4 >_{-1}^0.5 = nan
Newton iteration 5energy = nan<A u 5 , A u 5 >_{-1}^0.5 = nan
Newton iteration 6energy = nan<A u 6 , A u 6 >_{-1}^0.5 = nan
Newton iteration 7energy = nan<A u 7 , A u 7 >_{-1}^0.5 = nan
Newton iteration 8energy = nan<A u 8 , A u 8 >_{-1}^0.5 = nan
Newton iteration 9energy = nan<A u 9 , A u 9 >_{-1}^0.5 = nan
Newton iteration 10energy = nan<A u 10 , A u 10 >_{-1}^0.5 = nan
Newton iteration 11energy = nan<A u 11 , A u 11 >_{-1}^0.5 = nan
Newton iteration 12energy = nan<A u 12 , A u 12 >_{-1}^0.5 = nan
Newton iteration 13energy = nan<A u 13 , A u 13 >_{-1}^0.5 = nan
Newton iteration 14energy = nan<A u 14 , A u 14 >_{-1}^0.5 = nan
Newton iteration 15energy = nan<A u 15 , A u 15 >_{-1}^0.5 = nan
Newton iteration 16energy = nan<A u 16 , A u 16 >_{-1}^0.5 = nan
Newton iteration 17energy = nan<A u 17 , A u 17 >_{-1}^0.5 = nan
Newton iteration 18energy = nan<A u 18 , A u 18 >_{-1}^0.5 = nan
Newton iteration 19energy = nan<A u 19 , A u 19 >_{-1}^0.5 = nan
Newton iteration 20energy = nan<A u 20 , A u 20 >_{-1}^0.5 = nan
Newton iteration 21energy = nan<A u 21 , A u 21 >_{-1}^0.5 = nan
Newton iteration 22energy = nan<A u 22 , A u 22 >_{-1}^0.5 = nan
Newton iteration 23energy = nan<A u 23 , A u 23 >_{-1}^0.5 = nan
Newton iteration 24energy = nan<A u 24 , A u 24 >_{-1}^0.5 = nan
loadstep 39
Newton iteration 0energy = nan<A u 0 , A u 0 >_{-1}^0.5 = nan
Newton iteration 1energy = nan<A u 1 , A u 1 >_{-1}^0.5 = nan
Newton iteration 2energy = nan<A u 2 , A u 2 >_{-1}^0.5 = nan
Newton iteration 3energy = nan<A u 3 , A u 3 >_{-1}^0.5 = nan
Newton iteration 4energy = nan<A u 4 , A u 4 >_{-1}^0.5 = nan
Newton iteration 5energy = nan<A u 5 , A u 5 >_{-1}^0.5 = nan
Newton iteration 6energy = nan<A u 6 , A u 6 >_{-1}^0.5 = nan
Newton iteration 7energy = nan<A u 7 , A u 7 >_{-1}^0.5 = nan
Newton iteration 8energy = nan<A u 8 , A u 8 >_{-1}^0.5 = nan
Newton iteration 9energy = nan<A u 9 , A u 9 >_{-1}^0.5 = nan
Newton iteration 10energy = nan<A u 10 , A u 10 >_{-1}^0.5 = nan
Newton iteration 11energy = nan<A u 11 , A u 11 >_{-1}^0.5 = nan
Newton iteration 12energy = nan<A u 12 , A u 12 >_{-1}^0.5 = nan
Newton iteration 13energy = nan<A u 13 , A u 13 >_{-1}^0.5 = nan
Newton iteration 14energy = nan<A u 14 , A u 14 >_{-1}^0.5 = nan
Newton iteration 15energy = nan<A u 15 , A u 15 >_{-1}^0.5 = nan
Newton iteration 16energy = nan<A u 16 , A u 16 >_{-1}^0.5 = nan
Newton iteration 17energy = nan<A u 17 , A u 17 >_{-1}^0.5 = nan
Newton iteration 18energy = nan<A u 18 , A u 18 >_{-1}^0.5 = nan
Newton iteration 19energy = nan<A u 19 , A u 19 >_{-1}^0.5 = nan
Newton iteration 20energy = nan<A u 20 , A u 20 >_{-1}^0.5 = nan
Newton iteration 21energy = nan<A u 21 , A u 21 >_{-1}^0.5 = nan
Newton iteration 22energy = nan<A u 22 , A u 22 >_{-1}^0.5 = nan
Newton iteration 23energy = nan<A u 23 , A u 23 >_{-1}^0.5 = nan
Newton iteration 24energy = nan<A u 24 , A u 24 >_{-1}^0.5 = nan
loadstep 40
Newton iteration 0energy = nan<A u 0 , A u 0 >_{-1}^0.5 = nan
Newton iteration 1energy = nan<A u 1 , A u 1 >_{-1}^0.5 = nan
Newton iteration 2energy = nan<A u 2 , A u 2 >_{-1}^0.5 = nan
Newton iteration 3energy = nan<A u 3 , A u 3 >_{-1}^0.5 = nan
Newton iteration 4energy = nan<A u 4 , A u 4 >_{-1}^0.5 = nan
Newton iteration 5energy = nan<A u 5 , A u 5 >_{-1}^0.5 = nan
Newton iteration 6energy = nan<A u 6 , A u 6 >_{-1}^0.5 = nan
Newton iteration 7energy = nan<A u 7 , A u 7 >_{-1}^0.5 = nan
Newton iteration 8energy = nan<A u 8 , A u 8 >_{-1}^0.5 = nan
Newton iteration 9energy = nan<A u 9 , A u 9 >_{-1}^0.5 = nan
Newton iteration 10energy = nan<A u 10 , A u 10 >_{-1}^0.5 = nan
Newton iteration 11energy = nan<A u 11 , A u 11 >_{-1}^0.5 = nan
Newton iteration 12energy = nan<A u 12 , A u 12 >_{-1}^0.5 = nan
Newton iteration 13energy = nan<A u 13 , A u 13 >_{-1}^0.5 = nan
Newton iteration 14energy = nan<A u 14 , A u 14 >_{-1}^0.5 = nan
Newton iteration 15energy = nan<A u 15 , A u 15 >_{-1}^0.5 = nan
Newton iteration 16energy = nan<A u 16 , A u 16 >_{-1}^0.5 = nan
Newton iteration 17energy = nan<A u 17 , A u 17 >_{-1}^0.5 = nan
Newton iteration 18energy = nan<A u 18 , A u 18 >_{-1}^0.5 = nan
Newton iteration 19energy = nan<A u 19 , A u 19 >_{-1}^0.5 = nan
Newton iteration 20energy = nan<A u 20 , A u 20 >_{-1}^0.5 = nan
Newton iteration 21energy = nan<A u 21 , A u 21 >_{-1}^0.5 = nan
Newton iteration 22energy = nan<A u 22 , A u 22 >_{-1}^0.5 = nan
Newton iteration 23energy = nan<A u 23 , A u 23 >_{-1}^0.5 = nan
Newton iteration 24energy = nan<A u 24 , A u 24 >_{-1}^0.5 = nan
loadstep 41
Newton iteration 0energy = nan<A u 0 , A u 0 >_{-1}^0.5 = nan
Newton iteration 1energy = nan<A u 1 , A u 1 >_{-1}^0.5 = nan
Newton iteration 2energy = nan<A u 2 , A u 2 >_{-1}^0.5 = nan
Newton iteration 3energy = nan<A u 3 , A u 3 >_{-1}^0.5 = nan
Newton iteration 4energy = nan<A u 4 , A u 4 >_{-1}^0.5 = nan
Newton iteration 5energy = nan<A u 5 , A u 5 >_{-1}^0.5 = nan
Newton iteration 6energy = nan<A u 6 , A u 6 >_{-1}^0.5 = nan
Newton iteration 7energy = nan<A u 7 , A u 7 >_{-1}^0.5 = nan
Newton iteration 8energy = nan<A u 8 , A u 8 >_{-1}^0.5 = nan
Newton iteration 9energy = nan<A u 9 , A u 9 >_{-1}^0.5 = nan
Newton iteration 10energy = nan<A u 10 , A u 10 >_{-1}^0.5 = nan
Newton iteration 11energy = nan<A u 11 , A u 11 >_{-1}^0.5 = nan
Newton iteration 12energy = nan<A u 12 , A u 12 >_{-1}^0.5 = nan
Newton iteration 13energy = nan<A u 13 , A u 13 >_{-1}^0.5 = nan
Newton iteration 14energy = nan<A u 14 , A u 14 >_{-1}^0.5 = nan
Newton iteration 15energy = nan<A u 15 , A u 15 >_{-1}^0.5 = nan
Newton iteration 16energy = nan<A u 16 , A u 16 >_{-1}^0.5 = nan
Newton iteration 17energy = nan<A u 17 , A u 17 >_{-1}^0.5 = nan
Newton iteration 18energy = nan<A u 18 , A u 18 >_{-1}^0.5 = nan
Newton iteration 19energy = nan<A u 19 , A u 19 >_{-1}^0.5 = nan
Newton iteration 20energy = nan<A u 20 , A u 20 >_{-1}^0.5 = nan
Newton iteration 21energy = nan<A u 21 , A u 21 >_{-1}^0.5 = nan
Newton iteration 22energy = nan<A u 22 , A u 22 >_{-1}^0.5 = nan
Newton iteration 23energy = nan<A u 23 , A u 23 >_{-1}^0.5 = nan
Newton iteration 24energy = nan<A u 24 , A u 24 >_{-1}^0.5 = nan
loadstep 42
Newton iteration 0energy = nan<A u 0 , A u 0 >_{-1}^0.5 = nan
Newton iteration 1energy = nan<A u 1 , A u 1 >_{-1}^0.5 = nan
Newton iteration 2energy = nan<A u 2 , A u 2 >_{-1}^0.5 = nan
Newton iteration 3energy = nan<A u 3 , A u 3 >_{-1}^0.5 = nan
Newton iteration 4energy = nan<A u 4 , A u 4 >_{-1}^0.5 = nan
Newton iteration 5energy = nan<A u 5 , A u 5 >_{-1}^0.5 = nan
Newton iteration 6energy = nan<A u 6 , A u 6 >_{-1}^0.5 = nan
Newton iteration 7energy = nan<A u 7 , A u 7 >_{-1}^0.5 = nan
Newton iteration 8energy = nan<A u 8 , A u 8 >_{-1}^0.5 = nan
Newton iteration 9energy = nan<A u 9 , A u 9 >_{-1}^0.5 = nan
Newton iteration 10energy = nan<A u 10 , A u 10 >_{-1}^0.5 = nan
Newton iteration 11energy = nan<A u 11 , A u 11 >_{-1}^0.5 = nan
Newton iteration 12energy = nan<A u 12 , A u 12 >_{-1}^0.5 = nan
Newton iteration 13energy = nan<A u 13 , A u 13 >_{-1}^0.5 = nan
Newton iteration 14energy = nan<A u 14 , A u 14 >_{-1}^0.5 = nan
Newton iteration 15energy = nan<A u 15 , A u 15 >_{-1}^0.5 = nan
Newton iteration 16energy = nan<A u 16 , A u 16 >_{-1}^0.5 = nan
Newton iteration 17energy = nan<A u 17 , A u 17 >_{-1}^0.5 = nan
Newton iteration 18energy = nan<A u 18 , A u 18 >_{-1}^0.5 = nan
Newton iteration 19energy = nan<A u 19 , A u 19 >_{-1}^0.5 = nan
Newton iteration 20energy = nan<A u 20 , A u 20 >_{-1}^0.5 = nan
Newton iteration 21energy = nan<A u 21 , A u 21 >_{-1}^0.5 = nan
Newton iteration 22energy = nan<A u 22 , A u 22 >_{-1}^0.5 = nan
Newton iteration 23energy = nan<A u 23 , A u 23 >_{-1}^0.5 = nan
Newton iteration 24energy = nan<A u 24 , A u 24 >_{-1}^0.5 = nan
loadstep 43
Newton iteration 0energy = nan<A u 0 , A u 0 >_{-1}^0.5 = nan
Newton iteration 1energy = nan<A u 1 , A u 1 >_{-1}^0.5 = nan
Newton iteration 2energy = nan<A u 2 , A u 2 >_{-1}^0.5 = nan
Newton iteration 3energy = nan<A u 3 , A u 3 >_{-1}^0.5 = nan
Newton iteration 4energy = nan<A u 4 , A u 4 >_{-1}^0.5 = nan
Newton iteration 5energy = nan<A u 5 , A u 5 >_{-1}^0.5 = nan
Newton iteration 6energy = nan<A u 6 , A u 6 >_{-1}^0.5 = nan
Newton iteration 7energy = nan<A u 7 , A u 7 >_{-1}^0.5 = nan
Newton iteration 8energy = nan<A u 8 , A u 8 >_{-1}^0.5 = nan
Newton iteration 9energy = nan<A u 9 , A u 9 >_{-1}^0.5 = nan
Newton iteration 10energy = nan<A u 10 , A u 10 >_{-1}^0.5 = nan
Newton iteration 11energy = nan<A u 11 , A u 11 >_{-1}^0.5 = nan
Newton iteration 12energy = nan<A u 12 , A u 12 >_{-1}^0.5 = nan
Newton iteration 13energy = nan<A u 13 , A u 13 >_{-1}^0.5 = nan
Newton iteration 14energy = nan<A u 14 , A u 14 >_{-1}^0.5 = nan
Newton iteration 15energy = nan<A u 15 , A u 15 >_{-1}^0.5 = nan
Newton iteration 16energy = nan<A u 16 , A u 16 >_{-1}^0.5 = nan
Newton iteration 17energy = nan<A u 17 , A u 17 >_{-1}^0.5 = nan
Newton iteration 18energy = nan<A u 18 , A u 18 >_{-1}^0.5 = nan
Newton iteration 19energy = nan<A u 19 , A u 19 >_{-1}^0.5 = nan
Newton iteration 20energy = nan<A u 20 , A u 20 >_{-1}^0.5 = nan
Newton iteration 21energy = nan<A u 21 , A u 21 >_{-1}^0.5 = nan
Newton iteration 22energy = nan<A u 22 , A u 22 >_{-1}^0.5 = nan
Newton iteration 23energy = nan<A u 23 , A u 23 >_{-1}^0.5 = nan
Newton iteration 24energy = nan<A u 24 , A u 24 >_{-1}^0.5 = nan
loadstep 44
Newton iteration 0energy = nan<A u 0 , A u 0 >_{-1}^0.5 = nan
Newton iteration 1energy = nan<A u 1 , A u 1 >_{-1}^0.5 = nan
Newton iteration 2energy = nan<A u 2 , A u 2 >_{-1}^0.5 = nan
Newton iteration 3energy = nan<A u 3 , A u 3 >_{-1}^0.5 = nan
Newton iteration 4energy = nan<A u 4 , A u 4 >_{-1}^0.5 = nan
Newton iteration 5energy = nan<A u 5 , A u 5 >_{-1}^0.5 = nan
Newton iteration 6energy = nan<A u 6 , A u 6 >_{-1}^0.5 = nan
Newton iteration 7energy = nan<A u 7 , A u 7 >_{-1}^0.5 = nan
Newton iteration 8energy = nan<A u 8 , A u 8 >_{-1}^0.5 = nan
Newton iteration 9energy = nan<A u 9 , A u 9 >_{-1}^0.5 = nan
Newton iteration 10energy = nan<A u 10 , A u 10 >_{-1}^0.5 = nan
Newton iteration 11energy = nan<A u 11 , A u 11 >_{-1}^0.5 = nan
Newton iteration 12energy = nan<A u 12 , A u 12 >_{-1}^0.5 = nan
Newton iteration 13energy = nan<A u 13 , A u 13 >_{-1}^0.5 = nan
Newton iteration 14energy = nan<A u 14 , A u 14 >_{-1}^0.5 = nan
Newton iteration 15energy = nan<A u 15 , A u 15 >_{-1}^0.5 = nan
Newton iteration 16energy = nan<A u 16 , A u 16 >_{-1}^0.5 = nan
Newton iteration 17energy = nan<A u 17 , A u 17 >_{-1}^0.5 = nan
Newton iteration 18energy = nan<A u 18 , A u 18 >_{-1}^0.5 = nan
Newton iteration 19energy = nan<A u 19 , A u 19 >_{-1}^0.5 = nan
Newton iteration 20energy = nan<A u 20 , A u 20 >_{-1}^0.5 = nan
Newton iteration 21energy = nan<A u 21 , A u 21 >_{-1}^0.5 = nan
Newton iteration 22energy = nan<A u 22 , A u 22 >_{-1}^0.5 = nan
Newton iteration 23energy = nan<A u 23 , A u 23 >_{-1}^0.5 = nan
Newton iteration 24energy = nan<A u 24 , A u 24 >_{-1}^0.5 = nan
loadstep 45
Newton iteration 0energy = nan<A u 0 , A u 0 >_{-1}^0.5 = nan
Newton iteration 1energy = nan<A u 1 , A u 1 >_{-1}^0.5 = nan
Newton iteration 2energy = nan<A u 2 , A u 2 >_{-1}^0.5 = nan
Newton iteration 3energy = nan<A u 3 , A u 3 >_{-1}^0.5 = nan
Newton iteration 4energy = nan<A u 4 , A u 4 >_{-1}^0.5 = nan
Newton iteration 5energy = nan<A u 5 , A u 5 >_{-1}^0.5 = nan
Newton iteration 6energy = nan<A u 6 , A u 6 >_{-1}^0.5 = nan
Newton iteration 7energy = nan<A u 7 , A u 7 >_{-1}^0.5 = nan
Newton iteration 8energy = nan<A u 8 , A u 8 >_{-1}^0.5 = nan
Newton iteration 9energy = nan<A u 9 , A u 9 >_{-1}^0.5 = nan
Newton iteration 10energy = nan<A u 10 , A u 10 >_{-1}^0.5 = nan
Newton iteration 11energy = nan<A u 11 , A u 11 >_{-1}^0.5 = nan
Newton iteration 12energy = nan<A u 12 , A u 12 >_{-1}^0.5 = nan
Newton iteration 13energy = nan<A u 13 , A u 13 >_{-1}^0.5 = nan
Newton iteration 14energy = nan<A u 14 , A u 14 >_{-1}^0.5 = nan
Newton iteration 15energy = nan<A u 15 , A u 15 >_{-1}^0.5 = nan
Newton iteration 16energy = nan<A u 16 , A u 16 >_{-1}^0.5 = nan
Newton iteration 17energy = nan<A u 17 , A u 17 >_{-1}^0.5 = nan
Newton iteration 18energy = nan<A u 18 , A u 18 >_{-1}^0.5 = nan
Newton iteration 19energy = nan<A u 19 , A u 19 >_{-1}^0.5 = nan
Newton iteration 20energy = nan<A u 20 , A u 20 >_{-1}^0.5 = nan
Newton iteration 21energy = nan<A u 21 , A u 21 >_{-1}^0.5 = nan
Newton iteration 22energy = nan<A u 22 , A u 22 >_{-1}^0.5 = nan
Newton iteration 23energy = nan<A u 23 , A u 23 >_{-1}^0.5 = nan
Newton iteration 24energy = nan<A u 24 , A u 24 >_{-1}^0.5 = nan
loadstep 46
Newton iteration 0energy = nan<A u 0 , A u 0 >_{-1}^0.5 = nan
Newton iteration 1energy = nan<A u 1 , A u 1 >_{-1}^0.5 = nan
Newton iteration 2energy = nan<A u 2 , A u 2 >_{-1}^0.5 = nan
Newton iteration 3energy = nan<A u 3 , A u 3 >_{-1}^0.5 = nan
Newton iteration 4energy = nan<A u 4 , A u 4 >_{-1}^0.5 = nan
Newton iteration 5energy = nan<A u 5 , A u 5 >_{-1}^0.5 = nan
Newton iteration 6energy = nan<A u 6 , A u 6 >_{-1}^0.5 = nan
Newton iteration 7energy = nan<A u 7 , A u 7 >_{-1}^0.5 = nan
Newton iteration 8energy = nan<A u 8 , A u 8 >_{-1}^0.5 = nan
Newton iteration 9energy = nan<A u 9 , A u 9 >_{-1}^0.5 = nan
Newton iteration 10energy = nan<A u 10 , A u 10 >_{-1}^0.5 = nan
Newton iteration 11energy = nan<A u 11 , A u 11 >_{-1}^0.5 = nan
Newton iteration 12energy = nan<A u 12 , A u 12 >_{-1}^0.5 = nan
Newton iteration 13energy = nan<A u 13 , A u 13 >_{-1}^0.5 = nan
Newton iteration 14energy = nan<A u 14 , A u 14 >_{-1}^0.5 = nan
Newton iteration 15energy = nan<A u 15 , A u 15 >_{-1}^0.5 = nan
Newton iteration 16energy = nan<A u 16 , A u 16 >_{-1}^0.5 = nan
Newton iteration 17energy = nan<A u 17 , A u 17 >_{-1}^0.5 = nan
Newton iteration 18energy = nan<A u 18 , A u 18 >_{-1}^0.5 = nan
Newton iteration 19energy = nan<A u 19 , A u 19 >_{-1}^0.5 = nan
Newton iteration 20energy = nan<A u 20 , A u 20 >_{-1}^0.5 = nan
Newton iteration 21energy = nan<A u 21 , A u 21 >_{-1}^0.5 = nan
Newton iteration 22energy = nan<A u 22 , A u 22 >_{-1}^0.5 = nan
Newton iteration 23energy = nan<A u 23 , A u 23 >_{-1}^0.5 = nan
Newton iteration 24energy = nan<A u 24 , A u 24 >_{-1}^0.5 = nan
loadstep 47
Newton iteration 0energy = nan<A u 0 , A u 0 >_{-1}^0.5 = nan
Newton iteration 1energy = nan<A u 1 , A u 1 >_{-1}^0.5 = nan
Newton iteration 2energy = nan<A u 2 , A u 2 >_{-1}^0.5 = nan
Newton iteration 3energy = nan<A u 3 , A u 3 >_{-1}^0.5 = nan
Newton iteration 4energy = nan<A u 4 , A u 4 >_{-1}^0.5 = nan
Newton iteration 5energy = nan<A u 5 , A u 5 >_{-1}^0.5 = nan
Newton iteration 6energy = nan<A u 6 , A u 6 >_{-1}^0.5 = nan
Newton iteration 7energy = nan<A u 7 , A u 7 >_{-1}^0.5 = nan
Newton iteration 8energy = nan<A u 8 , A u 8 >_{-1}^0.5 = nan
Newton iteration 9energy = nan<A u 9 , A u 9 >_{-1}^0.5 = nan
Newton iteration 10energy = nan<A u 10 , A u 10 >_{-1}^0.5 = nan
Newton iteration 11energy = nan<A u 11 , A u 11 >_{-1}^0.5 = nan
Newton iteration 12energy = nan<A u 12 , A u 12 >_{-1}^0.5 = nan
Newton iteration 13energy = nan<A u 13 , A u 13 >_{-1}^0.5 = nan
Newton iteration 14energy = nan<A u 14 , A u 14 >_{-1}^0.5 = nan
Newton iteration 15energy = nan<A u 15 , A u 15 >_{-1}^0.5 = nan
Newton iteration 16energy = nan<A u 16 , A u 16 >_{-1}^0.5 = nan
Newton iteration 17energy = nan<A u 17 , A u 17 >_{-1}^0.5 = nan
Newton iteration 18energy = nan<A u 18 , A u 18 >_{-1}^0.5 = nan
Newton iteration 19energy = nan<A u 19 , A u 19 >_{-1}^0.5 = nan
Newton iteration 20energy = nan<A u 20 , A u 20 >_{-1}^0.5 = nan
Newton iteration 21energy = nan<A u 21 , A u 21 >_{-1}^0.5 = nan
Newton iteration 22energy = nan<A u 22 , A u 22 >_{-1}^0.5 = nan
Newton iteration 23energy = nan<A u 23 , A u 23 >_{-1}^0.5 = nan
Newton iteration 24energy = nan<A u 24 , A u 24 >_{-1}^0.5 = nan
loadstep 48
Newton iteration 0energy = nan<A u 0 , A u 0 >_{-1}^0.5 = nan
Newton iteration 1energy = nan<A u 1 , A u 1 >_{-1}^0.5 = nan
Newton iteration 2energy = nan<A u 2 , A u 2 >_{-1}^0.5 = nan
Newton iteration 3energy = nan<A u 3 , A u 3 >_{-1}^0.5 = nan
Newton iteration 4energy = nan<A u 4 , A u 4 >_{-1}^0.5 = nan
Newton iteration 5energy = nan<A u 5 , A u 5 >_{-1}^0.5 = nan
Newton iteration 6energy = nan<A u 6 , A u 6 >_{-1}^0.5 = nan
Newton iteration 7energy = nan<A u 7 , A u 7 >_{-1}^0.5 = nan
Newton iteration 8energy = nan<A u 8 , A u 8 >_{-1}^0.5 = nan
Newton iteration 9energy = nan<A u 9 , A u 9 >_{-1}^0.5 = nan
Newton iteration 10energy = nan<A u 10 , A u 10 >_{-1}^0.5 = nan
Newton iteration 11energy = nan<A u 11 , A u 11 >_{-1}^0.5 = nan
Newton iteration 12energy = nan<A u 12 , A u 12 >_{-1}^0.5 = nan
Newton iteration 13energy = nan<A u 13 , A u 13 >_{-1}^0.5 = nan
Newton iteration 14energy = nan<A u 14 , A u 14 >_{-1}^0.5 = nan
Newton iteration 15energy = nan<A u 15 , A u 15 >_{-1}^0.5 = nan
Newton iteration 16energy = nan<A u 16 , A u 16 >_{-1}^0.5 = nan
Newton iteration 17energy = nan<A u 17 , A u 17 >_{-1}^0.5 = nan
Newton iteration 18energy = nan<A u 18 , A u 18 >_{-1}^0.5 = nan
Newton iteration 19energy = nan<A u 19 , A u 19 >_{-1}^0.5 = nan
Newton iteration 20energy = nan<A u 20 , A u 20 >_{-1}^0.5 = nan
Newton iteration 21energy = nan<A u 21 , A u 21 >_{-1}^0.5 = nan
Newton iteration 22energy = nan<A u 22 , A u 22 >_{-1}^0.5 = nan
Newton iteration 23energy = nan<A u 23 , A u 23 >_{-1}^0.5 = nan
Newton iteration 24energy = nan<A u 24 , A u 24 >_{-1}^0.5 = nan
loadstep 49
Newton iteration 0energy = nan<A u 0 , A u 0 >_{-1}^0.5 = nan
Newton iteration 1energy = nan<A u 1 , A u 1 >_{-1}^0.5 = nan
Newton iteration 2energy = nan<A u 2 , A u 2 >_{-1}^0.5 = nan
Newton iteration 3energy = nan<A u 3 , A u 3 >_{-1}^0.5 = nan
Newton iteration 4energy = nan<A u 4 , A u 4 >_{-1}^0.5 = nan
Newton iteration 5energy = nan<A u 5 , A u 5 >_{-1}^0.5 = nan
Newton iteration 6energy = nan<A u 6 , A u 6 >_{-1}^0.5 = nan
Newton iteration 7energy = nan<A u 7 , A u 7 >_{-1}^0.5 = nan
Newton iteration 8energy = nan<A u 8 , A u 8 >_{-1}^0.5 = nan
Newton iteration 9energy = nan<A u 9 , A u 9 >_{-1}^0.5 = nan
Newton iteration 10energy = nan<A u 10 , A u 10 >_{-1}^0.5 = nan
Newton iteration 11energy = nan<A u 11 , A u 11 >_{-1}^0.5 = nan
Newton iteration 12energy = nan<A u 12 , A u 12 >_{-1}^0.5 = nan
Newton iteration 13energy = nan<A u 13 , A u 13 >_{-1}^0.5 = nan
Newton iteration 14energy = nan<A u 14 , A u 14 >_{-1}^0.5 = nan
Newton iteration 15energy = nan<A u 15 , A u 15 >_{-1}^0.5 = nan
Newton iteration 16energy = nan<A u 16 , A u 16 >_{-1}^0.5 = nan
Newton iteration 17energy = nan<A u 17 , A u 17 >_{-1}^0.5 = nan
Newton iteration 18energy = nan<A u 18 , A u 18 >_{-1}^0.5 = nan
Newton iteration 19energy = nan<A u 19 , A u 19 >_{-1}^0.5 = nan
Newton iteration 20energy = nan<A u 20 , A u 20 >_{-1}^0.5 = nan
Newton iteration 21energy = nan<A u 21 , A u 21 >_{-1}^0.5 = nan
Newton iteration 22energy = nan<A u 22 , A u 22 >_{-1}^0.5 = nan
Newton iteration 23energy = nan<A u 23 , A u 23 >_{-1}^0.5 = nan
Newton iteration 24energy = nan<A u 24 , A u 24 >_{-1}^0.5 = nan
Allen-Cahn equation¶
The Allen-Cahn equations describe the process of phase separation and is the (\(L^2\)) gradient-flow equation to the energy
\[E(v) = \int_{\Omega} \varepsilon \vert \nabla v \vert^2~+~v^2(1-v^2) ~ dx\]
i.e. the solution to the Allen-Cahn equation solves
\[\partial_t u = \frac{\delta E}{\delta u}\]
The quantity \(u\) is an indicator for a phase where \(-1\) refers to one phase and \(1\) to another phase.
The equation has two driving forces:
\(u\) is pulled into one of the two minima (\(-1\) and \(1\)) of the nonlinear term \(u^2(1-u^2)\) (separation of the phases)
the diffusion term scaled with \(\varepsilon\) enforces a smooth transition between the two phases. \(\varepsilon\) determines the size of the transition layer
We use the “SymbolicEnergy” feature to formulate the energy minimization problem and combine it with an implicit Euler discretization:
\[M u^{n+1} - M u^n = \Delta t \underbrace{\frac{\delta E}{\delta u}}_{=:A(u)} (u^{n+1})\]
which we can interpreted as a nonlinear minimization problem again with the energy
\[E^{IE}(v) = \int_{\Omega} \frac{\varepsilon}{2} \vert \nabla v \vert^2~+~v^2(1-v^2) + \frac{1}{2\Delta t} \vert v - u^n \vert^2 ~ dx\]
To solve the nonlinear equation at every time step we again rely on Newton’s method.
[11]:
from ngsolve import *
from netgen.geom2d import *
periodic = SplineGeometry()
pnts = [ (0,0), (1,0), (1,1), (0,1) ]
pnums = [periodic.AppendPoint(*p) for p in pnts]
lright = periodic.Append ( ["line", pnums[0], pnums[1]],bc="periodic")
btop = periodic.Append ( ["line", pnums[1], pnums[2]], bc="periodic")
periodic.Append ( ["line", pnums[3], pnums[2]], leftdomain=0, rightdomain=1, copy=lright, bc="periodic")
periodic.Append ( ["line", pnums[0], pnums[3]], leftdomain=0, rightdomain=1, copy=btop, bc="periodic")
mesh = Mesh (periodic.GenerateMesh(maxh=0.2))
V = Periodic(H1(mesh, order=4, dirichlet=[]))
u = V.TrialFunction()
eps = 4e-3
dt = 1e-1
gfu = GridFunction(V)
gfuold = GridFunction(V)
a = BilinearForm (V, symmetric=False)
a += Variation( (eps/2*grad(u)*grad(u) + ((1-u*u)*(1-u*u))
+ 0.5/dt*(u-gfuold)*(u-gfuold)) * dx)
[12]:
from math import pi
gfu = GridFunction(V)
gfu.Set(sin(2*pi*x))
#gfu.Set(sin(1e8*x)) #<- essentially a random function
Draw(gfu,mesh,"u")
SetVisualization (deformation=True)
t = 0
[13]:
for timestep in range(50):
gfuold.vec.data = gfu.vec
SolveNonlinearMinProblem(a,gfu)
Redraw()
t += dt
print("t = ", t)
Newton iteration 0energy = 0.41450192074963627<A u 0 , A u 0 >_{-1}^0.5 = 0.2809270463814763
Newton iteration 1energy = 0.3768978580277295<A u 1 , A u 1 >_{-1}^0.5 = 0.019174370471939695
Newton iteration 2energy = 0.37671342104048205<A u 2 , A u 2 >_{-1}^0.5 = 9.62893713459873e-05
Newton iteration 3energy = 0.3767134164045806<A u 3 , A u 3 >_{-1}^0.5 = 2.505456355595927e-09
Newton iteration 4energy = 0.3767134164045804<A u 4 , A u 4 >_{-1}^0.5 = 2.7184233888511976e-16
t = 0.1
Newton iteration 0energy = 0.34388128652385946<A u 0 , A u 0 >_{-1}^0.5 = 0.2360393474078293
Newton iteration 1energy = 0.317109443298749<A u 1 , A u 1 >_{-1}^0.5 = 0.013442999343193845
Newton iteration 2energy = 0.31701887934635165<A u 2 , A u 2 >_{-1}^0.5 = 4.660092984276528e-05
Newton iteration 3energy = 0.3170188782605195<A u 3 , A u 3 >_{-1}^0.5 = 5.76804809164461e-10
Newton iteration 4energy = 0.3170188782605194<A u 4 , A u 4 >_{-1}^0.5 = 3.1990285570820334e-16
t = 0.2
Newton iteration 0energy = 0.29549435651384603<A u 0 , A u 0 >_{-1}^0.5 = 0.18345925931027124
Newton iteration 1energy = 0.2791640303389011<A u 1 , A u 1 >_{-1}^0.5 = 0.00801465116310893
Newton iteration 2energy = 0.27913186993623207<A u 2 , A u 2 >_{-1}^0.5 = 1.6250158077614812e-05
Newton iteration 3energy = 0.27913186980419796<A u 3 , A u 3 >_{-1}^0.5 = 6.857117114914922e-11
Newton iteration 4energy = 0.27913186980419796<A u 4 , A u 4 >_{-1}^0.5 = 3.7531322332558395e-16
t = 0.30000000000000004
Newton iteration 0energy = 0.266927895079205<A u 0 , A u 0 >_{-1}^0.5 = 0.13349658930568908
Newton iteration 1energy = 0.2582050355110706<A u 1 , A u 1 >_{-1}^0.5 = 0.004188513106581714
Newton iteration 2energy = 0.2581962576290237<A u 2 , A u 2 >_{-1}^0.5 = 4.366756818932449e-06
Newton iteration 3energy = 0.25819625761948944<A u 3 , A u 3 >_{-1}^0.5 = 4.850627189892583e-12
Newton iteration 4energy = 0.25819625761948944<A u 4 , A u 4 >_{-1}^0.5 = 3.5334040570852337e-16
t = 0.4
Newton iteration 0energy = 0.2520431142368916<A u 0 , A u 0 >_{-1}^0.5 = 0.0923617548386862
Newton iteration 1energy = 0.24783887867786272<A u 1 , A u 1 >_{-1}^0.5 = 0.0019843467589302326
Newton iteration 2energy = 0.24783690922418358<A u 2 , A u 2 >_{-1}^0.5 = 9.685217012005846e-07
Newton iteration 3energy = 0.2478369092237145<A u 3 , A u 3 >_{-1}^0.5 = 2.347626218203708e-13
Newton iteration 4energy = 0.24783690922371454<A u 4 , A u 4 >_{-1}^0.5 = 3.744224708774473e-16
t = 0.5
Newton iteration 0energy = 0.24499322447001118<A u 0 , A u 0 >_{-1}^0.5 = 0.06164684190451998
Newton iteration 1energy = 0.24311100787346082<A u 1 , A u 1 >_{-1}^0.5 = 0.000877962350648967
Newton iteration 2energy = 0.24311062240965592<A u 2 , A u 2 >_{-1}^0.5 = 1.882207876194426e-07
Newton iteration 3energy = 0.2431106224096382<A u 3 , A u 3 >_{-1}^0.5 = 8.792617315901084e-15
t = 0.6
Newton iteration 0energy = 0.24187354814793974<A u 0 , A u 0 >_{-1}^0.5 = 0.04016635875651402
Newton iteration 1energy = 0.24107181028211977<A u 1 , A u 1 >_{-1}^0.5 = 0.0003715367480244102
Newton iteration 2energy = 0.24107174125819164<A u 2 , A u 2 >_{-1}^0.5 = 3.366894131842145e-08
Newton iteration 3energy = 0.24107174125819092<A u 3 , A u 3 >_{-1}^0.5 = 4.913943409229916e-16
t = 0.7
Newton iteration 0energy = 0.24055438873764529<A u 0 , A u 0 >_{-1}^0.5 = 0.0257820338592469
Newton iteration 1energy = 0.24022333391554057<A u 1 , A u 1 >_{-1}^0.5 = 0.00015331214259508882
Newton iteration 2energy = 0.24022332216294007<A u 2 , A u 2 >_{-1}^0.5 = 5.795968633379131e-09
Newton iteration 3energy = 0.24022332216294012<A u 3 , A u 3 >_{-1}^0.5 = 4.2179709338570137e-16
t = 0.7999999999999999
Newton iteration 0energy = 0.24001186185815887<A u 0 , A u 0 >_{-1}^0.5 = 0.016425332383409815
Newton iteration 1energy = 0.23987730469401844<A u 1 , A u 1 >_{-1}^0.5 = 6.283809116190501e-05
Newton iteration 2energy = 0.23987730271968472<A u 2 , A u 2 >_{-1}^0.5 = 1.0154930806997262e-09
Newton iteration 3energy = 0.2398773027196847<A u 3 , A u 3 >_{-1}^0.5 = 4.173970267623849e-16
t = 0.8999999999999999
Newton iteration 0energy = 0.23979158193635364<A u 0 , A u 0 >_{-1}^0.5 = 0.010464445156610182
Newton iteration 1energy = 0.23973691889059892<A u 1 , A u 1 >_{-1}^0.5 = 2.620620881243269e-05
Newton iteration 2energy = 0.23973691854721466<A u 2 , A u 2 >_{-1}^0.5 = 1.999460696404716e-10
Newton iteration 3energy = 0.23973691854721466<A u 3 , A u 3 >_{-1}^0.5 = 4.4248699202570095e-16
t = 0.9999999999999999
Newton iteration 0energy = 0.23970186094053098<A u 0 , A u 0 >_{-1}^0.5 = 0.006735415826965367
Newton iteration 1energy = 0.2396792028920519<A u 1 , A u 1 >_{-1}^0.5 = 1.1556763084958016e-05
Newton iteration 2energy = 0.23967920282527239<A u 2 , A u 2 >_{-1}^0.5 = 5.0994595020143205e-11
Newton iteration 3energy = 0.23967920282527236<A u 3 , A u 3 >_{-1}^0.5 = 4.767165104918626e-16
t = 1.0999999999999999
Newton iteration 0energy = 0.23966435830447438<A u 0 , A u 0 >_{-1}^0.5 = 0.004452542361404154
Newton iteration 1energy = 0.2396544535095988<A u 1 , A u 1 >_{-1}^0.5 = 5.694142646686867e-06
Newton iteration 2energy = 0.23965445349338702<A u 2 , A u 2 >_{-1}^0.5 = 1.7725053456860548e-11
Newton iteration 3energy = 0.23965445349338704<A u 3 , A u 3 >_{-1}^0.5 = 4.904553444933838e-16
t = 1.2
Newton iteration 0energy = 0.2396476714273789<A u 0 , A u 0 >_{-1}^0.5 = 0.003097908141818161
Newton iteration 1energy = 0.2396428758025269<A u 1 , A u 1 >_{-1}^0.5 = 3.2840536759943483e-06
Newton iteration 2energy = 0.23964287579713436<A u 2 , A u 2 >_{-1}^0.5 = 7.767160647134923e-12
Newton iteration 3energy = 0.23964287579713447<A u 3 , A u 3 >_{-1}^0.5 = 4.723813054601472e-16
t = 1.3
Newton iteration 0energy = 0.23963934988363733<A u 0 , A u 0 >_{-1}^0.5 = 0.002326255780764285
Newton iteration 1energy = 0.23963664545354402<A u 1 , A u 1 >_{-1}^0.5 = 2.2181767038573016e-06
Newton iteration 2energy = 0.23963664545108385<A u 2 , A u 2 >_{-1}^0.5 = 3.98362777332673e-12
Newton iteration 3energy = 0.23963664545108396<A u 3 , A u 3 >_{-1}^0.5 = 4.4351366529535486e-16
t = 1.4000000000000001
Newton iteration 0energy = 0.2396344792626691<A u 0 , A u 0 >_{-1}^0.5 = 0.001900681681945114
Newton iteration 1energy = 0.23963267365152283<A u 1 , A u 1 >_{-1}^0.5 = 1.6879462608945738e-06
Newton iteration 2energy = 0.23963267365009808<A u 2 , A u 2 >_{-1}^0.5 = 2.298264705049454e-12
Newton iteration 3energy = 0.23963267365009816<A u 3 , A u 3 >_{-1}^0.5 = 4.4197703441051045e-16
t = 1.5000000000000002
Newton iteration 0energy = 0.2396311154900333<A u 0 , A u 0 >_{-1}^0.5 = 0.0016629280504851128
Newton iteration 1energy = 0.23962973322336534<A u 1 , A u 1 >_{-1}^0.5 = 1.3855651852114338e-06
Newton iteration 2energy = 0.2396297332224055<A u 2 , A u 2 >_{-1}^0.5 = 1.458610570239534e-12
Newton iteration 3energy = 0.23962973322240552<A u 3 , A u 3 >_{-1}^0.5 = 4.38518424272123e-16
t = 1.6000000000000003
Newton iteration 0energy = 0.23962848049094196<A u 0 , A u 0 >_{-1}^0.5 = 0.001518763204050943
Newton iteration 1energy = 0.2396273274180991<A u 1 , A u 1 >_{-1}^0.5 = 1.1887244241115454e-06
Newton iteration 2energy = 0.23962732741739273<A u 2 , A u 2 >_{-1}^0.5 = 1.001096666960293e-12
Newton iteration 3energy = 0.23962732741739273<A u 3 , A u 3 >_{-1}^0.5 = 4.0840067389417865e-16
t = 1.7000000000000004
Newton iteration 0energy = 0.23962625398674178<A u 0 , A u 0 >_{-1}^0.5 = 0.0014193843873505012
Newton iteration 1energy = 0.2396252468233091<A u 1 , A u 1 >_{-1}^0.5 = 1.0456329094495533e-06
Newton iteration 2energy = 0.23962524682276243<A u 2 , A u 2 >_{-1}^0.5 = 7.301201343513507e-13
Newton iteration 3energy = 0.23962524682276243<A u 3 , A u 3 >_{-1}^0.5 = 4.29510672098558e-16
t = 1.8000000000000005
Newton iteration 0energy = 0.23962429644294758<A u 0 , A u 0 >_{-1}^0.5 = 0.0013419591604963835
Newton iteration 1energy = 0.23962339612718064<A u 1 , A u 1 >_{-1}^0.5 = 9.330532164733575e-07
Newton iteration 2energy = 0.23962339612674535<A u 2 , A u 2 >_{-1}^0.5 = 5.563593050694301e-13
Newton iteration 3energy = 0.23962339612674535<A u 3 , A u 3 >_{-1}^0.5 = 4.650852292595641e-16
t = 1.9000000000000006
Newton iteration 0energy = 0.23962254072403216<A u 0 , A u 0 >_{-1}^0.5 = 0.0012763226584944427
Newton iteration 1energy = 0.2396217263037926<A u 1 , A u 1 >_{-1}^0.5 = 8.39903365900313e-07
Newton iteration 2energy = 0.2396217263034399<A u 2 , A u 2 >_{-1}^0.5 = 4.3695229015325914e-13
Newton iteration 3energy = 0.23962172630343986<A u 3 , A u 3 >_{-1}^0.5 = 4.4327115109145793e-16
t = 2.0000000000000004
Newton iteration 0energy = 0.23962094959676447<A u 0 , A u 0 >_{-1}^0.5 = 0.001217971596191082
Newton iteration 1energy = 0.23962020792824282<A u 1 , A u 1 >_{-1}^0.5 = 7.604668959773798e-07
Newton iteration 2energy = 0.23962020792795358<A u 2 , A u 2 >_{-1}^0.5 = 3.5034624456400284e-13
Newton iteration 3energy = 0.23962020792795355<A u 3 , A u 3 >_{-1}^0.5 = 4.2207511307987207e-16
t = 2.1000000000000005
Newton iteration 0energy = 0.23961949892644918<A u 0 , A u 0 >_{-1}^0.5 = 0.001164835876370049
Newton iteration 1energy = 0.23961882055016764<A u 1 , A u 1 >_{-1}^0.5 = 6.914942047827968e-07
Newton iteration 2energy = 0.23961882054992856<A u 2 , A u 2 >_{-1}^0.5 = 2.8492711708547686e-13
Newton iteration 3energy = 0.2396188205499285<A u 3 , A u 3 >_{-1}^0.5 = 4.105118892875785e-16
t = 2.2000000000000006
Newton iteration 0energy = 0.23961817091741755<A u 0 , A u 0 >_{-1}^0.5 = 0.0011158806899804854
Newton iteration 1energy = 0.23961754835794166<A u 1 , A u 1 >_{-1}^0.5 = 6.309368258696765e-07
Newton iteration 2energy = 0.23961754835774274<A u 2 , A u 2 >_{-1}^0.5 = 2.341796123153525e-13
Newton iteration 3energy = 0.23961754835774274<A u 3 , A u 3 >_{-1}^0.5 = 4.061828856568734e-16
t = 2.3000000000000007
Newton iteration 0energy = 0.23961695130021168<A u 0 , A u 0 >_{-1}^0.5 = 0.0010705168654729913
Newton iteration 1energy = 0.23961637832545638<A u 1 , A u 1 >_{-1}^0.5 = 5.773776869402066e-07
Newton iteration 2energy = 0.23961637832528954<A u 2 , A u 2 >_{-1}^0.5 = 1.941085721921248e-13
Newton iteration 3energy = 0.2396163783252895<A u 3 , A u 3 >_{-1}^0.5 = 3.8189491824811617e-16
t = 2.400000000000001
Newton iteration 0energy = 0.2396158280858507<A u 0 , A u 0 >_{-1}^0.5 = 0.0010283540627867129
Newton iteration 1energy = 0.23961529935305312<A u 1 , A u 1 >_{-1}^0.5 = 5.297627567459455e-07
Newton iteration 2energy = 0.2396152993529128<A u 2 , A u 2 >_{-1}^0.5 = 1.6202721330818085e-13
Newton iteration 3energy = 0.23961529935291287<A u 3 , A u 3 >_{-1}^0.5 = 4.358270738139593e-16
t = 2.500000000000001
Newton iteration 0energy = 0.23961479095503696<A u 0 , A u 0 >_{-1}^0.5 = 0.0009890970729091886
Newton iteration 1energy = 0.23961430181779736<A u 1 , A u 1 >_{-1}^0.5 = 4.872666035761809e-07
Newton iteration 2energy = 0.23961430181767868<A u 2 , A u 2 >_{-1}^0.5 = 1.3596374979037828e-13
Newton iteration 3energy = 0.23961430181767868<A u 3 , A u 3 >_{-1}^0.5 = 4.251907858981771e-16
t = 2.600000000000001
Newton iteration 0energy = 0.2396138309145446<A u 0 , A u 0 >_{-1}^0.5 = 0.0009525015914669434
Newton iteration 1energy = 0.239613377301057<A u 1 , A u 1 >_{-1}^0.5 = 4.4921981870807446e-07
Newton iteration 2energy = 0.2396133773009562<A u 2 , A u 2 >_{-1}^0.5 = 1.1467946165587828e-13
Newton iteration 3energy = 0.23961337730095622<A u 3 , A u 3 >_{-1}^0.5 = 3.774188750075078e-16
t = 2.700000000000001
Newton iteration 0energy = 0.23961294007376727<A u 0 , A u 0 >_{-1}^0.5 = 0.000918354775793347
Newton iteration 1energy = 0.23961251839967906<A u 1 , A u 1 >_{-1}^0.5 = 4.1506656580398494e-07
Newton iteration 2energy = 0.23961251839959297<A u 2 , A u 2 >_{-1}^0.5 = 9.725202360614327e-14
t = 2.800000000000001
Newton iteration 0energy = 0.23961211148122993<A u 0 , A u 0 >_{-1}^0.5 = 0.0008864662194064388
Newton iteration 1energy = 0.239611718581682<A u 1 , A u 1 >_{-1}^0.5 = 3.8433776211053425e-07
Newton iteration 2energy = 0.23961171858160824<A u 2 , A u 2 >_{-1}^0.5 = 8.28029265589098e-14
t = 2.9000000000000012
Newton iteration 0energy = 0.23961133899534895<A u 0 , A u 0 >_{-1}^0.5 = 0.0008566633749454556
Newton iteration 1energy = 0.23961097206924228<A u 1 , A u 1 >_{-1}^0.5 = 3.566329547505521e-07
Newton iteration 2energy = 0.2396109720691786<A u 2 , A u 2 >_{-1}^0.5 = 7.081354970377196e-14
t = 3.0000000000000013
Newton iteration 0energy = 0.23961061717765034<A u 0 , A u 0 >_{-1}^0.5 = 0.0008287889263658277
Newton iteration 1energy = 0.23961027374068294<A u 1 , A u 1 >_{-1}^0.5 = 3.3160737848567906e-07
Newton iteration 2energy = 0.23961027374062788<A u 2 , A u 2 >_{-1}^0.5 = 6.079803952456371e-14
t = 3.1000000000000014
Newton iteration 0energy = 0.2396099412024393<A u 0 , A u 0 >_{-1}^0.5 = 0.0008026990506801483
Newton iteration 1energy = 0.2396096190469684<A u 1 , A u 1 >_{-1}^0.5 = 3.0896230212479397e-07
Newton iteration 2energy = 0.23960961904692069<A u 2 , A u 2 >_{-1}^0.5 = 5.241140706901035e-14
t = 3.2000000000000015
Newton iteration 0energy = 0.23960930677943382<A u 0 , A u 0 >_{-1}^0.5 = 0.0007782621114855629
Newton iteration 1energy = 0.23960900393990808<A u 1 , A u 1 >_{-1}^0.5 = 2.884375804236844e-07
Newton iteration 2energy = 0.23960900393986664<A u 2 , A u 2 >_{-1}^0.5 = 4.540118001573666e-14
t = 3.3000000000000016
Newton iteration 0energy = 0.23960871008705897<A u 0 , A u 0 >_{-1}^0.5 = 0.0007553575814414011
Newton iteration 1energy = 0.23960842481012087<A u 1 , A u 1 >_{-1}^0.5 = 2.6980575458671215e-07
Newton iteration 2energy = 0.23960842481008446<A u 2 , A u 2 >_{-1}^0.5 = 3.943996348918186e-14
t = 3.4000000000000017
Newton iteration 0energy = 0.23960814771471362<A u 0 , A u 0 >_{-1}^0.5 = 0.0007338751008247079
Newton iteration 1energy = 0.23960787843327264<A u 1 , A u 1 >_{-1}^0.5 = 2.5286728189389155e-07
Newton iteration 2energy = 0.23960787843324066<A u 2 , A u 2 >_{-1}^0.5 = 3.439428926348582e-14
t = 3.5000000000000018
Newton iteration 0energy = 0.23960761661268298<A u 0 , A u 0 >_{-1}^0.5 = 0.0007137136279155078
Newton iteration 1energy = 0.23960736192339585<A u 1 , A u 1 >_{-1}^0.5 = 2.374466157262674e-07
Newton iteration 2energy = 0.23960736192336773<A u 2 , A u 2 >_{-1}^0.5 = 3.0118209899550216e-14
t = 3.600000000000002
Newton iteration 0energy = 0.23960711404861648<A u 0 , A u 0 >_{-1}^0.5 = 0.0006947806591340658
Newton iteration 1energy = 0.2396068726922975<A u 1 , A u 1 >_{-1}^0.5 = 2.2338894228518974e-07
Newton iteration 2energy = 0.23960687269227246<A u 2 , A u 2 >_{-1}^0.5 = 2.646664649776906e-14
t = 3.700000000000002
Newton iteration 0energy = 0.23960663756965972<A u 0 , A u 0 >_{-1}^0.5 = 0.0006769915073045265
Newton iteration 1energy = 0.23960640841422343<A u 1 , A u 1 >_{-1}^0.5 = 2.105574362104978e-07
Newton iteration 2energy = 0.2396064084142012<A u 2 , A u 2 >_{-1}^0.5 = 2.336972902228967e-14
t = 3.800000000000002
Newton iteration 0energy = 0.23960618496947447<A u 0 , A u 0 >_{-1}^0.5 = 0.000660268631503423
Newton iteration 1energy = 0.23960596699506742<A u 1 , A u 1 >_{-1}^0.5 = 1.988309330887977e-07
Newton iteration 2energy = 0.23960596699504766<A u 2 , A u 2 >_{-1}^0.5 = 2.066429071905498e-14
t = 3.900000000000002
Newton iteration 0energy = 0.2396057542594827<A u 0 , A u 0 >_{-1}^0.5 = 0.0006445410144735523
Newton iteration 1energy = 0.23960554654551236<A u 1 , A u 1 >_{-1}^0.5 = 1.8810194235857957e-07
Newton iteration 2energy = 0.2396055465454947<A u 2 , A u 2 >_{-1}^0.5 = 1.8366811384808692e-14
t = 4.000000000000002
Newton iteration 0energy = 0.23960534364377462<A u 0 , A u 0 >_{-1}^0.5 = 0.000629743584836064
Newton iteration 1energy = 0.2396051453575799<A u 1 , A u 1 >_{-1}^0.5 = 1.782749416632698e-07
Newton iteration 2energy = 0.23960514535756394<A u 2 , A u 2 >_{-1}^0.5 = 1.6387693854311274e-14
t = 4.100000000000001
Newton iteration 0energy = 0.2396049514971884<A u 0 , A u 0 >_{-1}^0.5 = 0.0006158166819334044
Newton iteration 1energy = 0.23960476188413699<A u 1 , A u 1 >_{-1}^0.5 = 1.6926490630063148e-07
Newton iteration 2energy = 0.23960476188412266<A u 2 , A u 2 >_{-1}^0.5 = 1.4645667456707844e-14
t = 4.200000000000001
Newton iteration 0energy = 0.23960457634614388<A u 0 , A u 0 >_{-1}^0.5 = 0.0006027055613948682
Newton iteration 1energy = 0.23960439472096542<A u 1 , A u 1 >_{-1}^0.5 = 1.6099603651044563e-07
Newton iteration 2energy = 0.23960439472095246<A u 2 , A u 2 >_{-1}^0.5 = 1.3169049709341755e-14
t = 4.300000000000001
Newton iteration 0energy = 0.2396042168518652<A u 0 , A u 0 >_{-1}^0.5 = 0.0005903599396007538
Newton iteration 1energy = 0.23960404259105916<A u 1 , A u 1 >_{-1}^0.5 = 1.5340065262933946e-07
Newton iteration 2energy = 0.23960404259104737<A u 2 , A u 2 >_{-1}^0.5 = 1.1870290065198425e-14
t = 4.4
Newton iteration 0energy = 0.23960387179567857<A u 0 , A u 0 >_{-1}^0.5 = 0.000578733575219527
Newton iteration 1energy = 0.23960370433085457<A u 1 , A u 1 >_{-1}^0.5 = 1.464182337184534e-07
Newton iteration 2energy = 0.23960370433084374<A u 2 , A u 2 >_{-1}^0.5 = 1.0689611591568492e-14
t = 4.5
Newton iteration 0energy = 0.23960354006611326<A u 0 , A u 0 >_{-1}^0.5 = 0.0005677838859563624
Newton iteration 1energy = 0.23960337887814326<A u 1 , A u 1 >_{-1}^0.5 = 1.3999457887440088e-07
Newton iteration 2energy = 0.23960337887813346<A u 2 , A u 2 >_{-1}^0.5 = 9.775614352617596e-15
t = 4.6
Newton iteration 0energy = 0.23960322064757067<A u 0 , A u 0 >_{-1}^0.5 = 0.0005574715986083124
Newton iteration 1energy = 0.23960306526144648<A u 1 , A u 1 >_{-1}^0.5 = 1.340810746787933e-07
Newton iteration 2energy = 0.23960306526143751<A u 2 , A u 2 >_{-1}^0.5 = 8.915450208600919e-15
t = 4.699999999999999
Newton iteration 0energy = 0.23960291261035663<A u 0 , A u 0 >_{-1}^0.5 = 0.0005477604304871074
Newton iteration 1energy = 0.2396027625906619<A u 1 , A u 1 >_{-1}^0.5 = 1.286340534114446e-07
Newton iteration 2energy = 0.23960276259065363<A u 2 , A u 2 >_{-1}^0.5 = 8.152953746466837e-15
t = 4.799999999999999
Newton iteration 0energy = 0.2396026151019025<A u 0 , A u 0 >_{-1}^0.5 = 0.0005386168002517349
Newton iteration 1energy = 0.2396024700488198<A u 1 , A u 1 >_{-1}^0.5 = 1.236142308349331e-07
Newton iteration 2energy = 0.23960247004881194<A u 2 , A u 2 >_{-1}^0.5 = 7.477509263058163e-15
t = 4.899999999999999
Newton iteration 0energy = 0.23960232733902<A u 0 , A u 0 >_{-1}^0.5 = 0.0005300095661930656
Newton iteration 1energy = 0.2396021868848036<A u 1 , A u 1 >_{-1}^0.5 = 1.189862127140539e-07
Newton iteration 2energy = 0.23960218688479662<A u 2 , A u 2 >_{-1}^0.5 = 6.8851488394343004e-15
t = 4.999999999999998
Minimal energy extension (postscript in unit-2.1.3 )¶
\[u \in V^{ho,disc}, \quad u^{lo,cont} \in V^{lo,cont}, \quad \lambda \in V^{lo,disc},\]
[14]:
from netgen.geom2d import unit_square
mesh = Mesh(unit_square.GenerateMesh(maxh=0.1))
fes_ho = Discontinuous(H1(mesh, order=10))
fes_lo = H1(mesh, order=1, dirichlet=".*")
fes_lam = Discontinuous(H1(mesh, order=1))
fes = FESpace([fes_ho, fes_lo, fes_lam])
uho, ulo, lam = fes.TrialFunction()
\[\int_{\Omega} \frac12 \Vert \nabla u \Vert^2 - u + \sum_T \sum_{V \in V(T)} ((u-u^{lo})\cdot \lambda)|_{V} \longrightarrow \operatorname{min}!\]
[15]:
a = BilinearForm(fes)
a += Variation(0.5 * grad(uho)*grad(uho)*dx
- 1*uho*dx
+ (uho-ulo)*lam*dx(element_vb=BBND))
gfu = GridFunction(fes)
solvers.Newton(a=a, u=gfu)
Draw(gfu.components[0])
Newton iteration 0
err = 0.3937351156371161
Newton iteration 1
err = 3.874027236775335e-15
The minimization problem is solved by the solution of the PDE:
\[\int_{\Omega} \nabla u \cdot \nabla v = \int_{\Omega} 1 \cdot v \quad \forall ~ v \in V^{ho,disc}\]
under the constraint
\[u(v) = u^{lo}(v) \quad \text{ for all vertices } v \in V(T) \text{ for all } T.\]