operator/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