class ODECrawler: public ZGlass
ODE integrator using adaptive step functions from GSL. X corresponds to time (scalar). Y corresponds to state vector at a given time. Final state can be obtained from mY. Works in double precision. Trajectory is stored in ODEstorage object with at least mStoreDx between consequtive points. Maximum of mStoreMax intermediate points are stored. If mStoreMax = 0 no results are stored. If mStoreMax < 0 all points are stored (default). Start and end points (at x1 and x2) are always stored, unless mStoreMax=0. Example glass Moonraker. Macro moonraker.C.
Function Members (Methods)
public:
protected:
| virtual void | ZGlass::AdEndarkenment() |
| virtual void | ZGlass::AdEnlightenment() |
| virtual void | ZGlass::AdUnfoldment() |
| ZMIR* | ZGlass::assert_MIR_presence(const Exc_t& header, int what = 0) |
| ZGlass* | assert_odemaster(const Exc_t& eh) |
| virtual void | TObject::DoError(int level, const char* location, const char* fmt, va_list va) const |
| ZMIR* | ZGlass::get_MIR() |
| void | init_integration(Bool_t call_ode_start) |
| void | Integrate() |
| void | TObject::MakeZombie() |
| void | ZGlass::queen_check_in(ZGlass* l) |
| virtual void | ZGlass::reference_all() |
| virtual void | ZGlass::reference_links() |
| virtual Int_t | ZGlass::remove_references_to(ZGlass* lens) |
| void | ZGlass::set_link_or_die(ZGlass*& link, ZGlass* new_val, FID_t fid) |
| ZMIR* | ZGlass::suggest_MIR_presence(const Exc_t& header, int what = 0) |
| virtual void | ZGlass::unreference_all() |
| virtual void | ZGlass::unreference_links() |
| void | ZGlass::warn_caller(const TString& warning) |
private:
| void | _gsl_alloc() |
| void | _gsl_free() |
| void | _init() |
Data Members
public:
| enum StepFunc_e { | SF_RK2 | |
| SF_RK4 | ||
| SF_RKF45 | ||
| SF_RKCK | ||
| SF_RK8PD | ||
| SF_RK2Imp | ||
| SF_RK4Imp | ||
| SF_BSImp | ||
| SF_Gear1 | ||
| SF_Gear2 | ||
| }; | ||
| enum ZGlass::Bits_e { | kFixedNameBit | |
| kDyingBit | ||
| }; | ||
| enum ZGlass::MirComponents_e { | MC_Any | |
| MC_IsFlare | ||
| MC_IsBeam | ||
| MC_HasResultReq | ||
| MC_IsDetached | ||
| }; | ||
| enum TObject::EStatusBits { | kCanDelete | |
| kMustCleanup | ||
| kObjInCanvas | ||
| kIsReferenced | ||
| kHasUUID | ||
| kCannotPick | ||
| kNoContextMenu | ||
| kInvalidObject | ||
| }; | ||
| enum TObject::[unnamed] { | kIsOnHeap | |
| kNotDeleted | ||
| kZombie | ||
| kBitMask | ||
| kSingleKey | ||
| kOverwrite | ||
| kWriteDelete | ||
| }; |
protected:
| Bool_t | ZGlass::bAcceptRefs | X{GS} 7 BoolOut() |
| Bool_t | ZGlass::bMIRActive | X{GS} 7 BoolOut(-join=>1) |
| Bool_t | bStoreH | X{GS} 7 Bool() |
| Bool_t | ZGlass::bUseDispList | X{GST} 7 Bool(-join=>1) |
| Bool_t | ZGlass::bUseNameStack | X{GS} 7 Bool() |
| Double_t | mEpsAbs | X{GS} 7 Value(-range=>[0,1000], -join=>1) |
| Double_t | mEpsRel | X{GS} 7 Value(-range=>[0,1]) |
| UShort_t | ZGlass::mEyeRefCount | ! |
| Double_t | mFacDer | X{GS} 7 Value(-range=>[0,100]) |
| Double_t | mFacVal | X{GS} 7 Value(-range=>[0,100], -join=>1) |
| UShort_t | ZGlass::mFireRefCount | ! X{G} 7 ValOut(-width=>4) |
| UInt_t | ZGlass::mGlassBits | X{g} |
| ZLink<ZMirFilter> | ZGlass::mGuard | X{E} L{} |
| Double_t | mH | X{G} 7 ValOut(-join=>1) |
| Double_t | mH1 | X{GS} 7 Value(-join=>1) |
| Double_t | mHmin | X{GS} 7 Value() |
| Int_t | mMaxSteps | X{GS} 7 Value() |
| UShort_t | ZGlass::mMoonRefCount | ! X{G} 7 ValOut(-width=>4) |
| Int_t | mN | X{G} 7 ValOut() |
| TString | ZGlass::mName | X{RGE} 7 Textor() |
| ZLink<ZGlass> | mODEMaster | X{GS} L{a} |
| ODECrawler::StepFunc_e | mPrevStepFunc | ! |
| ZQueen* | ZGlass::mQueen | ! X{g} |
| GMutex | ZGlass::mReadMutex | ! |
| UShort_t | ZGlass::mRefCount | ! X{G} 7 ValOut(-width=>4, -join=>1) |
| Saturn* | ZGlass::mSaturn | ! X{g} |
| ID_t | ZGlass::mSaturnID | X{G} 7 ValOut(-range=>[0,MAX_ID,1], -width=>10) |
| TimeStamp_t | ZGlass::mStampReqTring | ! X{GS} TimeStamp of last change that requires retriangulation |
| Int_t | mStepChanged | ! X{G} 7 ValOut() |
| ODECrawler::StepFunc_e | mStepFunc | X{GS} 7 PhonyEnum() |
| Int_t | mStepOK | ! X{G} 7 ValOut(-join=>1) |
| ODEStorage* | mStorage | X{g} |
| Double_t | mStoreDx | X{GS} 7 Value(-range=>[0,1e12]) |
| Int_t | mStoreMax | X{GS} 7 Value(-join=>1) |
| Int_t | mStored | ! X{G} 7 ValOut(-join=>1) |
| UShort_t | ZGlass::mSunRefCount | ! X{G} 7 ValOut(-width=>4, -join=>1) |
| TimeStamp_t | ZGlass::mTimeStamp | ! X{GS} TimeStamp of last change |
| TString | ZGlass::mTitle | X{RGE} 7 Textor() |
| Double_t | mX1 | X{GS} 7 Value(-join=>1) |
| Double_t | mX2 | X{GS} 7 Value() |
| vector<Double_t> | mY | ! X{r} |
| set<NameChangeCB*>* | ZGlass::pspNameChangeCB | ! |
| set<RayAbsorber*>* | ZGlass::pspRayAbsorber | ! |
| static GledNS::LinkMemberInfo* | ZGlass::sap_Guard_lmi | |
| static GledNS::ClassInfo* | sap_ODECrawler_ci | |
| static GledNS::LinkMemberInfo* | sap_ODEMaster_lmi | |
| static GledNS::ClassInfo* | ZGlass::sap_ZGlass_ci |
private:
| Bool_t | m_crawling | ! |
| gsl_odeiv_control* | m_gsl_control | ! |
| gsl_odeiv_evolve* | m_gsl_evolve | ! |
| gsl_odeiv_step* | m_gsl_step | ! |
| gsl_odeiv_system* | m_gsl_system | ! |
| ODECrawlerMaster* | m_true_master | ! |
Class Charts
Function documentation
void _init()
ODECrawler(const Text_t* n = "ODECrawler", const Text_t* t = 0)
const gsl_odeiv_step_type* s_get_step_func(ODECrawler::StepFunc_e sf)
void _gsl_alloc()
void _gsl_free()
void AdvanceXLimits(Double_t delta_x = 0)
Advance x limits forward, so that X1 becomes X2 and new integration interval is delta_x. If delta_x is 0, length of the interval is preserved.
void SetStorage(ODEStorage* s)
void DetachStorage()
void init_integration(Bool_t call_ode_start)
Initialize storage arrays, clear integration state. If call_ode_start is true, call ODEStart() in master to get initial parameters. m_true_master must be set.
void Crawl(Bool_t call_ode_start = true)
Integrate with given master / parameters. If called via a MIR it is executed in a broadcasted detached thread. The object is locked for the duration of the execution.
void ChangeOrderInPlace(Int_t order)
ODECrawler(const Text_t* n = "ODECrawler", const Text_t* t = 0)