sparse matrix with data type as 'single'

Dear all:

I have a very large sparse matrix that I hope to reduce the size further more. I noticed that in full matrix, one could cast a matrix to a different data type, e.g., from double to single.

A = eye(100);
S = sparse(A); 
B = single(eye(100));
S = sparse(B);  % error: Undefined function 'sparse' for input arguments of type 'single'.
whos

I am wondering if there is any way to turn non-zero part of a sparse matrix into a data type of single so that I could save half of the space.

Thank you, and I look forward to hearing from you!

Best,

Long

Answer 1

[Walter Roberson]

No, that is not supported. sparse() is only supported for double and logical.

[Long Hong]

Thank you, Walter! Do you happen to know any alternative way to make the storage smaller for sparse matrices?

[Walter Roberson]

Not while using the defined sparse functions.

What operations are you doing with the matrices? For some restricted operations it might be practical to write replacement sparse routines. However if you are doing something like eigs() or using sparse \ (matrix division) then enough would have to be reimplemented to make it questionable as to whether it is worth doing.

[Long Hong]

Thank you for your reply, Walter!

I am using the large sparse matrix to do some basic operations, and two advanced ones - Incomplete Cholesky factorization (ichol) and use preconditioned conjugate gradients method to solve system of linear equations (pcg).

Could you please explain a bit more on the "replacement sparse routines"? I am quite new to the sparse world. Thank you!

[Walter Roberson]

https://www.researchgate.net/post/Can_I_perform_sparse_matrix_operations_in_Matlab_with_single_precision

[Long Hong]

Thank you, Walter! I will check this post out!