Teko_LSCPreconditionerFactory.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_LSCPreconditionerFactory_hpp__
#define __Teko_LSCPreconditionerFactory_hpp__
 
#include "Teko_BlockPreconditionerFactory.hpp"
#include "Teko_LSCStrategy.hpp"
 
namespace Teko {
namespace NS {  // Navier-Stokes specialization

class LSCPrecondState : public BlockPreconditionerState {
 public:
  LSCPrecondState() {}

  LinearOp invMass_;

  ModifiableLinearOp BQBt_;

  ModifiableLinearOp BHBt_;

  LinearOp BQBtmC_;
  InverseLinearOp invBQBtmC_;

  LinearOp BHBtmC_;
  InverseLinearOp invBHBtmC_;

  LinearOp aiD_;

  double gamma_;

  double alpha_;
};
 
class LSCPreconditionerFactory : public BlockPreconditionerFactory {
 public:


  LSCPreconditionerFactory(const LinearOp& invF, const LinearOp& invBQBtmC, const LinearOp& invD,
                           const LinearOp& invMass);

  LSCPreconditionerFactory(const LinearOp& invF, const LinearOp& invBQBtmC,
                           const LinearOp& invMass);

  LSCPreconditionerFactory(const Teuchos::RCP<LSCStrategy>& strategy);

  LSCPreconditionerFactory();

  virtual LinearOp buildPreconditionerOperator(BlockedLinearOp& blo,
                                               BlockPreconditionerState& state) const;

  virtual RCP<PreconditionerState> buildPreconditionerState() const {
    return rcp(new LSCPrecondState());
  }

  virtual Teuchos::RCP<Teuchos::ParameterList> getRequestedParameters() const;

  virtual bool updateRequestedParameters(const Teuchos::ParameterList& pl);
 
 protected:
  using Teko::BlockPreconditionerFactory::buildPreconditionerOperator;
 
  // Gimmie object
  Teuchos::RCP<LSCStrategy> invOpsStrategy_;
  bool isSymmetric_;

  virtual void initializeFromParameterList(const Teuchos::ParameterList& pl);
 
 public:
  static RCP<LSCStrategy> buildStrategy(const std::string& name,
                                        const Teuchos::ParameterList& settings,
                                        const RCP<const InverseLibrary>& invLib,
                                        const RCP<RequestHandler>& rh);

  static void addStrategy(const std::string& name, const RCP<Cloneable>& clone);
 
 private:
  static CloneFactory<LSCStrategy> strategyBuilder_;

  static void initializeStrategyBuilder();
};
 
}  // end namespace NS
}  // end namespace Teko
 
#endif