Mercurial > octave-nkf

comparison scripts/ode/ode45.m @ 20646:01586012300e

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
compute number of fevals correctly in ode45 (bug #46189) * scripts/ode/ode45.m: account for FSAL when counting FEvals
author Carlo de Falco <carlo.defalco@polimi.it>
date Mon, 12 Oct 2015 22:50:59 +0200
parents a22d8a2eb0e5
children 05c77bc1d204
comparison
equal deleted inserted replaced
20645:395140e53656 20646:01586012300e
406 ## Print additional information if option Stats is set 406 ## Print additional information if option Stats is set
407 if (strcmp (vodeoptions.Stats, "on")) 407 if (strcmp (vodeoptions.Stats, "on"))
408 vhavestats = true; 408 vhavestats = true;
409 vnsteps = solution.vcntloop-2; # vcntloop from 2..end 409 vnsteps = solution.vcntloop-2; # vcntloop from 2..end
410 vnfailed = (solution.vcntcycles-1)-(vnsteps)+1; # vcntcycl from 1..end 410 vnfailed = (solution.vcntcycles-1)-(vnsteps)+1; # vcntcycl from 1..end
411 vnfevals = 7*(solution.vcntcycles-1); # number of ode evaluations 411 vnfevals = 7 * (solution.vcntcycles-1) + 1; # number of ode evaluations
412 vndecomps = 0; # number of LU decompositions 412 vndecomps = 0; # number of LU decompositions
413 vnpds = 0; # number of partial derivatives 413 vnpds = 0; # number of partial derivatives
414 vnlinsols = 0; # no. of linear systems solutions 414 vnlinsols = 0; # no. of linear systems solutions
415 ## Print cost statistics if no output argument is given 415 ## Print cost statistics if no output argument is given
416 if (nargout == 0) 416 if (nargout == 0)