---

itsqr.h

Go to the documentation of this file.

// FILE itsqr.h 


#ifndef INCLUDED_ITSQR_H
#define INCLUDED_ITSQR_H


#include "itsbase.h"
#include "taylctrl.h"
#include "messgs.h"


typedef void (*QRFactor)( MATRIX & Q, const MATRIX & A );


class ITS_QR : public ITS_BASE
{       
  
public:                                                                 
  
  ITS_QR( int order = TAYLOR_CTRL::GetMaxOrder() );
  
  
  void Init( const PtrODENumeric ODE );  
  
  
  virtual void Compute( INTERVAL_VECTOR & Ytight, 
                        INTERVAL_VECTOR & LocErr,
                        const INTERVAL & h, int order,
                        PtrODENumeric ODE, 
                        PtrODESolver Solver );
  
  
  void SetQRFunction( QRFactor f ); 
  
  
  ~ITS_QR();
  

private:
  
  QRFactor QRFactorFunction;
  
  INTERVAL_MATRIX  AjInv;  
  
  MATRIX TmpMat1;
  
  void ResizeSpace( int n );
  
};



inline ITS_QR:: ITS_QR( int order ) : 
  ITS_BASE( order, MSG_ITS_QR_METHOD ) 
{ 
  assert( order>0 && order<=TAYLOR_CTRL::GetMaxOrder() );
}


inline void ITS_QR :: SetQRFunction( QRFactor f )
{
  assert( f!=NULL );
  QRFactorFunction = f;
}


inline ITS_QR :: ~ITS_QR() 
{ 
  ;
}


#endif 

---