matfile/denseOWrapper.hh

Go to the documentation of this file.

#ifndef DENSEOWRAPPER_HH_
#define DENSEOWRAPPER_HH_


using namespace matfile;

namespace concepts{

class DenseOWrapper{
  
  //Methoden die der Writer bzw. ArrayDenseWrapper erwartet
public: 
  
  
  typedef ArrayDenseWrapper Type; 
  
  template<class T>
  DenseOWrapper(concepts::Matrix<T>& matrix);
  
  virtual ~DenseOWrapper(){};
  
  bool isComplex() const{
  return isComplex_;
  }
   
  Real* realData(uint i, uint j) const{
    if(!isComplex_)
      return &((*matrix_r)(i,j));
    return &((*matrix_c)(i,j).real());
  }
  
  Real* complexData(uint i, uint j) const{
    return &((*matrix_c)(i,j).imag());
  }
  
  const MATtype realType() const {
    return realType_;
  }
  
  const MATtype imagType() const {
    return cmplxType_;
  }
  
  void size(miINT32_t& rows, miINT32_t& cols) const {
    rows = rows_;
    cols = cols_;
  }
    
  private:
    
    bool isComplex_;
    
    //Number of rows, colums
    miINT32_t rows_, cols_;

    //matrices
    concepts::DenseMatrix<Real>* matrix_r;
    concepts::DenseMatrix<Cmplx>* matrix_c;
   
    //Daten fuer matfile
    MATtype realType_, cmplxType_; 
};


template<class T>
DenseOWrapper::DenseOWrapper(concepts::Matrix< T >& matrix){
  
  //Setze standard Informationen
  isComplex_ = (NumTraits<T>::isComplex == 1);
  rows_ = matrix.dimX();
  cols_ = matrix.dimY();
  
  //Setze Informationen fuer matfile
  realType_ = NumTraits<T>::MatType;
  cmplxType_ = NumTraits<T>::MatType;
  
  //Setze Pointer auf die Matrix
  if(!isComplex_)
    matrix_r = dynamic_cast<DenseMatrix<Real>*>(&matrix);
  else
    matrix_c = dynamic_cast<DenseMatrix<Cmplx>*>(&matrix);
  }

//DenseOWrapper::~DenseOWrapper(){
//}

}
#endif