Thyra_ScaledAdjointLinearOpBase_def.hpp

00001 // @HEADER
00002 // ***********************************************************************
00003 // 
00004 //    Thyra: Interfaces and Support for Abstract Numerical Algorithms
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 THYRA_SCALED_ADJOINT_LINEAR_OP_BASE_HPP
00030 #define THYRA_SCALED_ADJOINT_LINEAR_OP_BASE_HPP
00031 
00032 #include "Thyra_ScaledAdjointLinearOpBase_decl.hpp"
00033 #include "Thyra_LinearOpBase.hpp"
00034 
00035 
00036 template<class Scalar>
00037 void Thyra::unwrap(
00038   const LinearOpBase<Scalar> &Op,
00039   Scalar *scalar,
00040   EOpTransp *transp,
00041   const LinearOpBase<Scalar>* *origOp
00042   )
00043 {
00044 #ifdef TEUCHOS_DEBUG
00045   TEST_FOR_EXCEPT( scalar==NULL );
00046   TEST_FOR_EXCEPT( transp==NULL );
00047   TEST_FOR_EXCEPT( origOp==NULL );
00048 #endif
00049   typedef Teuchos::ScalarTraits<Scalar>  ST;
00050   const ScaledAdjointLinearOpBase<Scalar>
00051     *saOp = dynamic_cast<const ScaledAdjointLinearOpBase<Scalar>*>(&Op);
00052   if(saOp) {
00053     *scalar = saOp->overallScalar();
00054     *transp = saOp->overallTransp();
00055     *origOp = &*saOp->getOrigOp();
00056   }
00057   else {
00058     *scalar = ST::one();
00059     *transp = NOTRANS;
00060     *origOp = &Op;
00061   }
00062 }
00063 
00064 
00065 template<class Scalar>
00066 void Thyra::unwrap(
00067   const RCP<const LinearOpBase<Scalar> > &Op,
00068   Scalar *scalar,
00069   EOpTransp *transp,
00070   RCP<const LinearOpBase<Scalar> > *origOp
00071   )
00072 {
00073 #ifdef TEUCHOS_DEBUG
00074   TEST_FOR_EXCEPT( scalar==NULL );
00075   TEST_FOR_EXCEPT( transp==NULL );
00076   TEST_FOR_EXCEPT( origOp==NULL );
00077 #endif
00078   typedef Teuchos::ScalarTraits<Scalar>  ST;
00079   RCP<const ScaledAdjointLinearOpBase<Scalar> >
00080     saOp = Teuchos::rcp_dynamic_cast<const ScaledAdjointLinearOpBase<Scalar> >(Op);
00081   if(saOp.get()) {
00082     *scalar = saOp->overallScalar();
00083     *transp = saOp->overallTransp();
00084     *origOp = saOp->getOrigOp();
00085   }
00086   else {
00087     *scalar = ST::one();
00088     *transp = NOTRANS;
00089     *origOp = Op;
00090   }
00091 }
00092 
00093 
00094 //
00095 // Explicit instant macro
00096 //
00097 
00098 #define THYRA_SCALED_ADJOINT_LINEAR_OP_BASE_INSTANT(SCALAR) \
00099   \
00100   template void unwrap( \
00101     const LinearOpBase<SCALAR > &Op, \
00102     SCALAR  *scalar, \
00103     EOpTransp *transp, \
00104     const LinearOpBase<SCALAR >* *origOp \
00105     ); \
00106    \
00107   template void unwrap( \
00108     const RCP<const LinearOpBase<SCALAR > > &Op, \
00109     SCALAR  *scalar, \
00110     EOpTransp *transp, \
00111     RCP<const LinearOpBase<SCALAR > > *origOp \
00112     );
00113 
00114 
00115 #endif  // THYRA_SCALED_ADJOINT_LINEAR_OP_BASE_HPP

Generated on Tue Jul 13 09:26:23 2010 for Thyra by  doxygen 1.4.7