All Questions
Tagged with blas sparse-matrix
24 questions
0
votes
1
answer
179
views
Fortran with Sparse BLAS not flushing memory
I have a subroutine that builds sparse matrices, and I need to call it several times. However, it seems that if I call this subroutine a lot of times (and/or if the sparse matrices are very large), ...
- 364
0
votes
1
answer
172
views
What is wrong with my sparse matrix-multiple vectors (SpMM) product function for CSR?
I have the following code for the sparse matrix-vector (SpMV) product in C assuming a CSR storage format:
void dcsrmv(SparseMatrixCSR *A, double *x, double *y) {
for (int i=0; i<A->m; i++) {
...
2
votes
0
answers
112
views
Sparse Blas on Mac OS Fortran
I want to use Sparse Blas operations on Mac OS in Fortran
Here is the programm
PROGRAM F95_EX
IMPLICIT NONE
INTEGER NMAX, NNZ
PARAMETER (NMAX=4, NNZ=6)
INTEGER i,n,a,istat
INTEGER,DIMENSION(:),...
- 139
-1
votes
2
answers
790
views
CUDA CSR Matrix-Matrix product transpose by itself
I have a very large, very sparse least-squares design matrix (A), which I would like to multiply by itself, as follows: N = A^T * A, where A & N are stored in CSR format. Obviously, A has more ...
- 97
1
vote
1
answer
711
views
number of operations for sparse*dense matrix multiplication
How many floating point operations does it take to multiply a CSR sparse x dense matrix or a dense x CSR sparse matrix using optimized sparse routines (like cuSparse or Eigen or Matlab).
In the limit ...
- 1,002
1
vote
0
answers
86
views
sparse-dense multiplication in matlab
Consider A*B in Matlab. A and B are matrices.
IF A and B are both dense it will call regular BLAS i.e. dGEMM.
What happens under the hood (which sparse library functions does Matlab call) if A and ...
- 1,002
4
votes
1
answer
506
views
Multicore Solving of Sparse Matrices in R
For an analysis I need to solve a large sparse system in the form Ax=b. A has the dimension of 3Mx3M and is positive definite.
For now, I use the 'Matrix' package which provides support for sparse ...
- 41
0
votes
1
answer
144
views
sparse BLAS solver used in a shared library doesn't work (returns '-1')
I have tried to program a solver-function in C (with Ubuntu in Code::Blocks). I want to call that solver from a Python script, so I have build it as a shared library and called it with 'ctypes' in ...
- 3
1
vote
2
answers
635
views
Sparse BLAS on OSX
I just discovered that SparseBLAS is included in the Accelerate Framework of OSX in OSX.
% ls -l /System/Library/Frameworks/Accelerate.framework/Versions/A/Frameworks/vecLib.framework/Versions/A
...
- 2,499
3
votes
2
answers
2k
views
Are sparse matrices typically stored in column major order or row major order?
A little background: I'm interested in doing some research on sparse matrix*vector multiplication.
I've been looking through this database of sparse matrices:
The University of Florida Sparse Matrix ...
- 3,667
0
votes
0
answers
273
views
Will intel mkl sparse BLAS coogemv internally change coo to csr storage?
There is sparse version of BLAS in intel MKL, for example doing multiplication of matrix and vector, we could use
mkl_?coogemv
Computes matrix-vector product of a sparse general matrix stored in
...
- 2,637
5
votes
3
answers
6k
views
Is sparse BLAS not included in BLAS?
I have a working LAPACK implementation and that, as far as I read, contains BLAS.
I want to use SPARSE BLAS and as far as I understand this website, SPARSE BLAS is part of BLAS.
But when I tried to ...
- 449
0
votes
1
answer
198
views
Sparse BLAS matrix-row vector product overwriting loop index
I am using the NIST Sparse BLAS v0.5 matrix-matrix multiplication routine, downloaded from http://math.nist.gov/~KRemington/fspblas/, to multiply a matrix one row at a time by a column vector. After ...
- 111
2
votes
0
answers
451
views
cuSPARSE csrmm with dense matrix in row-major format
I want to use the cuSPARSE csrmm function to multiply two matrices. The A matrix is sparse and the B matrix is dense. The dense matrix is in row-major format. Is there some nice way (trick) to ...
- 859
4
votes
1
answer
858
views
Does Scipy Sparse use the (sparse) BLAS library?
Numpy can use one of a number of BLAS libraries (eg. ATLAS, MKL, OpenBLAS etc.).
Does the scipy.sparse matrix module support the sparse BLAS library?
- 4,607