Reduced Wannier functions with Wannier#

In this tutorial we give a hands-on introduction to the reduced Wannier construction, which resolves the topological obstruction to exponentially localized Wannier functions by working in an enlarged (supercell) Hilbert subspace. Following the approach of Cole & Vanderbilt, Reduced Wannier representation for topological bands, Phys. Rev. B 111, 205139 (2025), we illustrate the method on the Haldane model and show how disentanglement plus maximal localization yields smooth, localized orbitals even when the primitive-cell bands carry nontrivial topology. This tutorial will make use of the Wannier class in pythtb, which implements the reduced Wannier construction as well as the maximal localization and disentanglement procedures.

See also

T. Cole and D. Vanderbilt, “Reduced Wannier representation for topological bands”, PRB 111, 205139 (2025). arXiv:2412.17084

from pythtb import Mesh, Wannier, WFArray
from pythtb.models import haldane
import numpy as np

Haldane model supercell construction#

We begin by constructing the Haldane model in the topological phase, in which case the non-zero Chern number of the occupied bands enforces a topological obstruction to constructing exponentially localized Wannier functions that respect the lattice symmetries.

# tight-binding parameters
delta = 1
t1 = 1
t2 = -0.4
prim_model = haldane(delta, t1, t2)

print(f"Chern number: {prim_model.chern_number((0, 1), (20, 20)):0.3f}")
Chern number: 1.000

To circumvent the obstruction, we use the procedure of constructing “reduced Wannier” functions that are localized in a smaller subspace of the original Hilbert space. Since the occupied space is one-dimensional, there exists no subspace. This is why we must construct a supercell, folding the bands back into the first Brillouin zone, to obtain more occupied bands to choose from. Here we will use a 2x2 supercell, which will give us 4 occupied bands to work with.

n_super_cell = 2
model = prim_model.make_supercell([[n_super_cell, 0], [0, n_super_cell]])
model.info(show=True, short=False)

Hide code cell output

----------------------------------------
       Tight-binding model report       
----------------------------------------
r-space dimension           = 2
k-space dimension           = 2
periodic directions         = [0, 1]
spinful                     = False
number of spin components   = 1
number of electronic states = 8
number of orbitals          = 8

Lattice vectors (Cartesian):
  # 0 ===> [ 2.000,  0.000]
  # 1 ===> [ 1.000,  1.732]
Volume of unit cell (Cartesian) = 3.464 [A^d]

Reciprocal lattice vectors (Cartesian):
  # 0 ===> [ 3.142, -1.814]
  # 1 ===> [ 0.000,  3.628]
Volume of reciprocal unit cell = 11.396 [A^-d]

Orbital vectors (Cartesian):
  # 0 ===> [ 0.500,  0.289]
  # 1 ===> [ 1.000,  0.577]
  # 2 ===> [ 1.000,  1.155]
  # 3 ===> [ 1.500,  1.443]
  # 4 ===> [ 1.500,  0.289]
  # 5 ===> [ 2.000,  0.577]
  # 6 ===> [ 2.000,  1.155]
  # 7 ===> [ 2.500,  1.443]

Orbital vectors (fractional):
  # 0 ===> [ 0.167,  0.167]
  # 1 ===> [ 0.333,  0.333]
  # 2 ===> [ 0.167,  0.667]
  # 3 ===> [ 0.333,  0.833]
  # 4 ===> [ 0.667,  0.167]
  # 5 ===> [ 0.833,  0.333]
  # 6 ===> [ 0.667,  0.667]
  # 7 ===> [ 0.833,  0.833]
----------------------------------------
Site energies:
  < 0 | H | 0 > = -1.000 
  < 1 | H | 1 > =  1.000 
  < 2 | H | 2 > = -1.000 
  < 3 | H | 3 > =  1.000 
  < 4 | H | 4 > = -1.000 
  < 5 | H | 5 > =  1.000 
  < 6 | H | 6 > = -1.000 
  < 7 | H | 7 > =  1.000 
Hoppings:
  < 0 | H | 1  + [ 0.0 ,  0.0 ] > = 1.0000+0.0000j
  < 0 | H | 5  + [-1.0 ,  0.0 ] > = 1.0000+0.0000j
  < 0 | H | 3  + [ 0.0 , -1.0 ] > = 1.0000+0.0000j
  < 0 | H | 4  + [ 0.0 ,  0.0 ] > = -0.0000-0.4000j
  < 0 | H | 6  + [-1.0 ,  0.0 ] > = -0.0000-0.4000j
  < 0 | H | 2  + [ 0.0 , -1.0 ] > = -0.0000-0.4000j
  < 1 | H | 5  + [-1.0 ,  0.0 ] > = -0.0000-0.4000j
  < 1 | H | 7  + [ 0.0 , -1.0 ] > = -0.0000-0.4000j
  < 1 | H | 3  + [ 0.0 ,  0.0 ] > = -0.0000-0.4000j
  < 2 | H | 3  + [ 0.0 ,  0.0 ] > = 1.0000+0.0000j
  < 2 | H | 7  + [-1.0 ,  0.0 ] > = 1.0000+0.0000j
  < 2 | H | 1  + [ 0.0 ,  0.0 ] > = 1.0000+0.0000j
  < 2 | H | 6  + [ 0.0 ,  0.0 ] > = -0.0000-0.4000j
  < 2 | H | 4  + [-1.0 ,  1.0 ] > = -0.0000-0.4000j
  < 2 | H | 0  + [ 0.0 ,  0.0 ] > = -0.0000-0.4000j
  < 3 | H | 7  + [-1.0 ,  0.0 ] > = -0.0000-0.4000j
  < 3 | H | 5  + [ 0.0 ,  0.0 ] > = -0.0000-0.4000j
  < 3 | H | 1  + [ 0.0 ,  1.0 ] > = -0.0000-0.4000j
  < 4 | H | 5  + [ 0.0 ,  0.0 ] > = 1.0000+0.0000j
  < 4 | H | 1  + [ 0.0 ,  0.0 ] > = 1.0000+0.0000j
  < 4 | H | 7  + [ 0.0 , -1.0 ] > = 1.0000+0.0000j
  < 4 | H | 0  + [ 1.0 ,  0.0 ] > = -0.0000-0.4000j
  < 4 | H | 2  + [ 0.0 ,  0.0 ] > = -0.0000-0.4000j
  < 4 | H | 6  + [ 0.0 , -1.0 ] > = -0.0000-0.4000j
  < 5 | H | 1  + [ 0.0 ,  0.0 ] > = -0.0000-0.4000j
  < 5 | H | 3  + [ 1.0 , -1.0 ] > = -0.0000-0.4000j
  < 5 | H | 7  + [ 0.0 ,  0.0 ] > = -0.0000-0.4000j
  < 6 | H | 7  + [ 0.0 ,  0.0 ] > = 1.0000+0.0000j
  < 6 | H | 3  + [ 0.0 ,  0.0 ] > = 1.0000+0.0000j
  < 6 | H | 5  + [ 0.0 ,  0.0 ] > = 1.0000+0.0000j
  < 6 | H | 2  + [ 1.0 ,  0.0 ] > = -0.0000-0.4000j
  < 6 | H | 0  + [ 0.0 ,  1.0 ] > = -0.0000-0.4000j
  < 6 | H | 4  + [ 0.0 ,  0.0 ] > = -0.0000-0.4000j
  < 7 | H | 3  + [ 0.0 ,  0.0 ] > = -0.0000-0.4000j
  < 7 | H | 1  + [ 1.0 ,  0.0 ] > = -0.0000-0.4000j
  < 7 | H | 5  + [ 0.0 ,  1.0 ] > = -0.0000-0.4000j
Hopping distances:
  | pos( 0 ) - pos( 1 ) + [ 0.0 ,  0.0 ] | =   0.577
  | pos( 0 ) - pos( 5 ) + [-1.0 ,  0.0 ] | =   0.577
  | pos( 0 ) - pos( 3 ) + [ 0.0 , -1.0 ] | =   0.577
  | pos( 0 ) - pos( 4 ) + [ 0.0 ,  0.0 ] | =   1.000
  | pos( 0 ) - pos( 6 ) + [-1.0 ,  0.0 ] | =   1.000
  | pos( 0 ) - pos( 2 ) + [ 0.0 , -1.0 ] | =   1.000
  | pos( 1 ) - pos( 5 ) + [-1.0 ,  0.0 ] | =   1.000
  | pos( 1 ) - pos( 7 ) + [ 0.0 , -1.0 ] | =   1.000
  | pos( 1 ) - pos( 3 ) + [ 0.0 ,  0.0 ] | =   1.000
  | pos( 2 ) - pos( 3 ) + [ 0.0 ,  0.0 ] | =   0.577
  | pos( 2 ) - pos( 7 ) + [-1.0 ,  0.0 ] | =   0.577
  | pos( 2 ) - pos( 1 ) + [ 0.0 ,  0.0 ] | =   0.577
  | pos( 2 ) - pos( 6 ) + [ 0.0 ,  0.0 ] | =   1.000
  | pos( 2 ) - pos( 4 ) + [-1.0 ,  1.0 ] | =   1.000
  | pos( 2 ) - pos( 0 ) + [ 0.0 ,  0.0 ] | =   1.000
  | pos( 3 ) - pos( 7 ) + [-1.0 ,  0.0 ] | =   1.000
  | pos( 3 ) - pos( 5 ) + [ 0.0 ,  0.0 ] | =   1.000
  | pos( 3 ) - pos( 1 ) + [ 0.0 ,  1.0 ] | =   1.000
  | pos( 4 ) - pos( 5 ) + [ 0.0 ,  0.0 ] | =   0.577
  | pos( 4 ) - pos( 1 ) + [ 0.0 ,  0.0 ] | =   0.577
  | pos( 4 ) - pos( 7 ) + [ 0.0 , -1.0 ] | =   0.577
  | pos( 4 ) - pos( 0 ) + [ 1.0 ,  0.0 ] | =   1.000
  | pos( 4 ) - pos( 2 ) + [ 0.0 ,  0.0 ] | =   1.000
  | pos( 4 ) - pos( 6 ) + [ 0.0 , -1.0 ] | =   1.000
  | pos( 5 ) - pos( 1 ) + [ 0.0 ,  0.0 ] | =   1.000
  | pos( 5 ) - pos( 3 ) + [ 1.0 , -1.0 ] | =   1.000
  | pos( 5 ) - pos( 7 ) + [ 0.0 ,  0.0 ] | =   1.000
  | pos( 6 ) - pos( 7 ) + [ 0.0 ,  0.0 ] | =   0.577
  | pos( 6 ) - pos( 3 ) + [ 0.0 ,  0.0 ] | =   0.577
  | pos( 6 ) - pos( 5 ) + [ 0.0 ,  0.0 ] | =   0.577
  | pos( 6 ) - pos( 2 ) + [ 1.0 ,  0.0 ] | =   1.000
  | pos( 6 ) - pos( 0 ) + [ 0.0 ,  1.0 ] | =   1.000
  | pos( 6 ) - pos( 4 ) + [ 0.0 ,  0.0 ] | =   1.000
  | pos( 7 ) - pos( 3 ) + [ 0.0 ,  0.0 ] | =   1.000
  | pos( 7 ) - pos( 1 ) + [ 1.0 ,  0.0 ] | =   1.000
  | pos( 7 ) - pos( 5 ) + [ 0.0 ,  1.0 ] | =   1.000

We construct the WFArray and diagonalize the model on a semi-full k-mesh. It is important that the mesh not include the endpoints \(k_i=1\), which correspond to the boundaries of the Brillouin zone. The Fourier transform requires a well-defined periodicity, which is disrupted by including these points. Therefore, we will use a k-mesh that spans the interior of the Brillouin zone, avoiding the boundaries. This is the default behavior of Mesh.build_grid.

.. note:: The Mesh.build_grid function automatically imposes periodic boundary conditions on the k-mesh. We can also explicitly state that the k-mesh is periodic, with the topology of a torus, we can use the Mesh.wind_bz function, specifying the mesh axis and k-component that is wrapped. Here, it is not necessary since the Mesh.build_grid function already does this for us. In other cases where we use a custom k-mesh, we may need to use Mesh.wind_bz to impose periodic boundary conditions.

nks = 20, 20  # number of k points along each dimension
mesh = Mesh(dim_k=2, axis_types=["k", "k"])
mesh.build_grid(shape=nks)
print(mesh)
Mesh Summary
========================================
Type: grid
Dimensionality: 2 k-dim(s) + 0 λ-dim(s)
Number of mesh points: 400
Full shape: (20, 20, 2)
k-axes: [Axis(type=k, name=k_0, size=20), Axis(type=k, name=k_1, size=20)]
λ-axes: []
Is a torus in k-space (all k-axes wind BZ): yes
Loops: (axis 0, comp 0, winds_bz=yes, closed=no), (axis 1, comp 1, winds_bz=yes, closed=no)

Now we pass this mesh to the WFArray constructor and solve the mesh.

wfa = WFArray(model.lattice, mesh)
wfa.solve_model(model)

We know that Wannierizing the full set of 4 occupied bands is obstructed by the topology of the band structure. We can try the next best thing and Wannierize a 3-dimensional subspace. To do this, we will choose a set of three trial wavefunctions centered on 3 of the low energy orbitals, where we would expect the localized Wannier functions of the trivial occupied bands to be located

n_orb = model.norb  # number of orbitals
n_occ = int(n_orb / 2)  # number of occupied bands (assume half-filling)

low_E_sites = np.arange(
    0, n_orb, 2
)  # low-energy sites defined to be indexed by even numbers
high_E_sites = np.arange(
    1, n_orb, 2
)  # high-energy sites defined to be indexed by odd numbers

omit_site = 6  # omitting one of the low energy sites
sites = list(np.setdiff1d(low_E_sites, [omit_site]))
tf_list = [
    [(orb, 1)] for orb in sites
]  # trial wavefunctions in form of [(orbital index, weight)]

n_tfs = len(tf_list)

print(f"Trial wavefunctions: {tf_list}")
print(f"# of Wannier functions: {n_tfs}")
print(f"# of occupied bands: {n_occ}")
print(f"Wannier fraction: {n_tfs / n_occ}")
Trial wavefunctions: [[(np.int64(0), 1)], [(np.int64(2), 1)], [(np.int64(4), 1)]]
# of Wannier functions: 3
# of occupied bands: 4
Wannier fraction: 0.75

Projection onto trial states#

Next, we initialize the Wannier object with the TBModel and WFArray objects. We initialize the Bloch-like states with project function which aligns the trial wavefunctions with the target bands specified by band_idxs.

WF = Wannier(wfa)

WF.project(tf_list, band_idxs=list(range(n_occ)))

This already gives us a set of Wannier functions that are exponentially localized, showing that this is a trivial subsapce of the obstructed manifold.

WF.info()
Wannier Function Report
===============================================================
WF 1: center = [0.50229186, 0.28999834]  Omega     = 0.27158248
WF 2: center = [1.00000000, 1.15205413]  Omega     = 0.27158248
WF 3: center = [1.49770814, 0.28999834]  Omega     = 0.27158248
Sum : center = [3.00000000, 1.73205081]  Omega tot = 0.81474745
---------------------------------------------------------------
Omega I   = 0.73887684
Omega D   = 0.01751802
Omega OD  = 0.05835260
Omega tot = 0.81474745

Disentanglement#

We can make these states even more localized with subspace selectio via the disentanglement procedure. This picks the subspace of the 4-band manfiold that minimizes the gauge-independent spread.

frozen_window = None  # frozen window in energy
outer_window = [-4, 0]  # outer window in energy

WF.disentangle(
    n_wfs=3,
    frozen_window=frozen_window,
    outer_window=outer_window,
    verbose=True,
    tf_speedup=True,
    max_iter=500,
    tol=1e-10,
)

Hide code cell output

Initial Omega_I: 0.7388768356161678
WARNING: All log messages before absl::InitializeLog() is called are written to STDERR
I0000 00:00:1779816470.412271    1211 port.cc:153] oneDNN custom operations are on. You may see slightly different numerical results due to floating-point round-off errors from different computation orders. To turn them off, set the environment variable `TF_ENABLE_ONEDNN_OPTS=0`.
I0000 00:00:1779816470.422218    1211 cudart_stub.cc:31] Could not find cuda drivers on your machine, GPU will not be used.
I0000 00:00:1779816471.315852    1211 cpu_feature_guard.cc:227] This TensorFlow binary is optimized to use available CPU instructions in performance-critical operations.
To enable the following instructions: AVX2 AVX512F AVX512_VNNI AVX512_BF16 FMA, in other operations, rebuild TensorFlow with the appropriate compiler flags.
WARNING: All log messages before absl::InitializeLog() is called are written to STDERR
I0000 00:00:1779816474.667384    1211 port.cc:153] oneDNN custom operations are on. You may see slightly different numerical results due to floating-point round-off errors from different computation orders. To turn them off, set the environment variable `TF_ENABLE_ONEDNN_OPTS=0`.
I0000 00:00:1779816474.672382    1211 cudart_stub.cc:31] Could not find cuda drivers on your machine, GPU will not be used.
E0000 00:00:1779816475.567401    1211 cuda_platform.cc:52] failed call to cuInit: INTERNAL: CUDA error: Failed call to cuInit: UNKNOWN ERROR (303)
iter    0 | Ω_I = 7.351166932e-01 | ΔΩ = -3.76014e-03
iter    1 | Ω_I = 7.318588717e-01 | ΔΩ = -3.25782e-03
iter    2 | Ω_I = 7.289963291e-01 | ΔΩ = -2.86254e-03
iter    3 | Ω_I = 7.264317378e-01 | ΔΩ = -2.56459e-03
iter    4 | Ω_I = 7.241199320e-01 | ΔΩ = -2.31181e-03
iter    5 | Ω_I = 7.220133072e-01 | ΔΩ = -2.10662e-03
iter    6 | Ω_I = 7.200724681e-01 | ΔΩ = -1.94084e-03
iter    7 | Ω_I = 7.182764535e-01 | ΔΩ = -1.79601e-03
iter    8 | Ω_I = 7.165979583e-01 | ΔΩ = -1.67850e-03
iter    9 | Ω_I = 7.150287030e-01 | ΔΩ = -1.56926e-03
iter   10 | Ω_I = 7.135503744e-01 | ΔΩ = -1.47833e-03
iter   11 | Ω_I = 7.121519414e-01 | ΔΩ = -1.39843e-03
iter   12 | Ω_I = 7.108224077e-01 | ΔΩ = -1.32953e-03
iter   13 | Ω_I = 7.095568950e-01 | ΔΩ = -1.26551e-03
iter   14 | Ω_I = 7.083470257e-01 | ΔΩ = -1.20987e-03
iter   15 | Ω_I = 7.071838033e-01 | ΔΩ = -1.16322e-03
iter   16 | Ω_I = 7.060706009e-01 | ΔΩ = -1.11320e-03
iter   17 | Ω_I = 7.049976651e-01 | ΔΩ = -1.07294e-03
iter   18 | Ω_I = 7.039620376e-01 | ΔΩ = -1.03563e-03
iter   19 | Ω_I = 7.029607689e-01 | ΔΩ = -1.00127e-03
iter   20 | Ω_I = 7.019978599e-01 | ΔΩ = -9.62909e-04
iter   21 | Ω_I = 7.010595038e-01 | ΔΩ = -9.38356e-04
iter   22 | Ω_I = 7.001495143e-01 | ΔΩ = -9.09990e-04
iter   23 | Ω_I = 6.992671529e-01 | ΔΩ = -8.82361e-04
iter   24 | Ω_I = 6.984133852e-01 | ΔΩ = -8.53768e-04
iter   25 | Ω_I = 6.975808188e-01 | ΔΩ = -8.32566e-04
iter   26 | Ω_I = 6.967677537e-01 | ΔΩ = -8.13065e-04
iter   27 | Ω_I = 6.959766306e-01 | ΔΩ = -7.91123e-04
iter   28 | Ω_I = 6.952110022e-01 | ΔΩ = -7.65628e-04
iter   29 | Ω_I = 6.944599578e-01 | ΔΩ = -7.51044e-04
iter   30 | Ω_I = 6.937276796e-01 | ΔΩ = -7.32278e-04
iter   31 | Ω_I = 6.930156823e-01 | ΔΩ = -7.11997e-04
iter   32 | Ω_I = 6.923188587e-01 | ΔΩ = -6.96824e-04
iter   33 | Ω_I = 6.916390343e-01 | ΔΩ = -6.79824e-04
iter   34 | Ω_I = 6.909733231e-01 | ΔΩ = -6.65711e-04
iter   35 | Ω_I = 6.903267742e-01 | ΔΩ = -6.46549e-04
iter   36 | Ω_I = 6.896922877e-01 | ΔΩ = -6.34486e-04
iter   37 | Ω_I = 6.890676022e-01 | ΔΩ = -6.24685e-04
iter   38 | Ω_I = 6.884623761e-01 | ΔΩ = -6.05226e-04
iter   39 | Ω_I = 6.878660844e-01 | ΔΩ = -5.96292e-04
iter   40 | Ω_I = 6.872874019e-01 | ΔΩ = -5.78683e-04
iter   41 | Ω_I = 6.867191895e-01 | ΔΩ = -5.68212e-04
iter   42 | Ω_I = 6.861625258e-01 | ΔΩ = -5.56664e-04
iter   43 | Ω_I = 6.856169607e-01 | ΔΩ = -5.45565e-04
iter   44 | Ω_I = 6.850849164e-01 | ΔΩ = -5.32044e-04
iter   45 | Ω_I = 6.845614562e-01 | ΔΩ = -5.23460e-04
iter   46 | Ω_I = 6.840538874e-01 | ΔΩ = -5.07569e-04
iter   47 | Ω_I = 6.835550108e-01 | ΔΩ = -4.98877e-04
iter   48 | Ω_I = 6.830635226e-01 | ΔΩ = -4.91488e-04
iter   49 | Ω_I = 6.825841277e-01 | ΔΩ = -4.79395e-04
iter   50 | Ω_I = 6.821178277e-01 | ΔΩ = -4.66300e-04
iter   51 | Ω_I = 6.816540239e-01 | ΔΩ = -4.63804e-04
iter   52 | Ω_I = 6.812038603e-01 | ΔΩ = -4.50164e-04
iter   53 | Ω_I = 6.807622755e-01 | ΔΩ = -4.41585e-04
iter   54 | Ω_I = 6.803287698e-01 | ΔΩ = -4.33506e-04
iter   55 | Ω_I = 6.799056789e-01 | ΔΩ = -4.23091e-04
iter   56 | Ω_I = 6.794873280e-01 | ΔΩ = -4.18351e-04
iter   57 | Ω_I = 6.790833581e-01 | ΔΩ = -4.03970e-04
iter   58 | Ω_I = 6.786817633e-01 | ΔΩ = -4.01595e-04
iter   59 | Ω_I = 6.782910015e-01 | ΔΩ = -3.90762e-04
iter   60 | Ω_I = 6.779054260e-01 | ΔΩ = -3.85576e-04
iter   61 | Ω_I = 6.775317437e-01 | ΔΩ = -3.73682e-04
iter   62 | Ω_I = 6.771633685e-01 | ΔΩ = -3.68375e-04
iter   63 | Ω_I = 6.768023003e-01 | ΔΩ = -3.61068e-04
iter   64 | Ω_I = 6.764512152e-01 | ΔΩ = -3.51085e-04
iter   65 | Ω_I = 6.761030724e-01 | ΔΩ = -3.48143e-04
iter   66 | Ω_I = 6.757666241e-01 | ΔΩ = -3.36448e-04
iter   67 | Ω_I = 6.754301743e-01 | ΔΩ = -3.36450e-04
iter   68 | Ω_I = 6.751017373e-01 | ΔΩ = -3.28437e-04
iter   69 | Ω_I = 6.747857396e-01 | ΔΩ = -3.15998e-04
iter   70 | Ω_I = 6.744699502e-01 | ΔΩ = -3.15789e-04
iter   71 | Ω_I = 6.741586216e-01 | ΔΩ = -3.11329e-04
iter   72 | Ω_I = 6.738575695e-01 | ΔΩ = -3.01052e-04
iter   73 | Ω_I = 6.735666825e-01 | ΔΩ = -2.90887e-04
iter   74 | Ω_I = 6.732713971e-01 | ΔΩ = -2.95285e-04
iter   75 | Ω_I = 6.729935405e-01 | ΔΩ = -2.77857e-04
iter   76 | Ω_I = 6.727116638e-01 | ΔΩ = -2.81877e-04
iter   77 | Ω_I = 6.724403879e-01 | ΔΩ = -2.71276e-04
iter   78 | Ω_I = 6.721693991e-01 | ΔΩ = -2.70989e-04
iter   79 | Ω_I = 6.719060678e-01 | ΔΩ = -2.63331e-04
iter   80 | Ω_I = 6.716521950e-01 | ΔΩ = -2.53873e-04
iter   81 | Ω_I = 6.713961083e-01 | ΔΩ = -2.56087e-04
iter   82 | Ω_I = 6.711513946e-01 | ΔΩ = -2.44714e-04
iter   83 | Ω_I = 6.709091140e-01 | ΔΩ = -2.42281e-04
iter   84 | Ω_I = 6.706725903e-01 | ΔΩ = -2.36524e-04
iter   85 | Ω_I = 6.704365720e-01 | ΔΩ = -2.36018e-04
iter   86 | Ω_I = 6.702060635e-01 | ΔΩ = -2.30509e-04
iter   87 | Ω_I = 6.699812283e-01 | ΔΩ = -2.24835e-04
iter   88 | Ω_I = 6.697635582e-01 | ΔΩ = -2.17670e-04
iter   89 | Ω_I = 6.695454575e-01 | ΔΩ = -2.18101e-04
iter   90 | Ω_I = 6.693326193e-01 | ΔΩ = -2.12838e-04
iter   91 | Ω_I = 6.691260166e-01 | ΔΩ = -2.06603e-04
iter   92 | Ω_I = 6.689240539e-01 | ΔΩ = -2.01963e-04
iter   93 | Ω_I = 6.687251747e-01 | ΔΩ = -1.98879e-04
iter   94 | Ω_I = 6.685286006e-01 | ΔΩ = -1.96574e-04
iter   95 | Ω_I = 6.683363847e-01 | ΔΩ = -1.92216e-04
iter   96 | Ω_I = 6.681490604e-01 | ΔΩ = -1.87324e-04
iter   97 | Ω_I = 6.679633203e-01 | ΔΩ = -1.85740e-04
iter   98 | Ω_I = 6.677838989e-01 | ΔΩ = -1.79421e-04
iter   99 | Ω_I = 6.676038806e-01 | ΔΩ = -1.80018e-04
iter  100 | Ω_I = 6.674312821e-01 | ΔΩ = -1.72599e-04
iter  101 | Ω_I = 6.672606475e-01 | ΔΩ = -1.70635e-04
iter  102 | Ω_I = 6.670937440e-01 | ΔΩ = -1.66903e-04
iter  103 | Ω_I = 6.669291534e-01 | ΔΩ = -1.64591e-04
iter  104 | Ω_I = 6.667673873e-01 | ΔΩ = -1.61766e-04
iter  105 | Ω_I = 6.666116967e-01 | ΔΩ = -1.55691e-04
iter  106 | Ω_I = 6.664543069e-01 | ΔΩ = -1.57390e-04
iter  107 | Ω_I = 6.663038570e-01 | ΔΩ = -1.50450e-04
iter  108 | Ω_I = 6.661512669e-01 | ΔΩ = -1.52590e-04
iter  109 | Ω_I = 6.660052260e-01 | ΔΩ = -1.46041e-04
iter  110 | Ω_I = 6.658650360e-01 | ΔΩ = -1.40190e-04
iter  111 | Ω_I = 6.657234845e-01 | ΔΩ = -1.41551e-04
iter  112 | Ω_I = 6.655846905e-01 | ΔΩ = -1.38794e-04
iter  113 | Ω_I = 6.654479489e-01 | ΔΩ = -1.36742e-04
iter  114 | Ω_I = 6.653163482e-01 | ΔΩ = -1.31601e-04
iter  115 | Ω_I = 6.651878872e-01 | ΔΩ = -1.28461e-04
iter  116 | Ω_I = 6.650559901e-01 | ΔΩ = -1.31897e-04
iter  117 | Ω_I = 6.649332955e-01 | ΔΩ = -1.22695e-04
iter  118 | Ω_I = 6.648075986e-01 | ΔΩ = -1.25697e-04
iter  119 | Ω_I = 6.646888379e-01 | ΔΩ = -1.18761e-04
iter  120 | Ω_I = 6.645712809e-01 | ΔΩ = -1.17557e-04
iter  121 | Ω_I = 6.644530230e-01 | ΔΩ = -1.18258e-04
iter  122 | Ω_I = 6.643395785e-01 | ΔΩ = -1.13445e-04
iter  123 | Ω_I = 6.642288027e-01 | ΔΩ = -1.10776e-04
iter  124 | Ω_I = 6.641188750e-01 | ΔΩ = -1.09928e-04
iter  125 | Ω_I = 6.640102889e-01 | ΔΩ = -1.08586e-04
iter  126 | Ω_I = 6.639058370e-01 | ΔΩ = -1.04452e-04
iter  127 | Ω_I = 6.638020325e-01 | ΔΩ = -1.03804e-04
iter  128 | Ω_I = 6.636998854e-01 | ΔΩ = -1.02147e-04
iter  129 | Ω_I = 6.635972009e-01 | ΔΩ = -1.02684e-04
iter  130 | Ω_I = 6.635012301e-01 | ΔΩ = -9.59708e-05
iter  131 | Ω_I = 6.634078942e-01 | ΔΩ = -9.33358e-05
iter  132 | Ω_I = 6.633102087e-01 | ΔΩ = -9.76856e-05
iter  133 | Ω_I = 6.632240272e-01 | ΔΩ = -8.61814e-05
iter  134 | Ω_I = 6.631295977e-01 | ΔΩ = -9.44295e-05
iter  135 | Ω_I = 6.630399383e-01 | ΔΩ = -8.96594e-05
iter  136 | Ω_I = 6.629521041e-01 | ΔΩ = -8.78342e-05
iter  137 | Ω_I = 6.628672228e-01 | ΔΩ = -8.48812e-05
iter  138 | Ω_I = 6.627844708e-01 | ΔΩ = -8.27520e-05
iter  139 | Ω_I = 6.627010649e-01 | ΔΩ = -8.34059e-05
iter  140 | Ω_I = 6.626215742e-01 | ΔΩ = -7.94907e-05
iter  141 | Ω_I = 6.625414908e-01 | ΔΩ = -8.00834e-05
iter  142 | Ω_I = 6.624616127e-01 | ΔΩ = -7.98781e-05
iter  143 | Ω_I = 6.623837754e-01 | ΔΩ = -7.78373e-05
iter  144 | Ω_I = 6.623123884e-01 | ΔΩ = -7.13870e-05
iter  145 | Ω_I = 6.622387271e-01 | ΔΩ = -7.36614e-05
iter  146 | Ω_I = 6.621660523e-01 | ΔΩ = -7.26747e-05
iter  147 | Ω_I = 6.620972684e-01 | ΔΩ = -6.87839e-05
iter  148 | Ω_I = 6.620227642e-01 | ΔΩ = -7.45042e-05
iter  149 | Ω_I = 6.619559584e-01 | ΔΩ = -6.68058e-05
iter  150 | Ω_I = 6.618912168e-01 | ΔΩ = -6.47416e-05
iter  151 | Ω_I = 6.618286513e-01 | ΔΩ = -6.25655e-05
iter  152 | Ω_I = 6.617582919e-01 | ΔΩ = -7.03593e-05
iter  153 | Ω_I = 6.616969772e-01 | ΔΩ = -6.13148e-05
iter  154 | Ω_I = 6.616333569e-01 | ΔΩ = -6.36202e-05
iter  155 | Ω_I = 6.615720330e-01 | ΔΩ = -6.13239e-05
iter  156 | Ω_I = 6.615121949e-01 | ΔΩ = -5.98381e-05
iter  157 | Ω_I = 6.614518339e-01 | ΔΩ = -6.03610e-05
iter  158 | Ω_I = 6.613947223e-01 | ΔΩ = -5.71115e-05
iter  159 | Ω_I = 6.613382156e-01 | ΔΩ = -5.65068e-05
iter  160 | Ω_I = 6.612833955e-01 | ΔΩ = -5.48201e-05
iter  161 | Ω_I = 6.612267631e-01 | ΔΩ = -5.66324e-05
iter  162 | Ω_I = 6.611766714e-01 | ΔΩ = -5.00917e-05
iter  163 | Ω_I = 6.611259528e-01 | ΔΩ = -5.07186e-05
iter  164 | Ω_I = 6.610725285e-01 | ΔΩ = -5.34243e-05
iter  165 | Ω_I = 6.610191963e-01 | ΔΩ = -5.33322e-05
iter  166 | Ω_I = 6.609709403e-01 | ΔΩ = -4.82560e-05
iter  167 | Ω_I = 6.609211569e-01 | ΔΩ = -4.97834e-05
iter  168 | Ω_I = 6.608711653e-01 | ΔΩ = -4.99916e-05
iter  169 | Ω_I = 6.608242078e-01 | ΔΩ = -4.69575e-05
iter  170 | Ω_I = 6.607769186e-01 | ΔΩ = -4.72892e-05
iter  171 | Ω_I = 6.607332276e-01 | ΔΩ = -4.36910e-05
iter  172 | Ω_I = 6.606898754e-01 | ΔΩ = -4.33522e-05
iter  173 | Ω_I = 6.606460474e-01 | ΔΩ = -4.38280e-05
iter  174 | Ω_I = 6.606040343e-01 | ΔΩ = -4.20131e-05
iter  175 | Ω_I = 6.605578384e-01 | ΔΩ = -4.61959e-05
iter  176 | Ω_I = 6.605179771e-01 | ΔΩ = -3.98613e-05
iter  177 | Ω_I = 6.604756561e-01 | ΔΩ = -4.23210e-05
iter  178 | Ω_I = 6.604369932e-01 | ΔΩ = -3.86629e-05
iter  179 | Ω_I = 6.603955477e-01 | ΔΩ = -4.14454e-05
iter  180 | Ω_I = 6.603582218e-01 | ΔΩ = -3.73260e-05
iter  181 | Ω_I = 6.603214865e-01 | ΔΩ = -3.67352e-05
iter  182 | Ω_I = 6.602783061e-01 | ΔΩ = -4.31805e-05
iter  183 | Ω_I = 6.602477239e-01 | ΔΩ = -3.05821e-05
iter  184 | Ω_I = 6.602096814e-01 | ΔΩ = -3.80425e-05
iter  185 | Ω_I = 6.601749593e-01 | ΔΩ = -3.47221e-05
iter  186 | Ω_I = 6.601403103e-01 | ΔΩ = -3.46490e-05
iter  187 | Ω_I = 6.601029956e-01 | ΔΩ = -3.73148e-05
iter  188 | Ω_I = 6.600744796e-01 | ΔΩ = -2.85160e-05
iter  189 | Ω_I = 6.600375358e-01 | ΔΩ = -3.69438e-05
iter  190 | Ω_I = 6.600065956e-01 | ΔΩ = -3.09403e-05
iter  191 | Ω_I = 6.599737813e-01 | ΔΩ = -3.28143e-05
iter  192 | Ω_I = 6.599476474e-01 | ΔΩ = -2.61339e-05
iter  193 | Ω_I = 6.599143257e-01 | ΔΩ = -3.33217e-05
iter  194 | Ω_I = 6.598847707e-01 | ΔΩ = -2.95549e-05
iter  195 | Ω_I = 6.598532022e-01 | ΔΩ = -3.15685e-05
iter  196 | Ω_I = 6.598236272e-01 | ΔΩ = -2.95751e-05
iter  197 | Ω_I = 6.597925711e-01 | ΔΩ = -3.10561e-05
iter  198 | Ω_I = 6.597662257e-01 | ΔΩ = -2.63454e-05
iter  199 | Ω_I = 6.597334997e-01 | ΔΩ = -3.27260e-05
iter  200 | Ω_I = 6.597115023e-01 | ΔΩ = -2.19974e-05
iter  201 | Ω_I = 6.596825033e-01 | ΔΩ = -2.89989e-05
iter  202 | Ω_I = 6.596560487e-01 | ΔΩ = -2.64546e-05
iter  203 | Ω_I = 6.596303867e-01 | ΔΩ = -2.56620e-05
iter  204 | Ω_I = 6.596021708e-01 | ΔΩ = -2.82159e-05
iter  205 | Ω_I = 6.595777553e-01 | ΔΩ = -2.44155e-05
iter  206 | Ω_I = 6.595535321e-01 | ΔΩ = -2.42232e-05
iter  207 | Ω_I = 6.595306676e-01 | ΔΩ = -2.28645e-05
iter  208 | Ω_I = 6.595064450e-01 | ΔΩ = -2.42226e-05
iter  209 | Ω_I = 6.594819346e-01 | ΔΩ = -2.45103e-05
iter  210 | Ω_I = 6.594597364e-01 | ΔΩ = -2.21983e-05
iter  211 | Ω_I = 6.594403932e-01 | ΔΩ = -1.93432e-05
iter  212 | Ω_I = 6.594166906e-01 | ΔΩ = -2.37026e-05
iter  213 | Ω_I = 6.593945124e-01 | ΔΩ = -2.21783e-05
iter  214 | Ω_I = 6.593756275e-01 | ΔΩ = -1.88848e-05
iter  215 | Ω_I = 6.593487618e-01 | ΔΩ = -2.68657e-05
iter  216 | Ω_I = 6.593281104e-01 | ΔΩ = -2.06513e-05
iter  217 | Ω_I = 6.593132007e-01 | ΔΩ = -1.49097e-05
iter  218 | Ω_I = 6.592905125e-01 | ΔΩ = -2.26882e-05
iter  219 | Ω_I = 6.592688767e-01 | ΔΩ = -2.16358e-05
iter  220 | Ω_I = 6.592547484e-01 | ΔΩ = -1.41283e-05
iter  221 | Ω_I = 6.592348250e-01 | ΔΩ = -1.99234e-05
iter  222 | Ω_I = 6.592132126e-01 | ΔΩ = -2.16124e-05
iter  223 | Ω_I = 6.591994046e-01 | ΔΩ = -1.38080e-05
iter  224 | Ω_I = 6.591750896e-01 | ΔΩ = -2.43149e-05
iter  225 | Ω_I = 6.591585226e-01 | ΔΩ = -1.65671e-05
iter  226 | Ω_I = 6.591445692e-01 | ΔΩ = -1.39534e-05
iter  227 | Ω_I = 6.591224128e-01 | ΔΩ = -2.21563e-05
iter  228 | Ω_I = 6.591071215e-01 | ΔΩ = -1.52913e-05
iter  229 | Ω_I = 6.590899119e-01 | ΔΩ = -1.72096e-05
iter  230 | Ω_I = 6.590729469e-01 | ΔΩ = -1.69649e-05
iter  231 | Ω_I = 6.590593874e-01 | ΔΩ = -1.35596e-05
iter  232 | Ω_I = 6.590412383e-01 | ΔΩ = -1.81491e-05
iter  233 | Ω_I = 6.590267245e-01 | ΔΩ = -1.45137e-05
iter  234 | Ω_I = 6.590099275e-01 | ΔΩ = -1.67970e-05
iter  235 | Ω_I = 6.589902795e-01 | ΔΩ = -1.96480e-05
iter  236 | Ω_I = 6.589800245e-01 | ΔΩ = -1.02550e-05
iter  237 | Ω_I = 6.589661444e-01 | ΔΩ = -1.38802e-05
iter  238 | Ω_I = 6.589491368e-01 | ΔΩ = -1.70076e-05
iter  239 | Ω_I = 6.589333051e-01 | ΔΩ = -1.58316e-05
iter  240 | Ω_I = 6.589187514e-01 | ΔΩ = -1.45537e-05
iter  241 | Ω_I = 6.589095753e-01 | ΔΩ = -9.17615e-06
iter  242 | Ω_I = 6.588955858e-01 | ΔΩ = -1.39895e-05
iter  243 | Ω_I = 6.588789578e-01 | ΔΩ = -1.66280e-05
iter  244 | Ω_I = 6.588665684e-01 | ΔΩ = -1.23894e-05
iter  245 | Ω_I = 6.588518065e-01 | ΔΩ = -1.47619e-05
iter  246 | Ω_I = 6.588434560e-01 | ΔΩ = -8.35057e-06
iter  247 | Ω_I = 6.588307225e-01 | ΔΩ = -1.27335e-05
iter  248 | Ω_I = 6.588143828e-01 | ΔΩ = -1.63396e-05
iter  249 | Ω_I = 6.588026264e-01 | ΔΩ = -1.17564e-05
iter  250 | Ω_I = 6.587907403e-01 | ΔΩ = -1.18861e-05
iter  251 | Ω_I = 6.587797136e-01 | ΔΩ = -1.10267e-05
iter  252 | Ω_I = 6.587671219e-01 | ΔΩ = -1.25916e-05
iter  253 | Ω_I = 6.587572683e-01 | ΔΩ = -9.85368e-06
iter  254 | Ω_I = 6.587476467e-01 | ΔΩ = -9.62154e-06
iter  255 | Ω_I = 6.587351223e-01 | ΔΩ = -1.25244e-05
iter  256 | Ω_I = 6.587236596e-01 | ΔΩ = -1.14627e-05
iter  257 | Ω_I = 6.587129779e-01 | ΔΩ = -1.06817e-05
iter  258 | Ω_I = 6.587008962e-01 | ΔΩ = -1.20817e-05
iter  259 | Ω_I = 6.586929626e-01 | ΔΩ = -7.93359e-06
iter  260 | Ω_I = 6.586823866e-01 | ΔΩ = -1.05760e-05
iter  261 | Ω_I = 6.586714048e-01 | ΔΩ = -1.09818e-05
iter  262 | Ω_I = 6.586634795e-01 | ΔΩ = -7.92524e-06
iter  263 | Ω_I = 6.586534088e-01 | ΔΩ = -1.00708e-05
iter  264 | Ω_I = 6.586440104e-01 | ΔΩ = -9.39839e-06
iter  265 | Ω_I = 6.586352885e-01 | ΔΩ = -8.72190e-06
iter  266 | Ω_I = 6.586200229e-01 | ΔΩ = -1.52655e-05
iter  267 | Ω_I = 6.586137644e-01 | ΔΩ = -6.25849e-06
iter  268 | Ω_I = 6.586028502e-01 | ΔΩ = -1.09142e-05
iter  269 | Ω_I = 6.585957724e-01 | ΔΩ = -7.07786e-06
iter  270 | Ω_I = 6.585847681e-01 | ΔΩ = -1.10043e-05
iter  271 | Ω_I = 6.585763662e-01 | ΔΩ = -8.40181e-06
iter  272 | Ω_I = 6.585674966e-01 | ΔΩ = -8.86967e-06
iter  273 | Ω_I = 6.585566632e-01 | ΔΩ = -1.08334e-05
iter  274 | Ω_I = 6.585485556e-01 | ΔΩ = -8.10760e-06
iter  275 | Ω_I = 6.585430333e-01 | ΔΩ = -5.52226e-06
iter  276 | Ω_I = 6.585337555e-01 | ΔΩ = -9.27783e-06
iter  277 | Ω_I = 6.585212934e-01 | ΔΩ = -1.24621e-05
iter  278 | Ω_I = 6.585164560e-01 | ΔΩ = -4.83740e-06
iter  279 | Ω_I = 6.585086461e-01 | ΔΩ = -7.80987e-06
iter  280 | Ω_I = 6.584994875e-01 | ΔΩ = -9.15864e-06
iter  281 | Ω_I = 6.584901508e-01 | ΔΩ = -9.33670e-06
iter  282 | Ω_I = 6.584879173e-01 | ΔΩ = -2.23351e-06
iter  283 | Ω_I = 6.584794414e-01 | ΔΩ = -8.47592e-06
iter  284 | Ω_I = 6.584727800e-01 | ΔΩ = -6.66137e-06
iter  285 | Ω_I = 6.584643014e-01 | ΔΩ = -8.47862e-06
iter  286 | Ω_I = 6.584590775e-01 | ΔΩ = -5.22384e-06
iter  287 | Ω_I = 6.584501586e-01 | ΔΩ = -8.91887e-06
iter  288 | Ω_I = 6.584446071e-01 | ΔΩ = -5.55157e-06
iter  289 | Ω_I = 6.584364625e-01 | ΔΩ = -8.14456e-06
iter  290 | Ω_I = 6.584286081e-01 | ΔΩ = -7.85437e-06
iter  291 | Ω_I = 6.584238953e-01 | ΔΩ = -4.71283e-06
iter  292 | Ω_I = 6.584182232e-01 | ΔΩ = -5.67208e-06
iter  293 | Ω_I = 6.584079885e-01 | ΔΩ = -1.02348e-05
iter  294 | Ω_I = 6.584006252e-01 | ΔΩ = -7.36330e-06
iter  295 | Ω_I = 6.583951526e-01 | ΔΩ = -5.47256e-06
iter  296 | Ω_I = 6.583896892e-01 | ΔΩ = -5.46335e-06
iter  297 | Ω_I = 6.583812488e-01 | ΔΩ = -8.44043e-06
iter  298 | Ω_I = 6.583775428e-01 | ΔΩ = -3.70598e-06
iter  299 | Ω_I = 6.583747619e-01 | ΔΩ = -2.78095e-06
iter  300 | Ω_I = 6.583676574e-01 | ΔΩ = -7.10444e-06
iter  301 | Ω_I = 6.583599135e-01 | ΔΩ = -7.74398e-06
iter  302 | Ω_I = 6.583521046e-01 | ΔΩ = -7.80889e-06
iter  303 | Ω_I = 6.583493122e-01 | ΔΩ = -2.79235e-06
iter  304 | Ω_I = 6.583398444e-01 | ΔΩ = -9.46788e-06
iter  305 | Ω_I = 6.583393483e-01 | ΔΩ = -4.96091e-07
iter  306 | Ω_I = 6.583308594e-01 | ΔΩ = -8.48890e-06
iter  307 | Ω_I = 6.583251343e-01 | ΔΩ = -5.72511e-06
iter  308 | Ω_I = 6.583194526e-01 | ΔΩ = -5.68168e-06
iter  309 | Ω_I = 6.583154758e-01 | ΔΩ = -3.97680e-06
iter  310 | Ω_I = 6.583110432e-01 | ΔΩ = -4.43257e-06
iter  311 | Ω_I = 6.583027282e-01 | ΔΩ = -8.31496e-06
iter  312 | Ω_I = 6.582995171e-01 | ΔΩ = -3.21115e-06
iter  313 | Ω_I = 6.582952406e-01 | ΔΩ = -4.27648e-06
iter  314 | Ω_I = 6.582889421e-01 | ΔΩ = -6.29849e-06
iter  315 | Ω_I = 6.582839137e-01 | ΔΩ = -5.02844e-06
iter  316 | Ω_I = 6.582793978e-01 | ΔΩ = -4.51583e-06
iter  317 | Ω_I = 6.582736517e-01 | ΔΩ = -5.74616e-06
iter  318 | Ω_I = 6.582694238e-01 | ΔΩ = -4.22785e-06
iter  319 | Ω_I = 6.582664620e-01 | ΔΩ = -2.96184e-06
iter  320 | Ω_I = 6.582594341e-01 | ΔΩ = -7.02792e-06
iter  321 | Ω_I = 6.582548180e-01 | ΔΩ = -4.61607e-06
iter  322 | Ω_I = 6.582512284e-01 | ΔΩ = -3.58964e-06
iter  323 | Ω_I = 6.582446067e-01 | ΔΩ = -6.62162e-06
iter  324 | Ω_I = 6.582422315e-01 | ΔΩ = -2.37526e-06
iter  325 | Ω_I = 6.582387705e-01 | ΔΩ = -3.46094e-06
iter  326 | Ω_I = 6.582359587e-01 | ΔΩ = -2.81180e-06
iter  327 | Ω_I = 6.582299553e-01 | ΔΩ = -6.00347e-06
iter  328 | Ω_I = 6.582288377e-01 | ΔΩ = -1.11756e-06
iter  329 | Ω_I = 6.582233581e-01 | ΔΩ = -5.47963e-06
iter  330 | Ω_I = 6.582189566e-01 | ΔΩ = -4.40148e-06
iter  331 | Ω_I = 6.582130775e-01 | ΔΩ = -5.87915e-06
iter  332 | Ω_I = 6.582098220e-01 | ΔΩ = -3.25544e-06
iter  333 | Ω_I = 6.582074082e-01 | ΔΩ = -2.41386e-06
iter  334 | Ω_I = 6.582024109e-01 | ΔΩ = -4.99722e-06
iter  335 | Ω_I = 6.582005591e-01 | ΔΩ = -1.85180e-06
iter  336 | Ω_I = 6.581965879e-01 | ΔΩ = -3.97122e-06
iter  337 | Ω_I = 6.581933033e-01 | ΔΩ = -3.28465e-06
iter  338 | Ω_I = 6.581881923e-01 | ΔΩ = -5.11098e-06
iter  339 | Ω_I = 6.581853188e-01 | ΔΩ = -2.87350e-06
iter  340 | Ω_I = 6.581811619e-01 | ΔΩ = -4.15684e-06
iter  341 | Ω_I = 6.581761564e-01 | ΔΩ = -5.00552e-06
iter  342 | Ω_I = 6.581752030e-01 | ΔΩ = -9.53447e-07
iter  343 | Ω_I = 6.581699533e-01 | ΔΩ = -5.24968e-06
iter  344 | Ω_I = 6.581660414e-01 | ΔΩ = -3.91190e-06
iter  345 | Ω_I = 6.581629771e-01 | ΔΩ = -3.06430e-06
iter  346 | Ω_I = 6.581573082e-01 | ΔΩ = -5.66894e-06
iter  347 | Ω_I = 6.581516145e-01 | ΔΩ = -5.69370e-06
iter  348 | Ω_I = 6.581523173e-01 | ΔΩ = 7.02803e-07
Warning: Ω_I is increasing. Reducing beta to 0.99.
iter  349 | Ω_I = 6.581475437e-01 | ΔΩ = -4.77354e-06
iter  350 | Ω_I = 6.581446307e-01 | ΔΩ = -2.91298e-06
iter  351 | Ω_I = 6.581419956e-01 | ΔΩ = -2.63510e-06
iter  352 | Ω_I = 6.581382541e-01 | ΔΩ = -3.74150e-06
iter  353 | Ω_I = 6.581367494e-01 | ΔΩ = -1.50476e-06
iter  354 | Ω_I = 6.581350756e-01 | ΔΩ = -1.67381e-06
iter  355 | Ω_I = 6.581293514e-01 | ΔΩ = -5.72416e-06
iter  356 | Ω_I = 6.581273438e-01 | ΔΩ = -2.00757e-06
iter  357 | Ω_I = 6.581251354e-01 | ΔΩ = -2.20846e-06
iter  358 | Ω_I = 6.581222858e-01 | ΔΩ = -2.84954e-06
iter  359 | Ω_I = 6.581186063e-01 | ΔΩ = -3.67951e-06
iter  360 | Ω_I = 6.581174050e-01 | ΔΩ = -1.20137e-06
iter  361 | Ω_I = 6.581126219e-01 | ΔΩ = -4.78310e-06
iter  362 | Ω_I = 6.581105755e-01 | ΔΩ = -2.04635e-06
iter  363 | Ω_I = 6.581068163e-01 | ΔΩ = -3.75923e-06
iter  364 | Ω_I = 6.581045937e-01 | ΔΩ = -2.22260e-06
iter  365 | Ω_I = 6.581014644e-01 | ΔΩ = -3.12924e-06
iter  366 | Ω_I = 6.580976768e-01 | ΔΩ = -3.78769e-06
iter  367 | Ω_I = 6.580980028e-01 | ΔΩ = 3.26094e-07
Warning: Ω_I is increasing. Reducing beta to 0.98.
iter  368 | Ω_I = 6.580945955e-01 | ΔΩ = -3.40731e-06
iter  369 | Ω_I = 6.580872062e-01 | ΔΩ = -7.38934e-06
iter  370 | Ω_I = 6.580871793e-01 | ΔΩ = -2.68831e-08
iter  371 | Ω_I = 6.580886804e-01 | ΔΩ = 1.50109e-06
Warning: Ω_I is increasing. Reducing beta to 0.97.
iter  372 | Ω_I = 6.580842222e-01 | ΔΩ = -4.45821e-06
iter  373 | Ω_I = 6.580793141e-01 | ΔΩ = -4.90809e-06
iter  374 | Ω_I = 6.580758372e-01 | ΔΩ = -3.47694e-06
iter  375 | Ω_I = 6.580765237e-01 | ΔΩ = 6.86529e-07
Warning: Ω_I is increasing. Reducing beta to 0.96.
iter  376 | Ω_I = 6.580722882e-01 | ΔΩ = -4.23545e-06
iter  377 | Ω_I = 6.580683063e-01 | ΔΩ = -3.98191e-06
iter  378 | Ω_I = 6.580681879e-01 | ΔΩ = -1.18455e-07
iter  379 | Ω_I = 6.580640650e-01 | ΔΩ = -4.12287e-06
iter  380 | Ω_I = 6.580636029e-01 | ΔΩ = -4.62146e-07
iter  381 | Ω_I = 6.580613121e-01 | ΔΩ = -2.29078e-06
iter  382 | Ω_I = 6.580578728e-01 | ΔΩ = -3.43925e-06
iter  383 | Ω_I = 6.580558289e-01 | ΔΩ = -2.04396e-06
iter  384 | Ω_I = 6.580550676e-01 | ΔΩ = -7.61294e-07
iter  385 | Ω_I = 6.580519957e-01 | ΔΩ = -3.07184e-06
iter  386 | Ω_I = 6.580503162e-01 | ΔΩ = -1.67950e-06
iter  387 | Ω_I = 6.580492152e-01 | ΔΩ = -1.10109e-06
iter  388 | Ω_I = 6.580480577e-01 | ΔΩ = -1.15745e-06
iter  389 | Ω_I = 6.580456081e-01 | ΔΩ = -2.44959e-06
iter  390 | Ω_I = 6.580421025e-01 | ΔΩ = -3.50563e-06
iter  391 | Ω_I = 6.580403382e-01 | ΔΩ = -1.76432e-06
iter  392 | Ω_I = 6.580366739e-01 | ΔΩ = -3.66428e-06
iter  393 | Ω_I = 6.580330963e-01 | ΔΩ = -3.57760e-06
iter  394 | Ω_I = 6.580339785e-01 | ΔΩ = 8.82262e-07
Warning: Ω_I is increasing. Reducing beta to 0.95.
iter  395 | Ω_I = 6.580320974e-01 | ΔΩ = -1.88117e-06
iter  396 | Ω_I = 6.580283203e-01 | ΔΩ = -3.77707e-06
iter  397 | Ω_I = 6.580272266e-01 | ΔΩ = -1.09371e-06
iter  398 | Ω_I = 6.580293680e-01 | ΔΩ = 2.14143e-06
Warning: Ω_I is increasing. Reducing beta to 0.94.
iter  399 | Ω_I = 6.580258044e-01 | ΔΩ = -3.56362e-06
iter  400 | Ω_I = 6.580219998e-01 | ΔΩ = -3.80459e-06
iter  401 | Ω_I = 6.580192676e-01 | ΔΩ = -2.73217e-06
iter  402 | Ω_I = 6.580174079e-01 | ΔΩ = -1.85971e-06
iter  403 | Ω_I = 6.580163456e-01 | ΔΩ = -1.06235e-06
iter  404 | Ω_I = 6.580141395e-01 | ΔΩ = -2.20611e-06
iter  405 | Ω_I = 6.580124204e-01 | ΔΩ = -1.71903e-06
iter  406 | Ω_I = 6.580125970e-01 | ΔΩ = 1.76569e-07
Warning: Ω_I is increasing. Reducing beta to 0.9299999999999999.
iter  407 | Ω_I = 6.580092953e-01 | ΔΩ = -3.30172e-06
iter  408 | Ω_I = 6.580084620e-01 | ΔΩ = -8.33272e-07
iter  409 | Ω_I = 6.580071689e-01 | ΔΩ = -1.29315e-06
iter  410 | Ω_I = 6.580045359e-01 | ΔΩ = -2.63296e-06
iter  411 | Ω_I = 6.580054582e-01 | ΔΩ = 9.22276e-07
Warning: Ω_I is increasing. Reducing beta to 0.9199999999999999.
iter  412 | Ω_I = 6.580010596e-01 | ΔΩ = -4.39857e-06
iter  413 | Ω_I = 6.580004705e-01 | ΔΩ = -5.89078e-07
iter  414 | Ω_I = 6.579994625e-01 | ΔΩ = -1.00806e-06
iter  415 | Ω_I = 6.579964433e-01 | ΔΩ = -3.01918e-06
iter  416 | Ω_I = 6.579941812e-01 | ΔΩ = -2.26207e-06
iter  417 | Ω_I = 6.579942199e-01 | ΔΩ = 3.86630e-08
Warning: Ω_I is increasing. Reducing beta to 0.9099999999999999.
iter  418 | Ω_I = 6.579926124e-01 | ΔΩ = -1.60746e-06
iter  419 | Ω_I = 6.579910473e-01 | ΔΩ = -1.56513e-06
iter  420 | Ω_I = 6.579906054e-01 | ΔΩ = -4.41926e-07
iter  421 | Ω_I = 6.579889593e-01 | ΔΩ = -1.64605e-06
iter  422 | Ω_I = 6.579868958e-01 | ΔΩ = -2.06350e-06
iter  423 | Ω_I = 6.579834481e-01 | ΔΩ = -3.44773e-06
iter  424 | Ω_I = 6.579838898e-01 | ΔΩ = 4.41742e-07
Warning: Ω_I is increasing. Reducing beta to 0.8999999999999999.
iter  425 | Ω_I = 6.579811459e-01 | ΔΩ = -2.74389e-06
iter  426 | Ω_I = 6.579772637e-01 | ΔΩ = -3.88219e-06
iter  427 | Ω_I = 6.579776190e-01 | ΔΩ = 3.55294e-07
Warning: Ω_I is increasing. Reducing beta to 0.8899999999999999.
iter  428 | Ω_I = 6.579754075e-01 | ΔΩ = -2.21153e-06
iter  429 | Ω_I = 6.579740492e-01 | ΔΩ = -1.35830e-06
iter  430 | Ω_I = 6.579749867e-01 | ΔΩ = 9.37444e-07
Warning: Ω_I is increasing. Reducing beta to 0.8799999999999999.
iter  431 | Ω_I = 6.579749001e-01 | ΔΩ = -8.65569e-08
iter  432 | Ω_I = 6.579710508e-01 | ΔΩ = -3.84928e-06
iter  433 | Ω_I = 6.579729035e-01 | ΔΩ = 1.85263e-06
Warning: Ω_I is increasing. Reducing beta to 0.8699999999999999.
iter  434 | Ω_I = 6.579695449e-01 | ΔΩ = -3.35854e-06
iter  435 | Ω_I = 6.579657480e-01 | ΔΩ = -3.79690e-06
iter  436 | Ω_I = 6.579674218e-01 | ΔΩ = 1.67378e-06
Warning: Ω_I is increasing. Reducing beta to 0.8599999999999999.
iter  437 | Ω_I = 6.579669736e-01 | ΔΩ = -4.48177e-07
iter  438 | Ω_I = 6.579632109e-01 | ΔΩ = -3.76273e-06
iter  439 | Ω_I = 6.579624715e-01 | ΔΩ = -7.39378e-07
iter  440 | Ω_I = 6.579629357e-01 | ΔΩ = 4.64174e-07
Warning: Ω_I is increasing. Reducing beta to 0.8499999999999999.
iter  441 | Ω_I = 6.579605401e-01 | ΔΩ = -2.39563e-06
iter  442 | Ω_I = 6.579593430e-01 | ΔΩ = -1.19702e-06
iter  443 | Ω_I = 6.579604502e-01 | ΔΩ = 1.10715e-06
Warning: Ω_I is increasing. Reducing beta to 0.8399999999999999.
iter  444 | Ω_I = 6.579589954e-01 | ΔΩ = -1.45481e-06
iter  445 | Ω_I = 6.579569119e-01 | ΔΩ = -2.08346e-06
iter  446 | Ω_I = 6.579562485e-01 | ΔΩ = -6.63363e-07
iter  447 | Ω_I = 6.579574962e-01 | ΔΩ = 1.24768e-06
Warning: Ω_I is increasing. Reducing beta to 0.8299999999999998.
iter  448 | Ω_I = 6.579526939e-01 | ΔΩ = -4.80236e-06
iter  449 | Ω_I = 6.579526356e-01 | ΔΩ = -5.83090e-08
iter  450 | Ω_I = 6.579521412e-01 | ΔΩ = -4.94361e-07
iter  451 | Ω_I = 6.579512099e-01 | ΔΩ = -9.31295e-07
iter  452 | Ω_I = 6.579479036e-01 | ΔΩ = -3.30627e-06
iter  453 | Ω_I = 6.579469894e-01 | ΔΩ = -9.14185e-07
iter  454 | Ω_I = 6.579478711e-01 | ΔΩ = 8.81635e-07
Warning: Ω_I is increasing. Reducing beta to 0.8199999999999998.
iter  455 | Ω_I = 6.579472793e-01 | ΔΩ = -5.91746e-07
iter  456 | Ω_I = 6.579457067e-01 | ΔΩ = -1.57259e-06
iter  457 | Ω_I = 6.579451580e-01 | ΔΩ = -5.48727e-07
iter  458 | Ω_I = 6.579429883e-01 | ΔΩ = -2.16968e-06
iter  459 | Ω_I = 6.579441387e-01 | ΔΩ = 1.15031e-06
Warning: Ω_I is increasing. Reducing beta to 0.8099999999999998.
iter  460 | Ω_I = 6.579427837e-01 | ΔΩ = -1.35497e-06
iter  461 | Ω_I = 6.579414132e-01 | ΔΩ = -1.37047e-06
iter  462 | Ω_I = 6.579414029e-01 | ΔΩ = -1.03181e-08
iter  463 | Ω_I = 6.579410796e-01 | ΔΩ = -3.23332e-07
iter  464 | Ω_I = 6.579388979e-01 | ΔΩ = -2.18170e-06
iter  465 | Ω_I = 6.579369629e-01 | ΔΩ = -1.93498e-06
iter  466 | Ω_I = 6.579354661e-01 | ΔΩ = -1.49675e-06
iter  467 | Ω_I = 6.579345072e-01 | ΔΩ = -9.58931e-07
iter  468 | Ω_I = 6.579345580e-01 | ΔΩ = 5.08433e-08
Warning: Ω_I is increasing. Reducing beta to 0.7999999999999998.
iter  469 | Ω_I = 6.579323715e-01 | ΔΩ = -2.18658e-06
iter  470 | Ω_I = 6.579323549e-01 | ΔΩ = -1.65452e-08
iter  471 | Ω_I = 6.579320477e-01 | ΔΩ = -3.07165e-07
iter  472 | Ω_I = 6.579292601e-01 | ΔΩ = -2.78765e-06
iter  473 | Ω_I = 6.579283158e-01 | ΔΩ = -9.44314e-07
iter  474 | Ω_I = 6.579297571e-01 | ΔΩ = 1.44133e-06
Warning: Ω_I is increasing. Reducing beta to 0.7899999999999998.
iter  475 | Ω_I = 6.579293338e-01 | ΔΩ = -4.23306e-07
iter  476 | Ω_I = 6.579292743e-01 | ΔΩ = -5.95444e-08
iter  477 | Ω_I = 6.579270819e-01 | ΔΩ = -2.19238e-06
iter  478 | Ω_I = 6.579259470e-01 | ΔΩ = -1.13488e-06
iter  479 | Ω_I = 6.579261709e-01 | ΔΩ = 2.23909e-07
Warning: Ω_I is increasing. Reducing beta to 0.7799999999999998.
iter  480 | Ω_I = 6.579219426e-01 | ΔΩ = -4.22833e-06
iter  481 | Ω_I = 6.579230637e-01 | ΔΩ = 1.12106e-06
Warning: Ω_I is increasing. Reducing beta to 0.7699999999999998.
iter  482 | Ω_I = 6.579204595e-01 | ΔΩ = -2.60416e-06
iter  483 | Ω_I = 6.579205340e-01 | ΔΩ = 7.44950e-08
Warning: Ω_I is increasing. Reducing beta to 0.7599999999999998.
iter  484 | Ω_I = 6.579201643e-01 | ΔΩ = -3.69733e-07
iter  485 | Ω_I = 6.579193915e-01 | ΔΩ = -7.72781e-07
iter  486 | Ω_I = 6.579201893e-01 | ΔΩ = 7.97843e-07
Warning: Ω_I is increasing. Reducing beta to 0.7499999999999998.
iter  487 | Ω_I = 6.579184851e-01 | ΔΩ = -1.70420e-06
iter  488 | Ω_I = 6.579152362e-01 | ΔΩ = -3.24889e-06
iter  489 | Ω_I = 6.579143134e-01 | ΔΩ = -9.22876e-07
iter  490 | Ω_I = 6.579141995e-01 | ΔΩ = -1.13873e-07
iter  491 | Ω_I = 6.579125211e-01 | ΔΩ = -1.67834e-06
iter  492 | Ω_I = 6.579129457e-01 | ΔΩ = 4.24545e-07
Warning: Ω_I is increasing. Reducing beta to 0.7399999999999998.
iter  493 | Ω_I = 6.579128766e-01 | ΔΩ = -6.90918e-08
iter  494 | Ω_I = 6.579100642e-01 | ΔΩ = -2.81242e-06
iter  495 | Ω_I = 6.579122456e-01 | ΔΩ = 2.18144e-06
Warning: Ω_I is increasing. Reducing beta to 0.7299999999999998.
iter  496 | Ω_I = 6.579106570e-01 | ΔΩ = -1.58859e-06
iter  497 | Ω_I = 6.579086571e-01 | ΔΩ = -1.99994e-06
iter  498 | Ω_I = 6.579088661e-01 | ΔΩ = 2.08980e-07
Warning: Ω_I is increasing. Reducing beta to 0.7199999999999998.
iter  499 | Ω_I = 6.579074860e-01 | ΔΩ = -1.38007e-06
WF.info()
Wannier Function Report
===================================================================
WF 1: center = [0.15002752, -0.60620957]  Omega     = 234.66083086
WF 2: center = [-0.40002465, -0.11548027]  Omega     = 235.62189606
WF 3: center = [0.04999661, -0.54848076]  Omega     = 222.02088753
Sum : center = [-0.20000052, -1.27017059]  Omega tot = 692.30361445
-------------------------------------------------------------------
Omega I   = 0.65789262
Omega D   = 686.33706926
Omega OD  = 5.30865257
Omega tot = 692.30361445

Maximal localization#

To obtain maximally localized Wannier functions, we follow this with another projection to initialize a smooth gauge, then maximal localization.

  • Note we must pass the flag tilde=True to indicate we are projecting the trial wavefunctions onto the tilde states and not the energy eigenstates

WF.project(use_tilde=True)
WF.info()
Wannier Function Report
===============================================================
WF 1: center = [0.48238924, 0.27850754]  Omega     = 0.23471731
WF 2: center = [1.00000001, 1.17503572]  Omega     = 0.23471689
WF 3: center = [1.51761076, 0.27850756]  Omega     = 0.23471889
Sum : center = [3.00000001, 1.73205081]  Omega tot = 0.70415309
---------------------------------------------------------------
Omega I   = 0.65789262
Omega D   = 0.00747840
Omega OD  = 0.03878207
Omega tot = 0.70415309
WF.maxloc(alpha=1 / 2, max_iter=1000, tol=1e-10, grad_min=1e-10, verbose=True)

Hide code cell output

iter    0 | Ω_tilde = 4.608330209e-02 | ΔΩ = -1.77166e-04 | ‖∇‖ = 2.31951e-05
iter    1 | Ω_tilde = 4.591826516e-02 | ΔΩ = -1.65037e-04 | ‖∇‖ = 2.30942e-05
iter    2 | Ω_tilde = 4.576414641e-02 | ΔΩ = -1.54119e-04 | ‖∇‖ = 2.29938e-05
iter    3 | Ω_tilde = 4.561991554e-02 | ΔΩ = -1.44231e-04 | ‖∇‖ = 2.28939e-05
iter    4 | Ω_tilde = 4.548468487e-02 | ΔΩ = -1.35231e-04 | ‖∇‖ = 2.27946e-05
iter    5 | Ω_tilde = 4.535768146e-02 | ΔΩ = -1.27003e-04 | ‖∇‖ = 2.26958e-05
iter    6 | Ω_tilde = 4.523822639e-02 | ΔΩ = -1.19455e-04 | ‖∇‖ = 2.25976e-05
iter    7 | Ω_tilde = 4.512571884e-02 | ΔΩ = -1.12508e-04 | ‖∇‖ = 2.25000e-05
iter    8 | Ω_tilde = 4.501962373e-02 | ΔΩ = -1.06095e-04 | ‖∇‖ = 2.24029e-05
iter    9 | Ω_tilde = 4.491946186e-02 | ΔΩ = -1.00162e-04 | ‖∇‖ = 2.23064e-05
iter   10 | Ω_tilde = 4.482480200e-02 | ΔΩ = -9.46599e-05 | ‖∇‖ = 2.22105e-05
iter   11 | Ω_tilde = 4.473525445e-02 | ΔΩ = -8.95476e-05 | ‖∇‖ = 2.21151e-05
iter   12 | Ω_tilde = 4.465046565e-02 | ΔΩ = -8.47888e-05 | ‖∇‖ = 2.20202e-05
iter   13 | Ω_tilde = 4.457011381e-02 | ΔΩ = -8.03518e-05 | ‖∇‖ = 2.19259e-05
iter   14 | Ω_tilde = 4.449390514e-02 | ΔΩ = -7.62087e-05 | ‖∇‖ = 2.18321e-05
iter   15 | Ω_tilde = 4.442157066e-02 | ΔΩ = -7.23345e-05 | ‖∇‖ = 2.17388e-05
iter   16 | Ω_tilde = 4.435286353e-02 | ΔΩ = -6.87071e-05 | ‖∇‖ = 2.16460e-05
iter   17 | Ω_tilde = 4.428755671e-02 | ΔΩ = -6.53068e-05 | ‖∇‖ = 2.15537e-05
iter   18 | Ω_tilde = 4.422544091e-02 | ΔΩ = -6.21158e-05 | ‖∇‖ = 2.14618e-05
iter   19 | Ω_tilde = 4.416632288e-02 | ΔΩ = -5.91180e-05 | ‖∇‖ = 2.13705e-05
iter   20 | Ω_tilde = 4.411002386e-02 | ΔΩ = -5.62990e-05 | ‖∇‖ = 2.12795e-05
iter   21 | Ω_tilde = 4.405637818e-02 | ΔΩ = -5.36457e-05 | ‖∇‖ = 2.11891e-05
iter   22 | Ω_tilde = 4.400523216e-02 | ΔΩ = -5.11460e-05 | ‖∇‖ = 2.10991e-05
iter   23 | Ω_tilde = 4.395644294e-02 | ΔΩ = -4.87892e-05 | ‖∇‖ = 2.10096e-05
iter   24 | Ω_tilde = 4.390987764e-02 | ΔΩ = -4.65653e-05 | ‖∇‖ = 2.09205e-05
iter   25 | Ω_tilde = 4.386541245e-02 | ΔΩ = -4.44652e-05 | ‖∇‖ = 2.08318e-05
iter   26 | Ω_tilde = 4.382293189e-02 | ΔΩ = -4.24806e-05 | ‖∇‖ = 2.07436e-05
iter   27 | Ω_tilde = 4.378232814e-02 | ΔΩ = -4.06037e-05 | ‖∇‖ = 2.06558e-05
iter   28 | Ω_tilde = 4.374350043e-02 | ΔΩ = -3.88277e-05 | ‖∇‖ = 2.05684e-05
iter   29 | Ω_tilde = 4.370635449e-02 | ΔΩ = -3.71459e-05 | ‖∇‖ = 2.04814e-05
iter   30 | Ω_tilde = 4.367080202e-02 | ΔΩ = -3.55525e-05 | ‖∇‖ = 2.03948e-05
iter   31 | Ω_tilde = 4.363676028e-02 | ΔΩ = -3.40417e-05 | ‖∇‖ = 2.03087e-05
iter   32 | Ω_tilde = 4.360415165e-02 | ΔΩ = -3.26086e-05 | ‖∇‖ = 2.02229e-05
iter   33 | Ω_tilde = 4.357290326e-02 | ΔΩ = -3.12484e-05 | ‖∇‖ = 2.01376e-05
iter   34 | Ω_tilde = 4.354294664e-02 | ΔΩ = -2.99566e-05 | ‖∇‖ = 2.00526e-05
iter   35 | Ω_tilde = 4.351421741e-02 | ΔΩ = -2.87292e-05 | ‖∇‖ = 1.99680e-05
iter   36 | Ω_tilde = 4.348665501e-02 | ΔΩ = -2.75624e-05 | ‖∇‖ = 1.98839e-05
iter   37 | Ω_tilde = 4.346020241e-02 | ΔΩ = -2.64526e-05 | ‖∇‖ = 1.98001e-05
iter   38 | Ω_tilde = 4.343480587e-02 | ΔΩ = -2.53965e-05 | ‖∇‖ = 1.97167e-05
iter   39 | Ω_tilde = 4.341041473e-02 | ΔΩ = -2.43911e-05 | ‖∇‖ = 1.96336e-05
iter   40 | Ω_tilde = 4.338698122e-02 | ΔΩ = -2.34335e-05 | ‖∇‖ = 1.95510e-05
iter   41 | Ω_tilde = 4.336446022e-02 | ΔΩ = -2.25210e-05 | ‖∇‖ = 1.94687e-05
iter   42 | Ω_tilde = 4.334280914e-02 | ΔΩ = -2.16511e-05 | ‖∇‖ = 1.93868e-05
iter   43 | Ω_tilde = 4.332198771e-02 | ΔΩ = -2.08214e-05 | ‖∇‖ = 1.93052e-05
iter   44 | Ω_tilde = 4.330195788e-02 | ΔΩ = -2.00298e-05 | ‖∇‖ = 1.92240e-05
iter   45 | Ω_tilde = 4.328268363e-02 | ΔΩ = -1.92742e-05 | ‖∇‖ = 1.91432e-05
iter   46 | Ω_tilde = 4.326413089e-02 | ΔΩ = -1.85527e-05 | ‖∇‖ = 1.90628e-05
iter   47 | Ω_tilde = 4.324626738e-02 | ΔΩ = -1.78635e-05 | ‖∇‖ = 1.89826e-05
iter   48 | Ω_tilde = 4.322906254e-02 | ΔΩ = -1.72048e-05 | ‖∇‖ = 1.89029e-05
iter   49 | Ω_tilde = 4.321248738e-02 | ΔΩ = -1.65752e-05 | ‖∇‖ = 1.88235e-05
iter   50 | Ω_tilde = 4.319651442e-02 | ΔΩ = -1.59730e-05 | ‖∇‖ = 1.87444e-05
iter   51 | Ω_tilde = 4.318111759e-02 | ΔΩ = -1.53968e-05 | ‖∇‖ = 1.86657e-05
iter   52 | Ω_tilde = 4.316627215e-02 | ΔΩ = -1.48454e-05 | ‖∇‖ = 1.85874e-05
iter   53 | Ω_tilde = 4.315195460e-02 | ΔΩ = -1.43176e-05 | ‖∇‖ = 1.85093e-05
iter   54 | Ω_tilde = 4.313814261e-02 | ΔΩ = -1.38120e-05 | ‖∇‖ = 1.84316e-05
iter   55 | Ω_tilde = 4.312481498e-02 | ΔΩ = -1.33276e-05 | ‖∇‖ = 1.83543e-05
iter   56 | Ω_tilde = 4.311195155e-02 | ΔΩ = -1.28634e-05 | ‖∇‖ = 1.82773e-05
iter   57 | Ω_tilde = 4.309953313e-02 | ΔΩ = -1.24184e-05 | ‖∇‖ = 1.82006e-05
iter   58 | Ω_tilde = 4.308754147e-02 | ΔΩ = -1.19917e-05 | ‖∇‖ = 1.81242e-05
iter   59 | Ω_tilde = 4.307595921e-02 | ΔΩ = -1.15823e-05 | ‖∇‖ = 1.80482e-05
iter   60 | Ω_tilde = 4.306476980e-02 | ΔΩ = -1.11894e-05 | ‖∇‖ = 1.79725e-05
iter   61 | Ω_tilde = 4.305395750e-02 | ΔΩ = -1.08123e-05 | ‖∇‖ = 1.78972e-05
iter   62 | Ω_tilde = 4.304350730e-02 | ΔΩ = -1.04502e-05 | ‖∇‖ = 1.78221e-05
iter   63 | Ω_tilde = 4.303340487e-02 | ΔΩ = -1.01024e-05 | ‖∇‖ = 1.77474e-05
iter   64 | Ω_tilde = 4.302363659e-02 | ΔΩ = -9.76828e-06 | ‖∇‖ = 1.76730e-05
iter   65 | Ω_tilde = 4.301418943e-02 | ΔΩ = -9.44716e-06 | ‖∇‖ = 1.75989e-05
iter   66 | Ω_tilde = 4.300505098e-02 | ΔΩ = -9.13845e-06 | ‖∇‖ = 1.75251e-05
iter   67 | Ω_tilde = 4.299620938e-02 | ΔΩ = -8.84160e-06 | ‖∇‖ = 1.74516e-05
iter   68 | Ω_tilde = 4.298765331e-02 | ΔΩ = -8.55607e-06 | ‖∇‖ = 1.73785e-05
iter   69 | Ω_tilde = 4.297937197e-02 | ΔΩ = -8.28135e-06 | ‖∇‖ = 1.73057e-05
iter   70 | Ω_tilde = 4.297135501e-02 | ΔΩ = -8.01695e-06 | ‖∇‖ = 1.72331e-05
iter   71 | Ω_tilde = 4.296359257e-02 | ΔΩ = -7.76244e-06 | ‖∇‖ = 1.71609e-05
iter   72 | Ω_tilde = 4.295607521e-02 | ΔΩ = -7.51736e-06 | ‖∇‖ = 1.70890e-05
iter   73 | Ω_tilde = 4.294879390e-02 | ΔΩ = -7.28131e-06 | ‖∇‖ = 1.70174e-05
iter   74 | Ω_tilde = 4.294174001e-02 | ΔΩ = -7.05390e-06 | ‖∇‖ = 1.69460e-05
iter   75 | Ω_tilde = 4.293490525e-02 | ΔΩ = -6.83475e-06 | ‖∇‖ = 1.68750e-05
iter   76 | Ω_tilde = 4.292828174e-02 | ΔΩ = -6.62352e-06 | ‖∇‖ = 1.68043e-05
iter   77 | Ω_tilde = 4.292186188e-02 | ΔΩ = -6.41986e-06 | ‖∇‖ = 1.67339e-05
iter   78 | Ω_tilde = 4.291563842e-02 | ΔΩ = -6.22346e-06 | ‖∇‖ = 1.66638e-05
iter   79 | Ω_tilde = 4.290960440e-02 | ΔΩ = -6.03401e-06 | ‖∇‖ = 1.65939e-05
iter   80 | Ω_tilde = 4.290375317e-02 | ΔΩ = -5.85123e-06 | ‖∇‖ = 1.65244e-05
iter   81 | Ω_tilde = 4.289807833e-02 | ΔΩ = -5.67484e-06 | ‖∇‖ = 1.64552e-05
iter   82 | Ω_tilde = 4.289257376e-02 | ΔΩ = -5.50457e-06 | ‖∇‖ = 1.63862e-05
iter   83 | Ω_tilde = 4.288723358e-02 | ΔΩ = -5.34018e-06 | ‖∇‖ = 1.63175e-05
iter   84 | Ω_tilde = 4.288205215e-02 | ΔΩ = -5.18143e-06 | ‖∇‖ = 1.62492e-05
iter   85 | Ω_tilde = 4.287702406e-02 | ΔΩ = -5.02809e-06 | ‖∇‖ = 1.61811e-05
iter   86 | Ω_tilde = 4.287214411e-02 | ΔΩ = -4.87995e-06 | ‖∇‖ = 1.61133e-05
iter   87 | Ω_tilde = 4.286740731e-02 | ΔΩ = -4.73680e-06 | ‖∇‖ = 1.60457e-05
iter   88 | Ω_tilde = 4.286280888e-02 | ΔΩ = -4.59844e-06 | ‖∇‖ = 1.59785e-05
iter   89 | Ω_tilde = 4.285834420e-02 | ΔΩ = -4.46468e-06 | ‖∇‖ = 1.59115e-05
iter   90 | Ω_tilde = 4.285400884e-02 | ΔΩ = -4.33535e-06 | ‖∇‖ = 1.58449e-05
iter   91 | Ω_tilde = 4.284979856e-02 | ΔΩ = -4.21028e-06 | ‖∇‖ = 1.57784e-05
iter   92 | Ω_tilde = 4.284570927e-02 | ΔΩ = -4.08929e-06 | ‖∇‖ = 1.57123e-05
iter   93 | Ω_tilde = 4.284173702e-02 | ΔΩ = -3.97225e-06 | ‖∇‖ = 1.56465e-05
iter   94 | Ω_tilde = 4.283787804e-02 | ΔΩ = -3.85898e-06 | ‖∇‖ = 1.55809e-05
iter   95 | Ω_tilde = 4.283412868e-02 | ΔΩ = -3.74936e-06 | ‖∇‖ = 1.55156e-05
iter   96 | Ω_tilde = 4.283048542e-02 | ΔΩ = -3.64325e-06 | ‖∇‖ = 1.54505e-05
iter   97 | Ω_tilde = 4.282694491e-02 | ΔΩ = -3.54051e-06 | ‖∇‖ = 1.53858e-05
iter   98 | Ω_tilde = 4.282350388e-02 | ΔΩ = -3.44103e-06 | ‖∇‖ = 1.53213e-05
iter   99 | Ω_tilde = 4.282015920e-02 | ΔΩ = -3.34468e-06 | ‖∇‖ = 1.52570e-05
iter  100 | Ω_tilde = 4.281690786e-02 | ΔΩ = -3.25134e-06 | ‖∇‖ = 1.51931e-05
iter  101 | Ω_tilde = 4.281374695e-02 | ΔΩ = -3.16091e-06 | ‖∇‖ = 1.51294e-05
iter  102 | Ω_tilde = 4.281067366e-02 | ΔΩ = -3.07329e-06 | ‖∇‖ = 1.50659e-05
iter  103 | Ω_tilde = 4.280768529e-02 | ΔΩ = -2.98837e-06 | ‖∇‖ = 1.50028e-05
iter  104 | Ω_tilde = 4.280477924e-02 | ΔΩ = -2.90605e-06 | ‖∇‖ = 1.49399e-05
iter  105 | Ω_tilde = 4.280195299e-02 | ΔΩ = -2.82625e-06 | ‖∇‖ = 1.48772e-05
iter  106 | Ω_tilde = 4.279920411e-02 | ΔΩ = -2.74888e-06 | ‖∇‖ = 1.48148e-05
iter  107 | Ω_tilde = 4.279653026e-02 | ΔΩ = -2.67384e-06 | ‖∇‖ = 1.47527e-05
iter  108 | Ω_tilde = 4.279392919e-02 | ΔΩ = -2.60107e-06 | ‖∇‖ = 1.46908e-05
iter  109 | Ω_tilde = 4.279139871e-02 | ΔΩ = -2.53048e-06 | ‖∇‖ = 1.46292e-05
iter  110 | Ω_tilde = 4.278893672e-02 | ΔΩ = -2.46199e-06 | ‖∇‖ = 1.45678e-05
iter  111 | Ω_tilde = 4.278654117e-02 | ΔΩ = -2.39554e-06 | ‖∇‖ = 1.45067e-05
iter  112 | Ω_tilde = 4.278421012e-02 | ΔΩ = -2.33106e-06 | ‖∇‖ = 1.44459e-05
iter  113 | Ω_tilde = 4.278194164e-02 | ΔΩ = -2.26847e-06 | ‖∇‖ = 1.43852e-05
iter  114 | Ω_tilde = 4.277973392e-02 | ΔΩ = -2.20773e-06 | ‖∇‖ = 1.43249e-05
iter  115 | Ω_tilde = 4.277758517e-02 | ΔΩ = -2.14875e-06 | ‖∇‖ = 1.42648e-05
iter  116 | Ω_tilde = 4.277549367e-02 | ΔΩ = -2.09149e-06 | ‖∇‖ = 1.42049e-05
iter  117 | Ω_tilde = 4.277345778e-02 | ΔΩ = -2.03589e-06 | ‖∇‖ = 1.41453e-05
iter  118 | Ω_tilde = 4.277147589e-02 | ΔΩ = -1.98190e-06 | ‖∇‖ = 1.40860e-05
iter  119 | Ω_tilde = 4.276954644e-02 | ΔΩ = -1.92945e-06 | ‖∇‖ = 1.40268e-05
iter  120 | Ω_tilde = 4.276766792e-02 | ΔΩ = -1.87851e-06 | ‖∇‖ = 1.39680e-05
iter  121 | Ω_tilde = 4.276583890e-02 | ΔΩ = -1.82902e-06 | ‖∇‖ = 1.39093e-05
iter  122 | Ω_tilde = 4.276405796e-02 | ΔΩ = -1.78094e-06 | ‖∇‖ = 1.38509e-05
iter  123 | Ω_tilde = 4.276232374e-02 | ΔΩ = -1.73422e-06 | ‖∇‖ = 1.37928e-05
iter  124 | Ω_tilde = 4.276063491e-02 | ΔΩ = -1.68882e-06 | ‖∇‖ = 1.37349e-05
iter  125 | Ω_tilde = 4.275899021e-02 | ΔΩ = -1.64470e-06 | ‖∇‖ = 1.36772e-05
iter  126 | Ω_tilde = 4.275738840e-02 | ΔΩ = -1.60181e-06 | ‖∇‖ = 1.36198e-05
iter  127 | Ω_tilde = 4.275582828e-02 | ΔΩ = -1.56012e-06 | ‖∇‖ = 1.35626e-05
iter  128 | Ω_tilde = 4.275430868e-02 | ΔΩ = -1.51960e-06 | ‖∇‖ = 1.35057e-05
iter  129 | Ω_tilde = 4.275282849e-02 | ΔΩ = -1.48019e-06 | ‖∇‖ = 1.34489e-05
iter  130 | Ω_tilde = 4.275138661e-02 | ΔΩ = -1.44188e-06 | ‖∇‖ = 1.33924e-05
iter  131 | Ω_tilde = 4.274998198e-02 | ΔΩ = -1.40462e-06 | ‖∇‖ = 1.33362e-05
iter  132 | Ω_tilde = 4.274861359e-02 | ΔΩ = -1.36839e-06 | ‖∇‖ = 1.32802e-05
iter  133 | Ω_tilde = 4.274728044e-02 | ΔΩ = -1.33316e-06 | ‖∇‖ = 1.32244e-05
iter  134 | Ω_tilde = 4.274598155e-02 | ΔΩ = -1.29888e-06 | ‖∇‖ = 1.31688e-05
iter  135 | Ω_tilde = 4.274471601e-02 | ΔΩ = -1.26554e-06 | ‖∇‖ = 1.31135e-05
iter  136 | Ω_tilde = 4.274348291e-02 | ΔΩ = -1.23311e-06 | ‖∇‖ = 1.30584e-05
iter  137 | Ω_tilde = 4.274228136e-02 | ΔΩ = -1.20155e-06 | ‖∇‖ = 1.30035e-05
iter  138 | Ω_tilde = 4.274111051e-02 | ΔΩ = -1.17085e-06 | ‖∇‖ = 1.29489e-05
iter  139 | Ω_tilde = 4.273996953e-02 | ΔΩ = -1.14098e-06 | ‖∇‖ = 1.28945e-05
iter  140 | Ω_tilde = 4.273885763e-02 | ΔΩ = -1.11191e-06 | ‖∇‖ = 1.28403e-05
iter  141 | Ω_tilde = 4.273777401e-02 | ΔΩ = -1.08361e-06 | ‖∇‖ = 1.27863e-05
iter  142 | Ω_tilde = 4.273671793e-02 | ΔΩ = -1.05608e-06 | ‖∇‖ = 1.27326e-05
iter  143 | Ω_tilde = 4.273568865e-02 | ΔΩ = -1.02928e-06 | ‖∇‖ = 1.26790e-05
iter  144 | Ω_tilde = 4.273468545e-02 | ΔΩ = -1.00320e-06 | ‖∇‖ = 1.26257e-05
iter  145 | Ω_tilde = 4.273370764e-02 | ΔΩ = -9.77807e-07 | ‖∇‖ = 1.25726e-05
iter  146 | Ω_tilde = 4.273275455e-02 | ΔΩ = -9.53089e-07 | ‖∇‖ = 1.25198e-05
iter  147 | Ω_tilde = 4.273182553e-02 | ΔΩ = -9.29025e-07 | ‖∇‖ = 1.24671e-05
iter  148 | Ω_tilde = 4.273091993e-02 | ΔΩ = -9.05595e-07 | ‖∇‖ = 1.24147e-05
iter  149 | Ω_tilde = 4.273003715e-02 | ΔΩ = -8.82783e-07 | ‖∇‖ = 1.23625e-05
iter  150 | Ω_tilde = 4.272917658e-02 | ΔΩ = -8.60571e-07 | ‖∇‖ = 1.23105e-05
iter  151 | Ω_tilde = 4.272833764e-02 | ΔΩ = -8.38941e-07 | ‖∇‖ = 1.22587e-05
iter  152 | Ω_tilde = 4.272751976e-02 | ΔΩ = -8.17878e-07 | ‖∇‖ = 1.22072e-05
iter  153 | Ω_tilde = 4.272672239e-02 | ΔΩ = -7.97365e-07 | ‖∇‖ = 1.21558e-05
iter  154 | Ω_tilde = 4.272594501e-02 | ΔΩ = -7.77387e-07 | ‖∇‖ = 1.21047e-05
iter  155 | Ω_tilde = 4.272518708e-02 | ΔΩ = -7.57929e-07 | ‖∇‖ = 1.20537e-05
iter  156 | Ω_tilde = 4.272444810e-02 | ΔΩ = -7.38976e-07 | ‖∇‖ = 1.20030e-05
iter  157 | Ω_tilde = 4.272372759e-02 | ΔΩ = -7.20516e-07 | ‖∇‖ = 1.19525e-05
iter  158 | Ω_tilde = 4.272302505e-02 | ΔΩ = -7.02533e-07 | ‖∇‖ = 1.19022e-05
iter  159 | Ω_tilde = 4.272234004e-02 | ΔΩ = -6.85015e-07 | ‖∇‖ = 1.18521e-05
iter  160 | Ω_tilde = 4.272167209e-02 | ΔΩ = -6.67949e-07 | ‖∇‖ = 1.18023e-05
iter  161 | Ω_tilde = 4.272102077e-02 | ΔΩ = -6.51323e-07 | ‖∇‖ = 1.17526e-05
iter  162 | Ω_tilde = 4.272038564e-02 | ΔΩ = -6.35124e-07 | ‖∇‖ = 1.17031e-05
iter  163 | Ω_tilde = 4.271976630e-02 | ΔΩ = -6.19342e-07 | ‖∇‖ = 1.16538e-05
iter  164 | Ω_tilde = 4.271916234e-02 | ΔΩ = -6.03964e-07 | ‖∇‖ = 1.16048e-05
iter  165 | Ω_tilde = 4.271857336e-02 | ΔΩ = -5.88980e-07 | ‖∇‖ = 1.15559e-05
iter  166 | Ω_tilde = 4.271799898e-02 | ΔΩ = -5.74380e-07 | ‖∇‖ = 1.15073e-05
iter  167 | Ω_tilde = 4.271743882e-02 | ΔΩ = -5.60152e-07 | ‖∇‖ = 1.14588e-05
iter  168 | Ω_tilde = 4.271689254e-02 | ΔΩ = -5.46286e-07 | ‖∇‖ = 1.14106e-05
iter  169 | Ω_tilde = 4.271635976e-02 | ΔΩ = -5.32774e-07 | ‖∇‖ = 1.13625e-05
iter  170 | Ω_tilde = 4.271584016e-02 | ΔΩ = -5.19606e-07 | ‖∇‖ = 1.13147e-05
iter  171 | Ω_tilde = 4.271533339e-02 | ΔΩ = -5.06772e-07 | ‖∇‖ = 1.12670e-05
iter  172 | Ω_tilde = 4.271483912e-02 | ΔΩ = -4.94263e-07 | ‖∇‖ = 1.12196e-05
iter  173 | Ω_tilde = 4.271435705e-02 | ΔΩ = -4.82072e-07 | ‖∇‖ = 1.11723e-05
iter  174 | Ω_tilde = 4.271388686e-02 | ΔΩ = -4.70189e-07 | ‖∇‖ = 1.11252e-05
iter  175 | Ω_tilde = 4.271342826e-02 | ΔΩ = -4.58606e-07 | ‖∇‖ = 1.10784e-05
iter  176 | Ω_tilde = 4.271298094e-02 | ΔΩ = -4.47315e-07 | ‖∇‖ = 1.10317e-05
iter  177 | Ω_tilde = 4.271254463e-02 | ΔΩ = -4.36310e-07 | ‖∇‖ = 1.09852e-05
iter  178 | Ω_tilde = 4.271211905e-02 | ΔΩ = -4.25581e-07 | ‖∇‖ = 1.09389e-05
iter  179 | Ω_tilde = 4.271170393e-02 | ΔΩ = -4.15123e-07 | ‖∇‖ = 1.08928e-05
iter  180 | Ω_tilde = 4.271129900e-02 | ΔΩ = -4.04927e-07 | ‖∇‖ = 1.08469e-05
iter  181 | Ω_tilde = 4.271090401e-02 | ΔΩ = -3.94987e-07 | ‖∇‖ = 1.08012e-05
iter  182 | Ω_tilde = 4.271051872e-02 | ΔΩ = -3.85297e-07 | ‖∇‖ = 1.07557e-05
iter  183 | Ω_tilde = 4.271014287e-02 | ΔΩ = -3.75850e-07 | ‖∇‖ = 1.07104e-05
iter  184 | Ω_tilde = 4.270977623e-02 | ΔΩ = -3.66639e-07 | ‖∇‖ = 1.06652e-05
iter  185 | Ω_tilde = 4.270941857e-02 | ΔΩ = -3.57658e-07 | ‖∇‖ = 1.06203e-05
iter  186 | Ω_tilde = 4.270906967e-02 | ΔΩ = -3.48902e-07 | ‖∇‖ = 1.05755e-05
iter  187 | Ω_tilde = 4.270872930e-02 | ΔΩ = -3.40364e-07 | ‖∇‖ = 1.05309e-05
iter  188 | Ω_tilde = 4.270839726e-02 | ΔΩ = -3.32040e-07 | ‖∇‖ = 1.04865e-05
iter  189 | Ω_tilde = 4.270807334e-02 | ΔΩ = -3.23922e-07 | ‖∇‖ = 1.04423e-05
iter  190 | Ω_tilde = 4.270775733e-02 | ΔΩ = -3.16007e-07 | ‖∇‖ = 1.03983e-05
iter  191 | Ω_tilde = 4.270744904e-02 | ΔΩ = -3.08289e-07 | ‖∇‖ = 1.03544e-05
iter  192 | Ω_tilde = 4.270714828e-02 | ΔΩ = -3.00763e-07 | ‖∇‖ = 1.03108e-05
iter  193 | Ω_tilde = 4.270685486e-02 | ΔΩ = -2.93424e-07 | ‖∇‖ = 1.02673e-05
iter  194 | Ω_tilde = 4.270656859e-02 | ΔΩ = -2.86267e-07 | ‖∇‖ = 1.02240e-05
iter  195 | Ω_tilde = 4.270628930e-02 | ΔΩ = -2.79287e-07 | ‖∇‖ = 1.01809e-05
iter  196 | Ω_tilde = 4.270601682e-02 | ΔΩ = -2.72481e-07 | ‖∇‖ = 1.01379e-05
iter  197 | Ω_tilde = 4.270575098e-02 | ΔΩ = -2.65843e-07 | ‖∇‖ = 1.00952e-05
iter  198 | Ω_tilde = 4.270549161e-02 | ΔΩ = -2.59369e-07 | ‖∇‖ = 1.00526e-05
iter  199 | Ω_tilde = 4.270523856e-02 | ΔΩ = -2.53055e-07 | ‖∇‖ = 1.00102e-05
iter  200 | Ω_tilde = 4.270499166e-02 | ΔΩ = -2.46898e-07 | ‖∇‖ = 9.96794e-06
iter  201 | Ω_tilde = 4.270475077e-02 | ΔΩ = -2.40892e-07 | ‖∇‖ = 9.92588e-06
iter  202 | Ω_tilde = 4.270451573e-02 | ΔΩ = -2.35035e-07 | ‖∇‖ = 9.88401e-06
iter  203 | Ω_tilde = 4.270428641e-02 | ΔΩ = -2.29322e-07 | ‖∇‖ = 9.84230e-06
iter  204 | Ω_tilde = 4.270406266e-02 | ΔΩ = -2.23750e-07 | ‖∇‖ = 9.80077e-06
iter  205 | Ω_tilde = 4.270384434e-02 | ΔΩ = -2.18316e-07 | ‖∇‖ = 9.75941e-06
iter  206 | Ω_tilde = 4.270363133e-02 | ΔΩ = -2.13015e-07 | ‖∇‖ = 9.71823e-06
iter  207 | Ω_tilde = 4.270342348e-02 | ΔΩ = -2.07845e-07 | ‖∇‖ = 9.67721e-06
iter  208 | Ω_tilde = 4.270322068e-02 | ΔΩ = -2.02801e-07 | ‖∇‖ = 9.63637e-06
iter  209 | Ω_tilde = 4.270302280e-02 | ΔΩ = -1.97882e-07 | ‖∇‖ = 9.59570e-06
iter  210 | Ω_tilde = 4.270282972e-02 | ΔΩ = -1.93084e-07 | ‖∇‖ = 9.55520e-06
iter  211 | Ω_tilde = 4.270264131e-02 | ΔΩ = -1.88403e-07 | ‖∇‖ = 9.51487e-06
iter  212 | Ω_tilde = 4.270245747e-02 | ΔΩ = -1.83837e-07 | ‖∇‖ = 9.47470e-06
iter  213 | Ω_tilde = 4.270227809e-02 | ΔΩ = -1.79384e-07 | ‖∇‖ = 9.43470e-06
iter  214 | Ω_tilde = 4.270210305e-02 | ΔΩ = -1.75039e-07 | ‖∇‖ = 9.39487e-06
iter  215 | Ω_tilde = 4.270193225e-02 | ΔΩ = -1.70801e-07 | ‖∇‖ = 9.35521e-06
iter  216 | Ω_tilde = 4.270176558e-02 | ΔΩ = -1.66667e-07 | ‖∇‖ = 9.31571e-06
iter  217 | Ω_tilde = 4.270160295e-02 | ΔΩ = -1.62634e-07 | ‖∇‖ = 9.27638e-06
iter  218 | Ω_tilde = 4.270144425e-02 | ΔΩ = -1.58699e-07 | ‖∇‖ = 9.23721e-06
iter  219 | Ω_tilde = 4.270128939e-02 | ΔΩ = -1.54861e-07 | ‖∇‖ = 9.19821e-06
iter  220 | Ω_tilde = 4.270113827e-02 | ΔΩ = -1.51117e-07 | ‖∇‖ = 9.15936e-06
iter  221 | Ω_tilde = 4.270099081e-02 | ΔΩ = -1.47464e-07 | ‖∇‖ = 9.12068e-06
iter  222 | Ω_tilde = 4.270084691e-02 | ΔΩ = -1.43900e-07 | ‖∇‖ = 9.08217e-06
iter  223 | Ω_tilde = 4.270070648e-02 | ΔΩ = -1.40424e-07 | ‖∇‖ = 9.04381e-06
iter  224 | Ω_tilde = 4.270056945e-02 | ΔΩ = -1.37032e-07 | ‖∇‖ = 9.00561e-06
iter  225 | Ω_tilde = 4.270043573e-02 | ΔΩ = -1.33723e-07 | ‖∇‖ = 8.96757e-06
iter  226 | Ω_tilde = 4.270030523e-02 | ΔΩ = -1.30495e-07 | ‖∇‖ = 8.92970e-06
iter  227 | Ω_tilde = 4.270017789e-02 | ΔΩ = -1.27345e-07 | ‖∇‖ = 8.89198e-06
iter  228 | Ω_tilde = 4.270005362e-02 | ΔΩ = -1.24272e-07 | ‖∇‖ = 8.85442e-06
iter  229 | Ω_tilde = 4.269993234e-02 | ΔΩ = -1.21274e-07 | ‖∇‖ = 8.81701e-06
iter  230 | Ω_tilde = 4.269981399e-02 | ΔΩ = -1.18349e-07 | ‖∇‖ = 8.77976e-06
iter  231 | Ω_tilde = 4.269969850e-02 | ΔΩ = -1.15495e-07 | ‖∇‖ = 8.74267e-06
iter  232 | Ω_tilde = 4.269958579e-02 | ΔΩ = -1.12711e-07 | ‖∇‖ = 8.70573e-06
iter  233 | Ω_tilde = 4.269947580e-02 | ΔΩ = -1.09994e-07 | ‖∇‖ = 8.66895e-06
iter  234 | Ω_tilde = 4.269936845e-02 | ΔΩ = -1.07343e-07 | ‖∇‖ = 8.63232e-06
iter  235 | Ω_tilde = 4.269926369e-02 | ΔΩ = -1.04757e-07 | ‖∇‖ = 8.59585e-06
iter  236 | Ω_tilde = 4.269916146e-02 | ΔΩ = -1.02234e-07 | ‖∇‖ = 8.55953e-06
iter  237 | Ω_tilde = 4.269906169e-02 | ΔΩ = -9.97715e-08 | ‖∇‖ = 8.52336e-06
iter  238 | Ω_tilde = 4.269896432e-02 | ΔΩ = -9.73691e-08 | ‖∇‖ = 8.48734e-06
iter  239 | Ω_tilde = 4.269886930e-02 | ΔΩ = -9.50251e-08 | ‖∇‖ = 8.45147e-06
iter  240 | Ω_tilde = 4.269877656e-02 | ΔΩ = -9.27379e-08 | ‖∇‖ = 8.41575e-06
iter  241 | Ω_tilde = 4.269868605e-02 | ΔΩ = -9.05063e-08 | ‖∇‖ = 8.38019e-06
iter  242 | Ω_tilde = 4.269859772e-02 | ΔΩ = -8.83288e-08 | ‖∇‖ = 8.34477e-06
iter  243 | Ω_tilde = 4.269851152e-02 | ΔΩ = -8.62040e-08 | ‖∇‖ = 8.30950e-06
iter  244 | Ω_tilde = 4.269842739e-02 | ΔΩ = -8.41308e-08 | ‖∇‖ = 8.27437e-06
iter  245 | Ω_tilde = 4.269834528e-02 | ΔΩ = -8.21077e-08 | ‖∇‖ = 8.23940e-06
iter  246 | Ω_tilde = 4.269826515e-02 | ΔΩ = -8.01337e-08 | ‖∇‖ = 8.20457e-06
iter  247 | Ω_tilde = 4.269818694e-02 | ΔΩ = -7.82075e-08 | ‖∇‖ = 8.16989e-06
iter  248 | Ω_tilde = 4.269811061e-02 | ΔΩ = -7.63278e-08 | ‖∇‖ = 8.13535e-06
iter  249 | Ω_tilde = 4.269803612e-02 | ΔΩ = -7.44938e-08 | ‖∇‖ = 8.10096e-06
iter  250 | Ω_tilde = 4.269796341e-02 | ΔΩ = -7.27041e-08 | ‖∇‖ = 8.06671e-06
iter  251 | Ω_tilde = 4.269789246e-02 | ΔΩ = -7.09576e-08 | ‖∇‖ = 8.03261e-06
iter  252 | Ω_tilde = 4.269782320e-02 | ΔΩ = -6.92534e-08 | ‖∇‖ = 7.99865e-06
iter  253 | Ω_tilde = 4.269775561e-02 | ΔΩ = -6.75904e-08 | ‖∇‖ = 7.96483e-06
iter  254 | Ω_tilde = 4.269768964e-02 | ΔΩ = -6.59677e-08 | ‖∇‖ = 7.93115e-06
iter  255 | Ω_tilde = 4.269762526e-02 | ΔΩ = -6.43841e-08 | ‖∇‖ = 7.89762e-06
iter  256 | Ω_tilde = 4.269756242e-02 | ΔΩ = -6.28387e-08 | ‖∇‖ = 7.86422e-06
iter  257 | Ω_tilde = 4.269750109e-02 | ΔΩ = -6.13307e-08 | ‖∇‖ = 7.83097e-06
iter  258 | Ω_tilde = 4.269744123e-02 | ΔΩ = -5.98592e-08 | ‖∇‖ = 7.79785e-06
iter  259 | Ω_tilde = 4.269738281e-02 | ΔΩ = -5.84232e-08 | ‖∇‖ = 7.76488e-06
iter  260 | Ω_tilde = 4.269732579e-02 | ΔΩ = -5.70218e-08 | ‖∇‖ = 7.73204e-06
iter  261 | Ω_tilde = 4.269727013e-02 | ΔΩ = -5.56542e-08 | ‖∇‖ = 7.69934e-06
iter  262 | Ω_tilde = 4.269721581e-02 | ΔΩ = -5.43196e-08 | ‖∇‖ = 7.66678e-06
iter  263 | Ω_tilde = 4.269716280e-02 | ΔΩ = -5.30172e-08 | ‖∇‖ = 7.63436e-06
iter  264 | Ω_tilde = 4.269711105e-02 | ΔΩ = -5.17462e-08 | ‖∇‖ = 7.60207e-06
iter  265 | Ω_tilde = 4.269706054e-02 | ΔΩ = -5.05060e-08 | ‖∇‖ = 7.56992e-06
iter  266 | Ω_tilde = 4.269701125e-02 | ΔΩ = -4.92955e-08 | ‖∇‖ = 7.53790e-06
iter  267 | Ω_tilde = 4.269696313e-02 | ΔΩ = -4.81143e-08 | ‖∇‖ = 7.50601e-06
iter  268 | Ω_tilde = 4.269691617e-02 | ΔΩ = -4.69615e-08 | ‖∇‖ = 7.47427e-06
iter  269 | Ω_tilde = 4.269687034e-02 | ΔΩ = -4.58364e-08 | ‖∇‖ = 7.44265e-06
iter  270 | Ω_tilde = 4.269682560e-02 | ΔΩ = -4.47385e-08 | ‖∇‖ = 7.41117e-06
iter  271 | Ω_tilde = 4.269678193e-02 | ΔΩ = -4.36670e-08 | ‖∇‖ = 7.37982e-06
iter  272 | Ω_tilde = 4.269673931e-02 | ΔΩ = -4.26213e-08 | ‖∇‖ = 7.34860e-06
iter  273 | Ω_tilde = 4.269669771e-02 | ΔΩ = -4.16008e-08 | ‖∇‖ = 7.31751e-06
iter  274 | Ω_tilde = 4.269665710e-02 | ΔΩ = -4.06049e-08 | ‖∇‖ = 7.28656e-06
iter  275 | Ω_tilde = 4.269661747e-02 | ΔΩ = -3.96329e-08 | ‖∇‖ = 7.25573e-06
iter  276 | Ω_tilde = 4.269657879e-02 | ΔΩ = -3.86843e-08 | ‖∇‖ = 7.22503e-06
iter  277 | Ω_tilde = 4.269654103e-02 | ΔΩ = -3.77586e-08 | ‖∇‖ = 7.19447e-06
iter  278 | Ω_tilde = 4.269650417e-02 | ΔΩ = -3.68550e-08 | ‖∇‖ = 7.16403e-06
iter  279 | Ω_tilde = 4.269646820e-02 | ΔΩ = -3.59732e-08 | ‖∇‖ = 7.13372e-06
iter  280 | Ω_tilde = 4.269643309e-02 | ΔΩ = -3.51127e-08 | ‖∇‖ = 7.10353e-06
iter  281 | Ω_tilde = 4.269639881e-02 | ΔΩ = -3.42727e-08 | ‖∇‖ = 7.07348e-06
iter  282 | Ω_tilde = 4.269636536e-02 | ΔΩ = -3.34531e-08 | ‖∇‖ = 7.04355e-06
iter  283 | Ω_tilde = 4.269633271e-02 | ΔΩ = -3.26530e-08 | ‖∇‖ = 7.01374e-06
iter  284 | Ω_tilde = 4.269630084e-02 | ΔΩ = -3.18722e-08 | ‖∇‖ = 6.98407e-06
iter  285 | Ω_tilde = 4.269626973e-02 | ΔΩ = -3.11102e-08 | ‖∇‖ = 6.95451e-06
iter  286 | Ω_tilde = 4.269623936e-02 | ΔΩ = -3.03665e-08 | ‖∇‖ = 6.92508e-06
iter  287 | Ω_tilde = 4.269620972e-02 | ΔΩ = -2.96406e-08 | ‖∇‖ = 6.89578e-06
iter  288 | Ω_tilde = 4.269618079e-02 | ΔΩ = -2.89322e-08 | ‖∇‖ = 6.86660e-06
iter  289 | Ω_tilde = 4.269615255e-02 | ΔΩ = -2.82407e-08 | ‖∇‖ = 6.83754e-06
iter  290 | Ω_tilde = 4.269612498e-02 | ΔΩ = -2.75659e-08 | ‖∇‖ = 6.80860e-06
iter  291 | Ω_tilde = 4.269609807e-02 | ΔΩ = -2.69073e-08 | ‖∇‖ = 6.77979e-06
iter  292 | Ω_tilde = 4.269607181e-02 | ΔΩ = -2.62644e-08 | ‖∇‖ = 6.75110e-06
iter  293 | Ω_tilde = 4.269604617e-02 | ΔΩ = -2.56370e-08 | ‖∇‖ = 6.72253e-06
iter  294 | Ω_tilde = 4.269602115e-02 | ΔΩ = -2.50246e-08 | ‖∇‖ = 6.69407e-06
iter  295 | Ω_tilde = 4.269599672e-02 | ΔΩ = -2.44270e-08 | ‖∇‖ = 6.66574e-06
iter  296 | Ω_tilde = 4.269597288e-02 | ΔΩ = -2.38436e-08 | ‖∇‖ = 6.63753e-06
iter  297 | Ω_tilde = 4.269594960e-02 | ΔΩ = -2.32743e-08 | ‖∇‖ = 6.60944e-06
iter  298 | Ω_tilde = 4.269592688e-02 | ΔΩ = -2.27186e-08 | ‖∇‖ = 6.58146e-06
iter  299 | Ω_tilde = 4.269590471e-02 | ΔΩ = -2.21762e-08 | ‖∇‖ = 6.55361e-06
iter  300 | Ω_tilde = 4.269588306e-02 | ΔΩ = -2.16468e-08 | ‖∇‖ = 6.52587e-06
iter  301 | Ω_tilde = 4.269586193e-02 | ΔΩ = -2.11301e-08 | ‖∇‖ = 6.49825e-06
iter  302 | Ω_tilde = 4.269584130e-02 | ΔΩ = -2.06258e-08 | ‖∇‖ = 6.47074e-06
iter  303 | Ω_tilde = 4.269582117e-02 | ΔΩ = -2.01336e-08 | ‖∇‖ = 6.44335e-06
iter  304 | Ω_tilde = 4.269580152e-02 | ΔΩ = -1.96532e-08 | ‖∇‖ = 6.41608e-06
iter  305 | Ω_tilde = 4.269578233e-02 | ΔΩ = -1.91842e-08 | ‖∇‖ = 6.38892e-06
iter  306 | Ω_tilde = 4.269576361e-02 | ΔΩ = -1.87265e-08 | ‖∇‖ = 6.36187e-06
iter  307 | Ω_tilde = 4.269574533e-02 | ΔΩ = -1.82798e-08 | ‖∇‖ = 6.33494e-06
iter  308 | Ω_tilde = 4.269572748e-02 | ΔΩ = -1.78437e-08 | ‖∇‖ = 6.30812e-06
iter  309 | Ω_tilde = 4.269571006e-02 | ΔΩ = -1.74181e-08 | ‖∇‖ = 6.28142e-06
iter  310 | Ω_tilde = 4.269569306e-02 | ΔΩ = -1.70027e-08 | ‖∇‖ = 6.25483e-06
iter  311 | Ω_tilde = 4.269567646e-02 | ΔΩ = -1.65972e-08 | ‖∇‖ = 6.22835e-06
iter  312 | Ω_tilde = 4.269566026e-02 | ΔΩ = -1.62015e-08 | ‖∇‖ = 6.20198e-06
iter  313 | Ω_tilde = 4.269564445e-02 | ΔΩ = -1.58152e-08 | ‖∇‖ = 6.17573e-06
iter  314 | Ω_tilde = 4.269562901e-02 | ΔΩ = -1.54381e-08 | ‖∇‖ = 6.14958e-06
iter  315 | Ω_tilde = 4.269561394e-02 | ΔΩ = -1.50701e-08 | ‖∇‖ = 6.12355e-06
iter  316 | Ω_tilde = 4.269559923e-02 | ΔΩ = -1.47109e-08 | ‖∇‖ = 6.09762e-06
iter  317 | Ω_tilde = 4.269558487e-02 | ΔΩ = -1.43602e-08 | ‖∇‖ = 6.07181e-06
iter  318 | Ω_tilde = 4.269557085e-02 | ΔΩ = -1.40180e-08 | ‖∇‖ = 6.04610e-06
iter  319 | Ω_tilde = 4.269555717e-02 | ΔΩ = -1.36839e-08 | ‖∇‖ = 6.02050e-06
iter  320 | Ω_tilde = 4.269554381e-02 | ΔΩ = -1.33579e-08 | ‖∇‖ = 5.99501e-06
iter  321 | Ω_tilde = 4.269553077e-02 | ΔΩ = -1.30396e-08 | ‖∇‖ = 5.96963e-06
iter  322 | Ω_tilde = 4.269551804e-02 | ΔΩ = -1.27289e-08 | ‖∇‖ = 5.94436e-06
iter  323 | Ω_tilde = 4.269550562e-02 | ΔΩ = -1.24257e-08 | ‖∇‖ = 5.91919e-06
iter  324 | Ω_tilde = 4.269549349e-02 | ΔΩ = -1.21296e-08 | ‖∇‖ = 5.89413e-06
iter  325 | Ω_tilde = 4.269548164e-02 | ΔΩ = -1.18408e-08 | ‖∇‖ = 5.86917e-06
iter  326 | Ω_tilde = 4.269547009e-02 | ΔΩ = -1.15587e-08 | ‖∇‖ = 5.84432e-06
iter  327 | Ω_tilde = 4.269545880e-02 | ΔΩ = -1.12834e-08 | ‖∇‖ = 5.81957e-06
iter  328 | Ω_tilde = 4.269544779e-02 | ΔΩ = -1.10147e-08 | ‖∇‖ = 5.79493e-06
iter  329 | Ω_tilde = 4.269543704e-02 | ΔΩ = -1.07524e-08 | ‖∇‖ = 5.77040e-06
iter  330 | Ω_tilde = 4.269542654e-02 | ΔΩ = -1.04964e-08 | ‖∇‖ = 5.74596e-06
iter  331 | Ω_tilde = 4.269541629e-02 | ΔΩ = -1.02465e-08 | ‖∇‖ = 5.72163e-06
iter  332 | Ω_tilde = 4.269540629e-02 | ΔΩ = -1.00025e-08 | ‖∇‖ = 5.69740e-06
iter  333 | Ω_tilde = 4.269539653e-02 | ΔΩ = -9.76444e-09 | ‖∇‖ = 5.67328e-06
iter  334 | Ω_tilde = 4.269538699e-02 | ΔΩ = -9.53197e-09 | ‖∇‖ = 5.64926e-06
iter  335 | Ω_tilde = 4.269537769e-02 | ΔΩ = -9.30512e-09 | ‖∇‖ = 5.62534e-06
iter  336 | Ω_tilde = 4.269536860e-02 | ΔΩ = -9.08368e-09 | ‖∇‖ = 5.60152e-06
iter  337 | Ω_tilde = 4.269535974e-02 | ΔΩ = -8.86747e-09 | ‖∇‖ = 5.57780e-06
iter  338 | Ω_tilde = 4.269535108e-02 | ΔΩ = -8.65651e-09 | ‖∇‖ = 5.55418e-06
iter  339 | Ω_tilde = 4.269534263e-02 | ΔΩ = -8.45044e-09 | ‖∇‖ = 5.53066e-06
iter  340 | Ω_tilde = 4.269533438e-02 | ΔΩ = -8.24941e-09 | ‖∇‖ = 5.50724e-06
iter  341 | Ω_tilde = 4.269532633e-02 | ΔΩ = -8.05313e-09 | ‖∇‖ = 5.48392e-06
iter  342 | Ω_tilde = 4.269531847e-02 | ΔΩ = -7.86159e-09 | ‖∇‖ = 5.46069e-06
iter  343 | Ω_tilde = 4.269531079e-02 | ΔΩ = -7.67453e-09 | ‖∇‖ = 5.43757e-06
iter  344 | Ω_tilde = 4.269530330e-02 | ΔΩ = -7.49195e-09 | ‖∇‖ = 5.41454e-06
iter  345 | Ω_tilde = 4.269529599e-02 | ΔΩ = -7.31375e-09 | ‖∇‖ = 5.39161e-06
iter  346 | Ω_tilde = 4.269528885e-02 | ΔΩ = -7.13984e-09 | ‖∇‖ = 5.36878e-06
iter  347 | Ω_tilde = 4.269528188e-02 | ΔΩ = -6.97003e-09 | ‖∇‖ = 5.34605e-06
iter  348 | Ω_tilde = 4.269527507e-02 | ΔΩ = -6.80423e-09 | ‖∇‖ = 5.32341e-06
iter  349 | Ω_tilde = 4.269526843e-02 | ΔΩ = -6.64246e-09 | ‖∇‖ = 5.30086e-06
iter  350 | Ω_tilde = 4.269526194e-02 | ΔΩ = -6.48451e-09 | ‖∇‖ = 5.27841e-06
iter  351 | Ω_tilde = 4.269525561e-02 | ΔΩ = -6.33034e-09 | ‖∇‖ = 5.25606e-06
iter  352 | Ω_tilde = 4.269524943e-02 | ΔΩ = -6.17984e-09 | ‖∇‖ = 5.23380e-06
iter  353 | Ω_tilde = 4.269524340e-02 | ΔΩ = -6.03292e-09 | ‖∇‖ = 5.21164e-06
iter  354 | Ω_tilde = 4.269523751e-02 | ΔΩ = -5.88951e-09 | ‖∇‖ = 5.18957e-06
iter  355 | Ω_tilde = 4.269523176e-02 | ΔΩ = -5.74947e-09 | ‖∇‖ = 5.16759e-06
iter  356 | Ω_tilde = 4.269522615e-02 | ΔΩ = -5.61283e-09 | ‖∇‖ = 5.14570e-06
iter  357 | Ω_tilde = 4.269522067e-02 | ΔΩ = -5.47946e-09 | ‖∇‖ = 5.12391e-06
iter  358 | Ω_tilde = 4.269521532e-02 | ΔΩ = -5.34919e-09 | ‖∇‖ = 5.10221e-06
iter  359 | Ω_tilde = 4.269521010e-02 | ΔΩ = -5.22209e-09 | ‖∇‖ = 5.08060e-06
iter  360 | Ω_tilde = 4.269520500e-02 | ΔΩ = -5.09798e-09 | ‖∇‖ = 5.05909e-06
iter  361 | Ω_tilde = 4.269520002e-02 | ΔΩ = -4.97683e-09 | ‖∇‖ = 5.03766e-06
iter  362 | Ω_tilde = 4.269519517e-02 | ΔΩ = -4.85857e-09 | ‖∇‖ = 5.01633e-06
iter  363 | Ω_tilde = 4.269519042e-02 | ΔΩ = -4.74315e-09 | ‖∇‖ = 4.99508e-06
iter  364 | Ω_tilde = 4.269518579e-02 | ΔΩ = -4.63048e-09 | ‖∇‖ = 4.97393e-06
iter  365 | Ω_tilde = 4.269518127e-02 | ΔΩ = -4.52047e-09 | ‖∇‖ = 4.95286e-06
iter  366 | Ω_tilde = 4.269517686e-02 | ΔΩ = -4.41307e-09 | ‖∇‖ = 4.93189e-06
iter  367 | Ω_tilde = 4.269517255e-02 | ΔΩ = -4.30822e-09 | ‖∇‖ = 4.91100e-06
iter  368 | Ω_tilde = 4.269516834e-02 | ΔΩ = -4.20591e-09 | ‖∇‖ = 4.89020e-06
iter  369 | Ω_tilde = 4.269516424e-02 | ΔΩ = -4.10598e-09 | ‖∇‖ = 4.86949e-06
iter  370 | Ω_tilde = 4.269516023e-02 | ΔΩ = -4.00853e-09 | ‖∇‖ = 4.84887e-06
iter  371 | Ω_tilde = 4.269515632e-02 | ΔΩ = -3.91329e-09 | ‖∇‖ = 4.82833e-06
iter  372 | Ω_tilde = 4.269515250e-02 | ΔΩ = -3.82038e-09 | ‖∇‖ = 4.80789e-06
iter  373 | Ω_tilde = 4.269514877e-02 | ΔΩ = -3.72963e-09 | ‖∇‖ = 4.78752e-06
iter  374 | Ω_tilde = 4.269514513e-02 | ΔΩ = -3.64110e-09 | ‖∇‖ = 4.76725e-06
iter  375 | Ω_tilde = 4.269514157e-02 | ΔΩ = -3.55461e-09 | ‖∇‖ = 4.74706e-06
iter  376 | Ω_tilde = 4.269513810e-02 | ΔΩ = -3.47023e-09 | ‖∇‖ = 4.72695e-06
iter  377 | Ω_tilde = 4.269513471e-02 | ΔΩ = -3.38784e-09 | ‖∇‖ = 4.70693e-06
iter  378 | Ω_tilde = 4.269513141e-02 | ΔΩ = -3.30744e-09 | ‖∇‖ = 4.68700e-06
iter  379 | Ω_tilde = 4.269512818e-02 | ΔΩ = -3.22890e-09 | ‖∇‖ = 4.66715e-06
iter  380 | Ω_tilde = 4.269512502e-02 | ΔΩ = -3.15224e-09 | ‖∇‖ = 4.64738e-06
iter  381 | Ω_tilde = 4.269512195e-02 | ΔΩ = -3.07741e-09 | ‖∇‖ = 4.62770e-06
iter  382 | Ω_tilde = 4.269511894e-02 | ΔΩ = -3.00440e-09 | ‖∇‖ = 4.60810e-06
iter  383 | Ω_tilde = 4.269511601e-02 | ΔΩ = -2.93305e-09 | ‖∇‖ = 4.58859e-06
iter  384 | Ω_tilde = 4.269511315e-02 | ΔΩ = -2.86346e-09 | ‖∇‖ = 4.56915e-06
iter  385 | Ω_tilde = 4.269511035e-02 | ΔΩ = -2.79550e-09 | ‖∇‖ = 4.54980e-06
iter  386 | Ω_tilde = 4.269510762e-02 | ΔΩ = -2.72917e-09 | ‖∇‖ = 4.53053e-06
iter  387 | Ω_tilde = 4.269510496e-02 | ΔΩ = -2.66443e-09 | ‖∇‖ = 4.51135e-06
iter  388 | Ω_tilde = 4.269510236e-02 | ΔΩ = -2.60119e-09 | ‖∇‖ = 4.49224e-06
iter  389 | Ω_tilde = 4.269509982e-02 | ΔΩ = -2.53945e-09 | ‖∇‖ = 4.47321e-06
iter  390 | Ω_tilde = 4.269509734e-02 | ΔΩ = -2.47921e-09 | ‖∇‖ = 4.45427e-06
iter  391 | Ω_tilde = 4.269509492e-02 | ΔΩ = -2.42039e-09 | ‖∇‖ = 4.43541e-06
iter  392 | Ω_tilde = 4.269509255e-02 | ΔΩ = -2.36294e-09 | ‖∇‖ = 4.41662e-06
iter  393 | Ω_tilde = 4.269509025e-02 | ΔΩ = -2.30698e-09 | ‖∇‖ = 4.39792e-06
iter  394 | Ω_tilde = 4.269508799e-02 | ΔΩ = -2.25213e-09 | ‖∇‖ = 4.37929e-06
iter  395 | Ω_tilde = 4.269508580e-02 | ΔΩ = -2.19879e-09 | ‖∇‖ = 4.36074e-06
iter  396 | Ω_tilde = 4.269508365e-02 | ΔΩ = -2.14659e-09 | ‖∇‖ = 4.34228e-06
iter  397 | Ω_tilde = 4.269508155e-02 | ΔΩ = -2.09570e-09 | ‖∇‖ = 4.32389e-06
iter  398 | Ω_tilde = 4.269507951e-02 | ΔΩ = -2.04603e-09 | ‖∇‖ = 4.30557e-06
iter  399 | Ω_tilde = 4.269507751e-02 | ΔΩ = -1.99744e-09 | ‖∇‖ = 4.28734e-06
iter  400 | Ω_tilde = 4.269507556e-02 | ΔΩ = -1.95008e-09 | ‖∇‖ = 4.26918e-06
iter  401 | Ω_tilde = 4.269507366e-02 | ΔΩ = -1.90386e-09 | ‖∇‖ = 4.25110e-06
iter  402 | Ω_tilde = 4.269507180e-02 | ΔΩ = -1.85866e-09 | ‖∇‖ = 4.23310e-06
iter  403 | Ω_tilde = 4.269506998e-02 | ΔΩ = -1.81466e-09 | ‖∇‖ = 4.21517e-06
iter  404 | Ω_tilde = 4.269506821e-02 | ΔΩ = -1.77158e-09 | ‖∇‖ = 4.19732e-06
iter  405 | Ω_tilde = 4.269506648e-02 | ΔΩ = -1.72958e-09 | ‖∇‖ = 4.17955e-06
iter  406 | Ω_tilde = 4.269506479e-02 | ΔΩ = -1.68862e-09 | ‖∇‖ = 4.16185e-06
iter  407 | Ω_tilde = 4.269506314e-02 | ΔΩ = -1.64851e-09 | ‖∇‖ = 4.14422e-06
iter  408 | Ω_tilde = 4.269506154e-02 | ΔΩ = -1.60946e-09 | ‖∇‖ = 4.12667e-06
iter  409 | Ω_tilde = 4.269505996e-02 | ΔΩ = -1.57131e-09 | ‖∇‖ = 4.10920e-06
iter  410 | Ω_tilde = 4.269505843e-02 | ΔΩ = -1.53404e-09 | ‖∇‖ = 4.09179e-06
iter  411 | Ω_tilde = 4.269505693e-02 | ΔΩ = -1.49771e-09 | ‖∇‖ = 4.07447e-06
iter  412 | Ω_tilde = 4.269505547e-02 | ΔΩ = -1.46218e-09 | ‖∇‖ = 4.05721e-06
iter  413 | Ω_tilde = 4.269505404e-02 | ΔΩ = -1.42756e-09 | ‖∇‖ = 4.04003e-06
iter  414 | Ω_tilde = 4.269505265e-02 | ΔΩ = -1.39369e-09 | ‖∇‖ = 4.02292e-06
iter  415 | Ω_tilde = 4.269505129e-02 | ΔΩ = -1.36068e-09 | ‖∇‖ = 4.00588e-06
iter  416 | Ω_tilde = 4.269504996e-02 | ΔΩ = -1.32845e-09 | ‖∇‖ = 3.98892e-06
iter  417 | Ω_tilde = 4.269504866e-02 | ΔΩ = -1.29691e-09 | ‖∇‖ = 3.97203e-06
iter  418 | Ω_tilde = 4.269504740e-02 | ΔΩ = -1.26620e-09 | ‖∇‖ = 3.95521e-06
iter  419 | Ω_tilde = 4.269504616e-02 | ΔΩ = -1.23618e-09 | ‖∇‖ = 3.93846e-06
iter  420 | Ω_tilde = 4.269504495e-02 | ΔΩ = -1.20693e-09 | ‖∇‖ = 3.92178e-06
iter  421 | Ω_tilde = 4.269504377e-02 | ΔΩ = -1.17829e-09 | ‖∇‖ = 3.90517e-06
iter  422 | Ω_tilde = 4.269504262e-02 | ΔΩ = -1.15039e-09 | ‖∇‖ = 3.88864e-06
iter  423 | Ω_tilde = 4.269504150e-02 | ΔΩ = -1.12311e-09 | ‖∇‖ = 3.87217e-06
iter  424 | Ω_tilde = 4.269504040e-02 | ΔΩ = -1.09651e-09 | ‖∇‖ = 3.85577e-06
iter  425 | Ω_tilde = 4.269503933e-02 | ΔΩ = -1.07057e-09 | ‖∇‖ = 3.83945e-06
iter  426 | Ω_tilde = 4.269503829e-02 | ΔΩ = -1.04520e-09 | ‖∇‖ = 3.82319e-06
iter  427 | Ω_tilde = 4.269503727e-02 | ΔΩ = -1.02041e-09 | ‖∇‖ = 3.80700e-06
iter  428 | Ω_tilde = 4.269503627e-02 | ΔΩ = -9.96245e-10 | ‖∇‖ = 3.79088e-06
iter  429 | Ω_tilde = 4.269503530e-02 | ΔΩ = -9.72655e-10 | ‖∇‖ = 3.77483e-06
iter  430 | Ω_tilde = 4.269503435e-02 | ΔΩ = -9.49594e-10 | ‖∇‖ = 3.75884e-06
iter  431 | Ω_tilde = 4.269503342e-02 | ΔΩ = -9.27109e-10 | ‖∇‖ = 3.74292e-06
iter  432 | Ω_tilde = 4.269503252e-02 | ΔΩ = -9.05177e-10 | ‖∇‖ = 3.72708e-06
iter  433 | Ω_tilde = 4.269503163e-02 | ΔΩ = -8.83706e-10 | ‖∇‖ = 3.71129e-06
iter  434 | Ω_tilde = 4.269503077e-02 | ΔΩ = -8.62788e-10 | ‖∇‖ = 3.69558e-06
iter  435 | Ω_tilde = 4.269502993e-02 | ΔΩ = -8.42307e-10 | ‖∇‖ = 3.67993e-06
iter  436 | Ω_tilde = 4.269502911e-02 | ΔΩ = -8.22402e-10 | ‖∇‖ = 3.66435e-06
iter  437 | Ω_tilde = 4.269502830e-02 | ΔΩ = -8.02935e-10 | ‖∇‖ = 3.64883e-06
iter  438 | Ω_tilde = 4.269502752e-02 | ΔΩ = -7.83929e-10 | ‖∇‖ = 3.63338e-06
iter  439 | Ω_tilde = 4.269502675e-02 | ΔΩ = -7.65338e-10 | ‖∇‖ = 3.61800e-06
iter  440 | Ω_tilde = 4.269502601e-02 | ΔΩ = -7.47207e-10 | ‖∇‖ = 3.60268e-06
iter  441 | Ω_tilde = 4.269502528e-02 | ΔΩ = -7.29537e-10 | ‖∇‖ = 3.58743e-06
iter  442 | Ω_tilde = 4.269502457e-02 | ΔΩ = -7.12259e-10 | ‖∇‖ = 3.57224e-06
iter  443 | Ω_tilde = 4.269502387e-02 | ΔΩ = -6.95418e-10 | ‖∇‖ = 3.55711e-06
iter  444 | Ω_tilde = 4.269502319e-02 | ΔΩ = -6.78877e-10 | ‖∇‖ = 3.54205e-06
iter  445 | Ω_tilde = 4.269502253e-02 | ΔΩ = -6.62889e-10 | ‖∇‖ = 3.52705e-06
iter  446 | Ω_tilde = 4.269502188e-02 | ΔΩ = -6.47177e-10 | ‖∇‖ = 3.51212e-06
iter  447 | Ω_tilde = 4.269502125e-02 | ΔΩ = -6.31811e-10 | ‖∇‖ = 3.49725e-06
iter  448 | Ω_tilde = 4.269502063e-02 | ΔΩ = -6.16905e-10 | ‖∇‖ = 3.48244e-06
iter  449 | Ω_tilde = 4.269502003e-02 | ΔΩ = -6.02253e-10 | ‖∇‖ = 3.46770e-06
iter  450 | Ω_tilde = 4.269501944e-02 | ΔΩ = -5.88016e-10 | ‖∇‖ = 3.45302e-06
iter  451 | Ω_tilde = 4.269501887e-02 | ΔΩ = -5.74101e-10 | ‖∇‖ = 3.43840e-06
iter  452 | Ω_tilde = 4.269501831e-02 | ΔΩ = -5.60509e-10 | ‖∇‖ = 3.42384e-06
iter  453 | Ω_tilde = 4.269501776e-02 | ΔΩ = -5.47262e-10 | ‖∇‖ = 3.40935e-06
iter  454 | Ω_tilde = 4.269501723e-02 | ΔΩ = -5.34292e-10 | ‖∇‖ = 3.39491e-06
iter  455 | Ω_tilde = 4.269501670e-02 | ΔΩ = -5.21667e-10 | ‖∇‖ = 3.38054e-06
iter  456 | Ω_tilde = 4.269501620e-02 | ΔΩ = -5.09273e-10 | ‖∇‖ = 3.36623e-06
iter  457 | Ω_tilde = 4.269501570e-02 | ΔΩ = -4.97270e-10 | ‖∇‖ = 3.35198e-06
iter  458 | Ω_tilde = 4.269501521e-02 | ΔΩ = -4.85429e-10 | ‖∇‖ = 3.33779e-06
iter  459 | Ω_tilde = 4.269501474e-02 | ΔΩ = -4.74025e-10 | ‖∇‖ = 3.32365e-06
iter  460 | Ω_tilde = 4.269501428e-02 | ΔΩ = -4.62806e-10 | ‖∇‖ = 3.30958e-06
iter  461 | Ω_tilde = 4.269501382e-02 | ΔΩ = -4.51771e-10 | ‖∇‖ = 3.29557e-06
iter  462 | Ω_tilde = 4.269501338e-02 | ΔΩ = -4.41127e-10 | ‖∇‖ = 3.28162e-06
iter  463 | Ω_tilde = 4.269501295e-02 | ΔΩ = -4.30714e-10 | ‖∇‖ = 3.26773e-06
iter  464 | Ω_tilde = 4.269501253e-02 | ΔΩ = -4.20485e-10 | ‖∇‖ = 3.25390e-06
iter  465 | Ω_tilde = 4.269501212e-02 | ΔΩ = -4.10556e-10 | ‖∇‖ = 3.24012e-06
iter  466 | Ω_tilde = 4.269501172e-02 | ΔΩ = -4.00857e-10 | ‖∇‖ = 3.22641e-06
iter  467 | Ω_tilde = 4.269501133e-02 | ΔΩ = -3.91297e-10 | ‖∇‖ = 3.21275e-06
iter  468 | Ω_tilde = 4.269501095e-02 | ΔΩ = -3.82127e-10 | ‖∇‖ = 3.19915e-06
iter  469 | Ω_tilde = 4.269501057e-02 | ΔΩ = -3.73097e-10 | ‖∇‖ = 3.18561e-06
iter  470 | Ω_tilde = 4.269501021e-02 | ΔΩ = -3.64204e-10 | ‖∇‖ = 3.17212e-06
iter  471 | Ω_tilde = 4.269500985e-02 | ΔΩ = -3.55611e-10 | ‖∇‖ = 3.15869e-06
iter  472 | Ω_tilde = 4.269500951e-02 | ΔΩ = -3.47179e-10 | ‖∇‖ = 3.14532e-06
iter  473 | Ω_tilde = 4.269500917e-02 | ΔΩ = -3.38978e-10 | ‖∇‖ = 3.13201e-06
iter  474 | Ω_tilde = 4.269500884e-02 | ΔΩ = -3.30961e-10 | ‖∇‖ = 3.11875e-06
iter  475 | Ω_tilde = 4.269500851e-02 | ΔΩ = -3.23151e-10 | ‖∇‖ = 3.10555e-06
iter  476 | Ω_tilde = 4.269500820e-02 | ΔΩ = -3.15525e-10 | ‖∇‖ = 3.09240e-06
iter  477 | Ω_tilde = 4.269500789e-02 | ΔΩ = -3.08015e-10 | ‖∇‖ = 3.07932e-06
iter  478 | Ω_tilde = 4.269500759e-02 | ΔΩ = -3.00758e-10 | ‖∇‖ = 3.06628e-06
iter  479 | Ω_tilde = 4.269500730e-02 | ΔΩ = -2.93593e-10 | ‖∇‖ = 3.05330e-06
iter  480 | Ω_tilde = 4.269500701e-02 | ΔΩ = -2.86705e-10 | ‖∇‖ = 3.04038e-06
iter  481 | Ω_tilde = 4.269500673e-02 | ΔΩ = -2.79886e-10 | ‖∇‖ = 3.02751e-06
iter  482 | Ω_tilde = 4.269500646e-02 | ΔΩ = -2.73320e-10 | ‖∇‖ = 3.01470e-06
iter  483 | Ω_tilde = 4.269500619e-02 | ΔΩ = -2.66823e-10 | ‖∇‖ = 3.00194e-06
iter  484 | Ω_tilde = 4.269500593e-02 | ΔΩ = -2.60511e-10 | ‖∇‖ = 2.98923e-06
iter  485 | Ω_tilde = 4.269500567e-02 | ΔΩ = -2.54337e-10 | ‖∇‖ = 2.97658e-06
iter  486 | Ω_tilde = 4.269500543e-02 | ΔΩ = -2.48347e-10 | ‖∇‖ = 2.96398e-06
iter  487 | Ω_tilde = 4.269500518e-02 | ΔΩ = -2.42449e-10 | ‖∇‖ = 2.95144e-06
iter  488 | Ω_tilde = 4.269500495e-02 | ΔΩ = -2.36782e-10 | ‖∇‖ = 2.93895e-06
iter  489 | Ω_tilde = 4.269500472e-02 | ΔΩ = -2.31115e-10 | ‖∇‖ = 2.92651e-06
iter  490 | Ω_tilde = 4.269500449e-02 | ΔΩ = -2.25655e-10 | ‖∇‖ = 2.91412e-06
iter  491 | Ω_tilde = 4.269500427e-02 | ΔΩ = -2.20333e-10 | ‖∇‖ = 2.90179e-06
iter  492 | Ω_tilde = 4.269500405e-02 | ΔΩ = -2.15127e-10 | ‖∇‖ = 2.88951e-06
iter  493 | Ω_tilde = 4.269500384e-02 | ΔΩ = -2.09989e-10 | ‖∇‖ = 2.87728e-06
iter  494 | Ω_tilde = 4.269500364e-02 | ΔΩ = -2.05128e-10 | ‖∇‖ = 2.86510e-06
iter  495 | Ω_tilde = 4.269500344e-02 | ΔΩ = -2.00221e-10 | ‖∇‖ = 2.85298e-06
iter  496 | Ω_tilde = 4.269500324e-02 | ΔΩ = -1.95452e-10 | ‖∇‖ = 2.84091e-06
iter  497 | Ω_tilde = 4.269500305e-02 | ΔΩ = -1.90868e-10 | ‖∇‖ = 2.82888e-06
iter  498 | Ω_tilde = 4.269500287e-02 | ΔΩ = -1.86329e-10 | ‖∇‖ = 2.81691e-06
iter  499 | Ω_tilde = 4.269500268e-02 | ΔΩ = -1.82021e-10 | ‖∇‖ = 2.80499e-06
iter  500 | Ω_tilde = 4.269500251e-02 | ΔΩ = -1.77598e-10 | ‖∇‖ = 2.79312e-06
iter  501 | Ω_tilde = 4.269500233e-02 | ΔΩ = -1.73451e-10 | ‖∇‖ = 2.78130e-06
iter  502 | Ω_tilde = 4.269500216e-02 | ΔΩ = -1.69351e-10 | ‖∇‖ = 2.76953e-06
iter  503 | Ω_tilde = 4.269500200e-02 | ΔΩ = -1.65319e-10 | ‖∇‖ = 2.75781e-06
iter  504 | Ω_tilde = 4.269500184e-02 | ΔΩ = -1.61472e-10 | ‖∇‖ = 2.74614e-06
iter  505 | Ω_tilde = 4.269500168e-02 | ΔΩ = -1.57601e-10 | ‖∇‖ = 2.73452e-06
iter  506 | Ω_tilde = 4.269500153e-02 | ΔΩ = -1.53915e-10 | ‖∇‖ = 2.72295e-06
iter  507 | Ω_tilde = 4.269500138e-02 | ΔΩ = -1.50229e-10 | ‖∇‖ = 2.71143e-06
iter  508 | Ω_tilde = 4.269500123e-02 | ΔΩ = -1.46751e-10 | ‖∇‖ = 2.69996e-06
iter  509 | Ω_tilde = 4.269500109e-02 | ΔΩ = -1.43226e-10 | ‖∇‖ = 2.68853e-06
iter  510 | Ω_tilde = 4.269500095e-02 | ΔΩ = -1.39839e-10 | ‖∇‖ = 2.67716e-06
iter  511 | Ω_tilde = 4.269500081e-02 | ΔΩ = -1.36545e-10 | ‖∇‖ = 2.66583e-06
iter  512 | Ω_tilde = 4.269500068e-02 | ΔΩ = -1.33320e-10 | ‖∇‖ = 2.65455e-06
iter  513 | Ω_tilde = 4.269500055e-02 | ΔΩ = -1.30163e-10 | ‖∇‖ = 2.64332e-06
iter  514 | Ω_tilde = 4.269500042e-02 | ΔΩ = -1.27076e-10 | ‖∇‖ = 2.63214e-06
iter  515 | Ω_tilde = 4.269500029e-02 | ΔΩ = -1.24104e-10 | ‖∇‖ = 2.62100e-06
iter  516 | Ω_tilde = 4.269500017e-02 | ΔΩ = -1.21133e-10 | ‖∇‖ = 2.60991e-06
iter  517 | Ω_tilde = 4.269500005e-02 | ΔΩ = -1.18322e-10 | ‖∇‖ = 2.59887e-06
iter  518 | Ω_tilde = 4.269499994e-02 | ΔΩ = -1.15511e-10 | ‖∇‖ = 2.58787e-06
iter  519 | Ω_tilde = 4.269499983e-02 | ΔΩ = -1.12770e-10 | ‖∇‖ = 2.57693e-06
iter  520 | Ω_tilde = 4.269499972e-02 | ΔΩ = -1.10098e-10 | ‖∇‖ = 2.56602e-06
iter  521 | Ω_tilde = 4.269499961e-02 | ΔΩ = -1.07494e-10 | ‖∇‖ = 2.55517e-06
iter  522 | Ω_tilde = 4.269499950e-02 | ΔΩ = -1.04983e-10 | ‖∇‖ = 2.54436e-06
iter  523 | Ω_tilde = 4.269499940e-02 | ΔΩ = -1.02472e-10 | ‖∇‖ = 2.53359e-06
iter  524 | Ω_tilde = 4.269499930e-02 | ΔΩ = -1.00053e-10 | ‖∇‖ = 2.52288e-06
iter  525 | Ω_tilde = 4.269499920e-02 | ΔΩ = -9.76801e-11 | ‖∇‖ = 2.51220e-06
iter  526 | Ω_tilde = 4.269499911e-02 | ΔΩ = -9.53994e-11 | ‖∇‖ = 2.50158e-06
iter  527 | Ω_tilde = 4.269499902e-02 | ΔΩ = -9.31417e-11 | ‖∇‖ = 2.49099e-06
iter  528 | Ω_tilde = 4.269499892e-02 | ΔΩ = -9.09070e-11 | ‖∇‖ = 2.48046e-06
iter  529 | Ω_tilde = 4.269499884e-02 | ΔΩ = -8.88106e-11 | ‖∇‖ = 2.46996e-06
iter  530 | Ω_tilde = 4.269499875e-02 | ΔΩ = -8.67142e-11 | ‖∇‖ = 2.45951e-06
iter  531 | Ω_tilde = 4.269499866e-02 | ΔΩ = -8.46177e-11 | ‖∇‖ = 2.44911e-06
iter  532 | Ω_tilde = 4.269499858e-02 | ΔΩ = -8.26365e-11 | ‖∇‖ = 2.43875e-06
iter  533 | Ω_tilde = 4.269499850e-02 | ΔΩ = -8.07013e-11 | ‖∇‖ = 2.42844e-06
iter  534 | Ω_tilde = 4.269499842e-02 | ΔΩ = -7.87661e-11 | ‖∇‖ = 2.41816e-06
iter  535 | Ω_tilde = 4.269499835e-02 | ΔΩ = -7.69231e-11 | ‖∇‖ = 2.40793e-06
iter  536 | Ω_tilde = 4.269499827e-02 | ΔΩ = -7.50801e-11 | ‖∇‖ = 2.39775e-06
iter  537 | Ω_tilde = 4.269499820e-02 | ΔΩ = -7.33523e-11 | ‖∇‖ = 2.38761e-06
iter  538 | Ω_tilde = 4.269499813e-02 | ΔΩ = -7.16014e-11 | ‖∇‖ = 2.37751e-06
iter  539 | Ω_tilde = 4.269499806e-02 | ΔΩ = -6.98966e-11 | ‖∇‖ = 2.36745e-06
iter  540 | Ω_tilde = 4.269499799e-02 | ΔΩ = -6.82839e-11 | ‖∇‖ = 2.35744e-06
iter  541 | Ω_tilde = 4.269499792e-02 | ΔΩ = -6.66252e-11 | ‖∇‖ = 2.34747e-06
iter  542 | Ω_tilde = 4.269499786e-02 | ΔΩ = -6.50587e-11 | ‖∇‖ = 2.33754e-06
iter  543 | Ω_tilde = 4.269499779e-02 | ΔΩ = -6.35382e-11 | ‖∇‖ = 2.32765e-06
iter  544 | Ω_tilde = 4.269499773e-02 | ΔΩ = -6.20177e-11 | ‖∇‖ = 2.31781e-06
iter  545 | Ω_tilde = 4.269499767e-02 | ΔΩ = -6.05663e-11 | ‖∇‖ = 2.30801e-06
iter  546 | Ω_tilde = 4.269499761e-02 | ΔΩ = -5.91380e-11 | ‖∇‖ = 2.29825e-06
iter  547 | Ω_tilde = 4.269499755e-02 | ΔΩ = -5.77557e-11 | ‖∇‖ = 2.28853e-06
iter  548 | Ω_tilde = 4.269499750e-02 | ΔΩ = -5.63734e-11 | ‖∇‖ = 2.27885e-06
iter  549 | Ω_tilde = 4.269499744e-02 | ΔΩ = -5.50142e-11 | ‖∇‖ = 2.26921e-06
iter  550 | Ω_tilde = 4.269499739e-02 | ΔΩ = -5.37702e-11 | ‖∇‖ = 2.25961e-06
iter  551 | Ω_tilde = 4.269499733e-02 | ΔΩ = -5.24570e-11 | ‖∇‖ = 2.25006e-06
iter  552 | Ω_tilde = 4.269499728e-02 | ΔΩ = -5.12360e-11 | ‖∇‖ = 2.24054e-06
iter  553 | Ω_tilde = 4.269499723e-02 | ΔΩ = -5.00380e-11 | ‖∇‖ = 2.23107e-06
iter  554 | Ω_tilde = 4.269499718e-02 | ΔΩ = -4.88401e-11 | ‖∇‖ = 2.22163e-06
iter  555 | Ω_tilde = 4.269499714e-02 | ΔΩ = -4.76651e-11 | ‖∇‖ = 2.21224e-06
iter  556 | Ω_tilde = 4.269499709e-02 | ΔΩ = -4.65824e-11 | ‖∇‖ = 2.20289e-06
iter  557 | Ω_tilde = 4.269499704e-02 | ΔΩ = -4.54766e-11 | ‖∇‖ = 2.19357e-06
iter  558 | Ω_tilde = 4.269499700e-02 | ΔΩ = -4.43938e-11 | ‖∇‖ = 2.18429e-06
iter  559 | Ω_tilde = 4.269499696e-02 | ΔΩ = -4.33340e-11 | ‖∇‖ = 2.17506e-06
iter  560 | Ω_tilde = 4.269499691e-02 | ΔΩ = -4.23434e-11 | ‖∇‖ = 2.16586e-06
iter  561 | Ω_tilde = 4.269499687e-02 | ΔΩ = -4.13298e-11 | ‖∇‖ = 2.15670e-06
iter  562 | Ω_tilde = 4.269499683e-02 | ΔΩ = -4.03161e-11 | ‖∇‖ = 2.14758e-06
iter  563 | Ω_tilde = 4.269499679e-02 | ΔΩ = -3.93946e-11 | ‖∇‖ = 2.13850e-06
iter  564 | Ω_tilde = 4.269499676e-02 | ΔΩ = -3.84731e-11 | ‖∇‖ = 2.12946e-06
iter  565 | Ω_tilde = 4.269499672e-02 | ΔΩ = -3.75285e-11 | ‖∇‖ = 2.12046e-06
iter  566 | Ω_tilde = 4.269499668e-02 | ΔΩ = -3.66992e-11 | ‖∇‖ = 2.11149e-06
iter  567 | Ω_tilde = 4.269499665e-02 | ΔΩ = -3.58007e-11 | ‖∇‖ = 2.10257e-06
iter  568 | Ω_tilde = 4.269499661e-02 | ΔΩ = -3.49483e-11 | ‖∇‖ = 2.09368e-06
iter  569 | Ω_tilde = 4.269499658e-02 | ΔΩ = -3.41650e-11 | ‖∇‖ = 2.08483e-06
iter  570 | Ω_tilde = 4.269499654e-02 | ΔΩ = -3.33357e-11 | ‖∇‖ = 2.07601e-06
iter  571 | Ω_tilde = 4.269499651e-02 | ΔΩ = -3.24833e-11 | ‖∇‖ = 2.06723e-06
iter  572 | Ω_tilde = 4.269499648e-02 | ΔΩ = -3.18152e-11 | ‖∇‖ = 2.05850e-06
iter  573 | Ω_tilde = 4.269499645e-02 | ΔΩ = -3.10319e-11 | ‖∇‖ = 2.04979e-06
iter  574 | Ω_tilde = 4.269499642e-02 | ΔΩ = -3.02716e-11 | ‖∇‖ = 2.04113e-06
iter  575 | Ω_tilde = 4.269499639e-02 | ΔΩ = -2.95575e-11 | ‖∇‖ = 2.03250e-06
iter  576 | Ω_tilde = 4.269499636e-02 | ΔΩ = -2.88663e-11 | ‖∇‖ = 2.02391e-06
iter  577 | Ω_tilde = 4.269499633e-02 | ΔΩ = -2.81982e-11 | ‖∇‖ = 2.01535e-06
iter  578 | Ω_tilde = 4.269499630e-02 | ΔΩ = -2.75301e-11 | ‖∇‖ = 2.00683e-06
iter  579 | Ω_tilde = 4.269499628e-02 | ΔΩ = -2.68851e-11 | ‖∇‖ = 1.99835e-06
iter  580 | Ω_tilde = 4.269499625e-02 | ΔΩ = -2.62170e-11 | ‖∇‖ = 1.98990e-06
iter  581 | Ω_tilde = 4.269499622e-02 | ΔΩ = -2.56641e-11 | ‖∇‖ = 1.98149e-06
iter  582 | Ω_tilde = 4.269499620e-02 | ΔΩ = -2.49729e-11 | ‖∇‖ = 1.97311e-06
iter  583 | Ω_tilde = 4.269499617e-02 | ΔΩ = -2.44661e-11 | ‖∇‖ = 1.96477e-06
iter  584 | Ω_tilde = 4.269499615e-02 | ΔΩ = -2.38211e-11 | ‖∇‖ = 1.95647e-06
iter  585 | Ω_tilde = 4.269499613e-02 | ΔΩ = -2.32912e-11 | ‖∇‖ = 1.94820e-06
iter  586 | Ω_tilde = 4.269499610e-02 | ΔΩ = -2.27613e-11 | ‖∇‖ = 1.93996e-06
iter  587 | Ω_tilde = 4.269499608e-02 | ΔΩ = -2.22314e-11 | ‖∇‖ = 1.93177e-06
iter  588 | Ω_tilde = 4.269499606e-02 | ΔΩ = -2.16325e-11 | ‖∇‖ = 1.92360e-06
iter  589 | Ω_tilde = 4.269499604e-02 | ΔΩ = -2.11487e-11 | ‖∇‖ = 1.91547e-06
iter  590 | Ω_tilde = 4.269499602e-02 | ΔΩ = -2.07340e-11 | ‖∇‖ = 1.90737e-06
iter  591 | Ω_tilde = 4.269499600e-02 | ΔΩ = -2.01580e-11 | ‖∇‖ = 1.89931e-06
iter  592 | Ω_tilde = 4.269499598e-02 | ΔΩ = -1.96743e-11 | ‖∇‖ = 1.89129e-06
iter  593 | Ω_tilde = 4.269499596e-02 | ΔΩ = -1.92826e-11 | ‖∇‖ = 1.88329e-06
iter  594 | Ω_tilde = 4.269499594e-02 | ΔΩ = -1.87527e-11 | ‖∇‖ = 1.87533e-06
iter  595 | Ω_tilde = 4.269499592e-02 | ΔΩ = -1.83381e-11 | ‖∇‖ = 1.86741e-06
iter  596 | Ω_tilde = 4.269499590e-02 | ΔΩ = -1.79464e-11 | ‖∇‖ = 1.85952e-06
iter  597 | Ω_tilde = 4.269499589e-02 | ΔΩ = -1.74626e-11 | ‖∇‖ = 1.85166e-06
iter  598 | Ω_tilde = 4.269499587e-02 | ΔΩ = -1.70479e-11 | ‖∇‖ = 1.84383e-06
iter  599 | Ω_tilde = 4.269499585e-02 | ΔΩ = -1.66793e-11 | ‖∇‖ = 1.83604e-06
iter  600 | Ω_tilde = 4.269499584e-02 | ΔΩ = -1.62647e-11 | ‖∇‖ = 1.82828e-06
iter  601 | Ω_tilde = 4.269499582e-02 | ΔΩ = -1.58730e-11 | ‖∇‖ = 1.82056e-06
iter  602 | Ω_tilde = 4.269499581e-02 | ΔΩ = -1.55735e-11 | ‖∇‖ = 1.81286e-06
iter  603 | Ω_tilde = 4.269499579e-02 | ΔΩ = -1.51589e-11 | ‖∇‖ = 1.80520e-06
iter  604 | Ω_tilde = 4.269499578e-02 | ΔΩ = -1.47672e-11 | ‖∇‖ = 1.79758e-06
iter  605 | Ω_tilde = 4.269499576e-02 | ΔΩ = -1.44908e-11 | ‖∇‖ = 1.78998e-06
iter  606 | Ω_tilde = 4.269499575e-02 | ΔΩ = -1.40530e-11 | ‖∇‖ = 1.78242e-06
iter  607 | Ω_tilde = 4.269499573e-02 | ΔΩ = -1.37996e-11 | ‖∇‖ = 1.77488e-06
iter  608 | Ω_tilde = 4.269499572e-02 | ΔΩ = -1.34310e-11 | ‖∇‖ = 1.76739e-06
iter  609 | Ω_tilde = 4.269499571e-02 | ΔΩ = -1.31315e-11 | ‖∇‖ = 1.75992e-06
iter  610 | Ω_tilde = 4.269499569e-02 | ΔΩ = -1.28090e-11 | ‖∇‖ = 1.75248e-06
iter  611 | Ω_tilde = 4.269499568e-02 | ΔΩ = -1.25326e-11 | ‖∇‖ = 1.74508e-06
iter  612 | Ω_tilde = 4.269499567e-02 | ΔΩ = -1.22100e-11 | ‖∇‖ = 1.73770e-06
iter  613 | Ω_tilde = 4.269499566e-02 | ΔΩ = -1.19566e-11 | ‖∇‖ = 1.73036e-06
iter  614 | Ω_tilde = 4.269499565e-02 | ΔΩ = -1.16110e-11 | ‖∇‖ = 1.72305e-06
iter  615 | Ω_tilde = 4.269499563e-02 | ΔΩ = -1.14037e-11 | ‖∇‖ = 1.71577e-06
iter  616 | Ω_tilde = 4.269499562e-02 | ΔΩ = -1.11042e-11 | ‖∇‖ = 1.70852e-06
iter  617 | Ω_tilde = 4.269499561e-02 | ΔΩ = -1.08508e-11 | ‖∇‖ = 1.70131e-06
iter  618 | Ω_tilde = 4.269499560e-02 | ΔΩ = -1.05974e-11 | ‖∇‖ = 1.69412e-06
iter  619 | Ω_tilde = 4.269499559e-02 | ΔΩ = -1.03440e-11 | ‖∇‖ = 1.68696e-06
iter  620 | Ω_tilde = 4.269499558e-02 | ΔΩ = -1.01136e-11 | ‖∇‖ = 1.67983e-06
iter  621 | Ω_tilde = 4.269499557e-02 | ΔΩ = -9.81409e-12 | ‖∇‖ = 1.67274e-06
iter  622 | Ω_tilde = 4.269499556e-02 | ΔΩ = -9.65282e-12 | ‖∇‖ = 1.66567e-06
iter  623 | Ω_tilde = 4.269499555e-02 | ΔΩ = -9.37637e-12 | ‖∇‖ = 1.65864e-06
iter  624 | Ω_tilde = 4.269499554e-02 | ΔΩ = -9.19207e-12 | ‖∇‖ = 1.65163e-06
iter  625 | Ω_tilde = 4.269499553e-02 | ΔΩ = -8.93866e-12 | ‖∇‖ = 1.64465e-06
iter  626 | Ω_tilde = 4.269499553e-02 | ΔΩ = -8.77739e-12 | ‖∇‖ = 1.63771e-06
iter  627 | Ω_tilde = 4.269499552e-02 | ΔΩ = -8.52397e-12 | ‖∇‖ = 1.63079e-06
iter  628 | Ω_tilde = 4.269499551e-02 | ΔΩ = -8.33968e-12 | ‖∇‖ = 1.62390e-06
iter  629 | Ω_tilde = 4.269499550e-02 | ΔΩ = -8.15537e-12 | ‖∇‖ = 1.61704e-06
iter  630 | Ω_tilde = 4.269499549e-02 | ΔΩ = -7.94803e-12 | ‖∇‖ = 1.61021e-06
iter  631 | Ω_tilde = 4.269499549e-02 | ΔΩ = -7.76373e-12 | ‖∇‖ = 1.60341e-06
iter  632 | Ω_tilde = 4.269499548e-02 | ΔΩ = -7.60246e-12 | ‖∇‖ = 1.59664e-06
iter  633 | Ω_tilde = 4.269499547e-02 | ΔΩ = -7.39513e-12 | ‖∇‖ = 1.58989e-06
iter  634 | Ω_tilde = 4.269499546e-02 | ΔΩ = -7.23386e-12 | ‖∇‖ = 1.58318e-06
iter  635 | Ω_tilde = 4.269499546e-02 | ΔΩ = -7.04956e-12 | ‖∇‖ = 1.57649e-06
iter  636 | Ω_tilde = 4.269499545e-02 | ΔΩ = -6.88830e-12 | ‖∇‖ = 1.56983e-06
iter  637 | Ω_tilde = 4.269499544e-02 | ΔΩ = -6.70399e-12 | ‖∇‖ = 1.56320e-06
iter  638 | Ω_tilde = 4.269499544e-02 | ΔΩ = -6.61184e-12 | ‖∇‖ = 1.55660e-06
iter  639 | Ω_tilde = 4.269499543e-02 | ΔΩ = -6.38146e-12 | ‖∇‖ = 1.55003e-06
iter  640 | Ω_tilde = 4.269499542e-02 | ΔΩ = -6.24324e-12 | ‖∇‖ = 1.54348e-06
iter  641 | Ω_tilde = 4.269499542e-02 | ΔΩ = -6.12805e-12 | ‖∇‖ = 1.53696e-06
iter  642 | Ω_tilde = 4.269499541e-02 | ΔΩ = -5.98982e-12 | ‖∇‖ = 1.53047e-06
iter  643 | Ω_tilde = 4.269499541e-02 | ΔΩ = -5.82855e-12 | ‖∇‖ = 1.52401e-06
iter  644 | Ω_tilde = 4.269499540e-02 | ΔΩ = -5.71337e-12 | ‖∇‖ = 1.51757e-06
iter  645 | Ω_tilde = 4.269499539e-02 | ΔΩ = -5.52907e-12 | ‖∇‖ = 1.51117e-06
iter  646 | Ω_tilde = 4.269499539e-02 | ΔΩ = -5.45994e-12 | ‖∇‖ = 1.50479e-06
iter  647 | Ω_tilde = 4.269499538e-02 | ΔΩ = -5.29869e-12 | ‖∇‖ = 1.49843e-06
iter  648 | Ω_tilde = 4.269499538e-02 | ΔΩ = -5.13742e-12 | ‖∇‖ = 1.49210e-06
iter  649 | Ω_tilde = 4.269499537e-02 | ΔΩ = -5.09135e-12 | ‖∇‖ = 1.48580e-06
iter  650 | Ω_tilde = 4.269499537e-02 | ΔΩ = -4.93008e-12 | ‖∇‖ = 1.47953e-06
iter  651 | Ω_tilde = 4.269499536e-02 | ΔΩ = -4.83793e-12 | ‖∇‖ = 1.47328e-06
iter  652 | Ω_tilde = 4.269499536e-02 | ΔΩ = -4.69971e-12 | ‖∇‖ = 1.46706e-06
iter  653 | Ω_tilde = 4.269499535e-02 | ΔΩ = -4.58451e-12 | ‖∇‖ = 1.46087e-06
iter  654 | Ω_tilde = 4.269499535e-02 | ΔΩ = -4.46933e-12 | ‖∇‖ = 1.45470e-06
iter  655 | Ω_tilde = 4.269499534e-02 | ΔΩ = -4.40022e-12 | ‖∇‖ = 1.44856e-06
iter  656 | Ω_tilde = 4.269499534e-02 | ΔΩ = -4.26199e-12 | ‖∇‖ = 1.44244e-06
iter  657 | Ω_tilde = 4.269499534e-02 | ΔΩ = -4.16984e-12 | ‖∇‖ = 1.43635e-06
iter  658 | Ω_tilde = 4.269499533e-02 | ΔΩ = -4.07768e-12 | ‖∇‖ = 1.43029e-06
iter  659 | Ω_tilde = 4.269499533e-02 | ΔΩ = -3.96250e-12 | ‖∇‖ = 1.42425e-06
iter  660 | Ω_tilde = 4.269499532e-02 | ΔΩ = -3.89339e-12 | ‖∇‖ = 1.41824e-06
iter  661 | Ω_tilde = 4.269499532e-02 | ΔΩ = -3.80123e-12 | ‖∇‖ = 1.41225e-06
iter  662 | Ω_tilde = 4.269499532e-02 | ΔΩ = -3.70908e-12 | ‖∇‖ = 1.40629e-06
iter  663 | Ω_tilde = 4.269499531e-02 | ΔΩ = -3.61693e-12 | ‖∇‖ = 1.40035e-06
iter  664 | Ω_tilde = 4.269499531e-02 | ΔΩ = -3.52478e-12 | ‖∇‖ = 1.39444e-06
iter  665 | Ω_tilde = 4.269499531e-02 | ΔΩ = -3.45567e-12 | ‖∇‖ = 1.38855e-06
iter  666 | Ω_tilde = 4.269499530e-02 | ΔΩ = -3.38655e-12 | ‖∇‖ = 1.38269e-06
iter  667 | Ω_tilde = 4.269499530e-02 | ΔΩ = -3.27136e-12 | ‖∇‖ = 1.37685e-06
iter  668 | Ω_tilde = 4.269499530e-02 | ΔΩ = -3.17921e-12 | ‖∇‖ = 1.37104e-06
iter  669 | Ω_tilde = 4.269499529e-02 | ΔΩ = -3.15618e-12 | ‖∇‖ = 1.36526e-06
iter  670 | Ω_tilde = 4.269499529e-02 | ΔΩ = -3.08706e-12 | ‖∇‖ = 1.35949e-06
iter  671 | Ω_tilde = 4.269499529e-02 | ΔΩ = -2.99491e-12 | ‖∇‖ = 1.35375e-06
iter  672 | Ω_tilde = 4.269499528e-02 | ΔΩ = -2.90276e-12 | ‖∇‖ = 1.34804e-06
iter  673 | Ω_tilde = 4.269499528e-02 | ΔΩ = -2.85668e-12 | ‖∇‖ = 1.34235e-06
iter  674 | Ω_tilde = 4.269499528e-02 | ΔΩ = -2.78758e-12 | ‖∇‖ = 1.33668e-06
iter  675 | Ω_tilde = 4.269499528e-02 | ΔΩ = -2.69541e-12 | ‖∇‖ = 1.33104e-06
iter  676 | Ω_tilde = 4.269499527e-02 | ΔΩ = -2.64935e-12 | ‖∇‖ = 1.32543e-06
iter  677 | Ω_tilde = 4.269499527e-02 | ΔΩ = -2.62631e-12 | ‖∇‖ = 1.31983e-06
iter  678 | Ω_tilde = 4.269499527e-02 | ΔΩ = -2.51112e-12 | ‖∇‖ = 1.31426e-06
iter  679 | Ω_tilde = 4.269499527e-02 | ΔΩ = -2.46504e-12 | ‖∇‖ = 1.30871e-06
iter  680 | Ω_tilde = 4.269499526e-02 | ΔΩ = -2.39592e-12 | ‖∇‖ = 1.30319e-06
iter  681 | Ω_tilde = 4.269499526e-02 | ΔΩ = -2.34986e-12 | ‖∇‖ = 1.29769e-06
iter  682 | Ω_tilde = 4.269499526e-02 | ΔΩ = -2.30378e-12 | ‖∇‖ = 1.29222e-06
iter  683 | Ω_tilde = 4.269499526e-02 | ΔΩ = -2.25770e-12 | ‖∇‖ = 1.28676e-06
iter  684 | Ω_tilde = 4.269499525e-02 | ΔΩ = -2.18859e-12 | ‖∇‖ = 1.28133e-06
iter  685 | Ω_tilde = 4.269499525e-02 | ΔΩ = -2.14251e-12 | ‖∇‖ = 1.27592e-06
iter  686 | Ω_tilde = 4.269499525e-02 | ΔΩ = -2.09644e-12 | ‖∇‖ = 1.27054e-06
iter  687 | Ω_tilde = 4.269499525e-02 | ΔΩ = -2.05036e-12 | ‖∇‖ = 1.26518e-06
iter  688 | Ω_tilde = 4.269499525e-02 | ΔΩ = -1.95821e-12 | ‖∇‖ = 1.25984e-06
iter  689 | Ω_tilde = 4.269499524e-02 | ΔΩ = -1.98125e-12 | ‖∇‖ = 1.25452e-06
iter  690 | Ω_tilde = 4.269499524e-02 | ΔΩ = -1.86606e-12 | ‖∇‖ = 1.24923e-06
iter  691 | Ω_tilde = 4.269499524e-02 | ΔΩ = -1.86606e-12 | ‖∇‖ = 1.24396e-06
iter  692 | Ω_tilde = 4.269499524e-02 | ΔΩ = -1.84302e-12 | ‖∇‖ = 1.23871e-06
iter  693 | Ω_tilde = 4.269499524e-02 | ΔΩ = -1.75087e-12 | ‖∇‖ = 1.23348e-06
iter  694 | Ω_tilde = 4.269499524e-02 | ΔΩ = -1.72783e-12 | ‖∇‖ = 1.22828e-06
iter  695 | Ω_tilde = 4.269499523e-02 | ΔΩ = -1.70480e-12 | ‖∇‖ = 1.22310e-06
iter  696 | Ω_tilde = 4.269499523e-02 | ΔΩ = -1.63568e-12 | ‖∇‖ = 1.21794e-06
iter  697 | Ω_tilde = 4.269499523e-02 | ΔΩ = -1.58960e-12 | ‖∇‖ = 1.21280e-06
iter  698 | Ω_tilde = 4.269499523e-02 | ΔΩ = -1.56657e-12 | ‖∇‖ = 1.20768e-06
iter  699 | Ω_tilde = 4.269499523e-02 | ΔΩ = -1.54353e-12 | ‖∇‖ = 1.20259e-06
iter  700 | Ω_tilde = 4.269499523e-02 | ΔΩ = -1.49745e-12 | ‖∇‖ = 1.19751e-06
iter  701 | Ω_tilde = 4.269499522e-02 | ΔΩ = -1.45138e-12 | ‖∇‖ = 1.19246e-06
iter  702 | Ω_tilde = 4.269499522e-02 | ΔΩ = -1.42834e-12 | ‖∇‖ = 1.18743e-06
iter  703 | Ω_tilde = 4.269499522e-02 | ΔΩ = -1.38227e-12 | ‖∇‖ = 1.18242e-06
iter  704 | Ω_tilde = 4.269499522e-02 | ΔΩ = -1.35923e-12 | ‖∇‖ = 1.17743e-06
iter  705 | Ω_tilde = 4.269499522e-02 | ΔΩ = -1.33619e-12 | ‖∇‖ = 1.17247e-06
iter  706 | Ω_tilde = 4.269499522e-02 | ΔΩ = -1.31316e-12 | ‖∇‖ = 1.16752e-06
iter  707 | Ω_tilde = 4.269499522e-02 | ΔΩ = -1.26708e-12 | ‖∇‖ = 1.16259e-06
iter  708 | Ω_tilde = 4.269499521e-02 | ΔΩ = -1.22100e-12 | ‖∇‖ = 1.15769e-06
iter  709 | Ω_tilde = 4.269499521e-02 | ΔΩ = -1.19797e-12 | ‖∇‖ = 1.15281e-06
iter  710 | Ω_tilde = 4.269499521e-02 | ΔΩ = -1.22100e-12 | ‖∇‖ = 1.14794e-06
iter  711 | Ω_tilde = 4.269499521e-02 | ΔΩ = -1.12885e-12 | ‖∇‖ = 1.14310e-06
iter  712 | Ω_tilde = 4.269499521e-02 | ΔΩ = -1.15189e-12 | ‖∇‖ = 1.13828e-06
iter  713 | Ω_tilde = 4.269499521e-02 | ΔΩ = -1.05974e-12 | ‖∇‖ = 1.13348e-06
iter  714 | Ω_tilde = 4.269499521e-02 | ΔΩ = -1.08277e-12 | ‖∇‖ = 1.12870e-06
iter  715 | Ω_tilde = 4.269499521e-02 | ΔΩ = -1.03670e-12 | ‖∇‖ = 1.12394e-06
iter  716 | Ω_tilde = 4.269499521e-02 | ΔΩ = -1.03671e-12 | ‖∇‖ = 1.11920e-06
iter  717 | Ω_tilde = 4.269499520e-02 | ΔΩ = -9.90617e-13 | ‖∇‖ = 1.11448e-06
iter  718 | Ω_tilde = 4.269499520e-02 | ΔΩ = -9.67587e-13 | ‖∇‖ = 1.10978e-06
iter  719 | Ω_tilde = 4.269499520e-02 | ΔΩ = -9.90624e-13 | ‖∇‖ = 1.10510e-06
iter  720 | Ω_tilde = 4.269499520e-02 | ΔΩ = -8.98476e-13 | ‖∇‖ = 1.10044e-06
iter  721 | Ω_tilde = 4.269499520e-02 | ΔΩ = -9.21513e-13 | ‖∇‖ = 1.09579e-06
iter  722 | Ω_tilde = 4.269499520e-02 | ΔΩ = -8.75432e-13 | ‖∇‖ = 1.09117e-06
iter  723 | Ω_tilde = 4.269499520e-02 | ΔΩ = -8.75439e-13 | ‖∇‖ = 1.08657e-06
iter  724 | Ω_tilde = 4.269499520e-02 | ΔΩ = -8.29357e-13 | ‖∇‖ = 1.08199e-06
iter  725 | Ω_tilde = 4.269499520e-02 | ΔΩ = -8.06320e-13 | ‖∇‖ = 1.07743e-06
iter  726 | Ω_tilde = 4.269499520e-02 | ΔΩ = -8.06327e-13 | ‖∇‖ = 1.07288e-06
iter  727 | Ω_tilde = 4.269499520e-02 | ΔΩ = -7.83283e-13 | ‖∇‖ = 1.06836e-06
iter  728 | Ω_tilde = 4.269499520e-02 | ΔΩ = -7.83283e-13 | ‖∇‖ = 1.06385e-06
iter  729 | Ω_tilde = 4.269499519e-02 | ΔΩ = -7.37209e-13 | ‖∇‖ = 1.05937e-06
iter  730 | Ω_tilde = 4.269499519e-02 | ΔΩ = -7.60246e-13 | ‖∇‖ = 1.05490e-06
iter  731 | Ω_tilde = 4.269499519e-02 | ΔΩ = -7.14172e-13 | ‖∇‖ = 1.05045e-06
iter  732 | Ω_tilde = 4.269499519e-02 | ΔΩ = -7.14172e-13 | ‖∇‖ = 1.04602e-06
iter  733 | Ω_tilde = 4.269499519e-02 | ΔΩ = -6.68091e-13 | ‖∇‖ = 1.04161e-06
iter  734 | Ω_tilde = 4.269499519e-02 | ΔΩ = -6.45060e-13 | ‖∇‖ = 1.03722e-06
iter  735 | Ω_tilde = 4.269499519e-02 | ΔΩ = -6.22016e-13 | ‖∇‖ = 1.03285e-06
iter  736 | Ω_tilde = 4.269499519e-02 | ΔΩ = -6.68098e-13 | ‖∇‖ = 1.02849e-06
iter  737 | Ω_tilde = 4.269499519e-02 | ΔΩ = -6.22023e-13 | ‖∇‖ = 1.02416e-06
iter  738 | Ω_tilde = 4.269499519e-02 | ΔΩ = -6.22016e-13 | ‖∇‖ = 1.01984e-06
iter  739 | Ω_tilde = 4.269499519e-02 | ΔΩ = -5.52905e-13 | ‖∇‖ = 1.01554e-06
iter  740 | Ω_tilde = 4.269499519e-02 | ΔΩ = -5.75949e-13 | ‖∇‖ = 1.01126e-06
iter  741 | Ω_tilde = 4.269499519e-02 | ΔΩ = -5.75942e-13 | ‖∇‖ = 1.00699e-06
iter  742 | Ω_tilde = 4.269499519e-02 | ΔΩ = -5.75942e-13 | ‖∇‖ = 1.00275e-06
iter  743 | Ω_tilde = 4.269499519e-02 | ΔΩ = -5.06831e-13 | ‖∇‖ = 9.98522e-07
iter  744 | Ω_tilde = 4.269499519e-02 | ΔΩ = -5.52905e-13 | ‖∇‖ = 9.94313e-07
iter  745 | Ω_tilde = 4.269499518e-02 | ΔΩ = -5.06831e-13 | ‖∇‖ = 9.90121e-07
iter  746 | Ω_tilde = 4.269499518e-02 | ΔΩ = -4.83794e-13 | ‖∇‖ = 9.85947e-07
iter  747 | Ω_tilde = 4.269499518e-02 | ΔΩ = -4.83794e-13 | ‖∇‖ = 9.81791e-07
iter  748 | Ω_tilde = 4.269499518e-02 | ΔΩ = -4.83794e-13 | ‖∇‖ = 9.77653e-07
iter  749 | Ω_tilde = 4.269499518e-02 | ΔΩ = -4.83794e-13 | ‖∇‖ = 9.73532e-07
iter  750 | Ω_tilde = 4.269499518e-02 | ΔΩ = -4.14682e-13 | ‖∇‖ = 9.69428e-07
iter  751 | Ω_tilde = 4.269499518e-02 | ΔΩ = -4.83794e-13 | ‖∇‖ = 9.65342e-07
iter  752 | Ω_tilde = 4.269499518e-02 | ΔΩ = -4.37719e-13 | ‖∇‖ = 9.61273e-07
iter  753 | Ω_tilde = 4.269499518e-02 | ΔΩ = -4.14675e-13 | ‖∇‖ = 9.57221e-07
iter  754 | Ω_tilde = 4.269499518e-02 | ΔΩ = -4.14682e-13 | ‖∇‖ = 9.53186e-07
iter  755 | Ω_tilde = 4.269499518e-02 | ΔΩ = -3.91645e-13 | ‖∇‖ = 9.49169e-07
iter  756 | Ω_tilde = 4.269499518e-02 | ΔΩ = -4.14675e-13 | ‖∇‖ = 9.45168e-07
iter  757 | Ω_tilde = 4.269499518e-02 | ΔΩ = -3.68608e-13 | ‖∇‖ = 9.41185e-07
iter  758 | Ω_tilde = 4.269499518e-02 | ΔΩ = -3.91638e-13 | ‖∇‖ = 9.37218e-07
iter  759 | Ω_tilde = 4.269499518e-02 | ΔΩ = -2.99490e-13 | ‖∇‖ = 9.33268e-07
iter  760 | Ω_tilde = 4.269499518e-02 | ΔΩ = -3.91645e-13 | ‖∇‖ = 9.29335e-07
iter  761 | Ω_tilde = 4.269499518e-02 | ΔΩ = -3.45564e-13 | ‖∇‖ = 9.25418e-07
iter  762 | Ω_tilde = 4.269499518e-02 | ΔΩ = -3.45571e-13 | ‖∇‖ = 9.21518e-07
iter  763 | Ω_tilde = 4.269499518e-02 | ΔΩ = -3.45564e-13 | ‖∇‖ = 9.17635e-07
iter  764 | Ω_tilde = 4.269499518e-02 | ΔΩ = -2.76452e-13 | ‖∇‖ = 9.13768e-07
iter  765 | Ω_tilde = 4.269499518e-02 | ΔΩ = -3.68608e-13 | ‖∇‖ = 9.09917e-07
iter  766 | Ω_tilde = 4.269499518e-02 | ΔΩ = -2.99490e-13 | ‖∇‖ = 9.06083e-07
iter  767 | Ω_tilde = 4.269499518e-02 | ΔΩ = -2.99490e-13 | ‖∇‖ = 9.02265e-07
iter  768 | Ω_tilde = 4.269499518e-02 | ΔΩ = -2.99490e-13 | ‖∇‖ = 8.98463e-07
iter  769 | Ω_tilde = 4.269499518e-02 | ΔΩ = -2.99497e-13 | ‖∇‖ = 8.94677e-07
iter  770 | Ω_tilde = 4.269499518e-02 | ΔΩ = -2.76452e-13 | ‖∇‖ = 8.90907e-07
iter  771 | Ω_tilde = 4.269499518e-02 | ΔΩ = -2.76452e-13 | ‖∇‖ = 8.87153e-07
iter  772 | Ω_tilde = 4.269499517e-02 | ΔΩ = -2.99490e-13 | ‖∇‖ = 8.83415e-07
iter  773 | Ω_tilde = 4.269499517e-02 | ΔΩ = -2.07341e-13 | ‖∇‖ = 8.79692e-07
iter  774 | Ω_tilde = 4.269499517e-02 | ΔΩ = -2.53415e-13 | ‖∇‖ = 8.75986e-07
iter  775 | Ω_tilde = 4.269499517e-02 | ΔΩ = -2.53415e-13 | ‖∇‖ = 8.72295e-07
iter  776 | Ω_tilde = 4.269499517e-02 | ΔΩ = -2.76452e-13 | ‖∇‖ = 8.68620e-07
iter  777 | Ω_tilde = 4.269499517e-02 | ΔΩ = -2.30378e-13 | ‖∇‖ = 8.64960e-07
iter  778 | Ω_tilde = 4.269499517e-02 | ΔΩ = -2.07341e-13 | ‖∇‖ = 8.61316e-07
iter  779 | Ω_tilde = 4.269499517e-02 | ΔΩ = -2.76452e-13 | ‖∇‖ = 8.57687e-07
iter  780 | Ω_tilde = 4.269499517e-02 | ΔΩ = -2.07341e-13 | ‖∇‖ = 8.54074e-07
iter  781 | Ω_tilde = 4.269499517e-02 | ΔΩ = -1.84297e-13 | ‖∇‖ = 8.50476e-07
iter  782 | Ω_tilde = 4.269499517e-02 | ΔΩ = -2.07341e-13 | ‖∇‖ = 8.46893e-07
iter  783 | Ω_tilde = 4.269499517e-02 | ΔΩ = -2.07341e-13 | ‖∇‖ = 8.43325e-07
iter  784 | Ω_tilde = 4.269499517e-02 | ΔΩ = -2.07341e-13 | ‖∇‖ = 8.39773e-07
iter  785 | Ω_tilde = 4.269499517e-02 | ΔΩ = -2.07341e-13 | ‖∇‖ = 8.36235e-07
iter  786 | Ω_tilde = 4.269499517e-02 | ΔΩ = -2.07341e-13 | ‖∇‖ = 8.32712e-07
iter  787 | Ω_tilde = 4.269499517e-02 | ΔΩ = -1.84297e-13 | ‖∇‖ = 8.29205e-07
iter  788 | Ω_tilde = 4.269499517e-02 | ΔΩ = -1.84304e-13 | ‖∇‖ = 8.25712e-07
iter  789 | Ω_tilde = 4.269499517e-02 | ΔΩ = -1.84304e-13 | ‖∇‖ = 8.22234e-07
iter  790 | Ω_tilde = 4.269499517e-02 | ΔΩ = -1.61260e-13 | ‖∇‖ = 8.18770e-07
iter  791 | Ω_tilde = 4.269499517e-02 | ΔΩ = -1.61267e-13 | ‖∇‖ = 8.15321e-07
iter  792 | Ω_tilde = 4.269499517e-02 | ΔΩ = -1.84304e-13 | ‖∇‖ = 8.11887e-07
iter  793 | Ω_tilde = 4.269499517e-02 | ΔΩ = -1.15186e-13 | ‖∇‖ = 8.08468e-07
iter  794 | Ω_tilde = 4.269499517e-02 | ΔΩ = -2.07341e-13 | ‖∇‖ = 8.05062e-07
iter  795 | Ω_tilde = 4.269499517e-02 | ΔΩ = -1.38230e-13 | ‖∇‖ = 8.01672e-07
iter  796 | Ω_tilde = 4.269499517e-02 | ΔΩ = -1.61260e-13 | ‖∇‖ = 7.98295e-07
iter  797 | Ω_tilde = 4.269499517e-02 | ΔΩ = -1.61267e-13 | ‖∇‖ = 7.94933e-07
iter  798 | Ω_tilde = 4.269499517e-02 | ΔΩ = -1.38230e-13 | ‖∇‖ = 7.91585e-07
iter  799 | Ω_tilde = 4.269499517e-02 | ΔΩ = -1.15186e-13 | ‖∇‖ = 7.88251e-07
iter  800 | Ω_tilde = 4.269499517e-02 | ΔΩ = -1.38230e-13 | ‖∇‖ = 7.84931e-07
iter  801 | Ω_tilde = 4.269499517e-02 | ΔΩ = -1.61260e-13 | ‖∇‖ = 7.81626e-07
iter  802 | Ω_tilde = 4.269499517e-02 | ΔΩ = -1.15193e-13 | ‖∇‖ = 7.78334e-07
iter  803 | Ω_tilde = 4.269499517e-02 | ΔΩ = -1.15186e-13 | ‖∇‖ = 7.75056e-07
iter  804 | Ω_tilde = 4.269499517e-02 | ΔΩ = -1.61267e-13 | ‖∇‖ = 7.71793e-07
iter  805 | Ω_tilde = 4.269499517e-02 | ΔΩ = -9.21485e-14 | ‖∇‖ = 7.68542e-07
iter  806 | Ω_tilde = 4.269499517e-02 | ΔΩ = -1.15193e-13 | ‖∇‖ = 7.65306e-07
iter  807 | Ω_tilde = 4.269499517e-02 | ΔΩ = -1.38223e-13 | ‖∇‖ = 7.62083e-07
iter  808 | Ω_tilde = 4.269499517e-02 | ΔΩ = -1.15193e-13 | ‖∇‖ = 7.58874e-07
iter  809 | Ω_tilde = 4.269499517e-02 | ΔΩ = -1.15186e-13 | ‖∇‖ = 7.55679e-07
iter  810 | Ω_tilde = 4.269499517e-02 | ΔΩ = -6.91114e-14 | ‖∇‖ = 7.52497e-07
iter  811 | Ω_tilde = 4.269499517e-02 | ΔΩ = -1.15193e-13 | ‖∇‖ = 7.49329e-07
iter  812 | Ω_tilde = 4.269499517e-02 | ΔΩ = -1.38223e-13 | ‖∇‖ = 7.46173e-07
iter  813 | Ω_tilde = 4.269499517e-02 | ΔΩ = -9.21554e-14 | ‖∇‖ = 7.43032e-07
iter  814 | Ω_tilde = 4.269499517e-02 | ΔΩ = -9.21485e-14 | ‖∇‖ = 7.39903e-07
iter  815 | Ω_tilde = 4.269499517e-02 | ΔΩ = -9.21554e-14 | ‖∇‖ = 7.36788e-07
iter  816 | Ω_tilde = 4.269499517e-02 | ΔΩ = -1.15186e-13 | ‖∇‖ = 7.33686e-07
iter  817 | Ω_tilde = 4.269499517e-02 | ΔΩ = -6.91114e-14 | ‖∇‖ = 7.30597e-07
iter  818 | Ω_tilde = 4.269499517e-02 | ΔΩ = -1.15193e-13 | ‖∇‖ = 7.27521e-07
iter  819 | Ω_tilde = 4.269499517e-02 | ΔΩ = -9.21485e-14 | ‖∇‖ = 7.24458e-07
iter  820 | Ω_tilde = 4.269499517e-02 | ΔΩ = -4.60743e-14 | ‖∇‖ = 7.21408e-07
iter  821 | Ω_tilde = 4.269499517e-02 | ΔΩ = -9.21554e-14 | ‖∇‖ = 7.18371e-07
iter  822 | Ω_tilde = 4.269499517e-02 | ΔΩ = -9.21485e-14 | ‖∇‖ = 7.15347e-07
iter  823 | Ω_tilde = 4.269499517e-02 | ΔΩ = -6.91114e-14 | ‖∇‖ = 7.12336e-07
iter  824 | Ω_tilde = 4.269499517e-02 | ΔΩ = -1.15193e-13 | ‖∇‖ = 7.09337e-07
iter  825 | Ω_tilde = 4.269499517e-02 | ΔΩ = -4.60743e-14 | ‖∇‖ = 7.06351e-07
iter  826 | Ω_tilde = 4.269499517e-02 | ΔΩ = -6.91114e-14 | ‖∇‖ = 7.03378e-07
iter  827 | Ω_tilde = 4.269499517e-02 | ΔΩ = -6.91114e-14 | ‖∇‖ = 7.00417e-07
iter  828 | Ω_tilde = 4.269499517e-02 | ΔΩ = -9.21554e-14 | ‖∇‖ = 6.97468e-07
iter  829 | Ω_tilde = 4.269499517e-02 | ΔΩ = -6.91114e-14 | ‖∇‖ = 6.94533e-07
iter  830 | Ω_tilde = 4.269499517e-02 | ΔΩ = -6.91114e-14 | ‖∇‖ = 6.91609e-07
iter  831 | Ω_tilde = 4.269499517e-02 | ΔΩ = -6.91183e-14 | ‖∇‖ = 6.88698e-07
iter  832 | Ω_tilde = 4.269499517e-02 | ΔΩ = -6.91114e-14 | ‖∇‖ = 6.85799e-07
iter  833 | Ω_tilde = 4.269499517e-02 | ΔΩ = -6.91114e-14 | ‖∇‖ = 6.82913e-07
iter  834 | Ω_tilde = 4.269499517e-02 | ΔΩ = -2.30371e-14 | ‖∇‖ = 6.80038e-07
iter  835 | Ω_tilde = 4.269499517e-02 | ΔΩ = -4.60812e-14 | ‖∇‖ = 6.77176e-07
iter  836 | Ω_tilde = 4.269499517e-02 | ΔΩ = -6.91114e-14 | ‖∇‖ = 6.74326e-07
iter  837 | Ω_tilde = 4.269499517e-02 | ΔΩ = -4.60743e-14 | ‖∇‖ = 6.71488e-07
iter  838 | Ω_tilde = 4.269499517e-02 | ΔΩ = -6.91114e-14 | ‖∇‖ = 6.68662e-07
iter  839 | Ω_tilde = 4.269499517e-02 | ΔΩ = -6.91183e-14 | ‖∇‖ = 6.65848e-07
iter  840 | Ω_tilde = 4.269499517e-02 | ΔΩ = -4.60743e-14 | ‖∇‖ = 6.63046e-07
iter  841 | Ω_tilde = 4.269499517e-02 | ΔΩ = -4.60743e-14 | ‖∇‖ = 6.60256e-07
iter  842 | Ω_tilde = 4.269499517e-02 | ΔΩ = -6.91114e-14 | ‖∇‖ = 6.57477e-07
iter  843 | Ω_tilde = 4.269499517e-02 | ΔΩ = -4.60743e-14 | ‖∇‖ = 6.54710e-07
iter  844 | Ω_tilde = 4.269499517e-02 | ΔΩ = -6.91183e-14 | ‖∇‖ = 6.51955e-07
iter  845 | Ω_tilde = 4.269499517e-02 | ΔΩ = -2.30371e-14 | ‖∇‖ = 6.49212e-07
iter  846 | Ω_tilde = 4.269499517e-02 | ΔΩ = -4.60743e-14 | ‖∇‖ = 6.46480e-07
iter  847 | Ω_tilde = 4.269499517e-02 | ΔΩ = -4.60743e-14 | ‖∇‖ = 6.43759e-07
iter  848 | Ω_tilde = 4.269499517e-02 | ΔΩ = -4.60743e-14 | ‖∇‖ = 6.41051e-07
iter  849 | Ω_tilde = 4.269499517e-02 | ΔΩ = -2.30371e-14 | ‖∇‖ = 6.38353e-07
iter  850 | Ω_tilde = 4.269499517e-02 | ΔΩ = -6.91183e-14 | ‖∇‖ = 6.35667e-07
iter  851 | Ω_tilde = 4.269499517e-02 | ΔΩ = -2.30371e-14 | ‖∇‖ = 6.32993e-07
iter  852 | Ω_tilde = 4.269499517e-02 | ΔΩ = -4.60743e-14 | ‖∇‖ = 6.30329e-07
iter  853 | Ω_tilde = 4.269499517e-02 | ΔΩ = -6.91114e-14 | ‖∇‖ = 6.27677e-07
iter  854 | Ω_tilde = 4.269499517e-02 | ΔΩ = -2.30371e-14 | ‖∇‖ = 6.25036e-07
iter  855 | Ω_tilde = 4.269499517e-02 | ΔΩ = -4.60812e-14 | ‖∇‖ = 6.22407e-07
iter  856 | Ω_tilde = 4.269499517e-02 | ΔΩ = -2.30371e-14 | ‖∇‖ = 6.19788e-07
iter  857 | Ω_tilde = 4.269499517e-02 | ΔΩ = -2.30371e-14 | ‖∇‖ = 6.17180e-07
iter  858 | Ω_tilde = 4.269499517e-02 | ΔΩ = -4.60743e-14 | ‖∇‖ = 6.14584e-07
iter  859 | Ω_tilde = 4.269499516e-02 | ΔΩ = -4.60743e-14 | ‖∇‖ = 6.11998e-07
iter  860 | Ω_tilde = 4.269499516e-02 | ΔΩ = -2.30371e-14 | ‖∇‖ = 6.09424e-07
iter  861 | Ω_tilde = 4.269499516e-02 | ΔΩ = -4.60743e-14 | ‖∇‖ = 6.06860e-07
iter  862 | Ω_tilde = 4.269499516e-02 | ΔΩ = -2.30371e-14 | ‖∇‖ = 6.04307e-07
iter  863 | Ω_tilde = 4.269499516e-02 | ΔΩ = -2.30441e-14 | ‖∇‖ = 6.01765e-07
iter  864 | Ω_tilde = 4.269499516e-02 | ΔΩ = -2.30371e-14 | ‖∇‖ = 5.99233e-07
iter  865 | Ω_tilde = 4.269499516e-02 | ΔΩ = -2.30371e-14 | ‖∇‖ = 5.96713e-07
iter  866 | Ω_tilde = 4.269499516e-02 | ΔΩ = -2.30371e-14 | ‖∇‖ = 5.94203e-07
iter  867 | Ω_tilde = 4.269499516e-02 | ΔΩ = -2.30371e-14 | ‖∇‖ = 5.91703e-07
iter  868 | Ω_tilde = 4.269499516e-02 | ΔΩ = -4.60743e-14 | ‖∇‖ = 5.89214e-07
iter  869 | Ω_tilde = 4.269499516e-02 | ΔΩ = -2.30371e-14 | ‖∇‖ = 5.86736e-07
iter  870 | Ω_tilde = 4.269499516e-02 | ΔΩ = -2.30371e-14 | ‖∇‖ = 5.84268e-07
iter  871 | Ω_tilde = 4.269499516e-02 | ΔΩ =  0.00000e+00 | ‖∇‖ = 5.81810e-07
iter  872 | Ω_tilde = 4.269499516e-02 | ΔΩ = -4.60743e-14 | ‖∇‖ = 5.79363e-07
iter  873 | Ω_tilde = 4.269499516e-02 | ΔΩ =  0.00000e+00 | ‖∇‖ = 5.76927e-07
iter  874 | Ω_tilde = 4.269499516e-02 | ΔΩ = -2.30441e-14 | ‖∇‖ = 5.74500e-07
iter  875 | Ω_tilde = 4.269499516e-02 | ΔΩ = -4.60743e-14 | ‖∇‖ = 5.72084e-07
iter  876 | Ω_tilde = 4.269499516e-02 | ΔΩ = -2.30371e-14 | ‖∇‖ = 5.69678e-07
iter  877 | Ω_tilde = 4.269499516e-02 | ΔΩ = -4.60743e-14 | ‖∇‖ = 5.67282e-07
iter  878 | Ω_tilde = 4.269499516e-02 | ΔΩ =  0.00000e+00 | ‖∇‖ = 5.64896e-07
iter  879 | Ω_tilde = 4.269499516e-02 | ΔΩ = -2.30371e-14 | ‖∇‖ = 5.62521e-07
iter  880 | Ω_tilde = 4.269499516e-02 | ΔΩ = -4.60743e-14 | ‖∇‖ = 5.60155e-07
iter  881 | Ω_tilde = 4.269499516e-02 | ΔΩ = -2.30371e-14 | ‖∇‖ = 5.57799e-07
iter  882 | Ω_tilde = 4.269499516e-02 | ΔΩ =  0.00000e+00 | ‖∇‖ = 5.55454e-07
iter  883 | Ω_tilde = 4.269499516e-02 | ΔΩ = -2.30371e-14 | ‖∇‖ = 5.53118e-07
iter  884 | Ω_tilde = 4.269499516e-02 | ΔΩ = -2.30441e-14 | ‖∇‖ = 5.50792e-07
iter  885 | Ω_tilde = 4.269499516e-02 | ΔΩ =  2.30441e-14 | ‖∇‖ = 5.48476e-07
iter  886 | Ω_tilde = 4.269499516e-02 | ΔΩ = -4.60812e-14 | ‖∇‖ = 5.46169e-07
iter  887 | Ω_tilde = 4.269499516e-02 | ΔΩ = -2.30371e-14 | ‖∇‖ = 5.43873e-07
iter  888 | Ω_tilde = 4.269499516e-02 | ΔΩ = -2.30371e-14 | ‖∇‖ = 5.41586e-07
iter  889 | Ω_tilde = 4.269499516e-02 | ΔΩ =  2.30371e-14 | ‖∇‖ = 5.39309e-07
iter  890 | Ω_tilde = 4.269499516e-02 | ΔΩ = -4.60743e-14 | ‖∇‖ = 5.37041e-07
iter  891 | Ω_tilde = 4.269499516e-02 | ΔΩ =  2.30371e-14 | ‖∇‖ = 5.34783e-07
iter  892 | Ω_tilde = 4.269499516e-02 | ΔΩ = -4.60743e-14 | ‖∇‖ = 5.32535e-07
iter  893 | Ω_tilde = 4.269499516e-02 | ΔΩ = -4.60743e-14 | ‖∇‖ = 5.30296e-07
iter  894 | Ω_tilde = 4.269499516e-02 | ΔΩ =  0.00000e+00 | ‖∇‖ = 5.28066e-07
iter  895 | Ω_tilde = 4.269499516e-02 | ΔΩ = -2.30371e-14 | ‖∇‖ = 5.25846e-07
iter  896 | Ω_tilde = 4.269499516e-02 | ΔΩ =  2.30371e-14 | ‖∇‖ = 5.23635e-07
iter  897 | Ω_tilde = 4.269499516e-02 | ΔΩ = -2.30371e-14 | ‖∇‖ = 5.21434e-07
iter  898 | Ω_tilde = 4.269499516e-02 | ΔΩ = -2.30371e-14 | ‖∇‖ = 5.19242e-07
iter  899 | Ω_tilde = 4.269499516e-02 | ΔΩ =  0.00000e+00 | ‖∇‖ = 5.17059e-07
iter  900 | Ω_tilde = 4.269499516e-02 | ΔΩ = -2.30371e-14 | ‖∇‖ = 5.14885e-07
iter  901 | Ω_tilde = 4.269499516e-02 | ΔΩ =  0.00000e+00 | ‖∇‖ = 5.12721e-07
iter  902 | Ω_tilde = 4.269499516e-02 | ΔΩ =  0.00000e+00 | ‖∇‖ = 5.10565e-07
iter  903 | Ω_tilde = 4.269499516e-02 | ΔΩ = -2.30441e-14 | ‖∇‖ = 5.08419e-07
iter  904 | Ω_tilde = 4.269499516e-02 | ΔΩ = -4.60743e-14 | ‖∇‖ = 5.06282e-07
iter  905 | Ω_tilde = 4.269499516e-02 | ΔΩ = -2.30371e-14 | ‖∇‖ = 5.04154e-07
iter  906 | Ω_tilde = 4.269499516e-02 | ΔΩ =  2.30371e-14 | ‖∇‖ = 5.02034e-07
iter  907 | Ω_tilde = 4.269499516e-02 | ΔΩ = -2.30371e-14 | ‖∇‖ = 4.99924e-07
iter  908 | Ω_tilde = 4.269499516e-02 | ΔΩ =  0.00000e+00 | ‖∇‖ = 4.97823e-07
iter  909 | Ω_tilde = 4.269499516e-02 | ΔΩ =  2.30371e-14 | ‖∇‖ = 4.95730e-07
iter  910 | Ω_tilde = 4.269499516e-02 | ΔΩ = -4.60743e-14 | ‖∇‖ = 4.93647e-07
iter  911 | Ω_tilde = 4.269499516e-02 | ΔΩ =  2.30371e-14 | ‖∇‖ = 4.91572e-07
iter  912 | Ω_tilde = 4.269499516e-02 | ΔΩ = -4.60743e-14 | ‖∇‖ = 4.89506e-07
iter  913 | Ω_tilde = 4.269499516e-02 | ΔΩ = -2.30371e-14 | ‖∇‖ = 4.87448e-07
iter  914 | Ω_tilde = 4.269499516e-02 | ΔΩ =  0.00000e+00 | ‖∇‖ = 4.85400e-07
iter  915 | Ω_tilde = 4.269499516e-02 | ΔΩ =  0.00000e+00 | ‖∇‖ = 4.83360e-07
iter  916 | Ω_tilde = 4.269499516e-02 | ΔΩ = -2.30371e-14 | ‖∇‖ = 4.81328e-07
iter  917 | Ω_tilde = 4.269499516e-02 | ΔΩ =  0.00000e+00 | ‖∇‖ = 4.79306e-07
iter  918 | Ω_tilde = 4.269499516e-02 | ΔΩ =  0.00000e+00 | ‖∇‖ = 4.77291e-07
iter  919 | Ω_tilde = 4.269499516e-02 | ΔΩ = -2.30371e-14 | ‖∇‖ = 4.75285e-07
iter  920 | Ω_tilde = 4.269499516e-02 | ΔΩ =  2.30371e-14 | ‖∇‖ = 4.73288e-07
iter  921 | Ω_tilde = 4.269499516e-02 | ΔΩ = -4.60743e-14 | ‖∇‖ = 4.71299e-07
iter  922 | Ω_tilde = 4.269499516e-02 | ΔΩ =  2.30371e-14 | ‖∇‖ = 4.69319e-07
iter  923 | Ω_tilde = 4.269499516e-02 | ΔΩ =  0.00000e+00 | ‖∇‖ = 4.67347e-07
iter  924 | Ω_tilde = 4.269499516e-02 | ΔΩ =  0.00000e+00 | ‖∇‖ = 4.65383e-07
iter  925 | Ω_tilde = 4.269499516e-02 | ΔΩ = -2.30371e-14 | ‖∇‖ = 4.63427e-07
iter  926 | Ω_tilde = 4.269499516e-02 | ΔΩ = -2.30371e-14 | ‖∇‖ = 4.61480e-07
iter  927 | Ω_tilde = 4.269499516e-02 | ΔΩ = -2.30441e-14 | ‖∇‖ = 4.59541e-07
iter  928 | Ω_tilde = 4.269499516e-02 | ΔΩ =  0.00000e+00 | ‖∇‖ = 4.57610e-07
iter  929 | Ω_tilde = 4.269499516e-02 | ΔΩ =  0.00000e+00 | ‖∇‖ = 4.55687e-07
iter  930 | Ω_tilde = 4.269499516e-02 | ΔΩ =  0.00000e+00 | ‖∇‖ = 4.53772e-07
iter  931 | Ω_tilde = 4.269499516e-02 | ΔΩ =  0.00000e+00 | ‖∇‖ = 4.51866e-07
iter  932 | Ω_tilde = 4.269499516e-02 | ΔΩ =  0.00000e+00 | ‖∇‖ = 4.49967e-07
iter  933 | Ω_tilde = 4.269499516e-02 | ΔΩ = -2.30371e-14 | ‖∇‖ = 4.48077e-07
iter  934 | Ω_tilde = 4.269499516e-02 | ΔΩ =  0.00000e+00 | ‖∇‖ = 4.46194e-07
iter  935 | Ω_tilde = 4.269499516e-02 | ΔΩ = -2.30371e-14 | ‖∇‖ = 4.44320e-07
iter  936 | Ω_tilde = 4.269499516e-02 | ΔΩ =  2.30371e-14 | ‖∇‖ = 4.42453e-07
iter  937 | Ω_tilde = 4.269499516e-02 | ΔΩ = -2.30371e-14 | ‖∇‖ = 4.40594e-07
iter  938 | Ω_tilde = 4.269499516e-02 | ΔΩ =  0.00000e+00 | ‖∇‖ = 4.38744e-07
iter  939 | Ω_tilde = 4.269499516e-02 | ΔΩ =  0.00000e+00 | ‖∇‖ = 4.36900e-07
iter  940 | Ω_tilde = 4.269499516e-02 | ΔΩ =  0.00000e+00 | ‖∇‖ = 4.35065e-07
iter  941 | Ω_tilde = 4.269499516e-02 | ΔΩ =  0.00000e+00 | ‖∇‖ = 4.33237e-07
iter  942 | Ω_tilde = 4.269499516e-02 | ΔΩ =  0.00000e+00 | ‖∇‖ = 4.31417e-07
iter  943 | Ω_tilde = 4.269499516e-02 | ΔΩ = -2.30371e-14 | ‖∇‖ = 4.29605e-07
iter  944 | Ω_tilde = 4.269499516e-02 | ΔΩ =  0.00000e+00 | ‖∇‖ = 4.27801e-07
iter  945 | Ω_tilde = 4.269499516e-02 | ΔΩ = -2.30371e-14 | ‖∇‖ = 4.26004e-07
iter  946 | Ω_tilde = 4.269499516e-02 | ΔΩ =  0.00000e+00 | ‖∇‖ = 4.24214e-07
iter  947 | Ω_tilde = 4.269499516e-02 | ΔΩ =  0.00000e+00 | ‖∇‖ = 4.22432e-07
iter  948 | Ω_tilde = 4.269499516e-02 | ΔΩ =  2.30371e-14 | ‖∇‖ = 4.20658e-07
iter  949 | Ω_tilde = 4.269499516e-02 | ΔΩ =  0.00000e+00 | ‖∇‖ = 4.18891e-07
iter  950 | Ω_tilde = 4.269499516e-02 | ΔΩ = -2.30371e-14 | ‖∇‖ = 4.17132e-07
iter  951 | Ω_tilde = 4.269499516e-02 | ΔΩ =  0.00000e+00 | ‖∇‖ = 4.15380e-07
iter  952 | Ω_tilde = 4.269499516e-02 | ΔΩ = -2.30371e-14 | ‖∇‖ = 4.13635e-07
iter  953 | Ω_tilde = 4.269499516e-02 | ΔΩ =  0.00000e+00 | ‖∇‖ = 4.11898e-07
iter  954 | Ω_tilde = 4.269499516e-02 | ΔΩ =  2.30371e-14 | ‖∇‖ = 4.10168e-07
iter  955 | Ω_tilde = 4.269499516e-02 | ΔΩ =  0.00000e+00 | ‖∇‖ = 4.08446e-07
iter  956 | Ω_tilde = 4.269499516e-02 | ΔΩ = -2.30371e-14 | ‖∇‖ = 4.06730e-07
iter  957 | Ω_tilde = 4.269499516e-02 | ΔΩ =  0.00000e+00 | ‖∇‖ = 4.05022e-07
iter  958 | Ω_tilde = 4.269499516e-02 | ΔΩ = -2.30371e-14 | ‖∇‖ = 4.03321e-07
iter  959 | Ω_tilde = 4.269499516e-02 | ΔΩ =  0.00000e+00 | ‖∇‖ = 4.01628e-07
iter  960 | Ω_tilde = 4.269499516e-02 | ΔΩ =  0.00000e+00 | ‖∇‖ = 3.99941e-07
iter  961 | Ω_tilde = 4.269499516e-02 | ΔΩ =  2.30371e-14 | ‖∇‖ = 3.98261e-07
iter  962 | Ω_tilde = 4.269499516e-02 | ΔΩ = -2.30371e-14 | ‖∇‖ = 3.96589e-07
iter  963 | Ω_tilde = 4.269499516e-02 | ΔΩ =  0.00000e+00 | ‖∇‖ = 3.94924e-07
iter  964 | Ω_tilde = 4.269499516e-02 | ΔΩ =  0.00000e+00 | ‖∇‖ = 3.93266e-07
iter  965 | Ω_tilde = 4.269499516e-02 | ΔΩ = -2.30371e-14 | ‖∇‖ = 3.91614e-07
iter  966 | Ω_tilde = 4.269499516e-02 | ΔΩ =  0.00000e+00 | ‖∇‖ = 3.89970e-07
iter  967 | Ω_tilde = 4.269499516e-02 | ΔΩ =  2.30371e-14 | ‖∇‖ = 3.88332e-07
iter  968 | Ω_tilde = 4.269499516e-02 | ΔΩ =  0.00000e+00 | ‖∇‖ = 3.86702e-07
iter  969 | Ω_tilde = 4.269499516e-02 | ΔΩ = -2.30371e-14 | ‖∇‖ = 3.85078e-07
iter  970 | Ω_tilde = 4.269499516e-02 | ΔΩ =  0.00000e+00 | ‖∇‖ = 3.83462e-07
iter  971 | Ω_tilde = 4.269499516e-02 | ΔΩ = -2.30371e-14 | ‖∇‖ = 3.81852e-07
iter  972 | Ω_tilde = 4.269499516e-02 | ΔΩ =  2.30371e-14 | ‖∇‖ = 3.80248e-07
iter  973 | Ω_tilde = 4.269499516e-02 | ΔΩ = -2.30371e-14 | ‖∇‖ = 3.78652e-07
iter  974 | Ω_tilde = 4.269499516e-02 | ΔΩ =  2.30371e-14 | ‖∇‖ = 3.77062e-07
iter  975 | Ω_tilde = 4.269499516e-02 | ΔΩ =  0.00000e+00 | ‖∇‖ = 3.75479e-07
iter  976 | Ω_tilde = 4.269499516e-02 | ΔΩ = -2.30371e-14 | ‖∇‖ = 3.73903e-07
iter  977 | Ω_tilde = 4.269499516e-02 | ΔΩ =  2.30371e-14 | ‖∇‖ = 3.72333e-07
iter  978 | Ω_tilde = 4.269499516e-02 | ΔΩ =  0.00000e+00 | ‖∇‖ = 3.70770e-07
iter  979 | Ω_tilde = 4.269499516e-02 | ΔΩ = -2.30371e-14 | ‖∇‖ = 3.69214e-07
iter  980 | Ω_tilde = 4.269499516e-02 | ΔΩ =  2.30371e-14 | ‖∇‖ = 3.67664e-07
iter  981 | Ω_tilde = 4.269499516e-02 | ΔΩ =  0.00000e+00 | ‖∇‖ = 3.66121e-07
iter  982 | Ω_tilde = 4.269499516e-02 | ΔΩ = -4.60743e-14 | ‖∇‖ = 3.64584e-07
iter  983 | Ω_tilde = 4.269499516e-02 | ΔΩ =  2.30371e-14 | ‖∇‖ = 3.63053e-07
iter  984 | Ω_tilde = 4.269499516e-02 | ΔΩ =  0.00000e+00 | ‖∇‖ = 3.61529e-07
iter  985 | Ω_tilde = 4.269499516e-02 | ΔΩ =  2.30371e-14 | ‖∇‖ = 3.60012e-07
iter  986 | Ω_tilde = 4.269499516e-02 | ΔΩ = -4.60743e-14 | ‖∇‖ = 3.58501e-07
iter  987 | Ω_tilde = 4.269499516e-02 | ΔΩ =  2.30371e-14 | ‖∇‖ = 3.56996e-07
iter  988 | Ω_tilde = 4.269499516e-02 | ΔΩ =  0.00000e+00 | ‖∇‖ = 3.55498e-07
iter  989 | Ω_tilde = 4.269499516e-02 | ΔΩ =  0.00000e+00 | ‖∇‖ = 3.54006e-07
iter  990 | Ω_tilde = 4.269499516e-02 | ΔΩ =  2.30371e-14 | ‖∇‖ = 3.52520e-07
iter  991 | Ω_tilde = 4.269499516e-02 | ΔΩ = -2.30371e-14 | ‖∇‖ = 3.51040e-07
iter  992 | Ω_tilde = 4.269499516e-02 | ΔΩ = -4.60812e-14 | ‖∇‖ = 3.49567e-07
iter  993 | Ω_tilde = 4.269499516e-02 | ΔΩ =  4.60812e-14 | ‖∇‖ = 3.48100e-07
iter  994 | Ω_tilde = 4.269499516e-02 | ΔΩ = -2.30371e-14 | ‖∇‖ = 3.46639e-07
iter  995 | Ω_tilde = 4.269499516e-02 | ΔΩ = -2.30441e-14 | ‖∇‖ = 3.45184e-07
iter  996 | Ω_tilde = 4.269499516e-02 | ΔΩ =  0.00000e+00 | ‖∇‖ = 3.43736e-07
iter  997 | Ω_tilde = 4.269499516e-02 | ΔΩ =  2.30441e-14 | ‖∇‖ = 3.42293e-07
iter  998 | Ω_tilde = 4.269499516e-02 | ΔΩ =  0.00000e+00 | ‖∇‖ = 3.40857e-07
iter  999 | Ω_tilde = 4.269499516e-02 | ΔΩ = -2.30441e-14 | ‖∇‖ = 3.39426e-07
WF.info()
Wannier Function Report
===============================================================
WF 1: center = [0.48350366, 0.27915095]  Omega     = 0.23352883
WF 2: center = [1.00000001, 1.17374889]  Omega     = 0.23352839
WF 3: center = [1.51649634, 0.27915097]  Omega     = 0.23353039
Sum : center = [3.00000001, 1.73205081]  Omega tot = 0.70058762
---------------------------------------------------------------
Omega I   = 0.65789262
Omega D   = 0.00682290
Omega OD  = 0.03587209
Omega tot = 0.70058762

Now the spreads have been minimized, and the Wannier functions are maximally localized. To help validate that the Wannier functions are indeed exponentially localized, we can plot the decay of each Wannier function’s weight away from its center with plot_decay. This will plot the absolute value of each Wannier function as a function of distance from its center on a logarithmic scale.

fig, ax = WF.plot_decay(0, show=True)
../_images/73c6fed9d18507dff58a161dc131e641cf1cd76e93797899e423c8198f65df7c.png
fig, ax = WF.plot_density(0, show=True)
../_images/d0bfb218afb8e9b4f715f4183ad9f711f95450e30966089a815d3c46332422e9.png

Note that we have effectively broken the primitive translational symmetry of the underlying lattice by choosing a subset of trial wavefunctions on three out of the four low energy sites in the supercell. We can see their positions using plot_centers

fig, ax = WF.plot_centers(
    color_home_cell=True, center_scale=15, legend=True, pmx=4, pmy=4, show=True
)
../_images/edd2afe14ebeabf03531f814dd34f856e52ff9e70cd33ddf695fcd397e96717a.png

Wannier interpolation#

We can view the Wannier interpolated bands by calling plot_interp_bands. We specify a set of high-symmetry k-points that defines the one-dimensional path along which the bands are plotted.

k_nodes = [
    [0, 0],
    [2 / 3, 1 / 3],
    [1 / 2, 1 / 2],
    [1 / 3, 2 / 3],
    [0, 0],
    [1 / 2, 1 / 2],
]
k_label = (r"$\Gamma $", r"$K$", r"$M$", r"$K^\prime$", r"$\Gamma $", r"$M$")
n_interp = 501
interp_energies = WF.interp_bands(k_nodes, n_interp=n_interp, ret_eigvecs=False)
fig, ax = model.plot_bands(
    k_nodes=k_nodes,
    nk=501,
    k_node_labels=k_label,
    proj_orb_idx=high_E_sites,
    cmap="plasma",
)

(k_vec, k_dist, k_node) = model.k_path(k_nodes, nk=n_interp, report=False)
ax.plot(k_dist, interp_energies, ls="--", c="lightgreen", lw=2, zorder=5, alpha=1)

# plot windows
if frozen_window is not None:
    ax.axhline(frozen_window[0], ls="--", c="b", label="frozen window")
    ax.axhline(frozen_window[1], ls="--", c="b")

ax.axhline(outer_window[0], ls=":", c="r", label="disentanglement window")
ax.axhline(outer_window[1], ls=":", c="r")
ax.legend()
<matplotlib.legend.Legend at 0x766f5991aa80>
../_images/bc8606baeb0a1f2578c24513f1804e18ca4e8eb6380ff8987d12611cc4d9ef9d.png