Function

sparse_outer_product_dense_double

Computes the outer product of the dense vector x and the sparse vector y, with both operands containing double-precision values.

Declaration

sparse_status sparse_outer_product_dense_double(sparse_dimension M, sparse_dimension N, sparse_dimension nz, double alpha, const double *restrict x, sparse_stride incx, const double *restrict y, const sparse_index *restrict indy, sparse_matrix_double *restrict C);

Parameters

M

The number of rows of x and the resulting matrix.

N

The number of columns of the resulting matrix. The number of nonzero values must be less than or equal to N.

nz

The number of nonzero values in the sparse vector y. Must be less than or equal to N.

alpha

Scalar multiplier of x.

x

Pointer to the dense vector x. Must be M number of elements. Negative strides are supported. Note, unlike dense BLAS routines, the pointer points to the last element when stride is negative.

incx

Increment between valid values in the dense vector x. Negative strides are supported.

y

Pointer to the dense storage for the values of the sparse vector y. The corresponding entry in indy holds the index of the value. Contains nz values.

indy

Pointer to the dense storage for the index values of the sparse vector y. The corresponding entry in y holds the values of the vector. Contains nz values.

C

Pointer to an uninitialized sparse matrix object. On success a newly allocated sparse matrix object is returned in this pointer. On error, this set to NULL.You are responsible for calling sparse_matrix_destroy on this matrix object.

Return Value

On success SPARSE_SUCCESS is returned an C is valid matrix object. The caller is responsible for cleaning up the sparse matrix object with sparse_matrix_destroy.

Will return SPARSE_ILLEGAL_PARAMETER if nz > N, and C will be unchanged.

Discussion

Compute the outer product of the dense vector x and the sparse vector y and return a new sparse matrix in the uninitialized pointer sparse matrix pointer C. C = alpha * x * y'. You are responsible for calling sparse_matrix_destroy on the returned matrix.The matrix object returned on success is a point wise based sparse matrix.

Indices in indx are always assumed to be stored in ascending order. Additionally, indices are assumed to be unique. The behavior of this function is undefined if either of these assumptions are not met.

All indices are 0 based (the first element of a pointer is ptr[0]).

See Also

Matrix-Vector Operations

sparse_matrix_vector_product_dense_double

Multiplies the dense vector x by the sparse matrix A and adds the result to the dense vector y, with all operands containing double-precision values.

sparse_matrix_vector_product_dense_float

Multiplies the dense vector x by the sparse matrix A and adds the result to the dense vector y, with all operands containing single-precision values.

sparse_vector_triangular_solve_dense_double

Solves the system of equations x = alpha * T⁻¹ * x for x where x is a dense vector and T is a triangular sparse matrix, with all operands containing double-precision values.

sparse_vector_triangular_solve_dense_float

Solves the system of equations x = alpha * T⁻¹ * x for x where x is a dense vector and T is a triangular sparse matrix, with all operands containing single-precision values.

sparse_outer_product_dense_float

Computes the outer product of the dense vector x and the sparse vector y, with both operands containing single-precision values.

sparse_permute_rows_double

Permutes the rows of the double-precision sparse matrix A based on the provided permutation array.

sparse_permute_rows_float

Permutes the rows of the single-precision sparse matrix A based on the provided permutation array.

sparse_permute_cols_double

Permutes the columns of the double-precision sparse matrix A based on the provided permutation array.

sparse_permute_cols_float

Permutes the columns of the single-precision sparse matrix A based on the provided permutation array.

sparse_elementwise_norm_double

Computes the specified element-wise norm of the double-precision sparse matrix A.

sparse_elementwise_norm_float

Computes the specified element-wise norm of the single-precision sparse matrix A.

sparse_operator_norm_double

Computes the specified operator norm of the double-precision sparse matrix A.

sparse_operator_norm_float

Computes the specified operator norm of the single-precision sparse matrix A.

sparse_matrix_trace_double

Computes the sum along the specified diagonal of the double-precision sparse matrix A.

sparse_matrix_trace_float

Computes the sum along the specified diagonal of the single-precision sparse matrix A.