00001 // @HEADER 00002 // *********************************************************************** 00003 // 00004 // Meros: Segregated Preconditioning Package 00005 // Copyright (2004) Sandia Corporation 00006 // 00007 // Under terms of Contract DE-AC04-94AL85000, there is a non-exclusive 00008 // license for use of this work by or on behalf of the U.S. Government. 00009 // 00010 // This library is free software; you can redistribute it and/or modify 00011 // it under the terms of the GNU Lesser General Public License as 00012 // published by the Free Software Foundation; either version 2.1 of the 00013 // License, or (at your option) any later version. 00014 // 00015 // This library is distributed in the hope that it will be useful, but 00016 // WITHOUT ANY WARRANTY; without even the implied warranty of 00017 // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU 00018 // Lesser General Public License for more details. 00019 // 00020 // You should have received a copy of the GNU Lesser General Public 00021 // License along with this library; if not, write to the Free Software 00022 // Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 00023 // USA 00024 // Questions? Contact Michael A. Heroux (maherou@sandia.gov) 00025 // 00026 // *********************************************************************** 00027 // @HEADER 00028 00029 #ifndef MEROS_INVERSE_OPERATOR_H 00030 #define MEROS_INVERSE_OPERATOR_H 00031 00032 #include "Thyra_VectorImpl.hpp" // need for LinOpDecl 00033 #include "Thyra_VectorSpaceImpl.hpp" // need for LinOpDecl 00034 #include "Thyra_LinearOperatorImpl.hpp" 00035 #include "Thyra_DefaultInverseLinearOp.hpp" 00036 #include "Meros_LinearSolveStrategy.hpp" 00037 00038 00039 namespace Meros 00040 { 00041 using namespace Teuchos; 00042 using namespace Thyra; 00043 00047 template <class Scalar> 00048 class InverseOperator 00049 : public Handleable<LinearOpBase<Scalar> >, 00050 public DefaultInverseLinearOp<Scalar> 00051 { 00052 public: 00054 InverseOperator(const LinearOperator<Scalar>& op, 00055 const LinearSolveStrategy<Scalar>& solver) 00056 : DefaultInverseLinearOp<Scalar>(solver.getLOWS(op)) 00057 { 00058 } 00059 00060 /* */ 00062 InverseOperator(const ConstLinearOperator<Scalar>& op, 00063 const LinearSolveStrategy<Scalar>& solver) 00064 : DefaultInverseLinearOp<Scalar>(solver.getLOWS(op)) 00065 { 00066 } 00067 00068 /* */ 00069 TEUCHOS_GET_RCP(LinearOpBase<Scalar> ); 00070 }; 00071 00072 00073 } // namespace Meros 00074 00075 #endif