daskr Module


Uses

  • module~~daskr~~UsesGraph module~daskr daskr module~daskr_kinds daskr_kinds module~daskr->module~daskr_kinds iso_fortran_env iso_fortran_env module~daskr_kinds->iso_fortran_env

Used by

  • module~~daskr~~UsedByGraph module~daskr daskr module~daskr_banpre daskr_banpre module~daskr_banpre->module~daskr module~daskr_ilupre daskr_ilupre module~daskr_ilupre->module~daskr proc~datv datv proc~datv->module~daskr proc~ddasid ddasid proc~ddasid->module~daskr proc~ddasik ddasik proc~ddasik->module~daskr proc~dfnrmd dfnrmd proc~dfnrmd->module~daskr proc~dfnrmk dfnrmk proc~dfnrmk->module~daskr proc~dlinsd dlinsd proc~dlinsd->module~daskr proc~dlinsk dlinsk proc~dlinsk->module~daskr proc~dmatd dmatd proc~dmatd->module~daskr proc~dnedd dnedd proc~dnedd->module~daskr proc~dnedk dnedk proc~dnedk->module~daskr proc~dnsd dnsd proc~dnsd->module~daskr proc~dnsid dnsid proc~dnsid->module~daskr proc~dnsik dnsik proc~dnsik->module~daskr proc~dnsk dnsk proc~dnsk->module~daskr proc~dslvd dslvd proc~dslvd->module~daskr proc~dslvk dslvk proc~dslvk->module~daskr proc~dspigm dspigm proc~dspigm->module~daskr proc~jac~2 jac proc~jac~2->module~daskr program~example_heat example_heat program~example_heat->module~daskr_banpre program~example_heatilu example_heatilu program~example_heatilu->module~daskr_ilupre

Variables

Type Visibility Attributes Name Initial
integer, public, parameter :: rwloc_cj = 1

cj: scalar used in forming the system Jacobian

integer, public, parameter :: rwloc_cjold = 2

cjold: scalar used in forming the system Jacobian on last step

integer, public, parameter :: rwloc_h = 3

h: step size for next step

integer, public, parameter :: rwloc_tn = 4

tn: current independent variable value

integer, public, parameter :: rwloc_hold = 7

hold: step size used on last step

integer, public, parameter :: iwloc_ml = 1

ml: lower bandwidth

integer, public, parameter :: iwloc_mu = 2

mu: upper bandwidth

integer, public, parameter :: iwloc_mtype = 4

mtype: method type

integer, public, parameter :: iwloc_k = 7

k: order of method for next step

integer, public, parameter :: iwloc_kold = 8

kold: order used on last step

integer, public, parameter :: iwloc_nst = 11

nst: number of steps taken

integer, public, parameter :: iwloc_nre = 12

nre: number of res calls

integer, public, parameter :: iwloc_nje = 13

nje: number of jac calls

integer, public, parameter :: iwloc_netf = 14

netf: number of error test failures

integer, public, parameter :: iwloc_ncfn = 15

ncfn: number of nonlinear convergence failures

integer, public, parameter :: iwloc_ncfl = 16

ncfl: number of linear iteration convergence failures

integer, public, parameter :: iwloc_leniw = 17

leniw: actual iwork length used

integer, public, parameter :: iwloc_lenrw = 18

lenrw: actual rwork length used

integer, public, parameter :: iwloc_nni = 19

nni: number of nonlinear iterations

integer, public, parameter :: iwloc_nli = 20

nli: number of linear (Krylov) iterations

integer, public, parameter :: iwloc_nps = 21

nps: number of psol calls

integer, public, parameter :: iwloc_npd = 22

npd: length of partial derivatives vector returned by jac

integer, public, parameter :: iwloc_miter = 23

miter: ??

integer, public, parameter :: iwloc_maxl = 24

maxl: ??

integer, public, parameter :: iwloc_kmp = 25

kmp: ??

integer, public, parameter :: iwloc_nrmax = 26

nrmax: ??

integer, public, parameter :: iwloc_lrwp = 29

lrwp: location of start of rwp in rwork.

integer, public, parameter :: iwloc_liwp = 30

liwp: location of start of iwp in iwork.

integer, public, parameter :: iwloc_kprint = 31

kprint: flag to turn on print

integer, public, parameter :: iwloc_mxnit = 32

mxnit: maximum number of Newton iterations

integer, public, parameter :: iwloc_mxnj = 33

mxnj: maximum number of Jacobian evaluations ?

integer, public, parameter :: iwloc_lsoff = 35

lsoff: flag to turn off linesearch

integer, public, parameter :: iwloc_nrtfn = 36

nrtfn: number of rt calls


Abstract Interfaces

abstract interface

  • public subroutine res_t(t, y, ydot, cj, delta, ires, rpar, ipar)

    Arguments

    Type IntentOptional Attributes Name
    real(kind=rk), intent(in) :: t
    real(kind=rk), intent(in) :: y(*)
    real(kind=rk), intent(in) :: ydot(*)
    real(kind=rk), intent(in) :: cj
    real(kind=rk), intent(out) :: delta(*)
    integer, intent(inout) :: ires
    real(kind=rk), intent(inout) :: rpar(*)
    integer, intent(inout) :: ipar(*)

abstract interface

  • public subroutine jac_t(t, y, ydot, pd, cj, rpar, ipar)

    Arguments

    Type IntentOptional Attributes Name
    real(kind=rk), intent(in) :: t
    real(kind=rk), intent(in) :: y(*)
    real(kind=rk), intent(in) :: ydot(*)
    real(kind=rk), intent(out) :: pd(*)
    real(kind=rk), intent(in) :: cj
    real(kind=rk), intent(inout) :: rpar(*)
    integer, intent(inout) :: ipar(*)

abstract interface

  • public subroutine jack_t(res, ires, neq, t, y, ydot, rewt, savr, wk, h, cj, rwp, iwp, ierr, rpar, ipar)

    Arguments

    Type IntentOptional Attributes Name
    procedure(res_t) :: res
    integer, intent(out) :: ires
    integer, intent(in) :: neq
    real(kind=rk), intent(in) :: t
    real(kind=rk), intent(inout) :: y(*)
    real(kind=rk), intent(inout) :: ydot(*)
    real(kind=rk), intent(in) :: rewt(*)
    real(kind=rk), intent(inout) :: savr(*)
    real(kind=rk), intent(inout) :: wk(*)
    real(kind=rk), intent(in) :: h
    real(kind=rk), intent(in) :: cj
    real(kind=rk), intent(inout) :: rwp(*)
    integer, intent(inout) :: iwp(*)
    integer, intent(inout) :: ierr
    real(kind=rk), intent(inout) :: rpar(*)
    integer, intent(inout) :: ipar(*)

abstract interface

  • public subroutine psol_t(neq, t, y, ydot, savr, wk, cj, wght, rwp, iwp, b, epslin, ierr, rpar, ipar)

    Arguments

    Type IntentOptional Attributes Name
    integer, intent(in) :: neq
    real(kind=rk), intent(in) :: t
    real(kind=rk), intent(in) :: y(*)
    real(kind=rk), intent(in) :: ydot(*)
    real(kind=rk), intent(in) :: savr(*)
    real(kind=rk), intent(inout) :: wk(*)
    real(kind=rk), intent(in) :: cj
    real(kind=rk), intent(in) :: wght(*)
    real(kind=rk), intent(inout) :: rwp(*)
    integer, intent(inout) :: iwp(*)
    real(kind=rk), intent(inout) :: b(*)
    real(kind=rk), intent(in) :: epslin
    integer, intent(inout) :: ierr
    real(kind=rk), intent(inout) :: rpar(*)
    integer, intent(inout) :: ipar(*)