5 SUBROUTINE datv (NEQ, Y, TN, YPRIME, SAVR, V, WGHT, YPTEM, RES,
6 * ires, psol, z, vtem, wp, iwp, cj, eplin, ier, nre, npsl,
78 INTEGER NEQ, IRES, IWP, IER, NRE, NPSL, IPAR
79 DOUBLE PRECISION Y, TN, YPRIME, SAVR, V, WGHT, YPTEM, Z, VTEM,
81 dimension y(*), yprime(*), savr(*), v(*), wght(*), yptem(*),
82 1 z(*), vtem(*), wp(*), iwp(*), rpar(*), ipar(*)
84 DOUBLE PRECISION EPLIN
92 10 vtem(i) = v(i)/wght(i)
99 yptem(i) = yprime(i) + vtem(i)*cj
100 20 z(i) = y(i) + vtem(i)
106 CALL res(tn,z,yptem,cj,vtem,ires,rpar,ipar)
108 IF (ires .LT. 0)
RETURN
114 70 z(i) = vtem(i) - savr(i)
118 CALL psol(neq, tn, y, yprime, savr, yptem, cj, wght, wp, iwp,
119 1 z, eplin, ier, rpar, ipar)
121 IF (ier .NE. 0)
RETURN
126 90 z(i) = z(i)*wght(i)
std::string dimension(void) const
subroutine datv(NEQ, Y, TN, YPRIME, SAVR, V, WGHT, YPTEM, RES, IRES, PSOL, Z, VTEM, WP, IWP, CJ, EPLIN, IER, NRE, NPSL, RPAR, IPAR)