1 //////////////////////////////////////////////////////////////////////// 2 // 3 // Copyright (C) 1996-2021 The Octave Project Developers 4 // 5 // See the file COPYRIGHT.md in the top-level directory of this 6 // distribution or <https://octave.org/copyright/>. 7 // 8 // This file is part of Octave. 9 // 10 // Octave is free software: you can redistribute it and/or modify it 11 // under the terms of the GNU General Public License as published by 12 // the Free Software Foundation, either version 3 of the License, or 13 // (at your option) any later version. 14 // 15 // Octave is distributed in the hope that it will be useful, but 16 // WITHOUT ANY WARRANTY; without even the implied warranty of 17 // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the 18 // GNU General Public License for more details. 19 // 20 // You should have received a copy of the GNU General Public License 21 // along with Octave; see the file COPYING. If not, see 22 // <https://www.gnu.org/licenses/>. 23 // 24 //////////////////////////////////////////////////////////////////////// 25 26 #if ! defined (octave_DASPK_h) 27 #define octave_DASPK_h 1 28 29 #include "octave-config.h" 30 31 #include <string> 32 33 #include "Array.h" 34 #include "DASPK-opts.h" 35 36 class Matrix; 37 38 class 39 OCTAVE_API 40 DASPK : public DAE, public DASPK_options 41 { 42 public: 43 DASPK(void)44 DASPK (void) 45 : DAE (), DASPK_options (), initialized (false), liw (0), lrw (0), 46 info (), iwork (), rwork (), abs_tol (), rel_tol () { } 47 DASPK(const ColumnVector & s,double tm,DAEFunc & f)48 DASPK (const ColumnVector& s, double tm, DAEFunc& f) 49 : DAE (s, tm, f), DASPK_options (), initialized (false), liw (0), 50 lrw (0), info (), iwork (), rwork (), abs_tol (), rel_tol () { } 51 DASPK(const ColumnVector & s,const ColumnVector & deriv,double tm,DAEFunc & f)52 DASPK (const ColumnVector& s, const ColumnVector& deriv, 53 double tm, DAEFunc& f) 54 : DAE (s, deriv, tm, f), DASPK_options (), initialized (false), 55 liw (0), lrw (0), info (), iwork (), rwork (), abs_tol (), 56 rel_tol () { } 57 58 ~DASPK (void) = default; 59 60 ColumnVector do_integrate (double t); 61 62 Matrix do_integrate (const ColumnVector& tout); 63 64 Matrix do_integrate (const ColumnVector& tout, const ColumnVector& tcrit); 65 66 Matrix integrate (const ColumnVector& tout, Matrix& xdot_out); 67 68 Matrix integrate (const ColumnVector& tout, Matrix& xdot_out, 69 const ColumnVector& tcrit); 70 71 std::string error_message (void) const; 72 73 private: 74 75 bool initialized; 76 77 octave_f77_int_type liw; 78 octave_f77_int_type lrw; 79 80 Array<octave_f77_int_type> info; 81 Array<octave_f77_int_type> iwork; 82 83 Array<double> rwork; 84 85 Array<double> abs_tol; 86 Array<double> rel_tol; 87 }; 88 89 #endif 90