Teko_StaticLSCStrategy.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_StaticLSCStrategy_hpp__
#define __Teko_StaticLSCStrategy_hpp__
 
#include "Teko_LSCStrategy.hpp"
 
namespace Teko {
namespace NS {
 
class LSCPrecondState;  // forward declaration
 
// constant, not very flexible strategy for driving LSCPreconditioenrFactory
class StaticLSCStrategy : public LSCStrategy {
 public:
  // Staiblized constructor
  StaticLSCStrategy(const LinearOp& invF, const LinearOp& invBQBtmC, const LinearOp& invD,
                    const LinearOp& invMass);
 
  // Stable constructor
  StaticLSCStrategy(const LinearOp& invF, const LinearOp& invBQBtmC, const LinearOp& invMass);

  virtual void buildState(BlockedLinearOp& /* A */, BlockPreconditionerState& /* state */) const {}

  virtual LinearOp getInvF(const BlockedLinearOp& /* A */,
                           BlockPreconditionerState& /* state */) const {
    return invF_;
  }

  virtual LinearOp getInvBQBt(const BlockedLinearOp& /* A */,
                              BlockPreconditionerState& /* state */) const {
    return invBQBtmC_;
  }

  virtual LinearOp getInvBHBt(const BlockedLinearOp& /* A */,
                              BlockPreconditionerState& /* state */) const {
    return invBQBtmC_;
  }

  virtual LinearOp getOuterStabilization(const BlockedLinearOp& /* A */,
                                         BlockPreconditionerState& /* state */) const {
    return invD_;
  }
 
  virtual LinearOp getInnerStabilization(const BlockedLinearOp& /* A */,
                                         BlockPreconditionerState& /* state */) const {
    return Teuchos::null;
  }

  virtual LinearOp getInvMass(const BlockedLinearOp& /* A */,
                              BlockPreconditionerState& /* state */) const {
    return invMass_;
  }

  virtual LinearOp getHScaling(const BlockedLinearOp& /* A */,
                               BlockPreconditionerState& /* state */) const {
    return invMass_;
  }

  virtual bool useFullLDU() const { return false; }

  virtual void setSymmetric(bool /* isSymmetric */) {}
 
 protected:
  // protected memebers
  LinearOp invF_;
  LinearOp invBQBtmC_;
  LinearOp invD_;
  LinearOp invMass_;
};
 
}  // end namespace NS
}  // end namespace Teko
 
#endif