---

odesol.h

Go to the documentation of this file.

// FILE odesol.h 


#ifndef INCLUDED_SOLUTION
#define INCLUDED_SOLUTION


#include <assert.h>
#include "vnodecnf.h"


class SOLUTION 
{

public:
  
  SOLUTION( int n );

  int GetDim() const;
  
  void SetTcur( double t );
  
  void SetTightEncl( const INTERVAL_VECTOR & Y );
  
  void SetInitEncl( const INTERVAL_VECTOR & Y );
  
  double GetTcur() const;
  
  double GetTprev() const; 
  
  const INTERVAL_VECTOR & GetTightEncl() const; 
  
  const INTERVAL_VECTOR & GetInitEncl() const; 

  SOLUTION & operator = ( const SOLUTION & S);
  
  virtual ~SOLUTION(); 

private:

  double tCur, tPrev;
  
  //point, tCur
  INTERVAL_VECTOR Ytight; 
  
  INTERVAL_VECTOR Yinit; 
  
};


typedef SOLUTION*  PtrSolution;  


inline SOLUTION :: SOLUTION( int n ) : tCur(0), tPrev(0)
{ 
  assert( n>0 ); 
  
  Resize( Ytight, n ); 
  Resize( Yinit, n ); 
}


inline void SOLUTION :: SetTcur( double t )  
{ 
  tPrev = tCur;
  tCur = t;
}

inline void SOLUTION :: SetTightEncl( const INTERVAL_VECTOR & Y ) 
{ 
  assert( GetDim()==Dimension(Y) );
  Ytight = Y; 
}


inline void SOLUTION :: SetInitEncl( const INTERVAL_VECTOR & Y ) 
{ 
  assert( GetDim()==Dimension(Y) );
  Yinit = Y; 
}


inline double SOLUTION :: GetTcur() const    
{ 
  return tCur; 
}


inline double SOLUTION :: GetTprev() const    
{ 
  return tPrev; 
}


inline const INTERVAL_VECTOR & SOLUTION :: GetTightEncl() const 
{ 
  assert(Dimension(Ytight)!=0);
  return Ytight; 
}


inline const INTERVAL_VECTOR & SOLUTION :: GetInitEncl() const 
{ 
  return Yinit; 
}


inline SOLUTION :: ~SOLUTION()
{
  ;
}


inline int SOLUTION :: GetDim() const
{
  return Dimension(Ytight);
}


#endif

---