5 DOUBLE PRECISION FUNCTION ddwnrm(NEQ,V,RWT,RPAR,IPAR)
20 IMPLICIT DOUBLE PRECISION(
a-
h,o-z)
26 IF(
abs(v(i)*rwt(i)) .GT. vmax) vmax =
abs(v(i)*rwt(i))
28 IF(vmax .LE. 0.0d0) go
to 30
31 20 sum = sum + ((v(i)*rwt(i))/vmax)**2
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
may be zero for pure relative error test tem the relative tolerance must be greater than or equal to
OCTAVE_EXPORT octave_value_list etc The functions then dimension(columns)
OCTAVE_EXPORT octave_value_list return the value of the option it must match the dimension of the state and the relative tolerance must also be a vector of the same length tem it must match the dimension of the state and the absolute tolerance must also be a vector of the same length The local error test applied at each integration step is xample roup abs(local error in x(i))<
double precision function ddwnrm(NEQ, V, RWT, RPAR, IPAR)
octave_value sqrt(void) const