/usr/include/CLHEP/GenericFunctions/ButcherTableau.hh is in libclhep-dev 2.1.4.1+dfsg-1.
This file is owned by root:root, with mode 0o644.
The actual contents of the file can be viewed below.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 | #ifndef _ButcherTableau_h_
#define _ButcherTableau_h_
// This class defines a Butcher Tableau, which completely specifies
// a Runge-Kutte integration scheme. Butcher Tableau are described
// in Numerical Methods for Ordinary Differential Equations, John
// Wiley & sons, West Sussex England.
//
// General form is :
//
// c|A
// ---
// |b^T
//
// where A is a matrix and b, c are column vectors.
//
// The Butcher Tableau Class presents itself as an empty structure
// that the user has to fill up. One can blithely fill write into
// any element of A, b, or c. Space is automatically allocated.
#include <vector>
#include <string>
namespace Genfun {
class ButcherTableau {
public:
// Constructor:
inline ButcherTableau(const std::string &name, unsigned int order);
// Returns the name:
inline const std::string & name() const;
// Returns the order:
inline unsigned int order() const;
// Returns the number of steps:
inline unsigned int nSteps() const;
// Write access to elements:
inline double & A(unsigned int i, unsigned int j);
inline double & b(unsigned int i);
inline double & c(unsigned int i);
// Read access to elements (inline for speed)
inline const double & A(unsigned int i, unsigned int j) const;
inline const double & b(unsigned int i) const;
inline const double & c(unsigned int i) const;
private:
std::vector< std::vector<double> > _A;
std::vector<double> _b;
std::vector<double> _c;
std::string _name;
unsigned int _order;
};
class EulerTableau: public ButcherTableau {
// Constructor:
public:
inline EulerTableau();
};
class MidpointTableau: public ButcherTableau {
// Constructor:
public:
inline MidpointTableau();
};
class TrapezoidTableau: public ButcherTableau {
// Constructor:
public:
inline TrapezoidTableau();
};
class RK31Tableau: public ButcherTableau {
// Constructor:
public:
inline RK31Tableau();
};
class RK32Tableau: public ButcherTableau {
// Constructor:
public:
inline RK32Tableau();
};
class ClassicalRungeKuttaTableau: public ButcherTableau {
// Constructor:
public:
inline ClassicalRungeKuttaTableau();
};
class ThreeEighthsRuleTableau: public ButcherTableau {
// Constructor:
public:
inline ThreeEighthsRuleTableau();
};
}
inline std::ostream & operator << (std::ostream & o, const Genfun::ButcherTableau & b);
#include "CLHEP/GenericFunctions/ButcherTableau.icc"
#endif
|