Parallel Colt 0.7.2

cern.colt.matrix.tfloat.impl
Class CCFloatMatrix2D

java.lang.Object
  extended by cern.colt.PersistentObject
      extended by cern.colt.matrix.AbstractMatrix
          extended by cern.colt.matrix.AbstractMatrix2D
              extended by cern.colt.matrix.tfloat.FloatMatrix2D
                  extended by cern.colt.matrix.tfloat.impl.WrapperFloatMatrix2D
                      extended by cern.colt.matrix.tfloat.impl.CCFloatMatrix2D
All Implemented Interfaces:
Serializable, Cloneable

public class CCFloatMatrix2D
extends WrapperFloatMatrix2D

Sparse row-compressed 2-d matrix holding float elements. First see the package summary and javadoc tree view to get the broad picture.

Implementation:

Internally uses the standard sparse row-compressed format, with two important differences that broaden the applicability of this storage format:


Note that this implementation is not synchronized.

Memory requirements:

Cells that

memory [bytes] = 4*rows + 12 * nonZeros.
Where nonZeros = cardinality() is the number of non-zero cells. Thus, a 1000 x 1000 matrix with 1000000 non-zero cells consumes 11.5 MB. The same 1000 x 1000 matrix with 1000 non-zero cells consumes 15 KB.

Time complexity:

Getting a cell value takes time O(log nzr) where nzr is the number of non-zeros of the touched row. This is usually quick, because typically there are only few nonzeros per row. So, in practice, get has expected constant time. Setting a cell value takes worst-case time O(nz) where nzr is the total number of non-zeros in the matrix. This can be extremely slow, but if you traverse coordinates properly (i.e. upwards), each write is done much quicker:

 // rather quick
 matrix.assign(0);
 for (int row = 0; row < rows; row++) {
     for (int column = 0; column < columns; column++) {
         if (someCondition)
             matrix.setQuick(row, column, someValue);
     }
 }
 
 // poor
 matrix.assign(0);
 for (int row = rows; --row >= 0;) {
     for (int column = columns; --column >= 0;) {
         if (someCondition)
             matrix.setQuick(row, column, someValue);
     }
 }
 
If for whatever reasons you can't iterate properly, consider to create an empty dense matrix, store your non-zeros in it, then call sparse.assign(dense). Under the circumstances, this is still rather quick.

Fast iteration over non-zeros can be done via forEachNonZero(cern.colt.function.tfloat.IntIntFloatFunction), which supplies your function with row, column and value of each nonzero. Although the internally implemented version is a bit more sophisticated, here is how a quite efficient user-level matrix-vector multiplication could look like:

 // Linear algebraic y = A * x
 A.forEachNonZero(new cern.colt.function.IntIntFloatFunction() {
     public float apply(int row, int column, float value) {
         y.setQuick(row, y.getQuick(row) + value * x.getQuick(column));
         return value;
     }
 });
 

Here is how a a quite efficient user-level combined scaling operation could look like:

 // Elementwise A = A + alpha*B
 B.forEachNonZero(new cern.colt.function.IntIntFloatFunction() {
     public float apply(int row, int column, float value) {
         A.setQuick(row, column, A.getQuick(row, column) + alpha * value);
         return value;
     }
 });
 
Method assign(FloatMatrix2D,cern.colt.function.tfloat.FloatFloatFunction) does just that if you supply FloatFunctions.plusMultSecond(float) as argument.

Version:
0.9, 04/14/2000
Author:
wolfgang.hoschek@cern.ch
See Also:
Serialized Form

Field Summary
 
Fields inherited from class cern.colt.PersistentObject
serialVersionUID
 
Constructor Summary
CCFloatMatrix2D(float[][] values)
          Constructs a matrix with a copy of the given values.
CCFloatMatrix2D(int rows, int columns)
          Constructs a matrix with a given number of rows and columns.
CCFloatMatrix2D(int rows, int columns, int nzmax)
          Constructs a matrix with a given number of rows and columns.
CCFloatMatrix2D(int rows, int columns, int[] columnPointers, IntArrayList rowIndexes, FloatArrayList values)
          Constructs a matrix with indexes and values given in compressed-column format.
CCFloatMatrix2D(int rows, int columns, IntArrayList rowIndexes, IntArrayList columnIndexes, float value)
          Constructs a matrix with indexes given in coordinate format and single value.
CCFloatMatrix2D(int rows, int columns, IntArrayList rowIndexes, IntArrayList columnIndexes, FloatArrayList values)
          Constructs a matrix with indexes and values given in coordinate format.
 
Method Summary
 FloatMatrix2D assign(float value)
          Sets all nonzero cells to the state specified by value.
 FloatMatrix2D assign(FloatFunction function)
          Assigns the result of a function to each nonzero cell;
 FloatMatrix2D assign(FloatMatrix2D source)
          Replaces all cell values of the receiver with the values of another matrix.
 FloatMatrix2D assign(FloatMatrix2D y, FloatFloatFunction function)
          Assigns the result of a function to each cell; x[row,col] = function(x[row,col],y[row,col]).
 int cardinality()
          Returns the number of cells having non-zero values; ignores tolerance.
 FloatMatrix2D forEachNonZero(IntIntFloatFunction function)
          Assigns the result of a function to each non-zero cell; x[row,col] = function(x[row,col]).
 int[] getColumnPointers()
           
 DenseFloatMatrix2D getFull()
           
 float getQuick(int row, int column)
          Returns the matrix cell value at coordinate [row,column].
 IntArrayList getRowindexes()
           
 FloatArrayList getValues()
           
 FloatMatrix2D like(int rows, int columns)
          Construct and returns a new empty matrix of the same dynamic type as the receiver, having the specified number of rows and columns.
 FloatMatrix1D like1D(int size)
          Construct and returns a new 1-d matrix of the corresponding dynamic type, entirely independent of the receiver.
 void setQuick(int row, int column, float value)
          Sets the matrix cell at coordinate [row,column] to the specified value.
 void trimToSize()
          Releases any superfluous internal memory.
 FloatMatrix1D zMult(FloatMatrix1D y, FloatMatrix1D z, float alpha, float beta, boolean transposeA)
          Linear algebraic matrix-vector multiplication; z = alpha * A * y + beta*z.
 FloatMatrix2D zMult(FloatMatrix2D B, FloatMatrix2D C, float alpha, float beta, boolean transposeA, boolean transposeB)
          Linear algebraic matrix-matrix multiplication; C = alpha * A x B + beta*C.
 
Methods inherited from class cern.colt.matrix.tfloat.impl.WrapperFloatMatrix2D
elements, equals, equals, setNcolumns, setNrows, vectorize, viewColumn, viewColumnFlip, viewDice, viewPart, viewRow, viewRowFlip, viewSelection, viewStrides
 
Methods inherited from class cern.colt.matrix.tfloat.FloatMatrix2D
aggregate, aggregate, aggregate, aggregate, assign, assign, assign, assign, assign, copy, get, getMaxLocation, getMinLocation, getNegativeValues, getNonZeros, getPositiveValues, like, normalize, set, toArray, toString, viewSelection, viewSelection, viewSorted, zAssign8Neighbors, zMult, zMult, zSum
 
Methods inherited from class cern.colt.matrix.AbstractMatrix2D
checkShape, checkShape, columns, columnStride, index, rows, rowStride, size, toStringShort
 
Methods inherited from class cern.colt.matrix.AbstractMatrix
ensureCapacity, isView
 
Methods inherited from class cern.colt.PersistentObject
clone
 
Methods inherited from class java.lang.Object
getClass, hashCode, notify, notifyAll, wait, wait, wait
 

Constructor Detail

CCFloatMatrix2D

public CCFloatMatrix2D(float[][] values)
Constructs a matrix with a copy of the given values. values is required to have the form values[row][column] and have exactly the same number of columns in every row.

The values are copied. So subsequent changes in values are not reflected in the matrix, and vice-versa.

Parameters:
values - The values to be filled into the new matrix.
Throws:
IllegalArgumentException - if for any 1 <= row < values.length: values[row].length != values[row-1].length .

CCFloatMatrix2D

public CCFloatMatrix2D(int rows,
                       int columns,
                       int[] columnPointers,
                       IntArrayList rowIndexes,
                       FloatArrayList values)
Constructs a matrix with indexes and values given in compressed-column format.

Parameters:
rows -
columns -
columnPointers -
rowIndexes -
values -

CCFloatMatrix2D

public CCFloatMatrix2D(int rows,
                       int columns)
Constructs a matrix with a given number of rows and columns. All entries are initially 0.

Parameters:
rows - the number of rows the matrix shall have.
columns - the number of columns the matrix shall have.
Throws:
IllegalArgumentException - if rows<0 || columns<0 || (float)columns*rows > Integer.MAX_VALUE .

CCFloatMatrix2D

public CCFloatMatrix2D(int rows,
                       int columns,
                       int nzmax)
Constructs a matrix with a given number of rows and columns. All entries are initially 0.

Parameters:
rows - the number of rows the matrix shall have.
columns - the number of columns the matrix shall have.
Throws:
IllegalArgumentException - if rows<0 || columns<0 || (float)columns*rows > Integer.MAX_VALUE .

CCFloatMatrix2D

public CCFloatMatrix2D(int rows,
                       int columns,
                       IntArrayList rowIndexes,
                       IntArrayList columnIndexes,
                       FloatArrayList values)
Constructs a matrix with indexes and values given in coordinate format.

Parameters:
rows -
columns -
rowIndexes -
columnIndexes -
values -

CCFloatMatrix2D

public CCFloatMatrix2D(int rows,
                       int columns,
                       IntArrayList rowIndexes,
                       IntArrayList columnIndexes,
                       float value)
Constructs a matrix with indexes given in coordinate format and single value.

Parameters:
rows -
columns -
rowIndexes -
columnIndexes -
value -
Method Detail

assign

public FloatMatrix2D assign(float value)
Sets all nonzero cells to the state specified by value.

Overrides:
assign in class FloatMatrix2D
Parameters:
value - the value to be filled into the cells.
Returns:
this (for convenience only).

assign

public FloatMatrix2D assign(FloatFunction function)
Assigns the result of a function to each nonzero cell;

Overrides:
assign in class FloatMatrix2D
Parameters:
function - a function object taking as argument the current cell's value.
Returns:
this (for convenience only).
See Also:
FloatFunctions

assign

public FloatMatrix2D assign(FloatMatrix2D source)
Replaces all cell values of the receiver with the values of another matrix. Both matrices must have the same number of rows and columns. If both matrices share the same cells (as is the case if they are views derived from the same matrix) and intersect in an ambiguous way, then replaces as if using an intermediate auxiliary deep copy of other.

Overrides:
assign in class FloatMatrix2D
Parameters:
source - the source matrix to copy from (may be identical to the receiver).
Returns:
this (for convenience only).
Throws:
IllegalArgumentException - if columns() != source.columns() || rows() != source.rows()

assign

public FloatMatrix2D assign(FloatMatrix2D y,
                            FloatFloatFunction function)
Description copied from class: FloatMatrix2D
Assigns the result of a function to each cell; x[row,col] = function(x[row,col],y[row,col]).

Example:

   // assign x[row,col] = x[row,col]<sup>y[row,col]</sup>
   m1 = 2 x 2 matrix 
   0 1 
   2 3
 
   m2 = 2 x 2 matrix 
   0 2 
   4 6
 
   m1.assign(m2, cern.jet.math.Functions.pow);
   -->
   m1 == 2 x 2 matrix
   1   1 
   16 729
 
 
For further examples, see the package doc.

Overrides:
assign in class FloatMatrix2D
Parameters:
y - the secondary matrix to operate on.
function - a function object taking as first argument the current cell's value of this, and as second argument the current cell's value of y,
Returns:
this (for convenience only).
See Also:
FloatFunctions

forEachNonZero

public FloatMatrix2D forEachNonZero(IntIntFloatFunction function)
Description copied from class: FloatMatrix2D
Assigns the result of a function to each non-zero cell; x[row,col] = function(x[row,col]). Use this method for fast special-purpose iteration. If you want to modify another matrix instead of this (i.e. work in read-only mode), simply return the input value unchanged. Parameters to function are as follows: first==row, second==column, third==nonZeroValue.

Overrides:
forEachNonZero in class FloatMatrix2D
Parameters:
function - a function object taking as argument the current non-zero cell's row, column and value.
Returns:
this (for convenience only).

getRowindexes

public IntArrayList getRowindexes()

getColumnPointers

public int[] getColumnPointers()

getValues

public FloatArrayList getValues()

getQuick

public float getQuick(int row,
                      int column)
Returns the matrix cell value at coordinate [row,column].

Provided with invalid parameters this method may return invalid objects without throwing any exception. You should only use this method when you are absolutely sure that the coordinate is within bounds. Precondition (unchecked): 0 <= column < columns() && 0 <= row < rows().

Overrides:
getQuick in class WrapperFloatMatrix2D
Parameters:
row - the index of the row-coordinate.
column - the index of the column-coordinate.
Returns:
the value at the specified coordinate.

like

public FloatMatrix2D like(int rows,
                          int columns)
Construct and returns a new empty matrix of the same dynamic type as the receiver, having the specified number of rows and columns. For example, if the receiver is an instance of type DenseFloatMatrix2D the new matrix must also be of type DenseFloatMatrix2D, if the receiver is an instance of type SparseFloatMatrix2D the new matrix must also be of type SparseFloatMatrix2D, etc. In general, the new matrix should have internal parametrization as similar as possible.

Overrides:
like in class WrapperFloatMatrix2D
Parameters:
rows - the number of rows the matrix shall have.
columns - the number of columns the matrix shall have.
Returns:
a new empty matrix of the same dynamic type.

like1D

public FloatMatrix1D like1D(int size)
Construct and returns a new 1-d matrix of the corresponding dynamic type, entirely independent of the receiver. For example, if the receiver is an instance of type DenseFloatMatrix2D the new matrix must be of type DenseFloatMatrix1D, if the receiver is an instance of type SparseFloatMatrix2D the new matrix must be of type SparseFloatMatrix1D, etc.

Overrides:
like1D in class WrapperFloatMatrix2D
Parameters:
size - the number of cells the matrix shall have.
Returns:
a new matrix of the corresponding dynamic type.

cardinality

public int cardinality()
Description copied from class: FloatMatrix2D
Returns the number of cells having non-zero values; ignores tolerance.

Overrides:
cardinality in class FloatMatrix2D
Returns:
cardinality

setQuick

public void setQuick(int row,
                     int column,
                     float value)
Sets the matrix cell at coordinate [row,column] to the specified value.

Provided with invalid parameters this method may access illegal indexes without throwing any exception. You should only use this method when you are absolutely sure that the coordinate is within bounds. Precondition (unchecked): 0 <= column < columns() && 0 <= row < rows().

Overrides:
setQuick in class WrapperFloatMatrix2D
Parameters:
row - the index of the row-coordinate.
column - the index of the column-coordinate.
value - the value to be filled into the specified cell.

getFull

public DenseFloatMatrix2D getFull()

trimToSize

public void trimToSize()
Description copied from class: AbstractMatrix
Releases any superfluous internal memory. An application can use this operation to minimize the storage of the receiver.

This default implementation does nothing. Override this method if necessary.

Overrides:
trimToSize in class AbstractMatrix

zMult

public FloatMatrix1D zMult(FloatMatrix1D y,
                           FloatMatrix1D z,
                           float alpha,
                           float beta,
                           boolean transposeA)
Description copied from class: FloatMatrix2D
Linear algebraic matrix-vector multiplication; z = alpha * A * y + beta*z. z[i] = alpha*Sum(A[i,j] * y[j]) + beta*z[i], i=0..A.rows()-1, j=0..y.size()-1 . Where A == this.
Note: Matrix shape conformance is checked after potential transpositions.

Overrides:
zMult in class FloatMatrix2D
Parameters:
y - the source vector.
z - the vector where results are to be stored. Set this parameter to null to indicate that a new result vector shall be constructed.
Returns:
z (for convenience only).

zMult

public FloatMatrix2D zMult(FloatMatrix2D B,
                           FloatMatrix2D C,
                           float alpha,
                           float beta,
                           boolean transposeA,
                           boolean transposeB)
Description copied from class: FloatMatrix2D
Linear algebraic matrix-matrix multiplication; C = alpha * A x B + beta*C. C[i,j] = alpha*Sum(A[i,k] * B[k,j]) + beta*C[i,j], k=0..n-1.
Matrix shapes: A(m x n), B(n x p), C(m x p).
Note: Matrix shape conformance is checked after potential transpositions.

Overrides:
zMult in class FloatMatrix2D
Parameters:
B - the second source matrix.
C - the matrix where results are to be stored. Set this parameter to null to indicate that a new result matrix shall be constructed.
Returns:
C (for convenience only).

Parallel Colt 0.7.2

Jump to the Parallel Colt Homepage