00001
00002
00003
00004
00005
00006
00007
00008
00009
00010
00011
00012
00013
00014
00015
00016
00017
00018
00019
00020
00021
00022
00023
00024
00025
00026
00027
00028
00029
00030
00031
00032
00033
00034
00035
00036 #ifndef SIBS_H
00037 #define SIBS_H
00038
00039 #include "ODESolver.H"
00040
00041
00042
00043 namespace Foam
00044 {
00045
00046
00047
00048
00049
00050 class SIBS
00051 :
00052 public ODESolver
00053 {
00054
00055
00056 static const label kMaxX_ = 7, iMaxX_ = kMaxX_ + 1;
00057 static const label nSeq_[iMaxX_];
00058
00059 static const scalar safe1, safe2, redMax, redMin, scaleMX;
00060
00061 mutable scalarField a_;
00062 mutable Matrix<scalar> alpha_;
00063 mutable Matrix<scalar> d_p_;
00064 mutable scalarField x_p_;
00065 mutable scalarField err_;
00066
00067 mutable scalarField yTemp_;
00068 mutable scalarField ySeq_;
00069 mutable scalarField yErr_;
00070 mutable scalarField dfdx_;
00071 mutable Matrix<scalar> dfdy_;
00072
00073 mutable label first_, kMax_, kOpt_;
00074 mutable scalar epsOld_, xNew_;
00075
00076
00077 void SIMPR
00078 (
00079 const ODE& ode,
00080 const scalar xStart,
00081 const scalarField& y,
00082 const scalarField& dydx,
00083 const scalarField& dfdx,
00084 const Matrix<scalar>& dfdy,
00085 const scalar deltaX,
00086 const label nSteps,
00087 scalarField& yEnd
00088 ) const;
00089
00090
00091 void polyExtrapolate
00092 (
00093 const label iest,
00094 const scalar xest,
00095 const scalarField& yest,
00096 scalarField& yz,
00097 scalarField& dy,
00098 scalarField& x_p,
00099 Matrix<scalar>& d_p
00100 ) const;
00101
00102
00103 public:
00104
00105
00106 TypeName("SIBS");
00107
00108
00109
00110
00111
00112 SIBS(const ODE& ode);
00113
00114
00115
00116
00117 void solve
00118 (
00119 const ODE& ode,
00120 scalar& x,
00121 scalarField& y,
00122 scalarField& dydx,
00123 const scalar eps,
00124 const scalarField& yScale,
00125 const scalar hTry,
00126 scalar& hDid,
00127 scalar& hNext
00128 ) const;
00129 };
00130
00131
00132
00133
00134 }
00135
00136
00137
00138 #endif
00139
00140