Mercurial > hg > octave-nkf
view libinterp/octave-value/ov-oncleanup.h @ 20811:a22d8a2eb0e5
fix adaptive strategy in ode solvers.
* script/ode/ode45.m: remove unused option OutputSave
* script/ode/private/integrate_adaptive.m: rewrite algorithm
to be more compatible.
* script/ode/private/runge_kutta_45_dorpri.m: use kahan summation
for time increment.
author | Carlo de Falco <carlo.defalco@polimi.it> |
---|---|
date | Sun, 11 Oct 2015 18:44:58 +0200 |
parents | 01895a6acdde |
children |
line wrap: on
line source
/* Copyright (C) 2010-2015 VZLU Prague This file is part of Octave. Octave is free software; you can redistribute it and/or modify it under the terms of the GNU General Public License as published by the Free Software Foundation; either version 3 of the License, or (at your option) any later version. Octave is distributed in the hope that it will be useful, but WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License for more details. You should have received a copy of the GNU General Public License along with Octave; see the file COPYING. If not, see <http://www.gnu.org/licenses/>. */ #include <iosfwd> #include "ov-base.h" #include "ov-struct.h" #include "ov.h" static void gripe_internal (void) { error ("onCleanup: internal error: cloning nonempty object"); } class octave_oncleanup : public octave_base_value { public: octave_oncleanup (void) : fcn () { } octave_oncleanup (const octave_value& fcn); octave_base_value *clone (void) const { if (fcn.is_defined ()) gripe_internal (); return empty_clone (); } octave_base_value *empty_clone (void) const { return new octave_oncleanup (); } ~octave_oncleanup (void); bool is_defined (void) const { return true; } bool is_object (void) const { return true; } // do we want this? octave_map map_value (void) const { return scalar_map_value (); } octave_scalar_map scalar_map_value (void) const; dim_vector dims (void) const { static dim_vector dv (1, 1); return dv; } bool save_ascii (std::ostream& os); bool load_ascii (std::istream& is); bool save_binary (std::ostream& os, bool& save_as_floats); bool load_binary (std::istream& is, bool swap, oct_mach_info::float_format fmt); bool save_hdf5 (octave_hdf5_id loc_id, const char *name, bool save_as_floats); bool load_hdf5 (octave_hdf5_id loc_id, const char *name); void print (std::ostream& os, bool pr_as_read_syntax = false); void print_raw (std::ostream& os, bool pr_as_read_syntax = false) const; private: DECLARE_OV_TYPEID_FUNCTIONS_AND_DATA protected: octave_value fcn; };