Mercurial > octave-nkf
view liboctave/DAE.h @ 100:197bd6c63baa
[project @ 1993-09-13 03:06:15 by jwe]
(syl): Use complex_matrix_value(), not ComplexMatrixLoad().
author | jwe |
---|---|
date | Mon, 13 Sep 1993 03:06:15 +0000 |
parents | 9a4c07481e61 |
children | 780cbbc57b7c |
line wrap: on
line source
// DAE.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 (_DAE_h) #define _DAE_h 1 #ifdef __GNUG__ #pragma interface #endif #include <iostream.h> #include "ODE.h" #include "DAEFunc.h" #include "Matrix.h" #include "f77-uscore.h" #ifndef Vector #define Vector ColumnVector #endif class DAE : public ODE, public DAEFunc { public: DAE (void); DAE (int); DAE (Vector& x, double time, DAEFunc& f); DAE (Vector& x, Vector& xdot, double time, DAEFunc& f); ~DAE (void); Vector deriv (void); virtual void initialize (Vector& x, double t); virtual void initialize (Vector& x, Vector& xdot, double t); Vector integrate (double t); Matrix integrate (const Vector& tout, Matrix& xdot_out); Matrix integrate (const Vector& tout, Matrix& xdot_out, const Vector& tcrit); protected: /* * Some of this is probably too closely related to DASSL, but hey, * this is just a first attempt... */ Vector xdot; private: int restart; int liw; int lrw; int idid; int *info; int *iwork; double *rwork; friend int ddassl_j (double *time, double *state, double *deriv, double *pd, double *cj, double *rpar, int *ipar); friend int ddassl_f (double *time, double *state, double *deriv, double *delta, int *ires, double *rpar, int *ipar); }; #endif