Procedures

ProcedureLocationProcedure TypeDescription
add_lvst dsparsk.f Subroutine
amudia dsparsk.f Subroutine
amux dsparsk.f Subroutine
aplb dsparsk.f Subroutine
aplb1 dsparsk.f Subroutine
aplsb dsparsk.f Subroutine
aplsca dsparsk.f Subroutine
atob dsparsk.f Subroutine
banja daskr_banpre Subroutine

This subroutine generates a banded preconditioner matrix that approximates the iteration matrix , where the DAE system is . The band matrix has half-bandwidths and . It is computed by making calls to the user's res routine and forming difference quotients, exactly as in the banded direct method option of DASKR. banja calls the LINPACK routine dgbfa to do an LU factorization of this matrix.
banps daskr_banpre Subroutine

This subroutine uses the factors produced by banja to solve linear systems for the banded preconditioner , given a vector . It calls the LINPACK routine dgbsl for this.
BFS dsparsk.f Subroutine
c1_average web_m Subroutine

This routine computes the spatial average value of .
cinit web_m Subroutine

This routine computes and loads the vectors of initial values.
coicsr dsparsk.f Subroutine
coocsr dsparsk.f Subroutine
cperm dsparsk.f Subroutine
csrbnd dsparsk.f Subroutine
csrcoo dsparsk.f Subroutine
csrcsc dsparsk.f Subroutine
csrcsc2 dsparsk.f Subroutine
csrdia dsparsk.f Subroutine
csrdns dsparsk.f Subroutine
D1MACH daux.f Function
DASKR daskr.f Subroutine
DATV daskr.f Subroutine
daxpy dlinpk.f Subroutine
dblstr dsparsk.f Subroutine
DCNST0 daskr.f Subroutine
DCNSTR daskr.f Subroutine
dcopy dlinpk.f Subroutine
DDASIC daskr.f Subroutine
DDASID daskr.f Subroutine
DDASIK daskr.f Subroutine
DDATRP daskr.f Subroutine
DDAWTS daskr.f Subroutine
ddot dlinpk.f Function
DDSTP daskr.f Subroutine
DDWNRM daskr.f Function
DFNRMD daskr.f Subroutine
DFNRMK daskr.f Subroutine
dgbfa dlinpk.f Subroutine
dgbsl dlinpk.f Subroutine
dgefa dlinpk.f Subroutine
dgesl dlinpk.f Subroutine
DGSET2 drbgpre.f Subroutine
DHELS daskr.f Subroutine
DHEQR daskr.f Subroutine
diamua dsparsk.f Subroutine
diapos dsparsk.f Subroutine
DINVWT daskr.f Subroutine
DJACILU dilupre.f Subroutine
DJCALC dilupre.f Subroutine
DJILU dilupre.f Subroutine
DJREORD dilupre.f Subroutine
DLINSD daskr.f Subroutine
DLINSK daskr.f Subroutine
DMATD daskr.f Subroutine
DMSET2 drbdpre.f Subroutine
DNEDD daskr.f Subroutine
DNEDK daskr.f Subroutine
dnrm2 dlinpk.f Function
DNSD daskr.f Subroutine
DNSID daskr.f Subroutine
DNSIK daskr.f Subroutine
DNSK daskr.f Subroutine
DORTH daskr.f Subroutine
dperm dsparsk.f Subroutine
DPSOLILU dilupre.f Subroutine
DRBDJA drbdpre.f Subroutine
DRBDPS drbdpre.f Subroutine
DRBGJA drbgpre.f Subroutine
DRBGPS drbgpre.f Subroutine
DRCHEK daskr.f Subroutine
DROOTS daskr.f Subroutine
dscal dlinpk.f Subroutine
DSLVD daskr.f Subroutine
DSLVK daskr.f Subroutine
DSPIGM daskr.f Subroutine
DSPSETUP dilupre.f Subroutine
DUMSUM daux.f Subroutine
dvperm dsparsk.f Subroutine
DYYPNW daskr.f Subroutine
f krdem1_m Subroutine

dy1/dt routine.
f web_m Subroutine

This routine computes the right-hand sides of all the equations and returns them in the array cprime. The interaction rates are computed by calls to rates, and these are saved in rpar for later use in preconditioning.
f krdem2_m Subroutine

dy/dt routine.
gauss_seidel web_m Subroutine

This routine provides the inverse of the spatial factor for a product preconditoner in an ns-species reaction-diffusion problem. It performs itmax Gauss-Seidel iterations to compute an approximation to , where , and represents the diffusion contributions to the Jacobian. The solution vector is returned in z.
getbwd dsparsk.f Subroutine
GSET1 drbgpre.f Subroutine
idamax dlinpk.f Function
ilut dsparsk.f Subroutine
ilutp dsparsk.f Subroutine
infdia dsparsk.f Subroutine
ivperm dsparsk.f Subroutine
IXSAV daux.f Function
jac krdem2_m Subroutine

Jacobian routine.
jacrs web_m Subroutine

This routine interfaces to subroutines DRBDJA or DRBGJA, depending on the flag jbg=ipar(2), to generate and preprocess the block-diagonal Jacobian corresponding to the reaction term .

Read more…
lusol dsparsk.f Subroutine
maskdeg dsparsk.f Function
out web_m Subroutine

This routine prints the values of the individual species densities at the current time t, to logical unit lout.
perphn dsparsk.f Subroutine
prtmt dsparsk.f Subroutine
psolrs web_m Subroutine

This routine applies the inverse of a product preconditioner matrix to the vector in the array b. Depending on the flag jpre, this involves a call to gs, for the inverse of the spatial factor, and/or a call to DRBDPS or DRBGPS for the inverse of the reaction-based factor (cj*I_d - dR/dy). The latter factor uses block-grouping (with a call to DRBGPS) if jbg == 1, and does not (with a call to DRBDPS) if jbg == 0. the flag jbg is passed as ipar(2). The array b is overwritten with the solution.
qsplit dsparsk.f Subroutine
rates web_m Subroutine

This routine computes one block of the rate term of the system , namely block (jx, jy), for use in preconditioning.
res heatilu_m Subroutine

User-supplied residuals subroutine. It computes the residuals for the 2D discretized heat equation, with zero boundary values.
res heat_m Subroutine

User-supplied residuals subroutine. It computes the residuals for the 2D discretized heat equation, with zero boundary values.
res krdem1_m Subroutine

Residuals routine.
res web_m Subroutine

This routine computes the residuals vector, using subroutine f for the right-hand sides.
res krdem2_m Subroutine

Residuals routine.
rnrms dsparsk.f Subroutine
roscal dsparsk.f Subroutine
rperm dsparsk.f Subroutine
rt heatilu_m Subroutine

Roots routine.

rt heat_m Subroutine

Roots routine.
rt krdem1_m Subroutine

Roots routine.
rt web_m Subroutine

Roots routine.
rt krdem2_m Subroutine

Roots routine.
rversp dsparsk.f Subroutine
setid web_m Subroutine

This routine sets the ID array in iwork, indicating which components are differential and which are algebraic.
setpar web_par Subroutine

This routine sets the basic problem parameters which are passed to the various routines via the module web_par.
stripes dsparsk.f Subroutine
to_string web_m Function
uinit heatilu_m Subroutine

This routine computes and loads the vector of initial values. The initial u values are given by the polynomial u = 16x(1-x)y(1-y). The initial uprime values are set to zero (DASKR corrects these during the first time step).
uinit heat_m Subroutine

This routine computes and loads the vector of initial values. The initial u values are given by the polynomial u = 16x(1-x)y(1-y). The initial uprime values are set to zero (DASKR corrects these during the first time step).
XERRWD daux.f Subroutine
XSETF daux.f Subroutine
XSETUN daux.f Subroutine
call~~graph~~CallGraph proc~add_lvst add_lvst proc~amudia amudia proc~amux amux proc~aplb aplb w w proc~aplb->w proc~aplb1 aplb1 proc~aplsb aplsb proc~aplsca aplsca diapos diapos proc~aplsca->diapos proc~atob atob proc~banja banja d1mach d1mach proc~banja->d1mach dgbfa dgbfa proc~banja->dgbfa proc~banps banps dgbsl dgbsl proc~banps->dgbsl proc~bfs BFS add_lvst add_lvst proc~bfs->add_lvst proc~c1_average c1_average proc~cinit cinit proc~f~2 f proc~cinit->proc~f~2 proc~coicsr coicsr proc~coocsr coocsr proc~cperm cperm proc~csrbnd csrbnd bd bd proc~csrbnd->bd getbwd getbwd proc~csrbnd->getbwd proc~csrcoo csrcoo proc~csrcsc csrcsc csrcsc2 csrcsc2 proc~csrcsc->csrcsc2 proc~csrcsc2 csrcsc2 proc~csrdia csrdia infdia infdia proc~csrdia->infdia proc~csrdns csrdns f f proc~csrdns->f ns ns proc~csrdns->ns proc~d1mach D1MACH dumsum dumsum proc~d1mach->dumsum proc~daskr DASKR proc~daskr->d1mach dcnst0 dcnst0 proc~daskr->dcnst0 ddasic ddasic proc~daskr->ddasic ddatrp ddatrp proc~daskr->ddatrp ddawts ddawts proc~daskr->ddawts ddstp ddstp proc~daskr->ddstp ddwnrm ddwnrm proc~daskr->ddwnrm dinvwt dinvwt proc~daskr->dinvwt drchek drchek proc~daskr->drchek xerrwd xerrwd proc~daskr->xerrwd proc~datv DATV proc~daxpy daxpy proc~dblstr dblstr bfs bfs proc~dblstr->bfs perphn perphn proc~dblstr->perphn stripes stripes proc~dblstr->stripes proc~dcnst0 DCNST0 proc~dcnstr DCNSTR proc~dcopy dcopy proc~ddasic DDASIC dcopy dcopy proc~ddasic->dcopy proc~ddasid DDASID dmatd dmatd proc~ddasid->dmatd dnsid dnsid proc~ddasid->dnsid proc~ddasik DDASIK proc~ddasik->dcopy dnsik dnsik proc~ddasik->dnsik proc~ddatrp DDATRP proc~ddawts DDAWTS proc~ddot ddot proc~ddstp DDSTP proc~ddstp->ddatrp proc~ddstp->ddwnrm proc~ddwnrm DDWNRM proc~dfnrmd DFNRMD proc~dfnrmd->ddwnrm dslvd dslvd proc~dfnrmd->dslvd proc~dfnrmk DFNRMK proc~dfnrmk->dcopy proc~dfnrmk->ddwnrm dscal dscal proc~dfnrmk->dscal proc~dgbfa dgbfa daxpy daxpy proc~dgbfa->daxpy proc~dgbfa->dscal proc~dgbsl dgbsl proc~dgbsl->daxpy proc~dgefa dgefa proc~dgefa->daxpy proc~dgefa->dscal proc~dgesl dgesl proc~dgesl->daxpy proc~dgset2 DGSET2 proc~dgset2->d1mach gset1 gset1 proc~dgset2->gset1 proc~dhels DHELS proc~dhels->daxpy proc~dheqr DHEQR proc~diamua diamua proc~diapos diapos proc~dinvwt DINVWT proc~djacilu DJACILU amudia amudia proc~djacilu->amudia djcalc djcalc proc~djacilu->djcalc djilu djilu proc~djacilu->djilu djreord djreord proc~djacilu->djreord dvperm dvperm proc~djacilu->dvperm prtmt prtmt proc~djacilu->prtmt roscal roscal proc~djacilu->roscal proc~djacilu->xerrwd proc~djcalc DJCALC coocsr coocsr proc~djcalc->coocsr proc~djcalc->xerrwd proc~djilu DJILU ilut ilut proc~djilu->ilut ilutp ilutp proc~djilu->ilutp proc~djilu->xerrwd proc~djreord DJREORD atob atob proc~djreord->atob proc~djreord->bfs dperm dperm proc~djreord->dperm rversp rversp proc~djreord->rversp proc~dlinsd DLINSD dcnstr dcnstr proc~dlinsd->dcnstr proc~dlinsd->dcopy dfnrmd dfnrmd proc~dlinsd->dfnrmd dyypnw dyypnw proc~dlinsd->dyypnw proc~dlinsd->xerrwd proc~dlinsk DLINSK proc~dlinsk->dcnstr proc~dlinsk->dcopy dfnrmk dfnrmk proc~dlinsk->dfnrmk proc~dlinsk->dyypnw proc~dlinsk->xerrwd proc~dmatd DMATD proc~dmatd->dgbfa dgefa dgefa proc~dmatd->dgefa proc~dmset2 DMSET2 proc~dmset2->d1mach proc~dnedd DNEDD proc~dnedd->ddwnrm proc~dnedd->dmatd dnsd dnsd proc~dnedd->dnsd proc~dnedk DNEDK proc~dnedk->ddwnrm dnsk dnsk proc~dnedk->dnsk proc~dnrm2 dnrm2 float float proc~dnrm2->float proc~dnsd DNSD proc~dnsd->ddwnrm proc~dnsd->dslvd proc~dnsid DNSID proc~dnsid->dcopy proc~dnsid->ddwnrm dlinsd dlinsd proc~dnsid->dlinsd proc~dnsid->dslvd proc~dnsik DNSIK proc~dnsik->dcopy proc~dnsik->ddwnrm proc~dnsik->dfnrmk dlinsk dlinsk proc~dnsik->dlinsk dslvk dslvk proc~dnsik->dslvk proc~dnsk DNSK proc~dnsk->ddwnrm proc~dnsk->dslvk proc~dorth DORTH proc~dorth->daxpy proc~dperm dperm cperm cperm proc~dperm->cperm rperm rperm proc~dperm->rperm proc~dpsolilu DPSOLILU lusol lusol proc~dpsolilu->lusol proc~drbdja DRBDJA proc~drbdja->dgefa proc~drbdps DRBDPS dgesl dgesl proc~drbdps->dgesl proc~drbgja DRBGJA proc~drbgja->dgefa jxr jxr proc~drbgja->jxr jyr jyr proc~drbgja->jyr proc~drbgps DRBGPS proc~drbgps->dgesl jigx jigx proc~drbgps->jigx jigy jigy proc~drbgps->jigy proc~drchek DRCHEK proc~drchek->dcopy proc~drchek->ddatrp droots droots proc~drchek->droots proc~droots DROOTS proc~droots->dcopy proc~dscal dscal proc~dslvd DSLVD proc~dslvd->dgbsl proc~dslvd->dgesl proc~dslvk DSLVK proc~dslvk->dcopy proc~dslvk->dscal dspigm dspigm proc~dslvk->dspigm proc~dspigm DSPIGM datv datv proc~dspigm->datv proc~dspigm->daxpy proc~dspigm->dcopy dhels dhels proc~dspigm->dhels dheqr dheqr proc~dspigm->dheqr dorth dorth proc~dspigm->dorth proc~dspigm->dscal proc~dspsetup DSPSETUP proc~dumsum DUMSUM proc~dvperm dvperm proc~dyypnw DYYPNW proc~f f proc~rates rates proc~f~2->proc~rates proc~f~3 f proc~gauss_seidel gauss_seidel proc~getbwd getbwd proc~gset1 GSET1 proc~idamax idamax proc~ilut ilut qsplit qsplit proc~ilut->qsplit proc~ilutp ilutp proc~ilutp->qsplit proc~infdia infdia proc~ivperm ivperm proc~ixsav IXSAV proc~jac jac proc~jacrs jacrs drbdja drbdja proc~jacrs->drbdja drbgja drbgja proc~jacrs->drbgja proc~lusol lusol proc~maskdeg maskdeg proc~out out proc~to_string to_string proc~out->proc~to_string proc~perphn perphn proc~perphn->bfs proc~prtmt prtmt 1h 1h proc~prtmt->1h alog10 alog10 proc~prtmt->alog10 proc~prtmt->f lseif lseif proc~prtmt->lseif proc~psolrs psolrs proc~psolrs->proc~gauss_seidel drbdps drbdps proc~psolrs->drbdps drbgps drbgps proc~psolrs->drbgps proc~qsplit qsplit proc~res res proc~res~2 res proc~res~3 res proc~res~3->proc~f proc~res~4 res proc~res~4->proc~f~2 proc~res~5 res proc~res~5->proc~f~3 proc~rnrms rnrms proc~roscal roscal diamua diamua proc~roscal->diamua rnrms rnrms proc~roscal->rnrms proc~rperm rperm proc~rt rt proc~rt~2 rt proc~rt~3 rt proc~rt~4 rt proc~rt~4->proc~c1_average proc~rt~5 rt proc~rversp rversp proc~setid setid proc~setpar setpar proc~stripes stripes proc~uinit uinit proc~uinit~2 uinit proc~xerrwd XERRWD proc~xsetf XSETF proc~xsetun XSETUN program~example_heat example_heat program~example_heat->proc~uinit~2 daskr daskr program~example_heat->daskr program~example_heatilu example_heatilu program~example_heatilu->proc~uinit program~example_heatilu->daskr dspsetup dspsetup program~example_heatilu->dspsetup program~example_web example_web program~example_web->proc~c1_average program~example_web->proc~cinit program~example_web->proc~out program~example_web->proc~setid program~example_web->proc~setpar program~example_web->daskr dgset2 dgset2 program~example_web->dgset2 dmset2 dmset2 program~example_web->dmset2 program~test_krdem1 test_krdem1 program~test_krdem1->daskr program~test_krdem2 test_krdem2 program~test_krdem2->daskr
Help