Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Jump to bottom

Nonlinear solver does not converge in some models with nullspace removal #5371

Open
jdannberg opened this issue Aug 23, 2023 · 0 comments
Open

Nonlinear solver does not converge in some models with nullspace removal #5371

jdannberg opened this issue Aug 23, 2023 · 0 comments
Labels
bug

Comments

@jdannberg
Copy link
Contributor

I have a very simple models setup (copied below) where the nonlinear solver does not converge (at all). The model has a prescribed horizontal stress at the top and bottom (traction boundary conditions), and periodic boundary conditions on the sides, and therefore needs nullspace removal. The viscosity is constant and the density is constant, so there is no nonlinearity. When I switch off the nullspace removal, the nonlinear solver converges without any problem (but this is of course not an option in more complex models). The solution is correct in both cases, as far as I can tell.

Specifically, the residual that does not decrease is the Stokes residual (the temperature is constant as well). I think something might be wrong with how we compute the residual when we use nullspace removal, because the nonlinear solver converges when I use the iterated advection and defect correction Stokes solver (which has its own, different residual computation). I found this quite confusing, because removing the nullspace should not change anything about the residual.

I tried a lot of things and could not figure out what exactly the problem was, but I wondered if the problem might be some interaction with setting the constraints (because it seems the degree of freedom in the bottom left corner is constrained in models with nullspace removal, so when I call current_constraints.set_zero and then current_constraints.distribute without solving in between, the velocity would always be zero on that degree of freedom).

set Dimension                              = 2
set End time                               = 20
set Output directory                       = output-constant-visc
set Use years in output instead of seconds = true
set Maximum first time step                = 10
set Maximum relative increase in time step = 1
set Nonlinear solver scheme                = iterated Advection and Stokes

set Surface pressure                       = 0
set Adiabatic surface temperature          = 1100

set Max nonlinear iterations = 50

subsection Geometry model
  set Model name = box

  subsection Box
    set X extent = 100000
    set Y extent = 100000
    set X periodic = true
  end
end

subsection Boundary temperature model
  set Fixed temperature boundary indicators   = top,bottom
end


subsection Boundary temperature model
  set List of model names = initial temperature
  subsection Initial temperature
    set Minimal temperature = 1100
  end
end



subsection Gravity model
  set Model name = vertical
  subsection Vertical
    set Magnitude = 0.0
  end
end


subsection Boundary traction model
  set Prescribed traction boundary indicators = top x: function, bottom x: function

  subsection Function
    set Variable names = x,y
    # We want to prescribe a tangential traction of 50 MPa
    # at the top boundary and -50 at the bottom boundary.
    set  Function constants = traction = 50e6, y_extent = 1e5
    set Function expression = if (y>y_extent/2, traction, -traction) ;  0
  end
end


subsection Boundary velocity model
  set Prescribed velocity boundary indicators = top y: function, bottom y: function

  subsection Function
    set Variable names = x,y
    set Function expression = 0;0
  end
end


subsection Initial temperature model
  set Model name = adiabatic

  subsection Adiabatic
    set Age top boundary layer      = 0
  end
end


subsection Material model
  set Model name = simpler

  subsection Simpler model
    set Reference density                = 3400
    set Thermal conductivity             = 0
    set Thermal expansion coefficient    = 0
    set Reference temperature            = 1100
    set Viscosity                        = 1e22
  end
end


subsection Mesh refinement
  set Initial adaptive refinement        = 0
  set Initial global refinement          = 1
  set Time steps between mesh refinement = 0
end


subsection Postprocess
  set List of postprocessors = visualization
  subsection Visualization
    set List of output variables = viscosity, density, shear stress, stress, stress second invariant
    set Time between graphical output = 0
  end
end


subsection Nullspace removal
  set Remove nullspace = net x translation
end

Here is the log output from that model run (with the current main):

-----------------------------------------------------------------------------
-- This is ASPECT, the Advanced Solver for Problems in Earth's ConvecTion.
--     . version 2.6.0-pre (main, 063741d02)
--     . using deal.II 9.5.0
--     .       with 32 bit indices and vectorization level 2 (256 bits)
--     . using Trilinos 13.2.0
--     . using p4est 2.3.2
--     . running in DEBUG mode
--     . running with 1 MPI process
-----------------------------------------------------------------------------

-----------------------------------------------------------------------------
-- For information on how to cite ASPECT, see:
--   https://aspect.geodynamics.org/citing.html?ver=2.6.0-pre&sha=063741d02&src=code
-----------------------------------------------------------------------------
Number of active cells: 4 (on 2 levels)
Number of degrees of freedom: 84 (50+9+25)

*** Timestep 0:  t=0 years, dt=0 years
   Solving temperature system... 0 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 1+0 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.0151e-16, 1
      Relative nonlinear residual (total system) after nonlinear iteration 1: 1

   Solving temperature system... 0 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 7+0 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.0151e-16, 0.928307
      Relative nonlinear residual (total system) after nonlinear iteration 2: 0.928307

   Solving temperature system... 0 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 7+0 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.0151e-16, 0.928307
      Relative nonlinear residual (total system) after nonlinear iteration 3: 0.928307

   Solving temperature system... 0 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 7+0 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.0151e-16, 0.928307
      Relative nonlinear residual (total system) after nonlinear iteration 4: 0.928307

   Solving temperature system... 0 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 7+0 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.0151e-16, 0.928307
      Relative nonlinear residual (total system) after nonlinear iteration 5: 0.928307

   Solving temperature system... 0 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 7+0 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.0151e-16, 0.928307
      Relative nonlinear residual (total system) after nonlinear iteration 6: 0.928307

   Solving temperature system... 0 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 7+0 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.0151e-16, 0.928307
      Relative nonlinear residual (total system) after nonlinear iteration 7: 0.928307

   Solving temperature system... 0 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 7+0 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.0151e-16, 0.928307
      Relative nonlinear residual (total system) after nonlinear iteration 8: 0.928307

   Solving temperature system... 0 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 7+0 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.0151e-16, 0.928307
      Relative nonlinear residual (total system) after nonlinear iteration 9: 0.928307

   Solving temperature system... 0 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 7+0 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.0151e-16, 0.928307
      Relative nonlinear residual (total system) after nonlinear iteration 10: 0.928307

   Solving temperature system... 0 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 7+0 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.0151e-16, 0.928307
      Relative nonlinear residual (total system) after nonlinear iteration 11: 0.928307

   Solving temperature system... 0 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 7+0 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.0151e-16, 0.928307
      Relative nonlinear residual (total system) after nonlinear iteration 12: 0.928307

   Solving temperature system... 0 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 7+0 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.0151e-16, 0.928307
      Relative nonlinear residual (total system) after nonlinear iteration 13: 0.928307

   Solving temperature system... 0 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 7+0 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.0151e-16, 0.928307
      Relative nonlinear residual (total system) after nonlinear iteration 14: 0.928307

   Solving temperature system... 0 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 7+0 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.0151e-16, 0.928307
      Relative nonlinear residual (total system) after nonlinear iteration 15: 0.928307

   Solving temperature system... 0 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 7+0 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.0151e-16, 0.928307
      Relative nonlinear residual (total system) after nonlinear iteration 16: 0.928307

   Solving temperature system... 0 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 7+0 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.0151e-16, 0.928307
      Relative nonlinear residual (total system) after nonlinear iteration 17: 0.928307

   Solving temperature system... 0 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 7+0 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.0151e-16, 0.928307
      Relative nonlinear residual (total system) after nonlinear iteration 18: 0.928307

   Solving temperature system... 0 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 7+0 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.0151e-16, 0.928307
      Relative nonlinear residual (total system) after nonlinear iteration 19: 0.928307

   Solving temperature system... 0 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 7+0 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.0151e-16, 0.928307
      Relative nonlinear residual (total system) after nonlinear iteration 20: 0.928307

   Solving temperature system... 0 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 7+0 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.0151e-16, 0.928307
      Relative nonlinear residual (total system) after nonlinear iteration 21: 0.928307

   Solving temperature system... 0 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 7+0 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.0151e-16, 0.928307
      Relative nonlinear residual (total system) after nonlinear iteration 22: 0.928307

   Solving temperature system... 0 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 7+0 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.0151e-16, 0.928307
      Relative nonlinear residual (total system) after nonlinear iteration 23: 0.928307

   Solving temperature system... 0 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 7+0 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.0151e-16, 0.928307
      Relative nonlinear residual (total system) after nonlinear iteration 24: 0.928307

   Solving temperature system... 0 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 7+0 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.0151e-16, 0.928307
      Relative nonlinear residual (total system) after nonlinear iteration 25: 0.928307

   Solving temperature system... 0 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 7+0 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.0151e-16, 0.928307
      Relative nonlinear residual (total system) after nonlinear iteration 26: 0.928307

   Solving temperature system... 0 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 7+0 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.0151e-16, 0.928307
      Relative nonlinear residual (total system) after nonlinear iteration 27: 0.928307

   Solving temperature system... 0 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 7+0 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.0151e-16, 0.928307
      Relative nonlinear residual (total system) after nonlinear iteration 28: 0.928307

   Solving temperature system... 0 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 7+0 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.0151e-16, 0.928307
      Relative nonlinear residual (total system) after nonlinear iteration 29: 0.928307

   Solving temperature system... 0 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 7+0 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.0151e-16, 0.928307
      Relative nonlinear residual (total system) after nonlinear iteration 30: 0.928307

   Solving temperature system... 0 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 7+0 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.0151e-16, 0.928307
      Relative nonlinear residual (total system) after nonlinear iteration 31: 0.928307

   Solving temperature system... 0 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 7+0 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.0151e-16, 0.928307
      Relative nonlinear residual (total system) after nonlinear iteration 32: 0.928307

   Solving temperature system... 0 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 7+0 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.0151e-16, 0.928307
      Relative nonlinear residual (total system) after nonlinear iteration 33: 0.928307

   Solving temperature system... 0 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 7+0 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.0151e-16, 0.928307
      Relative nonlinear residual (total system) after nonlinear iteration 34: 0.928307

   Solving temperature system... 0 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 7+0 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.0151e-16, 0.928307
      Relative nonlinear residual (total system) after nonlinear iteration 35: 0.928307

   Solving temperature system... 0 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 7+0 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.0151e-16, 0.928307
      Relative nonlinear residual (total system) after nonlinear iteration 36: 0.928307

   Solving temperature system... 0 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 7+0 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.0151e-16, 0.928307
      Relative nonlinear residual (total system) after nonlinear iteration 37: 0.928307

   Solving temperature system... 0 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 7+0 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.0151e-16, 0.928307
      Relative nonlinear residual (total system) after nonlinear iteration 38: 0.928307

   Solving temperature system... 0 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 7+0 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.0151e-16, 0.928307
      Relative nonlinear residual (total system) after nonlinear iteration 39: 0.928307

   Solving temperature system... 0 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 7+0 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.0151e-16, 0.928307
      Relative nonlinear residual (total system) after nonlinear iteration 40: 0.928307

   Solving temperature system... 0 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 7+0 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.0151e-16, 0.928307
      Relative nonlinear residual (total system) after nonlinear iteration 41: 0.928307

   Solving temperature system... 0 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 7+0 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.0151e-16, 0.928307
      Relative nonlinear residual (total system) after nonlinear iteration 42: 0.928307

   Solving temperature system... 0 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 7+0 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.0151e-16, 0.928307
      Relative nonlinear residual (total system) after nonlinear iteration 43: 0.928307

   Solving temperature system... 0 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 7+0 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.0151e-16, 0.928307
      Relative nonlinear residual (total system) after nonlinear iteration 44: 0.928307

   Solving temperature system... 0 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 7+0 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.0151e-16, 0.928307
      Relative nonlinear residual (total system) after nonlinear iteration 45: 0.928307

   Solving temperature system... 0 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 7+0 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.0151e-16, 0.928307
      Relative nonlinear residual (total system) after nonlinear iteration 46: 0.928307

   Solving temperature system... 0 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 7+0 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.0151e-16, 0.928307
      Relative nonlinear residual (total system) after nonlinear iteration 47: 0.928307

   Solving temperature system... 0 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 7+0 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.0151e-16, 0.928307
      Relative nonlinear residual (total system) after nonlinear iteration 48: 0.928307

   Solving temperature system... 0 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 7+0 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.0151e-16, 0.928307
      Relative nonlinear residual (total system) after nonlinear iteration 49: 0.928307

   Solving temperature system... 0 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 7+0 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.0151e-16, 0.928307
      Relative nonlinear residual (total system) after nonlinear iteration 50: 0.928307


   Postprocessing:
     Writing graphical output: output-constant-visc/solution/solution-00000

*** Timestep 1:  t=10 years, dt=10 years
   Solving temperature system... 0 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 7+0 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.10756e-16, 0.928307
      Relative nonlinear residual (total system) after nonlinear iteration 1: 0.928307

   Solving temperature system... 0 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 7+0 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.10756e-16, 0.928307
      Relative nonlinear residual (total system) after nonlinear iteration 2: 0.928307

   Solving temperature system... 0 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 7+0 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.10756e-16, 0.928307
      Relative nonlinear residual (total system) after nonlinear iteration 3: 0.928307

   Solving temperature system... 0 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 7+0 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.10756e-16, 0.928307
      Relative nonlinear residual (total system) after nonlinear iteration 4: 0.928307

   Solving temperature system... 0 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 7+0 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.10756e-16, 0.928307
      Relative nonlinear residual (total system) after nonlinear iteration 5: 0.928307

   Solving temperature system... 0 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 7+0 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.10756e-16, 0.928307
      Relative nonlinear residual (total system) after nonlinear iteration 6: 0.928307

   Solving temperature system... 0 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 7+0 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.10756e-16, 0.928307
      Relative nonlinear residual (total system) after nonlinear iteration 7: 0.928307

   Solving temperature system... 0 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 7+0 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.10756e-16, 0.928307
      Relative nonlinear residual (total system) after nonlinear iteration 8: 0.928307

   Solving temperature system... 0 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 7+0 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.10756e-16, 0.928307
      Relative nonlinear residual (total system) after nonlinear iteration 9: 0.928307

   Solving temperature system... 0 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 7+0 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.10756e-16, 0.928307
      Relative nonlinear residual (total system) after nonlinear iteration 10: 0.928307

   Solving temperature system... 0 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 7+0 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.10756e-16, 0.928307
      Relative nonlinear residual (total system) after nonlinear iteration 11: 0.928307

   Solving temperature system... 0 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 7+0 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.10756e-16, 0.928307
      Relative nonlinear residual (total system) after nonlinear iteration 12: 0.928307

   Solving temperature system... 0 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 7+0 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.10756e-16, 0.928307
      Relative nonlinear residual (total system) after nonlinear iteration 13: 0.928307

   Solving temperature system... 0 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 7+0 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.10756e-16, 0.928307
      Relative nonlinear residual (total system) after nonlinear iteration 14: 0.928307

   Solving temperature system... 0 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 7+0 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.10756e-16, 0.928307
      Relative nonlinear residual (total system) after nonlinear iteration 15: 0.928307

   Solving temperature system... 0 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 7+0 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.10756e-16, 0.928307
      Relative nonlinear residual (total system) after nonlinear iteration 16: 0.928307

   Solving temperature system... 0 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 7+0 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.10756e-16, 0.928307
      Relative nonlinear residual (total system) after nonlinear iteration 17: 0.928307

   Solving temperature system... 0 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 7+0 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.10756e-16, 0.928307
      Relative nonlinear residual (total system) after nonlinear iteration 18: 0.928307

   Solving temperature system... 0 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 7+0 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.10756e-16, 0.928307
      Relative nonlinear residual (total system) after nonlinear iteration 19: 0.928307

   Solving temperature system... 0 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 7+0 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.10756e-16, 0.928307
      Relative nonlinear residual (total system) after nonlinear iteration 20: 0.928307

   Solving temperature system... 0 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 7+0 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.10756e-16, 0.928307
      Relative nonlinear residual (total system) after nonlinear iteration 21: 0.928307

   Solving temperature system... 0 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 7+0 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.10756e-16, 0.928307
      Relative nonlinear residual (total system) after nonlinear iteration 22: 0.928307

   Solving temperature system... 0 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 7+0 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.10756e-16, 0.928307
      Relative nonlinear residual (total system) after nonlinear iteration 23: 0.928307

   Solving temperature system... 0 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 7+0 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.10756e-16, 0.928307
      Relative nonlinear residual (total system) after nonlinear iteration 24: 0.928307

   Solving temperature system... 0 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 7+0 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.10756e-16, 0.928307
      Relative nonlinear residual (total system) after nonlinear iteration 25: 0.928307

   Solving temperature system... 0 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 7+0 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.10756e-16, 0.928307
      Relative nonlinear residual (total system) after nonlinear iteration 26: 0.928307

   Solving temperature system... 0 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 7+0 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.10756e-16, 0.928307
      Relative nonlinear residual (total system) after nonlinear iteration 27: 0.928307

   Solving temperature system... 0 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 7+0 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.10756e-16, 0.928307
      Relative nonlinear residual (total system) after nonlinear iteration 28: 0.928307

   Solving temperature system... 0 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 7+0 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.10756e-16, 0.928307
      Relative nonlinear residual (total system) after nonlinear iteration 29: 0.928307

   Solving temperature system... 0 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 7+0 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.10756e-16, 0.928307
      Relative nonlinear residual (total system) after nonlinear iteration 30: 0.928307

   Solving temperature system... 0 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 7+0 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.10756e-16, 0.928307
      Relative nonlinear residual (total system) after nonlinear iteration 31: 0.928307

   Solving temperature system... 0 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 7+0 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.10756e-16, 0.928307
      Relative nonlinear residual (total system) after nonlinear iteration 32: 0.928307

   Solving temperature system... 0 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 7+0 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.10756e-16, 0.928307
      Relative nonlinear residual (total system) after nonlinear iteration 33: 0.928307

   Solving temperature system... 0 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 7+0 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.10756e-16, 0.928307
      Relative nonlinear residual (total system) after nonlinear iteration 34: 0.928307

   Solving temperature system... 0 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 7+0 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.10756e-16, 0.928307
      Relative nonlinear residual (total system) after nonlinear iteration 35: 0.928307

   Solving temperature system... 0 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 7+0 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.10756e-16, 0.928307
      Relative nonlinear residual (total system) after nonlinear iteration 36: 0.928307

   Solving temperature system... 0 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 7+0 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.10756e-16, 0.928307
      Relative nonlinear residual (total system) after nonlinear iteration 37: 0.928307

   Solving temperature system... 0 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 7+0 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.10756e-16, 0.928307
      Relative nonlinear residual (total system) after nonlinear iteration 38: 0.928307

   Solving temperature system... 0 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 7+0 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.10756e-16, 0.928307
      Relative nonlinear residual (total system) after nonlinear iteration 39: 0.928307

   Solving temperature system... 0 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 7+0 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.10756e-16, 0.928307
      Relative nonlinear residual (total system) after nonlinear iteration 40: 0.928307

   Solving temperature system... 0 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 7+0 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.10756e-16, 0.928307
      Relative nonlinear residual (total system) after nonlinear iteration 41: 0.928307

   Solving temperature system... 0 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 7+0 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.10756e-16, 0.928307
      Relative nonlinear residual (total system) after nonlinear iteration 42: 0.928307

   Solving temperature system... 0 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 7+0 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.10756e-16, 0.928307
      Relative nonlinear residual (total system) after nonlinear iteration 43: 0.928307

   Solving temperature system... 0 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 7+0 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.10756e-16, 0.928307
      Relative nonlinear residual (total system) after nonlinear iteration 44: 0.928307

   Solving temperature system... 0 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 7+0 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.10756e-16, 0.928307
      Relative nonlinear residual (total system) after nonlinear iteration 45: 0.928307

   Solving temperature system... 0 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 7+0 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.10756e-16, 0.928307
      Relative nonlinear residual (total system) after nonlinear iteration 46: 0.928307

   Solving temperature system... 0 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 7+0 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.10756e-16, 0.928307
      Relative nonlinear residual (total system) after nonlinear iteration 47: 0.928307

   Solving temperature system... 0 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 7+0 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.10756e-16, 0.928307
      Relative nonlinear residual (total system) after nonlinear iteration 48: 0.928307

   Solving temperature system... 0 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 7+0 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.10756e-16, 0.928307
      Relative nonlinear residual (total system) after nonlinear iteration 49: 0.928307

   Solving temperature system... 0 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 7+0 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.10756e-16, 0.928307
      Relative nonlinear residual (total system) after nonlinear iteration 50: 0.928307


   Postprocessing:
     Writing graphical output: output-constant-visc/solution/solution-00001

*** Timestep 2:  t=20 years, dt=10 years
   Solving temperature system... 0 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 7+0 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.2874e-16, 0.928307
      Relative nonlinear residual (total system) after nonlinear iteration 1: 0.928307

   Solving temperature system... 0 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 7+0 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.2874e-16, 0.928307
      Relative nonlinear residual (total system) after nonlinear iteration 2: 0.928307

   Solving temperature system... 0 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 7+0 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.2874e-16, 0.928307
      Relative nonlinear residual (total system) after nonlinear iteration 3: 0.928307

   Solving temperature system... 0 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 7+0 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.2874e-16, 0.928307
      Relative nonlinear residual (total system) after nonlinear iteration 4: 0.928307

   Solving temperature system... 0 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 7+0 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.2874e-16, 0.928307
      Relative nonlinear residual (total system) after nonlinear iteration 5: 0.928307

   Solving temperature system... 0 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 7+0 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.2874e-16, 0.928307
      Relative nonlinear residual (total system) after nonlinear iteration 6: 0.928307

   Solving temperature system... 0 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 7+0 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.2874e-16, 0.928307
      Relative nonlinear residual (total system) after nonlinear iteration 7: 0.928307

   Solving temperature system... 0 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 7+0 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.2874e-16, 0.928307
      Relative nonlinear residual (total system) after nonlinear iteration 8: 0.928307

   Solving temperature system... 0 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 7+0 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.2874e-16, 0.928307
      Relative nonlinear residual (total system) after nonlinear iteration 9: 0.928307

   Solving temperature system... 0 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 7+0 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.2874e-16, 0.928307
      Relative nonlinear residual (total system) after nonlinear iteration 10: 0.928307

   Solving temperature system... 0 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 7+0 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.2874e-16, 0.928307
      Relative nonlinear residual (total system) after nonlinear iteration 11: 0.928307

   Solving temperature system... 0 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 7+0 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.2874e-16, 0.928307
      Relative nonlinear residual (total system) after nonlinear iteration 12: 0.928307

   Solving temperature system... 0 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 7+0 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.2874e-16, 0.928307
      Relative nonlinear residual (total system) after nonlinear iteration 13: 0.928307

   Solving temperature system... 0 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 7+0 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.2874e-16, 0.928307
      Relative nonlinear residual (total system) after nonlinear iteration 14: 0.928307

   Solving temperature system... 0 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 7+0 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.2874e-16, 0.928307
      Relative nonlinear residual (total system) after nonlinear iteration 15: 0.928307

   Solving temperature system... 0 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 7+0 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.2874e-16, 0.928307
      Relative nonlinear residual (total system) after nonlinear iteration 16: 0.928307

   Solving temperature system... 0 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 7+0 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.2874e-16, 0.928307
      Relative nonlinear residual (total system) after nonlinear iteration 17: 0.928307

   Solving temperature system... 0 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 7+0 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.2874e-16, 0.928307
      Relative nonlinear residual (total system) after nonlinear iteration 18: 0.928307

   Solving temperature system... 0 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 7+0 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.2874e-16, 0.928307
      Relative nonlinear residual (total system) after nonlinear iteration 19: 0.928307

   Solving temperature system... 0 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 7+0 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.2874e-16, 0.928307
      Relative nonlinear residual (total system) after nonlinear iteration 20: 0.928307

   Solving temperature system... 0 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 7+0 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.2874e-16, 0.928307
      Relative nonlinear residual (total system) after nonlinear iteration 21: 0.928307

   Solving temperature system... 0 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 7+0 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.2874e-16, 0.928307
      Relative nonlinear residual (total system) after nonlinear iteration 22: 0.928307

   Solving temperature system... 0 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 7+0 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.2874e-16, 0.928307
      Relative nonlinear residual (total system) after nonlinear iteration 23: 0.928307

   Solving temperature system... 0 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 7+0 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.2874e-16, 0.928307
      Relative nonlinear residual (total system) after nonlinear iteration 24: 0.928307

   Solving temperature system... 0 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 7+0 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.2874e-16, 0.928307
      Relative nonlinear residual (total system) after nonlinear iteration 25: 0.928307

   Solving temperature system... 0 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 7+0 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.2874e-16, 0.928307
      Relative nonlinear residual (total system) after nonlinear iteration 26: 0.928307

   Solving temperature system... 0 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 7+0 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.2874e-16, 0.928307
      Relative nonlinear residual (total system) after nonlinear iteration 27: 0.928307

   Solving temperature system... 0 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 7+0 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.2874e-16, 0.928307
      Relative nonlinear residual (total system) after nonlinear iteration 28: 0.928307

   Solving temperature system... 0 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 7+0 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.2874e-16, 0.928307
      Relative nonlinear residual (total system) after nonlinear iteration 29: 0.928307

   Solving temperature system... 0 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 7+0 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.2874e-16, 0.928307
      Relative nonlinear residual (total system) after nonlinear iteration 30: 0.928307

   Solving temperature system... 0 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 7+0 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.2874e-16, 0.928307
      Relative nonlinear residual (total system) after nonlinear iteration 31: 0.928307

   Solving temperature system... 0 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 7+0 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.2874e-16, 0.928307
      Relative nonlinear residual (total system) after nonlinear iteration 32: 0.928307

   Solving temperature system... 0 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 7+0 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.2874e-16, 0.928307
      Relative nonlinear residual (total system) after nonlinear iteration 33: 0.928307

   Solving temperature system... 0 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 7+0 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.2874e-16, 0.928307
      Relative nonlinear residual (total system) after nonlinear iteration 34: 0.928307

   Solving temperature system... 0 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 7+0 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.2874e-16, 0.928307
      Relative nonlinear residual (total system) after nonlinear iteration 35: 0.928307

   Solving temperature system... 0 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 7+0 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.2874e-16, 0.928307
      Relative nonlinear residual (total system) after nonlinear iteration 36: 0.928307

   Solving temperature system... 0 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 7+0 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.2874e-16, 0.928307
      Relative nonlinear residual (total system) after nonlinear iteration 37: 0.928307

   Solving temperature system... 0 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 7+0 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.2874e-16, 0.928307
      Relative nonlinear residual (total system) after nonlinear iteration 38: 0.928307

   Solving temperature system... 0 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 7+0 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.2874e-16, 0.928307
      Relative nonlinear residual (total system) after nonlinear iteration 39: 0.928307

   Solving temperature system... 0 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 7+0 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.2874e-16, 0.928307
      Relative nonlinear residual (total system) after nonlinear iteration 40: 0.928307

   Solving temperature system... 0 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 7+0 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.2874e-16, 0.928307
      Relative nonlinear residual (total system) after nonlinear iteration 41: 0.928307

   Solving temperature system... 0 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 7+0 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.2874e-16, 0.928307
      Relative nonlinear residual (total system) after nonlinear iteration 42: 0.928307

   Solving temperature system... 0 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 7+0 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.2874e-16, 0.928307
      Relative nonlinear residual (total system) after nonlinear iteration 43: 0.928307

   Solving temperature system... 0 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 7+0 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.2874e-16, 0.928307
      Relative nonlinear residual (total system) after nonlinear iteration 44: 0.928307

   Solving temperature system... 0 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 7+0 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.2874e-16, 0.928307
      Relative nonlinear residual (total system) after nonlinear iteration 45: 0.928307

   Solving temperature system... 0 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 7+0 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.2874e-16, 0.928307
      Relative nonlinear residual (total system) after nonlinear iteration 46: 0.928307

   Solving temperature system... 0 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 7+0 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.2874e-16, 0.928307
      Relative nonlinear residual (total system) after nonlinear iteration 47: 0.928307

   Solving temperature system... 0 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 7+0 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.2874e-16, 0.928307
      Relative nonlinear residual (total system) after nonlinear iteration 48: 0.928307

   Solving temperature system... 0 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 7+0 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.2874e-16, 0.928307
      Relative nonlinear residual (total system) after nonlinear iteration 49: 0.928307

   Solving temperature system... 0 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 7+0 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.2874e-16, 0.928307
      Relative nonlinear residual (total system) after nonlinear iteration 50: 0.928307


   Postprocessing:
     Writing graphical output: output-constant-visc/solution/solution-00002

Termination requested by criterion: end time


+----------------------------------------------+------------+------------+
| Total wallclock time elapsed since start     |      3.88s |            |
|                                              |            |            |
| Section                          | no. calls |  wall time | % of total |
+----------------------------------+-----------+------------+------------+
| Assemble Stokes system           |       150 |     0.848s |        22% |
| Assemble temperature system      |       150 |     0.865s |        22% |
| Build Stokes preconditioner      |       150 |      1.01s |        26% |
| Build temperature preconditioner |       150 |    0.0139s |      0.36% |
| Initialization                   |         1 |     0.501s |        13% |
| Postprocessing                   |         3 |    0.0178s |      0.46% |
| Setup dof systems                |         1 |   0.00346s |         0% |
| Setup initial conditions         |         1 |   0.00882s |      0.23% |
| Setup matrices                   |         1 |   0.00351s |         0% |
| Solve Stokes system              |       150 |     0.512s |        13% |
| Solve temperature system         |       150 |    0.0474s |       1.2% |
+----------------------------------+-----------+------------+------------+

-- Total wallclock time elapsed including restarts: 4s
-----------------------------------------------------------------------------
-- For information on how to cite ASPECT, see:
--   https://aspect.geodynamics.org/citing.html?ver=2.6.0-pre&sha=063741d02&src=code
-----------------------------------------------------------------------------
@jdannberg jdannberg added the bug label Aug 23, 2023
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Assignees
No one assigned
Labels
bug
Projects
None yet
Milestone
No milestone
Development

No branches or pull requests
1 participant
@jdannberg