Thyra_LinearOpWithSolveBase.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_LINEAR_OP_WITH_SOLVE_BASE_HPP
00030 #define THYRA_LINEAR_OP_WITH_SOLVE_BASE_HPP
00031 
00032 #include "Thyra_LinearOpWithSolveBaseDecl.hpp"
00033 #include "Thyra_LinearOpBase.hpp"
00034 
00035 namespace Thyra {
00036 
00037 template <class RangeScalar, class DomainScalar>
00038 bool LinearOpWithSolveBase<RangeScalar,DomainScalar>::solveSupportsConj(EConj conj) const
00039 {
00040   return true;
00041 }
00042 
00043 template <class RangeScalar, class DomainScalar>
00044 bool LinearOpWithSolveBase<RangeScalar,DomainScalar>::solveTransposeSupportsConj(EConj conj) const
00045 {
00046   return false;
00047 }
00048 
00049 template <class RangeScalar, class DomainScalar>
00050 bool LinearOpWithSolveBase<RangeScalar,DomainScalar>::solveSupportsSolveMeasureType(EConj conj, const SolveMeasureType& solveMeasureType) const
00051 {
00052   return solveMeasureType.useDefault();
00053 }
00054 
00055 template <class RangeScalar, class DomainScalar>
00056 bool LinearOpWithSolveBase<RangeScalar,DomainScalar>::solveTransposeSupportsSolveMeasureType(EConj conj, const SolveMeasureType& solveMeasureType) const
00057 {
00058   return solveMeasureType.useDefault();
00059 }
00060 
00061 template <class RangeScalar, class DomainScalar>
00062 void LinearOpWithSolveBase<RangeScalar,DomainScalar>::solveTranspose(
00063   const EConj                                   conj
00064   ,const MultiVectorBase<DomainScalar>          &B
00065   ,MultiVectorBase<RangeScalar>                 *X
00066   ,const int                                    numBlocks
00067   ,const BlockSolveCriteria<PromotedScalar>     blockSolveCriteria[]
00068   ,SolveStatus<PromotedScalar>                  blockSolveStatus[]
00069   ) const
00070 {
00071   TEST_FOR_EXCEPTION(
00072     true,std::logic_error
00073     ,"LinearOpWithSolveBase<"<<Teuchos::ScalarTraits<RangeScalar>::name()<<","<<Teuchos::ScalarTraits<DomainScalar>::name()<<">::solveTranspose(...): "
00074     "Error, the concrete subclass described as { " << this->description() << " } "
00075     " with this->solveTransposeSupportsConj("<<toString(conj)<<")="<<this->solveTransposeSupportsConj(conj)
00076     << " did not override this function and does not support transposes."
00077     );
00078 }
00079 
00080 } // namespace Thyra
00081 
00082 #endif // THYRA_LINEAR_OP_WITH_SOLVE_BASE_HPP

Generated on Wed May 12 21:26:41 2010 for Thyra ANA Operator Solve Interfaces by  doxygen 1.4.7