transpose.hh

Go to the documentation of this file.

 
#ifndef transpose_hh
#define transpose_hh
 
#include "operator/matrix.hh"
 
namespace concepts {
 
  // ************************************************************* Transpose **
 
  template<class F>
  class Transpose : public Matrix<F> {
  public:
    typedef typename Realtype<F>::type r_type;
    typedef typename Cmplxtype<F>::type c_type;
 
    Transpose(Matrix<F>& a) : Matrix<F>(a.dimY(), a.dimX()), a_(a) { }
    virtual F operator()(const uint i, const uint j) const { return a_(j, i); }
    virtual F& operator()(const uint i, const uint j) { return a_(j, i); }
    virtual void operator()(const Function<r_type>& fncY,
                            Function<F>& fncX);
    virtual void operator()(const Function<c_type>& fncY,
                            Function<c_type>& fncX);
    template<class H, class I>
    void operator()(const Vector<H>& fncY, Vector<I>& fncX) {
      a_.transpMult(fncY, fncX);
    }
 
    virtual void transpMult(const Vector<r_type>& fncY, Vector<F>& fncX) {
      a_(fncY, fncX);
    }
    virtual void transpMult(const Vector<c_type>& fncY,
                            Vector<c_type>& fncX)
    {
      a_(fncY, fncX);
    }
  protected:
    std::ostream& info(std::ostream& os) const;
  private:
    Matrix<F>& a_;
  };
 
} // namespace concepts
 
#endif // transpose_hh