diff options
| author | Stanislaw Halik <sthalik@misaki.pl> | 2019-03-03 21:09:10 +0100 |
|---|---|---|
| committer | Stanislaw Halik <sthalik@misaki.pl> | 2019-03-03 21:10:13 +0100 |
| commit | f0238cfb6997c4acfc2bd200de7295f3fa36968f (patch) | |
| tree | b215183760e4f615b9c1dabc1f116383b72a1b55 /eigen/bench/BenchSparseUtil.h | |
| parent | 543edd372a5193d04b3de9f23c176ab439e51b31 (diff) | |
don't index Eigen
Diffstat (limited to 'eigen/bench/BenchSparseUtil.h')
| -rw-r--r-- | eigen/bench/BenchSparseUtil.h | 149 |
1 files changed, 0 insertions, 149 deletions
diff --git a/eigen/bench/BenchSparseUtil.h b/eigen/bench/BenchSparseUtil.h deleted file mode 100644 index 13981f6..0000000 --- a/eigen/bench/BenchSparseUtil.h +++ /dev/null @@ -1,149 +0,0 @@ - -#include <Eigen/Sparse> -#include <bench/BenchTimer.h> -#include <set> - -using namespace std; -using namespace Eigen; -using namespace Eigen; - -#ifndef SIZE -#define SIZE 1024 -#endif - -#ifndef DENSITY -#define DENSITY 0.01 -#endif - -#ifndef SCALAR -#define SCALAR double -#endif - -typedef SCALAR Scalar; -typedef Matrix<Scalar,Dynamic,Dynamic> DenseMatrix; -typedef Matrix<Scalar,Dynamic,1> DenseVector; -typedef SparseMatrix<Scalar> EigenSparseMatrix; - -void fillMatrix(float density, int rows, int cols, EigenSparseMatrix& dst) -{ - dst.reserve(double(rows)*cols*density); - for(int j = 0; j < cols; j++) - { - for(int i = 0; i < rows; i++) - { - Scalar v = (internal::random<float>(0,1) < density) ? internal::random<Scalar>() : 0; - if (v!=0) - dst.insert(i,j) = v; - } - } - dst.finalize(); -} - -void fillMatrix2(int nnzPerCol, int rows, int cols, EigenSparseMatrix& dst) -{ -// std::cout << "alloc " << nnzPerCol*cols << "\n"; - dst.reserve(nnzPerCol*cols); - for(int j = 0; j < cols; j++) - { - std::set<int> aux; - for(int i = 0; i < nnzPerCol; i++) - { - int k = internal::random<int>(0,rows-1); - while (aux.find(k)!=aux.end()) - k = internal::random<int>(0,rows-1); - aux.insert(k); - - dst.insert(k,j) = internal::random<Scalar>(); - } - } - dst.finalize(); -} - -void eiToDense(const EigenSparseMatrix& src, DenseMatrix& dst) -{ - dst.setZero(); - for (int j=0; j<src.cols(); ++j) - for (EigenSparseMatrix::InnerIterator it(src.derived(), j); it; ++it) - dst(it.index(),j) = it.value(); -} - -#ifndef NOGMM -#include "gmm/gmm.h" -typedef gmm::csc_matrix<Scalar> GmmSparse; -typedef gmm::col_matrix< gmm::wsvector<Scalar> > GmmDynSparse; -void eiToGmm(const EigenSparseMatrix& src, GmmSparse& dst) -{ - GmmDynSparse tmp(src.rows(), src.cols()); - for (int j=0; j<src.cols(); ++j) - for (EigenSparseMatrix::InnerIterator it(src.derived(), j); it; ++it) - tmp(it.index(),j) = it.value(); - gmm::copy(tmp, dst); -} -#endif - -#ifndef NOMTL -#include <boost/numeric/mtl/mtl.hpp> -typedef mtl::compressed2D<Scalar, mtl::matrix::parameters<mtl::tag::col_major> > MtlSparse; -typedef mtl::compressed2D<Scalar, mtl::matrix::parameters<mtl::tag::row_major> > MtlSparseRowMajor; -void eiToMtl(const EigenSparseMatrix& src, MtlSparse& dst) -{ - mtl::matrix::inserter<MtlSparse> ins(dst); - for (int j=0; j<src.cols(); ++j) - for (EigenSparseMatrix::InnerIterator it(src.derived(), j); it; ++it) - ins[it.index()][j] = it.value(); -} -#endif - -#ifdef CSPARSE -extern "C" { -#include "cs.h" -} -void eiToCSparse(const EigenSparseMatrix& src, cs* &dst) -{ - cs* aux = cs_spalloc (0, 0, 1, 1, 1); - for (int j=0; j<src.cols(); ++j) - for (EigenSparseMatrix::InnerIterator it(src.derived(), j); it; ++it) - if (!cs_entry(aux, it.index(), j, it.value())) - { - std::cout << "cs_entry error\n"; - exit(2); - } - dst = cs_compress(aux); -// cs_spfree(aux); -} -#endif // CSPARSE - -#ifndef NOUBLAS -#include <boost/numeric/ublas/vector.hpp> -#include <boost/numeric/ublas/matrix.hpp> -#include <boost/numeric/ublas/io.hpp> -#include <boost/numeric/ublas/triangular.hpp> -#include <boost/numeric/ublas/vector_sparse.hpp> -#include <boost/numeric/ublas/matrix_sparse.hpp> -#include <boost/numeric/ublas/vector_of_vector.hpp> -#include <boost/numeric/ublas/operation.hpp> - -typedef boost::numeric::ublas::compressed_matrix<Scalar,boost::numeric::ublas::column_major> UBlasSparse; - -void eiToUblas(const EigenSparseMatrix& src, UBlasSparse& dst) -{ - dst.resize(src.rows(), src.cols(), false); - for (int j=0; j<src.cols(); ++j) - for (EigenSparseMatrix::InnerIterator it(src.derived(), j); it; ++it) - dst(it.index(),j) = it.value(); -} - -template <typename EigenType, typename UblasType> -void eiToUblasVec(const EigenType& src, UblasType& dst) -{ - dst.resize(src.size()); - for (int j=0; j<src.size(); ++j) - dst[j] = src.coeff(j); -} -#endif - -#ifdef OSKI -extern "C" { -#include <oski/oski.h> -} -#endif |