Thyra_get_Epetra_Operator.hpp

// @HEADER
// *****************************************************************************
//    Thyra: Interfaces and Support for Abstract Numerical Algorithms
//
// Copyright 2004 NTESS and the Thyra contributors.
// SPDX-License-Identifier: BSD-3-Clause
// *****************************************************************************
// @HEADER

#ifndef THYRA_GET_EPETRA_OPERATOR_HPP
#define THYRA_GET_EPETRA_OPERATOR_HPP

#include "Thyra_EpetraTypes.hpp"
#include "Teuchos_TestForException.hpp"
#include "Teuchos_TypeNameTraits.hpp"
#include <stdexcept> // std::invalid_argument

namespace Thyra {

template<class Scalar>
Teuchos::RCP<Epetra_Operator>
get_Epetra_Operator (LinearOpBase<Scalar>& /* op */)
{
  TEUCHOS_TEST_FOR_EXCEPTION
    (true, std::invalid_argument, "Thyra::get_Epetra_Operator: This function "
     "only works if Scalar=double, because for Epetra objects, the only Scalar"
     " type is double.  Instead, Scalar = " <<
     Teuchos::TypeNameTraits<Scalar>::name () << ".");
}

template<>
Teuchos::RCP<Epetra_Operator>
get_Epetra_Operator( LinearOpBase<double> &op );

template<class Scalar>
Teuchos::RCP<const Epetra_Operator>
get_Epetra_Operator( const LinearOpBase<Scalar> & /* op */ )
{
  TEUCHOS_TEST_FOR_EXCEPTION
    (true, std::invalid_argument, "Thyra::get_Epetra_Operator: This function "
     "only works if Scalar=double, because for Epetra objects, the only Scalar"
     " type is double.  Instead, Scalar = " <<
     Teuchos::TypeNameTraits<Scalar>::name () << ".");
}

template<>
Teuchos::RCP<const Epetra_Operator>
get_Epetra_Operator( const LinearOpBase<double> &op );

} // namespace Thyra


#endif // THYRA_GET_EPETRA_OPERATOR_HPP

#if defined(Thyra_SHOW_DEPRECATED_WARNINGS)
#ifdef __GNUC__
#warning "The ThyraEpetraAdapters package is deprecated"
#endif
#endif