43 #ifndef THYRA_DEFAULT_SERIAL_DENSE_LINEAR_OP_WITH_SOLVE_DECL_HPP 44 #define THYRA_DEFAULT_SERIAL_DENSE_LINEAR_OP_WITH_SOLVE_DECL_HPP 47 #include "Thyra_LinearOpWithSolveBase.hpp" 48 #include "RTOpPack_LapackWrappers.hpp" 55 inline RTOpPack::ETransp convertToRTOpPackETransp(
const EOpTransp transp )
58 TEUCHOS_TEST_FOR_EXCEPT(transp ==
CONJ);
68 TEUCHOS_TEST_FOR_EXCEPT(
true);
90 template<
class Scalar>
106 RCP<const LinearOpBase<Scalar> >
getFwdOp()
const;
114 RCP<const VectorSpaceBase<Scalar> >
range()
const;
116 RCP<const VectorSpaceBase<Scalar> >
domain()
const;
161 RCP<const MultiVectorBase<Scalar> > M_;
168 static void factorize(
171 const Ptr<Array<int> > &ipiv
174 static void backsolve(
176 const ArrayView<const int> ipiv,
193 template<
class Scalar>
194 RCP<DefaultSerialDenseLinearOpWithSolve<Scalar> >
205 template<
class Scalar>
206 RCP<DefaultSerialDenseLinearOpWithSolve<Scalar> >
209 RCP<DefaultSerialDenseLinearOpWithSolve<Scalar> >
211 M_lows->initialize(M);
219 #endif // THYRA_DEFAULT_SERIAL_DENSE_LINEAR_OP_WITH_SOLVE_DECL_HPP RCP< DefaultSerialDenseLinearOpWithSolve< Scalar > > defaultSerialDenseLinearOpWithSolve()
Nonmember constructor.
Base class for all linear operators that can support a high-level solve operation.
EOpTransp
Enumeration for determining how a linear operator is applied. `*.
RCP< const LinearOpBase< Scalar > > getFwdOp() const
Simple concreate subclass of LinearOpWithSolveBase for serial dense matrices implemented using LAPACK...
Use the non-transposed operator.
DefaultSerialDenseLinearOpWithSolve()
Use the transposed operator with complex-conjugate clements (same as TRANS for real scalar types)...
Use the non-transposed operator with complex-conjugate elements (same as NOTRANS for real scalar type...
bool opSupportedImpl(EOpTransp M_trans) const
void initialize(const RCP< const MultiVectorBase< Scalar > > &M)
Use the transposed operator.
bool solveSupportsSolveMeasureTypeImpl(EOpTransp M_trans, const SolveMeasureType &solveMeasureType) const
bool solveSupportsImpl(EOpTransp M_trans) const
Interface for a collection of column vectors called a multi-vector.
RCP< const VectorSpaceBase< Scalar > > range() const
Simple struct for the return status from a solve.
void applyImpl(const EOpTransp M_trans, const MultiVectorBase< Scalar > &X, const Ptr< MultiVectorBase< Scalar > > &Y, const Scalar alpha, const Scalar beta) const
SolveStatus< Scalar > solveImpl(const EOpTransp transp, const MultiVectorBase< Scalar > &B, const Ptr< MultiVectorBase< Scalar > > &X, const Ptr< const SolveCriteria< Scalar > > solveCriteria) const
RCP< DefaultSerialDenseLinearOpWithSolve< Scalar > > defaultSerialDenseLinearOpWithSolve(const RCP< const MultiVectorBase< Scalar > > &M)
Nonmember constructor.
RCP< const VectorSpaceBase< Scalar > > domain() const
Simple struct that defines the requested solution criteria for a solve.