← Back to branch summary

~ubuntu-branches/debian/sid/octave3.0/sid

« back to all changes in this revision

Viewing changes to liboctave/base-de.h

  • Committer: Bazaar Package Importer
  • Author(s): Rafael Laboissiere
  • Date: 2007-12-23 16:04:15 UTC
  • Revision ID: james.westby@ubuntu.com-20071223160415-n4gk468dihy22e9v
Tags: upstream-3.0.0
Import upstream version 3.0.0

Show diffs side-by-side

added added

removed removed

Lines of Context:
liboctave/base-de.h
 
1
/*
 
2
 
 
3
Copyright (C) 1996, 1997, 2000, 2002, 2003, 2005, 2007 John W. Eaton
 
4
 
 
5
This file is part of Octave.
 
6
 
 
7
Octave is free software; you can redistribute it and/or modify it
 
8
under the terms of the GNU General Public License as published by the
 
9
Free Software Foundation; either version 3 of the License, or (at your
 
10
option) any later version.
 
11
 
 
12
Octave is distributed in the hope that it will be useful, but WITHOUT
 
13
ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
 
14
FITNESS FOR A PARTICULAR PURPOSE.  See the GNU General Public License
 
15
for more details.
 
16
 
 
17
You should have received a copy of the GNU General Public License
 
18
along with Octave; see the file COPYING.  If not, see
 
19
<http://www.gnu.org/licenses/>.
 
20
 
 
21
*/
 
22
 
 
23
#if !defined (octave_base_de_h)
 
24
#define octave_base_de_h 1
 
25
 
 
26
#include <string>
 
27
 
 
28
#include "dColVector.h"
 
29
#include "dMatrix.h"
 
30
 
 
31
class
 
32
base_diff_eqn
 
33
{
 
34
public:
 
35
 
 
36
  base_diff_eqn (void)
 
37
    : x (), t (0.0), stop_time (0.0), stop_time_set (false),
 
38
      restart (true), integration_error (false), istate (0) { } 
 
39
 
 
40
  base_diff_eqn (const ColumnVector& xx, double tt)
 
41
    : x (xx), t (tt), stop_time (0.0), stop_time_set (false),
 
42
      restart (true), integration_error (false), istate (0) { }
 
43
 
 
44
  base_diff_eqn (const base_diff_eqn& a)
 
45
    : x (a.x), t (a.t), stop_time (0.0), stop_time_set (false),
 
46
      restart (true), integration_error (false), istate (0) { }
 
47
 
 
48
  virtual ~base_diff_eqn (void) { }
 
49
 
 
50
  base_diff_eqn& operator = (const base_diff_eqn& a)
 
51
    {
 
52
      if (this != &a)
 
53
        {
 
54
          x = a.x;
 
55
          t = a.t;
 
56
          stop_time = a.stop_time;
 
57
          stop_time_set = a.stop_time_set;
 
58
          restart = a.restart;
 
59
          integration_error = a.integration_error;
 
60
          istate = a.istate;
 
61
        }
 
62
 
 
63
      return *this;
 
64
    }
 
65
 
 
66
  void initialize (const ColumnVector& x0, double t0)
 
67
    {
 
68
      x = x0;
 
69
      t = t0;
 
70
      integration_error = false;
 
71
      istate = 0;
 
72
      force_restart ();
 
73
    }
 
74
 
 
75
  octave_idx_type size (void) const { return x.capacity (); }
 
76
 
 
77
  ColumnVector state (void) const { return x; }
 
78
 
 
79
  double time (void) const { return t; }
 
80
 
 
81
  void set_stop_time (double tt)
 
82
    {
 
83
      stop_time_set = true;
 
84
      stop_time = tt;
 
85
      force_restart ();
 
86
    }
 
87
 
 
88
  void clear_stop_time (void)
 
89
    {
 
90
      stop_time_set = false;
 
91
      force_restart ();
 
92
    }
 
93
 
 
94
  virtual void force_restart (void) { restart = true; }
 
95
 
 
96
  bool integration_ok (void) const { return ! integration_error; }
 
97
 
 
98
  octave_idx_type integration_state (void) const { return istate; }
 
99
 
 
100
  virtual std::string error_message (void) const = 0;
 
101
 
 
102
protected:
 
103
 
 
104
  ColumnVector x;
 
105
 
 
106
  double t;
 
107
 
 
108
  double stop_time;
 
109
 
 
110
  bool stop_time_set;
 
111
 
 
112
  bool restart;
 
113
 
 
114
  bool integration_error;
 
115
 
 
116
  octave_idx_type istate;
 
117
};
 
118
 
 
119
#endif
 
120
 
 
121
/*
 
122
;;; Local Variables: ***
 
123
;;; mode: C++ ***
 
124
;;; End: ***
 
125
*/