About Site Map Submit Contact Us Log in | Create an account
Create an account Log In
Average Rating
User Rating:
Visitors Rating:
My rating:

Write review
  • Last update: 5 years ago
  • Total downloads: 21
  • Operating system: Linux, Mac OS X, Win All, BSD, Solaris
  • Publisher: Matt J (15 other programs)
See full specifications

scripts default iconRobust Sparse data types (Scripts) Publisher's description

The main file in this package, rsparse.m

The main file in this package, rsparse.m, creates a sparse array object that is a subclass of MATLAB's native sparse data type, but which is robust against certain bugs and errors (present as of R2009b).


The I/O syntax and options for rsparse() are identical to that of sparse().

Where sparse() would return a double sparse matrix, rsparse() will return a matrix of class RobustSparseDouble, a class defined in this package. Similarly, wherever sparse() would return a logical sparse matrix, rsparse() will return a matrix of user-defined class RobustSparseLogical.

EXAMPLE 1: Indexing robustness - this utilizes the sparse-sub-access package by Bruno Luong.

>> A=sparse(100000,100000,1), %native MATLAB sparse matrix

A =

(100000,100000) 1

>> B=rsparse(A), %convert to robust type


(100000,100000) 1

>> whos A B

Name Size Bytes Class Attributes

A 100000x100000 400016 double sparse
B 100000x100000 400072 RobustSparseDouble

>> A(logical(A))=2, %Wrong result due to indexing bug

A =

(65408,14101) 2
(100000,100000) 1

>> B(logical(B))=2, %RIGHT!!!


(100000,100000) 2

EXAMPLE 2: Robustness to mixed data type operations.

For some reason, MATLAB decides to issue an error when a mathematical operation involving sparse and non-double numeric data is attempted:


>> y=A*x, %ANNOYING....

??? Error using ==> mtimes
Sparse single arrays are not supported.

The RobustSparse types, however, will silently pre-convert the foreign data to doubles:


>>y=B*x, %BETTER!!!

y =



I'm hoping that TMW will fix these problems soon, rendering this package obsolete. Once they do, simply replace the code in rsparse.m with the following, and rsparse() will become equivalent to sparse()

function S=rsparse(varargin)



The package also includes "robustified" versions of other common sparfun toolbox routines: rspeye, rsprand, rsprandn, rprandsym,... They too will revert to their native MATLAB versions when the above replacement is made.


(1) Because this package uses various M-coded wrappers for built-in sparse routines, there is obviously some compromise in speed.

(2) I've overloaded all math operators (+ , - , .* , >=, <=, etc...) and some common functions like sum(), inv(), spfun(), etc... to return RobustSparse data types, where they would otherwise return ordinary sparse types.

For obvious reasons, however, I cannot do this for all functions out there capable of returning a sparse matrix (e.g., triu, tril, sin, cos, etc ... ), nor can I keep up with TMW's releases of new functions.

Therefore, you should bear in mind that B=triu(rsparse(A)) and similar will return an ordinary native sparse data type, B, with all of its vulnerabilities. You should take care to post convert such results to robust type, via B=rsparse(B). Or, you can add overloaded methods as you see fit, mimicking the others in RobustSparseDouble.m and RobustSparseLogical.m

System Requirements:

MATLAB 7.9 (2009b)
Program Release Status: Major Update
Program Install Support: Install and Uninstall

Robust Sparse data types (Scripts) Tags:

Click on a tag to find related softwares

Is Robust Sparse data types (Scripts) your software?

Manage your software

Most Popular

scripts default icon ASK, OOK, FSK, BPSK, QPSK, 8PSK modulation 1.1
ASK, OOK, FSK, BPSK, QPSK, 8PSK modulation contain several functions for digital modulation simulation
scripts default icon Simulink Communication Labs 1.1
Simulink Communication Labs allows you to learn communication systems in greater depth.
scripts default icon M-QAM modulation and demodulation 1.1
M-QAM modulation and demodulation is the QAM modulation and demodulation tech.
scripts default icon LZW Compression/Decompression 1.1
LZW Compression/Decompression - Updated LZW compressor and decompressor with reasonable performance
scripts default icon InSPIRE utility to plot a 2D displacement field (Scripts) 1.0
This program plots the deformation field (displace vectors) contained in vector.txt.