SparseDenseProduct.h
Go to the documentation of this file.
1 // This file is part of Eigen, a lightweight C++ template library
2 // for linear algebra.
3 //
4 // Copyright (C) 2008-2015 Gael Guennebaud <gael.guennebaud@inria.fr>
5 //
6 // This Source Code Form is subject to the terms of the Mozilla
7 // Public License v. 2.0. If a copy of the MPL was not distributed
8 // with this file, You can obtain one at http://mozilla.org/MPL/2.0/.
9 
10 #ifndef EIGEN_SPARSEDENSEPRODUCT_H
11 #define EIGEN_SPARSEDENSEPRODUCT_H
12 
13 namespace Eigen {
14 
15 namespace internal {
16 
17 template <> struct product_promote_storage_type<Sparse,Dense, OuterProduct> { typedef Sparse ret; };
18 template <> struct product_promote_storage_type<Dense,Sparse, OuterProduct> { typedef Sparse ret; };
19 
20 template<typename SparseLhsType, typename DenseRhsType, typename DenseResType,
21  typename AlphaType,
22  int LhsStorageOrder = ((SparseLhsType::Flags&RowMajorBit)==RowMajorBit) ? RowMajor : ColMajor,
23  bool ColPerCol = ((DenseRhsType::Flags&RowMajorBit)==0) || DenseRhsType::ColsAtCompileTime==1>
25 
26 template<typename SparseLhsType, typename DenseRhsType, typename DenseResType>
27 struct sparse_time_dense_product_impl<SparseLhsType,DenseRhsType,DenseResType, typename DenseResType::Scalar, RowMajor, true>
28 {
34  static void run(const SparseLhsType& lhs, const DenseRhsType& rhs, DenseResType& res, const typename Res::Scalar& alpha)
35  {
36  LhsEval lhsEval(lhs);
37 
38  Index n = lhs.outerSize();
39 #ifdef EIGEN_HAS_OPENMP
41  Index threads = Eigen::nbThreads();
42 #endif
43 
44  for(Index c=0; c<rhs.cols(); ++c)
45  {
46 #ifdef EIGEN_HAS_OPENMP
47  // This 20000 threshold has been found experimentally on 2D and 3D Poisson problems.
48  // It basically represents the minimal amount of work to be done to be worth it.
49  if(threads>1 && lhsEval.nonZerosEstimate() > 20000)
50  {
51  #pragma omp parallel for schedule(dynamic,(n+threads*4-1)/(threads*4)) num_threads(threads)
52  for(Index i=0; i<n; ++i)
53  processRow(lhsEval,rhs,res,alpha,i,c);
54  }
55  else
56 #endif
57  {
58  for(Index i=0; i<n; ++i)
59  processRow(lhsEval,rhs,res,alpha,i,c);
60  }
61  }
62  }
63 
64  static void processRow(const LhsEval& lhsEval, const DenseRhsType& rhs, DenseResType& res, const typename Res::Scalar& alpha, Index i, Index col)
65  {
66  typename Res::Scalar tmp(0);
67  for(LhsInnerIterator it(lhsEval,i); it ;++it)
68  tmp += it.value() * rhs.coeff(it.index(),col);
69  res.coeffRef(i,col) += alpha * tmp;
70  }
71 
72 };
73 
74 // FIXME: what is the purpose of the following specialization? Is it for the BlockedSparse format?
75 // -> let's disable it for now as it is conflicting with generic scalar*matrix and matrix*scalar operators
76 // template<typename T1, typename T2/*, int _Options, typename _StrideType*/>
77 // struct ScalarBinaryOpTraits<T1, Ref<T2/*, _Options, _StrideType*/> >
78 // {
79 // enum {
80 // Defined = 1
81 // };
82 // typedef typename CwiseUnaryOp<scalar_multiple2_op<T1, typename T2::Scalar>, T2>::PlainObject ReturnType;
83 // };
84 
85 template<typename SparseLhsType, typename DenseRhsType, typename DenseResType, typename AlphaType>
86 struct sparse_time_dense_product_impl<SparseLhsType,DenseRhsType,DenseResType, AlphaType, ColMajor, true>
87 {
92  typedef typename LhsEval::InnerIterator LhsInnerIterator;
93  static void run(const SparseLhsType& lhs, const DenseRhsType& rhs, DenseResType& res, const AlphaType& alpha)
94  {
95  LhsEval lhsEval(lhs);
96  for(Index c=0; c<rhs.cols(); ++c)
97  {
98  for(Index j=0; j<lhs.outerSize(); ++j)
99  {
100 // typename Res::Scalar rhs_j = alpha * rhs.coeff(j,c);
102  for(LhsInnerIterator it(lhsEval,j); it ;++it)
103  res.coeffRef(it.index(),c) += it.value() * rhs_j;
104  }
105  }
106  }
107 };
108 
109 template<typename SparseLhsType, typename DenseRhsType, typename DenseResType>
110 struct sparse_time_dense_product_impl<SparseLhsType,DenseRhsType,DenseResType, typename DenseResType::Scalar, RowMajor, false>
111 {
116  typedef typename LhsEval::InnerIterator LhsInnerIterator;
117  static void run(const SparseLhsType& lhs, const DenseRhsType& rhs, DenseResType& res, const typename Res::Scalar& alpha)
118  {
119  Index n = lhs.rows();
120  LhsEval lhsEval(lhs);
121 
122 #ifdef EIGEN_HAS_OPENMP
124  Index threads = Eigen::nbThreads();
125  // This 20000 threshold has been found experimentally on 2D and 3D Poisson problems.
126  // It basically represents the minimal amount of work to be done to be worth it.
127  if(threads>1 && lhsEval.nonZerosEstimate()*rhs.cols() > 20000)
128  {
129  #pragma omp parallel for schedule(dynamic,(n+threads*4-1)/(threads*4)) num_threads(threads)
130  for(Index i=0; i<n; ++i)
131  processRow(lhsEval,rhs,res,alpha,i);
132  }
133  else
134 #endif
135  {
136  for(Index i=0; i<n; ++i)
137  processRow(lhsEval, rhs, res, alpha, i);
138  }
139  }
140 
141  static void processRow(const LhsEval& lhsEval, const DenseRhsType& rhs, Res& res, const typename Res::Scalar& alpha, Index i)
142  {
143  typename Res::RowXpr res_i(res.row(i));
144  for(LhsInnerIterator it(lhsEval,i); it ;++it)
145  res_i += (alpha*it.value()) * rhs.row(it.index());
146  }
147 };
148 
149 template<typename SparseLhsType, typename DenseRhsType, typename DenseResType>
150 struct sparse_time_dense_product_impl<SparseLhsType,DenseRhsType,DenseResType, typename DenseResType::Scalar, ColMajor, false>
151 {
156  static void run(const SparseLhsType& lhs, const DenseRhsType& rhs, DenseResType& res, const typename Res::Scalar& alpha)
157  {
158  evaluator<Lhs> lhsEval(lhs);
159  for(Index j=0; j<lhs.outerSize(); ++j)
160  {
161  typename Rhs::ConstRowXpr rhs_j(rhs.row(j));
162  for(LhsInnerIterator it(lhsEval,j); it ;++it)
163  res.row(it.index()) += (alpha*it.value()) * rhs_j;
164  }
165  }
166 };
167 
168 template<typename SparseLhsType, typename DenseRhsType, typename DenseResType,typename AlphaType>
169 inline void sparse_time_dense_product(const SparseLhsType& lhs, const DenseRhsType& rhs, DenseResType& res, const AlphaType& alpha)
170 {
172 }
173 
174 } // end namespace internal
175 
176 namespace internal {
177 
178 template<typename Lhs, typename Rhs, int ProductType>
180  : generic_product_impl_base<Lhs,Rhs,generic_product_impl<Lhs,Rhs,SparseShape,DenseShape,ProductType> >
181 {
183 
184  template<typename Dest>
185  static void scaleAndAddTo(Dest& dst, const Lhs& lhs, const Rhs& rhs, const Scalar& alpha)
186  {
187  typedef typename nested_eval<Lhs,((Rhs::Flags&RowMajorBit)==0) ? 1 : Rhs::ColsAtCompileTime>::type LhsNested;
188  typedef typename nested_eval<Rhs,((Lhs::Flags&RowMajorBit)==0) ? 1 : Dynamic>::type RhsNested;
189  LhsNested lhsNested(lhs);
190  RhsNested rhsNested(rhs);
191  internal::sparse_time_dense_product(lhsNested, rhsNested, dst, alpha);
192  }
193 };
194 
195 template<typename Lhs, typename Rhs, int ProductType>
197  : generic_product_impl<Lhs, Rhs, SparseShape, DenseShape, ProductType>
198 {};
199 
200 template<typename Lhs, typename Rhs, int ProductType>
202  : generic_product_impl_base<Lhs,Rhs,generic_product_impl<Lhs,Rhs,DenseShape,SparseShape,ProductType> >
203 {
205 
206  template<typename Dst>
207  static void scaleAndAddTo(Dst& dst, const Lhs& lhs, const Rhs& rhs, const Scalar& alpha)
208  {
209  typedef typename nested_eval<Lhs,((Rhs::Flags&RowMajorBit)==0) ? Dynamic : 1>::type LhsNested;
210  typedef typename nested_eval<Rhs,((Lhs::Flags&RowMajorBit)==RowMajorBit) ? 1 : Lhs::RowsAtCompileTime>::type RhsNested;
211  LhsNested lhsNested(lhs);
212  RhsNested rhsNested(rhs);
213 
214  // transpose everything
215  Transpose<Dst> dstT(dst);
216  internal::sparse_time_dense_product(rhsNested.transpose(), lhsNested.transpose(), dstT, alpha);
217  }
218 };
219 
220 template<typename Lhs, typename Rhs, int ProductType>
222  : generic_product_impl<Lhs, Rhs, DenseShape, SparseShape, ProductType>
223 {};
224 
225 template<typename LhsT, typename RhsT, bool NeedToTranspose>
227 {
228 protected:
232 
233  // if the actual left-hand side is a dense vector,
234  // then build a sparse-view so that we can seamlessly iterate over it.
239 
243  typedef typename ProdXprType::Scalar Scalar;
244 
245 public:
246  enum {
247  Flags = NeedToTranspose ? RowMajorBit : 0,
249  };
250 
251  class InnerIterator : public LhsIterator
252  {
253  public:
255  : LhsIterator(xprEval.m_lhsXprImpl, 0),
256  m_outer(outer),
257  m_empty(false),
258  m_factor(get(xprEval.m_rhsXprImpl, outer, typename internal::traits<ActualRhs>::StorageKind() ))
259  {}
260 
262  EIGEN_STRONG_INLINE Index row() const { return NeedToTranspose ? m_outer : LhsIterator::index(); }
263  EIGEN_STRONG_INLINE Index col() const { return NeedToTranspose ? LhsIterator::index() : m_outer; }
264 
266  EIGEN_STRONG_INLINE operator bool() const { return LhsIterator::operator bool() && (!m_empty); }
267 
268  protected:
269  Scalar get(const RhsEval &rhs, Index outer, Dense = Dense()) const
270  {
271  return rhs.coeff(outer);
272  }
273 
275  {
276  typename RhsEval::InnerIterator it(rhs, outer);
277  if (it && it.index()==0 && it.value()!=Scalar(0))
278  return it.value();
279  m_empty = true;
280  return Scalar(0);
281  }
282 
284  bool m_empty;
286  };
287 
289  : m_lhs(lhs), m_lhsXprImpl(m_lhs), m_rhsXprImpl(rhs)
290  {
292  }
293 
294  // transpose case
296  : m_lhs(lhs), m_lhsXprImpl(m_lhs), m_rhsXprImpl(rhs)
297  {
299  }
300 
301 protected:
302  const LhsArg m_lhs;
305 };
306 
307 // sparse * dense outer product
308 template<typename Lhs, typename Rhs>
310  : sparse_dense_outer_product_evaluator<Lhs,Rhs, Lhs::IsRowMajor>
311 {
313 
315  typedef typename XprType::PlainObject PlainObject;
316 
317  explicit product_evaluator(const XprType& xpr)
318  : Base(xpr.lhs(), xpr.rhs())
319  {}
320 
321 };
322 
323 template<typename Lhs, typename Rhs>
325  : sparse_dense_outer_product_evaluator<Lhs,Rhs, Rhs::IsRowMajor>
326 {
328 
330  typedef typename XprType::PlainObject PlainObject;
331 
332  explicit product_evaluator(const XprType& xpr)
333  : Base(xpr.lhs(), xpr.rhs())
334  {}
335 
336 };
337 
338 } // end namespace internal
339 
340 } // end namespace Eigen
341 
342 #endif // EIGEN_SPARSEDENSEPRODUCT_H
Eigen::DefaultProduct
@ DefaultProduct
Definition: Constants.h:500
Eigen::internal::sparse_time_dense_product_impl< SparseLhsType, DenseRhsType, DenseResType, typename DenseResType::Scalar, ColMajor, false >::Res
internal::remove_all< DenseResType >::type Res
Definition: SparseDenseProduct.h:154
Eigen::internal::generic_product_impl< Lhs, Rhs, DenseShape, SparseShape, ProductType >::Scalar
Product< Lhs, Rhs >::Scalar Scalar
Definition: SparseDenseProduct.h:204
Eigen::internal::product_evaluator< Product< Lhs, Rhs, DefaultProduct >, OuterProduct, DenseShape, SparseShape >::PlainObject
XprType::PlainObject PlainObject
Definition: SparseDenseProduct.h:330
Eigen::internal::Lhs
@ Lhs
Definition: TensorContractionMapper.h:19
Eigen::internal::sparse_time_dense_product_impl< SparseLhsType, DenseRhsType, DenseResType, typename DenseResType::Scalar, ColMajor, false >::Lhs
internal::remove_all< SparseLhsType >::type Lhs
Definition: SparseDenseProduct.h:152
Eigen::HugeCost
const int HugeCost
Definition: Constants.h:44
Eigen::nbThreads
int nbThreads()
Definition: Parallelizer.h:63
Eigen
Namespace containing all symbols from the Eigen library.
Definition: jet.h:637
RowXpr
Block< Derived, 1, internal::traits< Derived >::ColsAtCompileTime, IsRowMajor > RowXpr
Definition: BlockMethods.h:17
col
m col(1)
Eigen::internal::product_promote_storage_type< Sparse, Dense, OuterProduct >::ret
Sparse ret
Definition: SparseDenseProduct.h:17
Eigen::internal::sparse_dense_outer_product_evaluator::InnerIterator::outer
EIGEN_STRONG_INLINE Index outer() const
Definition: SparseDenseProduct.h:261
alpha
RealScalar alpha
Definition: level1_cplx_impl.h:147
Eigen::DenseShape
Definition: Constants.h:528
Eigen::internal::SparseTriangularShape
Definition: SparseUtil.h:137
Eigen::Sparse
Definition: Constants.h:510
Eigen::internal::sparse_time_dense_product_impl< SparseLhsType, DenseRhsType, DenseResType, typename DenseResType::Scalar, RowMajor, false >::Res
internal::remove_all< DenseResType >::type Res
Definition: SparseDenseProduct.h:114
Eigen::internal::sparse_time_dense_product_impl< SparseLhsType, DenseRhsType, DenseResType, typename DenseResType::Scalar, ColMajor, false >::run
static void run(const SparseLhsType &lhs, const DenseRhsType &rhs, DenseResType &res, const typename Res::Scalar &alpha)
Definition: SparseDenseProduct.h:156
c
Scalar Scalar * c
Definition: benchVecAdd.cpp:17
Eigen::internal::product_evaluator< Product< Lhs, Rhs, DefaultProduct >, OuterProduct, DenseShape, SparseShape >::XprType
Product< Lhs, Rhs > XprType
Definition: SparseDenseProduct.h:329
Eigen::RowMajorBit
const unsigned int RowMajorBit
Definition: Constants.h:66
Eigen::internal::product_evaluator< Product< Lhs, Rhs, DefaultProduct >, OuterProduct, SparseShape, DenseShape >::PlainObject
XprType::PlainObject PlainObject
Definition: SparseDenseProduct.h:315
Eigen::internal::sparse_time_dense_product_impl< SparseLhsType, DenseRhsType, DenseResType, typename DenseResType::Scalar, ColMajor, false >::LhsInnerIterator
evaluator< Lhs >::InnerIterator LhsInnerIterator
Definition: SparseDenseProduct.h:155
Eigen::initParallel
void initParallel()
Definition: Parallelizer.h:53
Eigen::internal::sparse_dense_outer_product_evaluator::InnerIterator::get
Scalar get(const RhsEval &rhs, Index outer, Sparse=Sparse())
Definition: SparseDenseProduct.h:274
Eigen::internal::sparse_time_dense_product_impl< SparseLhsType, DenseRhsType, DenseResType, AlphaType, ColMajor, true >::run
static void run(const SparseLhsType &lhs, const DenseRhsType &rhs, DenseResType &res, const AlphaType &alpha)
Definition: SparseDenseProduct.h:93
type
Definition: pytypes.h:1525
Eigen::internal::generic_product_impl_base
Definition: ProductEvaluators.h:343
Eigen::internal::sparse_dense_outer_product_evaluator::Lhs1
conditional< NeedToTranspose, RhsT, LhsT >::type Lhs1
Definition: SparseDenseProduct.h:229
Eigen::internal::sparse_time_dense_product_impl< SparseLhsType, DenseRhsType, DenseResType, typename DenseResType::Scalar, RowMajor, true >::LhsInnerIterator
evaluator< Lhs >::InnerIterator LhsInnerIterator
Definition: SparseDenseProduct.h:32
Eigen::ScalarBinaryOpTraits
Determines whether the given binary operation of two numeric types is allowed and what the scalar ret...
Definition: XprHelper.h:801
Eigen::internal::nested_eval
Definition: XprHelper.h:457
Eigen::RowMajor
@ RowMajor
Definition: Constants.h:321
res
cout<< "Here is the matrix m:"<< endl<< m<< endl;Matrix< ptrdiff_t, 3, 1 > res
Definition: PartialRedux_count.cpp:3
Eigen::internal::generic_product_impl< Lhs, Rhs, SparseShape, DenseShape, ProductType >::Scalar
Product< Lhs, Rhs >::Scalar Scalar
Definition: SparseDenseProduct.h:182
Eigen::internal::sparse_dense_outer_product_evaluator::LhsIterator
evaluator< ActualLhs >::InnerIterator LhsIterator
Definition: SparseDenseProduct.h:242
Eigen::internal::product_evaluator< Product< Lhs, Rhs, DefaultProduct >, OuterProduct, SparseShape, DenseShape >::XprType
Product< Lhs, Rhs > XprType
Definition: SparseDenseProduct.h:314
Eigen::internal::sparse_time_dense_product_impl< SparseLhsType, DenseRhsType, DenseResType, typename DenseResType::Scalar, RowMajor, false >::LhsInnerIterator
LhsEval::InnerIterator LhsInnerIterator
Definition: SparseDenseProduct.h:116
Eigen::Transpose
Expression of the transpose of a matrix.
Definition: Transpose.h:52
Eigen::internal::product_evaluator< Product< Lhs, Rhs, DefaultProduct >, OuterProduct, SparseShape, DenseShape >::Base
sparse_dense_outer_product_evaluator< Lhs, Rhs, Lhs::IsRowMajor > Base
Definition: SparseDenseProduct.h:312
n
int n
Definition: BiCGSTAB_simple.cpp:1
Eigen::internal::generic_product_impl
Definition: ProductEvaluators.h:86
Eigen::internal::sparse_time_dense_product_impl< SparseLhsType, DenseRhsType, DenseResType, typename DenseResType::Scalar, RowMajor, false >::run
static void run(const SparseLhsType &lhs, const DenseRhsType &rhs, DenseResType &res, const typename Res::Scalar &alpha)
Definition: SparseDenseProduct.h:117
Eigen::internal::sparse_time_dense_product_impl< SparseLhsType, DenseRhsType, DenseResType, typename DenseResType::Scalar, RowMajor, false >::Rhs
internal::remove_all< DenseRhsType >::type Rhs
Definition: SparseDenseProduct.h:113
Eigen::internal::product_evaluator< Product< Lhs, Rhs, DefaultProduct >, OuterProduct, DenseShape, SparseShape >::product_evaluator
product_evaluator(const XprType &xpr)
Definition: SparseDenseProduct.h:332
Eigen::internal::generic_product_impl< Lhs, Rhs, DenseShape, SparseShape, ProductType >::scaleAndAddTo
static void scaleAndAddTo(Dst &dst, const Lhs &lhs, const Rhs &rhs, const Scalar &alpha)
Definition: SparseDenseProduct.h:207
Eigen::internal::true_type
Definition: Meta.h:96
Eigen::internal::sparse_dense_outer_product_evaluator::Flags
@ Flags
Definition: SparseDenseProduct.h:247
Eigen::internal::product_evaluator< Product< Lhs, Rhs, DefaultProduct >, OuterProduct, DenseShape, SparseShape >::Base
sparse_dense_outer_product_evaluator< Lhs, Rhs, Rhs::IsRowMajor > Base
Definition: SparseDenseProduct.h:327
j
std::ptrdiff_t j
Definition: tut_arithmetic_redux_minmax.cpp:2
Eigen::internal::sparse_time_dense_product
void sparse_time_dense_product(const SparseLhsType &lhs, const DenseRhsType &rhs, DenseResType &res, const AlphaType &alpha)
Definition: SparseDenseProduct.h:169
Eigen::internal::sparse_dense_outer_product_evaluator::m_lhsXprImpl
evaluator< ActualLhs > m_lhsXprImpl
Definition: SparseDenseProduct.h:303
Eigen::internal::sparse_dense_outer_product_evaluator
Definition: SparseDenseProduct.h:226
Eigen::internal::sparse_time_dense_product_impl
Definition: SparseDenseProduct.h:24
Eigen::internal::sparse_time_dense_product_impl< SparseLhsType, DenseRhsType, DenseResType, AlphaType, ColMajor, true >::LhsInnerIterator
LhsEval::InnerIterator LhsInnerIterator
Definition: SparseDenseProduct.h:92
Eigen::internal::get
Definition: CXX11Meta.h:115
Eigen::Dynamic
const int Dynamic
Definition: Constants.h:22
Eigen::SparseShape
Definition: Constants.h:537
Eigen::internal::sparse_dense_outer_product_evaluator::sparse_dense_outer_product_evaluator
sparse_dense_outer_product_evaluator(const ActualRhs &rhs, const Lhs1 &lhs)
Definition: SparseDenseProduct.h:295
EIGEN_STRONG_INLINE
#define EIGEN_STRONG_INLINE
Definition: Macros.h:917
Eigen::internal::sparse_time_dense_product_impl< SparseLhsType, DenseRhsType, DenseResType, typename DenseResType::Scalar, RowMajor, true >::run
static void run(const SparseLhsType &lhs, const DenseRhsType &rhs, DenseResType &res, const typename Res::Scalar &alpha)
Definition: SparseDenseProduct.h:34
Eigen::internal::sparse_dense_outer_product_evaluator::m_lhs
const LhsArg m_lhs
Definition: SparseDenseProduct.h:302
Eigen::internal::sparse_dense_outer_product_evaluator::sparse_dense_outer_product_evaluator
sparse_dense_outer_product_evaluator(const Lhs1 &lhs, const ActualRhs &rhs)
Definition: SparseDenseProduct.h:288
Eigen::internal::product_evaluator
Definition: ForwardDeclarations.h:164
Eigen::internal::sparse_time_dense_product_impl< SparseLhsType, DenseRhsType, DenseResType, typename DenseResType::Scalar, RowMajor, false >::Lhs
internal::remove_all< SparseLhsType >::type Lhs
Definition: SparseDenseProduct.h:112
Eigen::Product
Expression of the product of two arbitrary matrices or vectors.
Definition: Product.h:71
gtsam.examples.DogLegOptimizerExample.run
def run(args)
Definition: DogLegOptimizerExample.py:21
Eigen::internal::generic_product_impl< Lhs, Rhs, SparseShape, DenseShape, ProductType >::scaleAndAddTo
static void scaleAndAddTo(Dest &dst, const Lhs &lhs, const Rhs &rhs, const Scalar &alpha)
Definition: SparseDenseProduct.h:185
Eigen::internal::sparse_dense_outer_product_evaluator::InnerIterator
Definition: SparseDenseProduct.h:251
Eigen::Triplet< double >
Eigen::internal::product_promote_storage_type
Definition: XprHelper.h:578
Eigen::internal::evaluator
Definition: CoreEvaluators.h:90
Eigen::internal::sparse_dense_outer_product_evaluator::InnerIterator::m_empty
bool m_empty
Definition: SparseDenseProduct.h:284
Eigen::internal::sparse_time_dense_product_impl< SparseLhsType, DenseRhsType, DenseResType, typename DenseResType::Scalar, RowMajor, false >::LhsEval
evaluator< Lhs > LhsEval
Definition: SparseDenseProduct.h:115
Eigen::internal::sparse_time_dense_product_impl< SparseLhsType, DenseRhsType, DenseResType, AlphaType, ColMajor, true >::Lhs
internal::remove_all< SparseLhsType >::type Lhs
Definition: SparseDenseProduct.h:88
Eigen::internal::sparse_dense_outer_product_evaluator::LhsArg
conditional< is_same< typename internal::traits< Lhs1 >::StorageKind, Sparse >::value, Lhs1 const &, SparseView< Lhs1 > >::type LhsArg
Definition: SparseDenseProduct.h:238
Eigen::internal::sparse_time_dense_product_impl< SparseLhsType, DenseRhsType, DenseResType, typename DenseResType::Scalar, RowMajor, true >::Lhs
internal::remove_all< SparseLhsType >::type Lhs
Definition: SparseDenseProduct.h:29
Eigen::internal::sparse_time_dense_product_impl< SparseLhsType, DenseRhsType, DenseResType, AlphaType, ColMajor, true >::Res
internal::remove_all< DenseResType >::type Res
Definition: SparseDenseProduct.h:90
Eigen::internal::traits
Definition: ForwardDeclarations.h:17
gtsam::Sparse
Eigen::SparseMatrix< double > Sparse
Definition: AcceleratedPowerMethod.h:26
Eigen::internal::sparse_dense_outer_product_evaluator::InnerIterator::InnerIterator
InnerIterator(const sparse_dense_outer_product_evaluator &xprEval, Index outer)
Definition: SparseDenseProduct.h:254
Eigen::internal::sparse_dense_outer_product_evaluator::InnerIterator::value
EIGEN_STRONG_INLINE Scalar value() const
Definition: SparseDenseProduct.h:265
Eigen::internal::Rhs
@ Rhs
Definition: TensorContractionMapper.h:18
Eigen::internal::sparse_time_dense_product_impl< SparseLhsType, DenseRhsType, DenseResType, typename DenseResType::Scalar, RowMajor, true >::Rhs
internal::remove_all< DenseRhsType >::type Rhs
Definition: SparseDenseProduct.h:30
Eigen::internal::sparse_dense_outer_product_evaluator::ActualRhs
conditional< NeedToTranspose, LhsT, RhsT >::type ActualRhs
Definition: SparseDenseProduct.h:230
Eigen::internal::conditional
Definition: Meta.h:109
Eigen::internal::sparse_time_dense_product_impl< SparseLhsType, DenseRhsType, DenseResType, AlphaType, ColMajor, true >::LhsEval
evaluator< Lhs > LhsEval
Definition: SparseDenseProduct.h:91
Eigen::internal::sparse_time_dense_product_impl< SparseLhsType, DenseRhsType, DenseResType, typename DenseResType::Scalar, ColMajor, false >::Rhs
internal::remove_all< DenseRhsType >::type Rhs
Definition: SparseDenseProduct.h:153
Eigen::SparseView
Expression of a dense or sparse matrix with zero or too small values removed.
Definition: ForwardDeclarations.h:124
Eigen::internal::sparse_dense_outer_product_evaluator::RhsEval
evaluator< ActualRhs > RhsEval
Definition: SparseDenseProduct.h:241
Eigen::internal::sparse_time_dense_product_impl< SparseLhsType, DenseRhsType, DenseResType, typename DenseResType::Scalar, RowMajor, true >::processRow
static void processRow(const LhsEval &lhsEval, const DenseRhsType &rhs, DenseResType &res, const typename Res::Scalar &alpha, Index i, Index col)
Definition: SparseDenseProduct.h:64
Eigen::internal::sparse_time_dense_product_impl< SparseLhsType, DenseRhsType, DenseResType, AlphaType, ColMajor, true >::Rhs
internal::remove_all< DenseRhsType >::type Rhs
Definition: SparseDenseProduct.h:89
Eigen::internal::sparse_dense_outer_product_evaluator::InnerIterator::m_factor
Scalar m_factor
Definition: SparseDenseProduct.h:285
Eigen::internal::sparse_dense_outer_product_evaluator::Scalar
ProdXprType::Scalar Scalar
Definition: SparseDenseProduct.h:243
Eigen::internal::product_evaluator< Product< Lhs, Rhs, DefaultProduct >, OuterProduct, SparseShape, DenseShape >::product_evaluator
product_evaluator(const XprType &xpr)
Definition: SparseDenseProduct.h:317
Eigen::internal::sparse_time_dense_product_impl< SparseLhsType, DenseRhsType, DenseResType, typename DenseResType::Scalar, RowMajor, true >::Res
internal::remove_all< DenseResType >::type Res
Definition: SparseDenseProduct.h:31
Eigen::internal::sparse_dense_outer_product_evaluator::ActualLhs
conditional< is_same< typename internal::traits< Lhs1 >::StorageKind, Sparse >::value, Lhs1, SparseView< Lhs1 > >::type ActualLhs
Definition: SparseDenseProduct.h:236
internal
Definition: BandTriangularSolver.h:13
Eigen::internal::sparse_time_dense_product_impl< SparseLhsType, DenseRhsType, DenseResType, typename DenseResType::Scalar, RowMajor, false >::processRow
static void processRow(const LhsEval &lhsEval, const DenseRhsType &rhs, Res &res, const typename Res::Scalar &alpha, Index i)
Definition: SparseDenseProduct.h:141
Eigen::internal::product_promote_storage_type< Dense, Sparse, OuterProduct >::ret
Sparse ret
Definition: SparseDenseProduct.h:18
Eigen::internal::sparse_dense_outer_product_evaluator::InnerIterator::col
EIGEN_STRONG_INLINE Index col() const
Definition: SparseDenseProduct.h:263
Eigen::internal::sparse_dense_outer_product_evaluator::InnerIterator::get
Scalar get(const RhsEval &rhs, Index outer, Dense=Dense()) const
Definition: SparseDenseProduct.h:269
Eigen::ColMajor
@ ColMajor
Definition: Constants.h:319
EIGEN_INTERNAL_CHECK_COST_VALUE
#define EIGEN_INTERNAL_CHECK_COST_VALUE(C)
Definition: StaticAssert.h:218
Eigen::OuterProduct
@ OuterProduct
Definition: Constants.h:500
Eigen::internal::sparse_dense_outer_product_evaluator::InnerIterator::row
EIGEN_STRONG_INLINE Index row() const
Definition: SparseDenseProduct.h:262
Eigen::internal::sparse_dense_outer_product_evaluator::LhsEval
evaluator< ActualLhs > LhsEval
Definition: SparseDenseProduct.h:240
ConstRowXpr
const typedef Block< const Derived, 1, internal::traits< Derived >::ColsAtCompileTime, IsRowMajor > ConstRowXpr
Definition: BlockMethods.h:18
Eigen::internal::sparse_dense_outer_product_evaluator::m_rhsXprImpl
evaluator< ActualRhs > m_rhsXprImpl
Definition: SparseDenseProduct.h:304
Eigen::internal::sparse_dense_outer_product_evaluator::CoeffReadCost
@ CoeffReadCost
Definition: SparseDenseProduct.h:248
test_callbacks.value
value
Definition: test_callbacks.py:160
i
int i
Definition: BiCGSTAB_step_by_step.cpp:9
Eigen::internal::sparse_time_dense_product_impl< SparseLhsType, DenseRhsType, DenseResType, typename DenseResType::Scalar, RowMajor, true >::LhsEval
evaluator< Lhs > LhsEval
Definition: SparseDenseProduct.h:33
Eigen::internal::sparse_dense_outer_product_evaluator::ProdXprType
Product< LhsT, RhsT, DefaultProduct > ProdXprType
Definition: SparseDenseProduct.h:231
Eigen::internal::sparse_dense_outer_product_evaluator::InnerIterator::m_outer
Index m_outer
Definition: SparseDenseProduct.h:283
Scalar
SCALAR Scalar
Definition: bench_gemm.cpp:46
Eigen::Index
EIGEN_DEFAULT_DENSE_INDEX_TYPE Index
The Index type as used for the API.
Definition: Meta.h:74
Eigen::Dense
Definition: Constants.h:507


gtsam
Author(s):
autogenerated on Sat Nov 16 2024 04:04:30