5 SUBROUTINE dfnrmd (NEQ, Y, T, YPRIME, R, CJ, WT, RES, IRES,
6 * fnorm, wm, iwm, rpar, ipar)
36 IMPLICIT DOUBLE PRECISION (
a-
h,o-z)
44 CALL res(t,y,yprime,cj,r,ires,rpar,ipar)
45 IF (ires .LT. 0)
RETURN
49 CALL
dslvd(neq,r,wm,iwm)
53 fnorm =
ddwnrm(neq,r,wt,rpar,ipar)
calling an anonymous function involves an overhead quite comparable to the overhead of an m file function Passing a handle to a built in function is because the interpreter is not involved in the internal loop For a
subroutine dslvd(NEQ, DELTA, WM, IWM)
OCTAVE_EXPORT octave_value_list etc The functions then dimension(columns)
subroutine dfnrmd(NEQ, Y, T, YPRIME, R, CJ, WT, RES, IRES, FNORM, WM, IWM, RPAR, IPAR)
double precision function ddwnrm(NEQ, V, RWT, RPAR, IPAR)