Teko_BlockPreconditionerFactory.cpp

// @HEADER
// *****************************************************************************
//      Teko: A package for block and physics based preconditioning
//
// Copyright 2010 NTESS and the Teko contributors.
// SPDX-License-Identifier: BSD-3-Clause
// *****************************************************************************
// @HEADER
 
#include "Teko_BlockPreconditionerFactory.hpp"
 
#include "Teko_Preconditioner.hpp"
#include "Teko_InverseLibrary.hpp"
 
#include "Thyra_DefaultPreconditioner.hpp"
 
using namespace Thyra;
 
namespace Teko {

 
LinearOp BlockPreconditionerFactory::buildPreconditionerOperator(LinearOp &lo,
                                                                 PreconditionerState &state) const {
  // get the blocked linear operator
  RCP<LinearOpBase<double> > loA = Teuchos::rcp_const_cast<Thyra::LinearOpBase<double> >(lo);
  BlockedLinearOp A = Teuchos::rcp_dynamic_cast<Thyra::PhysicallyBlockedLinearOpBase<double> >(loA);
 
  state.setInitialized(false);
 
  return buildPreconditionerOperator(A, dynamic_cast<BlockPreconditionerState &>(state));
}

bool BlockPreconditionerFactory::isCompatible(
    const Thyra::LinearOpSourceBase<double> &fwdOpSrc) const {
  RCP<const Thyra::PhysicallyBlockedLinearOpBase<double> > A =
      Teuchos::rcp_dynamic_cast<const Thyra::PhysicallyBlockedLinearOpBase<double> >(
          fwdOpSrc.getOp());
  return A != Teuchos::null;
}
 
}  // end namespace Teko