Nonlinear Naghdi shell with arc-length continuation

This demo traces the load-displacement response of a thin elastic cylindrical roof panel subjected to a downward point load at its crown — a classical benchmark for geometric nonlinearity in shells Sze et al., 2004. The panel undergoes snap-through (and snap-back), so arc-length continuation is used. The mechanical model uses five-parameter Naghdi kinematics with partial selective reduced integration. The spatial formulation employs the tangential differential calculus (TDC) framework Hansbo et al., 2015Schöllhammer & Fries, 2019Neunteufel, 2021 (§7.2.1 of the latter).

Boundary conditions: the straight side edges are hinged (zero displacement, free rotation). The curved axial-end edges are free.

In particular, this demo emphasizes:

Naghdi shell kinematics with three Euler-angle director rotations,
through-thickness decomposition of Green-Lagrange strains and PK2 stresses into membrane, bending and shear contributions,
Crisfield Crisfield, 1981 arc-length continuation past a snap-through limit point including snap-back,
partial selective reduced integration to suppress membrane and shear locking,
a Lagrange-multiplier stabilisation of the unphysical drilling rotation,
validation against the digitised reference curve of Sze et al., 2004.


from mpi4py import MPI
from petsc4py import PETSc

import basix
import dolfinx
import ufl

import matplotlib.pyplot as plt
import mesh_opencylinder_gmshapi as mg
import numpy as np
import pyvista

import dolfiny

name = "tdc_shell_naghdi_cylindrical_roof"
comm = MPI.COMM_WORLD

Geometry and reference configuration¶

The benchmark geometry (Fig. 9a of Sze et al., 2004) is a thin cylindrical roof panel with radius $R = 2540$ , axial half-length $L = 254$ , half-opening angle $\theta = 0.1$ , and thickness $h = 6.35$ . Material parameters are $E = 3102.75$ and $\nu = 0.3$ . The reference load is $P = 1$ , with $\lambda$ the applied force.

The full panel spans $\theta \in [-0.1, 0.1]$ circumferentially and $y \in [0, 2L]$ axially, with the crown ridge along $\theta = 0$ . The point load acts at the midpoint of that ridge, topmid. Straight side edges (sides at $\theta = \pm 0.1$ ) are hinged. The arc-shaped axial ends (front, rear) are free.

Geometry of the cylindrical roof panel (Fig. 9a of Sze 2004). — Figure 1:Geometry of the cylindrical roof panel, reproduced from Sze *et al.*, 2004.

R = 2540  # cylinder radius
L = 254  # axial half-length
h = 6.35  # shell thickness
θ = 0.1  # half-opening angle
N = 11  # nodes per edge
p = 2  # polynomial order, physics
q = 2  # polynomial order, geometry
do_quads = False

gmsh_model, tdim = mg.mesh_opencylinder_gmshapi(
    name, Ly=L, R=R, θ=θ, nL=N, nR=N, order=q, do_quads=do_quads
)

mesh_data = dolfinx.io.gmsh.model_to_mesh(gmsh_model, comm, 0)
mesh = mesh_data.mesh
mts = mesh_data.cell_tags
gdim = mesh.geometry.dim

Info    : Meshing 1D...
Info    : [  0%] Meshing curve 1 (Circle)
Info    : [ 10%] Meshing curve 2 (Circle)
Info    : [ 20%] Meshing curve 3 (Circle)
Info    : [ 30%] Meshing curve 4 (Circle)
Info    : [ 40%] Meshing curve 5 (Circle)
Info    : [ 50%] Meshing curve 6 (Circle)
Info    : [ 60%] Meshing curve 7 (Line)
Info    : [ 60%] Meshing curve 8 (Line)
Info    : [ 70%] Meshing curve 9 (Line)
Info    : [ 80%] Meshing curve 10 (Line)
Info    : [ 90%] Meshing curve 11 (Line)
Info    : [100%] Meshing curve 12 (Line)
Info    : Done meshing 1D (Wall 0.000491684s, CPU 0.000881s)
Info    : Meshing 2D...
Info    : [  0%] Meshing surface 1 (Surface, Frontal-Delaunay)
Info    : [ 30%] Meshing surface 2 (Surface, Frontal-Delaunay)
Info    : [ 60%] Meshing surface 3 (Surface, Frontal-Delaunay)
Info    : [ 80%] Meshing surface 4 (Surface, Frontal-Delaunay)
Info    : Done meshing 2D (Wall 0.0168567s, CPU 0.017269s)
Info    : Meshing 3D...
Info    : Done meshing 3D (Wall 7.61263e-06s, CPU 7e-06s)
Info    : 528 nodes 1100 elements
Info    : Meshing order 2 (curvilinear on)...
Info    : [  0%] Meshing curve 1 order 2
Info    : [ 10%] Meshing curve 2 order 2
Info    : [ 20%] Meshing curve 3 order 2
Info    : [ 20%] Meshing curve 4 order 2
Info    : [ 30%] Meshing curve 5 order 2
Info    : [ 40%] Meshing curve 6 order 2
Info    : [ 40%] Meshing curve 7 order 2
Info    : [ 50%] Meshing curve 8 order 2
Info    : [ 60%] Meshing curve 9 order 2
Info    : [ 60%] Meshing curve 10 order 2
Info    : [ 70%] Meshing curve 11 order 2
Info    : [ 70%] Meshing curve 12 order 2
Info    : [ 80%] Meshing surface 1 order 2
Info    : [ 90%] Meshing surface 2 order 2
Info    : [ 90%] Meshing surface 3 order 2
Info    : [100%] Meshing surface 4 order 2
Info    : Surface mesh: worst distortion = 0.999724 (0 elements in ]0, 0.2]); worst gamma = 0.860592
Info    : Done meshing order 2 (Wall 0.15568s, CPU 0.154245s)

Naghdi shell kinematics via TDC¶

A material point in the shell is parametrised by a midsurface coordinate $x_0 \in \Omega$ and a through-thickness coordinate $\xi \in [-h/2, h/2]$ . Its placement in the reference and deformed configurations is

b_0(x_0, \xi) = x_0 + \xi\, d_0(x_0), \qquad b(x_0, \xi) = x_0 + u(x_0) + \xi\, d(x_0),

(1)

where $u$ is the midsurface displacement, $d_0 = n_0$ is the reference unit normal (the director) and $d$ is the deformed director. In the Naghdi model $d$ need not stay normal to the deformed midsurface: it is parametrised by three Euler-angle rotations stored in a vector field $r$ ,

d = R_x(r_0)\, R_y(r_1)\, R_z(r_2)\, d_0.

(2)

The third rotation component rotates the director about itself and carries no physical meaning. We suppress it by introducing a scalar Lagrange multiplier $c$ that enforces $r \cdot n_0 = 0$ weakly. The primary unknowns are therefore $u, r \in [H^1(\Omega)]^3$ and $c \in L^2(\Omega)$ , all discretised here with continuous Lagrange elements of degree $p = 2$ .

metadata = {"quadrature_degree": p * p}
metadata_prsi = {"quadrature_degree": p * (p - 1)}
dx = ufl.Measure("dx", domain=mesh, subdomain_data=mesh_data.cell_tags, metadata=metadata)
ds = ufl.Measure("ds", domain=mesh, subdomain_data=mesh_data.facet_tags, metadata=metadata)
dS = ufl.Measure("dS", domain=mesh, subdomain_data=mesh_data.facet_tags, metadata=metadata)
dP = ufl.Measure("dP", domain=mesh, subdomain_data=mesh_data.ridge_tags, metadata=metadata)

Ue = basix.ufl.element("Lagrange", mesh.basix_cell(), degree=p, shape=(gdim,))
Re = basix.ufl.element("Lagrange", mesh.basix_cell(), degree=p, shape=(gdim,))
Ce = basix.ufl.element("Lagrange", mesh.basix_cell(), degree=p, shape=())

Uf = dolfinx.fem.functionspace(mesh, Ue)
Rf = dolfinx.fem.functionspace(mesh, Re)
Cf = dolfinx.fem.functionspace(mesh, Ce)

u = dolfinx.fem.Function(Uf, name="u")
r = dolfinx.fem.Function(Rf, name="r")
c = dolfinx.fem.Function(Cf, name="c")
u_ = dolfinx.fem.Function(Uf, name="u_")
r_ = dolfinx.fem.Function(Rf, name="r_")

δm = ufl.TestFunctions(ufl.MixedFunctionSpace(Uf, Rf, Cf))
δu, δr, δc = δm

m = [u, r, c]

# Saint Venant-Kirchhoff parameters (plane stress)
t0 = dolfinx.fem.Constant(mesh, h)
E_val, nu = 3102.75, 0.3
matλ = dolfinx.fem.Constant(mesh, E_val * nu / (1 - nu**2))
matμ = dolfinx.fem.Constant(mesh, E_val / (2 * (1 + nu)))

# Hinged BCs on the straight side edges (θ=±θ): zero displacement, rotations free.
# The curved axial-end edges (front/rear) are unconstrained.
sides_dofs_Uf = dolfiny.mesh.locate_dofs_topological(
    Uf, mesh_data.facet_tags, mesh_data.physical_groups["sides"].tag
)

bcs = [dolfinx.fem.dirichletbc(u_, sides_dofs_Uf)]

Strain and stress decomposition¶

The configuration gradient $J = \nabla b$ and its reference counterpart $J_0 = \nabla b_0$ combine into the Green-Lagrange strain

E = \tfrac{1}{2}\,(J^\top J - J_0^\top J_0).

(3)

Because $J$ depends affinely on the through-thickness coordinate $\xi$ , the strain admits a clean separation into membrane, bending and shear parts via the tangent plane projector $P = I - n_0 \otimes n_0$ :

E_m = P\, E\big|_{\xi=0}\, P, \quad E_b = P\, \partial_\xi E\big|_{\xi=0}\, P, \quad E_s = E\big|_{\xi=0} - E_m.

(4)

The same projection applied to the second Piola-Kirchhoff stress $S = 2\mu E + \lambda\, \text{tr}(E)\, I$ yields the corresponding stress measures $S_m$ , $S_b$ , $S_s$ . Integration through the thickness gives the customary stress resultants

N = h\, S_m, \qquad Q = h\, S_s, \qquad M = \tfrac{h^3}{12}\, S_b.

(5)

I = ufl.Identity(mesh.geometry.dim)  # noqa: E741

We = basix.ufl.element("DG", mesh.basix_cell(), degree=q, shape=(gdim,))
W = dolfinx.fem.functionspace(mesh, We)
n0 = dolfinx.fem.Function(W, name="n0")
dolfiny.interpolation.interpolate(ufl.CellNormal(mesh), n0)

Sb_norm = dolfinx.fem.Function(Cf, name="Sb_norm")

P = I - ufl.outer(n0, n0)
x0 = ufl.SpatialCoordinate(mesh)

Ξe = basix.ufl.element("DG", mesh.basix_cell(), degree=q, shape=())
Ξ = dolfinx.fem.functionspace(mesh, Ξe)
ξ = dolfinx.fem.Function(Ξ, name="ξ")

# Reference and deformed placements
d0 = n0
b0 = x0 + ξ * d0

R0 = ufl.as_matrix(
    [[1, 0, 0], [0, ufl.cos(r[0]), -ufl.sin(r[0])], [0, ufl.sin(r[0]), ufl.cos(r[0])]]
)
R1 = ufl.as_matrix(
    [[ufl.cos(r[1]), 0, ufl.sin(r[1])], [0, 1, 0], [-ufl.sin(r[1]), 0, ufl.cos(r[1])]]
)
R2 = ufl.as_matrix(
    [[ufl.cos(r[2]), -ufl.sin(r[2]), 0], [ufl.sin(r[2]), ufl.cos(r[2]), 0], [0, 0, 1]]
)

d = (R0 * R1 * R2) * d0
b = x0 + u + ξ * d

# Configuration gradients (the explicit subtraction of d0⊗d0 makes J non-degenerate
# normal-to-surface so that derivatives in ξ are well defined).
J0 = ufl.grad(b0) - ufl.outer(d0, d0)
J0 = ufl.algorithms.apply_algebra_lowering.apply_algebra_lowering(J0)
J0 = ufl.algorithms.apply_derivatives.apply_derivatives(J0)
J0 = ufl.replace(J0, {ufl.grad(ξ): d0})

J = ufl.grad(b) - ufl.outer(d0, d0)
J = ufl.algorithms.apply_algebra_lowering.apply_algebra_lowering(J)
J = ufl.algorithms.apply_derivatives.apply_derivatives(J)
J = ufl.replace(J, {ufl.grad(ξ): d0})

# Green-Lagrange strain and its membrane/bending/shear projections
E = (J.T * J - J0.T * J0) / 2

Em = P * ufl.replace(E, {ξ: 0.0}) * P

Eb = ufl.diff(E, ξ)
Eb = ufl.algorithms.apply_algebra_lowering.apply_algebra_lowering(Eb)
Eb = ufl.algorithms.apply_derivatives.apply_derivatives(Eb)
Eb = P * ufl.replace(Eb, {ξ: 0.0}) * P

Es = ufl.replace(E, {ξ: 0.0}) - P * ufl.replace(E, {ξ: 0.0}) * P

δEm = ufl.derivative(Em, m, δm)
δEs = ufl.derivative(Es, m, δm)
δEb = ufl.derivative(Eb, m, δm)

# Saint Venant-Kirchhoff PK2 stress with the same membrane/bending/shear split
S = 2 * matμ * E + matλ * ufl.tr(E) * I

Sm = P * ufl.replace(S, {ξ: 0.0}) * P

Sb = ufl.diff(S, ξ)
Sb = ufl.algorithms.apply_algebra_lowering.apply_algebra_lowering(Sb)
Sb = ufl.algorithms.apply_derivatives.apply_derivatives(Sb)
Sb = P * ufl.replace(Sb, {ξ: 0.0}) * P

Ss = ufl.replace(S, {ξ: 0.0}) - P * ufl.replace(S, {ξ: 0.0}) * P

# Through-thickness integration → stress resultants
Sm *= t0
Ss *= t0
Sb *= t0**3 / 12

Variational form and locking treatment¶

Equilibrium is expressed as the principle of virtual work

\int_\Omega \big( \delta E_m : N + \delta E_s : Q + \delta E_b : M \big) \,\text{d}x \;=\; \lambda\, p_z\, \delta u_z \big|_{\text{topmid}},

(6)

where the right-hand side is the virtual work of a unit downward concentrated load $-e_z$ applied at the co-dimension-2 point topmid (crown centre), scaled by $\lambda$ . Since the reference load magnitude is 1, $\lambda$ is the applied force. Low-order shell elements are prone to membrane and transverse-shear locking. We address this with partial selective reduced integration Arnold & Brezzi, 1997, blending fully and reduced-integrated contributions via the element-wise weight $\alpha = h^2 / |J|$ , where $|J|$ is the cell Jacobian determinant. The drilling rotation $r \cdot n_0$ is also constrained weakly by the Lagrange multiplier $c$ , giving the augmented form

\delta c \,(r \cdot n_0)\, \mu + c \,(\delta r \cdot n_0)\, \mu.

(7)

P_max = 1.0  # benchmark reference load
λ = dolfinx.fem.Constant(mesh, 1.0)

# Element-wise weight for partial selective reduced integration (Arnold/Brezzi 1997)
Ae = basix.ufl.element("DG", mesh.basix_cell(), degree=0)
A = dolfinx.fem.functionspace(mesh, Ae)
α = dolfinx.fem.Function(A)
dolfiny.interpolation.interpolate(t0**2 / ufl.JacobianDeterminant(mesh), α)

f = (
    -ufl.inner(δEm, Sm) * α * dx
    - ufl.inner(δEm, Sm) * (1 - α) * dx(metadata=metadata_prsi)
    - ufl.inner(δEs, Ss) * α * dx
    - ufl.inner(δEs, Ss) * (1 - α) * dx(metadata=metadata_prsi)
    - ufl.inner(δEb, Sb) * dx
    + δc * ufl.inner(r, n0) * matμ * dx
    + c * ufl.inner(δr, n0) * matμ * dx
    + δc * dolfinx.fem.Constant(mesh, 0.0) * c * dx
    + λ
    * ufl.inner(δu, ufl.as_vector((0.0, 0.0, -P_max)))
    * dP(mesh_data.physical_groups["topmid"].tag)
)

F = ufl.extract_blocks(f)

Reference configuration¶

Before launching the continuation we visualise the un-deformed midsurface as a sanity check on the mesh and the inferred normal field $n_0$ . Bending stress resultants are zero in this state.

def plot_roof_pyvista(u, s, png, comm=MPI.COMM_WORLD):
    if comm.size > 1:
        return

    grid = pyvista.UnstructuredGrid(*dolfinx.plot.vtk_mesh(u.function_space))

    plotter = pyvista.Plotter(off_screen=True, theme=dolfiny.pyvista.theme)
    plotter.theme.font.fmt = "%1.3f"
    plotter.add_axes()
    plotter.enable_parallel_projection()

    grid.point_data["u"] = u.x.array.reshape(-1, 3)
    grid.point_data["stress"] = s.x.array

    grid_warped = grid.warp_by_vector("u", factor=1.0)

    levels = 5 if not grid.get_cell(0).is_linear else 0

    surf = plotter.add_mesh(
        grid_warped.extract_surface(nonlinear_subdivision=levels),
        scalars="stress",
        scalar_bar_args={"title": "Bending stress resultant [-]"},
        n_colors=10,
        color="white",
    )
    surf.mapper.scalar_range = [0.0, 300]

    plotter.add_mesh(
        grid_warped.separate_cells()
        .extract_surface(nonlinear_subdivision=levels)
        .extract_feature_edges(),
        style="wireframe",
        color="black",
        line_width=dolfiny.pyvista.pixels // 1000,
        render_lines_as_tubes=True,
    )

    plotter.view_xz()
    plotter.camera.azimuth += 45
    plotter.camera.elevation += 30
    # Frame on the reference geometry so the camera does not zoom with the deformation.
    plotter.camera.focal_point = grid.center
    plotter.camera.parallel_scale = 0.8 * max(
        grid.bounds[1] - grid.bounds[0],
        grid.bounds[3] - grid.bounds[2],
        grid.bounds[5] - grid.bounds[4],
    )
    plotter.show_axes()
    plotter.screenshot(png)
    plotter.close()
    plotter.deep_clean()


plot_roof_pyvista(u, Sb_norm, f"{name}_initial.png")

Figure 2:Reference configuration of the cylindrical roof panel.

Arc-length continuation¶

Around the snap-through limit point the load-displacement tangent becomes vertical: at fixed $\lambda$ no nearby equilibrium exists, and Newton iteration on the displacement alone diverges. The Crisfield method Crisfield, 1981 addresses this by promoting $\lambda$ to an unknown and appending a spherical arc-length constraint of prescribed step $\Delta s$ ,

\langle \Delta m,\, \Delta m \rangle + \psi^2\, \Delta\lambda^2 = \Delta s^2,

(8)

where $\langle \cdot,\, \cdot \rangle$ is an inner product on the discrete state $m = (u, r, c)$ and $\psi$ is a scaling factor. The augmented system is solved at every step by dolfiny’s Crisfield class, which wraps the underlying SNES nonlinear solver. We use a sparse direct LU (MUMPS) for the linearised systems. Step size control is adaptive: if a step fails — either because SNES does not converge or the arc-length quadratic has no real root — the step is retried with $\Delta s$ halved and a zero displacement predictor; on success the step size doubles back towards the prescribed maximum $\Delta s_{\max}$ .

with dolfiny.io.XDMFFile(comm, f"{name}.xdmf", "w") as ofile:
    if q <= 2:
        ofile.write_mesh_meshtags(mesh)

opts = PETSc.Options("roof")  # type: ignore[attr-defined]
opts["snes_type"] = "newtonls"
opts["snes_linesearch_type"] = "basic"
opts["snes_atol"] = 1.0e-14
opts["snes_rtol"] = 1.0e-08
# opts["snes_stol"] = 1.0e-08
opts["snes_max_it"] = 40
opts["ksp_type"] = "preonly"
opts["pc_type"] = "lu"
opts["pc_factor_mat_solver_type"] = "mumps"

u_step: list[np.ndarray] = []
λ_step: list[float] = []


def monitor(context=None):
    if comm.size > 1:
        return

    idx = dolfiny.mesh.locate_dofs_topological(
        Uf, mesh_data.ridge_tags, mesh_data.physical_groups["topmid"].tag
    )
    idx = dolfiny.function.unroll_dofs(idx, Uf.dofmap.bs)

    u_step.append(u.x.array[idx])
    λ_step.append(context.λ.value.item())


def block_inner(a1, a2):
    b1, b2 = [], []
    for mi in m:
        b1.append(dolfinx.fem.Function(mi.function_space, name=mi.name))
        b2.append(dolfinx.fem.Function(mi.function_space, name=mi.name))

    dolfinx.fem.petsc.assign(a1, b1)
    dolfinx.fem.petsc.assign(a2, b2)
    inner = 0.0
    for b1i, b2i in zip(b1, b2):
        inner += dolfiny.expression.assemble(ufl.inner(b1i, b2i), ufl.dx(mesh))
    return inner


def plot_load_displacement(u_step, λ_step, png):
    if comm.size > 1:
        return

    fig, ax1 = plt.subplots(figsize=(8, 6), dpi=400)
    ax1.set_title(
        f"Cylindrical roof $h={h}$, elems=[{(N - 1) * 2}, {(N - 1) * 2}]",
        fontsize=12,
    )
    ax1.set_xlabel("crown displacement $-u_z$ [-]", fontsize=12)
    ax1.set_ylabel("load $P$ [-]", fontsize=12)
    ax1.grid(linewidth=0.25)
    fig.tight_layout()

    sol_tdc = np.column_stack([λ_step, u_step])
    ax1.plot(
        -sol_tdc[:, 3],
        sol_tdc[:, 0],
        ls="-",
        lw=1.0,
        color="k",
        label="TDC (Naghdi, this work)",
    )

    # Overlay Sze (2004) Table 9d reference data (S4R, 24×24 mesh, P_max=3000).
    sze = np.loadtxt(f"{name}_sze2004.csv", delimiter=",", skiprows=1)
    ax1.plot(
        sze[:, 1],
        sze[:, 0],
        marker="o",
        ls="",
        color="tab:blue",
        label="Sze (2004), S4R 24×24",
    )

    ax1.legend()
    fig.tight_layout()
    fig.savefig(png)
    plt.close(fig)


problem = dolfiny.snesproblem.SNESProblem(F, m, bcs, prefix="roof")

continuation = dolfiny.continuation.Crisfield(problem, λ, inner=block_inner)
continuation.initialise(ds=0.1, λ=0.0, psi=1.0)
monitor(continuation)

for k in range(64):
    dolfiny.utils.pprint(f"\n*** Continuation step {k:d}")
    continuation.solve_step(ds=200)

    monitor(continuation)

    with dolfiny.io.XDMFFile(comm, f"{name}.xdmf", "a") as ofile:
        if q <= 2:
            ofile.write_function(u, k)
            dolfiny.projection.project(ufl.inner(Sb, Sb) ** 0.5, Sb_norm)
            ofile.write_function(Sb_norm, k)


*** Continuation step 0

# SNES iteration  0

# sub  0 [  6k] |x|=0.000e+00 |dx|=0.000e+00 |r|=1.000e-01 (u)

# sub  1 [  6k] |x|=0.000e+00 |dx|=0.000e+00 |r|=0.000e+00 (r)

# sub  2 [  2k] |x|=0.000e+00 |dx|=0.000e+00 |r|=0.000e+00 (c)

# all           |x|=0.000e+00 |dx|=0.000e+00 |r|=1.000e-01

# arc           |x|=1.000e-01 |dx|=0.000e+00 |r|=4.000e+04 (λ)

# SNES iteration  0, KSP iteration   0       |r|=1.000e-01

# SNES iteration  0, KSP iteration   1       |r|=7.882e-13

# SNES iteration  1

# sub  0 [  6k] |x|=1.328e-02 |dx|=1.328e-02 |r|=2.057e-05 (u)

# sub  1 [  6k] |x|=1.495e-04 |dx|=1.495e-04 |r|=2.061e-05 (r)

# sub  2 [  2k] |x|=8.824e-09 |dx|=8.824e-09 |r|=7.258e-13 (c)

# all           |x|=1.328e-02 |dx|=1.328e-02 |r|=2.912e-05

# SNES iteration  1, KSP iteration   0       |r|=1.000e+00

# SNES iteration  1, KSP iteration   1       |r|=5.177e-12

# arc           |x|=1.318e+02 |dx|=1.317e+02 |r|=1.237e-10 (λ)

# SNES iteration  1, KSP iteration   0       |r|=4.860e+01

# SNES iteration  1, KSP iteration   1       |r|=1.613e-10

# SNES iteration  2

# sub  0 [  6k] |x|=1.985e+01 |dx|=2.346e+00 |r|=6.305e-01 (u)

# sub  1 [  6k] |x|=2.228e-01 |dx|=2.619e-02 |r|=6.457e-01 (r)

# sub  2 [  2k] |x|=1.186e-05 |dx|=5.090e-07 |r|=1.498e-10 (c)

# all           |x|=1.985e+01 |dx|=2.346e+00 |r|=9.025e-01

# SNES iteration  2, KSP iteration   0       |r|=1.000e+00

# SNES iteration  2, KSP iteration   1       |r|=8.629e-12

# arc           |x|=1.178e+02 |dx|=1.407e+01 |r|=2.110e-10 (λ)

# SNES iteration  2, KSP iteration   0       |r|=7.173e-02

# SNES iteration  2, KSP iteration   1       |r|=2.365e-13

# SNES iteration  3

# sub  0 [  6k] |x|=1.751e+01 |dx|=9.574e-04 |r|=4.188e-07 (u)

# sub  1 [  6k] |x|=1.967e-01 |dx|=1.477e-05 |r|=2.371e-07 (r)

# sub  2 [  2k] |x|=1.057e-05 |dx|=8.352e-10 |r|=2.425e-13 (c)

# all           |x|=1.751e+01 |dx|=9.576e-04 |r|=4.813e-07

# SNES iteration  3, KSP iteration   0       |r|=1.000e+00

# SNES iteration  3, KSP iteration   1       |r|=6.757e-12

# arc           |x|=1.178e+02 |dx|=2.734e-03 |r|=3.274e-10 (λ)

# SNES iteration  3, KSP iteration   0       |r|=6.142e-07

# SNES iteration  3, KSP iteration   1       |r|=7.353e-19

# SNES iteration  4 success = CONVERGED_SNORM_RELATIVE

# sub  0 [  6k] |x|=1.751e+01 |dx|=5.714e-09 |r|=2.111e-09 (u)

# sub  1 [  6k] |x|=1.967e-01 |dx|=7.963e-11 |r|=5.736e-10 (r)

# sub  2 [  2k] |x|=1.057e-05 |dx|=4.415e-15 |r|=5.769e-14 (c)

# all           |x|=1.751e+01 |dx|=5.715e-09 |r|=2.187e-09

# arc           |x|=1.178e+02 |dx|=1.813e-08 |r|=4.075e-10 (λ)


*** Continuation step 1

# SNES iteration  0

# sub  0 [  6k] |x|=3.502e+01 |dx|=5.714e-09 |r|=7.028e+01 (u)

# sub  1 [  6k] |x|=3.933e-01 |dx|=7.963e-11 |r|=5.585e+01 (r)

# sub  2 [  2k] |x|=2.114e-05 |dx|=4.415e-15 |r|=1.315e-13 (c)

# all           |x|=3.502e+01 |dx|=5.715e-09 |r|=8.977e+01

# arc           |x|=2.355e+02 |dx|=0.000e+00 |r|=4.075e-10 (λ)

# SNES iteration  0, KSP iteration   0       |r|=8.977e+01

# SNES iteration  0, KSP iteration   1       |r|=3.148e-10

# SNES iteration  1

# sub  0 [  6k] |x|=4.004e+01 |dx|=5.053e+00 |r|=2.496e+00 (u)

# sub  1 [  6k] |x|=4.455e-01 |dx|=5.334e-02 |r|=2.321e+00 (r)

# sub  2 [  2k] |x|=2.162e-05 |dx|=1.065e-06 |r|=2.907e-10 (c)

# all           |x|=4.004e+01 |dx|=5.053e+00 |r|=3.409e+00

# SNES iteration  1, KSP iteration   0       |r|=1.000e+00

# SNES iteration  1, KSP iteration   1       |r|=9.445e-12

# arc           |x|=2.112e+02 |dx|=2.432e+01 |r|=2.037e-10 (λ)

# SNES iteration  1, KSP iteration   0       |r|=3.320e-01

# SNES iteration  1, KSP iteration   1       |r|=4.761e-13

# SNES iteration  2

# sub  0 [  6k] |x|=3.501e+01 |dx|=5.228e-03 |r|=1.681e-05 (u)

# sub  1 [  6k] |x|=3.906e-01 |dx|=9.087e-05 |r|=8.639e-06 (r)

# sub  2 [  2k] |x|=1.930e-05 |dx|=3.542e-09 |r|=4.713e-13 (c)

# all           |x|=3.501e+01 |dx|=5.228e-03 |r|=1.890e-05

# SNES iteration  2, KSP iteration   0       |r|=1.000e+00

# SNES iteration  2, KSP iteration   1       |r|=1.106e-11

# arc           |x|=2.112e+02 |dx|=1.278e-02 |r|=2.183e-11 (λ)

# SNES iteration  2, KSP iteration   0       |r|=1.333e-05

# SNES iteration  2, KSP iteration   1       |r|=4.983e-17

# SNES iteration  3 success = CONVERGED_FNORM_RELATIVE

# sub  0 [  6k] |x|=3.501e+01 |dx|=1.600e-07 |r|=2.138e-09 (u)

# sub  1 [  6k] |x|=3.906e-01 |dx|=2.205e-09 |r|=6.766e-10 (r)

# sub  2 [  2k] |x|=1.930e-05 |dx|=1.119e-13 |r|=1.114e-13 (c)

# all           |x|=3.501e+01 |dx|=1.600e-07 |r|=2.242e-09

# arc           |x|=2.112e+02 |dx|=1.188e-07 |r|=6.548e-11 (λ)


*** Continuation step 2

# SNES iteration  0

# sub  0 [  6k] |x|=5.251e+01 |dx|=1.600e-07 |r|=6.648e+01 (u)

# sub  1 [  6k] |x|=5.847e-01 |dx|=2.205e-09 |r|=4.296e+01 (r)

# sub  2 [  2k] |x|=2.804e-05 |dx|=1.119e-13 |r|=2.530e-13 (c)

# all           |x|=5.252e+01 |dx|=1.600e-07 |r|=7.915e+01

# arc           |x|=3.047e+02 |dx|=0.000e+00 |r|=6.548e-11 (λ)

# SNES iteration  0, KSP iteration   0       |r|=7.915e+01

# SNES iteration  0, KSP iteration   1       |r|=6.269e-10

# SNES iteration  1

# sub  0 [  6k] |x|=5.727e+01 |dx|=4.841e+00 |r|=1.820e+00 (u)

# sub  1 [  6k] |x|=6.298e-01 |dx|=4.694e-02 |r|=1.755e+00 (r)

# sub  2 [  2k] |x|=2.851e-05 |dx|=1.040e-06 |r|=6.158e-10 (c)

# all           |x|=5.728e+01 |dx|=4.841e+00 |r|=2.529e+00

# SNES iteration  1, KSP iteration   0       |r|=1.000e+00

# SNES iteration  1, KSP iteration   1       |r|=9.667e-12

# arc           |x|=2.859e+02 |dx|=1.876e+01 |r|=1.673e-10 (λ)

# SNES iteration  1, KSP iteration   0       |r|=3.965e-01

# SNES iteration  1, KSP iteration   1       |r|=1.246e-12

# SNES iteration  2

# sub  0 [  6k] |x|=5.248e+01 |dx|=1.271e-02 |r|=5.063e-05 (u)

# sub  1 [  6k] |x|=5.792e-01 |dx|=1.853e-04 |r|=2.226e-05 (r)

# sub  2 [  2k] |x|=2.663e-05 |dx|=4.627e-09 |r|=1.222e-12 (c)

# all           |x|=5.248e+01 |dx|=1.271e-02 |r|=5.530e-05

# SNES iteration  2, KSP iteration   0       |r|=1.000e+00

# SNES iteration  2, KSP iteration   1       |r|=2.361e-11

# arc           |x|=2.860e+02 |dx|=4.254e-02 |r|=1.819e-10 (λ)

# SNES iteration  2, KSP iteration   0       |r|=1.957e-05

# SNES iteration  2, KSP iteration   1       |r|=2.461e-16

# SNES iteration  3 success = CONVERGED_FNORM_RELATIVE

# sub  0 [  6k] |x|=5.249e+01 |dx|=3.411e-07 |r|=2.149e-09 (u)

# sub  1 [  6k] |x|=5.793e-01 |dx|=5.159e-09 |r|=7.711e-10 (r)

# sub  2 [  2k] |x|=2.663e-05 |dx|=1.770e-13 |r|=1.629e-13 (c)

# all           |x|=5.249e+01 |dx|=3.411e-07 |r|=2.283e-09

# arc           |x|=2.860e+02 |dx|=4.087e-07 |r|=2.401e-10 (λ)


*** Continuation step 3

# SNES iteration  0

# sub  0 [  6k] |x|=6.998e+01 |dx|=3.411e-07 |r|=6.087e+01 (u)

# sub  1 [  6k] |x|=7.681e-01 |dx|=5.159e-09 |r|=3.309e+01 (r)

# sub  2 [  2k] |x|=3.398e-05 |dx|=1.770e-13 |r|=3.783e-13 (c)

# all           |x|=6.998e+01 |dx|=3.411e-07 |r|=6.928e+01

# arc           |x|=3.607e+02 |dx|=0.000e+00 |r|=2.401e-10 (λ)

# SNES iteration  0, KSP iteration   0       |r|=6.928e+01

# SNES iteration  0, KSP iteration   1       |r|=7.539e-10

# SNES iteration  1

# sub  0 [  6k] |x|=7.424e+01 |dx|=4.412e+00 |r|=1.183e+00 (u)

# sub  1 [  6k] |x|=8.043e-01 |dx|=3.897e-02 |r|=1.244e+00 (r)

# sub  2 [  2k] |x|=3.444e-05 |dx|=1.019e-06 |r|=7.458e-10 (c)

# all           |x|=7.425e+01 |dx|=4.412e+00 |r|=1.717e+00

# SNES iteration  1, KSP iteration   0       |r|=1.000e+00

# SNES iteration  1, KSP iteration   1       |r|=1.911e-11

# arc           |x|=3.467e+02 |dx|=1.399e+01 |r|=3.420e-10 (λ)

# SNES iteration  1, KSP iteration   0       |r|=4.375e-01

# SNES iteration  1, KSP iteration   1       |r|=7.582e-13

# SNES iteration  2

# sub  0 [  6k] |x|=6.992e+01 |dx|=2.170e-02 |r|=1.017e-04 (u)

# sub  1 [  6k] |x|=7.605e-01 |dx|=2.904e-04 |r|=4.063e-05 (r)

# sub  2 [  2k] |x|=3.294e-05 |dx|=6.028e-09 |r|=5.745e-13 (c)

# all           |x|=6.992e+01 |dx|=2.171e-02 |r|=1.096e-04

# SNES iteration  2, KSP iteration   0       |r|=1.000e+00

# SNES iteration  2, KSP iteration   1       |r|=5.440e-11

# arc           |x|=3.468e+02 |dx|=6.395e-02 |r|=2.037e-10 (λ)

# SNES iteration  2, KSP iteration   0       |r|=2.628e-05

# SNES iteration  2, KSP iteration   1       |r|=6.936e-17

# SNES iteration  3 success = CONVERGED_FNORM_RELATIVE

# sub  0 [  6k] |x|=6.994e+01 |dx|=6.699e-07 |r|=2.136e-09 (u)

# sub  1 [  6k] |x|=7.607e-01 |dx|=1.045e-08 |r|=9.076e-10 (r)

# sub  2 [  2k] |x|=3.295e-05 |dx|=2.537e-13 |r|=2.150e-13 (c)

# all           |x|=6.994e+01 |dx|=6.700e-07 |r|=2.321e-09

# arc           |x|=3.468e+02 |dx|=1.301e-06 |r|=6.476e-10 (λ)


*** Continuation step 4

# SNES iteration  0

# sub  0 [  6k] |x|=8.739e+01 |dx|=6.699e-07 |r|=5.438e+01 (u)

# sub  1 [  6k] |x|=9.424e-01 |dx|=1.045e-08 |r|=2.626e+01 (r)

# sub  2 [  2k] |x|=3.930e-05 |dx|=2.537e-13 |r|=4.447e-13 (c)

# all           |x|=8.740e+01 |dx|=6.700e-07 |r|=6.038e+01

# arc           |x|=4.076e+02 |dx|=0.000e+00 |r|=6.476e-10 (λ)

# SNES iteration  0, KSP iteration   0       |r|=6.038e+01

# SNES iteration  0, KSP iteration   1       |r|=3.707e-10

# SNES iteration  1

# sub  0 [  6k] |x|=9.110e+01 |dx|=3.933e+00 |r|=7.347e-01 (u)

# sub  1 [  6k] |x|=9.700e-01 |dx|=3.165e-02 |r|=8.655e-01 (r)

# sub  2 [  2k] |x|=3.975e-05 |dx|=9.904e-07 |r|=3.578e-10 (c)

# all           |x|=9.111e+01 |dx|=3.933e+00 |r|=1.135e+00

# SNES iteration  1, KSP iteration   0       |r|=1.000e+00

# SNES iteration  1, KSP iteration   1       |r|=5.637e-11

# arc           |x|=3.973e+02 |dx|=1.037e+01 |r|=3.929e-10 (λ)

# SNES iteration  1, KSP iteration   0       |r|=4.443e-01

# SNES iteration  1, KSP iteration   1       |r|=2.051e-12

# SNES iteration  2

# sub  0 [  6k] |x|=8.730e+01 |dx|=2.944e-02 |r|=1.468e-04 (u)

# sub  1 [  6k] |x|=9.333e-01 |dx|=3.706e-04 |r|=5.732e-05 (r)

# sub  2 [  2k] |x|=3.856e-05 |dx|=7.295e-09 |r|=1.927e-12 (c)

# all           |x|=8.731e+01 |dx|=2.944e-02 |r|=1.576e-04

# SNES iteration  2, KSP iteration   0       |r|=1.000e+00

# SNES iteration  2, KSP iteration   1       |r|=4.505e-11

# arc           |x|=3.973e+02 |dx|=7.442e-02 |r|=1.965e-10 (λ)

# SNES iteration  2, KSP iteration   0       |r|=3.064e-05

# SNES iteration  2, KSP iteration   1       |r|=5.365e-17

# SNES iteration  3 success = CONVERGED_FNORM_RELATIVE

# sub  0 [  6k] |x|=8.733e+01 |dx|=1.070e-06 |r|=2.159e-09 (u)

# sub  1 [  6k] |x|=9.336e-01 |dx|=1.613e-08 |r|=1.063e-09 (r)

# sub  2 [  2k] |x|=3.857e-05 |dx|=3.244e-13 |r|=2.638e-13 (c)

# all           |x|=8.733e+01 |dx|=1.070e-06 |r|=2.406e-09

# arc           |x|=3.973e+02 |dx|=2.122e-06 |r|=6.257e-10 (λ)


*** Continuation step 5

# SNES iteration  0

# sub  0 [  6k] |x|=1.047e+02 |dx|=1.070e-06 |r|=4.779e+01 (u)

# sub  1 [  6k] |x|=1.107e+00 |dx|=1.613e-08 |r|=2.172e+01 (r)

# sub  2 [  2k] |x|=4.423e-05 |dx|=3.244e-13 |r|=5.829e-13 (c)

# all           |x|=1.048e+02 |dx|=1.070e-06 |r|=5.250e+01

# arc           |x|=4.479e+02 |dx|=0.000e+00 |r|=6.257e-10 (λ)

# SNES iteration  0, KSP iteration   0       |r|=5.250e+01

# SNES iteration  0, KSP iteration   1       |r|=2.525e-10

# SNES iteration  1

# sub  0 [  6k] |x|=1.080e+02 |dx|=3.550e+00 |r|=4.752e-01 (u)

# sub  1 [  6k] |x|=1.128e+00 |dx|=2.639e-02 |r|=6.301e-01 (r)

# sub  2 [  2k] |x|=4.467e-05 |dx|=9.571e-07 |r|=2.365e-10 (c)

# all           |x|=1.080e+02 |dx|=3.550e+00 |r|=7.892e-01

# SNES iteration  1, KSP iteration   0       |r|=1.000e+00

# SNES iteration  1, KSP iteration   1       |r|=6.090e-11

# arc           |x|=4.400e+02 |dx|=7.887e+00 |r|=5.020e-10 (λ)

# SNES iteration  1, KSP iteration   0       |r|=4.213e-01

# SNES iteration  1, KSP iteration   1       |r|=1.424e-12

# SNES iteration  2

# sub  0 [  6k] |x|=1.046e+02 |dx|=3.459e-02 |r|=1.686e-04 (u)

# sub  1 [  6k] |x|=1.097e+00 |dx|=4.128e-04 |r|=6.643e-05 (r)

# sub  2 [  2k] |x|=4.369e-05 |dx|=8.304e-09 |r|=1.177e-12 (c)

# all           |x|=1.046e+02 |dx|=3.460e-02 |r|=1.812e-04

# SNES iteration  2, KSP iteration   0       |r|=1.000e+00

# SNES iteration  2, KSP iteration   1       |r|=9.178e-11

# arc           |x|=4.401e+02 |dx|=7.521e-02 |r|=5.239e-10 (λ)

# SNES iteration  2, KSP iteration   0       |r|=3.206e-05

# SNES iteration  2, KSP iteration   1       |r|=1.645e-16

# SNES iteration  3 success = CONVERGED_FNORM_RELATIVE

# sub  0 [  6k] |x|=1.047e+02 |dx|=1.383e-06 |r|=2.098e-09 (u)

# sub  1 [  6k] |x|=1.097e+00 |dx|=1.994e-08 |r|=1.193e-09 (r)

# sub  2 [  2k] |x|=4.370e-05 |dx|=3.693e-13 |r|=3.177e-13 (c)

# all           |x|=1.047e+02 |dx|=1.383e-06 |r|=2.413e-09

# arc           |x|=4.401e+02 |dx|=2.490e-06 |r|=1.601e-10 (λ)


*** Continuation step 6

# SNES iteration  0

# sub  0 [  6k] |x|=1.220e+02 |dx|=1.383e-06 |r|=4.163e+01 (u)

# sub  1 [  6k] |x|=1.261e+00 |dx|=1.994e-08 |r|=1.855e+01 (r)

# sub  2 [  2k] |x|=4.889e-05 |dx|=3.693e-13 |r|=6.582e-13 (c)

# all           |x|=1.220e+02 |dx|=1.383e-06 |r|=4.557e+01

# arc           |x|=4.828e+02 |dx|=0.000e+00 |r|=1.601e-10 (λ)

# SNES iteration  0, KSP iteration   0       |r|=4.557e+01

# SNES iteration  0, KSP iteration   1       |r|=3.128e-10

# SNES iteration  1

# sub  0 [  6k] |x|=1.250e+02 |dx|=3.358e+00 |r|=3.530e-01 (u)

# sub  1 [  6k] |x|=1.279e+00 |dx|=2.359e-02 |r|=5.144e-01 (r)

# sub  2 [  2k] |x|=4.935e-05 |dx|=9.373e-07 |r|=3.016e-10 (c)

# all           |x|=1.250e+02 |dx|=3.358e+00 |r|=6.239e-01

# SNES iteration  1, KSP iteration   0       |r|=1.000e+00

# SNES iteration  1, KSP iteration   1       |r|=1.363e-10

# arc           |x|=4.764e+02 |dx|=6.332e+00 |r|=3.129e-10 (λ)

# SNES iteration  1, KSP iteration   0       |r|=3.813e-01

# SNES iteration  1, KSP iteration   1       |r|=3.278e-12

# SNES iteration  2

# sub  0 [  6k] |x|=1.219e+02 |dx|=3.725e-02 |r|=1.688e-04 (u)

# sub  1 [  6k] |x|=1.251e+00 |dx|=4.227e-04 |r|=6.716e-05 (r)

# sub  2 [  2k] |x|=4.848e-05 |dx|=9.071e-09 |r|=3.163e-12 (c)

# all           |x|=1.219e+02 |dx|=3.725e-02 |r|=1.817e-04

# SNES iteration  2, KSP iteration   0       |r|=1.000e+00

# SNES iteration  2, KSP iteration   1       |r|=8.329e-11

# arc           |x|=4.765e+02 |dx|=6.957e-02 |r|=8.731e-11 (λ)

# SNES iteration  2, KSP iteration   0       |r|=3.305e-05

# SNES iteration  2, KSP iteration   1       |r|=6.638e-17

# SNES iteration  3 success = CONVERGED_FNORM_RELATIVE

# sub  0 [  6k] |x|=1.219e+02 |dx|=1.562e-06 |r|=2.049e-09 (u)

# sub  1 [  6k] |x|=1.251e+00 |dx|=2.185e-08 |r|=1.351e-09 (r)

# sub  2 [  2k] |x|=4.849e-05 |dx|=4.268e-13 |r|=3.438e-13 (c)

# all           |x|=1.219e+02 |dx|=1.562e-06 |r|=2.454e-09

# arc           |x|=4.765e+02 |dx|=2.416e-06 |r|=3.274e-10 (λ)


*** Continuation step 7

# SNES iteration  0

# sub  0 [  6k] |x|=1.392e+02 |dx|=1.562e-06 |r|=3.609e+01 (u)

# sub  1 [  6k] |x|=1.406e+00 |dx|=2.185e-08 |r|=1.610e+01 (r)

# sub  2 [  2k] |x|=5.335e-05 |dx|=4.268e-13 |r|=7.399e-13 (c)

# all           |x|=1.392e+02 |dx|=1.562e-06 |r|=3.952e+01

# arc           |x|=5.130e+02 |dx|=0.000e+00 |r|=3.274e-10 (λ)

# SNES iteration  0, KSP iteration   0       |r|=3.952e+01

# SNES iteration  0, KSP iteration   1       |r|=4.368e-10

# SNES iteration  1

# sub  0 [  6k] |x|=1.422e+02 |dx|=3.418e+00 |r|=3.252e-01 (u)

# sub  1 [  6k] |x|=1.422e+00 |dx|=2.322e-02 |r|=4.937e-01 (r)

# sub  2 [  2k] |x|=5.385e-05 |dx|=9.556e-07 |r|=4.291e-10 (c)

# all           |x|=1.422e+02 |dx|=3.418e+00 |r|=5.912e-01

# SNES iteration  1, KSP iteration   0       |r|=1.000e+00

# SNES iteration  1, KSP iteration   1       |r|=1.441e-10

# arc           |x|=5.075e+02 |dx|=5.467e+00 |r|=2.183e-10 (λ)

# SNES iteration  1, KSP iteration   0       |r|=3.381e-01

# SNES iteration  1, KSP iteration   1       |r|=2.620e-12

# SNES iteration  2

# sub  0 [  6k] |x|=1.390e+02 |dx|=3.823e-02 |r|=1.574e-04 (u)

# sub  1 [  6k] |x|=1.396e+00 |dx|=4.130e-04 |r|=6.292e-05 (r)

# sub  2 [  2k] |x|=5.301e-05 |dx|=9.652e-09 |r|=2.469e-12 (c)

# all           |x|=1.390e+02 |dx|=3.823e-02 |r|=1.695e-04

# SNES iteration  2, KSP iteration   0       |r|=1.000e+00

# SNES iteration  2, KSP iteration   1       |r|=4.219e-11

# arc           |x|=5.076e+02 |dx|=6.068e-02 |r|=1.673e-10 (λ)

# SNES iteration  2, KSP iteration   0       |r|=3.550e-05

# SNES iteration  2, KSP iteration   1       |r|=8.897e-17

# SNES iteration  3 success = CONVERGED_FNORM_RELATIVE

# sub  0 [  6k] |x|=1.391e+02 |dx|=1.752e-06 |r|=2.102e-09 (u)

# sub  1 [  6k] |x|=1.396e+00 |dx|=2.428e-08 |r|=1.582e-09 (r)

# sub  2 [  2k] |x|=5.302e-05 |dx|=6.009e-13 |r|=3.868e-13 (c)

# all           |x|=1.391e+02 |dx|=1.753e-06 |r|=2.631e-09

# arc           |x|=5.076e+02 |dx|=2.110e-06 |r|=2.910e-11 (λ)


*** Continuation step 8

# SNES iteration  0

# sub  0 [  6k] |x|=1.563e+02 |dx|=1.752e-06 |r|=3.122e+01 (u)

# sub  1 [  6k] |x|=1.542e+00 |dx|=2.428e-08 |r|=1.406e+01 (r)

# sub  2 [  2k] |x|=5.763e-05 |dx|=6.009e-13 |r|=7.903e-13 (c)

# all           |x|=1.563e+02 |dx|=1.753e-06 |r|=3.424e+01

# arc           |x|=5.386e+02 |dx|=0.000e+00 |r|=2.910e-11 (λ)

# SNES iteration  0, KSP iteration   0       |r|=3.424e+01

# SNES iteration  0, KSP iteration   1       |r|=8.602e-10

# SNES iteration  1

# sub  0 [  6k] |x|=1.596e+02 |dx|=3.814e+00 |r|=3.833e-01 (u)

# sub  1 [  6k] |x|=1.559e+00 |dx|=2.549e-02 |r|=5.656e-01 (r)

# sub  2 [  2k] |x|=5.822e-05 |dx|=1.043e-06 |r|=8.554e-10 (c)

# all           |x|=1.597e+02 |dx|=3.814e+00 |r|=6.833e-01

# SNES iteration  1, KSP iteration   0       |r|=1.000e+00

# SNES iteration  1, KSP iteration   1       |r|=1.216e-10

# arc           |x|=5.335e+02 |dx|=5.091e+00 |r|=3.929e-10 (λ)

# SNES iteration  1, KSP iteration   0       |r|=3.026e-01

# SNES iteration  1, KSP iteration   1       |r|=3.026e-12

# SNES iteration  2

# sub  0 [  6k] |x|=1.561e+02 |dx|=3.862e-02 |r|=1.423e-04 (u)

# sub  1 [  6k] |x|=1.531e+00 |dx|=3.963e-04 |r|=5.773e-05 (r)

# sub  2 [  2k] |x|=5.732e-05 |dx|=1.018e-08 |r|=2.898e-12 (c)

# all           |x|=1.561e+02 |dx|=3.863e-02 |r|=1.535e-04

# SNES iteration  2, KSP iteration   0       |r|=1.000e+00

# SNES iteration  2, KSP iteration   1       |r|=1.201e-10

# arc           |x|=5.336e+02 |dx|=5.069e-02 |r|=7.130e-10 (λ)

# SNES iteration  2, KSP iteration   0       |r|=3.914e-05

# SNES iteration  2, KSP iteration   1       |r|=1.509e-15

# SNES iteration  3 success = CONVERGED_FNORM_RELATIVE

# sub  0 [  6k] |x|=1.561e+02 |dx|=5.527e-06 |r|=2.100e-09 (u)

# sub  1 [  6k] |x|=1.531e+00 |dx|=7.311e-08 |r|=1.694e-09 (r)

# sub  2 [  2k] |x|=5.733e-05 |dx|=2.882e-12 |r|=4.153e-13 (c)

# all           |x|=1.561e+02 |dx|=5.528e-06 |r|=2.698e-09

# arc           |x|=5.336e+02 |dx|=1.704e-06 |r|=3.638e-11 (λ)


*** Continuation step 9

# SNES iteration  0

# sub  0 [  6k] |x|=1.733e+02 |dx|=5.527e-06 |r|=2.701e+01 (u)

# sub  1 [  6k] |x|=1.668e+00 |dx|=7.311e-08 |r|=1.234e+01 (r)

# sub  2 [  2k] |x|=6.173e-05 |dx|=2.882e-12 |r|=8.702e-13 (c)

# all           |x|=1.733e+02 |dx|=5.528e-06 |r|=2.970e+01

# arc           |x|=5.596e+02 |dx|=0.000e+00 |r|=3.638e-11 (λ)

# SNES iteration  0, KSP iteration   0       |r|=2.970e+01

# SNES iteration  0, KSP iteration   1       |r|=1.178e-09

# SNES iteration  1

# sub  0 [  6k] |x|=1.775e+02 |dx|=4.741e+00 |r|=5.878e-01 (u)

# sub  1 [  6k] |x|=1.689e+00 |dx|=3.155e-02 |r|=7.826e-01 (r)

# sub  2 [  2k] |x|=6.248e-05 |dx|=1.257e-06 |r|=1.171e-09 (c)

# all           |x|=1.775e+02 |dx|=4.741e+00 |r|=9.788e-01

# SNES iteration  1, KSP iteration   0       |r|=1.000e+00

# SNES iteration  1, KSP iteration   1       |r|=1.515e-10

# arc           |x|=5.545e+02 |dx|=5.055e+00 |r|=1.091e-10 (λ)

# SNES iteration  1, KSP iteration   0       |r|=2.849e-01

# SNES iteration  1, KSP iteration   1       |r|=7.400e-12

# SNES iteration  2

# sub  0 [  6k] |x|=1.730e+02 |dx|=3.951e-02 |r|=1.275e-04 (u)

# sub  1 [  6k] |x|=1.657e+00 |dx|=3.809e-04 |r|=5.439e-05 (r)

# sub  2 [  2k] |x|=6.142e-05 |dx|=1.072e-08 |r|=7.341e-12 (c)

# all           |x|=1.730e+02 |dx|=3.952e-02 |r|=1.386e-04

# SNES iteration  2, KSP iteration   0       |r|=1.000e+00

# SNES iteration  2, KSP iteration   1       |r|=1.587e-10

# arc           |x|=5.546e+02 |dx|=4.078e-02 |r|=1.237e-10 (λ)

# SNES iteration  2, KSP iteration   0       |r|=4.218e-05

# SNES iteration  2, KSP iteration   1       |r|=1.085e-16

# SNES iteration  3 success = CONVERGED_FNORM_RELATIVE

# sub  0 [  6k] |x|=1.731e+02 |dx|=2.066e-06 |r|=2.124e-09 (u)

# sub  1 [  6k] |x|=1.657e+00 |dx|=3.042e-08 |r|=1.902e-09 (r)

# sub  2 [  2k] |x|=6.143e-05 |dx|=9.787e-13 |r|=4.428e-13 (c)

# all           |x|=1.731e+02 |dx|=2.066e-06 |r|=2.851e-09

# arc           |x|=5.546e+02 |dx|=1.334e-06 |r|=7.276e-12 (λ)


*** Continuation step 10

# SNES iteration  0

# sub  0 [  6k] |x|=1.901e+02 |dx|=2.066e-06 |r|=2.350e+01 (u)

# sub  1 [  6k] |x|=1.785e+00 |dx|=3.042e-08 |r|=1.102e+01 (r)

# sub  2 [  2k] |x|=6.561e-05 |dx|=9.787e-13 |r|=1.015e-12 (c)

# all           |x|=1.901e+02 |dx|=2.066e-06 |r|=2.595e+01

# arc           |x|=5.755e+02 |dx|=0.000e+00 |r|=7.276e-12 (λ)

# SNES iteration  0, KSP iteration   0       |r|=2.595e+01

# SNES iteration  0, KSP iteration   1       |r|=6.492e-10

# SNES iteration  1

# sub  0 [  6k] |x|=1.962e+02 |dx|=6.779e+00 |r|=1.228e+00 (u)

# sub  1 [  6k] |x|=1.817e+00 |dx|=4.526e-02 |r|=1.383e+00 (r)

# sub  2 [  2k] |x|=6.673e-05 |dx|=1.742e-06 |r|=6.267e-10 (c)

# all           |x|=1.962e+02 |dx|=6.779e+00 |r|=1.850e+00

# SNES iteration  1, KSP iteration   0       |r|=1.000e+00

# SNES iteration  1, KSP iteration   1       |r|=4.849e-10

# arc           |x|=5.703e+02 |dx|=5.260e+00 |r|=1.965e-10 (λ)

# SNES iteration  1, KSP iteration   0       |r|=2.954e-01

# SNES iteration  1, KSP iteration   1       |r|=8.371e-12

# SNES iteration  2

# sub  0 [  6k] |x|=1.898e+02 |dx|=4.299e-02 |r|=1.186e-04 (u)

# sub  1 [  6k] |x|=1.774e+00 |dx|=3.819e-04 |r|=5.570e-05 (r)

# sub  2 [  2k] |x|=6.528e-05 |dx|=1.184e-08 |r|=8.311e-12 (c)

# all           |x|=1.899e+02 |dx|=4.299e-02 |r|=1.310e-04

# SNES iteration  2, KSP iteration   0       |r|=1.000e+00

# SNES iteration  2, KSP iteration   1       |r|=1.721e-10

# arc           |x|=5.703e+02 |dx|=3.154e-02 |r|=2.183e-11 (λ)

# SNES iteration  2, KSP iteration   0       |r|=4.455e-05

# SNES iteration  2, KSP iteration   1       |r|=4.685e-16

# SNES iteration  3 success = CONVERGED_FNORM_RELATIVE

# sub  0 [  6k] |x|=1.899e+02 |dx|=2.345e-06 |r|=2.070e-09 (u)

# sub  1 [  6k] |x|=1.774e+00 |dx|=3.641e-08 |r|=1.972e-09 (r)

# sub  2 [  2k] |x|=6.529e-05 |dx|=1.245e-12 |r|=5.009e-13 (c)

# all           |x|=1.899e+02 |dx|=2.346e-06 |r|=2.859e-09

# arc           |x|=5.703e+02 |dx|=8.677e-07 |r|=5.093e-11 (λ)


*** Continuation step 11

# SNES iteration  0

# sub  0 [  6k] |x|=2.068e+02 |dx|=2.345e-06 |r|=2.074e+01 (u)

# sub  1 [  6k] |x|=1.892e+00 |dx|=3.641e-08 |r|=1.029e+01 (r)

# sub  2 [  2k] |x|=6.925e-05 |dx|=1.245e-12 |r|=1.057e-12 (c)

# all           |x|=2.068e+02 |dx|=2.346e-06 |r|=2.315e+01

# arc           |x|=5.861e+02 |dx|=0.000e+00 |r|=5.093e-11 (λ)

# SNES iteration  0, KSP iteration   0       |r|=2.315e+01

# SNES iteration  0, KSP iteration   1       |r|=2.155e-09

# SNES iteration  1

# sub  0 [  6k] |x|=2.182e+02 |dx|=1.242e+01 |r|=4.254e+00 (u)

# sub  1 [  6k] |x|=1.954e+00 |dx|=8.341e-02 |r|=3.868e+00 (r)

# sub  2 [  2k] |x|=7.135e-05 |dx|=3.110e-06 |r|=2.133e-09 (c)

# all           |x|=2.182e+02 |dx|=1.242e+01 |r|=5.749e+00

# SNES iteration  1, KSP iteration   0       |r|=1.000e+00

# SNES iteration  1, KSP iteration   1       |r|=1.997e-10

# arc           |x|=5.805e+02 |dx|=5.634e+00 |r|=9.313e-10 (λ)

# SNES iteration  1, KSP iteration   0       |r|=3.421e-01

# SNES iteration  1, KSP iteration   1       |r|=5.986e-12

# SNES iteration  2

# sub  0 [  6k] |x|=2.065e+02 |dx|=5.704e-02 |r|=1.417e-04 (u)

# sub  1 [  6k] |x|=1.881e+00 |dx|=4.497e-04 |r|=7.199e-05 (r)

# sub  2 [  2k] |x|=6.887e-05 |dx|=1.545e-08 |r|=5.873e-12 (c)

# all           |x|=2.065e+02 |dx|=5.704e-02 |r|=1.590e-04

# SNES iteration  2, KSP iteration   0       |r|=1.000e+00

# SNES iteration  2, KSP iteration   1       |r|=4.052e-10

# arc           |x|=5.805e+02 |dx|=2.329e-02 |r|=3.929e-10 (λ)

# SNES iteration  2, KSP iteration   0       |r|=4.669e-05

# SNES iteration  2, KSP iteration   1       |r|=2.993e-16

# SNES iteration  3 success = CONVERGED_FNORM_RELATIVE

# sub  0 [  6k] |x|=2.066e+02 |dx|=2.795e-06 |r|=2.125e-09 (u)

# sub  1 [  6k] |x|=1.881e+00 |dx|=4.581e-08 |r|=2.086e-09 (r)

# sub  2 [  2k] |x|=6.888e-05 |dx|=1.625e-12 |r|=5.268e-13 (c)

# all           |x|=2.066e+02 |dx|=2.795e-06 |r|=2.978e-09

# arc           |x|=5.805e+02 |dx|=2.124e-07 |r|=2.910e-11 (λ)


*** Continuation step 12

# SNES iteration  0

# sub  0 [  6k] |x|=2.234e+02 |dx|=2.795e-06 |r|=1.882e+01 (u)

# sub  1 [  6k] |x|=1.991e+00 |dx|=4.581e-08 |r|=1.031e+01 (r)

# sub  2 [  2k] |x|=7.258e-05 |dx|=1.625e-12 |r|=1.130e-12 (c)

# all           |x|=2.234e+02 |dx|=2.795e-06 |r|=2.146e+01

# arc           |x|=5.906e+02 |dx|=0.000e+00 |r|=2.910e-11 (λ)

# SNES iteration  0, KSP iteration   0       |r|=2.146e+01

# SNES iteration  0, KSP iteration   1       |r|=9.640e-09

# SNES iteration  1

# sub  0 [  6k] |x|=2.834e+02 |dx|=6.386e+01 |r|=1.150e+02 (u)

# sub  1 [  6k] |x|=2.334e+00 |dx|=4.306e-01 |r|=7.123e+01 (r)

# sub  2 [  2k] |x|=8.417e-05 |dx|=1.569e-05 |r|=9.517e-09 (c)

# all           |x|=2.834e+02 |dx|=6.386e+01 |r|=1.353e+02

# SNES iteration  1, KSP iteration   0       |r|=1.000e+00

# SNES iteration  1, KSP iteration   1       |r|=6.068e-09

# arc           |x|=5.845e+02 |dx|=6.111e+00 |r|=8.513e-09 (λ)

# SNES iteration  1, KSP iteration   0       |r|=4.244e-01

# SNES iteration  1, KSP iteration   1       |r|=6.209e-11

# SNES iteration  2

# sub  0 [  6k] |x|=2.227e+02 |dx|=3.370e-01 |r|=3.294e-03 (u)

# sub  1 [  6k] |x|=1.978e+00 |dx|=2.301e-03 |r|=2.138e-03 (r)

# sub  2 [  2k] |x|=7.208e-05 |dx|=8.408e-08 |r|=6.155e-11 (c)

# all           |x|=2.228e+02 |dx|=3.370e-01 |r|=3.927e-03

# SNES iteration  2, KSP iteration   0       |r|=1.000e+00

# SNES iteration  2, KSP iteration   1       |r|=1.846e-09

# arc           |x|=5.845e+02 |dx|=1.621e-02 |r|=1.019e-10 (λ)

# SNES iteration  2, KSP iteration   0       |r|=5.083e-05

# SNES iteration  2, KSP iteration   1       |r|=5.521e-15

# SNES iteration  3 success = CONVERGED_FNORM_RELATIVE

# sub  0 [  6k] |x|=2.231e+02 |dx|=1.745e-05 |r|=2.079e-09 (u)

# sub  1 [  6k] |x|=1.979e+00 |dx|=1.308e-07 |r|=2.424e-09 (r)

# sub  2 [  2k] |x|=7.214e-05 |dx|=4.407e-12 |r|=5.172e-13 (c)

# all           |x|=2.231e+02 |dx|=1.745e-05 |r|=3.194e-09

# arc           |x|=5.845e+02 |dx|=7.961e-07 |r|=8.731e-11 (λ)


*** Continuation step 13

# SNES iteration  0

# sub  0 [  6k] |x|=2.397e+02 |dx|=1.745e-05 |r|=1.789e+01 (u)

# sub  1 [  6k] |x|=2.080e+00 |dx|=1.308e-07 |r|=1.105e+01 (r)

# sub  2 [  2k] |x|=7.552e-05 |dx|=4.407e-12 |r|=1.151e-12 (c)

# all           |x|=2.397e+02 |dx|=1.745e-05 |r|=2.102e+01

# arc           |x|=5.886e+02 |dx|=0.000e+00 |r|=8.731e-11 (λ)

# SNES iteration  0, KSP iteration   0       |r|=2.102e+01

# SNES iteration  0, KSP iteration   1       |r|=6.464e-09

# SNES iteration  1

# sub  0 [  6k] |x|=2.188e+02 |dx|=2.255e+01 |r|=1.445e+01 (u)

# sub  1 [  6k] |x|=1.962e+00 |dx|=1.523e-01 |r|=9.897e+00 (r)

# sub  2 [  2k] |x|=7.168e-05 |dx|=5.481e-06 |r|=6.435e-09 (c)

# all           |x|=2.188e+02 |dx|=2.255e+01 |r|=1.751e+01

# SNES iteration  1, KSP iteration   0       |r|=1.000e+00

# SNES iteration  1, KSP iteration   1       |r|=9.869e-10

# arc           |x|=5.820e+02 |dx|=6.606e+00 |r|=1.746e-10 (λ)

# SNES iteration  1, KSP iteration   0       |r|=5.264e-01

# SNES iteration  1, KSP iteration   1       |r|=1.287e-11

# SNES iteration  2

# sub  0 [  6k] |x|=2.394e+02 |dx|=3.353e-02 |r|=9.347e-05 (u)

# sub  1 [  6k] |x|=2.068e+00 |dx|=2.830e-04 |r|=7.266e-05 (r)

# sub  2 [  2k] |x|=7.501e-05 |dx|=8.570e-09 |r|=1.285e-11 (c)

# all           |x|=2.394e+02 |dx|=3.353e-02 |r|=1.184e-04

# SNES iteration  2, KSP iteration   0       |r|=1.000e+00

# SNES iteration  2, KSP iteration   1       |r|=6.637e-10

# arc           |x|=5.820e+02 |dx|=1.029e-02 |r|=5.093e-10 (λ)

# SNES iteration  2, KSP iteration   0       |r|=6.185e-05

# SNES iteration  2, KSP iteration   1       |r|=6.404e-16

# SNES iteration  3 success = CONVERGED_FNORM_RELATIVE

# sub  0 [  6k] |x|=2.393e+02 |dx|=8.406e-06 |r|=2.132e-09 (u)

# sub  1 [  6k] |x|=2.068e+00 |dx|=9.162e-08 |r|=2.511e-09 (r)

# sub  2 [  2k] |x|=7.501e-05 |dx|=3.714e-12 |r|=5.435e-13 (c)

# all           |x|=2.393e+02 |dx|=8.407e-06 |r|=3.294e-09

# arc           |x|=5.820e+02 |dx|=2.407e-06 |r|=3.711e-10 (λ)


*** Continuation step 14

# SNES iteration  0

# sub  0 [  6k] |x|=2.558e+02 |dx|=8.406e-06 |r|=1.819e+01 (u)

# sub  1 [  6k] |x|=2.160e+00 |dx|=9.162e-08 |r|=1.214e+01 (r)

# sub  2 [  2k] |x|=7.800e-05 |dx|=3.714e-12 |r|=1.160e-12 (c)

# all           |x|=2.558e+02 |dx|=8.407e-06 |r|=2.187e+01

# arc           |x|=5.795e+02 |dx|=0.000e+00 |r|=3.711e-10 (λ)

# SNES iteration  0, KSP iteration   0       |r|=2.187e+01

# SNES iteration  0, KSP iteration   1       |r|=1.125e-09

# SNES iteration  1

# sub  0 [  6k] |x|=2.464e+02 |dx|=1.002e+01 |r|=2.841e+00 (u)

# sub  1 [  6k] |x|=2.106e+00 |dx|=6.763e-02 |r|=1.529e+00 (r)

# sub  2 [  2k] |x|=7.618e-05 |dx|=2.445e-06 |r|=1.096e-09 (c)

# all           |x|=2.464e+02 |dx|=1.002e+01 |r|=3.226e+00

# SNES iteration  1, KSP iteration   0       |r|=1.000e+00

# SNES iteration  1, KSP iteration   1       |r|=1.496e-10

# arc           |x|=5.725e+02 |dx|=6.979e+00 |r|=1.965e-10 (λ)

# SNES iteration  1, KSP iteration   0       |r|=6.082e-01

# SNES iteration  1, KSP iteration   1       |r|=2.345e-12

# SNES iteration  2

# sub  0 [  6k] |x|=2.554e+02 |dx|=1.577e-02 |r|=8.160e-05 (u)

# sub  1 [  6k] |x|=2.149e+00 |dx|=2.244e-04 |r|=6.325e-05 (r)

# sub  2 [  2k] |x|=7.740e-05 |dx|=7.460e-09 |r|=2.384e-12 (c)

# all           |x|=2.554e+02 |dx|=1.577e-02 |r|=1.032e-04

# SNES iteration  2, KSP iteration   0       |r|=1.000e+00

# SNES iteration  2, KSP iteration   1       |r|=4.194e-10

# arc           |x|=5.725e+02 |dx|=4.905e-03 |r|=7.276e-11 (λ)

# SNES iteration  2, KSP iteration   0       |r|=8.982e-05

# SNES iteration  2, KSP iteration   1       |r|=1.587e-15

# SNES iteration  3 success = CONVERGED_FNORM_RELATIVE

# sub  0 [  6k] |x|=2.554e+02 |dx|=9.605e-06 |r|=2.231e-09 (u)

# sub  1 [  6k] |x|=2.148e+00 |dx|=1.203e-07 |r|=2.644e-09 (r)

# sub  2 [  2k] |x|=7.740e-05 |dx|=4.967e-12 |r|=5.474e-13 (c)

# all           |x|=2.554e+02 |dx|=9.606e-06 |r|=3.459e-09

# arc           |x|=5.725e+02 |dx|=4.955e-06 |r|=2.328e-10 (λ)


*** Continuation step 15

# SNES iteration  0

# sub  0 [  6k] |x|=2.716e+02 |dx|=9.605e-06 |r|=2.031e+01 (u)

# sub  1 [  6k] |x|=2.232e+00 |dx|=1.203e-07 |r|=1.304e+01 (r)

# sub  2 [  2k] |x|=7.995e-05 |dx|=4.967e-12 |r|=1.203e-12 (c)

# all           |x|=2.716e+02 |dx|=9.606e-06 |r|=2.414e+01

# arc           |x|=5.630e+02 |dx|=0.000e+00 |r|=2.328e-10 (λ)

# SNES iteration  0, KSP iteration   0       |r|=2.414e+01

# SNES iteration  0, KSP iteration   1       |r|=1.851e-09

# SNES iteration  1

# sub  0 [  6k] |x|=2.654e+02 |dx|=6.493e+00 |r|=1.194e+00 (u)

# sub  1 [  6k] |x|=2.196e+00 |dx|=4.382e-02 |r|=5.236e-01 (r)

# sub  2 [  2k] |x|=7.867e-05 |dx|=1.632e-06 |r|=1.844e-09 (c)

# all           |x|=2.654e+02 |dx|=6.493e+00 |r|=1.303e+00

# SNES iteration  1, KSP iteration   0       |r|=1.000e+00

# SNES iteration  1, KSP iteration   1       |r|=3.274e-10

# arc           |x|=5.559e+02 |dx|=7.045e+00 |r|=9.459e-11 (λ)

# SNES iteration  1, KSP iteration   0       |r|=6.129e-01

# SNES iteration  1, KSP iteration   1       |r|=1.064e-12

# SNES iteration  2

# sub  0 [  6k] |x|=2.711e+02 |dx|=1.768e-02 |r|=1.267e-04 (u)

# sub  1 [  6k] |x|=2.221e+00 |dx|=2.722e-04 |r|=7.793e-05 (r)

# sub  2 [  2k] |x|=7.925e-05 |dx|=1.030e-08 |r|=1.101e-12 (c)

# all           |x|=2.711e+02 |dx|=1.768e-02 |r|=1.487e-04

# SNES iteration  2, KSP iteration   0       |r|=1.000e+00

# SNES iteration  2, KSP iteration   1       |r|=1.039e-10

# arc           |x|=5.559e+02 |dx|=7.932e-04 |r|=3.056e-10 (λ)

# SNES iteration  2, KSP iteration   0       |r|=1.511e-04

# SNES iteration  2, KSP iteration   1       |r|=1.177e-15

# SNES iteration  3 success = CONVERGED_FNORM_RELATIVE

# sub  0 [  6k] |x|=2.711e+02 |dx|=1.190e-05 |r|=2.123e-09 (u)

# sub  1 [  6k] |x|=2.221e+00 |dx|=1.588e-07 |r|=2.797e-09 (r)

# sub  2 [  2k] |x|=7.925e-05 |dx|=6.714e-12 |r|=5.266e-13 (c)

# all           |x|=2.711e+02 |dx|=1.190e-05 |r|=3.511e-09

# arc           |x|=5.559e+02 |dx|=7.780e-06 |r|=1.746e-10 (λ)


*** Continuation step 16

# SNES iteration  0

# sub  0 [  6k] |x|=2.871e+02 |dx|=1.190e-05 |r|=2.474e+01 (u)

# sub  1 [  6k] |x|=2.298e+00 |dx|=1.588e-07 |r|=1.326e+01 (r)

# sub  2 [  2k] |x|=8.128e-05 |dx|=6.714e-12 |r|=1.214e-12 (c)

# all           |x|=2.871e+02 |dx|=1.190e-05 |r|=2.807e+01

# arc           |x|=5.394e+02 |dx|=0.000e+00 |r|=1.746e-10 (λ)

# SNES iteration  0, KSP iteration   0       |r|=2.807e+01

# SNES iteration  0, KSP iteration   1       |r|=1.049e-09

# SNES iteration  1

# sub  0 [  6k] |x|=2.826e+02 |dx|=4.672e+00 |r|=6.397e-01 (u)

# sub  1 [  6k] |x|=2.271e+00 |dx|=3.170e-02 |r|=2.512e-01 (r)

# sub  2 [  2k] |x|=8.026e-05 |dx|=1.262e-06 |r|=1.044e-09 (c)

# all           |x|=2.826e+02 |dx|=4.672e+00 |r|=6.873e-01

# SNES iteration  1, KSP iteration   0       |r|=1.000e+00

# SNES iteration  1, KSP iteration   1       |r|=1.488e-10

# arc           |x|=5.327e+02 |dx|=6.653e+00 |r|=1.965e-10 (λ)

# SNES iteration  1, KSP iteration   0       |r|=5.130e-01

# SNES iteration  1, KSP iteration   1       |r|=5.969e-12

# SNES iteration  2

# sub  0 [  6k] |x|=2.866e+02 |dx|=2.262e-02 |r|=2.056e-04 (u)

# sub  1 [  6k] |x|=2.288e+00 |dx|=3.378e-04 |r|=8.253e-05 (r)

# sub  2 [  2k] |x|=8.052e-05 |dx|=1.361e-08 |r|=5.941e-12 (c)

# all           |x|=2.866e+02 |dx|=2.262e-02 |r|=2.215e-04

# SNES iteration  2, KSP iteration   0       |r|=1.000e+00

# SNES iteration  2, KSP iteration   1       |r|=3.023e-10

# arc           |x|=5.327e+02 |dx|=5.773e-03 |r|=2.838e-10 (λ)

# SNES iteration  2, KSP iteration   0       |r|=2.428e-04

# SNES iteration  2, KSP iteration   1       |r|=1.290e-14

# SNES iteration  3 success = CONVERGED_FNORM_RELATIVE

# sub  0 [  6k] |x|=2.866e+02 |dx|=4.202e-05 |r|=2.177e-09 (u)

# sub  1 [  6k] |x|=2.288e+00 |dx|=4.221e-07 |r|=2.921e-09 (r)

# sub  2 [  2k] |x|=8.052e-05 |dx|=1.483e-11 |r|=5.450e-13 (c)

# all           |x|=2.866e+02 |dx|=4.202e-05 |r|=3.643e-09

# arc           |x|=5.327e+02 |dx|=7.585e-06 |r|=1.164e-10 (λ)


*** Continuation step 17

# SNES iteration  0

# sub  0 [  6k] |x|=3.023e+02 |dx|=4.202e-05 |r|=3.082e+01 (u)

# sub  1 [  6k] |x|=2.360e+00 |dx|=4.221e-07 |r|=1.269e+01 (r)

# sub  2 [  2k] |x|=8.199e-05 |dx|=1.483e-11 |r|=1.162e-12 (c)

# all           |x|=3.023e+02 |dx|=4.202e-05 |r|=3.333e+01

# arc           |x|=5.095e+02 |dx|=0.000e+00 |r|=1.164e-10 (λ)

# SNES iteration  0, KSP iteration   0       |r|=3.333e+01

# SNES iteration  0, KSP iteration   1       |r|=5.341e-10

# SNES iteration  1

# sub  0 [  6k] |x|=2.990e+02 |dx|=3.434e+00 |r|=3.851e-01 (u)

# sub  1 [  6k] |x|=2.340e+00 |dx|=2.373e-02 |r|=1.707e-01 (r)

# sub  2 [  2k] |x|=8.111e-05 |dx|=1.057e-06 |r|=5.277e-10 (c)

# all           |x|=2.990e+02 |dx|=3.434e+00 |r|=4.212e-01

# SNES iteration  1, KSP iteration   0       |r|=1.000e+00

# SNES iteration  1, KSP iteration   1       |r|=2.729e-10

# arc           |x|=5.037e+02 |dx|=5.809e+00 |r|=7.276e-11 (λ)

# SNES iteration  1, KSP iteration   0       |r|=3.569e-01

# SNES iteration  1, KSP iteration   1       |r|=2.677e-12

# SNES iteration  2

# sub  0 [  6k] |x|=3.019e+02 |dx|=2.567e-02 |r|=2.560e-04 (u)

# sub  1 [  6k] |x|=2.352e+00 |dx|=3.693e-04 |r|=6.495e-05 (r)

# sub  2 [  2k] |x|=8.119e-05 |dx|=1.554e-08 |r|=2.591e-12 (c)

# all           |x|=3.019e+02 |dx|=2.567e-02 |r|=2.641e-04

# SNES iteration  2, KSP iteration   0       |r|=1.000e+00

# SNES iteration  2, KSP iteration   1       |r|=1.064e-10

# arc           |x|=5.037e+02 |dx|=6.400e-03 |r|=1.673e-10 (λ)

# SNES iteration  2, KSP iteration   0       |r|=2.886e-04

# SNES iteration  2, KSP iteration   1       |r|=9.674e-15

# SNES iteration  3 success = CONVERGED_FNORM_RELATIVE

# sub  0 [  6k] |x|=3.019e+02 |dx|=2.323e-05 |r|=2.228e-09 (u)

# sub  1 [  6k] |x|=2.352e+00 |dx|=2.342e-07 |r|=3.104e-09 (r)

# sub  2 [  2k] |x|=8.119e-05 |dx|=9.037e-12 |r|=4.724e-13 (c)

# all           |x|=3.019e+02 |dx|=2.323e-05 |r|=3.821e-09

# arc           |x|=5.037e+02 |dx|=3.611e-06 |r|=1.310e-10 (λ)


*** Continuation step 18

# SNES iteration  0

# sub  0 [  6k] |x|=3.174e+02 |dx|=2.323e-05 |r|=3.684e+01 (u)

# sub  1 [  6k] |x|=2.421e+00 |dx|=2.342e-07 |r|=1.169e+01 (r)

# sub  2 [  2k] |x|=8.209e-05 |dx|=9.037e-12 |r|=1.102e-12 (c)

# all           |x|=3.174e+02 |dx|=2.323e-05 |r|=3.865e+01

# arc           |x|=4.747e+02 |dx|=0.000e+00 |r|=1.310e-10 (λ)

# SNES iteration  0, KSP iteration   0       |r|=3.865e+01

# SNES iteration  0, KSP iteration   1       |r|=1.791e-10

# SNES iteration  1

# sub  0 [  6k] |x|=3.150e+02 |dx|=2.491e+00 |r|=2.539e-01 (u)

# sub  1 [  6k] |x|=2.406e+00 |dx|=1.799e-02 |r|=1.510e-01 (r)

# sub  2 [  2k] |x|=8.129e-05 |dx|=9.346e-07 |r|=1.673e-10 (c)

# all           |x|=3.150e+02 |dx|=2.491e+00 |r|=2.954e-01

# SNES iteration  1, KSP iteration   0       |r|=1.000e+00

# SNES iteration  1, KSP iteration   1       |r|=1.559e-10

# arc           |x|=4.700e+02 |dx|=4.703e+00 |r|=3.638e-11 (λ)

# SNES iteration  1, KSP iteration   0       |r|=2.207e-01

# SNES iteration  1, KSP iteration   1       |r|=1.886e-12

# SNES iteration  2

# sub  0 [  6k] |x|=3.170e+02 |dx|=1.952e-02 |r|=1.694e-04 (u)

# sub  1 [  6k] |x|=2.414e+00 |dx|=2.988e-04 |r|=3.434e-05 (r)

# sub  2 [  2k] |x|=8.127e-05 |dx|=1.381e-08 |r|=1.892e-12 (c)

# all           |x|=3.170e+02 |dx|=1.952e-02 |r|=1.728e-04

# SNES iteration  2, KSP iteration   0       |r|=1.000e+00

# SNES iteration  2, KSP iteration   1       |r|=3.621e-11

# arc           |x|=4.700e+02 |dx|=1.142e-03 |r|=1.455e-11 (λ)

# SNES iteration  2, KSP iteration   0       |r|=1.762e-04

# SNES iteration  2, KSP iteration   1       |r|=3.409e-16

# SNES iteration  3 success = CONVERGED_FNORM_RELATIVE

# sub  0 [  6k] |x|=3.170e+02 |dx|=3.162e-06 |r|=2.296e-09 (u)

# sub  1 [  6k] |x|=2.414e+00 |dx|=4.771e-08 |r|=3.208e-09 (r)

# sub  2 [  2k] |x|=8.127e-05 |dx|=2.372e-12 |r|=4.494e-13 (c)

# all           |x|=3.170e+02 |dx|=3.162e-06 |r|=3.945e-09

# arc           |x|=4.700e+02 |dx|=5.986e-07 |r|=4.366e-11 (λ)


*** Continuation step 19

# SNES iteration  0

# sub  0 [  6k] |x|=3.325e+02 |dx|=3.162e-06 |r|=4.117e+01 (u)

# sub  1 [  6k] |x|=2.483e+00 |dx|=4.771e-08 |r|=1.081e+01 (r)

# sub  2 [  2k] |x|=8.161e-05 |dx|=2.372e-12 |r|=9.386e-13 (c)

# all           |x|=3.325e+02 |dx|=3.162e-06 |r|=4.256e+01

# arc           |x|=4.362e+02 |dx|=0.000e+00 |r|=4.366e-11 (λ)

# SNES iteration  0, KSP iteration   0       |r|=4.256e+01

# SNES iteration  0, KSP iteration   1       |r|=2.112e-10

# SNES iteration  1

# sub  0 [  6k] |x|=3.308e+02 |dx|=1.774e+00 |r|=1.852e-01 (u)

# sub  1 [  6k] |x|=2.471e+00 |dx|=1.400e-02 |r|=1.456e-01 (r)

# sub  2 [  2k] |x|=8.087e-05 |dx|=8.596e-07 |r|=2.070e-10 (c)

# all           |x|=3.308e+02 |dx|=1.774e+00 |r|=2.355e-01

# SNES iteration  1, KSP iteration   0       |r|=1.000e+00

# SNES iteration  1, KSP iteration   1       |r|=6.146e-11

# arc           |x|=4.326e+02 |dx|=3.578e+00 |r|=1.528e-10 (λ)

# SNES iteration  1, KSP iteration   0       |r|=1.331e-01

# SNES iteration  1, KSP iteration   1       |r|=1.693e-12

# SNES iteration  2

# sub  0 [  6k] |x|=3.322e+02 |dx|=1.167e-02 |r|=6.212e-05 (u)

# sub  1 [  6k] |x|=2.477e+00 |dx|=1.845e-04 |r|=1.114e-05 (r)

# sub  2 [  2k] |x|=8.080e-05 |dx|=1.006e-08 |r|=1.714e-12 (c)

# all           |x|=3.322e+02 |dx|=1.167e-02 |r|=6.311e-05

# SNES iteration  2, KSP iteration   0       |r|=1.000e+00

# SNES iteration  2, KSP iteration   1       |r|=2.541e-10

# arc           |x|=4.326e+02 |dx|=6.080e-03 |r|=3.274e-10 (λ)

# SNES iteration  2, KSP iteration   0       |r|=5.347e-05

# SNES iteration  2, KSP iteration   1       |r|=6.927e-17

# SNES iteration  3 success = CONVERGED_FNORM_RELATIVE

# sub  0 [  6k] |x|=3.322e+02 |dx|=8.584e-07 |r|=2.305e-09 (u)

# sub  1 [  6k] |x|=2.477e+00 |dx|=1.351e-08 |r|=3.375e-09 (r)

# sub  2 [  2k] |x|=8.080e-05 |dx|=8.150e-13 |r|=4.198e-13 (c)

# all           |x|=3.322e+02 |dx|=8.585e-07 |r|=4.087e-09

# arc           |x|=4.326e+02 |dx|=7.902e-07 |r|=3.347e-10 (λ)


*** Continuation step 20

# SNES iteration  0

# sub  0 [  6k] |x|=3.477e+02 |dx|=8.584e-07 |r|=4.316e+01 (u)

# sub  1 [  6k] |x|=2.547e+00 |dx|=1.351e-08 |r|=1.035e+01 (r)

# sub  2 [  2k] |x|=8.061e-05 |dx|=8.150e-13 |r|=8.594e-13 (c)

# all           |x|=3.477e+02 |dx|=8.585e-07 |r|=4.438e+01

# arc           |x|=3.953e+02 |dx|=0.000e+00 |r|=3.347e-10 (λ)

# SNES iteration  0, KSP iteration   0       |r|=4.438e+01

# SNES iteration  0, KSP iteration   1       |r|=7.056e-11

# SNES iteration  1

# sub  0 [  6k] |x|=3.465e+02 |dx|=1.256e+00 |r|=1.490e-01 (u)

# sub  1 [  6k] |x|=2.538e+00 |dx|=1.149e-02 |r|=1.395e-01 (r)

# sub  2 [  2k] |x|=7.989e-05 |dx|=8.147e-07 |r|=6.398e-11 (c)

# all           |x|=3.465e+02 |dx|=1.256e+00 |r|=2.041e-01

# SNES iteration  1, KSP iteration   0       |r|=1.000e+00

# SNES iteration  1, KSP iteration   1       |r|=1.729e-10

# arc           |x|=3.927e+02 |dx|=2.602e+00 |r|=5.821e-11 (λ)

# SNES iteration  1, KSP iteration   0       |r|=9.200e-02

# SNES iteration  1, KSP iteration   1       |r|=7.806e-13

# SNES iteration  2

# sub  0 [  6k] |x|=3.475e+02 |dx|=6.006e-03 |r|=8.357e-06 (u)

# sub  1 [  6k] |x|=2.543e+00 |dx|=7.920e-05 |r|=3.497e-06 (r)

# sub  2 [  2k] |x|=7.981e-05 |dx|=6.711e-09 |r|=8.486e-13 (c)

# all           |x|=3.475e+02 |dx|=6.006e-03 |r|=9.059e-06

# SNES iteration  2, KSP iteration   0       |r|=1.000e+00

# SNES iteration  2, KSP iteration   1       |r|=6.032e-11

# arc           |x|=3.928e+02 |dx|=1.058e-02 |r|=1.091e-10 (λ)

# SNES iteration  2, KSP iteration   0       |r|=5.368e-06

# SNES iteration  2, KSP iteration   1       |r|=5.289e-17

# SNES iteration  3 success = CONVERGED_FNORM_RELATIVE

# sub  0 [  6k] |x|=3.474e+02 |dx|=1.893e-07 |r|=2.307e-09 (u)

# sub  1 [  6k] |x|=2.543e+00 |dx|=3.002e-09 |r|=3.434e-09 (r)

# sub  2 [  2k] |x|=7.981e-05 |dx|=3.344e-13 |r|=3.917e-13 (c)

# all           |x|=3.475e+02 |dx|=1.893e-07 |r|=4.137e-09

# arc           |x|=3.928e+02 |dx|=2.372e-07 |r|=1.310e-10 (λ)


*** Continuation step 21

# SNES iteration  0

# sub  0 [  6k] |x|=3.629e+02 |dx|=1.893e-07 |r|=4.312e+01 (u)

# sub  1 [  6k] |x|=2.614e+00 |dx|=3.002e-09 |r|=1.023e+01 (r)

# sub  2 [  2k] |x|=7.912e-05 |dx|=3.344e-13 |r|=8.001e-13 (c)

# all           |x|=3.630e+02 |dx|=1.893e-07 |r|=4.431e+01

# arc           |x|=3.529e+02 |dx|=0.000e+00 |r|=1.310e-10 (λ)

# SNES iteration  0, KSP iteration   0       |r|=4.431e+01

# SNES iteration  0, KSP iteration   1       |r|=8.117e-11

# SNES iteration  1

# sub  0 [  6k] |x|=3.621e+02 |dx|=9.094e-01 |r|=1.280e-01 (u)

# sub  1 [  6k] |x|=2.607e+00 |dx|=1.007e-02 |r|=1.287e-01 (r)

# sub  2 [  2k] |x|=7.842e-05 |dx|=7.900e-07 |r|=7.813e-11 (c)

# all           |x|=3.621e+02 |dx|=9.094e-01 |r|=1.815e-01

# SNES iteration  1, KSP iteration   0       |r|=1.000e+00

# SNES iteration  1, KSP iteration   1       |r|=5.687e-11

# arc           |x|=3.510e+02 |dx|=1.836e+00 |r|=2.910e-11 (λ)

# SNES iteration  1, KSP iteration   0       |r|=9.211e-02

# SNES iteration  1, KSP iteration   1       |r|=3.784e-12

# SNES iteration  2

# sub  0 [  6k] |x|=3.628e+02 |dx|=5.745e-03 |r|=4.964e-06 (u)

# sub  1 [  6k] |x|=2.610e+00 |dx|=6.408e-05 |r|=4.099e-06 (r)

# sub  2 [  2k] |x|=7.834e-05 |dx|=5.530e-09 |r|=3.754e-12 (c)

# all           |x|=3.628e+02 |dx|=5.745e-03 |r|=6.437e-06

# SNES iteration  2, KSP iteration   0       |r|=1.000e+00

# SNES iteration  2, KSP iteration   1       |r|=2.163e-10

# arc           |x|=3.510e+02 |dx|=1.080e-02 |r|=1.746e-10 (λ)

# SNES iteration  2, KSP iteration   0       |r|=1.014e-05

# SNES iteration  2, KSP iteration   1       |r|=3.057e-17

# SNES iteration  3 success = CONVERGED_FNORM_RELATIVE

# sub  0 [  6k] |x|=3.628e+02 |dx|=2.632e-07 |r|=2.316e-09 (u)

# sub  1 [  6k] |x|=2.610e+00 |dx|=4.182e-09 |r|=3.720e-09 (r)

# sub  2 [  2k] |x|=7.834e-05 |dx|=2.804e-13 |r|=3.417e-13 (c)

# all           |x|=3.628e+02 |dx|=2.632e-07 |r|=4.382e-09

# arc           |x|=3.510e+02 |dx|=1.882e-07 |r|=1.601e-10 (λ)


*** Continuation step 22

# SNES iteration  0

# sub  0 [  6k] |x|=3.784e+02 |dx|=2.632e-07 |r|=4.176e+01 (u)

# sub  1 [  6k] |x|=2.684e+00 |dx|=4.182e-09 |r|=1.026e+01 (r)

# sub  2 [  2k] |x|=7.719e-05 |dx|=2.804e-13 |r|=7.292e-13 (c)

# all           |x|=3.784e+02 |dx|=2.632e-07 |r|=4.300e+01

# arc           |x|=3.093e+02 |dx|=0.000e+00 |r|=1.601e-10 (λ)

# SNES iteration  0, KSP iteration   0       |r|=4.300e+01

# SNES iteration  0, KSP iteration   1       |r|=3.483e-11

# SNES iteration  1

# sub  0 [  6k] |x|=3.778e+02 |dx|=6.960e-01 |r|=1.136e-01 (u)

# sub  1 [  6k] |x|=2.678e+00 |dx|=9.327e-03 |r|=1.140e-01 (r)

# sub  2 [  2k] |x|=7.649e-05 |dx|=7.815e-07 |r|=2.989e-11 (c)

# all           |x|=3.778e+02 |dx|=6.960e-01 |r|=1.610e-01

# SNES iteration  1, KSP iteration   0       |r|=1.000e+00

# SNES iteration  1, KSP iteration   1       |r|=1.987e-10

# arc           |x|=3.080e+02 |dx|=1.269e+00 |r|=1.164e-10 (λ)

# SNES iteration  1, KSP iteration   0       |r|=1.049e-01

# SNES iteration  1, KSP iteration   1       |r|=9.742e-13

# SNES iteration  2

# sub  0 [  6k] |x|=3.783e+02 |dx|=6.510e-03 |r|=1.604e-05 (u)

# sub  1 [  6k] |x|=2.680e+00 |dx|=9.434e-05 |r|=6.780e-06 (r)

# sub  2 [  2k] |x|=7.643e-05 |dx|=5.708e-09 |r|=1.021e-12 (c)

# all           |x|=3.783e+02 |dx|=6.511e-03 |r|=1.741e-05

# SNES iteration  2, KSP iteration   0       |r|=1.000e+00

# SNES iteration  2, KSP iteration   1       |r|=1.430e-10

# arc           |x|=3.080e+02 |dx|=7.749e-03 |r|=4.220e-10 (λ)

# SNES iteration  2, KSP iteration   0       |r|=2.337e-05

# SNES iteration  2, KSP iteration   1       |r|=2.306e-16

# SNES iteration  3 success = CONVERGED_FNORM_RELATIVE

# sub  0 [  6k] |x|=3.783e+02 |dx|=4.786e-07 |r|=2.369e-09 (u)

# sub  1 [  6k] |x|=2.680e+00 |dx|=7.574e-09 |r|=3.875e-09 (r)

# sub  2 [  2k] |x|=7.643e-05 |dx|=3.706e-13 |r|=3.191e-13 (c)

# all           |x|=3.783e+02 |dx|=4.787e-07 |r|=4.542e-09

# arc           |x|=3.080e+02 |dx|=2.597e-07 |r|=6.548e-11 (λ)


*** Continuation step 23

# SNES iteration  0

# sub  0 [  6k] |x|=3.939e+02 |dx|=4.786e-07 |r|=3.984e+01 (u)

# sub  1 [  6k] |x|=2.756e+00 |dx|=7.574e-09 |r|=1.028e+01 (r)

# sub  2 [  2k] |x|=7.486e-05 |dx|=3.706e-13 |r|=6.899e-13 (c)

# all           |x|=3.939e+02 |dx|=4.787e-07 |r|=4.114e+01

# arc           |x|=2.651e+02 |dx|=0.000e+00 |r|=6.548e-11 (λ)

# SNES iteration  0, KSP iteration   0       |r|=4.114e+01

# SNES iteration  0, KSP iteration   1       |r|=1.897e-10

# SNES iteration  1

# sub  0 [  6k] |x|=3.935e+02 |dx|=5.763e-01 |r|=1.027e-01 (u)

# sub  1 [  6k] |x|=2.750e+00 |dx|=8.951e-03 |r|=9.866e-02 (r)

# sub  2 [  2k] |x|=7.416e-05 |dx|=7.901e-07 |r|=1.892e-10 (c)

# all           |x|=3.935e+02 |dx|=5.764e-01 |r|=1.424e-01

# SNES iteration  1, KSP iteration   0       |r|=1.000e+00

# SNES iteration  1, KSP iteration   1       |r|=7.527e-11

# arc           |x|=2.642e+02 |dx|=8.611e-01 |r|=4.366e-10 (λ)

# SNES iteration  1, KSP iteration   0       |r|=1.138e-01

# SNES iteration  1, KSP iteration   1       |r|=4.477e-13

# SNES iteration  2

# sub  0 [  6k] |x|=3.938e+02 |dx|=6.061e-03 |r|=2.135e-05 (u)

# sub  1 [  6k] |x|=2.752e+00 |dx|=1.039e-04 |r|=8.454e-06 (r)

# sub  2 [  2k] |x|=7.411e-05 |dx|=6.009e-09 |r|=5.054e-13 (c)

# all           |x|=3.938e+02 |dx|=6.062e-03 |r|=2.296e-05

# SNES iteration  2, KSP iteration   0       |r|=1.000e+00

# SNES iteration  2, KSP iteration   1       |r|=3.852e-11

# arc           |x|=2.642e+02 |dx|=3.069e-03 |r|=2.910e-11 (λ)

# SNES iteration  2, KSP iteration   0       |r|=2.545e-05

# SNES iteration  2, KSP iteration   1       |r|=1.051e-16

# SNES iteration  3 success = CONVERGED_FNORM_RELATIVE

# sub  0 [  6k] |x|=3.938e+02 |dx|=5.123e-07 |r|=2.391e-09 (u)

# sub  1 [  6k] |x|=2.752e+00 |dx|=8.410e-09 |r|=4.236e-09 (r)

# sub  2 [  2k] |x|=7.411e-05 |dx|=4.071e-13 |r|=2.755e-13 (c)

# all           |x|=3.938e+02 |dx|=5.123e-07 |r|=4.865e-09

# arc           |x|=2.642e+02 |dx|=1.843e-07 |r|=2.401e-10 (λ)


*** Continuation step 24

# SNES iteration  0

# sub  0 [  6k] |x|=4.096e+02 |dx|=5.123e-07 |r|=3.793e+01 (u)

# sub  1 [  6k] |x|=2.829e+00 |dx|=8.410e-09 |r|=1.025e+01 (r)

# sub  2 [  2k] |x|=7.218e-05 |dx|=4.071e-13 |r|=5.773e-13 (c)

# all           |x|=4.096e+02 |dx|=5.123e-07 |r|=3.929e+01

# arc           |x|=2.204e+02 |dx|=0.000e+00 |r|=2.401e-10 (λ)

# SNES iteration  0, KSP iteration   0       |r|=3.929e+01

# SNES iteration  0, KSP iteration   1       |r|=3.029e-10

# SNES iteration  1

# sub  0 [  6k] |x|=4.093e+02 |dx|=5.182e-01 |r|=9.516e-02 (u)

# sub  1 [  6k] |x|=2.824e+00 |dx|=8.819e-03 |r|=8.667e-02 (r)

# sub  2 [  2k] |x|=7.147e-05 |dx|=8.193e-07 |r|=3.026e-10 (c)

# all           |x|=4.093e+02 |dx|=5.183e-01 |r|=1.287e-01

# SNES iteration  1, KSP iteration   0       |r|=1.000e+00

# SNES iteration  1, KSP iteration   1       |r|=6.101e-11

# arc           |x|=2.198e+02 |dx|=5.649e-01 |r|=1.164e-10 (λ)

# SNES iteration  1, KSP iteration   0       |r|=1.189e-01

# SNES iteration  1, KSP iteration   1       |r|=4.476e-13

# SNES iteration  2

# sub  0 [  6k] |x|=4.095e+02 |dx|=5.001e-03 |r|=1.829e-05 (u)

# sub  1 [  6k] |x|=2.825e+00 |dx|=9.585e-05 |r|=8.507e-06 (r)

# sub  2 [  2k] |x|=7.143e-05 |dx|=6.573e-09 |r|=4.863e-13 (c)

# all           |x|=4.095e+02 |dx|=5.002e-03 |r|=2.017e-05

# SNES iteration  2, KSP iteration   0       |r|=1.000e+00

# SNES iteration  2, KSP iteration   1       |r|=2.714e-10

# arc           |x|=2.198e+02 |dx|=1.958e-03 |r|=2.619e-10 (λ)

# SNES iteration  2, KSP iteration   0       |r|=1.890e-05

# SNES iteration  2, KSP iteration   1       |r|=8.489e-17

# SNES iteration  3 success = CONVERGED_FNORM_RELATIVE

# sub  0 [  6k] |x|=4.095e+02 |dx|=4.050e-07 |r|=2.446e-09 (u)

# sub  1 [  6k] |x|=2.825e+00 |dx|=7.148e-09 |r|=4.247e-09 (r)

# sub  2 [  2k] |x|=7.143e-05 |dx|=4.293e-13 |r|=2.110e-13 (c)

# all           |x|=4.095e+02 |dx|=4.050e-07 |r|=4.901e-09

# arc           |x|=2.198e+02 |dx|=3.769e-08 |r|=4.293e-10 (λ)


*** Continuation step 25

# SNES iteration  0

# sub  0 [  6k] |x|=4.253e+02 |dx|=4.050e-07 |r|=3.645e+01 (u)

# sub  1 [  6k] |x|=2.905e+00 |dx|=7.148e-09 |r|=1.014e+01 (r)

# sub  2 [  2k] |x|=6.922e-05 |dx|=4.293e-13 |r|=4.732e-13 (c)

# all           |x|=4.253e+02 |dx|=4.050e-07 |r|=3.784e+01

# arc           |x|=1.754e+02 |dx|=0.000e+00 |r|=4.293e-10 (λ)

# SNES iteration  0, KSP iteration   0       |r|=3.784e+01

# SNES iteration  0, KSP iteration   1       |r|=3.555e-10

# SNES iteration  1

# sub  0 [  6k] |x|=4.250e+02 |dx|=5.021e-01 |r|=9.240e-02 (u)

# sub  1 [  6k] |x|=2.900e+00 |dx|=8.930e-03 |r|=8.288e-02 (r)

# sub  2 [  2k] |x|=6.849e-05 |dx|=8.747e-07 |r|=3.553e-10 (c)

# all           |x|=4.250e+02 |dx|=5.021e-01 |r|=1.241e-01

# SNES iteration  1, KSP iteration   0       |r|=1.000e+00

# SNES iteration  1, KSP iteration   1       |r|=9.118e-11

# arc           |x|=1.751e+02 |dx|=3.381e-01 |r|=3.711e-10 (λ)

# SNES iteration  1, KSP iteration   0       |r|=1.250e-01

# SNES iteration  1, KSP iteration   1       |r|=3.382e-13

# SNES iteration  2

# sub  0 [  6k] |x|=4.251e+02 |dx|=4.782e-03 |r|=1.176e-05 (u)

# sub  1 [  6k] |x|=2.900e+00 |dx|=8.277e-05 |r|=7.889e-06 (r)

# sub  2 [  2k] |x|=6.846e-05 |dx|=8.248e-09 |r|=3.699e-13 (c)

# all           |x|=4.251e+02 |dx|=4.783e-03 |r|=1.416e-05

# SNES iteration  2, KSP iteration   0       |r|=1.000e+00

# SNES iteration  2, KSP iteration   1       |r|=1.969e-10

# arc           |x|=1.750e+02 |dx|=6.554e-03 |r|=6.185e-10 (λ)

# SNES iteration  2, KSP iteration   0       |r|=1.208e-05

# SNES iteration  2, KSP iteration   1       |r|=4.700e-17

# SNES iteration  3 success = CONVERGED_FNORM_RELATIVE

# sub  0 [  6k] |x|=4.251e+02 |dx|=2.755e-07 |r|=2.476e-09 (u)

# sub  1 [  6k] |x|=2.900e+00 |dx|=5.518e-09 |r|=4.545e-09 (r)

# sub  2 [  2k] |x|=6.846e-05 |dx|=6.344e-13 |r|=1.705e-13 (c)

# all           |x|=4.251e+02 |dx|=2.755e-07 |r|=5.176e-09

# arc           |x|=1.750e+02 |dx|=2.188e-09 |r|=9.459e-11 (λ)


*** Continuation step 26

# SNES iteration  0

# sub  0 [  6k] |x|=4.410e+02 |dx|=2.755e-07 |r|=3.569e+01 (u)

# sub  1 [  6k] |x|=2.981e+00 |dx|=5.518e-09 |r|=9.987e+00 (r)

# sub  2 [  2k] |x|=6.608e-05 |dx|=6.344e-13 |r|=3.555e-13 (c)

# all           |x|=4.410e+02 |dx|=2.755e-07 |r|=3.706e+01

# arc           |x|=1.303e+02 |dx|=0.000e+00 |r|=9.459e-11 (λ)

# SNES iteration  0, KSP iteration   0       |r|=3.706e+01

# SNES iteration  0, KSP iteration   1       |r|=8.711e-11

# SNES iteration  1

# sub  0 [  6k] |x|=4.408e+02 |dx|=5.211e-01 |r|=9.668e-02 (u)

# sub  1 [  6k] |x|=2.977e+00 |dx|=9.360e-03 |r|=9.192e-02 (r)

# sub  2 [  2k] |x|=6.534e-05 |dx|=9.646e-07 |r|=8.626e-11 (c)

# all           |x|=4.408e+02 |dx|=5.212e-01 |r|=1.334e-01

# SNES iteration  1, KSP iteration   0       |r|=1.000e+00

# SNES iteration  1, KSP iteration   1       |r|=1.091e-10

# arc           |x|=1.302e+02 |dx|=1.430e-01 |r|=2.838e-10 (λ)

# SNES iteration  1, KSP iteration   0       |r|=1.369e-01

# SNES iteration  1, KSP iteration   1       |r|=1.061e-12

# SNES iteration  2

# sub  0 [  6k] |x|=4.408e+02 |dx|=6.803e-03 |r|=8.823e-06 (u)

# sub  1 [  6k] |x|=2.977e+00 |dx|=8.628e-05 |r|=7.846e-06 (r)

# sub  2 [  2k] |x|=6.533e-05 |dx|=1.187e-08 |r|=9.873e-13 (c)

# all           |x|=4.408e+02 |dx|=6.804e-03 |r|=1.181e-05

# SNES iteration  2, KSP iteration   0       |r|=1.000e+00

# SNES iteration  2, KSP iteration   1       |r|=1.005e-10

# arc           |x|=1.301e+02 |dx|=1.020e-02 |r|=6.548e-11 (λ)

# SNES iteration  2, KSP iteration   0       |r|=1.303e-05

# SNES iteration  2, KSP iteration   1       |r|=5.211e-16

# SNES iteration  3 success = CONVERGED_FNORM_RELATIVE

# sub  0 [  6k] |x|=4.408e+02 |dx|=5.144e-07 |r|=2.449e-09 (u)

# sub  1 [  6k] |x|=2.977e+00 |dx|=8.756e-09 |r|=4.418e-09 (r)

# sub  2 [  2k] |x|=6.532e-05 |dx|=1.262e-12 |r|=1.183e-13 (c)

# all           |x|=4.408e+02 |dx|=5.145e-07 |r|=5.051e-09

# arc           |x|=1.301e+02 |dx|=2.831e-07 |r|=2.547e-10 (λ)


*** Continuation step 27

# SNES iteration  0

# sub  0 [  6k] |x|=4.566e+02 |dx|=5.144e-07 |r|=3.586e+01 (u)

# sub  1 [  6k] |x|=3.058e+00 |dx|=8.756e-09 |r|=9.790e+00 (r)

# sub  2 [  2k] |x|=6.296e-05 |dx|=1.262e-12 |r|=2.573e-13 (c)

# all           |x|=4.566e+02 |dx|=5.145e-07 |r|=3.717e+01

# arc           |x|=8.523e+01 |dx|=0.000e+00 |r|=2.547e-10 (λ)

# SNES iteration  0, KSP iteration   0       |r|=3.717e+01

# SNES iteration  0, KSP iteration   1       |r|=1.716e-10

# SNES iteration  1

# sub  0 [  6k] |x|=4.564e+02 |dx|=5.781e-01 |r|=1.114e-01 (u)

# sub  1 [  6k] |x|=3.054e+00 |dx|=1.025e-02 |r|=1.168e-01 (r)

# sub  2 [  2k] |x|=6.223e-05 |dx|=1.101e-06 |r|=1.711e-10 (c)

# all           |x|=4.564e+02 |dx|=5.782e-01 |r|=1.614e-01

# SNES iteration  1, KSP iteration   0       |r|=1.000e+00

# SNES iteration  1, KSP iteration   1       |r|=2.198e-10

# arc           |x|=8.529e+01 |dx|=5.610e-02 |r|=8.004e-11 (λ)

# SNES iteration  1, KSP iteration   0       |r|=1.588e-01

# SNES iteration  1, KSP iteration   1       |r|=4.371e-12

# SNES iteration  2

# sub  0 [  6k] |x|=4.564e+02 |dx|=1.099e-02 |r|=1.606e-05 (u)

# sub  1 [  6k] |x|=3.054e+00 |dx|=1.294e-04 |r|=9.941e-06 (r)

# sub  2 [  2k] |x|=6.223e-05 |dx|=1.818e-08 |r|=1.112e-12 (c)

# all           |x|=4.564e+02 |dx|=1.100e-02 |r|=1.889e-05

# SNES iteration  2, KSP iteration   0       |r|=1.000e+00

# SNES iteration  2, KSP iteration   1       |r|=9.040e-11

# arc           |x|=8.528e+01 |dx|=1.228e-02 |r|=2.110e-10 (λ)

# SNES iteration  2, KSP iteration   0       |r|=2.512e-05

# SNES iteration  2, KSP iteration   1       |r|=9.071e-16

# SNES iteration  3 success = CONVERGED_FNORM_RELATIVE

# sub  0 [  6k] |x|=4.564e+02 |dx|=1.443e-06 |r|=2.466e-09 (u)

# sub  1 [  6k] |x|=3.054e+00 |dx|=2.252e-08 |r|=4.478e-09 (r)

# sub  2 [  2k] |x|=6.222e-05 |dx|=2.741e-12 |r|=7.010e-14 (c)

# all           |x|=4.564e+02 |dx|=1.443e-06 |r|=5.112e-09

# arc           |x|=8.528e+01 |dx|=1.303e-06 |r|=4.366e-11 (λ)


*** Continuation step 28

# SNES iteration  0

# sub  0 [  6k] |x|=4.722e+02 |dx|=1.443e-06 |r|=3.716e+01 (u)

# sub  1 [  6k] |x|=3.137e+00 |dx|=2.252e-08 |r|=9.569e+00 (r)

# sub  2 [  2k] |x|=6.018e-05 |dx|=2.741e-12 |r|=1.811e-13 (c)

# all           |x|=4.722e+02 |dx|=1.443e-06 |r|=3.837e+01

# arc           |x|=4.042e+01 |dx|=0.000e+00 |r|=4.366e-11 (λ)

# SNES iteration  0, KSP iteration   0       |r|=3.837e+01

# SNES iteration  0, KSP iteration   1       |r|=2.932e-10

# SNES iteration  1

# sub  0 [  6k] |x|=4.720e+02 |dx|=6.828e-01 |r|=1.426e-01 (u)

# sub  1 [  6k] |x|=3.133e+00 |dx|=1.180e-02 |r|=1.609e-01 (r)

# sub  2 [  2k] |x|=5.950e-05 |dx|=1.301e-06 |r|=2.926e-10 (c)

# all           |x|=4.720e+02 |dx|=6.829e-01 |r|=2.150e-01

# SNES iteration  1, KSP iteration   0       |r|=1.000e+00

# SNES iteration  1, KSP iteration   1       |r|=1.872e-10

# arc           |x|=4.072e+01 |dx|=2.994e-01 |r|=8.004e-11 (λ)

# SNES iteration  1, KSP iteration   0       |r|=1.945e-01

# SNES iteration  1, KSP iteration   1       |r|=1.720e-12

# SNES iteration  2

# sub  0 [  6k] |x|=4.719e+02 |dx|=1.756e-02 |r|=4.522e-05 (u)

# sub  1 [  6k] |x|=3.132e+00 |dx|=2.196e-04 |r|=1.712e-05 (r)

# sub  2 [  2k] |x|=5.951e-05 |dx|=2.856e-08 |r|=1.657e-12 (c)

# all           |x|=4.719e+02 |dx|=1.756e-02 |r|=4.836e-05

# SNES iteration  2, KSP iteration   0       |r|=1.000e+00

# SNES iteration  2, KSP iteration   1       |r|=1.061e-10

# arc           |x|=4.070e+01 |dx|=1.160e-02 |r|=7.276e-11 (λ)

# SNES iteration  2, KSP iteration   0       |r|=6.191e-05

# SNES iteration  2, KSP iteration   1       |r|=5.320e-15

# SNES iteration  3 success = CONVERGED_FNORM_RELATIVE

# sub  0 [  6k] |x|=4.719e+02 |dx|=5.523e-06 |r|=2.618e-09 (u)

# sub  1 [  6k] |x|=3.132e+00 |dx|=8.188e-08 |r|=4.914e-09 (r)

# sub  2 [  2k] |x|=5.951e-05 |dx|=7.957e-12 |r|=5.756e-14 (c)

# all           |x|=4.719e+02 |dx|=5.523e-06 |r|=5.568e-09

# arc           |x|=4.070e+01 |dx|=4.088e-06 |r|=1.455e-10 (λ)


*** Continuation step 29

# SNES iteration  0

# sub  0 [  6k] |x|=4.875e+02 |dx|=5.523e-06 |r|=3.982e+01 (u)

# sub  1 [  6k] |x|=3.216e+00 |dx|=8.188e-08 |r|=9.355e+00 (r)

# sub  2 [  2k] |x|=5.826e-05 |dx|=7.957e-12 |r|=1.602e-13 (c)

# all           |x|=4.875e+02 |dx|=5.523e-06 |r|=4.090e+01

# arc           |x|=3.870e+00 |dx|=0.000e+00 |r|=1.455e-10 (λ)

# SNES iteration  0, KSP iteration   0       |r|=4.090e+01

# SNES iteration  0, KSP iteration   1       |r|=2.029e-10

# SNES iteration  1

# sub  0 [  6k] |x|=4.874e+02 |dx|=8.526e-01 |r|=2.030e-01 (u)

# sub  1 [  6k] |x|=3.212e+00 |dx|=1.433e-02 |r|=2.327e-01 (r)

# sub  2 [  2k] |x|=5.773e-05 |dx|=1.593e-06 |r|=2.019e-10 (c)

# all           |x|=4.874e+02 |dx|=8.527e-01 |r|=3.088e-01

# SNES iteration  1, KSP iteration   0       |r|=1.000e+00

# SNES iteration  1, KSP iteration   1       |r|=9.129e-11

# arc           |x|=3.231e+00 |dx|=6.391e-01 |r|=4.366e-11 (λ)

# SNES iteration  1, KSP iteration   0       |r|=2.488e-01

# SNES iteration  1, KSP iteration   1       |r|=5.345e-11

# SNES iteration  2

# sub  0 [  6k] |x|=4.871e+02 |dx|=5.633e-02 |r|=6.522e-04 (u)

# sub  1 [  6k] |x|=3.211e+00 |dx|=7.989e-04 |r|=2.110e-04 (r)

# sub  2 [  2k] |x|=5.774e-05 |dx|=7.473e-08 |r|=5.343e-11 (c)

# all           |x|=4.872e+02 |dx|=5.634e-02 |r|=6.855e-04

# SNES iteration  2, KSP iteration   0       |r|=1.000e+00

# SNES iteration  2, KSP iteration   1       |r|=5.597e-11

# arc           |x|=3.236e+00 |dx|=5.227e-03 |r|=2.401e-10 (λ)

# SNES iteration  2, KSP iteration   0       |r|=7.052e-04

# SNES iteration  2, KSP iteration   1       |r|=4.055e-13

# SNES iteration  3 success = CONVERGED_FNORM_RELATIVE

# sub  0 [  6k] |x|=4.871e+02 |dx|=8.055e-04 |r|=1.167e-07 (u)

# sub  1 [  6k] |x|=3.211e+00 |dx|=1.165e-05 |r|=4.300e-08 (r)

# sub  2 [  2k] |x|=5.774e-05 |dx|=9.767e-10 |r|=4.167e-13 (c)

# all           |x|=4.872e+02 |dx|=8.055e-04 |r|=1.244e-07

# arc           |x|=3.236e+00 |dx|=2.264e-05 |r|=1.601e-10 (λ)


*** Continuation step 30

# SNES iteration  0

# sub  0 [  6k] |x|=5.026e+02 |dx|=8.055e-04 |r|=4.411e+01 (u)

# sub  1 [  6k] |x|=3.296e+00 |dx|=1.165e-05 |r|=9.224e+00 (r)

# sub  2 [  2k] |x|=5.801e-05 |dx|=9.767e-10 |r|=8.441e-13 (c)

# all           |x|=5.026e+02 |dx|=8.055e-04 |r|=4.506e+01

# arc           |x|=4.718e+01 |dx|=0.000e+00 |r|=1.601e-10 (λ)

# SNES iteration  0, KSP iteration   0       |r|=4.506e+01

# SNES iteration  0, KSP iteration   1       |r|=1.841e-10

# SNES iteration  1

# sub  0 [  6k] |x|=5.025e+02 |dx|=1.287e+00 |r|=4.619e-01 (u)

# sub  1 [  6k] |x|=3.293e+00 |dx|=2.052e-02 |r|=4.066e-01 (r)

# sub  2 [  2k] |x|=5.783e-05 |dx|=2.161e-06 |r|=1.817e-10 (c)

# all           |x|=5.025e+02 |dx|=1.287e+00 |r|=6.154e-01

# SNES iteration  1, KSP iteration   0       |r|=1.000e+00

# SNES iteration  1, KSP iteration   1       |r|=5.775e-11

# arc           |x|=4.601e+01 |dx|=1.165e+00 |r|=8.004e-11 (λ)

# SNES iteration  1, KSP iteration   0       |r|=4.583e-01

# SNES iteration  1, KSP iteration   1       |r|=1.900e-10

# SNES iteration  2

# sub  0 [  6k] |x|=5.021e+02 |dx|=3.932e-01 |r|=2.764e-02 (u)

# sub  1 [  6k] |x|=3.291e+00 |dx|=5.676e-03 |r|=9.419e-03 (r)

# sub  2 [  2k] |x|=5.780e-05 |dx|=4.811e-07 |r|=1.898e-10 (c)

# all           |x|=5.021e+02 |dx|=3.932e-01 |r|=2.920e-02

# SNES iteration  2, KSP iteration   0       |r|=1.000e+00

# SNES iteration  2, KSP iteration   1       |r|=2.099e-10

# arc           |x|=4.602e+01 |dx|=7.839e-03 |r|=2.256e-10 (λ)

# SNES iteration  2, KSP iteration   0       |r|=2.928e-02

# SNES iteration  2, KSP iteration   1       |r|=2.158e-12

# SNES iteration  3

# sub  0 [  6k] |x|=5.021e+02 |dx|=3.601e-03 |r|=2.540e-06 (u)

# sub  1 [  6k] |x|=3.291e+00 |dx|=5.250e-05 |r|=1.078e-06 (r)

# sub  2 [  2k] |x|=5.780e-05 |dx|=4.519e-09 |r|=2.170e-12 (c)

# all           |x|=5.021e+02 |dx|=3.602e-03 |r|=2.759e-06

# SNES iteration  3, KSP iteration   0       |r|=1.000e+00

# SNES iteration  3, KSP iteration   1       |r|=3.163e-10

# arc           |x|=4.602e+01 |dx|=4.259e-04 |r|=1.455e-10 (λ)

# SNES iteration  3, KSP iteration   0       |r|=2.698e-06

# SNES iteration  3, KSP iteration   1       |r|=1.019e-15

# SNES iteration  4 success = CONVERGED_FNORM_RELATIVE

# sub  0 [  6k] |x|=5.021e+02 |dx|=1.897e-06 |r|=2.619e-09 (u)

# sub  1 [  6k] |x|=3.291e+00 |dx|=2.740e-08 |r|=5.426e-09 (r)

# sub  2 [  2k] |x|=5.780e-05 |dx|=2.330e-12 |r|=1.645e-13 (c)

# all           |x|=5.021e+02 |dx|=1.897e-06 |r|=6.025e-09

# arc           |x|=4.602e+01 |dx|=1.506e-08 |r|=2.619e-10 (λ)


*** Continuation step 31

# SNES iteration  0

# sub  0 [  6k] |x|=5.172e+02 |dx|=1.897e-06 |r|=5.039e+01 (u)

# sub  1 [  6k] |x|=3.378e+00 |dx|=2.740e-08 |r|=9.387e+00 (r)

# sub  2 [  2k] |x|=6.050e-05 |dx|=2.330e-12 |r|=5.637e-13 (c)

# all           |x|=5.172e+02 |dx|=1.897e-06 |r|=5.126e+01

# arc           |x|=8.880e+01 |dx|=0.000e+00 |r|=2.619e-10 (λ)

# SNES iteration  0, KSP iteration   0       |r|=5.126e+01

# SNES iteration  0, KSP iteration   1       |r|=3.328e-10

# SNES iteration  1

# sub  0 [  6k] |x|=5.172e+02 |dx|=1.523e+00 |r|=5.715e-01 (u)

# sub  1 [  6k] |x|=3.376e+00 |dx|=2.440e-02 |r|=5.450e-01 (r)

# sub  2 [  2k] |x|=6.093e-05 |dx|=2.628e-06 |r|=3.308e-10 (c)

# all           |x|=5.173e+02 |dx|=1.523e+00 |r|=7.897e-01

# SNES iteration  1, KSP iteration   0       |r|=1.000e+00

# SNES iteration  1, KSP iteration   1       |r|=1.295e-10

# arc           |x|=8.690e+01 |dx|=1.905e+00 |r|=2.692e-10 (λ)

# SNES iteration  1, KSP iteration   0       |r|=4.466e-01

# SNES iteration  1, KSP iteration   1       |r|=2.694e-11

# SNES iteration  2

# sub  0 [  6k] |x|=5.165e+02 |dx|=9.883e-02 |r|=2.309e-03 (u)

# sub  1 [  6k] |x|=3.372e+00 |dx|=1.437e-03 |r|=5.677e-04 (r)

# sub  2 [  2k] |x|=6.076e-05 |dx|=1.462e-07 |r|=2.682e-11 (c)

# all           |x|=5.165e+02 |dx|=9.884e-02 |r|=2.378e-03

# SNES iteration  2, KSP iteration   0       |r|=1.000e+00

# SNES iteration  2, KSP iteration   1       |r|=1.008e-10

# arc           |x|=8.685e+01 |dx|=4.346e-02 |r|=1.237e-10 (λ)

# SNES iteration  2, KSP iteration   0       |r|=1.946e-03

# SNES iteration  2, KSP iteration   1       |r|=2.700e-13

# SNES iteration  3 success = CONVERGED_FNORM_RELATIVE

# sub  0 [  6k] |x|=5.165e+02 |dx|=2.730e-04 |r|=1.468e-08 (u)

# sub  1 [  6k] |x|=3.372e+00 |dx|=3.960e-06 |r|=7.120e-09 (r)

# sub  2 [  2k] |x|=6.076e-05 |dx|=3.404e-10 |r|=3.760e-13 (c)

# all           |x|=5.165e+02 |dx|=2.730e-04 |r|=1.631e-08

# arc           |x|=8.685e+01 |dx|=8.959e-05 |r|=2.328e-10 (λ)


*** Continuation step 32

# SNES iteration  0

# sub  0 [  6k] |x|=5.312e+02 |dx|=2.730e-04 |r|=5.900e+01 (u)

# sub  1 [  6k] |x|=3.461e+00 |dx|=3.960e-06 |r|=1.033e+01 (r)

# sub  2 [  2k] |x|=6.677e-05 |dx|=3.404e-10 |r|=8.108e-13 (c)

# all           |x|=5.312e+02 |dx|=2.730e-04 |r|=5.990e+01

# arc           |x|=1.277e+02 |dx|=0.000e+00 |r|=2.328e-10 (λ)

# SNES iteration  0, KSP iteration   0       |r|=5.990e+01

# SNES iteration  0, KSP iteration   1       |r|=6.254e-10

# SNES iteration  1

# sub  0 [  6k] |x|=5.315e+02 |dx|=2.056e+00 |r|=9.793e-01 (u)

# sub  1 [  6k] |x|=3.462e+00 |dx|=3.251e-02 |r|=8.433e-01 (r)

# sub  2 [  2k] |x|=6.814e-05 |dx|=3.406e-06 |r|=6.234e-10 (c)

# all           |x|=5.315e+02 |dx|=2.056e+00 |r|=1.292e+00

# SNES iteration  1, KSP iteration   0       |r|=1.000e+00

# SNES iteration  1, KSP iteration   1       |r|=5.745e-10

# arc           |x|=1.247e+02 |dx|=2.967e+00 |r|=5.966e-10 (λ)

# SNES iteration  1, KSP iteration   0       |r|=5.269e-01

# SNES iteration  1, KSP iteration   1       |r|=1.607e-11

# SNES iteration  2

# sub  0 [  6k] |x|=5.303e+02 |dx|=1.064e-01 |r|=2.663e-03 (u)

# sub  1 [  6k] |x|=3.455e+00 |dx|=1.589e-03 |r|=4.793e-04 (r)

# sub  2 [  2k] |x|=6.760e-05 |dx|=1.828e-07 |r|=1.588e-11 (c)

# all           |x|=5.303e+02 |dx|=1.064e-01 |r|=2.706e-03

# SNES iteration  2, KSP iteration   0       |r|=1.000e+00

# SNES iteration  2, KSP iteration   1       |r|=3.590e-10

# arc           |x|=1.246e+02 |dx|=9.912e-02 |r|=1.237e-10 (λ)

# SNES iteration  2, KSP iteration   0       |r|=1.348e-03

# SNES iteration  2, KSP iteration   1       |r|=1.943e-14

# SNES iteration  3 success = CONVERGED_FNORM_RELATIVE

# sub  0 [  6k] |x|=5.303e+02 |dx|=9.357e-05 |r|=3.443e-09 (u)

# sub  1 [  6k] |x|=3.454e+00 |dx|=1.406e-06 |r|=5.377e-09 (r)

# sub  2 [  2k] |x|=6.757e-05 |dx|=1.434e-10 |r|=3.401e-13 (c)

# all           |x|=5.303e+02 |dx|=9.358e-05 |r|=6.385e-09

# arc           |x|=1.246e+02 |dx|=1.434e-04 |r|=1.746e-10 (λ)


*** Continuation step 33

# SNES iteration  0

# sub  0 [  6k] |x|=5.443e+02 |dx|=9.357e-05 |r|=6.997e+01 (u)

# sub  1 [  6k] |x|=3.546e+00 |dx|=1.406e-06 |r|=1.291e+01 (r)

# sub  2 [  2k] |x|=7.741e-05 |dx|=1.434e-10 |r|=8.126e-13 (c)

# all           |x|=5.443e+02 |dx|=9.358e-05 |r|=7.115e+01

# arc           |x|=1.624e+02 |dx|=0.000e+00 |r|=1.746e-10 (λ)

# SNES iteration  0, KSP iteration   0       |r|=7.115e+01

# SNES iteration  0, KSP iteration   1       |r|=4.590e-10

# SNES iteration  1

# sub  0 [  6k] |x|=5.451e+02 |dx|=2.723e+00 |r|=1.661e+00 (u)

# sub  1 [  6k] |x|=3.551e+00 |dx|=4.247e-02 |r|=1.258e+00 (r)

# sub  2 [  2k] |x|=7.986e-05 |dx|=4.228e-06 |r|=4.546e-10 (c)

# all           |x|=5.451e+02 |dx|=2.723e+00 |r|=2.084e+00

# SNES iteration  1, KSP iteration   0       |r|=1.000e+00

# SNES iteration  1, KSP iteration   1       |r|=1.279e-10

# arc           |x|=1.582e+02 |dx|=4.144e+00 |r|=1.819e-10 (λ)

# SNES iteration  1, KSP iteration   0       |r|=5.748e-01

# SNES iteration  1, KSP iteration   1       |r|=3.775e-11

# SNES iteration  2

# sub  0 [  6k] |x|=5.433e+02 |dx|=1.468e-01 |r|=5.126e-03 (u)

# sub  1 [  6k] |x|=3.540e+00 |dx|=2.224e-03 |r|=9.424e-04 (r)

# sub  2 [  2k] |x|=7.862e-05 |dx|=2.518e-07 |r|=3.759e-11 (c)

# all           |x|=5.433e+02 |dx|=1.469e-01 |r|=5.212e-03

# SNES iteration  2, KSP iteration   0       |r|=1.000e+00

# SNES iteration  2, KSP iteration   1       |r|=1.550e-10

# arc           |x|=1.581e+02 |dx|=1.568e-01 |r|=2.547e-10 (λ)

# SNES iteration  2, KSP iteration   0       |r|=1.749e-03

# SNES iteration  2, KSP iteration   1       |r|=9.235e-14

# SNES iteration  3 success = CONVERGED_FNORM_RELATIVE

# sub  0 [  6k] |x|=5.433e+02 |dx|=1.348e-04 |r|=4.991e-09 (u)

# sub  1 [  6k] |x|=3.540e+00 |dx|=2.023e-06 |r|=5.833e-09 (r)

# sub  2 [  2k] |x|=7.855e-05 |dx|=2.030e-10 |r|=4.568e-13 (c)

# all           |x|=5.433e+02 |dx|=1.348e-04 |r|=7.677e-09

# arc           |x|=1.581e+02 |dx|=1.983e-04 |r|=2.401e-10 (λ)


*** Continuation step 34

# SNES iteration  0

# sub  0 [  6k] |x|=5.565e+02 |dx|=1.348e-04 |r|=8.212e+01 (u)

# sub  1 [  6k] |x|=3.636e+00 |dx|=2.023e-06 |r|=1.788e+01 (r)

# sub  2 [  2k] |x|=9.207e-05 |dx|=2.030e-10 |r|=9.477e-13 (c)

# all           |x|=5.566e+02 |dx|=1.348e-04 |r|=8.404e+01

# arc           |x|=1.916e+02 |dx|=0.000e+00 |r|=2.401e-10 (λ)

# SNES iteration  0, KSP iteration   0       |r|=8.404e+01

# SNES iteration  0, KSP iteration   1       |r|=1.177e-09

# SNES iteration  1

# sub  0 [  6k] |x|=5.578e+02 |dx|=3.123e+00 |r|=2.112e+00 (u)

# sub  1 [  6k] |x|=3.645e+00 |dx|=4.797e-02 |r|=1.564e+00 (r)

# sub  2 [  2k] |x|=9.497e-05 |dx|=4.291e-06 |r|=1.175e-09 (c)

# all           |x|=5.579e+02 |dx|=3.124e+00 |r|=2.627e+00

# SNES iteration  1, KSP iteration   0       |r|=1.000e+00

# SNES iteration  1, KSP iteration   1       |r|=2.181e-10

# arc           |x|=1.868e+02 |dx|=4.725e+00 |r|=4.875e-10 (λ)

# SNES iteration  1, KSP iteration   0       |r|=5.052e-01

# SNES iteration  1, KSP iteration   1       |r|=8.217e-11

# SNES iteration  2

# sub  0 [  6k] |x|=5.557e+02 |dx|=1.476e-01 |r|=5.050e-03 (u)

# sub  1 [  6k] |x|=3.630e+00 |dx|=2.256e-03 |r|=1.111e-03 (r)

# sub  2 [  2k] |x|=9.302e-05 |dx|=2.487e-07 |r|=8.194e-11 (c)

# all           |x|=5.557e+02 |dx|=1.476e-01 |r|=5.171e-03

# SNES iteration  2, KSP iteration   0       |r|=1.000e+00

# SNES iteration  2, KSP iteration   1       |r|=3.146e-10

# arc           |x|=1.867e+02 |dx|=1.516e-01 |r|=2.328e-10 (λ)

# SNES iteration  2, KSP iteration   0       |r|=1.266e-03

# SNES iteration  2, KSP iteration   1       |r|=3.308e-14

# SNES iteration  3 success = CONVERGED_FNORM_RELATIVE

# sub  0 [  6k] |x|=5.556e+02 |dx|=1.121e-04 |r|=3.793e-09 (u)

# sub  1 [  6k] |x|=3.630e+00 |dx|=1.676e-06 |r|=5.696e-09 (r)

# sub  2 [  2k] |x|=9.293e-05 |dx|=1.620e-10 |r|=5.617e-13 (c)

# all           |x|=5.556e+02 |dx|=1.121e-04 |r|=6.843e-09

# arc           |x|=1.867e+02 |dx|=1.516e-04 |r|=1.019e-10 (λ)


*** Continuation step 35

# SNES iteration  0

# sub  0 [  6k] |x|=5.682e+02 |dx|=1.121e-04 |r|=9.229e+01 (u)

# sub  1 [  6k] |x|=3.731e+00 |dx|=1.676e-06 |r|=2.511e+01 (r)

# sub  2 [  2k] |x|=1.092e-04 |dx|=1.620e-10 |r|=1.149e-12 (c)

# all           |x|=5.683e+02 |dx|=1.121e-04 |r|=9.564e+01

# arc           |x|=2.153e+02 |dx|=0.000e+00 |r|=1.019e-10 (λ)

# SNES iteration  0, KSP iteration   0       |r|=9.564e+01

# SNES iteration  0, KSP iteration   1       |r|=3.250e-09

# SNES iteration  1

# sub  0 [  6k] |x|=5.698e+02 |dx|=2.548e+00 |r|=1.339e+00 (u)

# sub  1 [  6k] |x|=3.740e+00 |dx|=3.788e-02 |r|=1.255e+00 (r)

# sub  2 [  2k] |x|=1.104e-04 |dx|=2.213e-06 |r|=3.249e-09 (c)

# all           |x|=5.698e+02 |dx|=2.548e+00 |r|=1.835e+00

# SNES iteration  1, KSP iteration   0       |r|=1.000e+00

# SNES iteration  1, KSP iteration   1       |r|=2.794e-10

# arc           |x|=2.117e+02 |dx|=3.620e+00 |r|=3.056e-10 (λ)

# SNES iteration  1, KSP iteration   0       |r|=4.387e-01

# SNES iteration  1, KSP iteration   1       |r|=2.850e-12

# SNES iteration  2

# sub  0 [  6k] |x|=5.679e+02 |dx|=5.849e-02 |r|=4.943e-04 (u)

# sub  1 [  6k] |x|=3.726e+00 |dx|=7.839e-04 |r|=2.036e-04 (r)

# sub  2 [  2k] |x|=1.088e-04 |dx|=6.712e-08 |r|=2.542e-12 (c)

# all           |x|=5.679e+02 |dx|=5.849e-02 |r|=5.346e-04

# SNES iteration  2, KSP iteration   0       |r|=1.000e+00

# SNES iteration  2, KSP iteration   1       |r|=2.048e-10

# arc           |x|=2.116e+02 |dx|=2.053e-02 |r|=2.547e-10 (λ)

# SNES iteration  2, KSP iteration   0       |r|=3.404e-04

# SNES iteration  2, KSP iteration   1       |r|=1.637e-14

# SNES iteration  3 success = CONVERGED_FNORM_RELATIVE

# sub  0 [  6k] |x|=5.679e+02 |dx|=4.985e-05 |r|=2.558e-09 (u)

# sub  1 [  6k] |x|=3.726e+00 |dx|=7.707e-07 |r|=5.613e-09 (r)

# sub  2 [  2k] |x|=1.088e-04 |dx|=6.388e-11 |r|=6.199e-13 (c)

# all           |x|=5.679e+02 |dx|=4.986e-05 |r|=6.168e-09

# arc           |x|=2.116e+02 |dx|=6.252e-05 |r|=2.328e-10 (λ)


*** Continuation step 36

# SNES iteration  0

# sub  0 [  6k] |x|=5.805e+02 |dx|=4.985e-05 |r|=9.594e+01 (u)

# sub  1 [  6k] |x|=3.834e+00 |dx|=7.707e-07 |r|=3.317e+01 (r)

# sub  2 [  2k] |x|=1.258e-04 |dx|=6.388e-11 |r|=1.429e-12 (c)

# all           |x|=5.805e+02 |dx|=4.986e-05 |r|=1.015e+02

# arc           |x|=2.366e+02 |dx|=0.000e+00 |r|=2.328e-10 (λ)

# SNES iteration  0, KSP iteration   0       |r|=1.015e+02

# SNES iteration  0, KSP iteration   1       |r|=3.613e-10

# SNES iteration  1

# sub  0 [  6k] |x|=5.814e+02 |dx|=2.060e+00 |r|=7.353e-01 (u)

# sub  1 [  6k] |x|=3.832e+00 |dx|=2.935e-02 |r|=6.440e-01 (r)

# sub  2 [  2k] |x|=1.234e-04 |dx|=2.549e-06 |r|=3.579e-10 (c)

# all           |x|=5.814e+02 |dx|=2.060e+00 |r|=9.774e-01

# SNES iteration  1, KSP iteration   0       |r|=1.000e+00

# SNES iteration  1, KSP iteration   1       |r|=9.503e-11

# arc           |x|=2.359e+02 |dx|=6.972e-01 |r|=3.129e-10 (λ)

# SNES iteration  1, KSP iteration   0       |r|=9.324e-01

# SNES iteration  1, KSP iteration   1       |r|=2.991e-11

# SNES iteration  2

# sub  0 [  6k] |x|=5.809e+02 |dx|=1.830e-01 |r|=6.124e-03 (u)

# sub  1 [  6k] |x|=3.829e+00 |dx|=2.131e-03 |r|=2.433e-03 (r)

# sub  2 [  2k] |x|=1.229e-04 |dx|=2.510e-07 |r|=2.962e-11 (c)

# all           |x|=5.809e+02 |dx|=1.830e-01 |r|=6.590e-03

# SNES iteration  2, KSP iteration   0       |r|=1.000e+00

# SNES iteration  2, KSP iteration   1       |r|=2.533e-10

# arc           |x|=2.361e+02 |dx|=1.932e-01 |r|=3.056e-10 (λ)

# SNES iteration  2, KSP iteration   0       |r|=2.369e-03

# SNES iteration  2, KSP iteration   1       |r|=9.721e-14

# SNES iteration  3 success = CONVERGED_FNORM_RELATIVE

# sub  0 [  6k] |x|=5.810e+02 |dx|=7.598e-04 |r|=5.143e-08 (u)

# sub  1 [  6k] |x|=3.829e+00 |dx|=7.550e-06 |r|=2.740e-08 (r)

# sub  2 [  2k] |x|=1.230e-04 |dx|=8.317e-10 |r|=7.782e-13 (c)

# all           |x|=5.810e+02 |dx|=7.598e-04 |r|=5.827e-08

# arc           |x|=2.361e+02 |dx|=1.127e-04 |r|=2.547e-10 (λ)


*** Continuation step 37

# SNES iteration  0

# sub  0 [  6k] |x|=5.944e+02 |dx|=7.598e-04 |r|=9.016e+01 (u)

# sub  1 [  6k] |x|=3.943e+00 |dx|=7.550e-06 |r|=3.927e+01 (r)

# sub  2 [  2k] |x|=1.381e-04 |dx|=8.317e-10 |r|=1.602e-12 (c)

# all           |x|=5.944e+02 |dx|=7.598e-04 |r|=9.834e+01

# arc           |x|=2.605e+02 |dx|=0.000e+00 |r|=2.547e-10 (λ)

# SNES iteration  0, KSP iteration   0       |r|=9.834e+01

# SNES iteration  0, KSP iteration   1       |r|=8.966e-10

# SNES iteration  1

# sub  0 [  6k] |x|=5.941e+02 |dx|=3.080e+00 |r|=1.963e+00 (u)

# sub  1 [  6k] |x|=3.928e+00 |dx|=4.569e-02 |r|=1.312e+00 (r)

# sub  2 [  2k] |x|=1.325e-04 |dx|=5.780e-06 |r|=8.934e-10 (c)

# all           |x|=5.941e+02 |dx|=3.081e+00 |r|=2.361e+00

# SNES iteration  1, KSP iteration   0       |r|=1.000e+00

# SNES iteration  1, KSP iteration   1       |r|=2.887e-10

# arc           |x|=2.625e+02 |dx|=1.966e+00 |r|=1.819e-10 (λ)

# SNES iteration  1, KSP iteration   0       |r|=2.057e+00

# SNES iteration  1, KSP iteration   1       |r|=3.885e-11

# SNES iteration  2

# sub  0 [  6k] |x|=5.950e+02 |dx|=2.187e-01 |r|=1.081e-02 (u)

# sub  1 [  6k] |x|=3.933e+00 |dx|=2.891e-03 |r|=4.228e-03 (r)

# sub  2 [  2k] |x|=1.324e-04 |dx|=3.839e-07 |r|=3.852e-11 (c)

# all           |x|=5.950e+02 |dx|=2.187e-01 |r|=1.161e-02

# SNES iteration  2, KSP iteration   0       |r|=1.000e+00

# SNES iteration  2, KSP iteration   1       |r|=6.938e-10

# arc           |x|=2.628e+02 |dx|=2.998e-01 |r|=2.037e-10 (λ)

# SNES iteration  2, KSP iteration   0       |r|=5.165e-03

# SNES iteration  2, KSP iteration   1       |r|=4.417e-14

# SNES iteration  3 success = CONVERGED_FNORM_RELATIVE

# sub  0 [  6k] |x|=5.951e+02 |dx|=4.906e-04 |r|=2.791e-08 (u)

# sub  1 [  6k] |x|=3.934e+00 |dx|=4.639e-06 |r|=1.926e-08 (r)

# sub  2 [  2k] |x|=1.325e-04 |dx|=5.785e-10 |r|=7.811e-13 (c)

# all           |x|=5.952e+02 |dx|=4.906e-04 |r|=3.391e-08

# arc           |x|=2.628e+02 |dx|=4.789e-04 |r|=1.310e-10 (λ)


*** Continuation step 38

# SNES iteration  0

# sub  0 [  6k] |x|=6.094e+02 |dx|=4.906e-04 |r|=8.010e+01 (u)

# sub  1 [  6k] |x|=4.048e+00 |dx|=4.639e-06 |r|=4.005e+01 (r)

# sub  2 [  2k] |x|=1.425e-04 |dx|=5.785e-10 |r|=1.719e-12 (c)

# all           |x|=6.095e+02 |dx|=4.906e-04 |r|=8.956e+01

# arc           |x|=2.895e+02 |dx|=0.000e+00 |r|=1.310e-10 (λ)

# SNES iteration  0, KSP iteration   0       |r|=8.956e+01

# SNES iteration  0, KSP iteration   1       |r|=2.779e-10

# SNES iteration  1

# sub  0 [  6k] |x|=6.084e+02 |dx|=3.360e+00 |r|=2.335e+00 (u)

# sub  1 [  6k] |x|=4.028e+00 |dx|=4.976e-02 |r|=1.505e+00 (r)

# sub  2 [  2k] |x|=1.360e-04 |dx|=6.803e-06 |r|=2.664e-10 (c)

# all           |x|=6.084e+02 |dx|=3.361e+00 |r|=2.778e+00

# SNES iteration  1, KSP iteration   0       |r|=1.000e+00

# SNES iteration  1, KSP iteration   1       |r|=3.218e-10

# arc           |x|=2.918e+02 |dx|=2.285e+00 |r|=2.037e-10 (λ)

# SNES iteration  1, KSP iteration   0       |r|=2.597e+00

# SNES iteration  1, KSP iteration   1       |r|=1.225e-10

# SNES iteration  2

# sub  0 [  6k] |x|=6.094e+02 |dx|=1.683e-01 |r|=5.846e-03 (u)

# sub  1 [  6k] |x|=4.033e+00 |dx|=2.084e-03 |r|=3.159e-03 (r)

# sub  2 [  2k] |x|=1.355e-04 |dx|=3.462e-07 |r|=1.224e-10 (c)

# all           |x|=6.094e+02 |dx|=1.684e-01 |r|=6.644e-03

# SNES iteration  2, KSP iteration   0       |r|=1.000e+00

# SNES iteration  2, KSP iteration   1       |r|=4.827e-10

# arc           |x|=2.920e+02 |dx|=2.156e-01 |r|=1.091e-10 (λ)

# SNES iteration  2, KSP iteration   0       |r|=5.219e-03

# SNES iteration  2, KSP iteration   1       |r|=9.550e-14

# SNES iteration  3 success = CONVERGED_FNORM_RELATIVE

# sub  0 [  6k] |x|=6.095e+02 |dx|=8.497e-04 |r|=1.023e-07 (u)

# sub  1 [  6k] |x|=4.034e+00 |dx|=1.098e-05 |r|=6.088e-08 (r)

# sub  2 [  2k] |x|=1.355e-04 |dx|=1.064e-09 |r|=7.788e-13 (c)

# all           |x|=6.095e+02 |dx|=8.498e-04 |r|=1.190e-07

# arc           |x|=2.920e+02 |dx|=5.217e-04 |r|=2.547e-10 (λ)


*** Continuation step 39

# SNES iteration  0

# sub  0 [  6k] |x|=6.240e+02 |dx|=8.497e-04 |r|=7.890e+01 (u)

# sub  1 [  6k] |x|=4.142e+00 |dx|=1.098e-05 |r|=3.551e+01 (r)

# sub  2 [  2k] |x|=1.388e-04 |dx|=1.064e-09 |r|=1.637e-12 (c)

# all           |x|=6.240e+02 |dx|=8.498e-04 |r|=8.652e+01

# arc           |x|=3.212e+02 |dx|=0.000e+00 |r|=2.547e-10 (λ)

# SNES iteration  0, KSP iteration   0       |r|=8.652e+01

# SNES iteration  0, KSP iteration   1       |r|=4.081e-10

# SNES iteration  1

# sub  0 [  6k] |x|=6.231e+02 |dx|=3.113e+00 |r|=1.910e+00 (u)

# sub  1 [  6k] |x|=4.124e+00 |dx|=4.685e-02 |r|=1.215e+00 (r)

# sub  2 [  2k] |x|=1.323e-04 |dx|=6.844e-06 |r|=4.014e-10 (c)

# all           |x|=6.231e+02 |dx|=3.114e+00 |r|=2.263e+00

# SNES iteration  1, KSP iteration   0       |r|=1.000e+00

# SNES iteration  1, KSP iteration   1       |r|=1.737e-10

# arc           |x|=3.219e+02 |dx|=7.335e-01 |r|=5.748e-10 (λ)

# SNES iteration  1, KSP iteration   0       |r|=2.351e+00

# SNES iteration  1, KSP iteration   1       |r|=3.630e-11

# SNES iteration  2

# sub  0 [  6k] |x|=6.233e+02 |dx|=1.436e-01 |r|=2.406e-03 (u)

# sub  1 [  6k] |x|=4.125e+00 |dx|=1.511e-03 |r|=1.940e-03 (r)

# sub  2 [  2k] |x|=1.318e-04 |dx|=2.651e-07 |r|=3.614e-11 (c)

# all           |x|=6.233e+02 |dx|=1.436e-01 |r|=3.090e-03

# SNES iteration  2, KSP iteration   0       |r|=1.000e+00

# SNES iteration  2, KSP iteration   1       |r|=1.284e-10

# arc           |x|=3.221e+02 |dx|=1.418e-01 |r|=4.366e-11 (λ)

# SNES iteration  2, KSP iteration   0       |r|=5.111e-03

# SNES iteration  2, KSP iteration   1       |r|=4.150e-14

# SNES iteration  3 success = CONVERGED_FNORM_RELATIVE

# sub  0 [  6k] |x|=6.233e+02 |dx|=4.440e-04 |r|=3.304e-08 (u)

# sub  1 [  6k] |x|=4.126e+00 |dx|=5.442e-06 |r|=1.673e-08 (r)

# sub  2 [  2k] |x|=1.318e-04 |dx|=4.101e-10 |r|=8.242e-13 (c)

# all           |x|=6.234e+02 |dx|=4.440e-04 |r|=3.703e-08

# arc           |x|=3.221e+02 |dx|=2.263e-04 |r|=3.201e-10 (λ)


*** Continuation step 40

# SNES iteration  0

# sub  0 [  6k] |x|=6.373e+02 |dx|=4.440e-04 |r|=9.004e+01 (u)

# sub  1 [  6k] |x|=4.227e+00 |dx|=5.442e-06 |r|=2.889e+01 (r)

# sub  2 [  2k] |x|=1.285e-04 |dx|=4.101e-10 |r|=1.734e-12 (c)

# all           |x|=6.374e+02 |dx|=4.440e-04 |r|=9.456e+01

# arc           |x|=3.522e+02 |dx|=0.000e+00 |r|=3.201e-10 (λ)

# SNES iteration  0, KSP iteration   0       |r|=9.456e+01

# SNES iteration  0, KSP iteration   1       |r|=6.382e-10

# SNES iteration  1

# sub  0 [  6k] |x|=6.381e+02 |dx|=3.698e+00 |r|=2.258e+00 (u)

# sub  1 [  6k] |x|=4.219e+00 |dx|=5.212e-02 |r|=1.299e+00 (r)

# sub  2 [  2k] |x|=1.220e-04 |dx|=6.990e-06 |r|=6.318e-10 (c)

# all           |x|=6.381e+02 |dx|=3.698e+00 |r|=2.605e+00

# SNES iteration  1, KSP iteration   0       |r|=1.000e+00

# SNES iteration  1, KSP iteration   1       |r|=2.595e-10

# arc           |x|=3.499e+02 |dx|=2.242e+00 |r|=3.274e-10 (λ)

# SNES iteration  1, KSP iteration   0       |r|=2.022e+00

# SNES iteration  1, KSP iteration   1       |r|=2.501e-11

# SNES iteration  2

# sub  0 [  6k] |x|=6.367e+02 |dx|=1.098e-01 |r|=2.403e-03 (u)

# sub  1 [  6k] |x|=4.211e+00 |dx|=1.387e-03 |r|=1.669e-03 (r)

# sub  2 [  2k] |x|=1.228e-04 |dx|=2.011e-07 |r|=2.486e-11 (c)

# all           |x|=6.367e+02 |dx|=1.098e-01 |r|=2.926e-03

# SNES iteration  2, KSP iteration   0       |r|=1.000e+00

# SNES iteration  2, KSP iteration   1       |r|=5.946e-10

# arc           |x|=3.501e+02 |dx|=1.114e-01 |r|=1.892e-10 (λ)

# SNES iteration  2, KSP iteration   0       |r|=2.244e-03

# SNES iteration  2, KSP iteration   1       |r|=7.013e-14

# SNES iteration  3 success = CONVERGED_FNORM_RELATIVE

# sub  0 [  6k] |x|=6.368e+02 |dx|=1.762e-04 |r|=6.021e-09 (u)

# sub  1 [  6k] |x|=4.211e+00 |dx|=2.101e-06 |r|=7.512e-09 (r)

# sub  2 [  2k] |x|=1.227e-04 |dx|=2.262e-10 |r|=7.629e-13 (c)

# all           |x|=6.368e+02 |dx|=1.762e-04 |r|=9.628e-09

# arc           |x|=3.501e+02 |dx|=1.916e-04 |r|=1.528e-10 (λ)


*** Continuation step 41

# SNES iteration  0

# sub  0 [  6k] |x|=6.504e+02 |dx|=1.762e-04 |r|=9.434e+01 (u)

# sub  1 [  6k] |x|=4.306e+00 |dx|=2.101e-06 |r|=2.560e+01 (r)

# sub  2 [  2k] |x|=1.143e-04 |dx|=2.262e-10 |r|=1.552e-12 (c)

# all           |x|=6.504e+02 |dx|=1.762e-04 |r|=9.776e+01

# arc           |x|=3.780e+02 |dx|=0.000e+00 |r|=1.528e-10 (λ)

# SNES iteration  0, KSP iteration   0       |r|=9.776e+01

# SNES iteration  0, KSP iteration   1       |r|=4.332e-10

# SNES iteration  1

# sub  0 [  6k] |x|=6.577e+02 |dx|=9.986e+00 |r|=6.242e+00 (u)

# sub  1 [  6k] |x|=4.327e+00 |dx|=9.589e-02 |r|=3.645e+00 (r)

# sub  2 [  2k] |x|=1.076e-04 |dx|=7.258e-06 |r|=3.582e-10 (c)

# all           |x|=6.577e+02 |dx|=9.987e+00 |r|=7.229e+00

# SNES iteration  1, KSP iteration   0       |r|=1.000e+00

# SNES iteration  1, KSP iteration   1       |r|=1.529e-10

# arc           |x|=3.714e+02 |dx|=6.613e+00 |r|=3.565e-10 (λ)

# SNES iteration  1, KSP iteration   0       |r|=3.871e+00

# SNES iteration  1, KSP iteration   1       |r|=4.558e-11

# SNES iteration  2

# sub  0 [  6k] |x|=6.511e+02 |dx|=4.465e-01 |r|=2.434e-02 (u)

# sub  1 [  6k] |x|=4.292e+00 |dx|=5.403e-03 |r|=1.293e-02 (r)

# sub  2 [  2k] |x|=1.118e-04 |dx|=2.860e-07 |r|=4.429e-11 (c)

# all           |x|=6.511e+02 |dx|=4.465e-01 |r|=2.756e-02

# SNES iteration  2, KSP iteration   0       |r|=1.000e+00

# SNES iteration  2, KSP iteration   1       |r|=1.887e-10

# arc           |x|=3.711e+02 |dx|=2.479e-01 |r|=7.276e-11 (λ)

# SNES iteration  2, KSP iteration   0       |r|=3.097e-02

# SNES iteration  2, KSP iteration   1       |r|=3.908e-13

# SNES iteration  3

# sub  0 [  6k] |x|=6.509e+02 |dx|=2.066e-03 |r|=1.224e-06 (u)

# sub  1 [  6k] |x|=4.291e+00 |dx|=3.573e-05 |r|=5.882e-07 (r)

# sub  2 [  2k] |x|=1.119e-04 |dx|=1.112e-09 |r|=8.359e-13 (c)

# all           |x|=6.509e+02 |dx|=2.066e-03 |r|=1.358e-06

# SNES iteration  3, KSP iteration   0       |r|=1.000e+00

# SNES iteration  3, KSP iteration   1       |r|=1.683e-10

# arc           |x|=3.711e+02 |dx|=1.048e-03 |r|=6.548e-10 (λ)

# SNES iteration  3, KSP iteration   0       |r|=1.092e-06

# SNES iteration  3, KSP iteration   1       |r|=2.396e-17

# SNES iteration  4 success = CONVERGED_FNORM_RELATIVE

# sub  0 [  6k] |x|=6.509e+02 |dx|=8.267e-08 |r|=2.487e-09 (u)

# sub  1 [  6k] |x|=4.291e+00 |dx|=1.462e-09 |r|=6.522e-09 (r)

# sub  2 [  2k] |x|=1.119e-04 |dx|=4.748e-14 |r|=7.450e-13 (c)

# all           |x|=6.509e+02 |dx|=8.268e-08 |r|=6.980e-09

# arc           |x|=3.711e+02 |dx|=4.148e-08 |r|=3.347e-10 (λ)


*** Continuation step 42

# SNES iteration  0

# sub  0 [  6k] |x|=6.652e+02 |dx|=8.267e-08 |r|=6.682e+01 (u)

# sub  1 [  6k] |x|=4.378e+00 |dx|=1.462e-09 |r|=2.547e+01 (r)

# sub  2 [  2k] |x|=1.017e-04 |dx|=4.748e-14 |r|=1.656e-12 (c)

# all           |x|=6.652e+02 |dx|=8.268e-08 |r|=7.151e+01

# arc           |x|=3.922e+02 |dx|=0.000e+00 |r|=3.347e-10 (λ)

# SNES iteration  0, KSP iteration   0       |r|=7.151e+01

# SNES iteration  0, KSP iteration   1       |r|=2.453e-09

# SNES iteration  1

# sub  0 [  6k] |x|=6.926e+02 |dx|=3.257e+01 |r|=4.380e+01 (u)

# sub  1 [  6k] |x|=4.492e+00 |dx|=2.729e-01 |r|=2.976e+01 (r)

# sub  2 [  2k] |x|=8.850e-05 |dx|=1.484e-05 |r|=2.313e-09 (c)

# all           |x|=6.926e+02 |dx|=3.257e+01 |r|=5.296e+01

# SNES iteration  1, KSP iteration   0       |r|=1.000e+00

# SNES iteration  1, KSP iteration   1       |r|=2.098e-09

# arc           |x|=3.833e+02 |dx|=8.984e+00 |r|=9.459e-11 (λ)

# SNES iteration  1, KSP iteration   0       |r|=5.320e+00

# SNES iteration  1, KSP iteration   1       |r|=1.693e-10

# SNES iteration  2

# sub  0 [  6k] |x|=6.671e+02 |dx|=1.732e+00 |r|=1.424e-01 (u)

# sub  1 [  6k] |x|=4.371e+00 |dx|=1.524e-02 |r|=9.237e-02 (r)

# sub  2 [  2k] |x|=1.015e-04 |dx|=8.501e-07 |r|=1.626e-10 (c)

# all           |x|=6.671e+02 |dx|=1.732e+00 |r|=1.698e-01

# SNES iteration  2, KSP iteration   0       |r|=1.000e+00

# SNES iteration  2, KSP iteration   1       |r|=2.843e-10

# arc           |x|=3.826e+02 |dx|=6.414e-01 |r|=2.838e-10 (λ)

# SNES iteration  2, KSP iteration   0       |r|=4.771e-02

# SNES iteration  2, KSP iteration   1       |r|=1.224e-12

# SNES iteration  3

# sub  0 [  6k] |x|=6.656e+02 |dx|=3.000e-03 |r|=2.918e-06 (u)

# sub  1 [  6k] |x|=4.365e+00 |dx|=5.290e-05 |r|=1.192e-06 (r)

# sub  2 [  2k] |x|=1.023e-04 |dx|=1.552e-09 |r|=1.421e-12 (c)

# all           |x|=6.657e+02 |dx|=3.001e-03 |r|=3.152e-06

# SNES iteration  3, KSP iteration   0       |r|=1.000e+00

# SNES iteration  3, KSP iteration   1       |r|=7.115e-10

# arc           |x|=3.826e+02 |dx|=1.762e-04 |r|=2.328e-10 (λ)

# SNES iteration  3, KSP iteration   0       |r|=3.159e-06

# SNES iteration  3, KSP iteration   1       |r|=5.030e-17

# SNES iteration  4 success = CONVERGED_FNORM_RELATIVE

# sub  0 [  6k] |x|=6.656e+02 |dx|=2.138e-07 |r|=2.486e-09 (u)

# sub  1 [  6k] |x|=4.365e+00 |dx|=3.594e-09 |r|=6.633e-09 (r)

# sub  2 [  2k] |x|=1.023e-04 |dx|=1.036e-13 |r|=6.817e-13 (c)

# all           |x|=6.657e+02 |dx|=2.138e-07 |r|=7.083e-09

# arc           |x|=3.826e+02 |dx|=2.511e-08 |r|=6.112e-10 (λ)


*** Continuation step 43

# SNES iteration  0

# sub  0 [  6k] |x|=6.805e+02 |dx|=2.138e-07 |r|=3.686e+01 (u)

# sub  1 [  6k] |x|=4.442e+00 |dx|=3.594e-09 |r|=2.173e+01 (r)

# sub  2 [  2k] |x|=9.307e-05 |dx|=1.036e-13 |r|=1.516e-12 (c)

# all           |x|=6.805e+02 |dx|=2.138e-07 |r|=4.279e+01

# arc           |x|=3.941e+02 |dx|=0.000e+00 |r|=6.112e-10 (λ)

# SNES iteration  0, KSP iteration   0       |r|=4.279e+01

# SNES iteration  0, KSP iteration   1       |r|=7.479e-09

# SNES iteration  1

# sub  0 [  6k] |x|=6.348e+02 |dx|=5.437e+01 |r|=1.125e+02 (u)

# sub  1 [  6k] |x|=4.254e+00 |dx|=4.492e-01 |r|=8.165e+01 (r)

# sub  2 [  2k] |x|=1.171e-04 |dx|=2.678e-05 |r|=7.352e-09 (c)

# all           |x|=6.349e+02 |dx|=5.437e+01 |r|=1.390e+02

# SNES iteration  1, KSP iteration   0       |r|=1.000e+00

# SNES iteration  1, KSP iteration   1       |r|=8.244e-10

# arc           |x|=3.855e+02 |dx|=8.593e+00 |r|=1.637e-09 (λ)

# SNES iteration  1, KSP iteration   0       |r|=2.773e+00

# SNES iteration  1, KSP iteration   1       |r|=1.332e-09

# SNES iteration  2

# sub  0 [  6k] |x|=6.773e+02 |dx|=3.592e+00 |r|=5.128e-01 (u)

# sub  1 [  6k] |x|=4.418e+00 |dx|=2.977e-02 |r|=3.976e-01 (r)

# sub  2 [  2k] |x|=9.563e-05 |dx|=1.914e-06 |r|=1.316e-09 (c)

# all           |x|=6.774e+02 |dx|=3.592e+00 |r|=6.489e-01

# SNES iteration  2, KSP iteration   0       |r|=1.000e+00

# SNES iteration  2, KSP iteration   1       |r|=3.609e-09

# arc           |x|=3.852e+02 |dx|=3.507e-01 |r|=1.019e-10 (λ)

# SNES iteration  2, KSP iteration   0       |r|=1.005e-02

# SNES iteration  2, KSP iteration   1       |r|=5.588e-13

# SNES iteration  3

# sub  0 [  6k] |x|=6.804e+02 |dx|=2.543e-03 |r|=3.883e-07 (u)

# sub  1 [  6k] |x|=4.432e+00 |dx|=2.222e-05 |r|=2.391e-07 (r)

# sub  2 [  2k] |x|=9.400e-05 |dx|=1.523e-09 |r|=8.452e-13 (c)

# all           |x|=6.804e+02 |dx|=2.543e-03 |r|=4.560e-07

# SNES iteration  3, KSP iteration   0       |r|=1.000e+00

# SNES iteration  3, KSP iteration   1       |r|=3.007e-09

# arc           |x|=3.852e+02 |dx|=2.394e-04 |r|=4.584e-10 (λ)

# SNES iteration  3, KSP iteration   0       |r|=1.444e-07

# SNES iteration  3, KSP iteration   1       |r|=7.460e-18

# SNES iteration  4 success = CONVERGED_FNORM_RELATIVE

# sub  0 [  6k] |x|=6.804e+02 |dx|=4.415e-08 |r|=2.476e-09 (u)

# sub  1 [  6k] |x|=4.432e+00 |dx|=3.745e-10 |r|=6.892e-09 (r)

# sub  2 [  2k] |x|=9.400e-05 |dx|=2.551e-14 |r|=6.034e-13 (c)

# all           |x|=6.804e+02 |dx|=4.416e-08 |r|=7.324e-09

# arc           |x|=3.852e+02 |dx|=4.184e-09 |r|=1.455e-10 (λ)


*** Continuation step 44

# SNES iteration  0

# sub  0 [  6k] |x|=6.952e+02 |dx|=4.415e-08 |r|=2.482e+01 (u)

# sub  1 [  6k] |x|=4.502e+00 |dx|=3.745e-10 |r|=1.989e+01 (r)

# sub  2 [  2k] |x|=8.596e-05 |dx|=2.551e-14 |r|=1.376e-12 (c)

# all           |x|=6.952e+02 |dx|=4.416e-08 |r|=3.180e+01

# arc           |x|=3.877e+02 |dx|=0.000e+00 |r|=1.455e-10 (λ)

# SNES iteration  0, KSP iteration   0       |r|=3.180e+01

# SNES iteration  0, KSP iteration   1       |r|=9.490e-10

# SNES iteration  1

# sub  0 [  6k] |x|=6.839e+02 |dx|=1.328e+01 |r|=6.858e+00 (u)

# sub  1 [  6k] |x|=4.449e+00 |dx|=1.102e-01 |r|=6.441e+00 (r)

# sub  2 [  2k] |x|=9.183e-05 |dx|=6.853e-06 |r|=8.386e-10 (c)

# all           |x|=6.839e+02 |dx|=1.328e+01 |r|=9.408e+00

# SNES iteration  1, KSP iteration   0       |r|=1.000e+00

# SNES iteration  1, KSP iteration   1       |r|=1.350e-10

# arc           |x|=3.795e+02 |dx|=8.173e+00 |r|=7.276e-12 (λ)

# SNES iteration  1, KSP iteration   0       |r|=1.267e+00

# SNES iteration  1, KSP iteration   1       |r|=2.870e-11

# SNES iteration  2

# sub  0 [  6k] |x|=6.947e+02 |dx|=2.202e-01 |r|=2.794e-03 (u)

# sub  1 [  6k] |x|=4.493e+00 |dx|=1.853e-03 |r|=1.888e-03 (r)

# sub  2 [  2k] |x|=8.656e-05 |dx|=1.299e-07 |r|=2.821e-11 (c)

# all           |x|=6.947e+02 |dx|=2.202e-01 |r|=3.372e-03

# SNES iteration  2, KSP iteration   0       |r|=1.000e+00

# SNES iteration  2, KSP iteration   1       |r|=8.444e-11

# arc           |x|=3.794e+02 |dx|=1.223e-01 |r|=2.110e-10 (λ)

# SNES iteration  2, KSP iteration   0       |r|=1.667e-03

# SNES iteration  2, KSP iteration   1       |r|=1.763e-14

# SNES iteration  3 success = CONVERGED_FNORM_RELATIVE

# sub  0 [  6k] |x|=6.948e+02 |dx|=1.339e-04 |r|=4.053e-09 (u)

# sub  1 [  6k] |x|=4.494e+00 |dx|=1.718e-06 |r|=7.132e-09 (r)

# sub  2 [  2k] |x|=8.647e-05 |dx|=8.948e-11 |r|=5.736e-13 (c)

# all           |x|=6.949e+02 |dx|=1.339e-04 |r|=8.203e-09

# arc           |x|=3.794e+02 |dx|=5.568e-05 |r|=1.091e-10 (λ)


*** Continuation step 45

# SNES iteration  0

# sub  0 [  6k] |x|=7.094e+02 |dx|=1.339e-04 |r|=2.470e+01 (u)

# sub  1 [  6k] |x|=4.560e+00 |dx|=1.718e-06 |r|=2.085e+01 (r)

# sub  2 [  2k] |x|=7.912e-05 |dx|=8.948e-11 |r|=1.279e-12 (c)

# all           |x|=7.095e+02 |dx|=1.339e-04 |r|=3.232e+01

# arc           |x|=3.736e+02 |dx|=0.000e+00 |r|=1.091e-10 (λ)

# SNES iteration  0, KSP iteration   0       |r|=3.232e+01

# SNES iteration  0, KSP iteration   1       |r|=2.092e-09

# SNES iteration  1

# sub  0 [  6k] |x|=7.029e+02 |dx|=7.626e+00 |r|=2.234e+00 (u)

# sub  1 [  6k] |x|=4.528e+00 |dx|=6.317e-02 |r|=2.323e+00 (r)

# sub  2 [  2k] |x|=8.231e-05 |dx|=3.809e-06 |r|=2.084e-09 (c)

# all           |x|=7.029e+02 |dx|=7.627e+00 |r|=3.223e+00

# SNES iteration  1, KSP iteration   0       |r|=1.000e+00

# SNES iteration  1, KSP iteration   1       |r|=1.887e-10

# arc           |x|=3.651e+02 |dx|=8.553e+00 |r|=0.000e+00 (λ)

# SNES iteration  1, KSP iteration   0       |r|=7.043e-01

# SNES iteration  1, KSP iteration   1       |r|=7.397e-12

# SNES iteration  2

# sub  0 [  6k] |x|=7.089e+02 |dx|=4.376e-02 |r|=2.818e-04 (u)

# sub  1 [  6k] |x|=4.553e+00 |dx|=5.293e-04 |r|=1.676e-04 (r)

# sub  2 [  2k] |x|=7.928e-05 |dx|=2.613e-08 |r|=7.325e-12 (c)

# all           |x|=7.090e+02 |dx|=4.376e-02 |r|=3.279e-04

# SNES iteration  2, KSP iteration   0       |r|=1.000e+00

# SNES iteration  2, KSP iteration   1       |r|=3.199e-10

# arc           |x|=3.650e+02 |dx|=3.144e-02 |r|=9.459e-11 (λ)

# SNES iteration  2, KSP iteration   0       |r|=3.806e-04

# SNES iteration  2, KSP iteration   1       |r|=1.333e-14

# SNES iteration  3 success = CONVERGED_FNORM_RELATIVE

# sub  0 [  6k] |x|=7.090e+02 |dx|=2.469e-05 |r|=2.712e-09 (u)

# sub  1 [  6k] |x|=4.553e+00 |dx|=3.401e-07 |r|=7.139e-09 (r)

# sub  2 [  2k] |x|=7.927e-05 |dx|=1.408e-11 |r|=5.185e-13 (c)

# all           |x|=7.090e+02 |dx|=2.469e-05 |r|=7.636e-09

# arc           |x|=3.650e+02 |dx|=1.507e-05 |r|=5.093e-11 (λ)


*** Continuation step 46

# SNES iteration  0

# sub  0 [  6k] |x|=7.232e+02 |dx|=2.469e-05 |r|=2.938e+01 (u)

# sub  1 [  6k] |x|=4.617e+00 |dx|=3.401e-07 |r|=2.383e+01 (r)

# sub  2 [  2k] |x|=7.221e-05 |dx|=1.408e-11 |r|=1.131e-12 (c)

# all           |x|=7.232e+02 |dx|=2.469e-05 |r|=3.783e+01

# arc           |x|=3.507e+02 |dx|=0.000e+00 |r|=5.093e-11 (λ)

# SNES iteration  0, KSP iteration   0       |r|=3.783e+01

# SNES iteration  0, KSP iteration   1       |r|=6.423e-10

# SNES iteration  1

# sub  0 [  6k] |x|=7.183e+02 |dx|=5.674e+00 |r|=1.276e+00 (u)

# sub  1 [  6k] |x|=4.592e+00 |dx|=4.738e-02 |r|=1.408e+00 (r)

# sub  2 [  2k] |x|=7.440e-05 |dx|=2.648e-06 |r|=6.261e-10 (c)

# all           |x|=7.183e+02 |dx|=5.674e+00 |r|=1.901e+00

# SNES iteration  1, KSP iteration   0       |r|=1.000e+00

# SNES iteration  1, KSP iteration   1       |r|=4.882e-11

# arc           |x|=3.409e+02 |dx|=9.746e+00 |r|=1.019e-10 (λ)

# SNES iteration  1, KSP iteration   0       |r|=4.779e-01

# SNES iteration  1, KSP iteration   1       |r|=5.896e-12

# SNES iteration  2

# sub  0 [  6k] |x|=7.227e+02 |dx|=1.922e-02 |r|=1.035e-04 (u)

# sub  1 [  6k] |x|=4.611e+00 |dx|=3.223e-04 |r|=5.058e-05 (r)

# sub  2 [  2k] |x|=7.208e-05 |dx|=8.213e-09 |r|=5.907e-12 (c)

# all           |x|=7.228e+02 |dx|=1.922e-02 |r|=1.152e-04

# SNES iteration  2, KSP iteration   0       |r|=1.000e+00

# SNES iteration  2, KSP iteration   1       |r|=7.400e-11

# arc           |x|=3.409e+02 |dx|=3.525e-03 |r|=1.819e-10 (λ)

# SNES iteration  2, KSP iteration   0       |r|=1.079e-04

# SNES iteration  2, KSP iteration   1       |r|=1.038e-15

# SNES iteration  3 success = CONVERGED_FNORM_RELATIVE

# sub  0 [  6k] |x|=7.227e+02 |dx|=5.959e-06 |r|=2.681e-09 (u)

# sub  1 [  6k] |x|=4.611e+00 |dx|=8.111e-08 |r|=7.331e-09 (r)

# sub  2 [  2k] |x|=7.208e-05 |dx|=2.735e-12 |r|=4.513e-13 (c)

# all           |x|=7.227e+02 |dx|=5.959e-06 |r|=7.806e-09

# arc           |x|=3.409e+02 |dx|=4.813e-06 |r|=6.548e-11 (λ)


*** Continuation step 47

# SNES iteration  0

# sub  0 [  6k] |x|=7.367e+02 |dx|=5.959e-06 |r|=3.512e+01 (u)

# sub  1 [  6k] |x|=4.674e+00 |dx|=8.111e-08 |r|=2.846e+01 (r)

# sub  2 [  2k] |x|=6.505e-05 |dx|=2.735e-12 |r|=9.064e-13 (c)

# all           |x|=7.367e+02 |dx|=5.959e-06 |r|=4.520e+01

# arc           |x|=3.168e+02 |dx|=0.000e+00 |r|=6.548e-11 (λ)

# SNES iteration  0, KSP iteration   0       |r|=4.520e+01

# SNES iteration  0, KSP iteration   1       |r|=1.616e-09

# SNES iteration  1

# sub  0 [  6k] |x|=7.326e+02 |dx|=4.777e+00 |r|=1.034e+00 (u)

# sub  1 [  6k] |x|=4.653e+00 |dx|=4.116e-02 |r|=1.190e+00 (r)

# sub  2 [  2k] |x|=6.674e-05 |dx|=2.082e-06 |r|=1.611e-09 (c)

# all           |x|=7.326e+02 |dx|=4.778e+00 |r|=1.577e+00

# SNES iteration  1, KSP iteration   0       |r|=1.000e+00

# SNES iteration  1, KSP iteration   1       |r|=4.519e-11

# arc           |x|=3.051e+02 |dx|=1.170e+01 |r|=1.019e-10 (λ)

# SNES iteration  1, KSP iteration   0       |r|=3.681e-01

# SNES iteration  1, KSP iteration   1       |r|=5.961e-12

# SNES iteration  2

# sub  0 [  6k] |x|=7.362e+02 |dx|=1.198e-02 |r|=4.433e-05 (u)

# sub  1 [  6k] |x|=4.670e+00 |dx|=1.984e-04 |r|=2.049e-05 (r)

# sub  2 [  2k] |x|=6.470e-05 |dx|=4.727e-09 |r|=5.962e-12 (c)

# all           |x|=7.362e+02 |dx|=1.198e-02 |r|=4.883e-05

# SNES iteration  2, KSP iteration   0       |r|=1.000e+00

# SNES iteration  2, KSP iteration   1       |r|=2.948e-11

# arc           |x|=3.052e+02 |dx|=1.414e-02 |r|=2.183e-11 (λ)

# SNES iteration  2, KSP iteration   0       |r|=3.204e-05

# SNES iteration  2, KSP iteration   1       |r|=6.689e-16

# SNES iteration  3 success = CONVERGED_FNORM_RELATIVE

# sub  0 [  6k] |x|=7.362e+02 |dx|=1.343e-06 |r|=2.712e-09 (u)

# sub  1 [  6k] |x|=4.670e+00 |dx|=1.757e-08 |r|=7.445e-09 (r)

# sub  2 [  2k] |x|=6.471e-05 |dx|=5.239e-13 |r|=3.957e-13 (c)

# all           |x|=7.362e+02 |dx|=1.343e-06 |r|=7.924e-09

# arc           |x|=3.052e+02 |dx|=1.399e-06 |r|=1.601e-10 (λ)


*** Continuation step 48

# SNES iteration  0

# sub  0 [  6k] |x|=7.499e+02 |dx|=1.343e-06 |r|=4.039e+01 (u)

# sub  1 [  6k] |x|=4.734e+00 |dx|=1.757e-08 |r|=3.462e+01 (r)

# sub  2 [  2k] |x|=5.752e-05 |dx|=5.239e-13 |r|=8.391e-13 (c)

# all           |x|=7.499e+02 |dx|=1.343e-06 |r|=5.320e+01

# arc           |x|=2.694e+02 |dx|=0.000e+00 |r|=1.601e-10 (λ)

# SNES iteration  0, KSP iteration   0       |r|=5.320e+01

# SNES iteration  0, KSP iteration   1       |r|=4.113e-10

# SNES iteration  1

# sub  0 [  6k] |x|=7.463e+02 |dx|=4.264e+00 |r|=1.006e+00 (u)

# sub  1 [  6k] |x|=4.716e+00 |dx|=3.872e-02 |r|=1.208e+00 (r)

# sub  2 [  2k] |x|=5.893e-05 |dx|=1.776e-06 |r|=3.981e-10 (c)

# all           |x|=7.463e+02 |dx|=4.265e+00 |r|=1.572e+00

# SNES iteration  1, KSP iteration   0       |r|=1.000e+00

# SNES iteration  1, KSP iteration   1       |r|=2.924e-11

# arc           |x|=2.551e+02 |dx|=1.429e+01 |r|=3.711e-10 (λ)

# SNES iteration  1, KSP iteration   0       |r|=2.997e-01

# SNES iteration  1, KSP iteration   1       |r|=1.739e-12

# SNES iteration  2

# sub  0 [  6k] |x|=7.496e+02 |dx|=6.074e-03 |r|=1.472e-05 (u)

# sub  1 [  6k] |x|=4.731e+00 |dx|=9.987e-05 |r|=8.565e-06 (r)

# sub  2 [  2k] |x|=5.703e-05 |dx|=2.763e-09 |r|=1.759e-12 (c)

# all           |x|=7.496e+02 |dx|=6.075e-03 |r|=1.703e-05

# SNES iteration  2, KSP iteration   0       |r|=1.000e+00

# SNES iteration  2, KSP iteration   1       |r|=5.738e-11

# arc           |x|=2.551e+02 |dx|=1.082e-02 |r|=7.276e-11 (λ)

# SNES iteration  2, KSP iteration   0       |r|=1.062e-05

# SNES iteration  2, KSP iteration   1       |r|=1.052e-16

# SNES iteration  3 success = CONVERGED_FNORM_RELATIVE

# sub  0 [  6k] |x|=7.496e+02 |dx|=2.298e-07 |r|=2.696e-09 (u)

# sub  1 [  6k] |x|=4.731e+00 |dx|=2.893e-09 |r|=7.236e-09 (r)

# sub  2 [  2k] |x|=5.703e-05 |dx|=1.144e-13 |r|=3.439e-13 (c)

# all           |x|=7.496e+02 |dx|=2.298e-07 |r|=7.722e-09

# arc           |x|=2.551e+02 |dx|=2.655e-07 |r|=4.366e-11 (λ)


*** Continuation step 49

# SNES iteration  0

# sub  0 [  6k] |x|=7.631e+02 |dx|=2.298e-07 |r|=4.446e+01 (u)

# sub  1 [  6k] |x|=4.798e+00 |dx|=2.893e-09 |r|=4.205e+01 (r)

# sub  2 [  2k] |x|=4.968e-05 |dx|=1.144e-13 |r|=6.903e-13 (c)

# all           |x|=7.631e+02 |dx|=2.298e-07 |r|=6.119e+01

# arc           |x|=2.050e+02 |dx|=0.000e+00 |r|=4.366e-11 (λ)

# SNES iteration  0, KSP iteration   0       |r|=6.119e+01

# SNES iteration  0, KSP iteration   1       |r|=2.147e-10

# SNES iteration  1

# sub  0 [  6k] |x|=7.599e+02 |dx|=3.892e+00 |r|=1.027e+00 (u)

# sub  1 [  6k] |x|=4.782e+00 |dx|=3.754e-02 |r|=1.296e+00 (r)

# sub  2 [  2k] |x|=5.091e-05 |dx|=1.606e-06 |r|=1.921e-10 (c)

# all           |x|=7.599e+02 |dx|=3.892e+00 |r|=1.654e+00

# SNES iteration  1, KSP iteration   0       |r|=1.000e+00

# SNES iteration  1, KSP iteration   1       |r|=2.253e-11

# arc           |x|=1.878e+02 |dx|=1.725e+01 |r|=5.239e-10 (λ)

# SNES iteration  1, KSP iteration   0       |r|=2.465e-01

# SNES iteration  1, KSP iteration   1       |r|=2.409e-13

# SNES iteration  2

# sub  0 [  6k] |x|=7.629e+02 |dx|=2.654e-03 |r|=4.948e-06 (u)

# sub  1 [  6k] |x|=4.797e+00 |dx|=5.081e-05 |r|=3.216e-06 (r)

# sub  2 [  2k] |x|=4.915e-05 |dx|=2.532e-09 |r|=3.591e-13 (c)

# all           |x|=7.629e+02 |dx|=2.655e-03 |r|=5.901e-06

# SNES iteration  2, KSP iteration   0       |r|=1.000e+00

# SNES iteration  2, KSP iteration   1       |r|=1.629e-11

# arc           |x|=1.878e+02 |dx|=1.405e-03 |r|=6.548e-11 (λ)

# SNES iteration  2, KSP iteration   0       |r|=6.165e-06

# SNES iteration  2, KSP iteration   1       |r|=8.227e-17

# SNES iteration  3 success = CONVERGED_FNORM_RELATIVE

# sub  0 [  6k] |x|=7.629e+02 |dx|=2.838e-07 |r|=2.771e-09 (u)

# sub  1 [  6k] |x|=4.797e+00 |dx|=3.735e-09 |r|=7.597e-09 (r)

# sub  2 [  2k] |x|=4.915e-05 |dx|=1.148e-13 |r|=2.715e-13 (c)

# all           |x|=7.629e+02 |dx|=2.838e-07 |r|=8.086e-09

# arc           |x|=1.878e+02 |dx|=6.441e-08 |r|=9.240e-10 (λ)


*** Continuation step 50

# SNES iteration  0

# sub  0 [  6k] |x|=7.764e+02 |dx|=2.838e-07 |r|=4.701e+01 (u)

# sub  1 [  6k] |x|=4.869e+00 |dx|=3.735e-09 |r|=5.028e+01 (r)

# sub  2 [  2k] |x|=4.197e-05 |dx|=1.148e-13 |r|=5.922e-13 (c)

# all           |x|=7.764e+02 |dx|=2.838e-07 |r|=6.883e+01

# arc           |x|=1.205e+02 |dx|=0.000e+00 |r|=9.240e-10 (λ)

# SNES iteration  0, KSP iteration   0       |r|=6.883e+01

# SNES iteration  0, KSP iteration   1       |r|=6.529e-10

# SNES iteration  1

# sub  0 [  6k] |x|=7.736e+02 |dx|=3.567e+00 |r|=1.019e+00 (u)

# sub  1 [  6k] |x|=4.854e+00 |dx|=3.634e-02 |r|=1.368e+00 (r)

# sub  2 [  2k] |x|=4.306e-05 |dx|=1.508e-06 |r|=6.470e-10 (c)

# all           |x|=7.736e+02 |dx|=3.567e+00 |r|=1.706e+00

# SNES iteration  1, KSP iteration   0       |r|=1.000e+00

# SNES iteration  1, KSP iteration   1       |r|=1.041e-11

# arc           |x|=1.002e+02 |dx|=2.023e+01 |r|=1.892e-10 (λ)

# SNES iteration  1, KSP iteration   0       |r|=2.028e-01

# SNES iteration  1, KSP iteration   1       |r|=4.823e-13

# SNES iteration  2

# sub  0 [  6k] |x|=7.764e+02 |dx|=3.560e-03 |r|=6.769e-06 (u)

# sub  1 [  6k] |x|=4.869e+00 |dx|=5.686e-05 |r|=1.834e-06 (r)

# sub  2 [  2k] |x|=4.158e-05 |dx|=2.832e-09 |r|=5.344e-13 (c)

# all           |x|=7.764e+02 |dx|=3.560e-03 |r|=7.013e-06

# SNES iteration  2, KSP iteration   0       |r|=1.000e+00

# SNES iteration  2, KSP iteration   1       |r|=2.054e-11

# arc           |x|=1.002e+02 |dx|=1.678e-02 |r|=3.711e-10 (λ)

# SNES iteration  2, KSP iteration   0       |r|=5.964e-06

# SNES iteration  2, KSP iteration   1       |r|=1.293e-17

# SNES iteration  3 success = CONVERGED_FNORM_RELATIVE

# sub  0 [  6k] |x|=7.764e+02 |dx|=6.282e-08 |r|=2.813e-09 (u)

# sub  1 [  6k] |x|=4.869e+00 |dx|=1.125e-09 |r|=7.687e-09 (r)

# sub  2 [  2k] |x|=4.158e-05 |dx|=4.998e-14 |r|=2.214e-13 (c)

# all           |x|=7.764e+02 |dx|=6.283e-08 |r|=8.185e-09

# arc           |x|=1.002e+02 |dx|=1.377e-07 |r|=1.892e-10 (λ)


*** Continuation step 51

# SNES iteration  0

# sub  0 [  6k] |x|=7.900e+02 |dx|=6.282e-08 |r|=4.817e+01 (u)

# sub  1 [  6k] |x|=4.947e+00 |dx|=1.125e-09 |r|=5.871e+01 (r)

# sub  2 [  2k] |x|=3.562e-05 |dx|=4.998e-14 |r|=4.799e-13 (c)

# all           |x|=7.901e+02 |dx|=6.283e-08 |r|=7.594e+01

# arc           |x|=1.265e+01 |dx|=0.000e+00 |r|=1.892e-10 (λ)

# SNES iteration  0, KSP iteration   0       |r|=7.594e+01

# SNES iteration  0, KSP iteration   1       |r|=1.715e-10

# SNES iteration  1

# sub  0 [  6k] |x|=7.876e+02 |dx|=3.260e+00 |r|=9.576e-01 (u)

# sub  1 [  6k] |x|=4.933e+00 |dx|=3.455e-02 |r|=1.382e+00 (r)

# sub  2 [  2k] |x|=3.652e-05 |dx|=1.442e-06 |r|=1.516e-10 (c)

# all           |x|=7.876e+02 |dx|=3.260e+00 |r|=1.681e+00

# SNES iteration  1, KSP iteration   0       |r|=1.000e+00

# SNES iteration  1, KSP iteration   1       |r|=1.329e-11

# arc           |x|=1.028e+01 |dx|=2.294e+01 |r|=1.455e-10 (λ)

# SNES iteration  1, KSP iteration   0       |r|=1.717e-01

# SNES iteration  1, KSP iteration   1       |r|=1.212e-12

# SNES iteration  2

# sub  0 [  6k] |x|=7.901e+02 |dx|=4.508e-03 |r|=7.762e-06 (u)

# sub  1 [  6k] |x|=4.948e+00 |dx|=6.544e-05 |r|=2.501e-06 (r)

# sub  2 [  2k] |x|=3.564e-05 |dx|=2.584e-09 |r|=1.219e-12 (c)

# all           |x|=7.902e+02 |dx|=4.509e-03 |r|=8.155e-06

# SNES iteration  2, KSP iteration   0       |r|=1.000e+00

# SNES iteration  2, KSP iteration   1       |r|=6.000e-12

# arc           |x|=1.031e+01 |dx|=2.895e-02 |r|=1.673e-10 (λ)

# SNES iteration  2, KSP iteration   0       |r|=5.423e-06

# SNES iteration  2, KSP iteration   1       |r|=4.517e-18

# SNES iteration  3 success = CONVERGED_FNORM_RELATIVE

# sub  0 [  6k] |x|=7.902e+02 |dx|=4.885e-08 |r|=2.795e-09 (u)

# sub  1 [  6k] |x|=4.948e+00 |dx|=9.981e-10 |r|=7.977e-09 (r)

# sub  2 [  2k] |x|=3.564e-05 |dx|=5.260e-14 |r|=1.764e-13 (c)

# all           |x|=7.902e+02 |dx|=4.886e-08 |r|=8.453e-09

# arc           |x|=1.031e+01 |dx|=2.155e-07 |r|=3.420e-10 (λ)


*** Continuation step 52

# SNES iteration  0

# sub  0 [  6k] |x|=8.040e+02 |dx|=4.885e-08 |r|=4.836e+01 (u)

# sub  1 [  6k] |x|=5.032e+00 |dx|=9.981e-10 |r|=6.681e+01 (r)

# sub  2 [  2k] |x|=3.313e-05 |dx|=5.260e-14 |r|=3.806e-13 (c)

# all           |x|=8.041e+02 |dx|=4.886e-08 |r|=8.248e+01

# arc           |x|=1.208e+02 |dx|=0.000e+00 |r|=3.420e-10 (λ)

# SNES iteration  0, KSP iteration   0       |r|=8.248e+01

# SNES iteration  0, KSP iteration   1       |r|=2.472e-10

# SNES iteration  1

# sub  0 [  6k] |x|=8.018e+02 |dx|=2.970e+00 |r|=8.546e-01 (u)

# sub  1 [  6k] |x|=5.018e+00 |dx|=3.221e-02 |r|=1.335e+00 (r)

# sub  2 [  2k] |x|=3.369e-05 |dx|=1.384e-06 |r|=2.367e-10 (c)

# all           |x|=8.018e+02 |dx|=2.970e+00 |r|=1.585e+00

# SNES iteration  1, KSP iteration   0       |r|=1.000e+00

# SNES iteration  1, KSP iteration   1       |r|=1.833e-11

# arc           |x|=1.461e+02 |dx|=2.525e+01 |r|=1.237e-10 (λ)

# SNES iteration  1, KSP iteration   0       |r|=1.520e-01

# SNES iteration  1, KSP iteration   1       |r|=3.177e-13

# SNES iteration  2

# sub  0 [  6k] |x|=8.042e+02 |dx|=4.391e-03 |r|=6.457e-06 (u)

# sub  1 [  6k] |x|=5.033e+00 |dx|=6.310e-05 |r|=2.935e-06 (r)

# sub  2 [  2k] |x|=3.381e-05 |dx|=2.216e-09 |r|=3.359e-13 (c)

# all           |x|=8.042e+02 |dx|=4.391e-03 |r|=7.093e-06

# SNES iteration  2, KSP iteration   0       |r|=1.000e+00

# SNES iteration  2, KSP iteration   1       |r|=6.865e-12

# arc           |x|=1.461e+02 |dx|=3.432e-02 |r|=1.819e-10 (λ)

# SNES iteration  2, KSP iteration   0       |r|=4.154e-06

# SNES iteration  2, KSP iteration   1       |r|=2.429e-18

# SNES iteration  3 success = CONVERGED_FNORM_RELATIVE

# sub  0 [  6k] |x|=8.042e+02 |dx|=4.359e-08 |r|=2.862e-09 (u)

# sub  1 [  6k] |x|=5.033e+00 |dx|=8.437e-10 |r|=8.702e-09 (r)

# sub  2 [  2k] |x|=3.381e-05 |dx|=4.882e-14 |r|=1.567e-13 (c)

# all           |x|=8.042e+02 |dx|=4.360e-08 |r|=9.161e-09

# arc           |x|=1.461e+02 |dx|=2.420e-07 |r|=2.838e-10 (λ)


*** Continuation step 53

# SNES iteration  0

# sub  0 [  6k] |x|=8.184e+02 |dx|=4.359e-08 |r|=4.806e+01 (u)

# sub  1 [  6k] |x|=5.123e+00 |dx|=8.437e-10 |r|=7.427e+01 (r)

# sub  2 [  2k] |x|=3.713e-05 |dx|=4.882e-14 |r|=3.337e-13 (c)

# all           |x|=8.184e+02 |dx|=4.360e-08 |r|=8.847e+01

# arc           |x|=2.819e+02 |dx|=0.000e+00 |r|=2.838e-10 (λ)

# SNES iteration  0, KSP iteration   0       |r|=8.847e+01

# SNES iteration  0, KSP iteration   1       |r|=9.208e-11

# SNES iteration  1

# sub  0 [  6k] |x|=8.164e+02 |dx|=2.702e+00 |r|=7.342e-01 (u)

# sub  1 [  6k] |x|=5.110e+00 |dx|=2.955e-02 |r|=1.246e+00 (r)

# sub  2 [  2k] |x|=3.727e-05 |dx|=1.321e-06 |r|=6.215e-11 (c)

# all           |x|=8.164e+02 |dx|=2.702e+00 |r|=1.446e+00

# SNES iteration  1, KSP iteration   0       |r|=1.000e+00

# SNES iteration  1, KSP iteration   1       |r|=2.945e-11

# arc           |x|=3.091e+02 |dx|=2.718e+01 |r|=5.093e-11 (λ)

# SNES iteration  1, KSP iteration   0       |r|=1.363e-01

# SNES iteration  1, KSP iteration   1       |r|=4.698e-13

# SNES iteration  2

# sub  0 [  6k] |x|=8.186e+02 |dx|=3.682e-03 |r|=4.451e-06 (u)

# sub  1 [  6k] |x|=5.125e+00 |dx|=5.432e-05 |r|=2.687e-06 (r)

# sub  2 [  2k] |x|=3.844e-05 |dx|=1.961e-09 |r|=4.735e-13 (c)

# all           |x|=8.186e+02 |dx|=3.683e-03 |r|=5.199e-06

# SNES iteration  2, KSP iteration   0       |r|=1.000e+00

# SNES iteration  2, KSP iteration   1       |r|=7.936e-12

# arc           |x|=3.091e+02 |dx|=3.303e-02 |r|=4.293e-10 (λ)

# SNES iteration  2, KSP iteration   0       |r|=2.842e-06

# SNES iteration  2, KSP iteration   1       |r|=1.197e-17

# SNES iteration  3 success = CONVERGED_FNORM_RELATIVE

# sub  0 [  6k] |x|=8.186e+02 |dx|=3.381e-08 |r|=2.874e-09 (u)

# sub  1 [  6k] |x|=5.125e+00 |dx|=6.130e-10 |r|=8.830e-09 (r)

# sub  2 [  2k] |x|=3.845e-05 |dx|=3.831e-14 |r|=1.364e-13 (c)

# all           |x|=8.186e+02 |dx|=3.382e-08 |r|=9.286e-09

# arc           |x|=3.091e+02 |dx|=2.321e-07 |r|=1.237e-10 (λ)


*** Continuation step 54

# SNES iteration  0

# sub  0 [  6k] |x|=8.331e+02 |dx|=3.381e-08 |r|=4.762e+01 (u)

# sub  1 [  6k] |x|=5.220e+00 |dx|=6.130e-10 |r|=8.103e+01 (r)

# sub  2 [  2k] |x|=4.764e-05 |dx|=3.831e-14 |r|=3.267e-13 (c)

# all           |x|=8.331e+02 |dx|=3.382e-08 |r|=9.399e+01

# arc           |x|=4.722e+02 |dx|=0.000e+00 |r|=1.237e-10 (λ)

# SNES iteration  0, KSP iteration   0       |r|=9.399e+01

# SNES iteration  0, KSP iteration   1       |r|=1.443e-10

# SNES iteration  1

# sub  0 [  6k] |x|=8.312e+02 |dx|=2.463e+00 |r|=6.179e-01 (u)

# sub  1 [  6k] |x|=5.208e+00 |dx|=2.687e-02 |r|=1.138e+00 (r)

# sub  2 [  2k] |x|=4.748e-05 |dx|=1.254e-06 |r|=1.280e-10 (c)

# all           |x|=8.313e+02 |dx|=2.463e+00 |r|=1.295e+00

# SNES iteration  1, KSP iteration   0       |r|=1.000e+00

# SNES iteration  1, KSP iteration   1       |r|=1.935e-11

# arc           |x|=5.010e+02 |dx|=2.884e+01 |r|=1.382e-10 (λ)

# SNES iteration  1, KSP iteration   0       |r|=1.212e-01

# SNES iteration  1, KSP iteration   1       |r|=2.679e-13

# SNES iteration  2

# sub  0 [  6k] |x|=8.333e+02 |dx|=2.871e-03 |r|=2.767e-06 (u)

# sub  1 [  6k] |x|=5.221e+00 |dx|=4.384e-05 |r|=2.083e-06 (r)

# sub  2 [  2k] |x|=4.932e-05 |dx|=1.722e-09 |r|=3.040e-13 (c)

# all           |x|=8.333e+02 |dx|=2.871e-03 |r|=3.464e-06

# SNES iteration  2, KSP iteration   0       |r|=1.000e+00

# SNES iteration  2, KSP iteration   1       |r|=2.504e-11

# arc           |x|=5.010e+02 |dx|=2.786e-02 |r|=1.528e-10 (λ)

# SNES iteration  2, KSP iteration   0       |r|=1.863e-06

# SNES iteration  2, KSP iteration   1       |r|=7.727e-18

# SNES iteration  3 success = CONVERGED_FNORM_RELATIVE

# sub  0 [  6k] |x|=8.333e+02 |dx|=2.382e-08 |r|=2.831e-09 (u)

# sub  1 [  6k] |x|=5.222e+00 |dx|=4.042e-10 |r|=9.071e-09 (r)

# sub  2 [  2k] |x|=4.933e-05 |dx|=2.721e-14 |r|=1.593e-13 (c)

# all           |x|=8.333e+02 |dx|=2.382e-08 |r|=9.502e-09

# arc           |x|=5.010e+02 |dx|=1.980e-07 |r|=2.910e-11 (λ)


*** Continuation step 55

# SNES iteration  0

# sub  0 [  6k] |x|=8.481e+02 |dx|=2.382e-08 |r|=4.720e+01 (u)

# sub  1 [  6k] |x|=5.322e+00 |dx|=4.042e-10 |r|=8.715e+01 (r)

# sub  2 [  2k] |x|=6.292e-05 |dx|=2.721e-14 |r|=3.360e-13 (c)

# all           |x|=8.481e+02 |dx|=2.382e-08 |r|=9.911e+01

# arc           |x|=6.929e+02 |dx|=0.000e+00 |r|=2.910e-11 (λ)

# SNES iteration  0, KSP iteration   0       |r|=9.911e+01

# SNES iteration  0, KSP iteration   1       |r|=2.191e-10

# SNES iteration  1

# sub  0 [  6k] |x|=8.464e+02 |dx|=2.255e+00 |r|=5.167e-01 (u)

# sub  1 [  6k] |x|=5.310e+00 |dx|=2.437e-02 |r|=1.028e+00 (r)

# sub  2 [  2k] |x|=6.259e-05 |dx|=1.184e-06 |r|=2.125e-10 (c)

# all           |x|=8.464e+02 |dx|=2.255e+00 |r|=1.150e+00

# SNES iteration  1, KSP iteration   0       |r|=1.000e+00

# SNES iteration  1, KSP iteration   1       |r|=1.152e-11

# arc           |x|=7.233e+02 |dx|=3.033e+01 |r|=3.783e-10 (λ)

# SNES iteration  1, KSP iteration   0       |r|=1.065e-01

# SNES iteration  1, KSP iteration   1       |r|=2.824e-13

# SNES iteration  2

# sub  0 [  6k] |x|=8.483e+02 |dx|=2.192e-03 |r|=1.642e-06 (u)

# sub  1 [  6k] |x|=5.323e+00 |dx|=3.425e-05 |r|=1.466e-06 (r)

# sub  2 [  2k] |x|=6.479e-05 |dx|=1.471e-09 |r|=3.311e-13 (c)

# all           |x|=8.483e+02 |dx|=2.192e-03 |r|=2.201e-06

# SNES iteration  2, KSP iteration   0       |r|=1.000e+00

# SNES iteration  2, KSP iteration   1       |r|=2.028e-11

# arc           |x|=7.233e+02 |dx|=2.160e-02 |r|=8.731e-11 (λ)

# SNES iteration  2, KSP iteration   0       |r|=1.228e-06

# SNES iteration  2, KSP iteration   1       |r|=2.239e-18

# SNES iteration  3 success = CONVERGED_FNORM_RELATIVE

# sub  0 [  6k] |x|=8.483e+02 |dx|=1.559e-08 |r|=2.861e-09 (u)

# sub  1 [  6k] |x|=5.323e+00 |dx|=2.514e-10 |r|=8.862e-09 (r)

# sub  2 [  2k] |x|=6.479e-05 |dx|=1.802e-14 |r|=1.836e-13 (c)

# all           |x|=8.483e+02 |dx|=1.560e-08 |r|=9.313e-09

# arc           |x|=7.233e+02 |dx|=1.518e-07 |r|=2.183e-10 (λ)


*** Continuation step 56

# SNES iteration  0

# sub  0 [  6k] |x|=8.634e+02 |dx|=1.559e-08 |r|=4.688e+01 (u)

# sub  1 [  6k] |x|=5.428e+00 |dx|=2.514e-10 |r|=9.275e+01 (r)

# sub  2 [  2k] |x|=8.166e-05 |dx|=1.802e-14 |r|=4.141e-13 (c)

# all           |x|=8.634e+02 |dx|=1.560e-08 |r|=1.039e+02

# arc           |x|=9.455e+02 |dx|=0.000e+00 |r|=2.183e-10 (λ)

# SNES iteration  0, KSP iteration   0       |r|=1.039e+02

# SNES iteration  0, KSP iteration   1       |r|=1.059e-10

# SNES iteration  1

# sub  0 [  6k] |x|=8.617e+02 |dx|=2.075e+00 |r|=4.338e-01 (u)

# sub  1 [  6k] |x|=5.416e+00 |dx|=2.215e-02 |r|=9.255e-01 (r)

# sub  2 [  2k] |x|=8.125e-05 |dx|=1.117e-06 |r|=9.081e-11 (c)

# all           |x|=8.617e+02 |dx|=2.075e+00 |r|=1.022e+00

# SNES iteration  1, KSP iteration   0       |r|=1.000e+00

# SNES iteration  1, KSP iteration   1       |r|=2.138e-11

# arc           |x|=9.773e+02 |dx|=3.173e+01 |r|=1.019e-10 (λ)

# SNES iteration  1, KSP iteration   0       |r|=9.308e-02

# SNES iteration  1, KSP iteration   1       |r|=1.239e-13

# SNES iteration  2

# sub  0 [  6k] |x|=8.635e+02 |dx|=1.689e-03 |r|=9.685e-07 (u)

# sub  1 [  6k] |x|=5.428e+00 |dx|=2.647e-05 |r|=9.816e-07 (r)

# sub  2 [  2k] |x|=8.363e-05 |dx|=1.255e-09 |r|=2.604e-13 (c)

# all           |x|=8.635e+02 |dx|=1.689e-03 |r|=1.379e-06

# SNES iteration  2, KSP iteration   0       |r|=1.000e+00

# SNES iteration  2, KSP iteration   1       |r|=8.012e-12

# arc           |x|=9.773e+02 |dx|=1.590e-02 |r|=9.459e-11 (λ)

# SNES iteration  2, KSP iteration   0       |r|=8.288e-07

# SNES iteration  2, KSP iteration   1       |r|=6.908e-19

# SNES iteration  3 success = CONVERGED_FNORM_RELATIVE

# sub  0 [  6k] |x|=8.635e+02 |dx|=9.677e-09 |r|=2.835e-09 (u)

# sub  1 [  6k] |x|=5.428e+00 |dx|=1.517e-10 |r|=9.054e-09 (r)

# sub  2 [  2k] |x|=8.363e-05 |dx|=1.131e-14 |r|=2.425e-13 (c)

# all           |x|=8.635e+02 |dx|=9.679e-09 |r|=9.488e-09

# arc           |x|=9.773e+02 |dx|=1.071e-07 |r|=3.929e-10 (λ)


*** Continuation step 57

# SNES iteration  0

# sub  0 [  6k] |x|=8.788e+02 |dx|=9.677e-09 |r|=4.664e+01 (u)

# sub  1 [  6k] |x|=5.537e+00 |dx|=1.517e-10 |r|=9.796e+01 (r)

# sub  2 [  2k] |x|=1.032e-04 |dx|=1.131e-14 |r|=4.912e-13 (c)

# all           |x|=8.789e+02 |dx|=9.679e-09 |r|=1.085e+02

# arc           |x|=1.231e+03 |dx|=0.000e+00 |r|=3.929e-10 (λ)

# SNES iteration  0, KSP iteration   0       |r|=1.085e+02

# SNES iteration  0, KSP iteration   1       |r|=2.634e-10

# SNES iteration  1

# sub  0 [  6k] |x|=8.773e+02 |dx|=1.922e+00 |r|=3.678e-01 (u)

# sub  1 [  6k] |x|=5.525e+00 |dx|=2.021e-02 |r|=8.347e-01 (r)

# sub  2 [  2k] |x|=1.028e-04 |dx|=1.054e-06 |r|=2.594e-10 (c)

# all           |x|=8.773e+02 |dx|=1.922e+00 |r|=9.122e-01

# SNES iteration  1, KSP iteration   0       |r|=1.000e+00

# SNES iteration  1, KSP iteration   1       |r|=8.193e-12

# arc           |x|=1.264e+03 |dx|=3.309e+01 |r|=2.910e-11 (λ)

# SNES iteration  1, KSP iteration   0       |r|=8.131e-02

# SNES iteration  1, KSP iteration   1       |r|=2.782e-13

# SNES iteration  2

# sub  0 [  6k] |x|=8.790e+02 |dx|=1.328e-03 |r|=5.843e-07 (u)

# sub  1 [  6k] |x|=5.537e+00 |dx|=2.052e-05 |r|=6.473e-07 (r)

# sub  2 [  2k] |x|=1.053e-04 |dx|=1.104e-09 |r|=4.003e-13 (c)

# all           |x|=8.790e+02 |dx|=1.328e-03 |r|=8.720e-07

# SNES iteration  2, KSP iteration   0       |r|=1.000e+00

# SNES iteration  2, KSP iteration   1       |r|=4.294e-12

# arc           |x|=1.264e+03 |dx|=1.138e-02 |r|=1.892e-10 (λ)

# SNES iteration  2, KSP iteration   0       |r|=5.734e-07

# SNES iteration  2, KSP iteration   1       |r|=2.037e-18

# SNES iteration  3 success = CONVERGED_FNORM_RELATIVE

# sub  0 [  6k] |x|=8.790e+02 |dx|=5.828e-09 |r|=2.826e-09 (u)

# sub  1 [  6k] |x|=5.537e+00 |dx|=9.089e-11 |r|=9.027e-09 (r)

# sub  2 [  2k] |x|=1.053e-04 |dx|=6.880e-15 |r|=2.800e-13 (c)

# all           |x|=8.790e+02 |dx|=5.829e-09 |r|=9.459e-09

# arc           |x|=1.264e+03 |dx|=7.197e-08 |r|=2.910e-10 (λ)


*** Continuation step 58

# SNES iteration  0

# sub  0 [  6k] |x|=8.945e+02 |dx|=5.828e-09 |r|=4.646e+01 (u)

# sub  1 [  6k] |x|=5.649e+00 |dx|=9.089e-11 |r|=1.029e+02 (r)

# sub  2 [  2k] |x|=1.274e-04 |dx|=6.880e-15 |r|=5.922e-13 (c)

# all           |x|=8.945e+02 |dx|=5.829e-09 |r|=1.129e+02

# arc           |x|=1.551e+03 |dx|=0.000e+00 |r|=2.910e-10 (λ)

# SNES iteration  0, KSP iteration   0       |r|=1.129e+02

# SNES iteration  0, KSP iteration   1       |r|=2.436e-10

# SNES iteration  1

# sub  0 [  6k] |x|=8.930e+02 |dx|=1.790e+00 |r|=3.158e-01 (u)

# sub  1 [  6k] |x|=5.638e+00 |dx|=1.855e-02 |r|=7.559e-01 (r)

# sub  2 [  2k] |x|=1.269e-04 |dx|=9.998e-07 |r|=2.394e-10 (c)

# all           |x|=8.930e+02 |dx|=1.790e+00 |r|=8.193e-01

# SNES iteration  1, KSP iteration   0       |r|=1.000e+00

# SNES iteration  1, KSP iteration   1       |r|=9.239e-12

# arc           |x|=1.586e+03 |dx|=3.443e+01 |r|=4.438e-10 (λ)

# SNES iteration  1, KSP iteration   0       |r|=7.130e-02

# SNES iteration  1, KSP iteration   1       |r|=1.693e-13

# SNES iteration  2 success = CONVERGED_FNORM_RELATIVE

# sub  0 [  6k] |x|=8.946e+02 |dx|=1.066e-03 |r|=3.656e-07 (u)

# sub  1 [  6k] |x|=5.649e+00 |dx|=1.611e-05 |r|=4.284e-07 (r)

# sub  2 [  2k] |x|=1.295e-04 |dx|=1.009e-09 |r|=3.536e-13 (c)

# all           |x|=8.946e+02 |dx|=1.066e-03 |r|=5.632e-07

# arc           |x|=1.586e+03 |dx|=8.047e-03 |r|=4.366e-10 (λ)


*** Continuation step 59

# SNES iteration  0

# sub  0 [  6k] |x|=9.103e+02 |dx|=1.066e-03 |r|=4.631e+01 (u)

# sub  1 [  6k] |x|=5.763e+00 |dx|=1.611e-05 |r|=1.076e+02 (r)

# sub  2 [  2k] |x|=1.539e-04 |dx|=1.009e-09 |r|=7.107e-13 (c)

# all           |x|=9.103e+02 |dx|=1.066e-03 |r|=1.172e+02

# arc           |x|=1.907e+03 |dx|=0.000e+00 |r|=4.366e-10 (λ)

# SNES iteration  0, KSP iteration   0       |r|=1.172e+02

# SNES iteration  0, KSP iteration   1       |r|=5.234e-10

# SNES iteration  1

# sub  0 [  6k] |x|=9.089e+02 |dx|=1.677e+00 |r|=2.749e-01 (u)

# sub  1 [  6k] |x|=5.752e+00 |dx|=1.712e-02 |r|=6.882e-01 (r)

# sub  2 [  2k] |x|=1.534e-04 |dx|=9.534e-07 |r|=5.216e-10 (c)

# all           |x|=9.089e+02 |dx|=1.677e+00 |r|=7.411e-01

# SNES iteration  1, KSP iteration   0       |r|=1.000e+00

# SNES iteration  1, KSP iteration   1       |r|=3.787e-12

# arc           |x|=1.943e+03 |dx|=3.577e+01 |r|=1.455e-10 (λ)

# SNES iteration  1, KSP iteration   0       |r|=6.290e-02

# SNES iteration  1, KSP iteration   1       |r|=2.827e-13

# SNES iteration  2 success = CONVERGED_FNORM_RELATIVE

# sub  0 [  6k] |x|=9.104e+02 |dx|=8.695e-04 |r|=2.381e-07 (u)

# sub  1 [  6k] |x|=5.763e+00 |dx|=1.285e-05 |r|=2.877e-07 (r)

# sub  2 [  2k] |x|=1.561e-04 |dx|=9.455e-10 |r|=4.756e-13 (c)

# all           |x|=9.104e+02 |dx|=8.696e-04 |r|=3.734e-07

# arc           |x|=1.943e+03 |dx|=5.698e-03 |r|=4.147e-10 (λ)


*** Continuation step 60

# SNES iteration  0

# sub  0 [  6k] |x|=9.262e+02 |dx|=8.695e-04 |r|=4.618e+01 (u)

# sub  1 [  6k] |x|=5.880e+00 |dx|=1.285e-05 |r|=1.122e+02 (r)

# sub  2 [  2k] |x|=1.829e-04 |dx|=9.455e-10 |r|=8.187e-13 (c)

# all           |x|=9.263e+02 |dx|=8.696e-04 |r|=1.213e+02

# arc           |x|=2.301e+03 |dx|=0.000e+00 |r|=4.147e-10 (λ)

# SNES iteration  0, KSP iteration   0       |r|=1.213e+02

# SNES iteration  0, KSP iteration   1       |r|=4.371e-10

# SNES iteration  1

# sub  0 [  6k] |x|=9.249e+02 |dx|=1.578e+00 |r|=2.423e-01 (u)

# sub  1 [  6k] |x|=5.869e+00 |dx|=1.589e-02 |r|=6.301e-01 (r)

# sub  2 [  2k] |x|=1.824e-04 |dx|=9.158e-07 |r|=4.353e-10 (c)

# all           |x|=9.249e+02 |dx|=1.578e+00 |r|=6.751e-01

# SNES iteration  1, KSP iteration   0       |r|=1.000e+00

# SNES iteration  1, KSP iteration   1       |r|=2.328e-12

# arc           |x|=2.338e+03 |dx|=3.710e+01 |r|=3.420e-10 (λ)

# SNES iteration  1, KSP iteration   0       |r|=5.586e-02

# SNES iteration  1, KSP iteration   1       |r|=6.674e-14

# SNES iteration  2 success = CONVERGED_FNORM_RELATIVE

# sub  0 [  6k] |x|=9.263e+02 |dx|=7.182e-04 |r|=1.606e-07 (u)

# sub  1 [  6k] |x|=5.880e+00 |dx|=1.043e-05 |r|=1.971e-07 (r)

# sub  2 [  2k] |x|=1.851e-04 |dx|=8.932e-10 |r|=4.235e-13 (c)

# all           |x|=9.263e+02 |dx|=7.182e-04 |r|=2.542e-07

# arc           |x|=2.338e+03 |dx|=4.074e-03 |r|=4.729e-10 (λ)


*** Continuation step 61

# SNES iteration  0

# sub  0 [  6k] |x|=9.423e+02 |dx|=7.182e-04 |r|=4.604e+01 (u)

# sub  1 [  6k] |x|=5.998e+00 |dx|=1.043e-05 |r|=1.167e+02 (r)

# sub  2 [  2k] |x|=2.142e-04 |dx|=8.932e-10 |r|=8.749e-13 (c)

# all           |x|=9.423e+02 |dx|=7.182e-04 |r|=1.254e+02

# arc           |x|=2.732e+03 |dx|=0.000e+00 |r|=4.729e-10 (λ)

# SNES iteration  0, KSP iteration   0       |r|=1.254e+02

# SNES iteration  0, KSP iteration   1       |r|=2.360e-10

# SNES iteration  1

# sub  0 [  6k] |x|=9.410e+02 |dx|=1.492e+00 |r|=2.162e-01 (u)

# sub  1 [  6k] |x|=5.988e+00 |dx|=1.483e-02 |r|=5.802e-01 (r)

# sub  2 [  2k] |x|=2.137e-04 |dx|=8.874e-07 |r|=2.334e-10 (c)

# all           |x|=9.410e+02 |dx|=1.492e+00 |r|=6.192e-01

# SNES iteration  1, KSP iteration   0       |r|=1.000e+00

# SNES iteration  1, KSP iteration   1       |r|=3.495e-12

# arc           |x|=2.770e+03 |dx|=3.843e+01 |r|=1.601e-10 (λ)

# SNES iteration  1, KSP iteration   0       |r|=4.995e-02

# SNES iteration  1, KSP iteration   1       |r|=6.528e-14

# SNES iteration  2 success = CONVERGED_FNORM_RELATIVE

# sub  0 [  6k] |x|=9.424e+02 |dx|=5.990e-04 |r|=1.115e-07 (u)

# sub  1 [  6k] |x|=5.998e+00 |dx|=8.611e-06 |r|=1.379e-07 (r)

# sub  2 [  2k] |x|=2.165e-04 |dx|=8.422e-10 |r|=4.815e-13 (c)

# all           |x|=9.424e+02 |dx|=5.990e-04 |r|=1.773e-07

# arc           |x|=2.770e+03 |dx|=2.961e-03 |r|=5.093e-11 (λ)


*** Continuation step 62

# SNES iteration  0

# sub  0 [  6k] |x|=9.585e+02 |dx|=5.990e-04 |r|=4.589e+01 (u)

# sub  1 [  6k] |x|=6.118e+00 |dx|=8.611e-06 |r|=1.211e+02 (r)

# sub  2 [  2k] |x|=2.480e-04 |dx|=8.422e-10 |r|=9.679e-13 (c)

# all           |x|=9.585e+02 |dx|=5.990e-04 |r|=1.295e+02

# arc           |x|=3.203e+03 |dx|=0.000e+00 |r|=5.093e-11 (λ)

# SNES iteration  0, KSP iteration   0       |r|=1.295e+02

# SNES iteration  0, KSP iteration   1       |r|=9.216e-11

# SNES iteration  1

# sub  0 [  6k] |x|=9.572e+02 |dx|=1.415e+00 |r|=1.950e-01 (u)

# sub  1 [  6k] |x|=6.108e+00 |dx|=1.391e-02 |r|=5.370e-01 (r)

# sub  2 [  2k] |x|=2.474e-04 |dx|=8.681e-07 |r|=8.528e-11 (c)

# all           |x|=9.572e+02 |dx|=1.415e+00 |r|=5.713e-01

# SNES iteration  1, KSP iteration   0       |r|=1.000e+00

# SNES iteration  1, KSP iteration   1       |r|=4.163e-12

# arc           |x|=3.243e+03 |dx|=3.977e+01 |r|=2.183e-11 (λ)

# SNES iteration  1, KSP iteration   0       |r|=4.494e-02

# SNES iteration  1, KSP iteration   1       |r|=4.924e-14

# SNES iteration  2 success = CONVERGED_FNORM_RELATIVE

# sub  0 [  6k] |x|=9.585e+02 |dx|=5.036e-04 |r|=7.919e-08 (u)

# sub  1 [  6k] |x|=6.118e+00 |dx|=7.208e-06 |r|=9.916e-08 (r)

# sub  2 [  2k] |x|=2.502e-04 |dx|=7.888e-10 |r|=5.437e-13 (c)

# all           |x|=9.585e+02 |dx|=5.037e-04 |r|=1.269e-07

# arc           |x|=3.243e+03 |dx|=2.197e-03 |r|=3.856e-10 (λ)


*** Continuation step 63

# SNES iteration  0

# sub  0 [  6k] |x|=9.747e+02 |dx|=5.036e-04 |r|=4.572e+01 (u)

# sub  1 [  6k] |x|=6.239e+00 |dx|=7.208e-06 |r|=1.254e+02 (r)

# sub  2 [  2k] |x|=2.841e-04 |dx|=7.888e-10 |r|=1.009e-12 (c)

# all           |x|=9.748e+02 |dx|=5.037e-04 |r|=1.335e+02

# arc           |x|=3.716e+03 |dx|=0.000e+00 |r|=3.856e-10 (λ)

# SNES iteration  0, KSP iteration   0       |r|=1.335e+02

# SNES iteration  0, KSP iteration   1       |r|=5.076e-10

# SNES iteration  1

# sub  0 [  6k] |x|=9.735e+02 |dx|=1.347e+00 |r|=1.775e-01 (u)

# sub  1 [  6k] |x|=6.230e+00 |dx|=1.310e-02 |r|=4.996e-01 (r)

# sub  2 [  2k] |x|=2.835e-04 |dx|=8.582e-07 |r|=5.059e-10 (c)

# all           |x|=9.735e+02 |dx|=1.348e+00 |r|=5.302e-01

# SNES iteration  1, KSP iteration   0       |r|=1.000e+00

# SNES iteration  1, KSP iteration   1       |r|=3.055e-12

# arc           |x|=3.757e+03 |dx|=4.110e+01 |r|=1.455e-11 (λ)

# SNES iteration  1, KSP iteration   0       |r|=4.066e-02

# SNES iteration  1, KSP iteration   1       |r|=4.984e-14

# SNES iteration  2 success = CONVERGED_FNORM_RELATIVE

# sub  0 [  6k] |x|=9.748e+02 |dx|=4.265e-04 |r|=5.738e-08 (u)

# sub  1 [  6k] |x|=6.239e+00 |dx|=6.105e-06 |r|=7.188e-08 (r)

# sub  2 [  2k] |x|=2.864e-04 |dx|=7.328e-10 |r|=5.996e-13 (c)

# all           |x|=9.748e+02 |dx|=4.265e-04 |r|=9.198e-08

# arc           |x|=3.757e+03 |dx|=1.667e-03 |r|=6.330e-10 (λ)

Deformed configuration and load-response curve¶

After the run the roof has snapped through and, for this geometry, snapped back. The bending stress resultant peaks along the curvature ridge at the crown. The load-displacement curve shows the characteristic limit-point softening, the load reversal during snap-back, and the subsequent post-buckling rise.

Final deformed cylindrical roof coloured by bending stress resultant magnitude. — Figure 3:Final deformed configuration coloured by bending stress resultant magnitude.

if comm.size == 1:
    sol_tdc = np.column_stack([λ_step, u_step])
    np.savetxt(f"{name}_tdc_solution.csv", sol_tdc, delimiter=",", header="λ,u", comments="")
    plot_load_displacement(u_step, λ_step, f"{name}.png")

Load factor versus crown displacement. — Figure 4:Load factor $\lambda$ versus downward crown displacement $-u_z$ traced by Crisfield arc-length continuation, showing the snap-through limit point, load reversal during snap-back, and post-buckling rise.

References¶

Sze, K. Y., Liu, X. H., & Lo, S. H. (2004). Popular benchmark problems for geometric nonlinear analysis of shells. Finite Elements in Analysis and Design, 40(11), 1551–1569. 10.1016/j.finel.2003.11.001
Hansbo, P., Larson, M. G., & Larsson, F. (2015). Tangential differential calculus and the finite element modeling of a large deformation elastic membrane problem. Computational Mechanics, 56(1), 87–95. 10.1007/s00466-015-1158-x
Schöllhammer, D., & Fries, T. P. (2019). Kirchhoff–Love shell theory based on tangential differential calculus. Computational Mechanics, 64, 113–131. 10.1007/s00466-018-1659-5
Neunteufel, M. (2021). Mixed finite element methods for nonlinear continuum mechanics and shells [Phdthesis, Technische Universität Wien]. 10.34726/hss.2021.85500
Crisfield, M. A. (1981). A fast incremental/iterative solution procedure that handles “snap-through.” Computers & Structures, 13(1–3), 55–62. 10.1016/0045-7949(81)90108-5
Arnold, D. N., & Brezzi, F. (1997). Locking-free finite element methods for shells. Mathematics of Computation, 66(217), 1–14. 10.1090/S0025-5718-97-00785-0