Sparse BLAS Routines

Sparse BLAS Level 1 Routines

?axpyi adds a scalar multiple of compressed sparse vector to a full-storage vector (see also ?axpy in Level 1 BLAS).
Fortran 77:
call saxpyi(nz, a, x, indx, y) call daxpyi(nz, a, x, indx, y) call caxpyi(nz, a, x, indx, y) call zaxpyi(nz, a, x, indx, y) Fortran 95:
call axpyi(x, indx, y [,a])

?doti computes the dot product of a compressed sparse real vector by a full-storage real vector (see also ?dot in Level 1 BLAS).
Fortran 77:
res = sdoti(nz, x, indx, y) res = ddoti(nz, x, indx, y) Fortran 95:
res = doti(x, indx, y)

?dotci computes the conjugated dot product of a compressed sparse complex vector with a full-storage complex vector (see also ?dotc in Level 1 BLAS).
Fortran 77:
res = cdotci(nz, x, indx, y) res = zdotci(nz, x, indx, y) Fortran 95:
res = dotci(x, indx, y)

?dotui computes the dot product of a compressed sparse complex vector by a full-storage complex vector, without conjugation (see also ?dotu in Level 1 BLAS).
Fortran 77:
res = cdotui(nz, x, indx, y) res = zdotui(nz, x, indx, y) Fortran 95:
res = dotui(x, indx, y)

?gthr gathers a full-storage sparse vector's elements into compressed form.
Fortran 77:
call sgthr(nz, y, x, indx) call dgthr(nz, y, x, indx) call cgthr(nz, y, x, indx) call zgthr(nz, y, x, indx) Fortran 95:
call gthr(x, indx, y)

?gthrz gathers a sparse vector's elements into compressed form, replacing them by zeros.
Fortran 77:
call sgthrz(nz, y, x, indx) call dgthrz(nz, y, x, indx) call cgthrz(nz, y, x, indx) call zgthrz(nz, y, x, indx) Fortran 95:
call gthrz(x, indx, y)

?roti applies Givens rotation to sparse vectors, one of which is in compressed form (see also ?rot in Level 1 BLAS).
Fortran 77:
call sroti(nz, x, indx, y, c, s) call droti(nz, x, indx, y, c, s) Fortran 95:
call roti(x, indx, y [,c] [,s])

?sctr converts compressed sparse vectors into full storage form.
Fortran 77:
call ssctr(nz, x, indx, y) call dsctr(nz, x, indx, y) call csctr(nz, x, indx, y) call zsctr(nz, x, indx, y) Fortran 95:
call sctr(x, indx, y)