SparseSymMatProd.h
Go to the documentation of this file.
1 // Copyright (C) 2016-2019 Yixuan Qiu <yixuan.qiu@cos.name>
2 //
3 // This Source Code Form is subject to the terms of the Mozilla
4 // Public License v. 2.0. If a copy of the MPL was not distributed
5 // with this file, You can obtain one at https://mozilla.org/MPL/2.0/.
6 
7 #ifndef SPARSE_SYM_MAT_PROD_H
8 #define SPARSE_SYM_MAT_PROD_H
9 
10 #include <Eigen/Core>
11 #include <Eigen/SparseCore>
12 
13 namespace Spectra {
14 
22 template <typename Scalar, int Uplo = Eigen::Lower, int Flags = 0, typename StorageIndex = int>
24 {
25 private:
32 
33  ConstGenericSparseMatrix m_mat;
34 
35 public:
43  SparseSymMatProd(ConstGenericSparseMatrix& mat) :
44  m_mat(mat)
45  {}
46 
50  Index rows() const { return m_mat.rows(); }
54  Index cols() const { return m_mat.cols(); }
55 
62  // y_out = A * x_in
63  void perform_op(const Scalar* x_in, Scalar* y_out) const
64  {
65  MapConstVec x(x_in, m_mat.cols());
66  MapVec y(y_out, m_mat.rows());
67  y.noalias() = m_mat.template selfadjointView<Uplo>() * x;
68  }
69 };
70 
71 } // namespace Spectra
72 
73 #endif // SPARSE_SYM_MAT_PROD_H
SCALAR Scalar
Definition: bench_gemm.cpp:46
Scalar * y
ConstGenericSparseMatrix m_mat
A versatible sparse matrix representation.
Definition: SparseMatrix.h:96
A matrix or vector expression mapping an existing array of data.
Definition: Map.h:94
EIGEN_DEFAULT_DENSE_INDEX_TYPE Index
The Index type as used for the API.
Definition: Meta.h:74
Eigen::Map< const Vector > MapConstVec
const Eigen::Ref< const SparseMatrix > ConstGenericSparseMatrix
Eigen::SparseMatrix< Scalar, Flags, StorageIndex > SparseMatrix
void perform_op(const Scalar *x_in, Scalar *y_out) const
set noclip points set clip one set noclip two set bar set border lt lw set xdata set ydata set zdata set x2data set y2data set boxwidth set dummy x
Eigen::Map< Vector > MapVec
SparseSymMatProd(ConstGenericSparseMatrix &mat)
Eigen::Matrix< Scalar, Eigen::Dynamic, 1 > Vector


gtsam
Author(s):
autogenerated on Tue Jul 4 2023 02:36:11