15template<
typename Decomposition,
typename RhsType,
typename StorageKind>
class SolveImpl;
32template<
typename Decomposition,
typename RhsType,
typename StorageKind>
struct solve_traits;
34template<
typename Decomposition,
typename RhsType>
35struct solve_traits<Decomposition,RhsType,Dense>
37 typedef typename make_proper_matrix_type<
typename RhsType::Scalar,
38 Decomposition::ColsAtCompileTime,
39 RhsType::ColsAtCompileTime,
40 RhsType::PlainObject::Options,
41 Decomposition::MaxColsAtCompileTime,
42 RhsType::MaxColsAtCompileTime>::type PlainObject;
45template<
typename Decomposition,
typename RhsType>
46struct traits<Solve<Decomposition, RhsType> >
47 : traits<typename solve_traits<Decomposition,RhsType,typename internal::traits<RhsType>::StorageKind>::PlainObject>
49 typedef typename solve_traits<Decomposition,RhsType,typename internal::traits<RhsType>::StorageKind>::PlainObject PlainObject;
50 typedef typename promote_index_type<typename Decomposition::StorageIndex, typename RhsType::StorageIndex>::type StorageIndex;
51 typedef traits<PlainObject> BaseTraits;
61template<
typename Decomposition,
typename RhsType>
62class Solve :
public SolveImpl<Decomposition,RhsType,typename internal::traits<RhsType>::StorageKind>
65 typedef typename internal::traits<Solve>::PlainObject PlainObject;
66 typedef typename internal::traits<Solve>::StorageIndex StorageIndex;
68 Solve(
const Decomposition &dec,
const RhsType &rhs)
69 : m_dec(dec), m_rhs(rhs)
72 EIGEN_DEVICE_FUNC EIGEN_CONSTEXPR
Index rows()
const EIGEN_NOEXCEPT {
return m_dec.cols(); }
73 EIGEN_DEVICE_FUNC EIGEN_CONSTEXPR
Index cols()
const EIGEN_NOEXCEPT {
return m_rhs.cols(); }
75 EIGEN_DEVICE_FUNC
const Decomposition& dec()
const {
return m_dec; }
76 EIGEN_DEVICE_FUNC
const RhsType& rhs()
const {
return m_rhs; }
79 const Decomposition &m_dec;
85template<
typename Decomposition,
typename RhsType>
86class SolveImpl<Decomposition,RhsType,
Dense>
87 :
public MatrixBase<Solve<Decomposition,RhsType> >
94 EIGEN_DENSE_PUBLIC_INTERFACE(Derived)
99 Scalar coeff(
Index i)
const;
103template<
typename Decomposition,
typename RhsType,
typename StorageKind>
104class SolveImpl :
public internal::generic_xpr_base<Solve<Decomposition,RhsType>, MatrixXpr, StorageKind>::type
107 typedef typename internal::generic_xpr_base<Solve<Decomposition,RhsType>, MatrixXpr, StorageKind>::type Base;
113template<
typename Decomposition,
typename RhsType>
114struct evaluator<Solve<Decomposition,RhsType> >
115 :
public evaluator<typename Solve<Decomposition,RhsType>::PlainObject>
117 typedef Solve<Decomposition,RhsType> SolveType;
118 typedef typename SolveType::PlainObject PlainObject;
119 typedef evaluator<PlainObject> Base;
123 EIGEN_DEVICE_FUNC
explicit evaluator(
const SolveType& solve)
124 : m_result(solve.rows(), solve.cols())
126 ::new (
static_cast<Base*
>(
this)) Base(m_result);
127 solve.dec()._solve_impl(solve.rhs(), m_result);
131 PlainObject m_result;
136template<typename DstXprType, typename DecType, typename RhsType, typename Scalar>
137struct Assignment<DstXprType, Solve<DecType,RhsType>, internal::assign_op<Scalar,Scalar>, Dense2Dense>
139 typedef Solve<DecType,RhsType> SrcXprType;
140 static void run(DstXprType &dst,
const SrcXprType &src,
const internal::assign_op<Scalar,Scalar> &)
142 Index dstRows = src.rows();
143 Index dstCols = src.cols();
144 if((dst.rows()!=dstRows) || (dst.cols()!=dstCols))
145 dst.resize(dstRows, dstCols);
147 src.dec()._solve_impl(src.rhs(), dst);
152template<
typename DstXprType,
typename DecType,
typename RhsType,
typename Scalar>
153struct Assignment<DstXprType, Solve<Transpose<const DecType>,RhsType>, internal::assign_op<Scalar,Scalar>, Dense2Dense>
155 typedef Solve<Transpose<const DecType>,RhsType> SrcXprType;
156 static void run(DstXprType &dst,
const SrcXprType &src,
const internal::assign_op<Scalar,Scalar> &)
158 Index dstRows = src.rows();
159 Index dstCols = src.cols();
160 if((dst.rows()!=dstRows) || (dst.cols()!=dstCols))
161 dst.resize(dstRows, dstCols);
163 src.dec().nestedExpression().template _solve_impl_transposed<false>(src.rhs(), dst);
168template<
typename DstXprType,
typename DecType,
typename RhsType,
typename Scalar>
169struct Assignment<DstXprType, Solve<CwiseUnaryOp<internal::scalar_conjugate_op<typename DecType::Scalar>, const Transpose<const DecType> >,RhsType>,
170 internal::assign_op<Scalar,Scalar>, Dense2Dense>
172 typedef Solve<CwiseUnaryOp<internal::scalar_conjugate_op<typename DecType::Scalar>,
const Transpose<const DecType> >,RhsType> SrcXprType;
173 static void run(DstXprType &dst,
const SrcXprType &src,
const internal::assign_op<Scalar,Scalar> &)
175 Index dstRows = src.rows();
176 Index dstCols = src.cols();
177 if((dst.rows()!=dstRows) || (dst.cols()!=dstCols))
178 dst.resize(dstRows, dstCols);
180 src.dec().nestedExpression().nestedExpression().template _solve_impl_transposed<true>(src.rhs(), dst);
Base class for all dense matrices, vectors, and expressions.
Definition: MatrixBase.h:50
Pseudo expression representing a solving operation.
Definition: Solve.h:63
const unsigned int EvalBeforeNestingBit
Definition: Constants.h:70
const unsigned int RowMajorBit
Definition: Constants.h:66
Namespace containing all symbols from the Eigen library.
Definition: Core:141
const int HugeCost
Definition: Constants.h:44
EIGEN_DEFAULT_DENSE_INDEX_TYPE Index
The Index type as used for the API.
Definition: Meta.h:74
Definition: Constants.h:507