Teko_ReorderedLinearOp.hpp

// @HEADER
// *****************************************************************************
//      Teko: A package for block and physics based preconditioning
//
// Copyright 2010 NTESS and the Teko contributors.
// SPDX-License-Identifier: BSD-3-Clause
// *****************************************************************************
// @HEADER

#ifndef __Teko_ReorderedLinearOp_hpp__
#define __Teko_ReorderedLinearOp_hpp__

#include "Teko_Utilities.hpp"
#include "Teko_ImplicitLinearOp.hpp"
#include "Teko_BlockedReordering.hpp"

namespace Teko {

class ReorderedLinearOp : public ImplicitLinearOp {
 public:
  ReorderedLinearOp(const Teuchos::RCP<const BlockReorderManager> &mgr,
                    const Teuchos::RCP<Thyra::LinearOpBase<double> > &blockedOp);

  Teuchos::RCP<const BlockReorderManager> getReorderManager() const { return mgr_; }

  Teko::ModifiableLinearOp getBlockedOp() const { return blockedOp_; }

  virtual VectorSpace range() const;

  virtual VectorSpace domain() const;

  virtual void implicitApply(const MultiVector &x, MultiVector &y, const double alpha = 1.0,
                             const double beta = 0.0) const;

  virtual void describe(Teuchos::FancyOStream &out_arg,
                        const Teuchos::EVerbosityLevel verbLevel) const;

 private:
  VectorSpace range_;
  VectorSpace domain_;
  Teuchos::RCP<const BlockReorderManager> mgr_;
  Teuchos::RCP<Thyra::LinearOpBase<double> > blockedOp_;
};

}  // end namespace Teko

#endif