MAGMA  1.5.0
Matrix Algebra for GPU and Multicore Architectures
 All Functions Groups
single-complex precision

Functions

magma_int_t magma_cgecsrmv (magma_trans_t transA, magma_int_t m, magma_int_t n, magmaFloatComplex alpha, magmaFloatComplex *d_val, magma_index_t *d_rowptr, magma_index_t *d_colind, magmaFloatComplex *d_x, magmaFloatComplex beta, magmaFloatComplex *d_y)
 This routine computes y = alpha * A * x + beta * y on the GPU. More...
 
magma_int_t magma_cgecsrmv_shift (magma_trans_t transA, magma_int_t m, magma_int_t n, magmaFloatComplex alpha, magmaFloatComplex lambda, magmaFloatComplex *d_val, magma_index_t *d_rowptr, magma_index_t *d_colind, magmaFloatComplex *d_x, magmaFloatComplex beta, int offset, int blocksize, magma_index_t *add_rows, magmaFloatComplex *d_y)
 This routine computes y = alpha * ( A -lambda I ) * x + beta * y on the GPU. More...
 
magma_int_t magma_cgeellmv (magma_trans_t transA, magma_int_t m, magma_int_t n, magma_int_t nnz_per_row, magmaFloatComplex alpha, magmaFloatComplex *d_val, magma_index_t *d_colind, magmaFloatComplex *d_x, magmaFloatComplex beta, magmaFloatComplex *d_y)
 This routine computes y = alpha * A * x + beta * y on the GPU. More...
 
magma_int_t magma_cgeellmv_shift (magma_trans_t transA, magma_int_t m, magma_int_t n, magma_int_t nnz_per_row, magmaFloatComplex alpha, magmaFloatComplex lambda, magmaFloatComplex *d_val, magma_index_t *d_colind, magmaFloatComplex *d_x, magmaFloatComplex beta, int offset, int blocksize, magma_index_t *add_rows, magmaFloatComplex *d_y)
 This routine computes y = alpha *( A - lambda I ) * x + beta * y on the GPU. More...
 
magma_int_t magma_cgeellrtmv (magma_trans_t transA, magma_int_t m, magma_int_t n, magma_int_t nnz_per_row, magmaFloatComplex alpha, magmaFloatComplex *d_val, magma_index_t *d_colind, magma_index_t *d_rowlength, magmaFloatComplex *d_x, magmaFloatComplex beta, magmaFloatComplex *d_y, magma_int_t alignment, magma_int_t blocksize)
 This routine computes y = alpha * A * x + beta * y on the GPU. More...
 
magma_int_t magma_cgeelltmv_shift (magma_trans_t transA, magma_int_t m, magma_int_t n, magma_int_t nnz_per_row, magmaFloatComplex alpha, magmaFloatComplex lambda, magmaFloatComplex *d_val, magma_index_t *d_colind, magmaFloatComplex *d_x, magmaFloatComplex beta, int offset, int blocksize, magma_index_t *add_rows, magmaFloatComplex *d_y)
 This routine computes y = alpha *( A - lambda I ) * x + beta * y on the GPU. More...
 
magma_int_t magma_cgesellpmv (magma_trans_t transA, magma_int_t m, magma_int_t n, magma_int_t blocksize, magma_int_t slices, magma_int_t alignment, magmaFloatComplex alpha, magmaFloatComplex *d_val, magma_index_t *d_colind, magma_index_t *d_rowptr, magmaFloatComplex *d_x, magmaFloatComplex beta, magmaFloatComplex *d_y)
 This routine computes y = alpha * A^t * x + beta * y on the GPU. More...
 
magma_int_t magma_cgesellcmv (magma_trans_t transA, magma_int_t m, magma_int_t n, magma_int_t blocksize, magma_int_t slices, magma_int_t alignment, magmaFloatComplex alpha, magmaFloatComplex *d_val, magma_index_t *d_colind, magma_index_t *d_rowptr, magmaFloatComplex *d_x, magmaFloatComplex beta, magmaFloatComplex *d_y)
 This routine computes y = alpha * A^t * x + beta * y on the GPU. More...
 
magma_int_t magma_cmdotc (int n, int k, magmaFloatComplex *v, magmaFloatComplex *r, magmaFloatComplex *d1, magmaFloatComplex *d2, magmaFloatComplex *skp)
 Computes the scalar product of a set of vectors v_i such that. More...
 
magma_int_t magma_cmgecsrmv (magma_trans_t transA, magma_int_t m, magma_int_t n, magma_int_t num_vecs, magmaFloatComplex alpha, magmaFloatComplex *d_val, magma_index_t *d_rowptr, magma_index_t *d_colind, magmaFloatComplex *d_x, magmaFloatComplex beta, magmaFloatComplex *d_y)
 This routine computes Y = alpha * A * X + beta * Y for X and Y sets of num_vec vectors on the GPU. More...
 
magma_int_t magma_cmgeellmv (magma_trans_t transA, magma_int_t m, magma_int_t n, magma_int_t num_vecs, magma_int_t nnz_per_row, magmaFloatComplex alpha, magmaFloatComplex *d_val, magma_index_t *d_colind, magmaFloatComplex *d_x, magmaFloatComplex beta, magmaFloatComplex *d_y)
 This routine computes Y = alpha * A * X + beta * Y for X and Y sets of num_vec vectors on the GPU. More...
 
magma_int_t magma_cmgeelltmv (magma_trans_t transA, magma_int_t m, magma_int_t n, magma_int_t num_vecs, magma_int_t nnz_per_row, magmaFloatComplex alpha, magmaFloatComplex *d_val, magma_index_t *d_colind, magmaFloatComplex *d_x, magmaFloatComplex beta, magmaFloatComplex *d_y)
 This routine computes Y = alpha * A * X + beta * Y for X and Y sets of num_vec vectors on the GPU. More...
 
magma_int_t magma_cmgesellpmv (magma_trans_t transA, magma_int_t m, magma_int_t n, magma_int_t num_vecs, magma_int_t blocksize, magma_int_t slices, magma_int_t alignment, magmaFloatComplex alpha, magmaFloatComplex *d_val, magma_index_t *d_colind, magma_index_t *d_rowptr, magmaFloatComplex *d_x, magmaFloatComplex beta, magmaFloatComplex *d_y)
 This routine computes Y = alpha * A^t * X + beta * Y on the GPU. More...
 

Detailed Description

Function Documentation

magma_int_t magma_cgecsrmv ( magma_trans_t  transA,
magma_int_t  m,
magma_int_t  n,
magmaFloatComplex  alpha,
magmaFloatComplex *  d_val,
magma_index_t *  d_rowptr,
magma_index_t *  d_colind,
magmaFloatComplex *  d_x,
magmaFloatComplex  beta,
magmaFloatComplex *  d_y 
)

This routine computes y = alpha * A * x + beta * y on the GPU.

The input format is CSR (val, row, col).

Parameters
transAmagma_trans_t transposition parameter for A
mmagma_int_t number of rows in A
nmagma_int_t number of columns in A
alphamagmaFloatComplex scalar multiplier
d_valmagmaFloatComplex* array containing values of A in CSR
d_rowptrmagma_int_t* rowpointer of A in CSR
d_colindmagma_int_t* columnindices of A in CSR
d_xmagmaFloatComplex* input vector x
betamagmaFloatComplex scalar multiplier
d_ymagmaFloatComplex* input/output vector y
magma_int_t magma_cgecsrmv_shift ( magma_trans_t  transA,
magma_int_t  m,
magma_int_t  n,
magmaFloatComplex  alpha,
magmaFloatComplex  lambda,
magmaFloatComplex *  d_val,
magma_index_t *  d_rowptr,
magma_index_t *  d_colind,
magmaFloatComplex *  d_x,
magmaFloatComplex  beta,
int  offset,
int  blocksize,
magma_index_t *  add_rows,
magmaFloatComplex *  d_y 
)

This routine computes y = alpha * ( A -lambda I ) * x + beta * y on the GPU.

It is a shifted version of the CSR-SpMV.

Parameters
transAmagma_trans_t transposition parameter for A
mmagma_int_t number of rows in A
nmagma_int_t number of columns in A
alphamagmaFloatComplex scalar multiplier
lambdamagmaFloatComplex scalar multiplier
d_valmagmaFloatComplex* array containing values of A in CSR
d_rowptrmagma_int_t* rowpointer of A in CSR
d_colindmagma_int_t* columnindices of A in CSR
d_xmagmaFloatComplex* input vector x
betamagmaFloatComplex scalar multiplier
offsetmagma_int_t in case not the main diagonal is scaled
blocksizemagma_int_t in case of processing multiple vectors
add_rowsmagma_int_t* in case the matrixpowerskernel is used
d_ymagmaFloatComplex* output vector y
magma_int_t magma_cgeellmv ( magma_trans_t  transA,
magma_int_t  m,
magma_int_t  n,
magma_int_t  nnz_per_row,
magmaFloatComplex  alpha,
magmaFloatComplex *  d_val,
magma_index_t *  d_colind,
magmaFloatComplex *  d_x,
magmaFloatComplex  beta,
magmaFloatComplex *  d_y 
)

This routine computes y = alpha * A * x + beta * y on the GPU.

Input format is ELLPACK.

Parameters
transAmagma_trans_t transposition parameter for A
mmagma_int_t number of rows in A
nmagma_int_t number of columns in A
nnz_per_rowmagma_int_t number of elements in the longest row
alphamagmaFloatComplex scalar multiplier
d_valmagmaFloatComplex* array containing values of A in ELLPACK
d_colindmagma_int_t* columnindices of A in ELLPACK
d_xmagmaFloatComplex* input vector x
betamagmaFloatComplex scalar multiplier
d_ymagmaFloatComplex* input/output vector y
magma_int_t magma_cgeellmv_shift ( magma_trans_t  transA,
magma_int_t  m,
magma_int_t  n,
magma_int_t  nnz_per_row,
magmaFloatComplex  alpha,
magmaFloatComplex  lambda,
magmaFloatComplex *  d_val,
magma_index_t *  d_colind,
magmaFloatComplex *  d_x,
magmaFloatComplex  beta,
int  offset,
int  blocksize,
magma_index_t *  add_rows,
magmaFloatComplex *  d_y 
)

This routine computes y = alpha *( A - lambda I ) * x + beta * y on the GPU.

Input format is ELLPACK. It is the shifted version of the ELLPACK SpMV.

Parameters
transAmagma_trans_t transposition parameter for A
mmagma_int_t number of rows in A
nmagma_int_t number of columns in A
nnz_per_rowmagma_int_t number of elements in the longest row
alphamagmaFloatComplex scalar multiplier
lambdamagmaFloatComplex scalar multiplier
d_valmagmaFloatComplex* array containing values of A in ELLPACK
d_colindmagma_int_t* columnindices of A in ELLPACK
d_xmagmaFloatComplex* input vector x
betamagmaFloatComplex scalar multiplier
offsetmagma_int_t in case not the main diagonal is scaled
blocksizemagma_int_t in case of processing multiple vectors
add_rowsmagma_int_t* in case the matrixpowerskernel is used
d_ymagmaFloatComplex* input/output vector y
magma_int_t magma_cgeellrtmv ( magma_trans_t  transA,
magma_int_t  m,
magma_int_t  n,
magma_int_t  nnz_per_row,
magmaFloatComplex  alpha,
magmaFloatComplex *  d_val,
magma_index_t *  d_colind,
magma_index_t *  d_rowlength,
magmaFloatComplex *  d_x,
magmaFloatComplex  beta,
magmaFloatComplex *  d_y,
magma_int_t  alignment,
magma_int_t  blocksize 
)

This routine computes y = alpha * A * x + beta * y on the GPU.

Input format is ELLRT. The ideas are taken from "Improving the performance of the sparse matrix vector product with GPUs", (CIT 2010), and modified to provide correct values.

Parameters
transAmagma_trans_t transposition parameter for A
mmagma_int_t number of rows
nmagma_int_t number of columns
nnz_per_rowmagma_int_t max number of nonzeros in a row
alphamagmaFloatComplex scalar alpha
d_valmagmaFloatComplex* val array
d_colindmagma_int_t* col indices
d_rowlengthmagma_int_t* number of elements in each row
d_xmagmaFloatComplex* input vector x
betamagmaFloatComplex scalar beta
d_ymagmaFloatComplex* output vector y
blocksizemagma_int_t threads per block
alignmentmagma_int_t threads assigned to each row
magma_int_t magma_cgeelltmv_shift ( magma_trans_t  transA,
magma_int_t  m,
magma_int_t  n,
magma_int_t  nnz_per_row,
magmaFloatComplex  alpha,
magmaFloatComplex  lambda,
magmaFloatComplex *  d_val,
magma_index_t *  d_colind,
magmaFloatComplex *  d_x,
magmaFloatComplex  beta,
int  offset,
int  blocksize,
magma_index_t *  add_rows,
magmaFloatComplex *  d_y 
)

This routine computes y = alpha *( A - lambda I ) * x + beta * y on the GPU.

Input format is ELL.

Parameters
transAmagma_trans_t transposition parameter for A
mmagma_int_t number of rows in A
nmagma_int_t number of columns in A
nnz_per_rowmagma_int_t number of elements in the longest row
alphamagmaFloatComplex scalar multiplier
lambdamagmaFloatComplex scalar multiplier
d_valmagmaFloatComplex* array containing values of A in ELL
d_colindmagma_int_t* columnindices of A in ELL
d_xmagmaFloatComplex* input vector x
betamagmaFloatComplex scalar multiplier
offsetmagma_int_t in case not the main diagonal is scaled
blocksizemagma_int_t in case of processing multiple vectors
add_rowsmagma_int_t* in case the matrixpowerskernel is used
d_ymagmaFloatComplex* input/output vector y
magma_int_t magma_cgesellcmv ( magma_trans_t  transA,
magma_int_t  m,
magma_int_t  n,
magma_int_t  blocksize,
magma_int_t  slices,
magma_int_t  alignment,
magmaFloatComplex  alpha,
magmaFloatComplex *  d_val,
magma_index_t *  d_colind,
magma_index_t *  d_rowptr,
magmaFloatComplex *  d_x,
magmaFloatComplex  beta,
magmaFloatComplex *  d_y 
)

This routine computes y = alpha * A^t * x + beta * y on the GPU.

Input format is SELLC/SELLP.

Parameters
transAmagma_trans_t transposition parameter for A
mmagma_int_t number of rows in A
nmagma_int_t number of columns in A
blocksizemagma_int_t number of rows in one ELL-slice
slicesmagma_int_t number of slices in matrix
alignmentmagma_int_t number of threads assigned to one row (=1)
alphamagmaFloatComplex scalar multiplier
d_valmagmaFloatComplex* array containing values of A in SELLC/P
d_colindmagma_int_t* columnindices of A in SELLC/P
d_rowptrmagma_int_t* rowpointer of SELLP
d_xmagmaFloatComplex* input vector x
betamagmaFloatComplex scalar multiplier
d_ymagmaFloatComplex* input/output vector y
magma_int_t magma_cgesellpmv ( magma_trans_t  transA,
magma_int_t  m,
magma_int_t  n,
magma_int_t  blocksize,
magma_int_t  slices,
magma_int_t  alignment,
magmaFloatComplex  alpha,
magmaFloatComplex *  d_val,
magma_index_t *  d_colind,
magma_index_t *  d_rowptr,
magmaFloatComplex *  d_x,
magmaFloatComplex  beta,
magmaFloatComplex *  d_y 
)

This routine computes y = alpha * A^t * x + beta * y on the GPU.

Input format is SELLP.

Parameters
transAmagma_trans_t transposition parameter for A
mmagma_int_t number of rows in A
nmagma_int_t number of columns in A
blocksizemagma_int_t number of rows in one ELL-slice
slicesmagma_int_t number of slices in matrix
alignmentmagma_int_t number of threads assigned to one row
alphamagmaFloatComplex scalar multiplier
d_valmagmaFloatComplex* array containing values of A in SELLP
d_colindmagma_int_t* columnindices of A in SELLP
d_rowptrmagma_int_t* rowpointer of SELLP
d_xmagmaFloatComplex* input vector x
betamagmaFloatComplex scalar multiplier
d_ymagmaFloatComplex* input/output vector y
magma_int_t magma_cmdotc ( int  n,
int  k,
magmaFloatComplex *  v,
magmaFloatComplex *  r,
magmaFloatComplex *  d1,
magmaFloatComplex *  d2,
magmaFloatComplex *  skp 
)

Computes the scalar product of a set of vectors v_i such that.

skp = ( <v_0,r>, <v_1,r>, .. )

Returns the vector skp.

Parameters
nint length of v_i and r
kint

vectors v_i

Parameters
vmagmaFloatComplex* v = (v_0 .. v_i.. v_k)
rmagmaFloatComplex* r
d1magmaFloatComplex* workspace
d2magmaFloatComplex* workspace
skpmagmaFloatComplex* vector[k] of scalar products (<v_i,r>...)
magma_int_t magma_cmgecsrmv ( magma_trans_t  transA,
magma_int_t  m,
magma_int_t  n,
magma_int_t  num_vecs,
magmaFloatComplex  alpha,
magmaFloatComplex *  d_val,
magma_index_t *  d_rowptr,
magma_index_t *  d_colind,
magmaFloatComplex *  d_x,
magmaFloatComplex  beta,
magmaFloatComplex *  d_y 
)

This routine computes Y = alpha * A * X + beta * Y for X and Y sets of num_vec vectors on the GPU.

Input format is CSR.

Parameters
transAmagma_trans_t transposition parameter for A
mmagma_int_t number of rows in A
nmagma_int_t number of columns in A
num_vecsmama_int_t number of vectors
alphamagmaFloatComplex scalar multiplier
d_valmagmaFloatComplex* array containing values of A in CSR
d_rowptrmagma_int_t* rowpointer of A in CSR
d_colindmagma_int_t* columnindices of A in CSR
d_xmagmaFloatComplex* input vector x
betamagmaFloatComplex scalar multiplier
d_ymagmaFloatComplex* input/output vector y
magma_int_t magma_cmgeellmv ( magma_trans_t  transA,
magma_int_t  m,
magma_int_t  n,
magma_int_t  num_vecs,
magma_int_t  nnz_per_row,
magmaFloatComplex  alpha,
magmaFloatComplex *  d_val,
magma_index_t *  d_colind,
magmaFloatComplex *  d_x,
magmaFloatComplex  beta,
magmaFloatComplex *  d_y 
)

This routine computes Y = alpha * A * X + beta * Y for X and Y sets of num_vec vectors on the GPU.

Input format is ELLPACK.

Parameters
transAmagma_trans_t transposition parameter for A
mmagma_int_t number of rows in A
nmagma_int_t number of columns in A
num_vecsmama_int_t number of vectors
nnz_per_rowmagma_int_t number of elements in the longest row
alphamagmaFloatComplex scalar multiplier
d_valmagmaFloatComplex* array containing values of A in ELLPACK
d_colindmagma_int_t* columnindices of A in ELLPACK
d_xmagmaFloatComplex* input vector x
betamagmaFloatComplex scalar multiplier
d_ymagmaFloatComplex* input/output vector y
magma_int_t magma_cmgeelltmv ( magma_trans_t  transA,
magma_int_t  m,
magma_int_t  n,
magma_int_t  num_vecs,
magma_int_t  nnz_per_row,
magmaFloatComplex  alpha,
magmaFloatComplex *  d_val,
magma_index_t *  d_colind,
magmaFloatComplex *  d_x,
magmaFloatComplex  beta,
magmaFloatComplex *  d_y 
)

This routine computes Y = alpha * A * X + beta * Y for X and Y sets of num_vec vectors on the GPU.

Input format is ELL.

Parameters
transAmagma_trans_t transposition parameter for A
mmagma_int_t number of rows in A
nmagma_int_t number of columns in A
num_vecsmama_int_t number of vectors
nnz_per_rowmagma_int_t number of elements in the longest row
alphamagmaFloatComplex scalar multiplier
d_valmagmaFloatComplex* array containing values of A in ELL
d_colindmagma_int_t* columnindices of A in ELL
d_xmagmaFloatComplex* input vector x
betamagmaFloatComplex scalar multiplier
d_ymagmaFloatComplex* input/output vector y
magma_int_t magma_cmgesellpmv ( magma_trans_t  transA,
magma_int_t  m,
magma_int_t  n,
magma_int_t  num_vecs,
magma_int_t  blocksize,
magma_int_t  slices,
magma_int_t  alignment,
magmaFloatComplex  alpha,
magmaFloatComplex *  d_val,
magma_index_t *  d_colind,
magma_index_t *  d_rowptr,
magmaFloatComplex *  d_x,
magmaFloatComplex  beta,
magmaFloatComplex *  d_y 
)

This routine computes Y = alpha * A^t * X + beta * Y on the GPU.

Input format is SELLP. Note, that the input format for X is row-major while the output format for Y is column major!

Parameters
transAmagma_trans_t transpose A?
mmagma_int_t number of rows in A
nmagma_int_t number of columns in A
num_vecsmagma_int_t number of columns in X and Y
blocksizemagma_int_t number of rows in one ELL-slice
slicesmagma_int_t number of slices in matrix
alignmentmagma_int_t number of threads assigned to one row
alphamagmaFloatComplex scalar multiplier
d_valmagmaFloatComplex* array containing values of A in SELLP
d_colindmagma_int_t* columnindices of A in SELLP
d_rowptrmagma_int_t* rowpointer of SELLP
d_xmagmaFloatComplex* input vector x
betamagmaFloatComplex scalar multiplier
d_ymagmaFloatComplex* input/output vector y