Mercurial > octave-nkf
diff liboctave/ODE.h @ 3:9a4c07481e61
[project @ 1993-08-08 01:20:23 by jwe]
Initial revision
author | jwe |
---|---|
date | Sun, 08 Aug 1993 01:21:46 +0000 |
parents | |
children | 780cbbc57b7c |
line wrap: on
line diff
--- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/liboctave/ODE.h Sun Aug 08 01:21:46 1993 +0000 @@ -0,0 +1,107 @@ +// ODE.h -*- C++ -*- +/* + +Copyright (C) 1992, 1993 John W. Eaton + +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 2, 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, write to the Free +Software Foundation, 675 Mass Ave, Cambridge, MA 02139, USA. + +*/ + +#if !defined (_ODE_h) +#define _ODE_h 1 + +#ifdef __GNUG__ +#pragma interface +#endif + +#include <iostream.h> +#include "ODEFunc.h" +#include "Matrix.h" + +class ODE : public ODEFunc +{ +public: + + ODE (void); + + ODE (int n); + + ODE (const ColumnVector& state, double time, const ODEFunc& f); + + virtual ~ODE (void); + + virtual int size (void) const; + virtual ColumnVector state (void) const; + virtual double time (void) const; + + virtual void force_restart (void); + virtual void initialize (const ColumnVector& x, double t); + virtual void set_stop_time (double t); + virtual void clear_stop_time (void); + + virtual ColumnVector integrate (double t); + + void integrate (int nsteps, double tstep, ostream& s); + + Matrix integrate (const ColumnVector& tout); + Matrix integrate (const ColumnVector& tout, const ColumnVector& tcrit); + +protected: + +/* + * Some of this is probably too closely related to LSODE, but hey, + * this is just a first attempt... + */ + + int n; + double t; + ColumnVector x; + + double absolute_tolerance; + double relative_tolerance; + + double stop_time; + int stop_time_set; + +private: + + int restart; + int method_flag; + int *iwork; + double *rwork; + int istate; + int itol; + int itask; + int iopt; + int liw; + int lrw; + + friend int lsode_f (int *neq, double *t, double *y, double *ydot); + + friend int lsode_j (int *neq, double *t, double *y, int *ml, int *mu, + double *pd, int *nrowpd); + +}; + +#endif + +/* +;;; Local Variables: *** +;;; mode: C++ *** +;;; page-delimiter: "^/\\*" *** +;;; End: *** +*/