Thyra Package Browser (Single Doxygen Collection) Version of the Day
Thyra_EpetraOperatorViewExtractorStd.cpp
Go to the documentation of this file.
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 #include "Thyra_EpetraOperatorViewExtractorStd.hpp"
00043 #include "Thyra_EpetraLinearOpBase.hpp"
00044 #include "Thyra_ScaledAdjointLinearOpBase.hpp"
00045 
00046 class Epetra_Operator;
00047 
00048 
00049 namespace Thyra {
00050 
00051 
00052 // Overridden from EpetraOperatorViewExtractorBase
00053 
00054 
00055 bool EpetraOperatorViewExtractorStd::isCompatible( const LinearOpBase<double> &fwdOp ) const
00056 {
00057   double wrappedScalar = 0.0;
00058   EOpTransp wrappedTransp = NOTRANS;
00059   const LinearOpBase<double> *wrappedFwdOp = NULL;
00060   ::Thyra::unwrap(fwdOp, &wrappedScalar, &wrappedTransp, &wrappedFwdOp);
00061   const EpetraLinearOpBase *eFwdOp = NULL;
00062   if( !(eFwdOp = dynamic_cast<const EpetraLinearOpBase*>(wrappedFwdOp)) )
00063     return false;
00064   return true;
00065 }
00066 
00067 
00068 void EpetraOperatorViewExtractorStd::getNonconstEpetraOpView(
00069   const RCP<LinearOpBase<double> > &fwdOp,
00070   const Ptr<RCP<Epetra_Operator> > &epetraOp,
00071   const Ptr<EOpTransp> &epetraOpTransp,
00072   const Ptr<EApplyEpetraOpAs> &epetraOpApplyAs,
00073   const Ptr<EAdjointEpetraOp> &epetraOpAdjointSupport,
00074     const Ptr<double> &epetraOpScalar
00075   ) const
00076 {
00077   TEUCHOS_TEST_FOR_EXCEPT(true);
00078   // ToDo: Implement once this is needed by just copying what is below and
00079   // removing the 'const' in the right places!
00080 }
00081 
00082 
00083 void EpetraOperatorViewExtractorStd::getEpetraOpView(
00084   const RCP<const LinearOpBase<double> > &fwdOp,
00085   const Ptr<RCP<const Epetra_Operator> > &epetraOp,
00086   const Ptr<EOpTransp> &epetraOpTransp,
00087   const Ptr<EApplyEpetraOpAs> &epetraOpApplyAs,
00088   const Ptr<EAdjointEpetraOp> &epetraOpAdjointSupport,
00089     const Ptr<double> &epetraOpScalar
00090   ) const
00091 {
00092   using Teuchos::outArg;
00093   double wrappedFwdOpScalar = 0.0;
00094   EOpTransp wrappedFwdOpTransp = NOTRANS;
00095   Teuchos::RCP<const LinearOpBase<double> > wrappedFwdOp; 
00096   unwrap(fwdOp,&wrappedFwdOpScalar, &wrappedFwdOpTransp, &wrappedFwdOp);
00097   Teuchos::RCP<const EpetraLinearOpBase> epetraFwdOp =
00098     Teuchos::rcp_dynamic_cast<const EpetraLinearOpBase>(wrappedFwdOp,true);
00099   EOpTransp epetra_epetraOpTransp;
00100   epetraFwdOp->getEpetraOpView(epetraOp, outArg(epetra_epetraOpTransp),
00101     epetraOpApplyAs, epetraOpAdjointSupport);
00102   *epetraOpTransp = trans_trans(real_trans(epetra_epetraOpTransp), wrappedFwdOpTransp);
00103   *epetraOpScalar = wrappedFwdOpScalar;
00104 }
00105 
00106 
00107 // ToDo: Refactor unwrap(...) to not use raw pointers!
00108 
00109 
00110 } // namespace Thyra
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Defines