From 9fcad57af8346bed0c802c652be408dd67d23e23 Mon Sep 17 00:00:00 2001 From: Matthew Fernandez Date: Wed, 30 Jun 2021 20:01:51 -0700 Subject: [PATCH] remove now unused SparseMatrix_multiply_dense2 --- lib/sparse/SparseMatrix.c | 34 ---------------------------------- 1 file changed, 34 deletions(-) diff --git a/lib/sparse/SparseMatrix.c b/lib/sparse/SparseMatrix.c index 0bf5cf316..366e18a9f 100644 --- a/lib/sparse/SparseMatrix.c +++ b/lib/sparse/SparseMatrix.c @@ -1055,40 +1055,6 @@ static void SparseMatrix_multiply_dense1(SparseMatrix A, real *v, real **res, in *res = u; -} - -static void SparseMatrix_multiply_dense2(SparseMatrix A, real *v, real **res, int dim, int transposed, int res_transposed){ - /* A v^T or A^T v^T where v a dense matrix of second dimension n or m. Real only for now. - transposed = FALSE: A*V^T, with A dimension m x n, V dimension dim x n, v[i*n+j] gives V[i,j]. Result of dimension m x dim - transposed = TRUE: A^T*V^T, with A dimension m x n, V dimension dim x m. v[i*m+j] gives V[i,j]. Result of dimension n x dim - */ - real *u, *rr; - int i, m, n; - assert(A->format == FORMAT_CSR); - assert(A->type == MATRIX_TYPE_REAL); - u = *res; - m = A->m; - n = A->n; - - if (!transposed){ - if (!u) u = MALLOC(sizeof(real)*((size_t)m)*((size_t) dim)); - for (i = 0; i < dim; i++){ - rr = &(u[m*i]); - SparseMatrix_multiply_vector(A, &(v[n*i]), &rr, transposed); - } - if (!res_transposed) dense_transpose(u, dim, m); - } else { - if (!u) u = MALLOC(sizeof(real)*((size_t)n)*((size_t)dim)); - for (i = 0; i < dim; i++){ - rr = &(u[n*i]); - SparseMatrix_multiply_vector(A, &(v[m*i]), &rr, transposed); - } - if (!res_transposed) dense_transpose(u, dim, n); - } - - *res = u; - - } void SparseMatrix_multiply_dense(SparseMatrix A, real *v, real **res, int res_transposed, int dim){ -- 2.50.1