timestep_stiff

Module

Description

$Id$ adapted from timestep_rkf, and from numerical recipe stiff algorithm ** AUTOMATIC CPARAM.INC GENERATION ************************ CPARAM logical, parameter :: lcourant_dt = .false. ***********************************************************

Quick access

Variables:

a21, a31, a32, advance_substeps, after_substep_sum_time, b1, b2, b3, b4, beta_hat, c21, c31, c32, c41, c42, c43, calc_pencils_sub_cycle, dt_alpha_ts, dt_beta_hat, dt_decrease, dt_increase, dt_major, e1, e2, e3, e4, errcon, errdf, errmax, errmaxs, f1, farr, farraymin, gam, grow, itter, maxtry, ode, ode_timestep_first, ode_timestep_second, pde_energy_only, rk3, rkck, rkck3, rkl_coeff, safety, shrnk, split_update, stiff, substeps, swap, update_after_substep

Routines:

initialize_timestep(), time_step()

Needed modules

Variables

  • timestep/a21 [real,private/parameter/optional/default=2.0]
  • timestep/a31 [real,private/parameter/optional/default=1.92]
  • timestep/a32 [real,private/parameter/optional/default=0.24]
  • timestep/after_substep_sum_time [real,public/optional/default=0.0]
  • timestep/b1 [real,private/parameter/optional/default=2.111111111111111]
  • timestep/b2 [real,private/parameter/optional/default=0.5]
  • timestep/b3 [real,private/parameter/optional/default=0.23148148148148148]
  • timestep/b4 [real,private/parameter/optional/default=1.1574074074074074]
  • timestep/c21 [real,private/parameter/optional/default=-8.0]
  • timestep/c31 [real,private/parameter/optional/default=14.88]
  • timestep/c32 [real,private/parameter/optional/default=2.4]
  • timestep/c41 [real,private/parameter/optional/default=-0.896]
  • timestep/c42 [real,private/parameter/optional/default=-0.432]
  • timestep/c43 [real,private/parameter/optional/default=-0.4]
  • timestep/dt_decrease [real,private/parameter/optional/default=-0.3333333333333333]
  • timestep/dt_increase [real,private/parameter/optional/default=-0.25]
  • timestep/e1 [real,private/parameter/optional/default=0.3148148148148148]
  • timestep/e2 [real,private/parameter/optional/default=0.19444444444444445]
  • timestep/e3 [real,private/parameter/optional/default=0.0]
  • timestep/e4 [real,private/parameter/optional/default=1.1574074074074074]
  • timestep/errcon [real,private/parameter/optional/default=0.1296]
  • timestep/farraymin (mvar) [real,private]
  • timestep/gam [real,private/parameter/optional/default=0.5]
  • timestep/grow [real,private/parameter/optional/default=1.5]
  • timestep/maxtry [integer,private/parameter/optional/default=40]
  • timestep/safety [real,private/parameter/optional/default=0.9]
  • timestep/shrnk [real,private/parameter/optional/default=0.5]
  • timestep/stiff [private]

Subroutines and functions

subroutine  timestep/initialize_timestep()
Call to:

warning(), fatal_error(), jacobn(), ludcmp(), pde(), lubksb()

subroutine  timestep/time_step(f, df, p)

Stiff solver, accurate to 2nd(?) order

28-aug-09/rplasson: coded

Parameters:

  • f (,,*,*) [real,contiguous]

  • df (,,*,*) [real,contiguous]

  • p [pencil_case]

Use :

messages

Call to:

warning(), fatal_error(), jacobn(), ludcmp(), pde(), lubksb()