octave: scripts/ode/ode45.m annotate

annotate scripts/ode/ode45.m @ 30564:796f54d4ddbf stable

update Octave Project Developers copyright for the new year In files that have the "Octave Project Developers" copyright notice, update for 2021. In all .txi and .texi files except gpl.txi and gpl.texi in the doc/liboctave and doc/interpreter directories, change the copyright to "Octave Project Developers", the same as used for other source files. Update copyright notices for 2022 (not done since 2019). For gpl.txi and gpl.texi, change the copyright notice to be "Free Software Foundation, Inc." and leave the date at 2007 only because this file only contains the text of the GPL, not anything created by the Octave Project Developers. Add Paul Thomas to contributors.in.

author	John W. Eaton <jwe@octave.org>
date	Tue, 28 Dec 2021 18:22:40 -0500
parents	7854d5752dd2
children	e1788b1a315f

rev	line source
27923 bd51beb6205e update formatting of copyright notices John W. Eaton <jwe@octave.org> parents: 27919 diff changeset	1 ########################################################################
bd51beb6205e update formatting of copyright notices John W. Eaton <jwe@octave.org> parents: 27919 diff changeset	2 ##
30564 796f54d4ddbf update Octave Project Developers copyright for the new year John W. Eaton <jwe@octave.org> parents: 29359 diff changeset	3 ## Copyright (C) 2006-2022 The Octave Project Developers
27918 b442ec6dda5c use centralized file for copyright info for individual contributors John W. Eaton <jwe@octave.org> parents: 26376 diff changeset	4 ##
27923 bd51beb6205e update formatting of copyright notices John W. Eaton <jwe@octave.org> parents: 27919 diff changeset	5 ## See the file COPYRIGHT.md in the top-level directory of this
bd51beb6205e update formatting of copyright notices John W. Eaton <jwe@octave.org> parents: 27919 diff changeset	6 ## distribution or <https://octave.org/copyright/>.
20533 fcb792acab9b Moving ode45, odeset, odeget, and levenshtein from odepkg to core. jcorno <jacopo.corno@gmail.com> parents: diff changeset	7 ##
fcb792acab9b Moving ode45, odeset, odeget, and levenshtein from odepkg to core. jcorno <jacopo.corno@gmail.com> parents: diff changeset	8 ## This file is part of Octave.
fcb792acab9b Moving ode45, odeset, odeget, and levenshtein from odepkg to core. jcorno <jacopo.corno@gmail.com> parents: diff changeset	9 ##
24534 194eb4bd202b maint: Update punctuation for GPL v3 license text. Rik <rik@octave.org> parents: 23220 diff changeset	10 ## Octave is free software: you can redistribute it and/or modify it
20533 fcb792acab9b Moving ode45, odeset, odeget, and levenshtein from odepkg to core. jcorno <jacopo.corno@gmail.com> parents: diff changeset	11 ## under the terms of the GNU General Public License as published by
24534 194eb4bd202b maint: Update punctuation for GPL v3 license text. Rik <rik@octave.org> parents: 23220 diff changeset	12 ## the Free Software Foundation, either version 3 of the License, or
22755 3a2b891d0b33 maint: Standardize Copyright formatting. Rik <rik@octave.org> parents: 22660 diff changeset	13 ## (at your option) any later version.
20533 fcb792acab9b Moving ode45, odeset, odeget, and levenshtein from odepkg to core. jcorno <jacopo.corno@gmail.com> parents: diff changeset	14 ##
fcb792acab9b Moving ode45, odeset, odeget, and levenshtein from odepkg to core. jcorno <jacopo.corno@gmail.com> parents: diff changeset	15 ## Octave is distributed in the hope that it will be useful, but
fcb792acab9b Moving ode45, odeset, odeget, and levenshtein from odepkg to core. jcorno <jacopo.corno@gmail.com> parents: diff changeset	16 ## WITHOUT ANY WARRANTY; without even the implied warranty of
22755 3a2b891d0b33 maint: Standardize Copyright formatting. Rik <rik@octave.org> parents: 22660 diff changeset	17 ## MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
3a2b891d0b33 maint: Standardize Copyright formatting. Rik <rik@octave.org> parents: 22660 diff changeset	18 ## GNU General Public License for more details.
20533 fcb792acab9b Moving ode45, odeset, odeget, and levenshtein from odepkg to core. jcorno <jacopo.corno@gmail.com> parents: diff changeset	19 ##
fcb792acab9b Moving ode45, odeset, odeget, and levenshtein from odepkg to core. jcorno <jacopo.corno@gmail.com> parents: diff changeset	20 ## You should have received a copy of the GNU General Public License
fcb792acab9b Moving ode45, odeset, odeget, and levenshtein from odepkg to core. jcorno <jacopo.corno@gmail.com> parents: diff changeset	21 ## along with Octave; see the file COPYING. If not, see
24534 194eb4bd202b maint: Update punctuation for GPL v3 license text. Rik <rik@octave.org> parents: 23220 diff changeset	22 ## <https://www.gnu.org/licenses/>.
27923 bd51beb6205e update formatting of copyright notices John W. Eaton <jwe@octave.org> parents: 27919 diff changeset	23 ##
bd51beb6205e update formatting of copyright notices John W. Eaton <jwe@octave.org> parents: 27919 diff changeset	24 ########################################################################
20533 fcb792acab9b Moving ode45, odeset, odeget, and levenshtein from odepkg to core. jcorno <jacopo.corno@gmail.com> parents: diff changeset	25
fcb792acab9b Moving ode45, odeset, odeget, and levenshtein from odepkg to core. jcorno <jacopo.corno@gmail.com> parents: diff changeset	26 ## -- texinfo --
20852 516bb87ea72e 2015 Code Sprint: remove class of function from docstring for all m-files. Rik <rik@octave.org> parents: 20719 diff changeset	27 ## @deftypefn {} {[@var{t}, @var{y}] =} ode45 (@var{fun}, @var{trange}, @var{init})
516bb87ea72e 2015 Code Sprint: remove class of function from docstring for all m-files. Rik <rik@octave.org> parents: 20719 diff changeset	28 ## @deftypefnx {} {[@var{t}, @var{y}] =} ode45 (@var{fun}, @var{trange}, @var{init}, @var{ode_opt})
516bb87ea72e 2015 Code Sprint: remove class of function from docstring for all m-files. Rik <rik@octave.org> parents: 20719 diff changeset	29 ## @deftypefnx {} {[@var{t}, @var{y}, @var{te}, @var{ye}, @var{ie}] =} ode45 (@dots{})
516bb87ea72e 2015 Code Sprint: remove class of function from docstring for all m-files. Rik <rik@octave.org> parents: 20719 diff changeset	30 ## @deftypefnx {} {@var{solution} =} ode45 (@dots{})
22939 54302b670139 doc: Document nargout=0 calling mode for ode23, ode45. Rik <rik@octave.org> parents: 22937 diff changeset	31 ## @deftypefnx {} {} ode45 (@dots{})
20533 fcb792acab9b Moving ode45, odeset, odeget, and levenshtein from odepkg to core. jcorno <jacopo.corno@gmail.com> parents: diff changeset	32 ##
20548 25623ef2ff4f doc: Rewrite docstrings for ode* family of functions. Rik <rik@octave.org> parents: 20543 diff changeset	33 ## Solve a set of non-stiff Ordinary Differential Equations (non-stiff ODEs)
20719 d903cccb8de8 doc: Periodic spellcheck of documentation. Rik <rik@octave.org> parents: 20634 diff changeset	34 ## with the well known explicit @nospell{Dormand-Prince} method of order 4.
20533 fcb792acab9b Moving ode45, odeset, odeget, and levenshtein from odepkg to core. jcorno <jacopo.corno@gmail.com> parents: diff changeset	35 ##
20634 80e630b37ba1 maint: Remove unnecessary 'v' prefix before variables in ODE m-files. Rik <rik@octave.org> parents: 20631 diff changeset	36 ## @var{fun} is a function handle, inline function, or string containing the
80e630b37ba1 maint: Remove unnecessary 'v' prefix before variables in ODE m-files. Rik <rik@octave.org> parents: 20631 diff changeset	37 ## name of the function that defines the ODE: @code{y' = f(t,y)}. The function
80e630b37ba1 maint: Remove unnecessary 'v' prefix before variables in ODE m-files. Rik <rik@octave.org> parents: 20631 diff changeset	38 ## must accept two inputs where the first is time @var{t} and the second is a
80e630b37ba1 maint: Remove unnecessary 'v' prefix before variables in ODE m-files. Rik <rik@octave.org> parents: 20631 diff changeset	39 ## column vector of unknowns @var{y}.
20533 fcb792acab9b Moving ode45, odeset, odeget, and levenshtein from odepkg to core. jcorno <jacopo.corno@gmail.com> parents: diff changeset	40 ##
20548 25623ef2ff4f doc: Rewrite docstrings for ode* family of functions. Rik <rik@octave.org> parents: 20543 diff changeset	41 ## @var{trange} specifies the time interval over which the ODE will be
20634 80e630b37ba1 maint: Remove unnecessary 'v' prefix before variables in ODE m-files. Rik <rik@octave.org> parents: 20631 diff changeset	42 ## evaluated. Typically, it is a two-element vector specifying the initial and
20548 25623ef2ff4f doc: Rewrite docstrings for ode* family of functions. Rik <rik@octave.org> parents: 20543 diff changeset	43 ## final times (@code{[tinit, tfinal]}). If there are more than two elements
25623ef2ff4f doc: Rewrite docstrings for ode* family of functions. Rik <rik@octave.org> parents: 20543 diff changeset	44 ## then the solution will also be evaluated at these intermediate time
22594 b8d525710075 cleanup ode solvers Carlo de Falco <carlo.defalco@polimi.it> parents: 22593 diff changeset	45 ## instances.
20533 fcb792acab9b Moving ode45, odeset, odeget, and levenshtein from odepkg to core. jcorno <jacopo.corno@gmail.com> parents: diff changeset	46 ##
20634 80e630b37ba1 maint: Remove unnecessary 'v' prefix before variables in ODE m-files. Rik <rik@octave.org> parents: 20631 diff changeset	47 ## By default, @code{ode45} uses an adaptive timestep with the
80e630b37ba1 maint: Remove unnecessary 'v' prefix before variables in ODE m-files. Rik <rik@octave.org> parents: 20631 diff changeset	48 ## @code{integrate_adaptive} algorithm. The tolerance for the timestep
25026 f886561f9696 doc: improve differential eqtn docs and mention ode15i/ode15s (bug #51965). Colin Macdonald <cbm@m.fsf.org> parents: 24534 diff changeset	49 ## computation may be changed by using the options @qcode{"RelTol"} and
f886561f9696 doc: improve differential eqtn docs and mention ode15i/ode15s (bug #51965). Colin Macdonald <cbm@m.fsf.org> parents: 24534 diff changeset	50 ## @qcode{"AbsTol"}.
20533 fcb792acab9b Moving ode45, odeset, odeget, and levenshtein from odepkg to core. jcorno <jacopo.corno@gmail.com> parents: diff changeset	51 ##
20634 80e630b37ba1 maint: Remove unnecessary 'v' prefix before variables in ODE m-files. Rik <rik@octave.org> parents: 20631 diff changeset	52 ## @var{init} contains the initial value for the unknowns. If it is a row
80e630b37ba1 maint: Remove unnecessary 'v' prefix before variables in ODE m-files. Rik <rik@octave.org> parents: 20631 diff changeset	53 ## vector then the solution @var{y} will be a matrix in which each column is
80e630b37ba1 maint: Remove unnecessary 'v' prefix before variables in ODE m-files. Rik <rik@octave.org> parents: 20631 diff changeset	54 ## the solution for the corresponding initial value in @var{init}.
20548 25623ef2ff4f doc: Rewrite docstrings for ode* family of functions. Rik <rik@octave.org> parents: 20543 diff changeset	55 ##
20634 80e630b37ba1 maint: Remove unnecessary 'v' prefix before variables in ODE m-files. Rik <rik@octave.org> parents: 20631 diff changeset	56 ## The optional fourth argument @var{ode_opt} specifies non-default options to
80e630b37ba1 maint: Remove unnecessary 'v' prefix before variables in ODE m-files. Rik <rik@octave.org> parents: 20631 diff changeset	57 ## the ODE solver. It is a structure generated by @code{odeset}.
80e630b37ba1 maint: Remove unnecessary 'v' prefix before variables in ODE m-files. Rik <rik@octave.org> parents: 20631 diff changeset	58 ##
80e630b37ba1 maint: Remove unnecessary 'v' prefix before variables in ODE m-files. Rik <rik@octave.org> parents: 20631 diff changeset	59 ## The function typically returns two outputs. Variable @var{t} is a
20548 25623ef2ff4f doc: Rewrite docstrings for ode* family of functions. Rik <rik@octave.org> parents: 20543 diff changeset	60 ## column vector and contains the times where the solution was found. The
25623ef2ff4f doc: Rewrite docstrings for ode* family of functions. Rik <rik@octave.org> parents: 20543 diff changeset	61 ## output @var{y} is a matrix in which each column refers to a different
25623ef2ff4f doc: Rewrite docstrings for ode* family of functions. Rik <rik@octave.org> parents: 20543 diff changeset	62 ## unknown of the problem and each row corresponds to a time in @var{t}.
20533 fcb792acab9b Moving ode45, odeset, odeget, and levenshtein from odepkg to core. jcorno <jacopo.corno@gmail.com> parents: diff changeset	63 ##
25026 f886561f9696 doc: improve differential eqtn docs and mention ode15i/ode15s (bug #51965). Colin Macdonald <cbm@m.fsf.org> parents: 24534 diff changeset	64 ## The output can also be returned as a structure @var{solution} which has a
f886561f9696 doc: improve differential eqtn docs and mention ode15i/ode15s (bug #51965). Colin Macdonald <cbm@m.fsf.org> parents: 24534 diff changeset	65 ## field @var{x} containing a row vector of times where the solution was
f886561f9696 doc: improve differential eqtn docs and mention ode15i/ode15s (bug #51965). Colin Macdonald <cbm@m.fsf.org> parents: 24534 diff changeset	66 ## evaluated and a field @var{y} containing the solution matrix such that each
f886561f9696 doc: improve differential eqtn docs and mention ode15i/ode15s (bug #51965). Colin Macdonald <cbm@m.fsf.org> parents: 24534 diff changeset	67 ## column corresponds to a time in @var{x}. Use
f886561f9696 doc: improve differential eqtn docs and mention ode15i/ode15s (bug #51965). Colin Macdonald <cbm@m.fsf.org> parents: 24534 diff changeset	68 ## @w{@code{fieldnames (@var{solution})}} to see the other fields and
21546 f7f97d7e9294 doc: Wrap m-file docstrings to 79 characters + newline (80 total). Rik <rik@octave.org> parents: 21443 diff changeset	69 ## additional information returned.
20634 80e630b37ba1 maint: Remove unnecessary 'v' prefix before variables in ODE m-files. Rik <rik@octave.org> parents: 20631 diff changeset	70 ##
28713 28d2511f2af2 maint: grammarcheck documentation ahead of 6.1 release. Rik <rik@octave.org> parents: 28045 diff changeset	71 ## If no output arguments are requested, and no @qcode{"OutputFcn"} is
28d2511f2af2 maint: grammarcheck documentation ahead of 6.1 release. Rik <rik@octave.org> parents: 28045 diff changeset	72 ## specified in @var{ode_opt}, then the @qcode{"OutputFcn"} is set to
28d2511f2af2 maint: grammarcheck documentation ahead of 6.1 release. Rik <rik@octave.org> parents: 28045 diff changeset	73 ## @code{odeplot} and the results of the solver are plotted immediately.
22939 54302b670139 doc: Document nargout=0 calling mode for ode23, ode45. Rik <rik@octave.org> parents: 22937 diff changeset	74 ##
25026 f886561f9696 doc: improve differential eqtn docs and mention ode15i/ode15s (bug #51965). Colin Macdonald <cbm@m.fsf.org> parents: 24534 diff changeset	75 ## If using the @qcode{"Events"} option then three additional outputs may be
f886561f9696 doc: improve differential eqtn docs and mention ode15i/ode15s (bug #51965). Colin Macdonald <cbm@m.fsf.org> parents: 24534 diff changeset	76 ## returned. @var{te} holds the time when an Event function returned a zero.
f886561f9696 doc: improve differential eqtn docs and mention ode15i/ode15s (bug #51965). Colin Macdonald <cbm@m.fsf.org> parents: 24534 diff changeset	77 ## @var{ye} holds the value of the solution at time @var{te}. @var{ie}
20634 80e630b37ba1 maint: Remove unnecessary 'v' prefix before variables in ODE m-files. Rik <rik@octave.org> parents: 20631 diff changeset	78 ## contains an index indicating which Event function was triggered in the case
80e630b37ba1 maint: Remove unnecessary 'v' prefix before variables in ODE m-files. Rik <rik@octave.org> parents: 20631 diff changeset	79 ## of multiple Event functions.
80e630b37ba1 maint: Remove unnecessary 'v' prefix before variables in ODE m-files. Rik <rik@octave.org> parents: 20631 diff changeset	80 ##
20719 d903cccb8de8 doc: Periodic spellcheck of documentation. Rik <rik@octave.org> parents: 20634 diff changeset	81 ## Example: Solve the @nospell{Van der Pol} equation
20548 25623ef2ff4f doc: Rewrite docstrings for ode* family of functions. Rik <rik@octave.org> parents: 20543 diff changeset	82 ##
20533 fcb792acab9b Moving ode45, odeset, odeget, and levenshtein from odepkg to core. jcorno <jacopo.corno@gmail.com> parents: diff changeset	83 ## @example
20548 25623ef2ff4f doc: Rewrite docstrings for ode* family of functions. Rik <rik@octave.org> parents: 20543 diff changeset	84 ## @group
20634 80e630b37ba1 maint: Remove unnecessary 'v' prefix before variables in ODE m-files. Rik <rik@octave.org> parents: 20631 diff changeset	85 ## fvdp = @@(@var{t},@var{y}) [@var{y}(2); (1 - @var{y}(1)^2) * @var{y}(2) - @var{y}(1)];
21443 acd6e203031d Alter BIST tests stop emitting warnings during runtests invocation. Rik <rik@octave.org> parents: 20908 diff changeset	86 ## [@var{t},@var{y}] = ode45 (fvdp, [0, 20], [2, 0]);
20548 25623ef2ff4f doc: Rewrite docstrings for ode* family of functions. Rik <rik@octave.org> parents: 20543 diff changeset	87 ## @end group
20533 fcb792acab9b Moving ode45, odeset, odeget, and levenshtein from odepkg to core. jcorno <jacopo.corno@gmail.com> parents: diff changeset	88 ## @end example
25026 f886561f9696 doc: improve differential eqtn docs and mention ode15i/ode15s (bug #51965). Colin Macdonald <cbm@m.fsf.org> parents: 24534 diff changeset	89 ## @seealso{odeset, odeget, ode23, ode15s}
20533 fcb792acab9b Moving ode45, odeset, odeget, and levenshtein from odepkg to core. jcorno <jacopo.corno@gmail.com> parents: diff changeset	90 ## @end deftypefn
fcb792acab9b Moving ode45, odeset, odeget, and levenshtein from odepkg to core. jcorno <jacopo.corno@gmail.com> parents: diff changeset	91
20634 80e630b37ba1 maint: Remove unnecessary 'v' prefix before variables in ODE m-files. Rik <rik@octave.org> parents: 20631 diff changeset	92 function varargout = ode45 (fun, trange, init, varargin)
20533 fcb792acab9b Moving ode45, odeset, odeget, and levenshtein from odepkg to core. jcorno <jacopo.corno@gmail.com> parents: diff changeset	93
fcb792acab9b Moving ode45, odeset, odeget, and levenshtein from odepkg to core. jcorno <jacopo.corno@gmail.com> parents: diff changeset	94 if (nargin < 3)
20553 45151de7423f maint: Clean up implementations of ode45.m, odeget.m, odeset.m. Rik <rik@octave.org> parents: 20549 diff changeset	95 print_usage ();
20533 fcb792acab9b Moving ode45, odeset, odeget, and levenshtein from odepkg to core. jcorno <jacopo.corno@gmail.com> parents: diff changeset	96 endif
20634 80e630b37ba1 maint: Remove unnecessary 'v' prefix before variables in ODE m-files. Rik <rik@octave.org> parents: 20631 diff changeset	97
28045 13dba3c069f8 Update input validation for odeXXX.m functions. Rik <rik@octave.org> parents: 28041 diff changeset	98 solver = "ode45";
22593 dba5074bdc79 simplify options management in ode solvers Carlo de Falco <carlo.defalco@polimi.it> parents: 22445 diff changeset	99 order = 5; # runge_kutta_45_dorpri uses local extrapolation
20631 00caf63edcdf maint: Remove obsolete ODE options from odeset/odeget/ode45. Rik <rik@octave.org> parents: 20620 diff changeset	100
20533 fcb792acab9b Moving ode45, odeset, odeget, and levenshtein from odepkg to core. jcorno <jacopo.corno@gmail.com> parents: diff changeset	101 if (nargin >= 4)
20549 e368ce72a844 maint: Use Octave coding conventions for ode* functions. Rik <rik@octave.org> parents: 20548 diff changeset	102 if (! isstruct (varargin{1}))
20634 80e630b37ba1 maint: Remove unnecessary 'v' prefix before variables in ODE m-files. Rik <rik@octave.org> parents: 20631 diff changeset	103 ## varargin{1:len} are parameters for fun
80e630b37ba1 maint: Remove unnecessary 'v' prefix before variables in ODE m-files. Rik <rik@octave.org> parents: 20631 diff changeset	104 odeopts = odeset ();
22593 dba5074bdc79 simplify options management in ode solvers Carlo de Falco <carlo.defalco@polimi.it> parents: 22445 diff changeset	105 funarguments = varargin;
22626 869c02fde46c Further clean-up of ode functions. Rik <rik@octave.org> parents: 22625 diff changeset	106 elseif (numel (varargin) > 1)
22593 dba5074bdc79 simplify options management in ode solvers Carlo de Falco <carlo.defalco@polimi.it> parents: 22445 diff changeset	107 ## varargin{1} is an ODE options structure opt
dba5074bdc79 simplify options management in ode solvers Carlo de Falco <carlo.defalco@polimi.it> parents: 22445 diff changeset	108 odeopts = varargin{1};
22626 869c02fde46c Further clean-up of ode functions. Rik <rik@octave.org> parents: 22625 diff changeset	109 funarguments = {varargin{2:numel (varargin)}};
28045 13dba3c069f8 Update input validation for odeXXX.m functions. Rik <rik@octave.org> parents: 28041 diff changeset	110 else
13dba3c069f8 Update input validation for odeXXX.m functions. Rik <rik@octave.org> parents: 28041 diff changeset	111 ## varargin{1} is an ODE options structure opt
22593 dba5074bdc79 simplify options management in ode solvers Carlo de Falco <carlo.defalco@polimi.it> parents: 22445 diff changeset	112 odeopts = varargin{1};
dba5074bdc79 simplify options management in ode solvers Carlo de Falco <carlo.defalco@polimi.it> parents: 22445 diff changeset	113 funarguments = {};
20533 fcb792acab9b Moving ode45, odeset, odeget, and levenshtein from odepkg to core. jcorno <jacopo.corno@gmail.com> parents: diff changeset	114 endif
20549 e368ce72a844 maint: Use Octave coding conventions for ode* functions. Rik <rik@octave.org> parents: 20548 diff changeset	115 else # nargin == 3
20634 80e630b37ba1 maint: Remove unnecessary 'v' prefix before variables in ODE m-files. Rik <rik@octave.org> parents: 20631 diff changeset	116 odeopts = odeset ();
22593 dba5074bdc79 simplify options management in ode solvers Carlo de Falco <carlo.defalco@polimi.it> parents: 22445 diff changeset	117 funarguments = {};
20533 fcb792acab9b Moving ode45, odeset, odeget, and levenshtein from odepkg to core. jcorno <jacopo.corno@gmail.com> parents: diff changeset	118 endif
fcb792acab9b Moving ode45, odeset, odeget, and levenshtein from odepkg to core. jcorno <jacopo.corno@gmail.com> parents: diff changeset	119
21443 acd6e203031d Alter BIST tests stop emitting warnings during runtests invocation. Rik <rik@octave.org> parents: 20908 diff changeset	120 if (! isnumeric (trange) \|\| ! isvector (trange))
20631 00caf63edcdf maint: Remove obsolete ODE options from odeset/odeget/ode45. Rik <rik@octave.org> parents: 20620 diff changeset	121 error ("Octave:invalid-input-arg",
20634 80e630b37ba1 maint: Remove unnecessary 'v' prefix before variables in ODE m-files. Rik <rik@octave.org> parents: 20631 diff changeset	122 "ode45: TRANGE must be a numeric vector");
20533 fcb792acab9b Moving ode45, odeset, odeget, and levenshtein from odepkg to core. jcorno <jacopo.corno@gmail.com> parents: diff changeset	123 endif
fcb792acab9b Moving ode45, odeset, odeget, and levenshtein from odepkg to core. jcorno <jacopo.corno@gmail.com> parents: diff changeset	124
22626 869c02fde46c Further clean-up of ode functions. Rik <rik@octave.org> parents: 22625 diff changeset	125 if (numel (trange) < 2)
20631 00caf63edcdf maint: Remove obsolete ODE options from odeset/odeget/ode45. Rik <rik@octave.org> parents: 20620 diff changeset	126 error ("Octave:invalid-input-arg",
20634 80e630b37ba1 maint: Remove unnecessary 'v' prefix before variables in ODE m-files. Rik <rik@octave.org> parents: 20631 diff changeset	127 "ode45: TRANGE must contain at least 2 elements");
80e630b37ba1 maint: Remove unnecessary 'v' prefix before variables in ODE m-files. Rik <rik@octave.org> parents: 20631 diff changeset	128 elseif (trange(1) == trange(2))
20631 00caf63edcdf maint: Remove obsolete ODE options from odeset/odeget/ode45. Rik <rik@octave.org> parents: 20620 diff changeset	129 error ("Octave:invalid-input-arg",
20634 80e630b37ba1 maint: Remove unnecessary 'v' prefix before variables in ODE m-files. Rik <rik@octave.org> parents: 20631 diff changeset	130 "ode45: invalid time span, TRANGE(1) == TRANGE(2)");
20533 fcb792acab9b Moving ode45, odeset, odeget, and levenshtein from odepkg to core. jcorno <jacopo.corno@gmail.com> parents: diff changeset	131 else
22593 dba5074bdc79 simplify options management in ode solvers Carlo de Falco <carlo.defalco@polimi.it> parents: 22445 diff changeset	132 direction = sign (trange(2) - trange(1));
20533 fcb792acab9b Moving ode45, odeset, odeget, and levenshtein from odepkg to core. jcorno <jacopo.corno@gmail.com> parents: diff changeset	133 endif
20634 80e630b37ba1 maint: Remove unnecessary 'v' prefix before variables in ODE m-files. Rik <rik@octave.org> parents: 20631 diff changeset	134 trange = trange(:);
20533 fcb792acab9b Moving ode45, odeset, odeget, and levenshtein from odepkg to core. jcorno <jacopo.corno@gmail.com> parents: diff changeset	135
21443 acd6e203031d Alter BIST tests stop emitting warnings during runtests invocation. Rik <rik@octave.org> parents: 20908 diff changeset	136 if (! isnumeric (init) \|\| ! isvector (init))
20631 00caf63edcdf maint: Remove obsolete ODE options from odeset/odeget/ode45. Rik <rik@octave.org> parents: 20620 diff changeset	137 error ("Octave:invalid-input-arg",
20634 80e630b37ba1 maint: Remove unnecessary 'v' prefix before variables in ODE m-files. Rik <rik@octave.org> parents: 20631 diff changeset	138 "ode45: INIT must be a numeric vector");
20533 fcb792acab9b Moving ode45, odeset, odeget, and levenshtein from odepkg to core. jcorno <jacopo.corno@gmail.com> parents: diff changeset	139 endif
20634 80e630b37ba1 maint: Remove unnecessary 'v' prefix before variables in ODE m-files. Rik <rik@octave.org> parents: 20631 diff changeset	140 init = init(:);
20533 fcb792acab9b Moving ode45, odeset, odeget, and levenshtein from odepkg to core. jcorno <jacopo.corno@gmail.com> parents: diff changeset	141
20634 80e630b37ba1 maint: Remove unnecessary 'v' prefix before variables in ODE m-files. Rik <rik@octave.org> parents: 20631 diff changeset	142 if (ischar (fun))
28041 5e44268dca6f Replace input validation relying on str2func with alternatives (bug #57351). Rik <rik@octave.org> parents: 27923 diff changeset	143 if (! exist (fun))
5e44268dca6f Replace input validation relying on str2func with alternatives (bug #57351). Rik <rik@octave.org> parents: 27923 diff changeset	144 error ("Octave:invalid-input-arg",
28045 13dba3c069f8 Update input validation for odeXXX.m functions. Rik <rik@octave.org> parents: 28041 diff changeset	145 ['ode45: function "' fun '" not found']);
28041 5e44268dca6f Replace input validation relying on str2func with alternatives (bug #57351). Rik <rik@octave.org> parents: 27923 diff changeset	146 endif
5e44268dca6f Replace input validation relying on str2func with alternatives (bug #57351). Rik <rik@octave.org> parents: 27923 diff changeset	147 fun = str2func (fun);
20613 27c091f4b66d allow first argument in ode45 to be a string Carlo de Falco <carlo.defalco@polimi.it> parents: 20612 diff changeset	148 endif
25803 23483673ba43 Use is_function_handle instead of isa (x, "function_handle"). Rik <rik@octave.org> parents: 25054 diff changeset	149 if (! is_function_handle (fun))
20631 00caf63edcdf maint: Remove obsolete ODE options from odeset/odeget/ode45. Rik <rik@octave.org> parents: 20620 diff changeset	150 error ("Octave:invalid-input-arg",
28045 13dba3c069f8 Update input validation for odeXXX.m functions. Rik <rik@octave.org> parents: 28041 diff changeset	151 "ode45: FUN must be a valid function handle");
20533 fcb792acab9b Moving ode45, odeset, odeget, and levenshtein from odepkg to core. jcorno <jacopo.corno@gmail.com> parents: diff changeset	152 endif
fcb792acab9b Moving ode45, odeset, odeget, and levenshtein from odepkg to core. jcorno <jacopo.corno@gmail.com> parents: diff changeset	153
20634 80e630b37ba1 maint: Remove unnecessary 'v' prefix before variables in ODE m-files. Rik <rik@octave.org> parents: 20631 diff changeset	154 ## Start preprocessing, have a look which options are set in odeopts,
20549 e368ce72a844 maint: Use Octave coding conventions for ode* functions. Rik <rik@octave.org> parents: 20548 diff changeset	155 ## check if an invalid or unused option is set
22626 869c02fde46c Further clean-up of ode functions. Rik <rik@octave.org> parents: 22625 diff changeset	156 [defaults, classes, attributes] = odedefaults (numel (init),
869c02fde46c Further clean-up of ode functions. Rik <rik@octave.org> parents: 22625 diff changeset	157 trange(1), trange(end));
22593 dba5074bdc79 simplify options management in ode solvers Carlo de Falco <carlo.defalco@polimi.it> parents: 22445 diff changeset	158
22660 7eb3f8ec1aed Remove FIXME comment about missing tests for ignored ODE options. Carlo de Falco <carlo.defalco@polimi.it> parents: 22659 diff changeset	159 ## FIXME: Refine is not correctly implemented yet
7eb3f8ec1aed Remove FIXME comment about missing tests for ignored ODE options. Carlo de Falco <carlo.defalco@polimi.it> parents: 22659 diff changeset	160 defaults = odeset (defaults, "Refine", 4);
7eb3f8ec1aed Remove FIXME comment about missing tests for ignored ODE options. Carlo de Falco <carlo.defalco@polimi.it> parents: 22659 diff changeset	161
7eb3f8ec1aed Remove FIXME comment about missing tests for ignored ODE options. Carlo de Falco <carlo.defalco@polimi.it> parents: 22659 diff changeset	162 persistent ode45_ignore_options = ...
7eb3f8ec1aed Remove FIXME comment about missing tests for ignored ODE options. Carlo de Falco <carlo.defalco@polimi.it> parents: 22659 diff changeset	163 {"BDF", "InitialSlope", "Jacobian", "JPattern",
7eb3f8ec1aed Remove FIXME comment about missing tests for ignored ODE options. Carlo de Falco <carlo.defalco@polimi.it> parents: 22659 diff changeset	164 "MassSingular", "MaxOrder", "MvPattern", "Vectorized"};
7eb3f8ec1aed Remove FIXME comment about missing tests for ignored ODE options. Carlo de Falco <carlo.defalco@polimi.it> parents: 22659 diff changeset	165
7eb3f8ec1aed Remove FIXME comment about missing tests for ignored ODE options. Carlo de Falco <carlo.defalco@polimi.it> parents: 22659 diff changeset	166 defaults = rmfield (defaults, ode45_ignore_options);
7eb3f8ec1aed Remove FIXME comment about missing tests for ignored ODE options. Carlo de Falco <carlo.defalco@polimi.it> parents: 22659 diff changeset	167 classes = rmfield (classes, ode45_ignore_options);
7eb3f8ec1aed Remove FIXME comment about missing tests for ignored ODE options. Carlo de Falco <carlo.defalco@polimi.it> parents: 22659 diff changeset	168 attributes = rmfield (attributes, ode45_ignore_options);
20533 fcb792acab9b Moving ode45, odeset, odeget, and levenshtein from odepkg to core. jcorno <jacopo.corno@gmail.com> parents: diff changeset	169
22626 869c02fde46c Further clean-up of ode functions. Rik <rik@octave.org> parents: 22625 diff changeset	170 odeopts = odemergeopts ("ode45", odeopts, defaults, classes, attributes);
20533 fcb792acab9b Moving ode45, odeset, odeget, and levenshtein from odepkg to core. jcorno <jacopo.corno@gmail.com> parents: diff changeset	171
22593 dba5074bdc79 simplify options management in ode solvers Carlo de Falco <carlo.defalco@polimi.it> parents: 22445 diff changeset	172 odeopts.funarguments = funarguments;
dba5074bdc79 simplify options management in ode solvers Carlo de Falco <carlo.defalco@polimi.it> parents: 22445 diff changeset	173 odeopts.direction = direction;
20634 80e630b37ba1 maint: Remove unnecessary 'v' prefix before variables in ODE m-files. Rik <rik@octave.org> parents: 20631 diff changeset	174
80e630b37ba1 maint: Remove unnecessary 'v' prefix before variables in ODE m-files. Rik <rik@octave.org> parents: 20631 diff changeset	175 if (! isempty (odeopts.NonNegative))
80e630b37ba1 maint: Remove unnecessary 'v' prefix before variables in ODE m-files. Rik <rik@octave.org> parents: 20631 diff changeset	176 if (isempty (odeopts.Mass))
80e630b37ba1 maint: Remove unnecessary 'v' prefix before variables in ODE m-files. Rik <rik@octave.org> parents: 20631 diff changeset	177 odeopts.havenonnegative = true;
80e630b37ba1 maint: Remove unnecessary 'v' prefix before variables in ODE m-files. Rik <rik@octave.org> parents: 20631 diff changeset	178 else
80e630b37ba1 maint: Remove unnecessary 'v' prefix before variables in ODE m-files. Rik <rik@octave.org> parents: 20631 diff changeset	179 odeopts.havenonnegative = false;
80e630b37ba1 maint: Remove unnecessary 'v' prefix before variables in ODE m-files. Rik <rik@octave.org> parents: 20631 diff changeset	180 warning ("Octave:invalid-input-arg",
22626 869c02fde46c Further clean-up of ode functions. Rik <rik@octave.org> parents: 22625 diff changeset	181 ['ode45: option "NonNegative" is ignored', ...
20634 80e630b37ba1 maint: Remove unnecessary 'v' prefix before variables in ODE m-files. Rik <rik@octave.org> parents: 20631 diff changeset	182 " when mass matrix is set\n"]);
80e630b37ba1 maint: Remove unnecessary 'v' prefix before variables in ODE m-files. Rik <rik@octave.org> parents: 20631 diff changeset	183 endif
80e630b37ba1 maint: Remove unnecessary 'v' prefix before variables in ODE m-files. Rik <rik@octave.org> parents: 20631 diff changeset	184 else
80e630b37ba1 maint: Remove unnecessary 'v' prefix before variables in ODE m-files. Rik <rik@octave.org> parents: 20631 diff changeset	185 odeopts.havenonnegative = false;
20533 fcb792acab9b Moving ode45, odeset, odeget, and levenshtein from odepkg to core. jcorno <jacopo.corno@gmail.com> parents: diff changeset	186 endif
fcb792acab9b Moving ode45, odeset, odeget, and levenshtein from odepkg to core. jcorno <jacopo.corno@gmail.com> parents: diff changeset	187
20634 80e630b37ba1 maint: Remove unnecessary 'v' prefix before variables in ODE m-files. Rik <rik@octave.org> parents: 20631 diff changeset	188 if (isempty (odeopts.OutputFcn) && nargout == 0)
80e630b37ba1 maint: Remove unnecessary 'v' prefix before variables in ODE m-files. Rik <rik@octave.org> parents: 20631 diff changeset	189 odeopts.OutputFcn = @odeplot;
80e630b37ba1 maint: Remove unnecessary 'v' prefix before variables in ODE m-files. Rik <rik@octave.org> parents: 20631 diff changeset	190 odeopts.haveoutputfunction = true;
80e630b37ba1 maint: Remove unnecessary 'v' prefix before variables in ODE m-files. Rik <rik@octave.org> parents: 20631 diff changeset	191 else
80e630b37ba1 maint: Remove unnecessary 'v' prefix before variables in ODE m-files. Rik <rik@octave.org> parents: 20631 diff changeset	192 odeopts.haveoutputfunction = ! isempty (odeopts.OutputFcn);
20533 fcb792acab9b Moving ode45, odeset, odeget, and levenshtein from odepkg to core. jcorno <jacopo.corno@gmail.com> parents: diff changeset	193 endif
fcb792acab9b Moving ode45, odeset, odeget, and levenshtein from odepkg to core. jcorno <jacopo.corno@gmail.com> parents: diff changeset	194
22593 dba5074bdc79 simplify options management in ode solvers Carlo de Falco <carlo.defalco@polimi.it> parents: 22445 diff changeset	195 if (isempty (odeopts.InitialStep))
dba5074bdc79 simplify options management in ode solvers Carlo de Falco <carlo.defalco@polimi.it> parents: 22445 diff changeset	196 odeopts.InitialStep = odeopts.direction * ...
22626 869c02fde46c Further clean-up of ode functions. Rik <rik@octave.org> parents: 22625 diff changeset	197 starting_stepsize (order, fun, trange(1), init,
869c02fde46c Further clean-up of ode functions. Rik <rik@octave.org> parents: 22625 diff changeset	198 odeopts.AbsTol, odeopts.RelTol,
22935 c9344df03da5 Allow case-insensitive option argument 'on' to ode solvers (bug #49918). Rik <rik@octave.org> parents: 22934 diff changeset	199 strcmpi (odeopts.NormControl, "on"),
22626 869c02fde46c Further clean-up of ode functions. Rik <rik@octave.org> parents: 22625 diff changeset	200 odeopts.funarguments);
22625 081a201b77c7 Clean up ode options implementation to follow Octave coding standards. Carlo de Falco <carlo.defalco@polimi.it> parents: 22605 diff changeset	201 endif
20533 fcb792acab9b Moving ode45, odeset, odeget, and levenshtein from odepkg to core. jcorno <jacopo.corno@gmail.com> parents: diff changeset	202
28045 13dba3c069f8 Update input validation for odeXXX.m functions. Rik <rik@octave.org> parents: 28041 diff changeset	203 if (! isempty (odeopts.Mass))
13dba3c069f8 Update input validation for odeXXX.m functions. Rik <rik@octave.org> parents: 28041 diff changeset	204 if (isnumeric (odeopts.Mass))
13dba3c069f8 Update input validation for odeXXX.m functions. Rik <rik@octave.org> parents: 28041 diff changeset	205 havemasshandle = false;
13dba3c069f8 Update input validation for odeXXX.m functions. Rik <rik@octave.org> parents: 28041 diff changeset	206 mass = odeopts.Mass; # constant mass
13dba3c069f8 Update input validation for odeXXX.m functions. Rik <rik@octave.org> parents: 28041 diff changeset	207 elseif (is_function_handle (odeopts.Mass))
13dba3c069f8 Update input validation for odeXXX.m functions. Rik <rik@octave.org> parents: 28041 diff changeset	208 havemasshandle = true; # mass defined by a function handle
13dba3c069f8 Update input validation for odeXXX.m functions. Rik <rik@octave.org> parents: 28041 diff changeset	209 else
13dba3c069f8 Update input validation for odeXXX.m functions. Rik <rik@octave.org> parents: 28041 diff changeset	210 error ("Octave:invalid-input-arg",
13dba3c069f8 Update input validation for odeXXX.m functions. Rik <rik@octave.org> parents: 28041 diff changeset	211 'ode45: "Mass" field must be a function handle or square matrix');
13dba3c069f8 Update input validation for odeXXX.m functions. Rik <rik@octave.org> parents: 28041 diff changeset	212 endif
13dba3c069f8 Update input validation for odeXXX.m functions. Rik <rik@octave.org> parents: 28041 diff changeset	213 else # no mass matrix - create a diag-matrix of ones for mass
20634 80e630b37ba1 maint: Remove unnecessary 'v' prefix before variables in ODE m-files. Rik <rik@octave.org> parents: 20631 diff changeset	214 havemasshandle = false; # mass = diag (ones (length (init), 1), 0);
20533 fcb792acab9b Moving ode45, odeset, odeget, and levenshtein from odepkg to core. jcorno <jacopo.corno@gmail.com> parents: diff changeset	215 endif
fcb792acab9b Moving ode45, odeset, odeget, and levenshtein from odepkg to core. jcorno <jacopo.corno@gmail.com> parents: diff changeset	216
21443 acd6e203031d Alter BIST tests stop emitting warnings during runtests invocation. Rik <rik@octave.org> parents: 20908 diff changeset	217 ## Starting the initialization of the core solver ode45
20634 80e630b37ba1 maint: Remove unnecessary 'v' prefix before variables in ODE m-files. Rik <rik@octave.org> parents: 20631 diff changeset	218
80e630b37ba1 maint: Remove unnecessary 'v' prefix before variables in ODE m-files. Rik <rik@octave.org> parents: 20631 diff changeset	219 if (havemasshandle) # Handle only the dynamic mass matrix,
22626 869c02fde46c Further clean-up of ode functions. Rik <rik@octave.org> parents: 22625 diff changeset	220 if (! strcmp (odeopts.MStateDependence, "none"))
869c02fde46c Further clean-up of ode functions. Rik <rik@octave.org> parents: 22625 diff changeset	221 ### constant mass matrices have already
20634 80e630b37ba1 maint: Remove unnecessary 'v' prefix before variables in ODE m-files. Rik <rik@octave.org> parents: 20631 diff changeset	222 mass = @(t,x) odeopts.Mass (t, x, odeopts.funarguments{:});
80e630b37ba1 maint: Remove unnecessary 'v' prefix before variables in ODE m-files. Rik <rik@octave.org> parents: 20631 diff changeset	223 fun = @(t,x) mass (t, x, odeopts.funarguments{:}) ...
22626 869c02fde46c Further clean-up of ode functions. Rik <rik@octave.org> parents: 22625 diff changeset	224 \ fun (t, x, odeopts.funarguments{:});
869c02fde46c Further clean-up of ode functions. Rik <rik@octave.org> parents: 22625 diff changeset	225 else
20634 80e630b37ba1 maint: Remove unnecessary 'v' prefix before variables in ODE m-files. Rik <rik@octave.org> parents: 20631 diff changeset	226 mass = @(t) odeopts.Mass (t, odeopts.funarguments{:});
80e630b37ba1 maint: Remove unnecessary 'v' prefix before variables in ODE m-files. Rik <rik@octave.org> parents: 20631 diff changeset	227 fun = @(t,x) mass (t, odeopts.funarguments{:}) ...
22626 869c02fde46c Further clean-up of ode functions. Rik <rik@octave.org> parents: 22625 diff changeset	228 \ fun (t, x, odeopts.funarguments{:});
20533 fcb792acab9b Moving ode45, odeset, odeget, and levenshtein from odepkg to core. jcorno <jacopo.corno@gmail.com> parents: diff changeset	229 endif
fcb792acab9b Moving ode45, odeset, odeget, and levenshtein from odepkg to core. jcorno <jacopo.corno@gmail.com> parents: diff changeset	230 endif
fcb792acab9b Moving ode45, odeset, odeget, and levenshtein from odepkg to core. jcorno <jacopo.corno@gmail.com> parents: diff changeset	231
22659 04fc7e9c5f96 Disable Refine option of ode solvers for certain input/output combinations (partial fix bug #49408). Rik <rik@octave.org> parents: 22658 diff changeset	232 if (nargout == 1)
04fc7e9c5f96 Disable Refine option of ode solvers for certain input/output combinations (partial fix bug #49408). Rik <rik@octave.org> parents: 22658 diff changeset	233 ## Single output requires auto-selected intermediate times,
04fc7e9c5f96 Disable Refine option of ode solvers for certain input/output combinations (partial fix bug #49408). Rik <rik@octave.org> parents: 22658 diff changeset	234 ## which is obtained by NOT specifying specific solution times.
04fc7e9c5f96 Disable Refine option of ode solvers for certain input/output combinations (partial fix bug #49408). Rik <rik@octave.org> parents: 22658 diff changeset	235 trange = [trange(1); trange(end)];
04fc7e9c5f96 Disable Refine option of ode solvers for certain input/output combinations (partial fix bug #49408). Rik <rik@octave.org> parents: 22658 diff changeset	236 odeopts.Refine = []; # disable Refine when single output requested
04fc7e9c5f96 Disable Refine option of ode solvers for certain input/output combinations (partial fix bug #49408). Rik <rik@octave.org> parents: 22658 diff changeset	237 elseif (numel (trange) > 2)
04fc7e9c5f96 Disable Refine option of ode solvers for certain input/output combinations (partial fix bug #49408). Rik <rik@octave.org> parents: 22658 diff changeset	238 odeopts.Refine = []; # disable Refine when specific times requested
22658 56d7d423aff9 Return intermediate integration times when single output argument used in ode solver (bug #49417). Rik <rik@octave.org> parents: 22655 diff changeset	239 endif
56d7d423aff9 Return intermediate integration times when single output argument used in ode solver (bug #49417). Rik <rik@octave.org> parents: 22655 diff changeset	240
22593 dba5074bdc79 simplify options management in ode solvers Carlo de Falco <carlo.defalco@polimi.it> parents: 22445 diff changeset	241 solution = integrate_adaptive (@runge_kutta_45_dorpri,
dba5074bdc79 simplify options management in ode solvers Carlo de Falco <carlo.defalco@polimi.it> parents: 22445 diff changeset	242 order, fun, trange, init, odeopts);
22625 081a201b77c7 Clean up ode options implementation to follow Octave coding standards. Carlo de Falco <carlo.defalco@polimi.it> parents: 22605 diff changeset	243
20533 fcb792acab9b Moving ode45, odeset, odeget, and levenshtein from odepkg to core. jcorno <jacopo.corno@gmail.com> parents: diff changeset	244 ## Postprocessing, do whatever when terminating integration algorithm
20634 80e630b37ba1 maint: Remove unnecessary 'v' prefix before variables in ODE m-files. Rik <rik@octave.org> parents: 20631 diff changeset	245 if (odeopts.haveoutputfunction) # Cleanup plotter
22639 7efa2d0e22c9 More Matlab-compatible implementation of OutputFcn. Clean up odeplot.m. Rik <rik@octave.org> parents: 22626 diff changeset	246 feval (odeopts.OutputFcn, [], [], "done", odeopts.funarguments{:});
20533 fcb792acab9b Moving ode45, odeset, odeget, and levenshtein from odepkg to core. jcorno <jacopo.corno@gmail.com> parents: diff changeset	247 endif
22593 dba5074bdc79 simplify options management in ode solvers Carlo de Falco <carlo.defalco@polimi.it> parents: 22445 diff changeset	248 if (! isempty (odeopts.Events)) # Cleanup event function handling
22626 869c02fde46c Further clean-up of ode functions. Rik <rik@octave.org> parents: 22625 diff changeset	249 ode_event_handler (odeopts.Events, solution.t(end),
22655 6b134d294d61 ode solvers: use ordinary transpose instead of Hermitian conjugate (bug #49410). Carlo de Falco <carlo.defalco@polimi.it> parents: 22654 diff changeset	250 solution.x(end,:).', "done", odeopts.funarguments{:});
20533 fcb792acab9b Moving ode45, odeset, odeget, and levenshtein from odepkg to core. jcorno <jacopo.corno@gmail.com> parents: diff changeset	251 endif
fcb792acab9b Moving ode45, odeset, odeget, and levenshtein from odepkg to core. jcorno <jacopo.corno@gmail.com> parents: diff changeset	252
fcb792acab9b Moving ode45, odeset, odeget, and levenshtein from odepkg to core. jcorno <jacopo.corno@gmail.com> parents: diff changeset	253 ## Print additional information if option Stats is set
22935 c9344df03da5 Allow case-insensitive option argument 'on' to ode solvers (bug #49918). Rik <rik@octave.org> parents: 22934 diff changeset	254 if (strcmpi (odeopts.Stats, "on"))
21955 0f3e875d9078 Fix statistics on solution solving for ode45, ode23 (bug #48243). jcorno <jacopo.corno@gmail.com> parents: 21634 diff changeset	255 nsteps = solution.cntloop; # cntloop from 2..end
0f3e875d9078 Fix statistics on solution solving for ode45, ode23 (bug #48243). jcorno <jacopo.corno@gmail.com> parents: 21634 diff changeset	256 nfailed = solution.cntcycles - nsteps; # cntcycl from 1..end
0f3e875d9078 Fix statistics on solution solving for ode45, ode23 (bug #48243). jcorno <jacopo.corno@gmail.com> parents: 21634 diff changeset	257 nfevals = 6 * solution.cntcycles + 1; # number of ode evaluations
20634 80e630b37ba1 maint: Remove unnecessary 'v' prefix before variables in ODE m-files. Rik <rik@octave.org> parents: 20631 diff changeset	258 ndecomps = 0; # number of LU decompositions
80e630b37ba1 maint: Remove unnecessary 'v' prefix before variables in ODE m-files. Rik <rik@octave.org> parents: 20631 diff changeset	259 npds = 0; # number of partial derivatives
80e630b37ba1 maint: Remove unnecessary 'v' prefix before variables in ODE m-files. Rik <rik@octave.org> parents: 20631 diff changeset	260 nlinsols = 0; # no. of linear systems solutions
22934 bb452f84a299 Fix printing of solution stats in ode23/ode45 solvers (bug #49918). Carlo de Falco <carlo.defalco@polimi.it> parents: 22755 diff changeset	261
bb452f84a299 Fix printing of solution stats in ode23/ode45 solvers (bug #49918). Carlo de Falco <carlo.defalco@polimi.it> parents: 22755 diff changeset	262 printf ("Number of successful steps: %d\n", nsteps);
bb452f84a299 Fix printing of solution stats in ode23/ode45 solvers (bug #49918). Carlo de Falco <carlo.defalco@polimi.it> parents: 22755 diff changeset	263 printf ("Number of failed attempts: %d\n", nfailed);
bb452f84a299 Fix printing of solution stats in ode23/ode45 solvers (bug #49918). Carlo de Falco <carlo.defalco@polimi.it> parents: 22755 diff changeset	264 printf ("Number of function calls: %d\n", nfevals);
20533 fcb792acab9b Moving ode45, odeset, odeget, and levenshtein from odepkg to core. jcorno <jacopo.corno@gmail.com> parents: diff changeset	265 endif
fcb792acab9b Moving ode45, odeset, odeget, and levenshtein from odepkg to core. jcorno <jacopo.corno@gmail.com> parents: diff changeset	266
20553 45151de7423f maint: Clean up implementations of ode45.m, odeget.m, odeset.m. Rik <rik@octave.org> parents: 20549 diff changeset	267 if (nargout == 2)
45151de7423f maint: Clean up implementations of ode45.m, odeget.m, odeset.m. Rik <rik@octave.org> parents: 20549 diff changeset	268 varargout{1} = solution.t; # Time stamps are first output argument
45151de7423f maint: Clean up implementations of ode45.m, odeget.m, odeset.m. Rik <rik@octave.org> parents: 20549 diff changeset	269 varargout{2} = solution.x; # Results are second output argument
45151de7423f maint: Clean up implementations of ode45.m, odeget.m, odeset.m. Rik <rik@octave.org> parents: 20549 diff changeset	270 elseif (nargout == 1)
22934 bb452f84a299 Fix printing of solution stats in ode23/ode45 solvers (bug #49918). Carlo de Falco <carlo.defalco@polimi.it> parents: 22755 diff changeset	271 varargout{1}.x = solution.t.'; # Time stamps saved in field x (row vector)
bb452f84a299 Fix printing of solution stats in ode23/ode45 solvers (bug #49918). Carlo de Falco <carlo.defalco@polimi.it> parents: 22755 diff changeset	272 varargout{1}.y = solution.x.'; # Results are saved in field y (row vector)
20634 80e630b37ba1 maint: Remove unnecessary 'v' prefix before variables in ODE m-files. Rik <rik@octave.org> parents: 20631 diff changeset	273 varargout{1}.solver = solver; # Solver name is saved in field solver
22593 dba5074bdc79 simplify options management in ode solvers Carlo de Falco <carlo.defalco@polimi.it> parents: 22445 diff changeset	274 if (! isempty (odeopts.Events))
20634 80e630b37ba1 maint: Remove unnecessary 'v' prefix before variables in ODE m-files. Rik <rik@octave.org> parents: 20631 diff changeset	275 varargout{1}.xe = solution.event{3}; # Time info when an event occurred
80e630b37ba1 maint: Remove unnecessary 'v' prefix before variables in ODE m-files. Rik <rik@octave.org> parents: 20631 diff changeset	276 varargout{1}.ye = solution.event{4}; # Results when an event occurred
22955 3c72c72233e3 ode23.m, ode45.m: Switch output creation order of struct fields to match Matlab. Rik <rik@octave.org> parents: 22940 diff changeset	277 varargout{1}.ie = solution.event{2}; # Index info which event occurred
20533 fcb792acab9b Moving ode45, odeset, odeget, and levenshtein from odepkg to core. jcorno <jacopo.corno@gmail.com> parents: diff changeset	278 endif
22935 c9344df03da5 Allow case-insensitive option argument 'on' to ode solvers (bug #49918). Rik <rik@octave.org> parents: 22934 diff changeset	279 if (strcmpi (odeopts.Stats, "on"))
20553 45151de7423f maint: Clean up implementations of ode45.m, odeget.m, odeset.m. Rik <rik@octave.org> parents: 20549 diff changeset	280 varargout{1}.stats = struct ();
20634 80e630b37ba1 maint: Remove unnecessary 'v' prefix before variables in ODE m-files. Rik <rik@octave.org> parents: 20631 diff changeset	281 varargout{1}.stats.nsteps = nsteps;
80e630b37ba1 maint: Remove unnecessary 'v' prefix before variables in ODE m-files. Rik <rik@octave.org> parents: 20631 diff changeset	282 varargout{1}.stats.nfailed = nfailed;
80e630b37ba1 maint: Remove unnecessary 'v' prefix before variables in ODE m-files. Rik <rik@octave.org> parents: 20631 diff changeset	283 varargout{1}.stats.nfevals = nfevals;
80e630b37ba1 maint: Remove unnecessary 'v' prefix before variables in ODE m-files. Rik <rik@octave.org> parents: 20631 diff changeset	284 varargout{1}.stats.npds = npds;
80e630b37ba1 maint: Remove unnecessary 'v' prefix before variables in ODE m-files. Rik <rik@octave.org> parents: 20631 diff changeset	285 varargout{1}.stats.ndecomps = ndecomps;
80e630b37ba1 maint: Remove unnecessary 'v' prefix before variables in ODE m-files. Rik <rik@octave.org> parents: 20631 diff changeset	286 varargout{1}.stats.nlinsols = nlinsols;
20533 fcb792acab9b Moving ode45, odeset, odeget, and levenshtein from odepkg to core. jcorno <jacopo.corno@gmail.com> parents: diff changeset	287 endif
22937 2b2ac017c901 Fix undefined return argument for more than 2 outputs from ode solver (bug #49890). Rik <rik@octave.org> parents: 22660 diff changeset	288 elseif (nargout > 2)
20553 45151de7423f maint: Clean up implementations of ode45.m, odeget.m, odeset.m. Rik <rik@octave.org> parents: 20549 diff changeset	289 varargout = cell (1,5);
45151de7423f maint: Clean up implementations of ode45.m, odeget.m, odeset.m. Rik <rik@octave.org> parents: 20549 diff changeset	290 varargout{1} = solution.t;
45151de7423f maint: Clean up implementations of ode45.m, odeget.m, odeset.m. Rik <rik@octave.org> parents: 20549 diff changeset	291 varargout{2} = solution.x;
22593 dba5074bdc79 simplify options management in ode solvers Carlo de Falco <carlo.defalco@polimi.it> parents: 22445 diff changeset	292 if (! isempty (odeopts.Events))
20634 80e630b37ba1 maint: Remove unnecessary 'v' prefix before variables in ODE m-files. Rik <rik@octave.org> parents: 20631 diff changeset	293 varargout{3} = solution.event{3}; # Time info when an event occurred
80e630b37ba1 maint: Remove unnecessary 'v' prefix before variables in ODE m-files. Rik <rik@octave.org> parents: 20631 diff changeset	294 varargout{4} = solution.event{4}; # Results when an event occurred
80e630b37ba1 maint: Remove unnecessary 'v' prefix before variables in ODE m-files. Rik <rik@octave.org> parents: 20631 diff changeset	295 varargout{5} = solution.event{2}; # Index info which event occurred
20533 fcb792acab9b Moving ode45, odeset, odeget, and levenshtein from odepkg to core. jcorno <jacopo.corno@gmail.com> parents: diff changeset	296 endif
fcb792acab9b Moving ode45, odeset, odeget, and levenshtein from odepkg to core. jcorno <jacopo.corno@gmail.com> parents: diff changeset	297 endif
fcb792acab9b Moving ode45, odeset, odeget, and levenshtein from odepkg to core. jcorno <jacopo.corno@gmail.com> parents: diff changeset	298
fcb792acab9b Moving ode45, odeset, odeget, and levenshtein from odepkg to core. jcorno <jacopo.corno@gmail.com> parents: diff changeset	299 endfunction
fcb792acab9b Moving ode45, odeset, odeget, and levenshtein from odepkg to core. jcorno <jacopo.corno@gmail.com> parents: diff changeset	300
20549 e368ce72a844 maint: Use Octave coding conventions for ode* functions. Rik <rik@octave.org> parents: 20548 diff changeset	301
21599 f29d68e24c5a ode23.m, ode45.m: More reformatting of demos to Octave coding standards. Rik <rik@octave.org> parents: 21596 diff changeset	302 %!demo
f29d68e24c5a ode23.m, ode45.m: More reformatting of demos to Octave coding standards. Rik <rik@octave.org> parents: 21596 diff changeset	303 %! ## Demonstrate convergence order for ode45
f29d68e24c5a ode23.m, ode45.m: More reformatting of demos to Octave coding standards. Rik <rik@octave.org> parents: 21596 diff changeset	304 %! tol = 1e-5 ./ 10.^[0:8];
f29d68e24c5a ode23.m, ode45.m: More reformatting of demos to Octave coding standards. Rik <rik@octave.org> parents: 21596 diff changeset	305 %! for i = 1 : numel (tol)
f29d68e24c5a ode23.m, ode45.m: More reformatting of demos to Octave coding standards. Rik <rik@octave.org> parents: 21596 diff changeset	306 %! opt = odeset ("RelTol", tol(i), "AbsTol", realmin);
f29d68e24c5a ode23.m, ode45.m: More reformatting of demos to Octave coding standards. Rik <rik@octave.org> parents: 21596 diff changeset	307 %! [t, y] = ode45 (@(t, y) -y, [0, 1], 1, opt);
f29d68e24c5a ode23.m, ode45.m: More reformatting of demos to Octave coding standards. Rik <rik@octave.org> parents: 21596 diff changeset	308 %! h(i) = 1 / (numel (t) - 1);
f29d68e24c5a ode23.m, ode45.m: More reformatting of demos to Octave coding standards. Rik <rik@octave.org> parents: 21596 diff changeset	309 %! err(i) = norm (y .* exp (t) - 1, Inf);
f29d68e24c5a ode23.m, ode45.m: More reformatting of demos to Octave coding standards. Rik <rik@octave.org> parents: 21596 diff changeset	310 %! endfor
f29d68e24c5a ode23.m, ode45.m: More reformatting of demos to Octave coding standards. Rik <rik@octave.org> parents: 21596 diff changeset	311 %!
f29d68e24c5a ode23.m, ode45.m: More reformatting of demos to Octave coding standards. Rik <rik@octave.org> parents: 21596 diff changeset	312 %! ## Estimate order visually
f29d68e24c5a ode23.m, ode45.m: More reformatting of demos to Octave coding standards. Rik <rik@octave.org> parents: 21596 diff changeset	313 %! loglog (h, tol, "-ob",
22445 cbac595d23ac Fix %!demo blocks for ode23 and ode45 Mike Miller <mtmiller@octave.org> parents: 22323 diff changeset	314 %! h, err, "-b",
21599 f29d68e24c5a ode23.m, ode45.m: More reformatting of demos to Octave coding standards. Rik <rik@octave.org> parents: 21596 diff changeset	315 %! h, (h/h(end)) .^ 4 .* tol(end), "k--",
22626 869c02fde46c Further clean-up of ode functions. Rik <rik@octave.org> parents: 22625 diff changeset	316 %! h, (h/h(end)) .^ 5 .* tol(end), "k-");
21599 f29d68e24c5a ode23.m, ode45.m: More reformatting of demos to Octave coding standards. Rik <rik@octave.org> parents: 21596 diff changeset	317 %! axis tight
f29d68e24c5a ode23.m, ode45.m: More reformatting of demos to Octave coding standards. Rik <rik@octave.org> parents: 21596 diff changeset	318 %! xlabel ("h");
f29d68e24c5a ode23.m, ode45.m: More reformatting of demos to Octave coding standards. Rik <rik@octave.org> parents: 21596 diff changeset	319 %! ylabel ("err(h)");
f29d68e24c5a ode23.m, ode45.m: More reformatting of demos to Octave coding standards. Rik <rik@octave.org> parents: 21596 diff changeset	320 %! title ("Convergence plot for ode45");
f29d68e24c5a ode23.m, ode45.m: More reformatting of demos to Octave coding standards. Rik <rik@octave.org> parents: 21596 diff changeset	321 %! legend ("imposed tolerance", "ode45 (relative) error",
f29d68e24c5a ode23.m, ode45.m: More reformatting of demos to Octave coding standards. Rik <rik@octave.org> parents: 21596 diff changeset	322 %! "order 4", "order 5", "location", "northwest");
22647 cf8ec7abb4ae Clean up demo and test blocks for ode functions. Rik <rik@octave.org> parents: 22644 diff changeset	323 %!
cf8ec7abb4ae Clean up demo and test blocks for ode functions. Rik <rik@octave.org> parents: 22644 diff changeset	324 %! ## Estimate order numerically
cf8ec7abb4ae Clean up demo and test blocks for ode functions. Rik <rik@octave.org> parents: 22644 diff changeset	325 %! p = diff (log (err)) ./ diff (log (h))
21599 f29d68e24c5a ode23.m, ode45.m: More reformatting of demos to Octave coding standards. Rik <rik@octave.org> parents: 21596 diff changeset	326
22934 bb452f84a299 Fix printing of solution stats in ode23/ode45 solvers (bug #49918). Carlo de Falco <carlo.defalco@polimi.it> parents: 22755 diff changeset	327 ## We are using the Van der Pol equation for all tests.
bb452f84a299 Fix printing of solution stats in ode23/ode45 solvers (bug #49918). Carlo de Falco <carlo.defalco@polimi.it> parents: 22755 diff changeset	328 ## Further tests also define a reference solution (computed at high accuracy)
28563 5a07c798eb08 avoid function call input or output argument number mismatch John W. Eaton <jwe@octave.org> parents: 28045 diff changeset	329 %!function ydot = fpol (t, y, varargin) # The Van der Pol ODE
22626 869c02fde46c Further clean-up of ode functions. Rik <rik@octave.org> parents: 22625 diff changeset	330 %! ydot = [y(2); (1 - y(1)^2) * y(2) - y(1)];
20533 fcb792acab9b Moving ode45, odeset, odeget, and levenshtein from odepkg to core. jcorno <jacopo.corno@gmail.com> parents: diff changeset	331 %!endfunction
21443 acd6e203031d Alter BIST tests stop emitting warnings during runtests invocation. Rik <rik@octave.org> parents: 20908 diff changeset	332 %!function ref = fref () # The computed reference solution
22626 869c02fde46c Further clean-up of ode functions. Rik <rik@octave.org> parents: 22625 diff changeset	333 %! ref = [0.32331666704577, -1.83297456798624];
20533 fcb792acab9b Moving ode45, odeset, odeget, and levenshtein from odepkg to core. jcorno <jacopo.corno@gmail.com> parents: diff changeset	334 %!endfunction
20634 80e630b37ba1 maint: Remove unnecessary 'v' prefix before variables in ODE m-files. Rik <rik@octave.org> parents: 20631 diff changeset	335 %!function [val, trm, dir] = feve (t, y, varargin)
28563 5a07c798eb08 avoid function call input or output argument number mismatch John W. Eaton <jwe@octave.org> parents: 28045 diff changeset	336 %! val = fpol (t, y, varargin{:}); # We use the derivatives
22934 bb452f84a299 Fix printing of solution stats in ode23/ode45 solvers (bug #49918). Carlo de Falco <carlo.defalco@polimi.it> parents: 22755 diff changeset	337 %! trm = zeros (2,1); # that's why component 2
bb452f84a299 Fix printing of solution stats in ode23/ode45 solvers (bug #49918). Carlo de Falco <carlo.defalco@polimi.it> parents: 22755 diff changeset	338 %! dir = ones (2,1); # does not seem to be exact
20549 e368ce72a844 maint: Use Octave coding conventions for ode* functions. Rik <rik@octave.org> parents: 20548 diff changeset	339 %!endfunction
20634 80e630b37ba1 maint: Remove unnecessary 'v' prefix before variables in ODE m-files. Rik <rik@octave.org> parents: 20631 diff changeset	340 %!function [val, trm, dir] = fevn (t, y, varargin)
28563 5a07c798eb08 avoid function call input or output argument number mismatch John W. Eaton <jwe@octave.org> parents: 28045 diff changeset	341 %! val = fpol (t, y, varargin{:}); # We use the derivatives
22934 bb452f84a299 Fix printing of solution stats in ode23/ode45 solvers (bug #49918). Carlo de Falco <carlo.defalco@polimi.it> parents: 22755 diff changeset	342 %! trm = ones (2,1); # that's why component 2
bb452f84a299 Fix printing of solution stats in ode23/ode45 solvers (bug #49918). Carlo de Falco <carlo.defalco@polimi.it> parents: 22755 diff changeset	343 %! dir = ones (2,1); # does not seem to be exact
20549 e368ce72a844 maint: Use Octave coding conventions for ode* functions. Rik <rik@octave.org> parents: 20548 diff changeset	344 %!endfunction
20634 80e630b37ba1 maint: Remove unnecessary 'v' prefix before variables in ODE m-files. Rik <rik@octave.org> parents: 20631 diff changeset	345 %!function mas = fmas (t, y, varargin)
22934 bb452f84a299 Fix printing of solution stats in ode23/ode45 solvers (bug #49918). Carlo de Falco <carlo.defalco@polimi.it> parents: 22755 diff changeset	346 %! mas = [1, 0; 0, 1]; # Dummy mass matrix for tests
20549 e368ce72a844 maint: Use Octave coding conventions for ode* functions. Rik <rik@octave.org> parents: 20548 diff changeset	347 %!endfunction
20634 80e630b37ba1 maint: Remove unnecessary 'v' prefix before variables in ODE m-files. Rik <rik@octave.org> parents: 20631 diff changeset	348 %!function mas = fmsa (t, y, varargin)
22934 bb452f84a299 Fix printing of solution stats in ode23/ode45 solvers (bug #49918). Carlo de Falco <carlo.defalco@polimi.it> parents: 22755 diff changeset	349 %! mas = sparse ([1, 0; 0, 1]); # A sparse dummy matrix
20549 e368ce72a844 maint: Use Octave coding conventions for ode* functions. Rik <rik@octave.org> parents: 20548 diff changeset	350 %!endfunction
20634 80e630b37ba1 maint: Remove unnecessary 'v' prefix before variables in ODE m-files. Rik <rik@octave.org> parents: 20631 diff changeset	351 %!function out = fout (t, y, flag, varargin)
22639 7efa2d0e22c9 More Matlab-compatible implementation of OutputFcn. Clean up odeplot.m. Rik <rik@octave.org> parents: 22626 diff changeset	352 %! out = false;
7efa2d0e22c9 More Matlab-compatible implementation of OutputFcn. Clean up odeplot.m. Rik <rik@octave.org> parents: 22626 diff changeset	353 %! if (strcmp (flag, "init"))
7efa2d0e22c9 More Matlab-compatible implementation of OutputFcn. Clean up odeplot.m. Rik <rik@octave.org> parents: 22626 diff changeset	354 %! if (! isequal (size (t), [2, 1]))
7efa2d0e22c9 More Matlab-compatible implementation of OutputFcn. Clean up odeplot.m. Rik <rik@octave.org> parents: 22626 diff changeset	355 %! error ('fout: step "init"');
7efa2d0e22c9 More Matlab-compatible implementation of OutputFcn. Clean up odeplot.m. Rik <rik@octave.org> parents: 22626 diff changeset	356 %! endif
22626 869c02fde46c Further clean-up of ode functions. Rik <rik@octave.org> parents: 22625 diff changeset	357 %! elseif (isempty (flag))
22639 7efa2d0e22c9 More Matlab-compatible implementation of OutputFcn. Clean up odeplot.m. Rik <rik@octave.org> parents: 22626 diff changeset	358 %! if (! isequal (size (t), [1, 1]))
7efa2d0e22c9 More Matlab-compatible implementation of OutputFcn. Clean up odeplot.m. Rik <rik@octave.org> parents: 22626 diff changeset	359 %! error ('fout: step "calc"');
7efa2d0e22c9 More Matlab-compatible implementation of OutputFcn. Clean up odeplot.m. Rik <rik@octave.org> parents: 22626 diff changeset	360 %! endif
7efa2d0e22c9 More Matlab-compatible implementation of OutputFcn. Clean up odeplot.m. Rik <rik@octave.org> parents: 22626 diff changeset	361 %! elseif (strcmp (flag, "done"))
7efa2d0e22c9 More Matlab-compatible implementation of OutputFcn. Clean up odeplot.m. Rik <rik@octave.org> parents: 22626 diff changeset	362 %! if (! isempty (t))
7efa2d0e22c9 More Matlab-compatible implementation of OutputFcn. Clean up odeplot.m. Rik <rik@octave.org> parents: 22626 diff changeset	363 %! warning ('fout: step "done"');
7efa2d0e22c9 More Matlab-compatible implementation of OutputFcn. Clean up odeplot.m. Rik <rik@octave.org> parents: 22626 diff changeset	364 %! endif
22626 869c02fde46c Further clean-up of ode functions. Rik <rik@octave.org> parents: 22625 diff changeset	365 %! else
22639 7efa2d0e22c9 More Matlab-compatible implementation of OutputFcn. Clean up odeplot.m. Rik <rik@octave.org> parents: 22626 diff changeset	366 %! error ("fout: invalid flag <%s>", flag);
22626 869c02fde46c Further clean-up of ode functions. Rik <rik@octave.org> parents: 22625 diff changeset	367 %! endif
20549 e368ce72a844 maint: Use Octave coding conventions for ode* functions. Rik <rik@octave.org> parents: 20548 diff changeset	368 %!endfunction
20533 fcb792acab9b Moving ode45, odeset, odeget, and levenshtein from odepkg to core. jcorno <jacopo.corno@gmail.com> parents: diff changeset	369 %!
20549 e368ce72a844 maint: Use Octave coding conventions for ode* functions. Rik <rik@octave.org> parents: 20548 diff changeset	370 %!test # two output arguments
20634 80e630b37ba1 maint: Remove unnecessary 'v' prefix before variables in ODE m-files. Rik <rik@octave.org> parents: 20631 diff changeset	371 %! [t, y] = ode45 (@fpol, [0 2], [2 0]);
80e630b37ba1 maint: Remove unnecessary 'v' prefix before variables in ODE m-files. Rik <rik@octave.org> parents: 20631 diff changeset	372 %! assert ([t(end), y(end,:)], [2, fref], 1e-2);
20549 e368ce72a844 maint: Use Octave coding conventions for ode* functions. Rik <rik@octave.org> parents: 20548 diff changeset	373 %!test # not too many steps
20634 80e630b37ba1 maint: Remove unnecessary 'v' prefix before variables in ODE m-files. Rik <rik@octave.org> parents: 20631 diff changeset	374 %! [t, y] = ode45 (@fpol, [0 2], [2 0]);
80e630b37ba1 maint: Remove unnecessary 'v' prefix before variables in ODE m-files. Rik <rik@octave.org> parents: 20631 diff changeset	375 %! assert (size (t) < 20);
20549 e368ce72a844 maint: Use Octave coding conventions for ode* functions. Rik <rik@octave.org> parents: 20548 diff changeset	376 %!test # anonymous function instead of real function
22644 23941beb44cd Cleanup variable names, obsolete ChangeLog info in ode23, ode45. Rik <rik@octave.org> parents: 22643 diff changeset	377 %! fvdp = @(t,y) [y(2); (1 - y(1)^2) * y(2) - y(1)];
23941beb44cd Cleanup variable names, obsolete ChangeLog info in ode23, ode45. Rik <rik@octave.org> parents: 22643 diff changeset	378 %! [t, y] = ode45 (fvdp, [0 2], [2 0]);
20634 80e630b37ba1 maint: Remove unnecessary 'v' prefix before variables in ODE m-files. Rik <rik@octave.org> parents: 20631 diff changeset	379 %! assert ([t(end), y(end,:)], [2, fref], 1e-2);
20613 27c091f4b66d allow first argument in ode45 to be a string Carlo de Falco <carlo.defalco@polimi.it> parents: 20612 diff changeset	380 %!test # string instead of function
20634 80e630b37ba1 maint: Remove unnecessary 'v' prefix before variables in ODE m-files. Rik <rik@octave.org> parents: 20631 diff changeset	381 %! [t, y] = ode45 ("fpol", [0 2], [2 0]);
80e630b37ba1 maint: Remove unnecessary 'v' prefix before variables in ODE m-files. Rik <rik@octave.org> parents: 20631 diff changeset	382 %! assert ([t(end), y(end,:)], [2, fref], 1e-2);
20549 e368ce72a844 maint: Use Octave coding conventions for ode* functions. Rik <rik@octave.org> parents: 20548 diff changeset	383 %!test # extra input arguments passed through
20634 80e630b37ba1 maint: Remove unnecessary 'v' prefix before variables in ODE m-files. Rik <rik@octave.org> parents: 20631 diff changeset	384 %! [t, y] = ode45 (@fpol, [0 2], [2 0], 12, 13, "KL");
80e630b37ba1 maint: Remove unnecessary 'v' prefix before variables in ODE m-files. Rik <rik@octave.org> parents: 20631 diff changeset	385 %! assert ([t(end), y(end,:)], [2, fref], 1e-2);
20631 00caf63edcdf maint: Remove obsolete ODE options from odeset/odeget/ode45. Rik <rik@octave.org> parents: 20620 diff changeset	386 %!test # empty ODEOPT structure but extra input arguments
20634 80e630b37ba1 maint: Remove unnecessary 'v' prefix before variables in ODE m-files. Rik <rik@octave.org> parents: 20631 diff changeset	387 %! opt = odeset;
80e630b37ba1 maint: Remove unnecessary 'v' prefix before variables in ODE m-files. Rik <rik@octave.org> parents: 20631 diff changeset	388 %! [t, y] = ode45 (@fpol, [0 2], [2 0], opt, 12, 13, "KL");
80e630b37ba1 maint: Remove unnecessary 'v' prefix before variables in ODE m-files. Rik <rik@octave.org> parents: 20631 diff changeset	389 %! assert ([t(end), y(end,:)], [2, fref], 1e-2);
20549 e368ce72a844 maint: Use Octave coding conventions for ode* functions. Rik <rik@octave.org> parents: 20548 diff changeset	390 %!test # Solve another anonymous function below zero
e368ce72a844 maint: Use Octave coding conventions for ode* functions. Rik <rik@octave.org> parents: 20548 diff changeset	391 %! vref = [0, 14.77810590694212];
20634 80e630b37ba1 maint: Remove unnecessary 'v' prefix before variables in ODE m-files. Rik <rik@octave.org> parents: 20631 diff changeset	392 %! [t, y] = ode45 (@(t,y) y, [-2 0], 2);
80e630b37ba1 maint: Remove unnecessary 'v' prefix before variables in ODE m-files. Rik <rik@octave.org> parents: 20631 diff changeset	393 %! assert ([t(end), y(end,:)], vref, 1e-1);
20549 e368ce72a844 maint: Use Octave coding conventions for ode* functions. Rik <rik@octave.org> parents: 20548 diff changeset	394 %!test # InitialStep option
20634 80e630b37ba1 maint: Remove unnecessary 'v' prefix before variables in ODE m-files. Rik <rik@octave.org> parents: 20631 diff changeset	395 %! opt = odeset ("InitialStep", 1e-8);
80e630b37ba1 maint: Remove unnecessary 'v' prefix before variables in ODE m-files. Rik <rik@octave.org> parents: 20631 diff changeset	396 %! [t, y] = ode45 (@fpol, [0 0.2], [2 0], opt);
80e630b37ba1 maint: Remove unnecessary 'v' prefix before variables in ODE m-files. Rik <rik@octave.org> parents: 20631 diff changeset	397 %! assert ([t(2)-t(1)], [1e-8], 1e-9);
20549 e368ce72a844 maint: Use Octave coding conventions for ode* functions. Rik <rik@octave.org> parents: 20548 diff changeset	398 %!test # MaxStep option
20634 80e630b37ba1 maint: Remove unnecessary 'v' prefix before variables in ODE m-files. Rik <rik@octave.org> parents: 20631 diff changeset	399 %! opt = odeset ("MaxStep", 1e-3);
80e630b37ba1 maint: Remove unnecessary 'v' prefix before variables in ODE m-files. Rik <rik@octave.org> parents: 20631 diff changeset	400 %! sol = ode45 (@fpol, [0 0.2], [2 0], opt);
80e630b37ba1 maint: Remove unnecessary 'v' prefix before variables in ODE m-files. Rik <rik@octave.org> parents: 20631 diff changeset	401 %! assert ([sol.x(5)-sol.x(4)], [1e-3], 1e-3);
22659 04fc7e9c5f96 Disable Refine option of ode solvers for certain input/output combinations (partial fix bug #49408). Rik <rik@octave.org> parents: 22658 diff changeset	402 %!test # Solve with intermediate step
04fc7e9c5f96 Disable Refine option of ode solvers for certain input/output combinations (partial fix bug #49408). Rik <rik@octave.org> parents: 22658 diff changeset	403 %! [t, y] = ode45 (@fpol, [0 1 2], [2 0]);
28912 0de38a6ef693 maint: Use Octave convention of space after function name in scripts dir. Rik <rik@octave.org> parents: 28896 diff changeset	404 %! assert (any ((t-1) == 0));
22659 04fc7e9c5f96 Disable Refine option of ode solvers for certain input/output combinations (partial fix bug #49408). Rik <rik@octave.org> parents: 22658 diff changeset	405 %! assert ([t(end), y(end,:)], [2, fref], 1e-3);
20549 e368ce72a844 maint: Use Octave coding conventions for ode* functions. Rik <rik@octave.org> parents: 20548 diff changeset	406 %!test # Solve in backward direction starting at t=0
e368ce72a844 maint: Use Octave coding conventions for ode* functions. Rik <rik@octave.org> parents: 20548 diff changeset	407 %! vref = [-1.205364552835178, 0.951542399860817];
20634 80e630b37ba1 maint: Remove unnecessary 'v' prefix before variables in ODE m-files. Rik <rik@octave.org> parents: 20631 diff changeset	408 %! sol = ode45 (@fpol, [0 -2], [2 0]);
22654 bc61ed076549 Change orientation of output fields in struct returned from ode solvers (bug #49402). Sebastian Schöps <sebastian@schoeps.org> parents: 22647 diff changeset	409 %! assert ([sol.x(end); sol.y(:,end)], [-2; vref'], 1e-2);
20549 e368ce72a844 maint: Use Octave coding conventions for ode* functions. Rik <rik@octave.org> parents: 20548 diff changeset	410 %!test # Solve in backward direction starting at t=2
e368ce72a844 maint: Use Octave coding conventions for ode* functions. Rik <rik@octave.org> parents: 20548 diff changeset	411 %! vref = [-1.205364552835178, 0.951542399860817];
20634 80e630b37ba1 maint: Remove unnecessary 'v' prefix before variables in ODE m-files. Rik <rik@octave.org> parents: 20631 diff changeset	412 %! sol = ode45 (@fpol, [2 -2], fref);
22654 bc61ed076549 Change orientation of output fields in struct returned from ode solvers (bug #49402). Sebastian Schöps <sebastian@schoeps.org> parents: 22647 diff changeset	413 %! assert ([sol.x(end); sol.y(:,end)], [-2; vref'], 1e-2);
22659 04fc7e9c5f96 Disable Refine option of ode solvers for certain input/output combinations (partial fix bug #49408). Rik <rik@octave.org> parents: 22658 diff changeset	414 %!test # Solve in backward direction starting at t=2, with intermediate step
20549 e368ce72a844 maint: Use Octave coding conventions for ode* functions. Rik <rik@octave.org> parents: 20548 diff changeset	415 %! vref = [-1.205364552835178, 0.951542399860817];
22659 04fc7e9c5f96 Disable Refine option of ode solvers for certain input/output combinations (partial fix bug #49408). Rik <rik@octave.org> parents: 22658 diff changeset	416 %! [t, y] = ode45 (@fpol, [2 0 -2], fref);
28912 0de38a6ef693 maint: Use Octave convention of space after function name in scripts dir. Rik <rik@octave.org> parents: 28896 diff changeset	417 %! idx = find (y < 0, 1, "first") - 1;
22659 04fc7e9c5f96 Disable Refine option of ode solvers for certain input/output combinations (partial fix bug #49408). Rik <rik@octave.org> parents: 22658 diff changeset	418 %! assert ([t(idx), y(idx,:)], [0,2,0], 1e-2);
04fc7e9c5f96 Disable Refine option of ode solvers for certain input/output combinations (partial fix bug #49408). Rik <rik@octave.org> parents: 22658 diff changeset	419 %! assert ([t(end), y(end,:)], [-2, vref], 1e-2);
20549 e368ce72a844 maint: Use Octave coding conventions for ode* functions. Rik <rik@octave.org> parents: 20548 diff changeset	420 %!test # Solve another anonymous function in backward direction
e368ce72a844 maint: Use Octave coding conventions for ode* functions. Rik <rik@octave.org> parents: 20548 diff changeset	421 %! vref = [-1, 0.367879437558975];
20634 80e630b37ba1 maint: Remove unnecessary 'v' prefix before variables in ODE m-files. Rik <rik@octave.org> parents: 20631 diff changeset	422 %! sol = ode45 (@(t,y) y, [0 -1], 1);
22654 bc61ed076549 Change orientation of output fields in struct returned from ode solvers (bug #49402). Sebastian Schöps <sebastian@schoeps.org> parents: 22647 diff changeset	423 %! assert ([sol.x(end); sol.y(:,end)], vref', 1e-3);
20549 e368ce72a844 maint: Use Octave coding conventions for ode* functions. Rik <rik@octave.org> parents: 20548 diff changeset	424 %!test # Solve another anonymous function below zero
e368ce72a844 maint: Use Octave coding conventions for ode* functions. Rik <rik@octave.org> parents: 20548 diff changeset	425 %! vref = [0, 14.77810590694212];
20634 80e630b37ba1 maint: Remove unnecessary 'v' prefix before variables in ODE m-files. Rik <rik@octave.org> parents: 20631 diff changeset	426 %! sol = ode45 (@(t,y) y, [-2 0], 2);
22654 bc61ed076549 Change orientation of output fields in struct returned from ode solvers (bug #49402). Sebastian Schöps <sebastian@schoeps.org> parents: 22647 diff changeset	427 %! assert ([sol.x(end); sol.y(:,end)], vref', 1e-3);
20549 e368ce72a844 maint: Use Octave coding conventions for ode* functions. Rik <rik@octave.org> parents: 20548 diff changeset	428 %!test # Solve in backward direction starting at t=0 with MaxStep option
e368ce72a844 maint: Use Octave coding conventions for ode* functions. Rik <rik@octave.org> parents: 20548 diff changeset	429 %! vref = [-1.205364552835178, 0.951542399860817];
20634 80e630b37ba1 maint: Remove unnecessary 'v' prefix before variables in ODE m-files. Rik <rik@octave.org> parents: 20631 diff changeset	430 %! opt = odeset ("MaxStep", 1e-3);
80e630b37ba1 maint: Remove unnecessary 'v' prefix before variables in ODE m-files. Rik <rik@octave.org> parents: 20631 diff changeset	431 %! sol = ode45 (@fpol, [0 -2], [2 0], opt);
80e630b37ba1 maint: Remove unnecessary 'v' prefix before variables in ODE m-files. Rik <rik@octave.org> parents: 20631 diff changeset	432 %! assert ([abs(sol.x(8)-sol.x(7))], [1e-3], 1e-3);
22654 bc61ed076549 Change orientation of output fields in struct returned from ode solvers (bug #49402). Sebastian Schöps <sebastian@schoeps.org> parents: 22647 diff changeset	433 %! assert ([sol.x(end); sol.y(:,end)], [-2; vref'], 1e-3);
20549 e368ce72a844 maint: Use Octave coding conventions for ode* functions. Rik <rik@octave.org> parents: 20548 diff changeset	434 %!test # AbsTol option
20634 80e630b37ba1 maint: Remove unnecessary 'v' prefix before variables in ODE m-files. Rik <rik@octave.org> parents: 20631 diff changeset	435 %! opt = odeset ("AbsTol", 1e-5);
80e630b37ba1 maint: Remove unnecessary 'v' prefix before variables in ODE m-files. Rik <rik@octave.org> parents: 20631 diff changeset	436 %! sol = ode45 (@fpol, [0 2], [2 0], opt);
22654 bc61ed076549 Change orientation of output fields in struct returned from ode solvers (bug #49402). Sebastian Schöps <sebastian@schoeps.org> parents: 22647 diff changeset	437 %! assert ([sol.x(end); sol.y(:,end)], [2; fref'], 1e-3);
20549 e368ce72a844 maint: Use Octave coding conventions for ode* functions. Rik <rik@octave.org> parents: 20548 diff changeset	438 %!test # AbsTol and RelTol option
20634 80e630b37ba1 maint: Remove unnecessary 'v' prefix before variables in ODE m-files. Rik <rik@octave.org> parents: 20631 diff changeset	439 %! opt = odeset ("AbsTol", 1e-8, "RelTol", 1e-8);
80e630b37ba1 maint: Remove unnecessary 'v' prefix before variables in ODE m-files. Rik <rik@octave.org> parents: 20631 diff changeset	440 %! sol = ode45 (@fpol, [0 2], [2 0], opt);
22654 bc61ed076549 Change orientation of output fields in struct returned from ode solvers (bug #49402). Sebastian Schöps <sebastian@schoeps.org> parents: 22647 diff changeset	441 %! assert ([sol.x(end); sol.y(:,end)], [2; fref'], 1e-3);
20549 e368ce72a844 maint: Use Octave coding conventions for ode* functions. Rik <rik@octave.org> parents: 20548 diff changeset	442 %!test # RelTol and NormControl option -- higher accuracy
20634 80e630b37ba1 maint: Remove unnecessary 'v' prefix before variables in ODE m-files. Rik <rik@octave.org> parents: 20631 diff changeset	443 %! opt = odeset ("RelTol", 1e-8, "NormControl", "on");
80e630b37ba1 maint: Remove unnecessary 'v' prefix before variables in ODE m-files. Rik <rik@octave.org> parents: 20631 diff changeset	444 %! sol = ode45 (@fpol, [0 2], [2 0], opt);
22654 bc61ed076549 Change orientation of output fields in struct returned from ode solvers (bug #49402). Sebastian Schöps <sebastian@schoeps.org> parents: 22647 diff changeset	445 %! assert ([sol.x(end); sol.y(:,end)], [2; fref'], 1e-5);
20549 e368ce72a844 maint: Use Octave coding conventions for ode* functions. Rik <rik@octave.org> parents: 20548 diff changeset	446 %!test # Keeps initial values while integrating
20634 80e630b37ba1 maint: Remove unnecessary 'v' prefix before variables in ODE m-files. Rik <rik@octave.org> parents: 20631 diff changeset	447 %! opt = odeset ("NonNegative", 2);
80e630b37ba1 maint: Remove unnecessary 'v' prefix before variables in ODE m-files. Rik <rik@octave.org> parents: 20631 diff changeset	448 %! sol = ode45 (@fpol, [0 2], [2 0], opt);
22654 bc61ed076549 Change orientation of output fields in struct returned from ode solvers (bug #49402). Sebastian Schöps <sebastian@schoeps.org> parents: 22647 diff changeset	449 %! assert ([sol.x(end); sol.y(:,end)], [2; 2; 0], 0.5);
20549 e368ce72a844 maint: Use Octave coding conventions for ode* functions. Rik <rik@octave.org> parents: 20548 diff changeset	450 %!test # Details of OutputSel and Refine can't be tested
20634 80e630b37ba1 maint: Remove unnecessary 'v' prefix before variables in ODE m-files. Rik <rik@octave.org> parents: 20631 diff changeset	451 %! opt = odeset ("OutputFcn", @fout, "OutputSel", 1, "Refine", 5);
80e630b37ba1 maint: Remove unnecessary 'v' prefix before variables in ODE m-files. Rik <rik@octave.org> parents: 20631 diff changeset	452 %! sol = ode45 (@fpol, [0 2], [2 0], opt);
80e630b37ba1 maint: Remove unnecessary 'v' prefix before variables in ODE m-files. Rik <rik@octave.org> parents: 20631 diff changeset	453 %!test # Stats must add further elements in sol
80e630b37ba1 maint: Remove unnecessary 'v' prefix before variables in ODE m-files. Rik <rik@octave.org> parents: 20631 diff changeset	454 %! opt = odeset ("Stats", "on");
22934 bb452f84a299 Fix printing of solution stats in ode23/ode45 solvers (bug #49918). Carlo de Falco <carlo.defalco@polimi.it> parents: 22755 diff changeset	455 %! stat_str = evalc ("sol = ode45 (@fpol, [0 2], [2 0], opt);");
bb452f84a299 Fix printing of solution stats in ode23/ode45 solvers (bug #49918). Carlo de Falco <carlo.defalco@polimi.it> parents: 22755 diff changeset	456 %! assert (strncmp (stat_str, "Number of successful steps:", 27));
20634 80e630b37ba1 maint: Remove unnecessary 'v' prefix before variables in ODE m-files. Rik <rik@octave.org> parents: 20631 diff changeset	457 %! assert (isfield (sol, "stats"));
80e630b37ba1 maint: Remove unnecessary 'v' prefix before variables in ODE m-files. Rik <rik@octave.org> parents: 20631 diff changeset	458 %! assert (isfield (sol.stats, "nsteps"));
80e630b37ba1 maint: Remove unnecessary 'v' prefix before variables in ODE m-files. Rik <rik@octave.org> parents: 20631 diff changeset	459 %!test # Events option add further elements in sol
80e630b37ba1 maint: Remove unnecessary 'v' prefix before variables in ODE m-files. Rik <rik@octave.org> parents: 20631 diff changeset	460 %! opt = odeset ("Events", @feve);
80e630b37ba1 maint: Remove unnecessary 'v' prefix before variables in ODE m-files. Rik <rik@octave.org> parents: 20631 diff changeset	461 %! sol = ode45 (@fpol, [0 10], [2 0], opt);
80e630b37ba1 maint: Remove unnecessary 'v' prefix before variables in ODE m-files. Rik <rik@octave.org> parents: 20631 diff changeset	462 %! assert (isfield (sol, "ie"));
80e630b37ba1 maint: Remove unnecessary 'v' prefix before variables in ODE m-files. Rik <rik@octave.org> parents: 20631 diff changeset	463 %! assert (sol.ie(1), 2);
80e630b37ba1 maint: Remove unnecessary 'v' prefix before variables in ODE m-files. Rik <rik@octave.org> parents: 20631 diff changeset	464 %! assert (isfield (sol, "xe"));
80e630b37ba1 maint: Remove unnecessary 'v' prefix before variables in ODE m-files. Rik <rik@octave.org> parents: 20631 diff changeset	465 %! assert (isfield (sol, "ye"));
20549 e368ce72a844 maint: Use Octave coding conventions for ode* functions. Rik <rik@octave.org> parents: 20548 diff changeset	466 %!test # Events option, now stop integration
21443 acd6e203031d Alter BIST tests stop emitting warnings during runtests invocation. Rik <rik@octave.org> parents: 20908 diff changeset	467 %! warning ("off", "integrate_adaptive:unexpected_termination", "local");
20634 80e630b37ba1 maint: Remove unnecessary 'v' prefix before variables in ODE m-files. Rik <rik@octave.org> parents: 20631 diff changeset	468 %! opt = odeset ("Events", @fevn, "NormControl", "on");
80e630b37ba1 maint: Remove unnecessary 'v' prefix before variables in ODE m-files. Rik <rik@octave.org> parents: 20631 diff changeset	469 %! sol = ode45 (@fpol, [0 10], [2 0], opt);
80e630b37ba1 maint: Remove unnecessary 'v' prefix before variables in ODE m-files. Rik <rik@octave.org> parents: 20631 diff changeset	470 %! assert ([sol.ie, sol.xe, sol.ye],
20553 45151de7423f maint: Clean up implementations of ode45.m, odeget.m, odeset.m. Rik <rik@octave.org> parents: 20549 diff changeset	471 %! [2.0, 2.496110, -0.830550, -2.677589], 6e-1);
20549 e368ce72a844 maint: Use Octave coding conventions for ode* functions. Rik <rik@octave.org> parents: 20548 diff changeset	472 %!test # Events option, five output arguments
21443 acd6e203031d Alter BIST tests stop emitting warnings during runtests invocation. Rik <rik@octave.org> parents: 20908 diff changeset	473 %! warning ("off", "integrate_adaptive:unexpected_termination", "local");
20634 80e630b37ba1 maint: Remove unnecessary 'v' prefix before variables in ODE m-files. Rik <rik@octave.org> parents: 20631 diff changeset	474 %! opt = odeset ("Events", @fevn, "NormControl", "on");
80e630b37ba1 maint: Remove unnecessary 'v' prefix before variables in ODE m-files. Rik <rik@octave.org> parents: 20631 diff changeset	475 %! [t, y, vxe, ye, vie] = ode45 (@fpol, [0 10], [2 0], opt);
80e630b37ba1 maint: Remove unnecessary 'v' prefix before variables in ODE m-files. Rik <rik@octave.org> parents: 20631 diff changeset	476 %! assert ([vie, vxe, ye],
20553 45151de7423f maint: Clean up implementations of ode45.m, odeget.m, odeset.m. Rik <rik@octave.org> parents: 20549 diff changeset	477 %! [2.0, 2.496110, -0.830550, -2.677589], 6e-1);
20549 e368ce72a844 maint: Use Octave coding conventions for ode* functions. Rik <rik@octave.org> parents: 20548 diff changeset	478 %!test # Mass option as function
20634 80e630b37ba1 maint: Remove unnecessary 'v' prefix before variables in ODE m-files. Rik <rik@octave.org> parents: 20631 diff changeset	479 %! opt = odeset ("Mass", @fmas);
80e630b37ba1 maint: Remove unnecessary 'v' prefix before variables in ODE m-files. Rik <rik@octave.org> parents: 20631 diff changeset	480 %! sol = ode45 (@fpol, [0 2], [2 0], opt);
22654 bc61ed076549 Change orientation of output fields in struct returned from ode solvers (bug #49402). Sebastian Schöps <sebastian@schoeps.org> parents: 22647 diff changeset	481 %! assert ([sol.x(end); sol.y(:,end)], [2; fref'], 1e-3);
20549 e368ce72a844 maint: Use Octave coding conventions for ode* functions. Rik <rik@octave.org> parents: 20548 diff changeset	482 %!test # Mass option as matrix
20634 80e630b37ba1 maint: Remove unnecessary 'v' prefix before variables in ODE m-files. Rik <rik@octave.org> parents: 20631 diff changeset	483 %! opt = odeset ("Mass", eye (2,2));
80e630b37ba1 maint: Remove unnecessary 'v' prefix before variables in ODE m-files. Rik <rik@octave.org> parents: 20631 diff changeset	484 %! sol = ode45 (@fpol, [0 2], [2 0], opt);
22654 bc61ed076549 Change orientation of output fields in struct returned from ode solvers (bug #49402). Sebastian Schöps <sebastian@schoeps.org> parents: 22647 diff changeset	485 %! assert ([sol.x(end); sol.y(:,end)], [2; fref'], 1e-3);
20549 e368ce72a844 maint: Use Octave coding conventions for ode* functions. Rik <rik@octave.org> parents: 20548 diff changeset	486 %!test # Mass option as sparse matrix
20634 80e630b37ba1 maint: Remove unnecessary 'v' prefix before variables in ODE m-files. Rik <rik@octave.org> parents: 20631 diff changeset	487 %! opt = odeset ("Mass", sparse (eye (2,2)));
80e630b37ba1 maint: Remove unnecessary 'v' prefix before variables in ODE m-files. Rik <rik@octave.org> parents: 20631 diff changeset	488 %! sol = ode45 (@fpol, [0 2], [2 0], opt);
22654 bc61ed076549 Change orientation of output fields in struct returned from ode solvers (bug #49402). Sebastian Schöps <sebastian@schoeps.org> parents: 22647 diff changeset	489 %! assert ([sol.x(end); sol.y(:,end)], [2; fref'], 1e-3);
20549 e368ce72a844 maint: Use Octave coding conventions for ode* functions. Rik <rik@octave.org> parents: 20548 diff changeset	490 %!test # Mass option as function and sparse matrix
20634 80e630b37ba1 maint: Remove unnecessary 'v' prefix before variables in ODE m-files. Rik <rik@octave.org> parents: 20631 diff changeset	491 %! opt = odeset ("Mass", @fmsa);
80e630b37ba1 maint: Remove unnecessary 'v' prefix before variables in ODE m-files. Rik <rik@octave.org> parents: 20631 diff changeset	492 %! sol = ode45 (@fpol, [0 2], [2 0], opt);
22654 bc61ed076549 Change orientation of output fields in struct returned from ode solvers (bug #49402). Sebastian Schöps <sebastian@schoeps.org> parents: 22647 diff changeset	493 %! assert ([sol.x(end); sol.y(:,end)], [2; fref'], 1e-3);
20549 e368ce72a844 maint: Use Octave coding conventions for ode* functions. Rik <rik@octave.org> parents: 20548 diff changeset	494 %!test # Mass option as function and MStateDependence
20634 80e630b37ba1 maint: Remove unnecessary 'v' prefix before variables in ODE m-files. Rik <rik@octave.org> parents: 20631 diff changeset	495 %! opt = odeset ("Mass", @fmas, "MStateDependence", "strong");
80e630b37ba1 maint: Remove unnecessary 'v' prefix before variables in ODE m-files. Rik <rik@octave.org> parents: 20631 diff changeset	496 %! sol = ode45 (@fpol, [0 2], [2 0], opt);
22654 bc61ed076549 Change orientation of output fields in struct returned from ode solvers (bug #49402). Sebastian Schöps <sebastian@schoeps.org> parents: 22647 diff changeset	497 %! assert ([sol.x(end); sol.y(:,end)], [2; fref'], 1e-3);
21443 acd6e203031d Alter BIST tests stop emitting warnings during runtests invocation. Rik <rik@octave.org> parents: 20908 diff changeset	498
22660 7eb3f8ec1aed Remove FIXME comment about missing tests for ignored ODE options. Carlo de Falco <carlo.defalco@polimi.it> parents: 22659 diff changeset	499 ## Note: The following options have no effect on this solver
7eb3f8ec1aed Remove FIXME comment about missing tests for ignored ODE options. Carlo de Falco <carlo.defalco@polimi.it> parents: 22659 diff changeset	500 ## therefore it makes no sense to test them here:
7eb3f8ec1aed Remove FIXME comment about missing tests for ignored ODE options. Carlo de Falco <carlo.defalco@polimi.it> parents: 22659 diff changeset	501 ##
7eb3f8ec1aed Remove FIXME comment about missing tests for ignored ODE options. Carlo de Falco <carlo.defalco@polimi.it> parents: 22659 diff changeset	502 ## "BDF"
7eb3f8ec1aed Remove FIXME comment about missing tests for ignored ODE options. Carlo de Falco <carlo.defalco@polimi.it> parents: 22659 diff changeset	503 ## "InitialSlope"
7eb3f8ec1aed Remove FIXME comment about missing tests for ignored ODE options. Carlo de Falco <carlo.defalco@polimi.it> parents: 22659 diff changeset	504 ## "JPattern"
7eb3f8ec1aed Remove FIXME comment about missing tests for ignored ODE options. Carlo de Falco <carlo.defalco@polimi.it> parents: 22659 diff changeset	505 ## "Jacobian"
7eb3f8ec1aed Remove FIXME comment about missing tests for ignored ODE options. Carlo de Falco <carlo.defalco@polimi.it> parents: 22659 diff changeset	506 ## "MassSingular"
7eb3f8ec1aed Remove FIXME comment about missing tests for ignored ODE options. Carlo de Falco <carlo.defalco@polimi.it> parents: 22659 diff changeset	507 ## "MaxOrder"
7eb3f8ec1aed Remove FIXME comment about missing tests for ignored ODE options. Carlo de Falco <carlo.defalco@polimi.it> parents: 22659 diff changeset	508 ## "MvPattern"
7eb3f8ec1aed Remove FIXME comment about missing tests for ignored ODE options. Carlo de Falco <carlo.defalco@polimi.it> parents: 22659 diff changeset	509 ## "Vectorized"
22626 869c02fde46c Further clean-up of ode functions. Rik <rik@octave.org> parents: 22625 diff changeset	510
22655 6b134d294d61 ode solvers: use ordinary transpose instead of Hermitian conjugate (bug #49410). Carlo de Falco <carlo.defalco@polimi.it> parents: 22654 diff changeset	511 %!test # Check that imaginary part of solution does not get inverted
6b134d294d61 ode solvers: use ordinary transpose instead of Hermitian conjugate (bug #49410). Carlo de Falco <carlo.defalco@polimi.it> parents: 22654 diff changeset	512 %! sol = ode45 (@(x,y) 1, [0 1], 1i);
28929 9e43deb9bfc3 maint: Use semicolon after assert statement inside %!test blocks. Rik <rik@octave.org> parents: 28912 diff changeset	513 %! assert (imag (sol.y), ones (size (sol.y)));
22655 6b134d294d61 ode solvers: use ordinary transpose instead of Hermitian conjugate (bug #49410). Carlo de Falco <carlo.defalco@polimi.it> parents: 22654 diff changeset	514 %! [x, y] = ode45 (@(x,y) 1, [0 1], 1i);
28929 9e43deb9bfc3 maint: Use semicolon after assert statement inside %!test blocks. Rik <rik@octave.org> parents: 28912 diff changeset	515 %! assert (imag (y), ones (size (y)));
22655 6b134d294d61 ode solvers: use ordinary transpose instead of Hermitian conjugate (bug #49410). Carlo de Falco <carlo.defalco@polimi.it> parents: 22654 diff changeset	516
28896 90fea9cc9caa test: Add expected error message <Invalid call> to BIST tests for nargin. Rik <rik@octave.org> parents: 28892 diff changeset	517 %!error <Invalid call> ode45 ()
90fea9cc9caa test: Add expected error message <Invalid call> to BIST tests for nargin. Rik <rik@octave.org> parents: 28892 diff changeset	518 %!error <Invalid call> ode45 (1)
90fea9cc9caa test: Add expected error message <Invalid call> to BIST tests for nargin. Rik <rik@octave.org> parents: 28892 diff changeset	519 %!error <Invalid call> ode45 (1,2)
22626 869c02fde46c Further clean-up of ode functions. Rik <rik@octave.org> parents: 22625 diff changeset	520 %!error <TRANGE must be a numeric> ode45 (@fpol, {[0 25]}, [3 15 1])
869c02fde46c Further clean-up of ode functions. Rik <rik@octave.org> parents: 22625 diff changeset	521 %!error <TRANGE must be a .* vector> ode45 (@fpol, [0 25; 25 0], [3 15 1])
869c02fde46c Further clean-up of ode functions. Rik <rik@octave.org> parents: 22625 diff changeset	522 %!error <TRANGE must contain at least 2 elements> ode45 (@fpol, [1], [3 15 1])
869c02fde46c Further clean-up of ode functions. Rik <rik@octave.org> parents: 22625 diff changeset	523 %!error <invalid time span> ode45 (@fpol, [1 1], [3 15 1])
869c02fde46c Further clean-up of ode functions. Rik <rik@octave.org> parents: 22625 diff changeset	524 %!error <INIT must be a numeric> ode45 (@fpol, [0 25], {[3 15 1]})
869c02fde46c Further clean-up of ode functions. Rik <rik@octave.org> parents: 22625 diff changeset	525 %!error <INIT must be a .* vector> ode45 (@fpol, [0 25], [3 15 1; 3 15 1])
869c02fde46c Further clean-up of ode functions. Rik <rik@octave.org> parents: 22625 diff changeset	526 %!error <FUN must be a valid function handle> ode45 (1, [0 25], [3 15 1])