Stack Overflow
Stack Overflow

Questions tagged [suitesparse]

SuiteSparse is a suite of sparse matrix algorithms, including UMFPACK, CHOLMOD, SPQR, KLU, BTF, ordering methods (AMD, CAMD, COLAMD, and CCOLAMD), CSparse, CXSparse, UFget, spqr_rank, Factorize, SSMULT, SFMULT, and many other packages.

SuiteSparse is a suite of sparse matrix algorithms, including:

  • UMFPACK: multifrontal LU factorization. Appears as LU and x=A\b in MATLAB.
  • CHOLMOD: supernodal Cholesky. Appears as CHOL and x=A\b in MATLAB. Used in Google Ceres. Now with CUDA acceleration, in collaboration with NVIDIA.
  • SPQR: multifrontal QR. Appears as QR and x=A\b in MATLAB. CUDA acceleration just released (October 10, 2014, SuiteSparse 4.4.0), and submitted to ACM Trans. Math. Software
  • KLU and BTF: sparse LU factorization, well-suited for circuit simulation. Appears in Xyce by Sandia, and many commercial circuit simulation packages.
  • Ordering methods (AMD, CAMD, COLAMD, and CCOLAMD). AMD and COLAMD appear in MATLAB.
  • CSparse and CXSparse: a concise sparse Cholesky factorization package for the book, "direct methods for sparse linear systems" from Tim Davis published by SIAM.
  • UFget: a MATLAB interface for the UF Sparse Matrix Collection
  • spqr_rank: a MATLAB package for reliable sparse rank detection, null set bases, pseudoinverse solutions, and basic solutions.
  • Factorize: an object-oriented solver for MATLAB (a reusable backslash).
  • SSMULT and SFMULT: sparse matrix multiplication. Appears as C=A*B in MATLAB.
  • ... and many other packages.
50 questions
0
votes
0 answers

Using scikit-sparse's cholmod in cython wrapper?

I'm trying to write a Cython wrapper to interface C code from Python. The C library makes use of Suitesparse's CHOLMOD, so I figured installing scikit-sparse (which uses a cholmod.pyx that contains everything I need) would be a neat way. However, I…
lz129
  • 21
  • 4
0
votes
1 answer

How to find (Q, R ) from SuiteSparseQR_factorization object?

In C++ interface of SuiteSparse, I can use SuiteSparseQR_factorization *QR; QR = SuiteSparseQR_factorize(A) ; to calculate QR decomposition of matrix A so that I can reuse QR for further calculation. But I wonder can I get the real Q,R…
ask_q2015
  • 1
0
votes
1 answer

CHOLMOD sparse dense multiplication issue

I am testing out cholmod to solve a sparse matrix system. As b of the Ax=b equation I want to use something that will lead to a known result, so I choose to declare as xe a vector of ones and b=A*xe. I tried some different approaches to do the…
David
  • 369
  • 3
  • 14
0
votes
1 answer

Is it possible to use CHOLMOD with matrices already in device memory?

I would like to use CHOLMOD's GPU acceleration, and have found several simple examples on how to use the library for Cholesky decomposition. However all of the examples provide the matrices to CHOLMOD in host memory, and allow it to copy them to the…
nitronoid
  • 1,459
  • 11
  • 26
0
votes
1 answer

error LNK2001: unresolved external symbol Visaul C++ SuiteSparse lib

I used Visual C++ 2015. I received the following errors when I tried to build my g2o sample code. Error LNK2001 unresolved external symbol __imp_fflush Error LNK2001 unresolved external symbol __imp_fflush Error LNK2001 unresolved external…
Shengren Huang
  • 1
  • 1
0
votes
0 answers

Configure ceres-solver with CMake and can't find SuiteSparse config header

I'm trying to use CMake to configure Ceres to build on my Windows machine. I have Suitesparse v4.4.4 built on my machine and have pointed to all the appropriate directories and libraries. So I believe I have SuiteSparse properly included in the…
gboy
  • 633
  • 1
  • 6
  • 13
0
votes
0 answers

Is there a way to get/set the factorization of large sparse matrices within Eigen solvers?

I have a program where the system matrix within the sparse solver (using the Eigen library) is often the same, and could be precomputed once and loaded as a file. Ideally I would want to serialize the whole solver, but there seems to be no way to do…
sordid
  • 177
  • 6
0
votes
1 answer

Example of compress column format for rank-deficient matrices

It is the first time I deal with column-compress storage (CCS) format to store matrices. After googling a bit, if I am right, in a matrix having n nonzero elements the CCS is as follows: -we define a vector A_v of dimensions n x 1 storing the n…
venom
  • 173
  • 1
  • 6
0
votes
0 answers

Using CMake to link to a TPL that was compiled with NVCC

How do I properly link my project to a TPL (CHOLMOD) that employs CUDA? I want to link my CMake project to CHOLMOD (libcholmod.so). CHOLMOD was compiled separately with -GPU_BLAS, which is simply a flag that enables GPU functionality. The CHOLMOD…
trickleboast
  • 49
  • 7
0
votes
1 answer

Indexing into CHOLMOD dense vector array

I have a cholmod_dense data structure: cholmod_dense* ex = cholmod_l_solve(CHOLMOD_A, L, B, &com); and I want to extract the values and copy them to another variable. This means I need to index into the double array and copy values over. for (int…
trickleboast
  • 49
  • 7
0
votes
0 answers

Creating a dense matrix in C++ that can be solved using CHOLMOD

I'm working on a project that is trying to optimize the time it takes to solve a dense matrix. Right now, matrices are being processed within C++ and as it stands there are no plans to change that. Looking at the CHOLMOD user manual and examples…
Alex
  • 2,145
  • 6
  • 36
  • 72
0
votes
3 answers

Brew + suite-sparse421

I have installed brew on OSX and I have installed successfully several things. However when I try to install suite-sparse421 (hard constraint on the version): brew tap homebrew/versions brew install suite-sparse421 I get the following error: Error:…
dim_tz
  • 1,451
  • 5
  • 20
  • 37
0
votes
1 answer

Failed to use cs_qrsol from CXSparse to solve x=A\b in C++ when A matrix is large

I am trying to solve a linear equation system x = A\b, using CXSparse library by Tim Davis (http://faculty.cse.tamu.edu/davis/suitesparse.html). I develop my C++ program (with OpenCV) using MS Visual Studio 2012 on Windows 7 x64. I have an A matrix,…
Sonia
  • 195
  • 2
  • 11
0
votes
0 answers

Eigen SuiteSparse (Eigen::SPQR) mismatch num rows

i am using Eigen::SPQR to solve a sparse matrix system Here is the minimal example which would give run time error # include # include # include # include /// test eigen SPQR …
P. Preet
  • 15
  • 6
0
votes
1 answer

Matlab choose linear direct sparse solver

I need to solve a sparse linear system in Matalab using the QR solver of SuiteSparse. How can I choose it? I know that, depending on the matrix structure, the QR solver is used when calling the "\" operator. But I haven't find any way to directly…
Marco Agnese
  • 349
  • 4
  • 15
1 2
3
4