MATH_GlobalMultiLevelsPolynomialToeplitzMatrix.h

Go to the documentation of this file.

#ifndef MATH_GlobalMultiLevelsPolynomialToeplitzMatrix_H
#define MATH_GlobalMultiLevelsPolynomialToeplitzMatrix_H

#include "MATH_MultiLevelsFFTToeplitzMatrix.h"



DEFINE_SPTR(MATH_GlobalMultiLevelsPolynomialToeplitzMatrix);

class MATH_GlobalMultiLevelsPolynomialToeplitzMatrix : public MATH_MultiLevelsFFTToeplitzMatrix  {
    SP_OBJECT(MATH_GlobalMultiLevelsPolynomialToeplitzMatrix);
    
public:
    
    
    
private:
  
    
    
protected: 
    
    // METHOD
    
    // CONSTRUCTORS 
  
    MATH_GlobalMultiLevelsPolynomialToeplitzMatrix(void);
    
    
    
    // DESTRUCTORS 
    
  
    virtual ~MATH_GlobalMultiLevelsPolynomialToeplitzMatrix(void);   
    
    
public:
    // NEW
    inline static SP::MATH_GlobalMultiLevelsPolynomialToeplitzMatrix New() {
        SP::MATH_GlobalMultiLevelsPolynomialToeplitzMatrix ret(new MATH_GlobalMultiLevelsPolynomialToeplitzMatrix(),
                                                  MATH_GlobalMultiLevelsPolynomialToeplitzMatrix::Delete()); 
        ret->setThis(ret);
        return ret;
    };

    virtual SP::MATH_MultiLevelsToeplitzMatrix NewInstance() const {
        return New();
    }

 
    virtual void copy(const MATH_MultiLevelsToeplitzMatrix& mat) {
        MATH_MultiLevelsToeplitzMatrix::copy(mat);
    }
    

    // SET methods
 
public:

    // GET methods 

   
  

  

    // diagonalization methods
    //=========================
protected:
    virtual void buildSpectralVectorProjector(FFTW_FFT& proj) const {
        MATH_MultiLevelsFFTToeplitzMatrix::buildSpectralVectorProjector(proj);
    }
    virtual void buildSpectralSpace();
    virtual void recoverCanonicalSpace();
    
    
    virtual void buildColumnLeadingDimensions();
    
    virtual void buildRowLeadingDimensions();
    
    //product in spectral space methods
    //==================================
protected:
    virtual void vectorProductInSpectralSpace(const FFTW_FFT& A,FFTW_FFT& Z) const {
        Z.multiply(A);
    }
    
    virtual void projectIntoSpectralSpace(const tUIndex& N, const  tReal* X) const;
  
    virtual void recoverFromSpectralSpace(const tUIndex& N,tReal* X) const;
   
    
};
#endif