4

Is there native support for Compressed Sparse Column (CSC) format for sparse matrices, like importing and manipulating them?

rcollyer
  • 33,976
  • 7
  • 92
  • 191
my account_ram
  • 2,158
  • 11
  • 25

3 Answers3

9

Mathematica has, I think, very efficient SparseArray functionality. SarseArrays are supported by a wide variety of functions, e.g LinearSolve, Eigensystem... The SparseArray is stored in CSR (Look for the SparseArray entry). If you are interested here is a description of the internals. For the Import/Export you could either make use of the HarwellBoeing format (as pointed out by Szabolcs) or make use of the positions and values and then create the SparseArray using:

positions = {{1, 1}, {2, 2}, {3, 3}, {1, 3}};
values = {1, 2, 3, 4};
sa = SparseArray[positions -> values]

To extract positions and values from SparseArrays use:

sa["NonzeroPositions"]
(*{{1, 1}, {1, 3}, {2, 2}, {3, 3}}*)
sa["NonzeroValues"]
(*{1, 4, 2, 3}*)

and then export those. "Properties" has more info.

Here is a link where examples of efficient use of SparseArrays was discussed and it shows the addition of values into the same position. But there are many more cases where SparseArrays are useful. In fact, I think, they are quite underestimated.

  • 1
    +1. I was not aware of all these Mathgroup discussions. Another possibly relevant link in this context: http://mathematica.stackexchange.com/questions/276/efficient-way-to-combine-sparsearray-objects/287#287 – Leonid Shifrin May 11 '12 at 09:00
6

The first Google hit for "compressed sparse column" is this page, which explains that it's also called the Harwell-Boeing format. If this is what you need, this format is supported by Mathematica.

Szabolcs
  • 234,956
  • 30
  • 623
  • 1,263
1

You can view the full list of supported formats. You can also view the list of formats by category.

You can evaluate $ImportFormats or $ExportFormats for a simple list of formats available.

Mr.Wizard
  • 271,378
  • 34
  • 587
  • 1,371