00001
00002
00003
00004
00005
00006
00007
00008
00009
00010
00011
00012
00013
00014
00015
00016
00017
00018
00019
00020
00021
00022
00023 #if !defined (octave_LSODE_h)
00024 #define octave_LSODE_h 1
00025
00026 #if defined (__GNUG__) && defined (USE_PRAGMA_INTERFACE_IMPLEMENTATION)
00027 #pragma interface
00028 #endif
00029
00030 #include <cfloat>
00031 #include <cmath>
00032
00033 #include "LSODE-opts.h"
00034
00035 class
00036 LSODE : public ODE, public LSODE_options
00037 {
00038 public:
00039
00040 LSODE (void) : ODE (), LSODE_options (), initialized (false) { }
00041
00042 LSODE (const ColumnVector& s, double tm, const ODEFunc& f)
00043 : ODE (s, tm, f), LSODE_options (), initialized (false) { }
00044
00045 ~LSODE (void) { }
00046
00047 ColumnVector do_integrate (double t);
00048
00049 Matrix do_integrate (const ColumnVector& tout);
00050
00051 Matrix do_integrate (const ColumnVector& tout, const ColumnVector& tcrit);
00052
00053 std::string error_message (void) const;
00054
00055 private:
00056
00057 bool initialized;
00058
00059 int method_flag;
00060 int maxord;
00061 int itask;
00062 int iopt;
00063 int itol;
00064
00065 int liw;
00066 int lrw;
00067
00068 Array<int> iwork;
00069 Array<double> rwork;
00070
00071 double rel_tol;
00072
00073 Array<double> abs_tol;
00074
00075 double *px;
00076 double *pabs_tol;
00077 int *piwork;
00078 double *prwork;
00079 };
00080
00081 #endif
00082
00083
00084
00085
00086
00087