00001
00002
00003
00004
00005
00006
00007
00008
00009
00010
00011
00012
00013
00014
00015
00016
00017
00018
00019
00020
00021
00022
00023
00024
00025
00026
00027
00028
00029 #ifndef THYRA_SINGLE_SCALAR_LINEAR_OP_WITH_SOLVE_BASE_DECL_HPP
00030 #define THYRA_SINGLE_SCALAR_LINEAR_OP_WITH_SOLVE_BASE_DECL_HPP
00031
00032 #include "Thyra_LinearOpWithSolveBaseDecl.hpp"
00033 #include "Thyra_SingleScalarLinearOpBaseDecl.hpp"
00034
00035 namespace Thyra {
00036
00038 template <class Scalar>
00039 class SingleScalarLinearOpWithSolveBase
00040 : virtual public LinearOpWithSolveBase<Scalar>
00041 , virtual protected SingleScalarLinearOpBase<Scalar>
00042 {
00043 public:
00044
00046 SingleScalarLinearOpBase<Scalar>::apply;
00047
00050
00052 bool solveSupportsConj(EConj conj) const;
00054 bool solveTransposeSupportsConj(EConj conj) const;
00056 bool solveSupportsSolveMeasureType(EConj conj, const SolveMeasureType& solveMeasureType) const;
00058 bool solveTransposeSupportsSolveMeasureType(EConj conj, const SolveMeasureType& solveMeasureType) const;
00060 void solve(
00061 const EConj conj
00062 ,const MultiVectorBase<Scalar> &B
00063 ,MultiVectorBase<Scalar> *X
00064 ,const int numBlocks
00065 ,const BlockSolveCriteria<Scalar> blockSolveCriteria[]
00066 ,SolveStatus<Scalar> blockSolveStatus[]
00067 ) const;
00069 void solveTranspose(
00070 const EConj conj
00071 ,const MultiVectorBase<Scalar> &B
00072 ,MultiVectorBase<Scalar> *X
00073 ,const int numBlocks
00074 ,const BlockSolveCriteria<Scalar> blockSolveCriteria[]
00075 ,SolveStatus<Scalar> blockSolveStatus[]
00076 ) const;
00077
00079
00080 protected:
00081
00084
00086 virtual bool solveSupportsTrans(ETransp M_trans) const = 0;
00088 virtual bool solveSupportsSolveMeasureType(ETransp M_trans, const SolveMeasureType& solveMeasureType) const = 0;
00090 virtual void solve(
00091 const ETransp M_trans
00092 ,const MultiVectorBase<Scalar> &B
00093 ,MultiVectorBase<Scalar> *X
00094 ,const int numBlocks
00095 ,const BlockSolveCriteria<Scalar> blockSolveCriteria[]
00096 ,SolveStatus<Scalar> blockSolveStatus[]
00097 ) const = 0;
00098
00100
00101 };
00102
00103 }
00104
00105 #endif // THYRA_SINGLE_SCALAR_LINEAR_OP_WITH_SOLVE_BASE_DECL_HPP