Teko_LU2x2InverseOp.hpp

Go to the documentation of this file.

// @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_LU2x2InverseOp_hpp__
#define __Teko_LU2x2InverseOp_hpp__

#include "Teko_Utilities.hpp"
#include "Teko_BlockImplicitLinearOp.hpp"

namespace Teko {

class LU2x2InverseOp : public BlockImplicitLinearOp {
 public:
  LU2x2InverseOp(const BlockedLinearOp &A, const LinearOp &invA00, const LinearOp &invS);

  LU2x2InverseOp(const BlockedLinearOp &A, const LinearOp &hatInvA00, const LinearOp &tildeInvA00,
                 const LinearOp &invS);



  virtual VectorSpace range() const { return productRange_; }

  virtual VectorSpace domain() const { return productDomain_; }

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

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

 protected:
  using BlockImplicitLinearOp::implicitApply;

  // fundamental operators to use
  const BlockedLinearOp A_;     
  const LinearOp hatInvA00_;    
  const LinearOp tildeInvA00_;  
  const LinearOp invS_;         

  // some blocks of A
  const LinearOp A10_;  
  const LinearOp A01_;  

  Teuchos::RCP<const Thyra::ProductVectorSpaceBase<double> >
      productRange_;  
  Teuchos::RCP<const Thyra::ProductVectorSpaceBase<double> >
      productDomain_;  

 private:
  // hide me!
  LU2x2InverseOp();
  LU2x2InverseOp(const LU2x2InverseOp &);
};

inline LinearOp createLU2x2InverseOp(BlockedLinearOp &A, LinearOp &invA00, LinearOp &invS) {
  return Teuchos::rcp(new LU2x2InverseOp(A, invA00, invS));
}

inline LinearOp createLU2x2InverseOp(BlockedLinearOp &A, LinearOp &invA00, LinearOp &invS,
                                     const std::string &str) {
  Teuchos::RCP<Thyra::LinearOpBase<double> > result =
      Teuchos::rcp(new LU2x2InverseOp(A, invA00, invS));
  result->setObjectLabel(str);

  return result;
}

inline LinearOp createLU2x2InverseOp(BlockedLinearOp &A, LinearOp &hatInvA00, LinearOp &tildeInvA00,
                                     LinearOp &invS) {
  return Teuchos::rcp(new LU2x2InverseOp(A, hatInvA00, tildeInvA00, invS));
}

inline LinearOp createLU2x2InverseOp(BlockedLinearOp &A, LinearOp &hatInvA00, LinearOp &tildeInvA00,
                                     LinearOp &invS, const std::string &str) {
  Teuchos::RCP<Thyra::LinearOpBase<double> > result =
      Teuchos::rcp(new LU2x2InverseOp(A, hatInvA00, tildeInvA00, invS));
  result->setObjectLabel(str);

  return result;
}

}  // end namespace Teko

#endif