blockOperator.hh

Go to the documentation of this file.

 
#ifndef blockOperator_hh
#define blockOperator_hh
 
#include <string>
#include <memory>
#include "basics/typedefs.hh"
#include "operator/compositions.hh"
#include "vectorial/vectorial.hh"
 
namespace vectorial {
 
  // ********************************************************* BlockOperator **
 
  template<class F>
  class BlockOperator : public concepts::Operator<F> {
  public:
    BlockOperator(concepts::Operator<F>* A = 0,
                  concepts::Operator<F>* B = 0,
                  concepts::Operator<F>* C = 0,
                  concepts::Operator<F>* D = 0);
    BlockOperator(uint dimX, uint dimY,
                  concepts::Operator<F>* A = 0,
                  concepts::Operator<F>* B = 0,
                  concepts::Operator<F>* C = 0,
                  concepts::Operator<F>* D = 0);
 
    virtual void operator()(const concepts::Function<F>& fncY,
                            concepts::Function<F>& fncX);
    void operator()(const concepts::Vector<F>& fncY,
                    concepts::Vector<F>& fncX);
  protected:
    virtual std::ostream& info(std::ostream& os) const;
  private:
    uint bdimX_[2], bdimY_[2];
    concepts::Operator<F>* A_;
    concepts::Operator<F>* B_;
    concepts::Operator<F>* C_;
    concepts::Operator<F>* D_;
  };
 
} // namespace vectorial
 
#endif // blockOperator_hh