// $Header: /cvs/gled-1.2/Numerica/Glasses/ODECrawler.h,v 1.4 2005/01/19 15:36:21 matevz Exp $

#ifndef Numerica_ODECrawler_H
#define Numerica_ODECrawler_H

#include <TVectorF.h>
#include <TVectorD.h>
#include <Glasses/Eventor.h>

class ODECrawlerMaster {
public:
  virtual ~ODECrawlerMaster() {}
  virtual UInt_t ODEOrder() = 0;
  virtual void ODEDerivatives(const Double_t x, const TVectorD& y, TVectorD& d) = 0;
  virtual void ODEStart(TVectorD& v, Double_t& x1, Double_t& x2) = 0;

  ClassDef(ODECrawlerMaster,1)
};

class ODECrawler : public Eventor {

private:
  Double_t hTINY, hSAFETY, hPGROW, hPSHRNK, hERRCON;	//!
  ODECrawlerMaster*	hTrueMaster;			//!

  void _init();

protected:
  ZGlass*	mODEMaster;	// X{gS} L{}

  UInt_t	mGuessesOK;	//! X{g}  7 ValOut(-join=>1)
  UInt_t	mGuessesBad;	//! X{g}  7 ValOut()
  UInt_t	mStored;	//! X{g}  7 ValOut(-join=>1)
  UInt_t	mMaxSteps;	//  X{gS} 7 Value()
  UInt_t	mStoreMax;	//  X{gS} 7 Value(-join=>1)
  Double_t	mStoreDx;	//  X{gS} 7 Value()
  TVectorF*	mXStored;	//! X{g}
  TVectorF*	mYStored;	//! X{g} ROOT now allows streaming ... check syntax
  
  TVectorD	mY;		//! X{r}
  UInt_t	mN;		// X{g}  7 ValOut(-join=>1)
  Double_t	mAcc;		// X{gS} 7 Value(-range=>[0,1])
  Double_t	mX1;		// X{gS} 7 Value(-join=>1)
  Double_t	mX2;		// X{gS} 7 Value()
  Double_t	mH1;		// X{gS} 7 Value(-join=>1)
  Double_t	mHmin;		// X{gS} 7 Value()
  
  Int_t		Rkqs(TVectorD& y, TVectorD& dydx, Double_t& x, Double_t htry,
		     TVectorD& yscal, Double_t& hdid, Double_t& hnext);
  void 		Rkck(TVectorD& y, TVectorD& dydx, Double_t x, Double_t h,
		     TVectorD& yout, TVectorD& yerr);
  void		Crawl();

public:
  ODECrawler(Text_t* n="ODECrawler", Text_t* t=0) :
    Eventor(n,t), mODEMaster(0) { _init(); }
  ODECrawler(ZGlass* m, const Text_t* n="ODECrawler", const Text_t* t=0) :
    Eventor(n,t), mODEMaster(m) { _init(); }
  virtual ~ODECrawler();

  void	DetachStored() { mXStored = mYStored = 0; }

  // virtuals
  virtual Operator::Arg* PreDance(Operator::Arg* op_arg=0);
  virtual void Operate(Operator::Arg* op_arg) throw(Operator::Exception);


#include "ODECrawler.h7"
  ClassDef(ODECrawler, 1)
}; // endclass ODECrawler

GlassIODef(ODECrawler);

#endif