Incompressible Navier-Stokes with heat conduction and convection

Download this notebook

Download meshes here

# Video generated with Paraview from the output files
from IPython.display import Video
Video("visualizations/navier_stokes_temperature.mp4", embed=True, width=800, height=400)

import os

import jax
import jax.numpy as jnp
import meshio
from flax.core import FrozenDict

from autopdex import utility, spaces, seeder, models, dae

jax.config.update("jax_enable_x64", True)

Load the mesh that was generated with GMSH

mesh = meshio.read("meshes\navier_stokes_mesh_v.msh")
coords_v = jnp.asarray(mesh.points[:, :2])
cells_v = jnp.asarray(mesh.cells_dict["triangle6"])

mesh = meshio.read("meshes\navier_stokes_mesh_p.msh")
coords_p = jnp.asarray(mesh.points[:, :2])
cells_p = jnp.asarray(mesh.cells_dict["triangle"])

coords_T = coords_v
cells_T = cells_v

node_coordinates = {
  '1velocity': coords_v,
  '2pressure': coords_p,
  '3temperature': coords_T,
}

Boundary conditions:

• Parabolic inflow at x = 0,

• No-slip on the cylinder and top/bottom walls.

• Inflow (x = 0): T = 0 (Dirichlet)

• Cylinder: T ramps from 0 to 20 over the first 3.5 seconds (Dirichlet)

• Top and bottom walls: insulated (Neumann, natural condition)

• p = 0 at lower right corner,

def on_left(x):
  return jnp.isclose(x[:, 0], 0.0)
def on_cylinder(x):
  return jnp.isclose((x[:, 0] - 0.2)**2 + (x[:, 1] - 0.2)**2, 0.05**2)
def on_lower_right(x):
  return jnp.isclose(x, jnp.array([2.2, 0.])).all(axis=1)
def on_top_bottom_circle(x):
  return (jnp.isclose(x[:, 1], 0.0) + jnp.isclose(x[:, 1], 0.41) + jnp.isclose(
      (x[:, 0] - 0.2)**2 + (x[:, 1] - 0.2)**2, 0.05**2))

selection_v_left = utility.dof_select(on_left(coords_v), jnp.asarray([True, True]))
selection_v_remaining = utility.dof_select(on_top_bottom_circle(coords_v), jnp.asarray([True, True]))
selection_p = utility.dof_select(on_lower_right(coords_p), True)
selection_temp_left = utility.dof_select(on_left(coords_T), True)
selection_temp_cylinder = utility.dof_select(on_cylinder(coords_T), True)

dirichlet_dofs = {
  '1velocity': selection_v_left + selection_v_remaining,
  '2pressure': selection_p,
  '3temperature': selection_temp_left + selection_temp_cylinder,
}

# Set initial Dirichlet conditions
dirichlet_conditions = utility.dict_zeros_like(dirichlet_dofs, dtype=jnp.float64)

# Time-dependent Dirichlet condition for the velocity inflow
def get_dirichlet_conditions(t):

  def condition(y_coor):
    u_max = 1.5
    u = 4 * u_max * y_coor * (0.41 - y_coor) / 0.41**2

    # 3.5 seconds sinusoidal ramp-up, then constant
    cond_list = [t <= 3.5, t > 3.5]
    func_list = [u * jnp.sin(jnp.pi / 7 * t), u]
    return jnp.asarray([jnp.piecewise(t, cond_list, func_list), 0.])

  new_values = jax.vmap(condition)(coords_v[:, 1])
  return jnp.where(selection_v_left, new_values, 0.0)

# This function is called before each time step to update the boundary conditions
def pre_step_updates(t, settings):
  settings['current time'] = t

  # Update velocity boundary conditions at time t
  new_dirichlet = get_dirichlet_conditions(t)
  settings['dirichlet conditions']['1velocity'] = new_dirichlet

  # Update temperature boundary conditions at time t
  # Cylinder temperature: ramp from 0 to 20 within the first 3.5 seconds
  cylinder_temp = jnp.where(t <= 3.5, 20 * jnp.sin(jnp.pi / 7 * t), 20.0)
  new_temp_bc = jnp.zeros_like(dirichlet_dofs['3temperature'], dtype=jnp.float64)
  new_temp_bc = new_temp_bc.at[selection_temp_cylinder].set(cylinder_temp)
  settings['dirichlet conditions']['3temperature'] = new_temp_bc

  return settings

Weak formulation: Coupled equations (Navier–Stokes and energy balance)

def weak_form_integrand(x_int, trial_ansatz, test_ansatz, settings, static_settings, elem_number, set):
  t = settings['current time']

  # Pressure (p) and test function (q)
  p = trial_ansatz['2pressure'](x_int, t)
  q = test_ansatz['2pressure'](x_int)

  # Velocity (u) and its derivatives
  u = trial_ansatz['1velocity'](x_int, t)
  grad_u, du_dt = jax.jacfwd(trial_ansatz['1velocity'], (0, 1))(x_int, t)

  # Test function (v) and its derivatives
  v = test_ansatz['1velocity'](x_int)
  grad_v = jax.jacfwd(test_ansatz['1velocity'])(x_int)

  # Temperature (T) and test function (φ)
  T_val = trial_ansatz['3temperature'](x_int, t)
  grad_T, dT_dt = jax.jacfwd(trial_ansatz['3temperature'], (0, 1))(x_int, t)
  phi, grad_phi = jax.value_and_grad(test_ansatz['3temperature'])(x_int)

  # Divergence of u and v
  div_u = grad_u[0, 0] + grad_u[1, 1]
  div_v = grad_v[0, 0] + grad_v[1, 1]

  # Temperature-dependent viscosity:
  # Base viscosity and temperature coefficient
  nu0 = 0.001
  beta = 0.05
  nu = nu0 * (1 + beta * T_val)

  # Weak form of the momentum equation
  weak_momentum = (du_dt @ v + jnp.einsum('j,ij->i', u, grad_u) @ v + nu * jnp.einsum('ij,ij->', grad_u, grad_v) -
                   p * div_v)

  # Weak form of the continuity equation
  weak_continuity = -q * div_u

  # Thermal diffusivity parameter (α)
  alpha = 0.001

  # Weak form of the energy (temperature) equation
  weak_energy = (dT_dt * phi + (u @ grad_T) * phi + alpha * jnp.dot(grad_T, grad_phi))

  return weak_momentum + weak_continuity + weak_energy

Set-up the element and time stepping procedure

# Ansatz functions for the fields
ansatz_fun = {
  '1velocity': spaces.fem_iso_line_tri_tet,  # quadratic (6 nodes per triangle)
  '2pressure': spaces.fem_iso_line_tri_tet,  # linear (3 nodes per triangle)
  '3temperature': spaces.fem_iso_line_tri_tet,
}

# Numerical integration
ref_int_coor, ref_int_weights = seeder.int_pts_ref_tri(order=4)

# Set-up the 'user residual' based on a weak-in-space formulation of the PDEs.
user_residual = models.mixed_reference_domain_residual_time(
    weak_form_integrand,
    ansatz_fun,
    ref_int_coor,
    ref_int_weights,
    mapping_key='1velocity',
)

# The settings are analougs to static analyses. Additionally, the 'current time' has to be initialized.
settings = {
  'connectivity': ({
    '1velocity': cells_v,
    '2pressure': cells_p,
    '3temperature': cells_T,
  },),
  'node coordinates': node_coordinates,
  'dirichlet dofs': dirichlet_dofs,
  'dirichlet conditions': dirichlet_conditions,
  'current time': 0.0,
}

# The static settings are also passed to the PDE functions.
# When used through the time stepping manager, also the 'dae' and 'time integrators' have to be specified.
static_settings = FrozenDict({
  'assembling mode': ('user residual',),
  'solution structure': ('nodal imposition',),
  'model': (user_residual,),
  'solver type': 'newton',
  'solver backend': 'pardiso',
  'solver': 'lu',
  'verbose': 1,
  'dae': 'call pde',    # Here we need the keyword 'call pde' in order to call the assembling routines for PDE specified by the 'model'.
  'time integrators': {
    '1velocity': dae.BackwardDiffFormula(2),
    '2pressure': dae.BackwardEuler(),     # The time derivative is not used, but an integrator has to be specified,
                                          #   that is compatible with the stages of the other fields
    '3temperature': dae.AdamsMoulton(1),
  },
})

# Define the postprocessing policy and pass the function for pre-step updates
manager = dae.TimeSteppingManager(static_settings,
                              save_policy=dae.SaveAllPolicy(),
                              pre_step_updates=pre_step_updates)

# Start with zero velocity and temperature
initial_conditions = {
  '1velocity': jnp.zeros((coords_v.shape[0], 2)),
  '2pressure': jnp.zeros((coords_p.shape[0], )),
  '3temperature': jnp.zeros((coords_T.shape[0], )),
}

Run the time stepping loop

t_final = 10
num_time_steps = 500
result = manager.run(initial_conditions, t_final / num_time_steps, t_final, num_time_steps, settings)

Iteration 0, Residual norm: 0.005942067640679757
Iteration 1, Residual norm: 3.083083564174187e-06
Iteration 2, Residual norm: 2.0176470260160768e-11

Iteration 0, Residual norm: 0.011434755869271162
Iteration 1, Residual norm: 3.1790372706762684e-06
Iteration 2, Residual norm: 2.0711606733126938e-11

Iteration 0, Residual norm: 0.00937001504683964
Iteration 1, Residual norm: 3.246223061661343e-06
Iteration 2, Residual norm: 1.92156329092327e-11

Iteration 0, Residual norm: 0.009560302298064465
Iteration 1, Residual norm: 3.295802526767916e-06
Iteration 2, Residual norm: 1.7878385789079422e-11

Iteration 0, Residual norm: 0.009423944133330217
Iteration 1, Residual norm: 3.3406106410982776e-06
Iteration 2, Residual norm: 1.698135650782057e-11

Iteration 0, Residual norm: 0.009582405843599125
Iteration 1, Residual norm: 3.3840118939866265e-06
Iteration 2, Residual norm: 1.637922174459858e-11

Iteration 0, Residual norm: 0.009465234837007126
Iteration 1, Residual norm: 3.426353714204454e-06
Iteration 2, Residual norm: 1.5932583358916965e-11

Iteration 0, Residual norm: 0.009602299192699582
Iteration 1, Residual norm: 3.468362844678119e-06
Iteration 2, Residual norm: 1.5557872845498258e-11

Iteration 0, Residual norm: 0.009499080365139594
Iteration 1, Residual norm: 3.5102368864519675e-06
Iteration 2, Residual norm: 1.5205549167374453e-11

Iteration 0, Residual norm: 0.009620353203828313
Iteration 1, Residual norm: 3.552088555828836e-06
Iteration 2, Residual norm: 1.4842556544569422e-11

Iteration 0, Residual norm: 0.009528201096226157
Iteration 1, Residual norm: 3.5938953190744814e-06
Iteration 2, Residual norm: 1.4461062135875258e-11

Iteration 0, Residual norm: 0.009637216139047883
Iteration 1, Residual norm: 3.6355556594700234e-06
Iteration 2, Residual norm: 1.4063105460082403e-11

Iteration 0, Residual norm: 0.00955440796003213
Iteration 1, Residual norm: 3.676912928783038e-06
Iteration 2, Residual norm: 1.3667600297644581e-11

Iteration 0, Residual norm: 0.009653783663757264
Iteration 1, Residual norm: 3.717852416037243e-06
Iteration 2, Residual norm: 1.3292711104811013e-11

Iteration 0, Residual norm: 0.00957934393395053
Iteration 1, Residual norm: 3.7583062542686667e-06
Iteration 2, Residual norm: 1.2963083034060339e-11

Iteration 0, Residual norm: 0.009671249915480505
Iteration 1, Residual norm: 3.7983949200369786e-06
Iteration 2, Residual norm: 1.269324136153106e-11
Progress: 3%, Time: 3.20e-01, accepted step: True, dt: 2.00e-02, iterations: 2

Iteration 0, Residual norm: 0.00960474587143688
Iteration 1, Residual norm: 3.83838760873752e-06
Iteration 2, Residual norm: 1.2497832708448826e-11

Iteration 0, Residual norm: 0.009691132926760711
Iteration 1, Residual norm: 3.8788152469234815e-06
Iteration 2, Residual norm: 1.237727361752124e-11

Iteration 0, Residual norm: 0.009632561071612
Iteration 1, Residual norm: 3.92035324777269e-06
Iteration 2, Residual norm: 1.2331945978711878e-11

Iteration 0, Residual norm: 0.00971528990396227
Iteration 1, Residual norm: 3.963846361816647e-06
Iteration 2, Residual norm: 1.2349682534325384e-11

Iteration 0, Residual norm: 0.009665009917747517
Iteration 1, Residual norm: 4.010111573409758e-06
Iteration 2, Residual norm: 1.2423276177819488e-11

Iteration 0, Residual norm: 0.009745935846707767
Iteration 1, Residual norm: 4.059901970588883e-06
Iteration 2, Residual norm: 1.2537270070646944e-11

Iteration 0, Residual norm: 0.009704637181814512
Iteration 1, Residual norm: 4.113720221745619e-06
Iteration 2, Residual norm: 1.2686320438629631e-11

Iteration 0, Residual norm: 0.009785674493879626
Iteration 1, Residual norm: 4.1718416805773e-06
Iteration 2, Residual norm: 1.2857663330101479e-11

Iteration 0, Residual norm: 0.009754357593915252
Iteration 1, Residual norm: 4.234252867465582e-06
Iteration 2, Residual norm: 1.3048919594612684e-11

Iteration 0, Residual norm: 0.009837515775075953
Iteration 1, Residual norm: 4.300787119960269e-06
Iteration 2, Residual norm: 1.324526147379282e-11

Iteration 0, Residual norm: 0.009817452893787333
Iteration 1, Residual norm: 4.371142151636899e-06
Iteration 2, Residual norm: 1.343952902087998e-11

Iteration 0, Residual norm: 0.00990481640549892
Iteration 1, Residual norm: 4.444996243502149e-06
Iteration 2, Residual norm: 1.3607753238641657e-11

Iteration 0, Residual norm: 0.009897456004102254
Iteration 1, Residual norm: 4.5219439429365535e-06
Iteration 2, Residual norm: 1.3736196603738689e-11

Iteration 0, Residual norm: 0.009991089827118578
Iteration 1, Residual norm: 4.601495641212458e-06
Iteration 2, Residual norm: 1.3798173592002715e-11

Iteration 0, Residual norm: 0.009997896874198064
Iteration 1, Residual norm: 4.682955518380378e-06
Iteration 2, Residual norm: 1.3786396431552196e-11

Iteration 0, Residual norm: 0.01009969175804651
Iteration 1, Residual norm: 4.765446326449205e-06
Iteration 2, Residual norm: 1.3688862735853463e-11

Iteration 0, Residual norm: 0.010121947048339082
Iteration 1, Residual norm: 4.847903140722904e-06
Iteration 2, Residual norm: 1.3518250126869838e-11

Iteration 0, Residual norm: 0.010233436345046323
Iteration 1, Residual norm: 4.929226312768863e-06
Iteration 2, Residual norm: 1.3283840147077372e-11

Iteration 0, Residual norm: 0.010272028470130937
Iteration 1, Residual norm: 5.008363496805256e-06
Iteration 2, Residual norm: 1.3011141812166564e-11

Iteration 0, Residual norm: 0.010394216586094615
Iteration 1, Residual norm: 5.084462923391846e-06
Iteration 2, Residual norm: 1.2716459189980691e-11

Iteration 0, Residual norm: 0.010449480402856419
Iteration 1, Residual norm: 5.156885916204653e-06
Iteration 2, Residual norm: 1.2427409369645346e-11

Iteration 0, Residual norm: 0.01058275633170828
Iteration 1, Residual norm: 5.225226238494636e-06
Iteration 2, Residual norm: 1.2172839280113758e-11

Iteration 0, Residual norm: 0.010654448553998363
Iteration 1, Residual norm: 5.289159451465758e-06
Iteration 2, Residual norm: 1.1989931590559946e-11

Iteration 0, Residual norm: 0.010798666765228974
Iteration 1, Residual norm: 5.348323348432126e-06
Iteration 2, Residual norm: 1.1895460950803366e-11
Progress: 8%, Time: 8.00e-01, accepted step: True, dt: 2.00e-02, iterations: 2

Iteration 0, Residual norm: 0.010886113823408667
Iteration 1, Residual norm: 5.402161769000609e-06
Iteration 2, Residual norm: 1.1883288869593874e-11

Iteration 0, Residual norm: 0.011040779117876648
Iteration 1, Residual norm: 5.4499452899271745e-06
Iteration 2, Residual norm: 1.1925092000482914e-11

Iteration 0, Residual norm: 0.01114300432376483
Iteration 1, Residual norm: 5.490793801373466e-06
Iteration 2, Residual norm: 1.1995428923897384e-11

Iteration 0, Residual norm: 0.011307253316715862
Iteration 1, Residual norm: 5.523767178013513e-06
Iteration 2, Residual norm: 1.206690259379299e-11

Iteration 0, Residual norm: 0.011422737682424572
Iteration 1, Residual norm: 5.547999262600308e-06
Iteration 2, Residual norm: 1.211493021371268e-11

Iteration 0, Residual norm: 0.011594884319932285
Iteration 1, Residual norm: 5.563207010436393e-06
Iteration 2, Residual norm: 1.2115021523593506e-11

Iteration 0, Residual norm: 0.011721003351061417
Iteration 1, Residual norm: 5.570462293890364e-06
Iteration 2, Residual norm: 1.2052206915900782e-11

Iteration 0, Residual norm: 0.011898095524228948
Iteration 1, Residual norm: 5.5725212547122275e-06
Iteration 2, Residual norm: 1.1899047538793475e-11

Iteration 0, Residual norm: 0.01203105733111836
Iteration 1, Residual norm: 5.572653672153008e-06
Iteration 2, Residual norm: 1.1638519602561489e-11

Iteration 0, Residual norm: 0.012209383935132557
Iteration 1, Residual norm: 5.5724160347919755e-06
Iteration 2, Residual norm: 1.1298503272998808e-11

Iteration 0, Residual norm: 0.012345245329472964
Iteration 1, Residual norm: 5.570907583891123e-06
Iteration 2, Residual norm: 1.0941918307232853e-11

Iteration 0, Residual norm: 0.012521511345133083
Iteration 1, Residual norm: 5.567251160847854e-06
Iteration 2, Residual norm: 1.0604234657973843e-11

Iteration 0, Residual norm: 0.012657064247259311
Iteration 1, Residual norm: 5.564149915069339e-06
Iteration 2, Residual norm: 1.0296085750263001e-11

Iteration 0, Residual norm: 0.012828707174882211
Iteration 1, Residual norm: 5.5684136351986575e-06
Iteration 2, Residual norm: 1.0026345172256312e-11

Iteration 0, Residual norm: 0.012961389556542977
Iteration 1, Residual norm: 5.58736986828839e-06
Iteration 2, Residual norm: 9.797202588807951e-12

Iteration 0, Residual norm: 0.01312651386589791
Iteration 1, Residual norm: 5.623865148725294e-06
Iteration 2, Residual norm: 9.602721267053303e-12

Iteration 0, Residual norm: 0.013254756684964682
Iteration 1, Residual norm: 5.674136461640454e-06
Iteration 2, Residual norm: 9.441939823616639e-12

Iteration 0, Residual norm: 0.013412858496113619
Iteration 1, Residual norm: 5.730717816823645e-06
Iteration 2, Residual norm: 9.288251675941507e-12

Iteration 0, Residual norm: 0.013536750655814214
Iteration 1, Residual norm: 5.787727518778484e-06
Iteration 2, Residual norm: 9.096978077320147e-12

Iteration 0, Residual norm: 0.013688772879634627
Iteration 1, Residual norm: 5.843251169851383e-06
Iteration 2, Residual norm: 8.855389658624439e-12

Iteration 0, Residual norm: 0.013809308655788293
Iteration 1, Residual norm: 5.89785190503211e-06
Iteration 2, Residual norm: 8.58633899728263e-12

Iteration 0, Residual norm: 0.013956365885706605
Iteration 1, Residual norm: 5.952071180412619e-06
Iteration 2, Residual norm: 8.308163404681622e-12

Iteration 0, Residual norm: 0.014074195758535576
Iteration 1, Residual norm: 6.004781224075708e-06
Iteration 2, Residual norm: 8.024842452062071e-12

Iteration 0, Residual norm: 0.014216795626150148
Iteration 1, Residual norm: 6.0527246209766074e-06
Iteration 2, Residual norm: 7.732555219754151e-12

Iteration 0, Residual norm: 0.014332014284792655
Iteration 1, Residual norm: 6.091800765016758e-06
Iteration 2, Residual norm: 7.430186080547348e-12
Progress: 13%, Time: 1.30e+00, accepted step: True, dt: 2.00e-02, iterations: 2

Iteration 0, Residual norm: 0.0144702196322019
Iteration 1, Residual norm: 6.119048994369564e-06
Iteration 2, Residual norm: 7.119309208984339e-12

Iteration 0, Residual norm: 0.014582564703022122
Iteration 1, Residual norm: 6.133644851821215e-06
Iteration 2, Residual norm: 6.802921875127348e-12

Iteration 0, Residual norm: 0.014716046064072197
Iteration 1, Residual norm: 6.1365767110078145e-06
Iteration 2, Residual norm: 6.485641230972342e-12

Iteration 0, Residual norm: 0.014824877477481386
Iteration 1, Residual norm: 6.129833767609812e-06
Iteration 2, Residual norm: 6.174668469685792e-12

Iteration 0, Residual norm: 0.014953066322762996
Iteration 1, Residual norm: 6.115157766630786e-06
Iteration 2, Residual norm: 5.873046223412052e-12

Iteration 0, Residual norm: 0.01505782975487296
Iteration 1, Residual norm: 6.092887594694906e-06
Iteration 2, Residual norm: 5.58133320321569e-12

Iteration 0, Residual norm: 0.015180555626277454
Iteration 1, Residual norm: 6.062303124127311e-06
Iteration 2, Residual norm: 5.304078261643279e-12

Iteration 0, Residual norm: 0.015281216249177835
Iteration 1, Residual norm: 6.0235507549061864e-06
Iteration 2, Residual norm: 5.046193446477015e-12

Iteration 0, Residual norm: 0.01539862657012962
Iteration 1, Residual norm: 5.978505350357472e-06
Iteration 2, Residual norm: 4.8057841227735134e-12

Iteration 0, Residual norm: 0.0154951685109001
Iteration 1, Residual norm: 5.929378020369113e-06
Iteration 2, Residual norm: 4.580422554703962e-12

Iteration 0, Residual norm: 0.015607289206815506
Iteration 1, Residual norm: 5.877121072286494e-06
Iteration 2, Residual norm: 4.3686304855966e-12

Iteration 0, Residual norm: 0.01569972458973703
Iteration 1, Residual norm: 5.821593720346689e-06
Iteration 2, Residual norm: 4.170062475720763e-12

Iteration 0, Residual norm: 0.01580680010195569
Iteration 1, Residual norm: 5.76262948073344e-06
Iteration 2, Residual norm: 3.984917474251317e-12

Iteration 0, Residual norm: 0.015895369281852628
Iteration 1, Residual norm: 5.7007222037733865e-06
Iteration 2, Residual norm: 3.8136403549830965e-12

Iteration 0, Residual norm: 0.015997638665400275
Iteration 1, Residual norm: 5.636793557017735e-06
Iteration 2, Residual norm: 3.656366580752949e-12

Iteration 0, Residual norm: 0.01608232975678911
Iteration 1, Residual norm: 5.571502173805968e-06
Iteration 2, Residual norm: 3.5131317918182154e-12

Iteration 0, Residual norm: 0.016179684568465733
Iteration 1, Residual norm: 5.50491399012325e-06
Iteration 2, Residual norm: 3.383245531031223e-12

Iteration 0, Residual norm: 0.016260221733429572
Iteration 1, Residual norm: 5.436970939806803e-06
Iteration 2, Residual norm: 3.266330593402433e-12

Iteration 0, Residual norm: 0.01635235195890206
Iteration 1, Residual norm: 5.368118904230099e-06
Iteration 2, Residual norm: 3.1605159558454366e-12

Iteration 0, Residual norm: 0.01642818024733699
Iteration 1, Residual norm: 5.299468034585617e-06
Iteration 2, Residual norm: 3.0629161602258784e-12

Iteration 0, Residual norm: 0.016514318298478152
Iteration 1, Residual norm: 5.2323169619908505e-06
Iteration 2, Residual norm: 2.9706999196783265e-12

Iteration 0, Residual norm: 0.016584381456525507
Iteration 1, Residual norm: 5.167522211077226e-06
Iteration 2, Residual norm: 2.8815401482648137e-12

Iteration 0, Residual norm: 0.016663503311427588
Iteration 1, Residual norm: 5.10501111547951e-06
Iteration 2, Residual norm: 2.792769383462624e-12

Iteration 0, Residual norm: 0.016726999261580004
Iteration 1, Residual norm: 5.04379799271334e-06
Iteration 2, Residual norm: 2.701623502766624e-12

Iteration 0, Residual norm: 0.01679885094085868
Iteration 1, Residual norm: 4.9826554143748095e-06
Iteration 2, Residual norm: 2.6065036925869467e-12
Progress: 18%, Time: 1.80e+00, accepted step: True, dt: 2.00e-02, iterations: 2

Iteration 0, Residual norm: 0.01685598839636612
Iteration 1, Residual norm: 4.921150804033531e-06
Iteration 2, Residual norm: 2.507698866303901e-12

Iteration 0, Residual norm: 0.016921156858846504
Iteration 1, Residual norm: 4.859789441074075e-06
Iteration 2, Residual norm: 2.4058765710584427e-12

Iteration 0, Residual norm: 0.016972584505898114
Iteration 1, Residual norm: 4.799028017878025e-06
Iteration 2, Residual norm: 2.3027778079894867e-12

Iteration 0, Residual norm: 0.01703173330351161
Iteration 1, Residual norm: 4.738948699063895e-06
Iteration 2, Residual norm: 2.2017729381850226e-12

Iteration 0, Residual norm: 0.017078016393583067
Iteration 1, Residual norm: 4.680113629335209e-06
Iteration 2, Residual norm: 2.107025172599876e-12

Iteration 0, Residual norm: 0.017131716480749976
Iteration 1, Residual norm: 4.623377334540193e-06
Iteration 2, Residual norm: 2.021954618222771e-12

Iteration 0, Residual norm: 0.01717343895524096
Iteration 1, Residual norm: 4.568820688593562e-06
Iteration 2, Residual norm: 1.9488624781697994e-12

Iteration 0, Residual norm: 0.017222457783721362
Iteration 1, Residual norm: 4.515953273397569e-06
Iteration 2, Residual norm: 1.8885092025504926e-12

Iteration 0, Residual norm: 0.017260508698196775
Iteration 1, Residual norm: 4.464607282359142e-06
Iteration 2, Residual norm: 1.839628349602606e-12

Iteration 0, Residual norm: 0.017305774001879503
Iteration 1, Residual norm: 4.41461865987398e-06
Iteration 2, Residual norm: 1.7987056620208892e-12

Iteration 0, Residual norm: 0.0173408483830963
Iteration 1, Residual norm: 4.365502608491893e-06
Iteration 2, Residual norm: 1.7611334221560803e-12

Iteration 0, Residual norm: 0.01738269668081021
Iteration 1, Residual norm: 4.317271352067131e-06
Iteration 2, Residual norm: 1.7223135298310179e-12

Iteration 0, Residual norm: 0.017414595071628555
Iteration 1, Residual norm: 4.270650874723819e-06
Iteration 2, Residual norm: 1.678563618437388e-12

Iteration 0, Residual norm: 0.017452318112167003
Iteration 1, Residual norm: 4.226121730371338e-06
Iteration 2, Residual norm: 1.627433422980858e-12

Iteration 0, Residual norm: 0.0174798764738006
Iteration 1, Residual norm: 4.183809665136766e-06
Iteration 2, Residual norm: 1.5683314324910052e-12

Iteration 0, Residual norm: 0.01751209144660161
Iteration 1, Residual norm: 4.143829594303126e-06
Iteration 2, Residual norm: 1.5027520417741028e-12

Iteration 0, Residual norm: 0.017533881184986
Iteration 1, Residual norm: 4.105664434877991e-06
Iteration 2, Residual norm: 1.4342740827010887e-12

Iteration 0, Residual norm: 0.017559422376621375
Iteration 1, Residual norm: 4.06793946935722e-06
Iteration 2, Residual norm: 1.3678660738618183e-12

Iteration 0, Residual norm: 0.017574730829008924
Iteration 1, Residual norm: 4.029523533867059e-06
Iteration 2, Residual norm: 1.309087242740065e-12

Iteration 0, Residual norm: 0.01759347257352215
Iteration 1, Residual norm: 3.990105115614938e-06
Iteration 2, Residual norm: 1.2626565938260789e-12

Iteration 0, Residual norm: 0.017602621019615454
Iteration 1, Residual norm: 3.9498500265200645e-06
Iteration 2, Residual norm: 1.230918340494261e-12

Iteration 0, Residual norm: 0.01761519704544317
Iteration 1, Residual norm: 3.90921714912368e-06
Iteration 2, Residual norm: 1.2124777765343599e-12

Iteration 0, Residual norm: 0.017618934847158197
Iteration 1, Residual norm: 3.8688043817895506e-06
Iteration 2, Residual norm: 1.2023135485836216e-12

Iteration 0, Residual norm: 0.017626175293842153
Iteration 1, Residual norm: 3.828938441826432e-06
Iteration 2, Residual norm: 1.193408561728101e-12

Iteration 0, Residual norm: 0.01762538768437859
Iteration 1, Residual norm: 3.789539853051448e-06
Iteration 2, Residual norm: 1.1790876469912113e-12
Progress: 23%, Time: 2.30e+00, accepted step: True, dt: 2.00e-02, iterations: 2

Iteration 0, Residual norm: 0.017628302280460677
Iteration 1, Residual norm: 3.7502938773440775e-06
Iteration 2, Residual norm: 1.154878348353578e-12

Iteration 0, Residual norm: 0.01762406337944352
Iteration 1, Residual norm: 3.7109277854515315e-06
Iteration 2, Residual norm: 1.1196193872214694e-12

Iteration 0, Residual norm: 0.017623724230459804
Iteration 1, Residual norm: 3.6713935429100587e-06
Iteration 2, Residual norm: 1.0756026658868473e-12

Iteration 0, Residual norm: 0.017616960265243384
Iteration 1, Residual norm: 3.631915089867095e-06
Iteration 2, Residual norm: 1.0280499084099378e-12

Iteration 0, Residual norm: 0.01761409828608765
Iteration 1, Residual norm: 3.592961758195121e-06
Iteration 2, Residual norm: 9.839976898021857e-13

Iteration 0, Residual norm: 0.01760528296557941
Iteration 1, Residual norm: 3.5551942574706895e-06
Iteration 2, Residual norm: 9.506357696799764e-13

Iteration 0, Residual norm: 0.017600139858338963
Iteration 1, Residual norm: 3.519309658128731e-06
Iteration 2, Residual norm: 9.329692688837601e-13

Iteration 0, Residual norm: 0.01758928042577693
Iteration 1, Residual norm: 3.4858620472329426e-06
Iteration 2, Residual norm: 9.321501485418342e-13

Iteration 0, Residual norm: 0.0175816866213454
Iteration 1, Residual norm: 3.4551175985326036e-06
Iteration 2, Residual norm: 9.453660901784062e-13

Iteration 0, Residual norm: 0.01756844069502321
Iteration 1, Residual norm: 3.426995260113752e-06
Iteration 2, Residual norm: 9.666291753957187e-13

Iteration 0, Residual norm: 0.01755792078738339
Iteration 1, Residual norm: 3.4010310943568025e-06
Iteration 2, Residual norm: 9.882576723308703e-13

Iteration 0, Residual norm: 0.017541662503287016
Iteration 1, Residual norm: 3.3763624119592713e-06
Iteration 2, Residual norm: 1.0029721318823203e-12

Iteration 0, Residual norm: 0.017527453720901007
Iteration 1, Residual norm: 3.3518111657182373e-06
Iteration 2, Residual norm: 1.0057252576395454e-12

Iteration 0, Residual norm: 0.017507222876677074
Iteration 1, Residual norm: 3.3261655297264377e-06
Iteration 2, Residual norm: 9.949025891036272e-13

Iteration 0, Residual norm: 0.017488103228787896
Iteration 1, Residual norm: 3.2985099414547925e-06
Iteration 2, Residual norm: 9.727980957580197e-13

Iteration 0, Residual norm: 0.017462250563870012
Iteration 1, Residual norm: 3.2684181401944657e-06
Iteration 2, Residual norm: 9.452517589429219e-13

Iteration 0, Residual norm: 0.017435969507987172
Iteration 1, Residual norm: 3.235969425265437e-06
Iteration 2, Residual norm: 9.199675882352656e-13

Iteration 0, Residual norm: 0.017401394843274375
Iteration 1, Residual norm: 3.201752493560684e-06
Iteration 2, Residual norm: 9.038877950554409e-13

Iteration 0, Residual norm: 0.01736382852355292
Iteration 1, Residual norm: 3.1669201656856997e-06
Iteration 2, Residual norm: 9.010234410944237e-13

Iteration 0, Residual norm: 0.017315333050104513
Iteration 1, Residual norm: 3.133094281740209e-06
Iteration 2, Residual norm: 9.112894725882214e-13

Iteration 0, Residual norm: 0.017260567244834255
Iteration 1, Residual norm: 3.10177624393113e-06
Iteration 2, Residual norm: 9.304602418960033e-13

Iteration 0, Residual norm: 0.017192322406893147
Iteration 1, Residual norm: 3.073411285480623e-06
Iteration 2, Residual norm: 9.52084447649878e-13

Iteration 0, Residual norm: 0.017115797582088
Iteration 1, Residual norm: 3.0469177060366224e-06
Iteration 2, Residual norm: 9.705455857675589e-13

Iteration 0, Residual norm: 0.01702545267857513
Iteration 1, Residual norm: 3.0204052027876576e-06
Iteration 2, Residual norm: 9.823030544015768e-13

Iteration 0, Residual norm: 0.016927159434321732
Iteration 1, Residual norm: 2.9926694010687393e-06
Iteration 2, Residual norm: 9.868350482085837e-13
Progress: 28%, Time: 2.80e+00, accepted step: True, dt: 2.00e-02, iterations: 2

Iteration 0, Residual norm: 0.016815946009930215
Iteration 1, Residual norm: 2.964180474751981e-06
Iteration 2, Residual norm: 9.869868453953798e-13

Iteration 0, Residual norm: 0.016696469062462898
Iteration 1, Residual norm: 2.9367881930421283e-06
Iteration 2, Residual norm: 9.863013742895608e-13

Iteration 0, Residual norm: 0.016562419581627748
Iteration 1, Residual norm: 2.9127742019607173e-06
Iteration 2, Residual norm: 9.875801620825408e-13

Iteration 0, Residual norm: 0.016416220235118885
Iteration 1, Residual norm: 2.8943897374064218e-06
Iteration 2, Residual norm: 9.944354428891489e-13

Iteration 0, Residual norm: 0.01625034601036815
Iteration 1, Residual norm: 2.8840003646336123e-06
Iteration 2, Residual norm: 1.0105133795525218e-12

Iteration 0, Residual norm: 0.016066225544480114
Iteration 1, Residual norm: 2.8838668726125854e-06
Iteration 2, Residual norm: 1.0378998421613039e-12

Iteration 0, Residual norm: 0.01585719507639812
Iteration 1, Residual norm: 2.895168374515539e-06
Iteration 2, Residual norm: 1.0751425385129615e-12

Iteration 0, Residual norm: 0.015626557747124164
Iteration 1, Residual norm: 2.9172554443880675e-06
Iteration 2, Residual norm: 1.117159521311074e-12

Iteration 0, Residual norm: 0.015371757212634035
Iteration 1, Residual norm: 2.9481385132070183e-06
Iteration 2, Residual norm: 1.1583129214965866e-12

Iteration 0, Residual norm: 0.015100870242376795
Iteration 1, Residual norm: 2.985789868434706e-06
Iteration 2, Residual norm: 1.195377762580895e-12

Iteration 0, Residual norm: 0.014816682033249344
Iteration 1, Residual norm: 3.0291206491668915e-06
Iteration 2, Residual norm: 1.2289072438608683e-12

Iteration 0, Residual norm: 0.014530713118974966
Iteration 1, Residual norm: 3.0781739265116127e-06
Iteration 2, Residual norm: 1.2626387671619359e-12

Iteration 0, Residual norm: 0.014247124059358442
Iteration 1, Residual norm: 3.133817406425776e-06
Iteration 2, Residual norm: 1.3042550937816335e-12

Iteration 0, Residual norm: 0.01397554349650886
Iteration 1, Residual norm: 3.197196194299907e-06
Iteration 2, Residual norm: 1.3618323897584458e-12

Iteration 0, Residual norm: 0.013716683089690334
Iteration 1, Residual norm: 3.2690311381601693e-06
Iteration 2, Residual norm: 1.4382591512582353e-12

Iteration 0, Residual norm: 0.013475537495217285
Iteration 1, Residual norm: 3.348973007165095e-06
Iteration 2, Residual norm: 1.5294699974343737e-12

Iteration 0, Residual norm: 0.013249249210883829
Iteration 1, Residual norm: 3.435481843791e-06
Iteration 2, Residual norm: 1.6258427969544135e-12

Iteration 0, Residual norm: 0.013040003636687506
Iteration 1, Residual norm: 3.526709456423416e-06
Iteration 2, Residual norm: 1.7186005050051114e-12

Iteration 0, Residual norm: 0.012844020711395067
Iteration 1, Residual norm: 3.6221442691708465e-06
Iteration 2, Residual norm: 1.8028207794114744e-12

Iteration 0, Residual norm: 0.012662999218641759
Iteration 1, Residual norm: 3.7236078522651386e-06
Iteration 2, Residual norm: 1.8827836935754234e-12

Iteration 0, Residual norm: 0.012493722396608
Iteration 1, Residual norm: 3.834370673280844e-06
Iteration 2, Residual norm: 1.9690961171997534e-12

Iteration 0, Residual norm: 0.012337939247003786
Iteration 1, Residual norm: 3.9572859196703415e-06
Iteration 2, Residual norm: 2.078911165676858e-12

Iteration 0, Residual norm: 0.012192931600363922
Iteration 1, Residual norm: 4.094075012654467e-06
Iteration 2, Residual norm: 2.2299876806222838e-12

Iteration 0, Residual norm: 0.012060176728580497
Iteration 1, Residual norm: 4.245849005159069e-06
Iteration 2, Residual norm: 2.4241606061389698e-12

Iteration 0, Residual norm: 0.01193721075218783
Iteration 1, Residual norm: 4.412734675574865e-06
Iteration 2, Residual norm: 2.6602844384502122e-12
Progress: 33%, Time: 3.30e+00, accepted step: True, dt: 2.00e-02, iterations: 2

Iteration 0, Residual norm: 0.011825134216728093
Iteration 1, Residual norm: 4.592490643698649e-06
Iteration 2, Residual norm: 2.924446696906542e-12

Iteration 0, Residual norm: 0.011721742482435371
Iteration 1, Residual norm: 4.781150829797111e-06
Iteration 2, Residual norm: 3.201410489735488e-12

Iteration 0, Residual norm: 0.011627786243579966
Iteration 1, Residual norm: 4.976366565823727e-06
Iteration 2, Residual norm: 3.469130182071291e-12

Iteration 0, Residual norm: 0.011541447176083899
Iteration 1, Residual norm: 5.1794431775803465e-06
Iteration 2, Residual norm: 3.711529579197743e-12

Iteration 0, Residual norm: 0.011463490863466791
Iteration 1, Residual norm: 5.393696133783385e-06
Iteration 2, Residual norm: 3.949553632085775e-12

Iteration 0, Residual norm: 0.011393102905695595
Iteration 1, Residual norm: 5.6221560534618725e-06
Iteration 2, Residual norm: 4.20331224928758e-12

Iteration 0, Residual norm: 0.01133164875561391
Iteration 1, Residual norm: 5.867728548125711e-06
Iteration 2, Residual norm: 4.526700463734386e-12

Iteration 0, Residual norm: 0.011279435024540983
Iteration 1, Residual norm: 6.134191540038184e-06
Iteration 2, Residual norm: 4.973413259430721e-12

Iteration 0, Residual norm: 0.011237910543142698
Iteration 1, Residual norm: 6.424871286186398e-06
Iteration 2, Residual norm: 5.555651717170105e-12

Iteration 0, Residual norm: 0.011207636821634738
Iteration 1, Residual norm: 6.7393521870504905e-06
Iteration 2, Residual norm: 6.262584857496041e-12

Iteration 0, Residual norm: 0.011189488943621808
Iteration 1, Residual norm: 7.075689588253995e-06
Iteration 2, Residual norm: 7.0347092105813396e-12

Iteration 0, Residual norm: 0.011185366424589071
Iteration 1, Residual norm: 7.426091771363348e-06
Iteration 2, Residual norm: 7.812724077912937e-12

Iteration 0, Residual norm: 0.01119708424950805
Iteration 1, Residual norm: 7.78413163983329e-06
Iteration 2, Residual norm: 8.550481558022462e-12

Iteration 0, Residual norm: 0.011225050419519777
Iteration 1, Residual norm: 8.154382624633783e-06
Iteration 2, Residual norm: 9.233657943404621e-12

Iteration 0, Residual norm: 0.011269467019630924
Iteration 1, Residual norm: 8.54547255811871e-06
Iteration 2, Residual norm: 9.874364884208138e-12

Iteration 0, Residual norm: 0.011329986992973613
Iteration 1, Residual norm: 8.964087425395654e-06
Iteration 2, Residual norm: 1.0558144320776246e-11

Iteration 0, Residual norm: 0.011405688667194252
Iteration 1, Residual norm: 9.413967490781106e-06
Iteration 2, Residual norm: 1.1443557942209775e-11

Iteration 0, Residual norm: 0.011495255947286296
Iteration 1, Residual norm: 9.899197989651361e-06
Iteration 2, Residual norm: 1.2604534381530919e-11

Iteration 0, Residual norm: 0.011597304941988087
Iteration 1, Residual norm: 1.042316096702374e-05
Iteration 2, Residual norm: 1.4136162098184106e-11

Iteration 0, Residual norm: 0.01171062537256893
Iteration 1, Residual norm: 1.0982187761923521e-05
Iteration 2, Residual norm: 1.596468030494911e-11

Iteration 0, Residual norm: 0.011834194446131291
Iteration 1, Residual norm: 1.1565670634384724e-05
Iteration 2, Residual norm: 1.8001255366178885e-11

Iteration 0, Residual norm: 0.011967026235608679
Iteration 1, Residual norm: 1.2165917853842216e-05
Iteration 2, Residual norm: 2.0118344860861208e-11

Iteration 0, Residual norm: 0.01210818652783199
Iteration 1, Residual norm: 1.2783616499771244e-05
Iteration 2, Residual norm: 2.2125732578393836e-11

Iteration 0, Residual norm: 0.01225713701380524
Iteration 1, Residual norm: 1.342259955295757e-05
Iteration 2, Residual norm: 2.4079520053457636e-11

Iteration 0, Residual norm: 0.012414196288248566
Iteration 1, Residual norm: 1.4085789779281099e-05
Iteration 2, Residual norm: 2.5943279489892967e-11
Progress: 38%, Time: 3.80e+00, accepted step: True, dt: 2.00e-02, iterations: 2

Iteration 0, Residual norm: 0.012580660930023177
Iteration 1, Residual norm: 1.4779822622565314e-05
Iteration 2, Residual norm: 2.7948049689866132e-11

Iteration 0, Residual norm: 0.01275844932338966
Iteration 1, Residual norm: 1.551963611116504e-05
Iteration 2, Residual norm: 3.0453568082557246e-11

Iteration 0, Residual norm: 0.012949540582011641
Iteration 1, Residual norm: 1.6323226166713852e-05
Iteration 2, Residual norm: 3.370733305513556e-11

Iteration 0, Residual norm: 0.013155627370846757
Iteration 1, Residual norm: 1.7199737740623287e-05
Iteration 2, Residual norm: 3.7929410792020855e-11

Iteration 0, Residual norm: 0.01337806935909628
Iteration 1, Residual norm: 1.8141305735031687e-05
Iteration 2, Residual norm: 4.2958903784516723e-11

Iteration 0, Residual norm: 0.01361792092583375
Iteration 1, Residual norm: 1.912716740531642e-05
Iteration 2, Residual norm: 4.8545143244966376e-11

Iteration 0, Residual norm: 0.013875752469379839
Iteration 1, Residual norm: 2.0139867696167983e-05
Iteration 2, Residual norm: 5.434048190458787e-11

Iteration 0, Residual norm: 0.014151295061724713
Iteration 1, Residual norm: 2.1178787455710874e-05
Iteration 2, Residual norm: 5.99820859285564e-11

Iteration 0, Residual norm: 0.014443168317675241
Iteration 1, Residual norm: 2.2256357651672796e-05
Iteration 2, Residual norm: 6.531989788435823e-11

Iteration 0, Residual norm: 0.014748993563938948
Iteration 1, Residual norm: 2.338343994697977e-05
Iteration 2, Residual norm: 7.036184451450491e-11

Iteration 0, Residual norm: 0.015066026388357034
Iteration 1, Residual norm: 2.4564822500294995e-05
Iteration 2, Residual norm: 7.579783557599924e-11

Iteration 0, Residual norm: 0.0153922617909474
Iteration 1, Residual norm: 2.5811359085873856e-05
Iteration 2, Residual norm: 8.217520785675647e-11

Iteration 0, Residual norm: 0.01572758750529765
Iteration 1, Residual norm: 2.7145385678378806e-05
Iteration 2, Residual norm: 9.031990376914701e-11

Iteration 0, Residual norm: 0.016074194025131314
Iteration 1, Residual norm: 2.857925855718493e-05
Iteration 2, Residual norm: 1.0081008660330972e-10

Iteration 0, Residual norm: 0.016435510210178594
Iteration 1, Residual norm: 3.009774130425789e-05
Iteration 2, Residual norm: 1.1318871300928597e-10

Iteration 0, Residual norm: 0.01681393193555064
Iteration 1, Residual norm: 3.16804623646427e-05
Iteration 2, Residual norm: 1.272287333803341e-10

Iteration 0, Residual norm: 0.017208889259078295
Iteration 1, Residual norm: 3.332670508604946e-05
Iteration 2, Residual norm: 1.416072395636253e-10

Iteration 0, Residual norm: 0.017616921171861023
Iteration 1, Residual norm: 3.504147697624766e-05
Iteration 2, Residual norm: 1.5616505738026537e-10

Iteration 0, Residual norm: 0.018033785802637547
Iteration 1, Residual norm: 3.681873752379783e-05
Iteration 2, Residual norm: 1.7032338695069503e-10

Iteration 0, Residual norm: 0.01845698322070806
Iteration 1, Residual norm: 3.8656629092457276e-05
Iteration 2, Residual norm: 1.8442469084952966e-10

Iteration 0, Residual norm: 0.018887083614162893
Iteration 1, Residual norm: 4.058030494292954e-05
Iteration 2, Residual norm: 1.999879596525905e-10

Iteration 0, Residual norm: 0.019327504251413368
Iteration 1, Residual norm: 4.2635243240179694e-05
Iteration 2, Residual norm: 2.1842517793894534e-10

Iteration 0, Residual norm: 0.019783301941262783
Iteration 1, Residual norm: 4.485854926445236e-05
Iteration 2, Residual norm: 2.4174581822797643e-10

Iteration 0, Residual norm: 0.02025966965895382
Iteration 1, Residual norm: 4.725613887465512e-05
Iteration 2, Residual norm: 2.701917893805504e-10

Iteration 0, Residual norm: 0.020760508056045932
Iteration 1, Residual norm: 4.979756713506663e-05
Iteration 2, Residual norm: 3.0328607621875963e-10
Progress: 43%, Time: 4.30e+00, accepted step: True, dt: 2.00e-02, iterations: 2

Iteration 0, Residual norm: 0.02128729425477691
Iteration 1, Residual norm: 5.24437205915819e-05
Iteration 2, Residual norm: 3.3965438975168046e-10

Iteration 0, Residual norm: 0.02183850544940066
Iteration 1, Residual norm: 5.518826367820016e-05
Iteration 2, Residual norm: 3.7719899173288103e-10

Iteration 0, Residual norm: 0.022409910399355647
Iteration 1, Residual norm: 5.806134845352005e-05
Iteration 2, Residual norm: 4.150381603227028e-10

Iteration 0, Residual norm: 0.02299579003956119
Iteration 1, Residual norm: 6.109262701548514e-05
Iteration 2, Residual norm: 4.5246421699433806e-10

Iteration 0, Residual norm: 0.023590830093829895
Iteration 1, Residual norm: 6.428634774804226e-05
Iteration 2, Residual norm: 4.91978408750079e-10

Iteration 0, Residual norm: 0.02419238837427635
Iteration 1, Residual norm: 6.76501239080419e-05
Iteration 2, Residual norm: 5.380235072435863e-10

Iteration 0, Residual norm: 0.024802645850082985
Iteration 1, Residual norm: 7.123757232395583e-05
Iteration 2, Residual norm: 5.932115782878317e-10

Iteration 0, Residual norm: 0.02542942963289629
Iteration 1, Residual norm: 7.510425565165533e-05
Iteration 2, Residual norm: 6.652412336894293e-10

Iteration 0, Residual norm: 0.026083905762222304
Iteration 1, Residual norm: 7.92117039554372e-05
Iteration 2, Residual norm: 7.509788447421741e-10

Iteration 0, Residual norm: 0.026774511108568856
Iteration 1, Residual norm: 8.348702710312799e-05
Iteration 2, Residual norm: 8.492260634771324e-10

Iteration 0, Residual norm: 0.027500458116066297
Iteration 1, Residual norm: 8.795487211041185e-05
Iteration 2, Residual norm: 9.542066011969578e-10

Iteration 0, Residual norm: 0.028250858759405358
Iteration 1, Residual norm: 9.266686256495015e-05
Iteration 2, Residual norm: 1.0615958990900919e-09

Iteration 0, Residual norm: 0.029011416689549593
Iteration 1, Residual norm: 9.75906554124078e-05
Iteration 2, Residual norm: 1.1691730682591428e-09

Iteration 0, Residual norm: 0.029773226153906967
Iteration 1, Residual norm: 0.00010265647232539463
Iteration 2, Residual norm: 1.2745565243771753e-09

Iteration 0, Residual norm: 0.030537146941532103
Iteration 1, Residual norm: 0.0001078703658933681
Iteration 2, Residual norm: 1.3849294422110646e-09

Iteration 0, Residual norm: 0.03131207417052306
Iteration 1, Residual norm: 0.00011333882225290065
Iteration 2, Residual norm: 1.5058254465359692e-09

Iteration 0, Residual norm: 0.032109662915588845
Iteration 1, Residual norm: 0.00011917133052033773
Iteration 2, Residual norm: 1.6518058305038131e-09

Iteration 0, Residual norm: 0.032939493561351944
Iteration 1, Residual norm: 0.00012542861535038425
Iteration 2, Residual norm: 1.8269449164051249e-09

Iteration 0, Residual norm: 0.033806638270763245
Iteration 1, Residual norm: 0.00013207983650447312
Iteration 2, Residual norm: 2.0326523324147034e-09

Iteration 0, Residual norm: 0.03471116656062808
Iteration 1, Residual norm: 0.0001390022918969188
Iteration 2, Residual norm: 2.2649602392068796e-09

Iteration 0, Residual norm: 0.03564841876399231
Iteration 1, Residual norm: 0.00014612721537933193
Iteration 2, Residual norm: 2.5110609756500495e-09

Iteration 0, Residual norm: 0.036610032923312434
Iteration 1, Residual norm: 0.00015350821931135535
Iteration 2, Residual norm: 2.764082303793718e-09

Iteration 0, Residual norm: 0.03758678228713414
Iteration 1, Residual norm: 0.00016124029828149875
Iteration 2, Residual norm: 3.01947006613034e-09

Iteration 0, Residual norm: 0.03857160550491331
Iteration 1, Residual norm: 0.00016938593391806488
Iteration 2, Residual norm: 3.2782609206702026e-09

Iteration 0, Residual norm: 0.039561459748635994
Iteration 1, Residual norm: 0.0001779651051213299
Iteration 2, Residual norm: 3.572834057949781e-09
Progress: 48%, Time: 4.80e+00, accepted step: True, dt: 2.00e-02, iterations: 2

Iteration 0, Residual norm: 0.040558934891965825
Iteration 1, Residual norm: 0.00018704940602565727
Iteration 2, Residual norm: 3.912740139698075e-09

Iteration 0, Residual norm: 0.04157430862366844
Iteration 1, Residual norm: 0.00019682937266447148
Iteration 2, Residual norm: 4.341928620052009e-09

Iteration 0, Residual norm: 0.04262685555998007
Iteration 1, Residual norm: 0.000207285615639636
Iteration 2, Residual norm: 4.879362172007381e-09

Iteration 0, Residual norm: 0.043738120822728224
Iteration 1, Residual norm: 0.00021808413882364332
Iteration 2, Residual norm: 5.497603638199546e-09

Iteration 0, Residual norm: 0.04491452828480908
Iteration 1, Residual norm: 0.000229293037109641
Iteration 2, Residual norm: 6.176204342710208e-09

Iteration 0, Residual norm: 0.046135957418458155
Iteration 1, Residual norm: 0.00024126123012458214
Iteration 2, Residual norm: 6.8851308595718045e-09

Iteration 0, Residual norm: 0.047368679810316616
Iteration 1, Residual norm: 0.0002539022716415044
Iteration 2, Residual norm: 7.619118096149719e-09

Iteration 0, Residual norm: 0.048590741614408045
Iteration 1, Residual norm: 0.0002669270628117447
Iteration 2, Residual norm: 8.31992124694394e-09

Iteration 0, Residual norm: 0.04980402731173367
Iteration 1, Residual norm: 0.00028017777833182087
Iteration 2, Residual norm: 9.069525039430552e-09

Iteration 0, Residual norm: 0.05102925789818551
Iteration 1, Residual norm: 0.0002939142917779828
Iteration 2, Residual norm: 9.832959752619917e-09

Iteration 0, Residual norm: 0.05229051243016246
Iteration 1, Residual norm: 0.00030849851890354664
Iteration 2, Residual norm: 1.0713815067865509e-08
Iteration 3, Residual norm: 4.3289891400389517e-16

Iteration 0, Residual norm: 0.053602631689226216
Iteration 1, Residual norm: 0.0003240568665624462
Iteration 2, Residual norm: 1.1779641054054977e-08
Iteration 3, Residual norm: 4.2476427096642545e-16

Iteration 0, Residual norm: 0.05496970999210133
Iteration 1, Residual norm: 0.00034067095959037174
Iteration 2, Residual norm: 1.2975665769964e-08
Iteration 3, Residual norm: 4.334999943870725e-16

Iteration 0, Residual norm: 0.05638862006662401
Iteration 1, Residual norm: 0.00035805008711967405
Iteration 2, Residual norm: 1.4384412884463663e-08
Iteration 3, Residual norm: 4.383564422457012e-16

Iteration 0, Residual norm: 0.057852668916840135
Iteration 1, Residual norm: 0.00037581415151641707
Iteration 2, Residual norm: 1.5876176330655374e-08
Iteration 3, Residual norm: 4.3733538184654023e-16

Iteration 0, Residual norm: 0.05934906670099847
Iteration 1, Residual norm: 0.00039410023415725374
Iteration 2, Residual norm: 1.742090837868644e-08
Iteration 3, Residual norm: 4.3305024968746757e-16

Iteration 0, Residual norm: 0.060863507099493085
Iteration 1, Residual norm: 0.00041310939805114905
Iteration 2, Residual norm: 1.902241035124936e-08
Iteration 3, Residual norm: 4.3888789391952715e-16

Iteration 0, Residual norm: 0.062389419493451795
Iteration 1, Residual norm: 0.0004330317639994646
Iteration 2, Residual norm: 2.0617692365663452e-08
Iteration 3, Residual norm: 4.3692477418690094e-16

Iteration 0, Residual norm: 0.06392582100242201
Iteration 1, Residual norm: 0.00045416675883442763
Iteration 2, Residual norm: 2.239392845629146e-08
Iteration 3, Residual norm: 4.397549705570144e-16

Iteration 0, Residual norm: 0.06547574602027195
Iteration 1, Residual norm: 0.0004765253519858039
Iteration 2, Residual norm: 2.4466546300836272e-08
Iteration 3, Residual norm: 4.4136806699524327e-16

Iteration 0, Residual norm: 0.06704314365766355
Iteration 1, Residual norm: 0.0005004287464724877
Iteration 2, Residual norm: 2.6960805023652978e-08
Iteration 3, Residual norm: 4.4233395647604623e-16

Iteration 0, Residual norm: 0.0686431144521709
Iteration 1, Residual norm: 0.0005264146833190612
Iteration 2, Residual norm: 3.014402509770777e-08
Iteration 3, Residual norm: 4.52423187641808e-16

Iteration 0, Residual norm: 0.07031857987128293
Iteration 1, Residual norm: 0.0005532169906283604
Iteration 2, Residual norm: 3.393209529807678e-08
Iteration 3, Residual norm: 4.641432119793577e-16

Iteration 0, Residual norm: 0.07211173561454129
Iteration 1, Residual norm: 0.0005800664018968385
Iteration 2, Residual norm: 3.813569677444444e-08
Iteration 3, Residual norm: 4.836306693401289e-16

Iteration 0, Residual norm: 0.07400563269602009
Iteration 1, Residual norm: 0.00060913330252568
Iteration 2, Residual norm: 4.2269436238431555e-08
Iteration 3, Residual norm: 5.084798914969919e-16
Progress: 53%, Time: 5.30e+00, accepted step: True, dt: 2.00e-02, iterations: 3

Iteration 0, Residual norm: 0.07592481972362394
Iteration 1, Residual norm: 0.0006403326264158266
Iteration 2, Residual norm: 4.70779520932972e-08
Iteration 3, Residual norm: 5.461594671912185e-16

Iteration 0, Residual norm: 0.07781391672707046
Iteration 1, Residual norm: 0.0006723877947819507
Iteration 2, Residual norm: 5.118812088906975e-08
Iteration 3, Residual norm: 5.926397997347642e-16

Iteration 0, Residual norm: 0.07966956151503315
Iteration 1, Residual norm: 0.0007049841920479372
Iteration 2, Residual norm: 5.574521830765021e-08
Iteration 3, Residual norm: 6.295785732981425e-16

Iteration 0, Residual norm: 0.08152802965223892
Iteration 1, Residual norm: 0.0007379528915187003
Iteration 2, Residual norm: 6.064378872133386e-08
Iteration 3, Residual norm: 6.613139172468367e-16

Iteration 0, Residual norm: 0.08344017101909151
Iteration 1, Residual norm: 0.0007730911915680547
Iteration 2, Residual norm: 6.525779854362631e-08
Iteration 3, Residual norm: 6.918985640852047e-16

Iteration 0, Residual norm: 0.08542882458753265
Iteration 1, Residual norm: 0.0008106084671508423
Iteration 2, Residual norm: 7.189075595231563e-08
Iteration 3, Residual norm: 7.514764954497104e-16

Iteration 0, Residual norm: 0.08749384626662526
Iteration 1, Residual norm: 0.0008501997283938865
Iteration 2, Residual norm: 7.843778521820882e-08
Iteration 3, Residual norm: 8.731652911325792e-16

Iteration 0, Residual norm: 0.08961978911617884
Iteration 1, Residual norm: 0.0008924078429546365
Iteration 2, Residual norm: 8.64114448273145e-08
Iteration 3, Residual norm: 1.0229078254085801e-15

Iteration 0, Residual norm: 0.09180483411829132
Iteration 1, Residual norm: 0.0009350048286497834
Iteration 2, Residual norm: 9.517657852203972e-08
Iteration 3, Residual norm: 1.273942919217698e-15

Iteration 0, Residual norm: 0.0940424092188872
Iteration 1, Residual norm: 0.0009782562260071566
Iteration 2, Residual norm: 1.0354484711223436e-07
Iteration 3, Residual norm: 1.5692910541231721e-15

Iteration 0, Residual norm: 0.09629738444514979
Iteration 1, Residual norm: 0.0010235506070195481
Iteration 2, Residual norm: 1.1285786810725556e-07
Iteration 3, Residual norm: 1.814790072399762e-15

Iteration 0, Residual norm: 0.09856582044799794
Iteration 1, Residual norm: 0.0010700738277997928
Iteration 2, Residual norm: 1.217976597161618e-07
Iteration 3, Residual norm: 2.0676868499219904e-15

Iteration 0, Residual norm: 0.10085017239804431
Iteration 1, Residual norm: 0.001119752497155214
Iteration 2, Residual norm: 1.3184296606757966e-07
Iteration 3, Residual norm: 2.330995793989161e-15

Iteration 0, Residual norm: 0.10316015544303471
Iteration 1, Residual norm: 0.0011729017937792696
Iteration 2, Residual norm: 1.437306062774068e-07
Iteration 3, Residual norm: 2.5756658725040175e-15

Iteration 0, Residual norm: 0.10550474882300064
Iteration 1, Residual norm: 0.0012290722393878467
Iteration 2, Residual norm: 1.5778065164847303e-07
Iteration 3, Residual norm: 3.02017150717256e-15

Iteration 0, Residual norm: 0.10784872300913424
Iteration 1, Residual norm: 0.0012902117185359912
Iteration 2, Residual norm: 1.7576650961435557e-07
Iteration 3, Residual norm: 3.7528942882219384e-15

Iteration 0, Residual norm: 0.11021892762528943
Iteration 1, Residual norm: 0.0013562796660354828
Iteration 2, Residual norm: 1.9580801702454856e-07
Iteration 3, Residual norm: 5.129795261007401e-15

Iteration 0, Residual norm: 0.1127413678455037
Iteration 1, Residual norm: 0.0014182458652435125
Iteration 2, Residual norm: 2.2173851579231492e-07
Iteration 3, Residual norm: 6.843347570056464e-15

Iteration 0, Residual norm: 0.11549966889927814
Iteration 1, Residual norm: 0.001482042404586796
Iteration 2, Residual norm: 2.4153451539861957e-07
Iteration 3, Residual norm: 8.8428318028817e-15

Iteration 0, Residual norm: 0.11835002810693804
Iteration 1, Residual norm: 0.0015561872045532247
Iteration 2, Residual norm: 2.6781068961287765e-07
Iteration 3, Residual norm: 1.090610304618478e-14

Iteration 0, Residual norm: 0.12113517075838888
Iteration 1, Residual norm: 0.0016291412658276446
Iteration 2, Residual norm: 2.9132275242944837e-07
Iteration 3, Residual norm: 1.2788908304938477e-14

Iteration 0, Residual norm: 0.1238434984856812
Iteration 1, Residual norm: 0.0017045790855059694
Iteration 2, Residual norm: 3.096184515683209e-07
Iteration 3, Residual norm: 1.4246897707094935e-14

Iteration 0, Residual norm: 0.12649412069405486
Iteration 1, Residual norm: 0.001780509240707367
Iteration 2, Residual norm: 3.4101893276611823e-07
Iteration 3, Residual norm: 1.577544081907163e-14

Iteration 0, Residual norm: 0.12919845177668096
Iteration 1, Residual norm: 0.0018578496922507985
Iteration 2, Residual norm: 3.5997320584035896e-07
Iteration 3, Residual norm: 1.664949766072625e-14

Iteration 0, Residual norm: 0.13201320829806568
Iteration 1, Residual norm: 0.0019447447366512532
Iteration 2, Residual norm: 3.9578034792025683e-07
Iteration 3, Residual norm: 1.8761644913770783e-14
Progress: 58%, Time: 5.80e+00, accepted step: True, dt: 2.00e-02, iterations: 3

Iteration 0, Residual norm: 0.1349556980043562
Iteration 1, Residual norm: 0.0020322110205930883
Iteration 2, Residual norm: 4.3123097811859507e-07
Iteration 3, Residual norm: 2.2435480051625948e-14

Iteration 0, Residual norm: 0.13794714124641805
Iteration 1, Residual norm: 0.0021266705918535924
Iteration 2, Residual norm: 4.6786048561746666e-07
Iteration 3, Residual norm: 2.8123948258306868e-14

Iteration 0, Residual norm: 0.14095948612875142
Iteration 1, Residual norm: 0.002223196949148833
Iteration 2, Residual norm: 5.1517578042688e-07
Iteration 3, Residual norm: 3.510912474435011e-14

Iteration 0, Residual norm: 0.14407939737168393
Iteration 1, Residual norm: 0.002315770172292853
Iteration 2, Residual norm: 5.524942577791033e-07
Iteration 3, Residual norm: 4.241440400455387e-14

Iteration 0, Residual norm: 0.14722478874509645
Iteration 1, Residual norm: 0.0024150399650592564
Iteration 2, Residual norm: 5.971079781507471e-07
Iteration 3, Residual norm: 4.916828908700894e-14

Iteration 0, Residual norm: 0.1503568544586179
Iteration 1, Residual norm: 0.0025162949876921213
Iteration 2, Residual norm: 6.396845812477558e-07
Iteration 3, Residual norm: 5.584134311915774e-14

Iteration 0, Residual norm: 0.15351777619339954
Iteration 1, Residual norm: 0.002624138134015225
Iteration 2, Residual norm: 6.861296840321782e-07
Iteration 3, Residual norm: 6.100131103174626e-14

Iteration 0, Residual norm: 0.1566792051245617
Iteration 1, Residual norm: 0.00273945741920986
Iteration 2, Residual norm: 7.484365365662368e-07
Iteration 3, Residual norm: 6.554641335266859e-14

Iteration 0, Residual norm: 0.15992084018529856
Iteration 1, Residual norm: 0.0028635329401617035
Iteration 2, Residual norm: 8.135620620970467e-07
Iteration 3, Residual norm: 7.758320391318119e-14

Iteration 0, Residual norm: 0.16322357113493383
Iteration 1, Residual norm: 0.002991431522009201
Iteration 2, Residual norm: 9.081445472792433e-07
Iteration 3, Residual norm: 9.888933523944899e-14

Iteration 0, Residual norm: 0.1664016533990776
Iteration 1, Residual norm: 0.0031331520541291995
Iteration 2, Residual norm: 9.86126548382084e-07
Iteration 3, Residual norm: 1.283278628397747e-13

Iteration 0, Residual norm: 0.1695316829499612
Iteration 1, Residual norm: 0.0032770784218021745
Iteration 2, Residual norm: 1.1081486459847217e-06
Iteration 3, Residual norm: 1.7229166823837112e-13

Iteration 0, Residual norm: 0.17305351171418712
Iteration 1, Residual norm: 0.003390379808675136
Iteration 2, Residual norm: 1.1893581836268737e-06
Iteration 3, Residual norm: 1.9970981106178172e-13

Iteration 0, Residual norm: 0.17693494097125753
Iteration 1, Residual norm: 0.0035361203448379154
Iteration 2, Residual norm: 1.2676236065359578e-06
Iteration 3, Residual norm: 2.408042989065267e-13

Iteration 0, Residual norm: 0.1807316548469977
Iteration 1, Residual norm: 0.003687923821771033
Iteration 2, Residual norm: 1.3846815900919535e-06
Iteration 3, Residual norm: 2.6303722493813884e-13

Iteration 0, Residual norm: 0.18439656425761597
Iteration 1, Residual norm: 0.0038320846473342287
Iteration 2, Residual norm: 1.4121042839702995e-06
Iteration 3, Residual norm: 2.716704348933714e-13

Iteration 0, Residual norm: 0.18786921658131606
Iteration 1, Residual norm: 0.003989185013250466
Iteration 2, Residual norm: 1.5513119462948992e-06
Iteration 3, Residual norm: 2.9667754286932296e-13

Iteration 0, Residual norm: 0.1913320784201204
Iteration 1, Residual norm: 0.004140265595940841
Iteration 2, Residual norm: 1.6274581136472307e-06
Iteration 3, Residual norm: 3.003824903861435e-13

Iteration 0, Residual norm: 0.19486212314762527
Iteration 1, Residual norm: 0.004310668391858129
Iteration 2, Residual norm: 1.759143251145653e-06
Iteration 3, Residual norm: 3.535198530960069e-13

Iteration 0, Residual norm: 0.19864862937258587
Iteration 1, Residual norm: 0.004478201430099936
Iteration 2, Residual norm: 1.9186637010323565e-06
Iteration 3, Residual norm: 4.354620777533477e-13

Iteration 0, Residual norm: 0.2025592635599618
Iteration 1, Residual norm: 0.004650654128683909
Iteration 2, Residual norm: 2.048057221256353e-06
Iteration 3, Residual norm: 5.396009212797871e-13

Iteration 0, Residual norm: 0.20631404317478239
Iteration 1, Residual norm: 0.0048292842442083156
Iteration 2, Residual norm: 2.2121641417259284e-06
Iteration 3, Residual norm: 6.326152560962114e-13

Iteration 0, Residual norm: 0.2100332249890166
Iteration 1, Residual norm: 0.004996695930355576
Iteration 2, Residual norm: 2.33216687375628e-06
Iteration 3, Residual norm: 7.225723335605405e-13

Iteration 0, Residual norm: 0.21400054720971554
Iteration 1, Residual norm: 0.00516584672243713
Iteration 2, Residual norm: 2.4665632182694823e-06
Iteration 3, Residual norm: 8.051596632101657e-13

Iteration 0, Residual norm: 0.21786686996454535
Iteration 1, Residual norm: 0.005352105216974939
Iteration 2, Residual norm: 2.616029007527971e-06
Iteration 3, Residual norm: 8.371131803960541e-13
Progress: 63%, Time: 6.30e+00, accepted step: True, dt: 2.00e-02, iterations: 3

Iteration 0, Residual norm: 0.22165151028818741
Iteration 1, Residual norm: 0.005547368573975562
Iteration 2, Residual norm: 2.762946468906511e-06
Iteration 3, Residual norm: 8.518086118935559e-13

Iteration 0, Residual norm: 0.22559084123251208
Iteration 1, Residual norm: 0.005762412350153078
Iteration 2, Residual norm: 3.0039154404329086e-06
Iteration 3, Residual norm: 9.48484835802586e-13

Iteration 0, Residual norm: 0.22936878975642758
Iteration 1, Residual norm: 0.005978110372150826
Iteration 2, Residual norm: 3.229453362246148e-06
Iteration 3, Residual norm: 1.189775792176669e-12

Iteration 0, Residual norm: 0.2333618575599783
Iteration 1, Residual norm: 0.006205371404988143
Iteration 2, Residual norm: 3.6037005814208456e-06
Iteration 3, Residual norm: 1.5807539493199848e-12

Iteration 0, Residual norm: 0.23751972102213997
Iteration 1, Residual norm: 0.006404103447340503
Iteration 2, Residual norm: 3.7732719177556123e-06
Iteration 3, Residual norm: 1.7960600090234709e-12

Iteration 0, Residual norm: 0.24105251087072177
Iteration 1, Residual norm: 0.006649534731873857
Iteration 2, Residual norm: 4.051242031810062e-06
Iteration 3, Residual norm: 2.2868453414780524e-12

Iteration 0, Residual norm: 0.24461025343360707
Iteration 1, Residual norm: 0.006834773996769487
Iteration 2, Residual norm: 4.272660825432189e-06
Iteration 3, Residual norm: 2.3725093663997e-12

Iteration 0, Residual norm: 0.24889306441600292
Iteration 1, Residual norm: 0.007011834817402481
Iteration 2, Residual norm: 4.293689858741998e-06
Iteration 3, Residual norm: 2.4966022717941874e-12

Iteration 0, Residual norm: 0.25343235224669247
Iteration 1, Residual norm: 0.007237248773539902
Iteration 2, Residual norm: 4.587072878107576e-06
Iteration 3, Residual norm: 2.4580535618669233e-12

Iteration 0, Residual norm: 0.25778880138845583
Iteration 1, Residual norm: 0.007471688783668808
Iteration 2, Residual norm: 4.636323035598774e-06
Iteration 3, Residual norm: 2.343177529277834e-12

Iteration 0, Residual norm: 0.26173779332232977
Iteration 1, Residual norm: 0.007713982586737289
Iteration 2, Residual norm: 4.967125675178639e-06
Iteration 3, Residual norm: 2.5555088328002772e-12

Iteration 0, Residual norm: 0.26547694634872476
Iteration 1, Residual norm: 0.007951622373198783
Iteration 2, Residual norm: 5.277345869754092e-06
Iteration 3, Residual norm: 2.9648661658591383e-12

Iteration 0, Residual norm: 0.26932658032385426
Iteration 1, Residual norm: 0.00818909249547238
Iteration 2, Residual norm: 5.6683139551170035e-06
Iteration 3, Residual norm: 3.757699834512841e-12

Iteration 0, Residual norm: 0.2732910704068165
Iteration 1, Residual norm: 0.008429542024833885
Iteration 2, Residual norm: 6.0458294887691415e-06
Iteration 3, Residual norm: 4.509215770369642e-12

Iteration 0, Residual norm: 0.2773806792888133
Iteration 1, Residual norm: 0.008624214916316847
Iteration 2, Residual norm: 6.329901199619995e-06
Iteration 3, Residual norm: 5.0862793624062726e-12

Iteration 0, Residual norm: 0.28162609578404446
Iteration 1, Residual norm: 0.008839847593282959
Iteration 2, Residual norm: 6.521042173416841e-06
Iteration 3, Residual norm: 5.439140266549999e-12

Iteration 0, Residual norm: 0.2852893793641211
Iteration 1, Residual norm: 0.009047328289603246
Iteration 2, Residual norm: 6.764730104979649e-06
Iteration 3, Residual norm: 5.541363050981962e-12

Iteration 0, Residual norm: 0.28901388514478543
Iteration 1, Residual norm: 0.00927621368909049
Iteration 2, Residual norm: 6.9649474203487434e-06
Iteration 3, Residual norm: 5.4011339396888494e-12

Iteration 0, Residual norm: 0.293198592057713
Iteration 1, Residual norm: 0.009510976801662735
Iteration 2, Residual norm: 7.308498618263931e-06
Iteration 3, Residual norm: 5.215262367477034e-12

Iteration 0, Residual norm: 0.2969328392976585
Iteration 1, Residual norm: 0.009786859840348008
Iteration 2, Residual norm: 7.648421498446916e-06
Iteration 3, Residual norm: 6.135882751708958e-12

Iteration 0, Residual norm: 0.30056489139043546
Iteration 1, Residual norm: 0.01010081522046462
Iteration 2, Residual norm: 8.523554744289794e-06
Iteration 3, Residual norm: 7.877338598661376e-12

Iteration 0, Residual norm: 0.3045885478657812
Iteration 1, Residual norm: 0.010313329222886152
Iteration 2, Residual norm: 8.825077242285172e-06
Iteration 3, Residual norm: 9.405110253737699e-12

Iteration 0, Residual norm: 0.3079913402517118
Iteration 1, Residual norm: 0.010572478789730963
Iteration 2, Residual norm: 9.514999165672975e-06
Iteration 3, Residual norm: 1.1698974587375116e-11

Iteration 0, Residual norm: 0.31189540691522355
Iteration 1, Residual norm: 0.010703720216669087
Iteration 2, Residual norm: 9.627803756179659e-06
Iteration 3, Residual norm: 1.1478385893348331e-11

Iteration 0, Residual norm: 0.3158488009947228
Iteration 1, Residual norm: 0.010877680868763902
Iteration 2, Residual norm: 9.726335221005283e-06
Iteration 3, Residual norm: 1.2234368752811363e-11
Progress: 68%, Time: 6.80e+00, accepted step: True, dt: 2.00e-02, iterations: 3

Iteration 0, Residual norm: 0.31856722250681463
Iteration 1, Residual norm: 0.011068673517697513
Iteration 2, Residual norm: 9.73352778274683e-06
Iteration 3, Residual norm: 1.081421627866839e-11

Iteration 0, Residual norm: 0.3218703713828704
Iteration 1, Residual norm: 0.01124458727114867
Iteration 2, Residual norm: 9.653447103933216e-06
Iteration 3, Residual norm: 1.0275826720270932e-11

Iteration 0, Residual norm: 0.32610738557912006
Iteration 1, Residual norm: 0.01149166238360537
Iteration 2, Residual norm: 1.0051144555466854e-05
Iteration 3, Residual norm: 9.901114765114911e-12

Iteration 0, Residual norm: 0.3300357326684362
Iteration 1, Residual norm: 0.011757696842316027
Iteration 2, Residual norm: 1.0634491321770542e-05
Iteration 3, Residual norm: 1.12212947458934e-11

Iteration 0, Residual norm: 0.3336317586285278
Iteration 1, Residual norm: 0.012049783783848238
Iteration 2, Residual norm: 1.1407264764880565e-05
Iteration 3, Residual norm: 1.4723585996560018e-11

Iteration 0, Residual norm: 0.3363667344349594
Iteration 1, Residual norm: 0.012274585371849127
Iteration 2, Residual norm: 1.2113550752219698e-05
Iteration 3, Residual norm: 1.8043848459495675e-11

Iteration 0, Residual norm: 0.33942228485374176
Iteration 1, Residual norm: 0.012452850609492826
Iteration 2, Residual norm: 1.274697494467838e-05
Iteration 3, Residual norm: 2.050325121757295e-11

Iteration 0, Residual norm: 0.3428372777709831
Iteration 1, Residual norm: 0.012591075306553077
Iteration 2, Residual norm: 1.2939463393970722e-05
Iteration 3, Residual norm: 2.06822733492599e-11

Iteration 0, Residual norm: 0.3457116736905306
Iteration 1, Residual norm: 0.012697002209522905
Iteration 2, Residual norm: 1.2715917181859286e-05
Iteration 3, Residual norm: 2.0197705916731133e-11

Iteration 0, Residual norm: 0.34882596605054494
Iteration 1, Residual norm: 0.012854335284155714
Iteration 2, Residual norm: 1.2815148388186778e-05
Iteration 3, Residual norm: 1.8152228852992073e-11

Iteration 0, Residual norm: 0.3521012387789246
Iteration 1, Residual norm: 0.013055698852099236
Iteration 2, Residual norm: 1.3145902366786079e-05
Iteration 3, Residual norm: 1.6867995758609745e-11

Iteration 0, Residual norm: 0.35473178801437694
Iteration 1, Residual norm: 0.013298999128994191
Iteration 2, Residual norm: 1.3662241965809047e-05
Iteration 3, Residual norm: 1.839620539140458e-11

Iteration 0, Residual norm: 0.3574923132619936
Iteration 1, Residual norm: 0.013591746088471966
Iteration 2, Residual norm: 1.4886274203637111e-05
Iteration 3, Residual norm: 2.3107378724075755e-11

Iteration 0, Residual norm: 0.36077237975695703
Iteration 1, Residual norm: 0.013831735293655842
Iteration 2, Residual norm: 1.5385490778971986e-05
Iteration 3, Residual norm: 2.9655920339158e-11

Iteration 0, Residual norm: 0.362720978307318
Iteration 1, Residual norm: 0.014050892144507795
Iteration 2, Residual norm: 1.674316127821887e-05
Iteration 3, Residual norm: 3.4765315656848825e-11

Iteration 0, Residual norm: 0.365127430540182
Iteration 1, Residual norm: 0.014118537572481564
Iteration 2, Residual norm: 1.6588038474196882e-05
Iteration 3, Residual norm: 3.62930486784337e-11

Iteration 0, Residual norm: 0.36802821711245315
Iteration 1, Residual norm: 0.014151761307619023
Iteration 2, Residual norm: 1.672344663235211e-05
Iteration 3, Residual norm: 3.56762008014248e-11

Iteration 0, Residual norm: 0.3698796104258863
Iteration 1, Residual norm: 0.014141613105648175
Iteration 2, Residual norm: 1.5920953212357726e-05
Iteration 3, Residual norm: 2.996117111876189e-11

Iteration 0, Residual norm: 0.3724577124647899
Iteration 1, Residual norm: 0.014095649876895761
Iteration 2, Residual norm: 1.5648513892205163e-05
Iteration 3, Residual norm: 2.656791203851484e-11

Iteration 0, Residual norm: 0.3746770208936881
Iteration 1, Residual norm: 0.014298620826131186
Iteration 2, Residual norm: 1.5620056431416326e-05
Iteration 3, Residual norm: 2.3659117336665372e-11

Iteration 0, Residual norm: 0.37613521630207813
Iteration 1, Residual norm: 0.01449857461646841
Iteration 2, Residual norm: 1.6422764334789697e-05
Iteration 3, Residual norm: 2.7430078956740215e-11

Iteration 0, Residual norm: 0.3788694548099462
Iteration 1, Residual norm: 0.014792900874341666
Iteration 2, Residual norm: 1.785423004155051e-05
Iteration 3, Residual norm: 3.562793263520637e-11

Iteration 0, Residual norm: 0.3815656792553009
Iteration 1, Residual norm: 0.015024804136537787
Iteration 2, Residual norm: 1.915921475091844e-05
Iteration 3, Residual norm: 4.612301646031793e-11

Iteration 0, Residual norm: 0.38349095709443615
Iteration 1, Residual norm: 0.015158211891688558
Iteration 2, Residual norm: 1.9693887354938544e-05
Iteration 3, Residual norm: 5.127598425995025e-11

Iteration 0, Residual norm: 0.3846366121661186
Iteration 1, Residual norm: 0.015215509887329828
Iteration 2, Residual norm: 2.0036732501276102e-05
Iteration 3, Residual norm: 5.293479953053061e-11
Progress: 73%, Time: 7.30e+00, accepted step: True, dt: 2.00e-02, iterations: 3

Iteration 0, Residual norm: 0.38566994704689206
Iteration 1, Residual norm: 0.015211738328468899
Iteration 2, Residual norm: 1.982215902306665e-05
Iteration 3, Residual norm: 4.8205831732904936e-11

Iteration 0, Residual norm: 0.38753719432006456
Iteration 1, Residual norm: 0.015148696443028076
Iteration 2, Residual norm: 1.8770293530007905e-05
Iteration 3, Residual norm: 4.1937457365394394e-11

Iteration 0, Residual norm: 0.38891076900711596
Iteration 1, Residual norm: 0.01522382179377092
Iteration 2, Residual norm: 1.8298148697908346e-05
Iteration 3, Residual norm: 3.5330347207032123e-11

Iteration 0, Residual norm: 0.39045416904546304
Iteration 1, Residual norm: 0.015427569718121991
Iteration 2, Residual norm: 1.936542772217189e-05
Iteration 3, Residual norm: 3.68860480986284e-11

Iteration 0, Residual norm: 0.3925158272116454
Iteration 1, Residual norm: 0.01572210222257112
Iteration 2, Residual norm: 2.055496235047993e-05
Iteration 3, Residual norm: 4.491950835336485e-11

Iteration 0, Residual norm: 0.39429520752093733
Iteration 1, Residual norm: 0.01596496232969542
Iteration 2, Residual norm: 2.2112941028408166e-05
Iteration 3, Residual norm: 6.269517448705137e-11

Iteration 0, Residual norm: 0.39562167916606966
Iteration 1, Residual norm: 0.01617271221389305
Iteration 2, Residual norm: 2.3658716131436557e-05
Iteration 3, Residual norm: 7.188238331158199e-11

Iteration 0, Residual norm: 0.397191179806052
Iteration 1, Residual norm: 0.016217889696848307
Iteration 2, Residual norm: 2.3652660839255513e-05
Iteration 3, Residual norm: 7.924259648584783e-11

Iteration 0, Residual norm: 0.3981650766548322
Iteration 1, Residual norm: 0.016129520506515685
Iteration 2, Residual norm: 2.326814340252073e-05
Iteration 3, Residual norm: 7.090665742183671e-11

Iteration 0, Residual norm: 0.3984467025848995
Iteration 1, Residual norm: 0.015972287130225343
Iteration 2, Residual norm: 2.1908398649503813e-05
Iteration 3, Residual norm: 6.385437654045895e-11

Iteration 0, Residual norm: 0.3997833956886327
Iteration 1, Residual norm: 0.01581942431160715
Iteration 2, Residual norm: 2.092550540791034e-05
Iteration 3, Residual norm: 5.191744607217354e-11

Iteration 0, Residual norm: 0.40059953240816076
Iteration 1, Residual norm: 0.015820169595062567
Iteration 2, Residual norm: 2.0404818438456462e-05
Iteration 3, Residual norm: 4.269703637749672e-11

Iteration 0, Residual norm: 0.40131293903087295
Iteration 1, Residual norm: 0.015893979353619178
Iteration 2, Residual norm: 2.11252355333461e-05
Iteration 3, Residual norm: 4.587018719075715e-11

Iteration 0, Residual norm: 0.4021992691833975
Iteration 1, Residual norm: 0.01616307944013734
Iteration 2, Residual norm: 2.2894369037718808e-05
Iteration 3, Residual norm: 6.283862042877661e-11

Iteration 0, Residual norm: 0.4032218592820377
Iteration 1, Residual norm: 0.016449712648825176
Iteration 2, Residual norm: 2.475944949248552e-05
Iteration 3, Residual norm: 8.146118598116732e-11

Iteration 0, Residual norm: 0.40416080472043864
Iteration 1, Residual norm: 0.016617561693799354
Iteration 2, Residual norm: 2.662665119825395e-05
Iteration 3, Residual norm: 9.489483852757098e-11

Iteration 0, Residual norm: 0.40520865919914545
Iteration 1, Residual norm: 0.01670402455408219
Iteration 2, Residual norm: 2.645348340136628e-05
Iteration 3, Residual norm: 1.0048593024738737e-10

Iteration 0, Residual norm: 0.4057363603469226
Iteration 1, Residual norm: 0.01649660629606952
Iteration 2, Residual norm: 2.5131798626485168e-05
Iteration 3, Residual norm: 8.963193672579445e-11

Iteration 0, Residual norm: 0.4057680025453868
Iteration 1, Residual norm: 0.016417422158540553
Iteration 2, Residual norm: 2.4025111191897176e-05
Iteration 3, Residual norm: 7.798085893811094e-11

Iteration 0, Residual norm: 0.40572363819531343
Iteration 1, Residual norm: 0.01631177116701185
Iteration 2, Residual norm: 2.3312178366783427e-05
Iteration 3, Residual norm: 5.83697898097787e-11

Iteration 0, Residual norm: 0.4060808625037903
Iteration 1, Residual norm: 0.016412401901058814
Iteration 2, Residual norm: 2.2548697356599077e-05
Iteration 3, Residual norm: 5.431822733633203e-11

Iteration 0, Residual norm: 0.4065703173159433
Iteration 1, Residual norm: 0.016618034217965243
Iteration 2, Residual norm: 2.3827133992731177e-05
Iteration 3, Residual norm: 6.354850974441053e-11

Iteration 0, Residual norm: 0.40752816311595746
Iteration 1, Residual norm: 0.017003664086637678
Iteration 2, Residual norm: 2.727481475456398e-05
Iteration 3, Residual norm: 9.591051932837998e-11

Iteration 0, Residual norm: 0.40839592564492616
Iteration 1, Residual norm: 0.01728133243226889
Iteration 2, Residual norm: 2.822125271630404e-05
Iteration 3, Residual norm: 1.1685057783572002e-10

Iteration 0, Residual norm: 0.4087909020839234
Iteration 1, Residual norm: 0.017372922197447443
Iteration 2, Residual norm: 2.98846160006334e-05
Iteration 3, Residual norm: 1.304101453026035e-10
Progress: 78%, Time: 7.80e+00, accepted step: True, dt: 2.00e-02, iterations: 3

Iteration 0, Residual norm: 0.40944536693308187
Iteration 1, Residual norm: 0.017239646399709794
Iteration 2, Residual norm: 2.8720820456757373e-05
Iteration 3, Residual norm: 1.2565976743418334e-10

Iteration 0, Residual norm: 0.4099670392209871
Iteration 1, Residual norm: 0.01704898613037333
Iteration 2, Residual norm: 2.7465688495866576e-05
Iteration 3, Residual norm: 1.1283224281204688e-10

Iteration 0, Residual norm: 0.4104344999490957
Iteration 1, Residual norm: 0.016761844694985355
Iteration 2, Residual norm: 2.5332710405395033e-05
Iteration 3, Residual norm: 8.668449520582284e-11

Iteration 0, Residual norm: 0.4102927642133465
Iteration 1, Residual norm: 0.016588497581987945
Iteration 2, Residual norm: 2.4044087684132584e-05
Iteration 3, Residual norm: 6.678033830177554e-11

Iteration 0, Residual norm: 0.4104137621172633
Iteration 1, Residual norm: 0.016581576243537097
Iteration 2, Residual norm: 2.3748317180749554e-05
Iteration 3, Residual norm: 5.912891458364579e-11

Iteration 0, Residual norm: 0.41120653387211303
Iteration 1, Residual norm: 0.016823042450419583
Iteration 2, Residual norm: 2.5570942220570745e-05
Iteration 3, Residual norm: 7.447008091910134e-11

Iteration 0, Residual norm: 0.4114912112171119
Iteration 1, Residual norm: 0.01708614144295133
Iteration 2, Residual norm: 2.8447526125158015e-05
Iteration 3, Residual norm: 1.0812187058094565e-10

Iteration 0, Residual norm: 0.4115874391238847
Iteration 1, Residual norm: 0.017429437365670878
Iteration 2, Residual norm: 3.043488179255112e-05
Iteration 3, Residual norm: 1.4383454088876058e-10

Iteration 0, Residual norm: 0.41203820803181374
Iteration 1, Residual norm: 0.01745783989117173
Iteration 2, Residual norm: 3.0998250895941207e-05
Iteration 3, Residual norm: 1.4833406315113367e-10

Iteration 0, Residual norm: 0.41212736086454704
Iteration 1, Residual norm: 0.017413752154974985
Iteration 2, Residual norm: 3.0706072027215694e-05
Iteration 3, Residual norm: 1.5070550914269238e-10

Iteration 0, Residual norm: 0.4120465262837663
Iteration 1, Residual norm: 0.017216093522522912
Iteration 2, Residual norm: 2.9377783150836624e-05
Iteration 3, Residual norm: 1.163341060460791e-10

Iteration 0, Residual norm: 0.4118907232270895
Iteration 1, Residual norm: 0.01694869710630076
Iteration 2, Residual norm: 2.629113475528651e-05
Iteration 3, Residual norm: 1.0013473727601188e-10

Iteration 0, Residual norm: 0.4115630854271427
Iteration 1, Residual norm: 0.016737509390415753
Iteration 2, Residual norm: 2.5065270994137417e-05
Iteration 3, Residual norm: 7.119179224542511e-11

Iteration 0, Residual norm: 0.41171210960064075
Iteration 1, Residual norm: 0.016946572633617627
Iteration 2, Residual norm: 2.5668024982972324e-05
Iteration 3, Residual norm: 7.820020345037334e-11

Iteration 0, Residual norm: 0.41181607617825955
Iteration 1, Residual norm: 0.017316494413922162
Iteration 2, Residual norm: 2.7190795237236873e-05
Iteration 3, Residual norm: 1.0219575927181364e-10

Iteration 0, Residual norm: 0.4123239041451812
Iteration 1, Residual norm: 0.01764909247749405
Iteration 2, Residual norm: 3.1826604171003235e-05
Iteration 3, Residual norm: 1.4828553995076029e-10

Iteration 0, Residual norm: 0.41320143132111353
Iteration 1, Residual norm: 0.017948902609989277
Iteration 2, Residual norm: 3.3660852876344655e-05
Iteration 3, Residual norm: 1.7809771209009572e-10

Iteration 0, Residual norm: 0.4134551686164413
Iteration 1, Residual norm: 0.018006650308748557
Iteration 2, Residual norm: 3.317666315840499e-05
Iteration 3, Residual norm: 1.8778656519682345e-10

Iteration 0, Residual norm: 0.4138552598377585
Iteration 1, Residual norm: 0.01774398824033158
Iteration 2, Residual norm: 3.225859527115808e-05
Iteration 3, Residual norm: 1.608869670217203e-10

Iteration 0, Residual norm: 0.41414145490720394
Iteration 1, Residual norm: 0.01737712240907362
Iteration 2, Residual norm: 2.9490993212351977e-05
Iteration 3, Residual norm: 1.35750297978993e-10

Iteration 0, Residual norm: 0.4146399964449795
Iteration 1, Residual norm: 0.01710679773730588
Iteration 2, Residual norm: 2.7083324030337245e-05
Iteration 3, Residual norm: 1.0109826315636842e-10

Iteration 0, Residual norm: 0.4151472763910619
Iteration 1, Residual norm: 0.01689848295451137
Iteration 2, Residual norm: 2.5583503943118743e-05
Iteration 3, Residual norm: 7.716476185146177e-11

Iteration 0, Residual norm: 0.41502117480303075
Iteration 1, Residual norm: 0.016987243455447575
Iteration 2, Residual norm: 2.6276439565350736e-05
Iteration 3, Residual norm: 8.050910142408828e-11

Iteration 0, Residual norm: 0.4153563272045062
Iteration 1, Residual norm: 0.017257087198348665
Iteration 2, Residual norm: 2.875772027136696e-05
Iteration 3, Residual norm: 1.1757378497198869e-10

Iteration 0, Residual norm: 0.41618662503441545
Iteration 1, Residual norm: 0.017662479717493875
Iteration 2, Residual norm: 3.192004740177442e-05
Iteration 3, Residual norm: 1.607418338744861e-10
Progress: 83%, Time: 8.30e+00, accepted step: True, dt: 2.00e-02, iterations: 3

Iteration 0, Residual norm: 0.41669492059614677
Iteration 1, Residual norm: 0.01786142302941864
Iteration 2, Residual norm: 3.464775809482919e-05
Iteration 3, Residual norm: 1.9119939783959655e-10

Iteration 0, Residual norm: 0.41670381931134226
Iteration 1, Residual norm: 0.017983506972546998
Iteration 2, Residual norm: 3.5256249323714746e-05
Iteration 3, Residual norm: 1.999691361489253e-10

Iteration 0, Residual norm: 0.4169949449333513
Iteration 1, Residual norm: 0.01773662533933171
Iteration 2, Residual norm: 3.286450583227901e-05
Iteration 3, Residual norm: 1.8608138411909428e-10

Iteration 0, Residual norm: 0.4167046777068769
Iteration 1, Residual norm: 0.01743045039976055
Iteration 2, Residual norm: 2.9962014605499046e-05
Iteration 3, Residual norm: 1.4616369659525365e-10

Iteration 0, Residual norm: 0.4165400256711665
Iteration 1, Residual norm: 0.017148240016162455
Iteration 2, Residual norm: 2.7991287562473867e-05
Iteration 3, Residual norm: 1.1149870816261615e-10

Iteration 0, Residual norm: 0.41634511148608067
Iteration 1, Residual norm: 0.017047914994805454
Iteration 2, Residual norm: 2.6328753265978224e-05
Iteration 3, Residual norm: 8.170190633297555e-11

Iteration 0, Residual norm: 0.41632721561995617
Iteration 1, Residual norm: 0.01702008267790791
Iteration 2, Residual norm: 2.713709646085188e-05
Iteration 3, Residual norm: 9.091585397880879e-11

Iteration 0, Residual norm: 0.4161985025770361
Iteration 1, Residual norm: 0.017581133543019982
Iteration 2, Residual norm: 3.0017607180876016e-05
Iteration 3, Residual norm: 1.4625407034498679e-10

Iteration 0, Residual norm: 0.4166706670421025
Iteration 1, Residual norm: 0.018044025661028296
Iteration 2, Residual norm: 3.324189865915553e-05
Iteration 3, Residual norm: 1.9341503338243857e-10

Iteration 0, Residual norm: 0.4174892205117695
Iteration 1, Residual norm: 0.01820388932817171
Iteration 2, Residual norm: 3.7780146155880686e-05
Iteration 3, Residual norm: 2.2874722495496904e-10

Iteration 0, Residual norm: 0.41775204045740905
Iteration 1, Residual norm: 0.018203997563116733
Iteration 2, Residual norm: 3.540139939129065e-05
Iteration 3, Residual norm: 2.3523339357941666e-10

Iteration 0, Residual norm: 0.4177589244490357
Iteration 1, Residual norm: 0.01799247640804579
Iteration 2, Residual norm: 3.396927261091286e-05
Iteration 3, Residual norm: 2.0289965166999394e-10

Iteration 0, Residual norm: 0.4183992303692509
Iteration 1, Residual norm: 0.01759031733883261
Iteration 2, Residual norm: 3.064451547014079e-05
Iteration 3, Residual norm: 1.5560732647910388e-10

Iteration 0, Residual norm: 0.41885511363547845
Iteration 1, Residual norm: 0.017195911587761726
Iteration 2, Residual norm: 2.8245055914874467e-05
Iteration 3, Residual norm: 1.1648766776978382e-10

Iteration 0, Residual norm: 0.41956624330142916
Iteration 1, Residual norm: 0.017001926853824987
Iteration 2, Residual norm: 2.6769765346779148e-05
Iteration 3, Residual norm: 8.705581596926934e-11

Iteration 0, Residual norm: 0.42022057584390293
Iteration 1, Residual norm: 0.017179590697032624
Iteration 2, Residual norm: 2.7503159275411556e-05
Iteration 3, Residual norm: 9.732467865242452e-11

Iteration 0, Residual norm: 0.41996550699223256
Iteration 1, Residual norm: 0.017446087373421217
Iteration 2, Residual norm: 3.0765947216297945e-05
Iteration 3, Residual norm: 1.462304482037592e-10

Iteration 0, Residual norm: 0.42063375294204125
Iteration 1, Residual norm: 0.01796149675362013
Iteration 2, Residual norm: 3.5048384032828e-05
Iteration 3, Residual norm: 2.033691465752993e-10

Iteration 0, Residual norm: 0.4210260466378244
Iteration 1, Residual norm: 0.01809796230789082
Iteration 2, Residual norm: 3.686728828107792e-05
Iteration 3, Residual norm: 2.385652542329931e-10

Iteration 0, Residual norm: 0.42147139118710464
Iteration 1, Residual norm: 0.018172321900783762
Iteration 2, Residual norm: 3.5939244869163746e-05
Iteration 3, Residual norm: 2.498782014565684e-10

Iteration 0, Residual norm: 0.4214460897955869
Iteration 1, Residual norm: 0.01787304233328903
Iteration 2, Residual norm: 3.4723696130962444e-05
Iteration 3, Residual norm: 2.0591669543652033e-10

Iteration 0, Residual norm: 0.42134253683835593
Iteration 1, Residual norm: 0.01759695522171612
Iteration 2, Residual norm: 3.183763079490126e-05
Iteration 3, Residual norm: 1.7026426613482295e-10

Iteration 0, Residual norm: 0.4210921811865449
Iteration 1, Residual norm: 0.01723592758688025
Iteration 2, Residual norm: 2.9053990143556278e-05
Iteration 3, Residual norm: 1.164638237199717e-10

Iteration 0, Residual norm: 0.42057819863794105
Iteration 1, Residual norm: 0.017041033723668324
Iteration 2, Residual norm: 2.649339179968828e-05
Iteration 3, Residual norm: 9.484707158804887e-11

Iteration 0, Residual norm: 0.42063253516513144
Iteration 1, Residual norm: 0.017086246985513396
Iteration 2, Residual norm: 2.7693480604851645e-05
Iteration 3, Residual norm: 9.741518998635991e-11
Progress: 88%, Time: 8.80e+00, accepted step: True, dt: 2.00e-02, iterations: 3

Iteration 0, Residual norm: 0.42093734443448555
Iteration 1, Residual norm: 0.01756101332063506
Iteration 2, Residual norm: 3.144672573457703e-05
Iteration 3, Residual norm: 1.8357202148210453e-10

Iteration 0, Residual norm: 0.42067609957178886
Iteration 1, Residual norm: 0.0181166002965844
Iteration 2, Residual norm: 3.4094800616718215e-05
Iteration 3, Residual norm: 2.4796255071052103e-10

Iteration 0, Residual norm: 0.4210023995246885
Iteration 1, Residual norm: 0.018390222488794596
Iteration 2, Residual norm: 3.951974988471211e-05
Iteration 3, Residual norm: 2.83513739170702e-10

Iteration 0, Residual norm: 0.421778695176264
Iteration 1, Residual norm: 0.01835132538157983
Iteration 2, Residual norm: 3.7931998919551974e-05
Iteration 3, Residual norm: 2.681141908466191e-10

Iteration 0, Residual norm: 0.42139676835313533
Iteration 1, Residual norm: 0.018052923295085694
Iteration 2, Residual norm: 3.499725445684995e-05
Iteration 3, Residual norm: 2.3651380607538845e-10

Iteration 0, Residual norm: 0.42174311613750703
Iteration 1, Residual norm: 0.017604986071537675
Iteration 2, Residual norm: 3.1800486355468845e-05
Iteration 3, Residual norm: 1.735171361426053e-10

Iteration 0, Residual norm: 0.4219847885837681
Iteration 1, Residual norm: 0.017260377257949502
Iteration 2, Residual norm: 2.867299980279622e-05
Iteration 3, Residual norm: 1.2312540752900835e-10

Iteration 0, Residual norm: 0.42273515540546264
Iteration 1, Residual norm: 0.017101472473301217
Iteration 2, Residual norm: 2.6857707427371766e-05
Iteration 3, Residual norm: 9.035370452021041e-11

Iteration 0, Residual norm: 0.4237801051802359
Iteration 1, Residual norm: 0.017138760032933822
Iteration 2, Residual norm: 2.8136282187132693e-05
Iteration 3, Residual norm: 1.0095885752115298e-10

Iteration 0, Residual norm: 0.4235630009439897
Iteration 1, Residual norm: 0.01757907768748693
Iteration 2, Residual norm: 3.202474980269318e-05
Iteration 3, Residual norm: 1.6846480255304481e-10

Iteration 0, Residual norm: 0.4238580885553052
Iteration 1, Residual norm: 0.018004903657576666
Iteration 2, Residual norm: 3.537089535104391e-05
Iteration 3, Residual norm: 2.459862993086491e-10

Iteration 0, Residual norm: 0.42410813518641866
Iteration 1, Residual norm: 0.018261581029300304
Iteration 2, Residual norm: 3.825751839996686e-05
Iteration 3, Residual norm: 2.7263482627685223e-10

Iteration 0, Residual norm: 0.4245580928369263
Iteration 1, Residual norm: 0.01825442552643422
Iteration 2, Residual norm: 3.885470668783002e-05
Iteration 3, Residual norm: 2.870252775223323e-10

Iteration 0, Residual norm: 0.4247170443210456
Iteration 1, Residual norm: 0.018018048935821575
Iteration 2, Residual norm: 3.642467359082057e-05
Iteration 3, Residual norm: 2.236495093551401e-10

Iteration 0, Residual norm: 0.42451106449001613
Iteration 1, Residual norm: 0.017622168528767114
Iteration 2, Residual norm: 3.18225215857668e-05
Iteration 3, Residual norm: 1.958635918202975e-10

Iteration 0, Residual norm: 0.42401896004703377
Iteration 1, Residual norm: 0.017199156160372348
Iteration 2, Residual norm: 2.8990368766564812e-05
Iteration 3, Residual norm: 1.2614283416063776e-10

Iteration 0, Residual norm: 0.4238167622120945
Iteration 1, Residual norm: 0.017096633106383062
Iteration 2, Residual norm: 2.7227471445097613e-05
Iteration 3, Residual norm: 1.0233641773115833e-10

Iteration 0, Residual norm: 0.423434830758013
Iteration 1, Residual norm: 0.01712571918773703
Iteration 2, Residual norm: 2.807774010070889e-05
Iteration 3, Residual norm: 1.0186307330556102e-10

Iteration 0, Residual norm: 0.42357998360099636
Iteration 1, Residual norm: 0.01741744991510641
Iteration 2, Residual norm: 3.136670002600557e-05
Iteration 3, Residual norm: 1.7619440254005413e-10

Iteration 0, Residual norm: 0.42322283379729747
Iteration 1, Residual norm: 0.018095129413680563
Iteration 2, Residual norm: 3.582288203510395e-05
Iteration 3, Residual norm: 2.6728399918392285e-10

Iteration 0, Residual norm: 0.42344802369063456
Iteration 1, Residual norm: 0.01848305560946816
Iteration 2, Residual norm: 3.8482583938376394e-05
Iteration 3, Residual norm: 3.077423186909387e-10

Iteration 0, Residual norm: 0.42394620395546834
Iteration 1, Residual norm: 0.018369082187567293
Iteration 2, Residual norm: 4.033822305457379e-05
Iteration 3, Residual norm: 2.9077570082833885e-10

Iteration 0, Residual norm: 0.42381887945223645
Iteration 1, Residual norm: 0.018093477216108513
Iteration 2, Residual norm: 3.589868028615305e-05
Iteration 3, Residual norm: 2.5909749034016815e-10

Iteration 0, Residual norm: 0.4237132825436826
Iteration 1, Residual norm: 0.01767976103868201
Iteration 2, Residual norm: 3.252962595415096e-05
Iteration 3, Residual norm: 1.909188398288026e-10

Iteration 0, Residual norm: 0.4239369762045043
Iteration 1, Residual norm: 0.01729633741772239
Iteration 2, Residual norm: 2.892691860494848e-05
Iteration 3, Residual norm: 1.3477710956473782e-10
Progress: 93%, Time: 9.30e+00, accepted step: True, dt: 2.00e-02, iterations: 3

Iteration 0, Residual norm: 0.4243428045991893
Iteration 1, Residual norm: 0.01707897344677241
Iteration 2, Residual norm: 2.7418920196119438e-05
Iteration 3, Residual norm: 1.0011185916489635e-10

Iteration 0, Residual norm: 0.4254141266703298
Iteration 1, Residual norm: 0.017119732652096233
Iteration 2, Residual norm: 2.8175993462394212e-05
Iteration 3, Residual norm: 1.0671956912343411e-10

Iteration 0, Residual norm: 0.42587172697799863
Iteration 1, Residual norm: 0.0175725058579319
Iteration 2, Residual norm: 3.160212683914248e-05
Iteration 3, Residual norm: 1.7808229973317441e-10

Iteration 0, Residual norm: 0.4255537565835712
Iteration 1, Residual norm: 0.01795170965704576
Iteration 2, Residual norm: 3.645922531667446e-05
Iteration 3, Residual norm: 2.552388754726858e-10

Iteration 0, Residual norm: 0.4258201535415503
Iteration 1, Residual norm: 0.018369371731933485
Iteration 2, Residual norm: 4.0242656305981375e-05
Iteration 3, Residual norm: 3.0306833907654406e-10

Iteration 0, Residual norm: 0.4258973615612214
Iteration 1, Residual norm: 0.018295936398596044
Iteration 2, Residual norm: 3.9415946252359523e-05
Iteration 3, Residual norm: 3.073389054219753e-10

Iteration 0, Residual norm: 0.42631233206848873
Iteration 1, Residual norm: 0.018068918491444735
Iteration 2, Residual norm: 3.616834910240143e-05
Iteration 3, Residual norm: 2.6633896552091727e-10

Iteration 0, Residual norm: 0.4262210860520482
Iteration 1, Residual norm: 0.017630567848249273
Iteration 2, Residual norm: 3.314400181681357e-05
Iteration 3, Residual norm: 2.0917821758079431e-10

Iteration 0, Residual norm: 0.4258112878471475
Iteration 1, Residual norm: 0.017308223495916465
Iteration 2, Residual norm: 2.985575471993474e-05
Iteration 3, Residual norm: 1.415420342340747e-10

Iteration 0, Residual norm: 0.4255469588972253
Iteration 1, Residual norm: 0.01706400736641649
Iteration 2, Residual norm: 2.799349882633769e-05
Iteration 3, Residual norm: 1.0486581682824197e-10

Iteration 0, Residual norm: 0.4247844038534671
Iteration 1, Residual norm: 0.017072876247755834
Iteration 2, Residual norm: 2.7560635579345143e-05
Iteration 3, Residual norm: 1.0800337044064268e-10

Iteration 0, Residual norm: 0.42489082320098864
Iteration 1, Residual norm: 0.0173351574090911
Iteration 2, Residual norm: 3.119187682693767e-05
Iteration 3, Residual norm: 1.6065129389707063e-10

Iteration 0, Residual norm: 0.42494403811172815
Iteration 1, Residual norm: 0.01794685835629224
Iteration 2, Residual norm: 3.623460551501993e-05
Iteration 3, Residual norm: 2.9587320606556713e-10

Iteration 0, Residual norm: 0.4245690063796409
Iteration 1, Residual norm: 0.018394938251536763
Iteration 2, Residual norm: 3.764325765587482e-05
Iteration 3, Residual norm: 3.3986816649338886e-10

Iteration 0, Residual norm: 0.4246168206956899
Iteration 1, Residual norm: 0.018432127649213526
Iteration 2, Residual norm: 4.1716881666903654e-05
Iteration 3, Residual norm: 3.3368786663316344e-10

Iteration 0, Residual norm: 0.42509982393272144
Iteration 1, Residual norm: 0.01822290089362166
Iteration 2, Residual norm: 3.727750960370478e-05
Iteration 3, Residual norm: 2.841421948946759e-10

Iteration 0, Residual norm: 0.4247032002608616
Iteration 1, Residual norm: 0.01774746052243213
Iteration 2, Residual norm: 3.3524922596765166e-05
Iteration 3, Residual norm: 2.2463990927446399e-10

Iteration 0, Residual norm: 0.4248787762798086
Iteration 1, Residual norm: 0.017271777932571043
Iteration 2, Residual norm: 2.9962194267327104e-05
Iteration 3, Residual norm: 1.564203311145294e-10

Iteration 0, Residual norm: 0.42498658368865366
Iteration 1, Residual norm: 0.017078331898525638
Iteration 2, Residual norm: 2.774726830759237e-05
Iteration 3, Residual norm: 1.0845560740337895e-10

Iteration 0, Residual norm: 0.4258475039039725
Iteration 1, Residual norm: 0.01709889426708523
Iteration 2, Residual norm: 2.7516384134032294e-05
Iteration 3, Residual norm: 1.0298470493970621e-10

Iteration 0, Residual norm: 0.4269170343393192
Iteration 1, Residual norm: 0.017440992272868723
Iteration 2, Residual norm: 3.133597089492727e-05
Iteration 3, Residual norm: 1.6423286059011267e-10

Iteration 0, Residual norm: 0.42645521527021174
Iteration 1, Residual norm: 0.017929677523024958
Iteration 2, Residual norm: 3.6893462388842486e-05
Iteration 3, Residual norm: 2.5495999069570405e-10

Iteration 0, Residual norm: 0.42662314118543737
Iteration 1, Residual norm: 0.01834314473128481
Iteration 2, Residual norm: 3.92080356918143e-05
Iteration 3, Residual norm: 3.3461727961877216e-10

Iteration 0, Residual norm: 0.4264418184040563
Iteration 1, Residual norm: 0.018334071590209278
Iteration 2, Residual norm: 3.9826232200902204e-05
Iteration 3, Residual norm: 3.231740367851326e-10

Iteration 0, Residual norm: 0.4267871512276016
Iteration 1, Residual norm: 0.018157268339797288
Iteration 2, Residual norm: 3.8481828806413826e-05
Iteration 3, Residual norm: 3.0655220574963335e-10
Progress: 98%, Time: 9.80e+00, accepted step: True, dt: 2.00e-02, iterations: 3

Iteration 0, Residual norm: 0.4269889464938852
Iteration 1, Residual norm: 0.01775517039612255
Iteration 2, Residual norm: 3.5154722209351746e-05
Iteration 3, Residual norm: 2.1352904149585835e-10

Iteration 0, Residual norm: 0.4266299941396428
Iteration 1, Residual norm: 0.017362885041909915
Iteration 2, Residual norm: 3.01293361652026e-05
Iteration 3, Residual norm: 1.6996168140680265e-10

Iteration 0, Residual norm: 0.42614661765642986
Iteration 1, Residual norm: 0.016986471799577383
Iteration 2, Residual norm: 2.7877214566973985e-05
Iteration 3, Residual norm: 1.0895402198804907e-10

Iteration 0, Residual norm: 0.425777115853739
Iteration 1, Residual norm: 0.017063893878073392
Iteration 2, Residual norm: 2.719539537360905e-05
Iteration 3, Residual norm: 1.1020699101857853e-10

Iteration 0, Residual norm: 0.42547001812925295
Iteration 1, Residual norm: 0.017256991253760325
Iteration 2, Residual norm: 3.081449032900793e-05
Iteration 3, Residual norm: 1.5052922934999215e-10

Iteration 0, Residual norm: 0.4255913908404619
Iteration 1, Residual norm: 0.017741379551166506
Iteration 2, Residual norm: 3.535949591015188e-05
Iteration 3, Residual norm: 2.724639010204767e-10

Iteration 0, Residual norm: 0.4249685214518685
Iteration 1, Residual norm: 0.018306637724405248
Iteration 2, Residual norm: 3.868210146840206e-05
Iteration 3, Residual norm: 3.490612596662382e-10

Iteration 0, Residual norm: 0.4250377122869643
Iteration 1, Residual norm: 0.0185229747928201
Iteration 2, Residual norm: 4.103915123435495e-05
Iteration 3, Residual norm: 3.7066783606339863e-10

Iteration 0, Residual norm: 0.42547423713985194
Iteration 1, Residual norm: 0.018257663296971503
Iteration 2, Residual norm: 3.93216630606771e-05
Iteration 3, Residual norm: 3.047449010625403e-10

Iteration 0, Residual norm: 0.4250873255862509
Iteration 1, Residual norm: 0.017814957943792287
Iteration 2, Residual norm: 3.471221081819317e-05
Iteration 3, Residual norm: 2.44932417622294e-10

Postprocessing

output_dir = "./output_navier_stokes"
os.makedirs(output_dir, exist_ok=True)

history = result.history
timesteps = history.t
velocity_history = history.q['1velocity']
pressure_history = history.q['2pressure']
temperature_history = history.q['3temperature']

# Export one VTK file for each time step
for i, t_val in enumerate(timesteps):
  # Export velocity and temperature
  point_data = {
      'velocity': jnp.pad(velocity_history[i], ((0, 0), (0, 1)), constant_values=0),
      'temperature': temperature_history[i],
  }
  meshio.Mesh(
      points=jnp.pad(coords_v, ((0, 0), (0, 1)), constant_values=0),
      cells={
          'triangle6': cells_v
      },
      point_data=point_data,
  ).write(f"{output_dir}/navier_stokes_t{i}.vtk")

  # Export pressure
  point_data = {
      'pressure': pressure_history[i],
  }
  meshio.Mesh(
      points=jnp.pad(coords_p, ((0, 0), (0, 1)), constant_values=0),
      cells={
          'triangle': cells_p
      },
      point_data=point_data,
  ).write(f"{output_dir}/navier_stokes_p_t{i}.vtk")