remove unused SparseMatrix_multiply_by_scaler

diff --git a/lib/sparse/SparseMatrix.c b/lib/sparse/SparseMatrix.c
index 34830f809219df7ae201cfae46a271044f1bf466..850afeedf0ad026d944a48ff2ae80c6ff305e99c 100644 (file)
--- a/lib/sparse/SparseMatrix.c
+++ b/lib/sparse/SparseMatrix.c
@@ -1347,52 +1347,6 @@ SparseMatrix SparseMatrix_scaled_by_vector(SparseMatrix A, real *v, int apply_to
   return A;
 
 }
-SparseMatrix SparseMatrix_multiply_by_scaler(SparseMatrix A, real s){
-  /* A scaled by a number */
-  int i, j, *ia, m;
-  real *a, *b = NULL;
-  int *ai;
-  assert(A->format == FORMAT_CSR);
-
-  switch (A->type){
-  case MATRIX_TYPE_INTEGER:
-    b = MALLOC(sizeof(real)*A->nz);
-    ai = (int*) A->a;
-    for (i = 0; i < A->nz; i++) b[i] = ai[i];
-    FREE(A->a);
-    A->a = b;
-    A->type = MATRIX_TYPE_REAL;
-    // fall through
-  case MATRIX_TYPE_REAL:
-    a = (real*) A->a;
-    ia = A->ia;
-    m = A->m;
-    for (i = 0; i < m; i++){
-      for (j = ia[i]; j < ia[i+1]; j++){
-       a[j] *= s;
-      }
-    }
-    break;
-  case MATRIX_TYPE_COMPLEX:
-    a = (real*) A->a;
-    ia = A->ia;
-    m = A->m;
-    for (i = 0; i < m; i++){
-      for (j = ia[i]; j < ia[i+1]; j++){
-       a[2*j] *= s;
-       a[2*j+1] *= s;
-      }
-    }
-
-    break;
-  default:
-    fprintf(stderr,"warning: scaling of matrix this type is not supported\n");
-  }
-
-  return A;
-
-}
-
 
 SparseMatrix SparseMatrix_multiply(SparseMatrix A, SparseMatrix B){
   int m;

diff --git a/lib/sparse/SparseMatrix.h b/lib/sparse/SparseMatrix.h
index 466e3a684123e338690c9ccd71bf4d3f4e938500..2cead463c69cd5fb57f9dec42a4383c87f355b66 100644 (file)
--- a/lib/sparse/SparseMatrix.h
+++ b/lib/sparse/SparseMatrix.h
@@ -97,7 +97,6 @@ int SparseMatrix_has_diagonal(SparseMatrix A);
 SparseMatrix SparseMatrix_normalize_by_row(SparseMatrix A);/* divide by max of each row */
 SparseMatrix SparseMatrix_crop(SparseMatrix A, real epsilon);/*remove any entry <= epsilon*/
 SparseMatrix SparseMatrix_scaled_by_vector(SparseMatrix A, real *v, int apply_to_row);
-SparseMatrix SparseMatrix_multiply_by_scaler(SparseMatrix A, real s);
 SparseMatrix SparseMatrix_make_undirected(SparseMatrix A);/* make it strictly low diag only, and set flag to undirected */
 int SparseMatrix_connectedQ(SparseMatrix A);
 real SparseMatrix_pseudo_diameter_only(SparseMatrix A);