public abstract class DoubleMatrix3D extends AbstractMatrix3D
A matrix has a number of slices, rows and columns, which are assigned upon instance construction - The matrix's size is then slices()*rows()*columns(). Elements are accessed via [slice,row,column] coordinates. Legal coordinates range from [0,0,0] to [slices()-1,rows()-1,columns()-1]. Any attempt to access an element at a coordinate slice<0 || slice>=slices() || row<0 || row>=rows() || column<0 || column>=column() will throw an IndexOutOfBoundsException.
Note that this implementation is not synchronized.
Modifier and Type | Method and Description |
---|---|
double |
aggregate(DoubleDoubleFunction aggr,
DoubleFunction f)
Applies a function to each cell and aggregates the results.
|
double |
aggregate(DoubleDoubleFunction aggr,
DoubleFunction f,
DoubleProcedure cond)
Applies a function to each cell that satisfies a condition and aggregates
the results.
|
double |
aggregate(DoubleDoubleFunction aggr,
DoubleFunction f,
IntArrayList sliceList,
IntArrayList rowList,
IntArrayList columnList)
Applies a function to all cells with a given indexes and aggregates the
results.
|
double |
aggregate(DoubleMatrix3D other,
DoubleDoubleFunction aggr,
DoubleDoubleFunction f)
Applies a function to each corresponding cell of two matrices and
aggregates the results.
|
DoubleMatrix3D |
assign(double value)
Sets all cells to the state specified by value.
|
DoubleMatrix3D |
assign(double[] values)
Sets all cells to the state specified by values.
|
DoubleMatrix3D |
assign(double[][][] values)
Sets all cells to the state specified by values.
|
DoubleMatrix3D |
assign(DoubleFunction function)
Assigns the result of a function to each cell;
x[slice,row,col] = function(x[slice,row,col]).
|
DoubleMatrix3D |
assign(DoubleMatrix3D other)
Replaces all cell values of the receiver with the values of another
matrix.
|
DoubleMatrix3D |
assign(DoubleMatrix3D y,
DoubleDoubleFunction function)
Assigns the result of a function to each cell;
x[row,col] = function(x[row,col],y[row,col]).
|
DoubleMatrix3D |
assign(DoubleMatrix3D y,
DoubleDoubleFunction function,
IntArrayList sliceList,
IntArrayList rowList,
IntArrayList columnList)
Assigns the result of a function to all cells with a given indexes
|
DoubleMatrix3D |
assign(DoubleProcedure cond,
double value)
Assigns a value to all cells that satisfy a condition.
|
DoubleMatrix3D |
assign(DoubleProcedure cond,
DoubleFunction f)
Assigns the result of a function to all cells that satisfy a condition.
|
int |
cardinality()
Returns the number of cells having non-zero values; ignores tolerance.
|
DoubleMatrix3D |
copy()
Constructs and returns a deep copy of the receiver.
|
abstract Object |
elements()
Returns the elements of this matrix.
|
boolean |
equals(double value)
Returns whether all cells are equal to the given value.
|
boolean |
equals(Object obj)
Compares this object against the specified object.
|
double |
get(int slice,
int row,
int column)
Returns the matrix cell value at coordinate [slice,row,column].
|
double[] |
getMaxLocation()
Return maximum value of this matrix together with its location
|
double[] |
getMinLocation()
Returns minimum value of this matrix together with its location
|
void |
getNegativeValues(IntArrayList sliceList,
IntArrayList rowList,
IntArrayList columnList,
DoubleArrayList valueList)
Fills the coordinates and values of cells having negative values into the
specified lists.
|
void |
getNonZeros(IntArrayList sliceList,
IntArrayList rowList,
IntArrayList columnList,
DoubleArrayList valueList)
Fills the coordinates and values of cells having non-zero values into the
specified lists.
|
void |
getPositiveValues(IntArrayList sliceList,
IntArrayList rowList,
IntArrayList columnList,
DoubleArrayList valueList)
Fills the coordinates and values of cells having positive values into the
specified lists.
|
abstract double |
getQuick(int slice,
int row,
int column)
Returns the matrix cell value at coordinate [slice,row,column].
|
DoubleMatrix3D |
like()
Construct and returns a new empty matrix of the same dynamic type
as the receiver, having the same number of slices, rows and columns.
|
abstract DoubleMatrix3D |
like(int slices,
int rows,
int columns)
Construct and returns a new empty matrix of the same dynamic type
as the receiver, having the specified number of slices, rows and columns.
|
abstract DoubleMatrix2D |
like2D(int rows,
int columns)
Construct and returns a new 2-d matrix of the corresponding dynamic
type, sharing the same cells.
|
void |
normalize()
Normalizes this matrix, i.e.
|
void |
set(int slice,
int row,
int column,
double value)
Sets the matrix cell at coordinate [slice,row,column] to the
specified value.
|
abstract void |
setQuick(int slice,
int row,
int column,
double value)
Sets the matrix cell at coordinate [slice,row,column] to the
specified value.
|
double[][][] |
toArray()
Constructs and returns a 2-dimensional array containing the cell values.
|
String |
toString()
Returns a string representation using default formatting.
|
abstract DoubleMatrix1D |
vectorize()
Returns a vector obtained by stacking the columns of each slice of the
matrix on top of one another.
|
DoubleMatrix2D |
viewColumn(int column)
Constructs and returns a new 2-dimensional slice view representing
the slices and rows of the given column.
|
DoubleMatrix3D |
viewColumnFlip()
Constructs and returns a new flip view along the column axis.
|
DoubleMatrix3D |
viewDice(int axis0,
int axis1,
int axis2)
Constructs and returns a new dice view; Swaps dimensions (axes);
Example: 3 x 4 x 5 matrix --> 4 x 3 x 5 matrix.
|
DoubleMatrix3D |
viewPart(int slice,
int row,
int column,
int depth,
int height,
int width)
Constructs and returns a new sub-range view that is a
depth x height x width sub matrix starting at
[slice,row,column]; Equivalent to
view().part(slice,row,column,depth,height,width); Provided for
convenience only.
|
DoubleMatrix2D |
viewRow(int row)
Constructs and returns a new 2-dimensional slice view representing
the slices and columns of the given row.
|
DoubleMatrix3D |
viewRowFlip()
Constructs and returns a new flip view along the row axis.
|
DoubleMatrix3D |
viewSelection(DoubleMatrix2DProcedure condition)
Constructs and returns a new selection view that is a matrix
holding all slices matching the given condition.
|
DoubleMatrix3D |
viewSelection(int[] sliceIndexes,
int[] rowIndexes,
int[] columnIndexes)
Constructs and returns a new selection view that is a matrix
holding the indicated cells.
|
DoubleMatrix2D |
viewSlice(int slice)
Constructs and returns a new 2-dimensional slice view representing
the rows and columns of the given slice.
|
DoubleMatrix3D |
viewSliceFlip()
Constructs and returns a new flip view along the slice axis.
|
DoubleMatrix3D |
viewSorted(int row,
int column)
Sorts the matrix slices into ascending order, according to the natural
ordering of the matrix values in the given [row,column]
position.
|
DoubleMatrix3D |
viewStrides(int sliceStride,
int rowStride,
int columnStride)
Constructs and returns a new stride view which is a sub matrix
consisting of every i-th cell.
|
void |
zAssign27Neighbors(DoubleMatrix3D B,
Double27Function function)
27 neighbor stencil transformation.
|
double |
zSum()
Returns the sum of all cells; Sum( x[i,j,k] ).
|
checkShape, checkShape, columns, columnStride, index, rows, rowStride, size, slices, sliceStride, toStringShort
ensureCapacity, isView, trimToSize
clone
public double aggregate(DoubleDoubleFunction aggr, DoubleFunction f)
Example:
cern.jet.math.Functions F = cern.jet.math.Functions.functions; 2 x 2 x 2 matrix 0 1 2 3 4 5 6 7 // Sum( x[slice,row,col]*x[slice,row,col] ) matrix.aggregate(F.plus,F.square); --> 140For further examples, see the package doc.
aggr
- an aggregation function taking as first argument the current
aggregation and as second argument the transformed current
cell value.f
- a function transforming the current cell value.DoubleFunctions
public double aggregate(DoubleDoubleFunction aggr, DoubleFunction f, DoubleProcedure cond)
aggr
- an aggregation function taking as first argument the current
aggregation and as second argument the transformed current
cell value.f
- a function transforming the current cell value.cond
- a condition.DoubleFunctions
public double aggregate(DoubleDoubleFunction aggr, DoubleFunction f, IntArrayList sliceList, IntArrayList rowList, IntArrayList columnList)
aggr
- an aggregation function taking as first argument the current
aggregation and as second argument the transformed current
cell value.f
- a function transforming the current cell value.sliceList
- slice indexes.rowList
- row indexes.columnList
- column indexes.DoubleFunctions
public double aggregate(DoubleMatrix3D other, DoubleDoubleFunction aggr, DoubleDoubleFunction f)
Example:
cern.jet.math.Functions F = cern.jet.math.Functions.functions; x = 2 x 2 x 2 matrix 0 1 2 3 4 5 6 7 y = 2 x 2 x 2 matrix 0 1 2 3 4 5 6 7 // Sum( x[slice,row,col] * y[slice,row,col] ) x.aggregate(y, F.plus, F.mult); --> 140 // Sum( (x[slice,row,col] + y[slice,row,col])ˆ2 ) x.aggregate(y, F.plus, F.chain(F.square,F.plus)); --> 560For further examples, see the package doc.
aggr
- an aggregation function taking as first argument the current
aggregation and as second argument the transformed current
cell values.f
- a function transforming the current cell values.IllegalArgumentException
- if
slices() != other.slices() || rows() != other.rows() || columns() != other.columns()DoubleFunctions
public DoubleMatrix3D assign(DoubleFunction function)
Example:
matrix = 1 x 2 x 2 matrix 0.5 1.5 2.5 3.5 // change each cell to its sine matrix.assign(cern.jet.math.Functions.sin); --> 1 x 2 x 2 matrix 0.479426 0.997495 0.598472 -0.350783For further examples, see the package doc.
function
- a function object taking as argument the current cell's value.DoubleFunctions
public DoubleMatrix3D assign(DoubleProcedure cond, DoubleFunction f)
cond
- a condition.f
- a function object.DoubleFunctions
public DoubleMatrix3D assign(DoubleProcedure cond, double value)
cond
- a condition.value
- a value.public DoubleMatrix3D assign(double value)
value
- the value to be filled into the cells.public DoubleMatrix3D assign(double[] values)
The values are copied. So subsequent changes in values are not reflected in the matrix, and vice-versa.
values
- the values to be filled into the cells.IllegalArgumentException
- if values.length != slices()*rows()*columns()public DoubleMatrix3D assign(double[][][] values)
The values are copied. So subsequent changes in values are not reflected in the matrix, and vice-versa.
values
- the values to be filled into the cells.IllegalArgumentException
- if
values.length != slices() || for any 0 <= slice < slices(): values[slice].length != rows()
.IllegalArgumentException
- if
for any 0 <= column < columns(): values[slice][row].length != columns()
.public DoubleMatrix3D assign(DoubleMatrix3D other)
other
- the source matrix to copy from (may be identical to the
receiver).IllegalArgumentException
- if
slices() != other.slices() || rows() != other.rows() || columns() != other.columns()public DoubleMatrix3D assign(DoubleMatrix3D y, DoubleDoubleFunction function)
Example:
// assign x[row,col] = x[row,col]<sup>y[row,col]</sup> m1 = 1 x 2 x 2 matrix 0 1 2 3 m2 = 1 x 2 x 2 matrix 0 2 4 6 m1.assign(m2, cern.jet.math.Functions.pow); --> m1 == 1 x 2 x 2 matrix 1 1 16 729For further examples, see the package doc.
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,IllegalArgumentException
- if
slices() != other.slices() || rows() != other.rows() || columns() != other.columns()DoubleFunctions
public DoubleMatrix3D assign(DoubleMatrix3D y, DoubleDoubleFunction function, IntArrayList sliceList, IntArrayList rowList, IntArrayList columnList)
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, *sliceList
- slice indexes.rowList
- row indexes.columnList
- column indexes.IllegalArgumentException
- if
slices() != other.slices() || rows() != other.rows() || columns() != other.columns()DoubleFunctions
public int cardinality()
public DoubleMatrix3D copy()
Note that the returned matrix is an independent deep copy. The returned matrix is not backed by this matrix, so changes in the returned matrix are not reflected in this matrix, and vice-versa.
public abstract Object elements()
public boolean equals(double value)
value
- the value to test against.public boolean equals(Object obj)
true
if and only if the argument is not null
and is at least a DoubleMatrix3D
object that has the same
number of slices, rows and columns as the receiver and has exactly the
same values at the same coordinates.public double get(int slice, int row, int column)
slice
- the index of the slice-coordinate.row
- the index of the row-coordinate.column
- the index of the column-coordinate.IndexOutOfBoundsException
- if
slice<0 || slice>=slices() || row<0 || row>=rows() || column<0 || column>=column()
.public double[] getMaxLocation()
public double[] getMinLocation()
public void getNegativeValues(IntArrayList sliceList, IntArrayList rowList, IntArrayList columnList, DoubleArrayList valueList)
sliceList
- the list to be filled with slice indexes, can have any size.rowList
- the list to be filled with row indexes, can have any size.columnList
- the list to be filled with column indexes, can have any size.valueList
- the list to be filled with values, can have any size.public void getNonZeros(IntArrayList sliceList, IntArrayList rowList, IntArrayList columnList, DoubleArrayList valueList)
In general, fill order is unspecified. This implementation fill
like:
for (slice = 0..slices-1) for (row = 0..rows-1) for (column = 0..colums-1) do ...
. However, subclasses are free to us any other order, even an order that
may change over time as cell values are changed. (Of course, result lists
indexes are guaranteed to correspond to the same cell). For an example,
see
getNonZeros(IntArrayList,IntArrayList,IntArrayList,DoubleArrayList)
.
sliceList
- the list to be filled with slice indexes, can have any size.rowList
- the list to be filled with row indexes, can have any size.columnList
- the list to be filled with column indexes, can have any size.valueList
- the list to be filled with values, can have any size.public void getPositiveValues(IntArrayList sliceList, IntArrayList rowList, IntArrayList columnList, DoubleArrayList valueList)
sliceList
- the list to be filled with slice indexes, can have any size.rowList
- the list to be filled with row indexes, can have any size.columnList
- the list to be filled with column indexes, can have any size.valueList
- the list to be filled with values, can have any size.public abstract double getQuick(int slice, int row, int 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): slice<0 || slice>=slices() || row<0 || row>=rows() || column<0 || column>=column().
slice
- the index of the slice-coordinate.row
- the index of the row-coordinate.column
- the index of the column-coordinate.public DoubleMatrix3D like()
public abstract DoubleMatrix3D like(int slices, int rows, int columns)
slices
- the number of slices the matrix shall have.rows
- the number of rows the matrix shall have.columns
- the number of columns the matrix shall have.public abstract DoubleMatrix2D like2D(int rows, int columns)
rows
- the number of rows the matrix shall have.columns
- the number of columns the matrix shall have.public void normalize()
public void set(int slice, int row, int column, double value)
slice
- the index of the slice-coordinate.row
- the index of the row-coordinate.column
- the index of the column-coordinate.value
- the value to be filled into the specified cell.IndexOutOfBoundsException
- if
row<0 || row>=rows() || slice<0 || slice>=slices() || column<0 || column>=column()
.public abstract void setQuick(int slice, int row, int column, double 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): slice<0 || slice>=slices() || row<0 || row>=rows() || column<0 || column>=column().
slice
- the index of the slice-coordinate.row
- the index of the row-coordinate.column
- the index of the column-coordinate.value
- the value to be filled into the specified cell.public double[][][] toArray()
The values are copied. So subsequent changes in values are not reflected in the matrix, and vice-versa.
public String toString()
toString
in class Object
DoubleFormatter
public abstract DoubleMatrix1D vectorize()
public DoubleMatrix2D viewColumn(int column)
To obtain a slice view on subranges, construct a sub-ranging view ( view().part(...)), then apply this method to the sub-range view. To obtain 1-dimensional views, apply this method, then apply another slice view (methods viewColumn, viewRow) on the intermediate 2-dimensional view. To obtain 1-dimensional views on subranges, apply both steps.
column
- the index of the column to fix.IndexOutOfBoundsException
- if column < 0 || column >= columns().viewSlice(int)
,
viewRow(int)
public DoubleMatrix3D viewColumnFlip()
viewSliceFlip()
,
viewRowFlip()
public DoubleMatrix3D viewDice(int axis0, int axis1, int axis2)
axis0
- the axis that shall become axis 0 (legal values 0..2).axis1
- the axis that shall become axis 1 (legal values 0..2).axis2
- the axis that shall become axis 2 (legal values 0..2).IllegalArgumentException
- if some of the parameters are equal or not in range 0..2.public DoubleMatrix3D viewPart(int slice, int row, int column, int depth, int height, int width)
slice
- The index of the slice-coordinate.row
- The index of the row-coordinate.column
- The index of the column-coordinate.depth
- The depth of the box.height
- The height of the box.width
- The width of the box.IndexOutOfBoundsException
- if
slice<0 || depth<0 || slice+depth>slices() || row<0 || height<0 || row+height>rows() || column<0 || width<0 || column+width>columns()public DoubleMatrix2D viewRow(int row)
To obtain a slice view on subranges, construct a sub-ranging view ( view().part(...)), then apply this method to the sub-range view. To obtain 1-dimensional views, apply this method, then apply another slice view (methods viewColumn, viewRow) on the intermediate 2-dimensional view. To obtain 1-dimensional views on subranges, apply both steps.
row
- the index of the row to fix.IndexOutOfBoundsException
- if row < 0 || row >= row().viewSlice(int)
,
viewColumn(int)
public DoubleMatrix3D viewRowFlip()
viewSliceFlip()
,
viewColumnFlip()
public DoubleMatrix3D viewSelection(DoubleMatrix2DProcedure condition)
Example:
// extract and view all slices which have an aggregate sum > 1000 matrix.viewSelection(new DoubleMatrix2DProcedure() { public final boolean apply(DoubleMatrix2D m) { return m.zSum > 1000; } });For further examples, see the package doc. The returned view is backed by this matrix, so changes in the returned view are reflected in this matrix, and vice-versa.
condition
- The condition to be matched.public DoubleMatrix3D viewSelection(int[] sliceIndexes, int[] rowIndexes, int[] columnIndexes)
DoubleMatrix2D.viewSelection(int[],int[])
.
Note that modifying the index arguments after this call has returned has no effect on the view. The returned view is backed by this matrix, so changes in the returned view are reflected in this matrix, and vice-versa.
sliceIndexes
- The slices of the cells that shall be visible in the new view.
To indicate that all slices shall be visible, simply
set this parameter to null.rowIndexes
- The rows of the cells that shall be visible in the new view.
To indicate that all rows shall be visible, simply set
this parameter to null.columnIndexes
- The columns of the cells that shall be visible in the new
view. To indicate that all columns shall be visible,
simply set this parameter to null.IndexOutOfBoundsException
- if !(0 <= sliceIndexes[i] < slices()) for any
i=0..sliceIndexes.length()-1.IndexOutOfBoundsException
- if !(0 <= rowIndexes[i] < rows()) for any
i=0..rowIndexes.length()-1.IndexOutOfBoundsException
- if !(0 <= columnIndexes[i] < columns()) for any
i=0..columnIndexes.length()-1.public DoubleMatrix2D viewSlice(int slice)
To obtain a slice view on subranges, construct a sub-ranging view ( view().part(...)), then apply this method to the sub-range view. To obtain 1-dimensional views, apply this method, then apply another slice view (methods viewColumn, viewRow) on the intermediate 2-dimensional view. To obtain 1-dimensional views on subranges, apply both steps.
slice
- the index of the slice to fix.IndexOutOfBoundsException
- if slice < 0 || slice >= slices().viewRow(int)
,
viewColumn(int)
public DoubleMatrix3D viewSliceFlip()
viewRowFlip()
,
viewColumnFlip()
public DoubleMatrix3D viewSorted(int row, int column)
DoubleSorting.sort(DoubleMatrix3D,int,int)
. For more advanced sorting functionality, see
DoubleSorting
.IndexOutOfBoundsException
- if
row < 0 || row >= rows() || column < 0 || column >= columns()
.public DoubleMatrix3D viewStrides(int sliceStride, int rowStride, int columnStride)
sliceStride
- the slice step factor.rowStride
- the row step factor.columnStride
- the column step factor.IndexOutOfBoundsException
- if sliceStride<=0 || rowStride<=0 || columnStride<=0
.public void zAssign27Neighbors(DoubleMatrix3D B, Double27Function function)
B[k,i,j] = function.apply( A[k-1,i-1,j-1], A[k-1,i-1,j], A[k-1,i-1,j+1], A[k-1,i, j-1], A[k-1,i, j], A[k-1,i, j+1], A[k-1,i+1,j-1], A[k-1,i+1,j], A[k-1,i+1,j+1], A[k ,i-1,j-1], A[k ,i-1,j], A[k ,i-1,j+1], A[k ,i, j-1], A[k ,i, j], A[k ,i, j+1], A[k ,i+1,j-1], A[k ,i+1,j], A[k ,i+1,j+1], A[k+1,i-1,j-1], A[k+1,i-1,j], A[k+1,i-1,j+1], A[k+1,i, j-1], A[k+1,i, j], A[k+1,i, j+1], A[k+1,i+1,j-1], A[k+1,i+1,j], A[k+1,i+1,j+1] ) x x x - - x x x - - - - x o x - - x o x - - - - x x x - - x x x ... - x x x - - - - - - - - - x o x - - - - - - - - - x x xMake sure that cells of this and B do not overlap. In case of overlapping views, behaviour is unspecified.
Example:
final double alpha = 0.25; final double beta = 0.75; cern.colt.function.Double27Function f = new cern.colt.function.Double27Function() { public final double apply( double a000, double a001, double a002, double a010, double a011, double a012, double a020, double a021, double a022, double a100, double a101, double a102, double a110, double a111, double a112, double a120, double a121, double a122, double a200, double a201, double a202, double a210, double a211, double a212, double a220, double a221, double a222) { return beta*a111 + alpha*(a000 + ... + a222); } }; A.zAssign27Neighbors(B,f);
B
- the matrix to hold the results.function
- the function to be applied to the 27 cells.NullPointerException
- if function==null.IllegalArgumentException
- if
rows() != B.rows() || columns() != B.columns() || slices() != B.slices()
.public double zSum()
Jump to the Parallel Colt Homepage