Thyra Version of the Day
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 // Redistribution and use in source and binary forms, with or without
00011 // modification, are permitted provided that the following conditions are
00012 // met:
00013 //
00014 // 1. Redistributions of source code must retain the above copyright
00015 // notice, this list of conditions and the following disclaimer.
00016 //
00017 // 2. Redistributions in binary form must reproduce the above copyright
00018 // notice, this list of conditions and the following disclaimer in the
00019 // documentation and/or other materials provided with the distribution.
00020 //
00021 // 3. Neither the name of the Corporation nor the names of the
00022 // contributors may be used to endorse or promote products derived from
00023 // this software without specific prior written permission.
00024 //
00025 // THIS SOFTWARE IS PROVIDED BY SANDIA CORPORATION "AS IS" AND ANY
00026 // EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
00027 // IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR
00028 // PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL SANDIA CORPORATION OR THE
00029 // CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL,
00030 // EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO,
00031 // PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR
00032 // PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF
00033 // LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING
00034 // NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF THIS
00035 // SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
00036 //
00037 // Questions? Contact Roscoe A. Bartlett (bartlettra@ornl.gov) 
00038 // 
00039 // ***********************************************************************
00040 // @HEADER
00041 
00042 #ifndef THYRA_SCALED_ADJOINT_LINEAR_OP_BASE_HPP
00043 #define THYRA_SCALED_ADJOINT_LINEAR_OP_BASE_HPP
00044 
00045 #include "Thyra_ScaledAdjointLinearOpBase_decl.hpp"
00046 #include "Thyra_LinearOpBase.hpp"
00047 
00048 
00049 template<class Scalar>
00050 void Thyra::unwrap(
00051   const LinearOpBase<Scalar> &Op,
00052   Scalar *scalar,
00053   EOpTransp *transp,
00054   const LinearOpBase<Scalar>* *origOp
00055   )
00056 {
00057 #ifdef TEUCHOS_DEBUG
00058   TEUCHOS_TEST_FOR_EXCEPT( scalar==NULL );
00059   TEUCHOS_TEST_FOR_EXCEPT( transp==NULL );
00060   TEUCHOS_TEST_FOR_EXCEPT( origOp==NULL );
00061 #endif
00062   typedef Teuchos::ScalarTraits<Scalar>  ST;
00063   const ScaledAdjointLinearOpBase<Scalar>
00064     *saOp = dynamic_cast<const ScaledAdjointLinearOpBase<Scalar>*>(&Op);
00065   if(saOp) {
00066     *scalar = saOp->overallScalar();
00067     *transp = saOp->overallTransp();
00068     *origOp = &*saOp->getOrigOp();
00069   }
00070   else {
00071     *scalar = ST::one();
00072     *transp = NOTRANS;
00073     *origOp = &Op;
00074   }
00075 }
00076 
00077 
00078 template<class Scalar>
00079 void Thyra::unwrap(
00080   const RCP<const LinearOpBase<Scalar> > &Op,
00081   Scalar *scalar,
00082   EOpTransp *transp,
00083   RCP<const LinearOpBase<Scalar> > *origOp
00084   )
00085 {
00086 #ifdef TEUCHOS_DEBUG
00087   TEUCHOS_TEST_FOR_EXCEPT( scalar==NULL );
00088   TEUCHOS_TEST_FOR_EXCEPT( transp==NULL );
00089   TEUCHOS_TEST_FOR_EXCEPT( origOp==NULL );
00090 #endif
00091   typedef Teuchos::ScalarTraits<Scalar>  ST;
00092   RCP<const ScaledAdjointLinearOpBase<Scalar> >
00093     saOp = Teuchos::rcp_dynamic_cast<const ScaledAdjointLinearOpBase<Scalar> >(Op);
00094   if(saOp.get()) {
00095     *scalar = saOp->overallScalar();
00096     *transp = saOp->overallTransp();
00097     *origOp = saOp->getOrigOp();
00098   }
00099   else {
00100     *scalar = ST::one();
00101     *transp = NOTRANS;
00102     *origOp = Op;
00103   }
00104 }
00105 
00106 
00107 //
00108 // Explicit instant macro
00109 //
00110 
00111 #define THYRA_SCALED_ADJOINT_LINEAR_OP_BASE_INSTANT(SCALAR) \
00112   \
00113   template void unwrap( \
00114     const LinearOpBase<SCALAR > &Op, \
00115     SCALAR  *scalar, \
00116     EOpTransp *transp, \
00117     const LinearOpBase<SCALAR >* *origOp \
00118     ); \
00119    \
00120   template void unwrap( \
00121     const RCP<const LinearOpBase<SCALAR > > &Op, \
00122     SCALAR  *scalar, \
00123     EOpTransp *transp, \
00124     RCP<const LinearOpBase<SCALAR > > *origOp \
00125     );
00126 
00127 
00128 #endif  // THYRA_SCALED_ADJOINT_LINEAR_OP_BASE_HPP
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Defines