ARTS
2.4.0(git:4fb77825)
|
The MatrixView class. More...
#include <matpackI.h>
Public Types | |
typedef Iterator2D | iterator |
Public Types inherited from ConstMatrixView | |
typedef ConstIterator2D | const_iterator |
Public Member Functions | |
constexpr | MatrixView (const MatrixView &)=default |
Numeric & | operator() (Index r, Index c) |
Plain index operator. More... | |
Numeric & | get (Index r, Index c) |
Get element implementation without assertions. More... | |
MatrixView | operator() (const Range &r, const Range &c) |
Index operator for subrange. More... | |
VectorView | operator() (const Range &r, Index c) |
Index operator returning a column as an object of type VectorView. More... | |
VectorView | operator() (Index r, const Range &c) |
Index operator returning a row as an object of type VectorView. More... | |
Iterator2D | begin () |
** Return const iterator to first row. Has to be redefined here, since it is More... | |
Iterator2D | end () |
Return iterator behind last row. More... | |
MatrixView & | operator= (const ConstMatrixView &v) |
Assignment operator. More... | |
MatrixView & | operator= (const MatrixView &v) |
Assignment from MatrixView to MatrixView. More... | |
MatrixView & | operator= (const Matrix &v) |
Assignment from a Matrix. More... | |
MatrixView & | operator= (const ConstVectorView &v) |
Assignment from a vector. More... | |
MatrixView & | operator= (Numeric x) |
Assigning a scalar to a MatrixView will set all elements to this value. More... | |
MatrixView & | operator*= (Numeric x) |
Multiplication by scalar. More... | |
MatrixView & | operator/= (Numeric x) |
Division by scalar. More... | |
MatrixView & | operator+= (Numeric x) |
Addition of scalar. More... | |
MatrixView & | operator-= (Numeric x) |
Subtraction of scalar. More... | |
MatrixView & | operator*= (const ConstMatrixView &x) |
Element-vise multiplication by another Matrix. More... | |
MatrixView & | operator/= (const ConstMatrixView &x) |
Element-vise division by another Matrix. More... | |
MatrixView & | operator+= (const ConstMatrixView &x) |
Element-vise addition of another Matrix. More... | |
MatrixView & | operator-= (const ConstMatrixView &x) |
Element-vise subtraction of another Matrix. More... | |
MatrixView & | operator*= (const ConstVectorView &x) |
Element-vise multiplication by a Vector (acting like a 1-column Matrix). More... | |
MatrixView & | operator/= (const ConstVectorView &x) |
Element-vise division by a Vector (acting like a 1-column Matrix). More... | |
MatrixView & | operator+= (const ConstVectorView &x) |
Element-vise addition of a Vector (acting like a 1-column Matrix). More... | |
MatrixView & | operator-= (const ConstVectorView &x) |
Element-vise subtraction of a Vector (acting like a 1-column Matrix). More... | |
const Numeric * | get_c_array () const |
Conversion to plain C-array. More... | |
Numeric * | get_c_array () |
Conversion to plain C-array. More... | |
virtual | ~MatrixView ()=default |
Destructor. More... | |
ConstIterator2D | begin () const |
Return const iterator to first row. More... | |
ConstIterator2D | end () const |
Return const iterator behind last row. More... | |
Numeric | get (Index r, Index c) const |
Get element implementation without assertions. More... | |
Public Member Functions inherited from ConstMatrixView | |
constexpr | ConstMatrixView (const ConstMatrixView &)=default |
constexpr | ConstMatrixView (ConstMatrixView &&)=default |
ConstMatrixView & | operator= (const ConstMatrixView &)=default |
ConstMatrixView & | operator= (ConstMatrixView &&)=default |
bool | empty () const |
Returns true if variable size is zero. More... | |
Index | nrows () const |
Returns the number of rows. More... | |
Index | ncols () const |
Returns the number of columns. More... | |
Numeric | operator() (Index r, Index c) const |
Plain const index operator. More... | |
Numeric | get (Index r, Index c) const |
Get element implementation without assertions. More... | |
ConstMatrixView | operator() (const Range &r, const Range &c) const |
Const index operator for subrange. More... | |
ConstVectorView | operator() (const Range &r, Index c) const |
Const index operator returning a column as an object of type ConstVectorView. More... | |
ConstVectorView | operator() (Index r, const Range &c) const |
Const index operator returning a row as an object of type ConstVectorView. More... | |
ConstIterator2D | begin () const |
Return const iterator to first row. More... | |
ConstIterator2D | end () const |
Return const iterator behind last row. More... | |
ConstVectorView | diagonal () const |
Matrix diagonal as vector. More... | |
virtual | ~ConstMatrixView ()=default |
Destructor. More... | |
Protected Member Functions | |
MatrixView ()=default | |
MatrixView (Numeric *data, const Range &r, const Range &c) | |
Explicit constructor. More... | |
MatrixView (Numeric *data, const Range &pr, const Range &pc, const Range &nr, const Range &nc) | |
Recursive constructor. More... | |
Protected Member Functions inherited from ConstMatrixView | |
ConstMatrixView ()=default | |
ConstMatrixView (Numeric *data, const Range &r, const Range &c) | |
Explicit constructor. More... | |
ConstMatrixView (Numeric *data, const Range &pr, const Range &pc, const Range &nr, const Range &nc) | |
Recursive constructor. More... | |
Friends | |
class | VectorView |
class | Iterator3D |
class | Tensor3View |
class | Tensor4View |
class | Tensor5View |
class | Tensor6View |
class | Tensor7View |
class | ComplexMatrixView |
ConstMatrixView | transpose (ConstMatrixView m) |
Const version of transpose. More... | |
MatrixView | transpose (MatrixView m) |
Returns the transpose. More... | |
void | mult (MatrixView, const ConstMatrixView &, const ConstMatrixView &) |
Matrix-Matrix Multiplication. More... | |
Additional Inherited Members | |
Protected Attributes inherited from ConstMatrixView | |
Range | mrr {0, 0, 1} |
The row range of mdata that is actually used. More... | |
Range | mcr {0, 0, 1} |
The column range of mdata that is actually used. More... | |
Numeric * | mdata {nullptr} |
Pointer to the plain C array that holds the data. More... | |
The MatrixView class.
This contains the main implementation of a Matrix. It defines the concepts of MatrixView. Plus additionally the recursive subrange operator, which makes it possible to create a MatrixView from a subrange of a MatrixView.
The class Matrix is just a special case of a MatrixView which also allocates storage.
Definition at line 1093 of file matpackI.h.
typedef Iterator2D MatrixView::iterator |
Definition at line 1104 of file matpackI.h.
|
constexprdefault |
|
virtualdefault |
Destructor.
|
protecteddefault |
Explicit constructor.
This one is used by Matrix to initialize its own MatrixView part. The row range rr must have a stride to account for the length of one row.
Definition at line 882 of file matpackI.cc.
|
protected |
Recursive constructor.
This is used to construct SubMatrices from SubMatrices. That means that the new ranges have to be interpreted relative to the original ranges.
The new ranges may contain -1 for the extent which acts as a joker. However, the used Range constructor converts this to an explicit range, consistent with the original Range.
Definition at line 901 of file matpackI.cc.
Iterator2D MatrixView::begin | ( | ) |
** Return const iterator to first row. Has to be redefined here, since it is
** Return const iterator behind last row. */ Return iterator to first row.
Definition at line 618 of file matpackI.cc.
References ConstMatrixView::mcr, ConstMatrixView::mdata, ConstMatrixView::mrr, Range::mstart, Range::mstride, and VectorView.
Referenced by copy(), Matrix::Matrix(), mult_general(), operator*=(), operator+=(), operator-=(), operator/=(), and operator=().
ConstIterator2D ConstMatrixView::begin |
Return const iterator to first row.
Definition at line 1019 of file matpackI.cc.
Iterator2D MatrixView::end | ( | ) |
Return iterator behind last row.
Definition at line 623 of file matpackI.cc.
References ConstMatrixView::mcr, ConstMatrixView::mdata, Range::mextent, ConstMatrixView::mrr, Range::mstart, Range::mstride, and VectorView.
Referenced by copy(), mult_general(), operator*=(), operator+=(), operator-=(), operator/=(), and operator=().
ConstIterator2D ConstMatrixView::end |
Return const iterator behind last row.
Definition at line 1020 of file matpackI.cc.
Get element implementation without assertions.
Definition at line 1118 of file matpackI.h.
|
inline |
Get element implementation without assertions.
Definition at line 1009 of file matpackI.h.
Numeric * MatrixView::get_c_array | ( | ) |
Conversion to plain C-array.
This function returns a pointer to the raw data. It fails if the MatrixView is not pointing to the beginning of a Matrix or the stride is not 1 because the caller expects to get a C array with continuous data.
Definition at line 750 of file matpackI.cc.
References ConstMatrixView::mcr, ConstMatrixView::mdata, Range::mextent, ConstMatrixView::mrr, Range::mstart, and Range::mstride.
const Numeric * MatrixView::get_c_array | ( | ) | const |
Conversion to plain C-array.
This function returns a pointer to the raw data. It fails if the MatrixView is not pointing to the beginning of a Matrix or the stride is not 1 because the caller expects to get a C array with continuous data.
Definition at line 735 of file matpackI.cc.
References ConstMatrixView::mcr, ConstMatrixView::mdata, Range::mextent, ConstMatrixView::mrr, Range::mstart, and Range::mstride.
Referenced by get_covariance_matrix_block(), get_variable_value(), nca_get_data_Matrix(), nca_put_var_Matrix(), nca_read_from_file(), nca_write_to_file(), Reduce(), InteractiveWorkspace::set_matrix_variable(), and test38().
MatrixView MatrixView::operator() | ( | const Range & | r, |
const Range & | c | ||
) |
Index operator for subrange.
We have to also account for the case, that *this is already a subrange of a Matrix. This allows correct recursive behavior.
Definition at line 576 of file matpackI.cc.
VectorView MatrixView::operator() | ( | const Range & | r, |
Index | c | ||
) |
Index operator returning a column as an object of type VectorView.
r | A range of rows. |
c | Index of selected column |
Definition at line 584 of file matpackI.cc.
VectorView MatrixView::operator() | ( | Index | r, |
const Range & | c | ||
) |
Index operator returning a row as an object of type VectorView.
r | Index of selected row. |
c | Range of columns |
Definition at line 596 of file matpackI.cc.
Plain index operator.
Definition at line 1108 of file matpackI.h.
MatrixView & MatrixView::operator*= | ( | const ConstMatrixView & | x | ) |
Element-vise multiplication by another Matrix.
Definition at line 760 of file matpackI.cc.
References begin(), ConstVectorView::begin(), end(), ConstMatrixView::ncols(), ConstMatrixView::nrows(), and ARTS::Var::x().
MatrixView & MatrixView::operator*= | ( | const ConstVectorView & | x | ) |
Element-vise multiplication by a Vector (acting like a 1-column Matrix).
Definition at line 824 of file matpackI.cc.
References begin(), end(), ConstMatrixView::ncols(), ConstMatrixView::nrows(), and ARTS::Var::x().
MatrixView & MatrixView::operator*= | ( | Numeric | x | ) |
Multiplication by scalar.
Definition at line 690 of file matpackI.cc.
MatrixView & MatrixView::operator+= | ( | const ConstMatrixView & | x | ) |
Element-vise addition of another Matrix.
Definition at line 792 of file matpackI.cc.
References begin(), ConstVectorView::begin(), end(), ConstMatrixView::ncols(), ConstMatrixView::nrows(), and ARTS::Var::x().
MatrixView & MatrixView::operator+= | ( | const ConstVectorView & | x | ) |
Element-vise addition of a Vector (acting like a 1-column Matrix).
Definition at line 852 of file matpackI.cc.
References begin(), end(), ConstMatrixView::ncols(), ConstMatrixView::nrows(), and ARTS::Var::x().
MatrixView & MatrixView::operator+= | ( | Numeric | x | ) |
MatrixView & MatrixView::operator-= | ( | const ConstMatrixView & | x | ) |
Element-vise subtraction of another Matrix.
Definition at line 808 of file matpackI.cc.
References begin(), ConstVectorView::begin(), end(), ConstMatrixView::ncols(), ConstMatrixView::nrows(), and ARTS::Var::x().
MatrixView & MatrixView::operator-= | ( | const ConstVectorView & | x | ) |
Element-vise subtraction of a Vector (acting like a 1-column Matrix).
Definition at line 866 of file matpackI.cc.
References begin(), end(), ConstMatrixView::ncols(), ConstMatrixView::nrows(), and ARTS::Var::x().
MatrixView & MatrixView::operator-= | ( | Numeric | x | ) |
MatrixView & MatrixView::operator/= | ( | const ConstMatrixView & | x | ) |
Element-vise division by another Matrix.
Definition at line 776 of file matpackI.cc.
References begin(), ConstVectorView::begin(), end(), ConstMatrixView::ncols(), ConstMatrixView::nrows(), and ARTS::Var::x().
MatrixView & MatrixView::operator/= | ( | const ConstVectorView & | x | ) |
Element-vise division by a Vector (acting like a 1-column Matrix).
Definition at line 838 of file matpackI.cc.
References begin(), end(), ConstMatrixView::ncols(), ConstMatrixView::nrows(), and ARTS::Var::x().
MatrixView & MatrixView::operator/= | ( | Numeric | x | ) |
MatrixView & MatrixView::operator= | ( | const ConstMatrixView & | m | ) |
Assignment operator.
This copies the data from another MatrixView to this MatrixView. Dimensions must agree! Resizing would destroy the selection that we might have done in this MatrixView by setting its range.
Definition at line 633 of file matpackI.cc.
References begin(), ConstMatrixView::begin(), copy(), ConstMatrixView::end(), ConstMatrixView::mcr, Range::mextent, and ConstMatrixView::mrr.
MatrixView & MatrixView::operator= | ( | const ConstVectorView & | v | ) |
Assignment from a vector.
This copies the data from a VectorView to this MatrixView. Dimensions must agree! Resizing would destroy the selection that we might have done in this MatrixView by setting its range.
Definition at line 672 of file matpackI.cc.
References begin(), ConstMatrixView::begin(), copy(), ConstMatrixView::end(), ConstMatrixView::mcr, Range::mextent, ConstMatrixView::mrr, and ConstVectorView::nelem().
MatrixView & MatrixView::operator= | ( | const Matrix & | m | ) |
Assignment from a Matrix.
This must exist to overide the automatically generated assignment operators, which don't copy the contents!
Definition at line 659 of file matpackI.cc.
References begin(), copy(), end(), ConstMatrixView::mcr, Range::mextent, and ConstMatrixView::mrr.
MatrixView & MatrixView::operator= | ( | const MatrixView & | m | ) |
Assignment from MatrixView to MatrixView.
This is a tricky one. The problem is that since MatrixView is derived from ConstMatrixView, a default = operator is generated by the compiler, which does not do what we want. So we need this one to override the default.
Definition at line 647 of file matpackI.cc.
References begin(), copy(), end(), ConstMatrixView::mcr, Range::mextent, and ConstMatrixView::mrr.
MatrixView & MatrixView::operator= | ( | Numeric | x | ) |
Assigning a scalar to a MatrixView will set all elements to this value.
Definition at line 684 of file matpackI.cc.
References begin(), copy(), end(), and ARTS::Var::x().
|
friend |
Definition at line 1169 of file matpackI.h.
|
friend |
Definition at line 1163 of file matpackI.h.
|
friend |
Matrix-Matrix Multiplication.
Performs the matrix multiplication A = B * C. The dimensions must match, i.e. A must be a m times n matrix, B a m times k matrix and C a k times c matrix. No memory reallocation takes place, only the data is copied. Using this function on overlapping MatrixViews belonging to the same Matrix will lead to unpredictable results. In particular, this means that A and B must not be the same matrix!
If the memory layout allows it, the multiplication is performed using BLAS' _dgemm, which leads to a significant speed up of the operation. To be compatible with BLAS the matrix views A, B and C must satisfy:
That means that A and B can be ConstMatrixView objects corresponding to transposed/non-transposed submatrices of a matrix, that are continuous along their first/second dimension. C must correspond to a non-transposed submatrix of a matrix, that is continuous along its second dimension.
[in,out] | A | The matrix A, that will hold the result of the multiplication. |
[in] | B | The matrix B |
[in] | C | The matrix C |
Definition at line 1242 of file matpackI.cc.
|
friend |
Definition at line 1164 of file matpackI.h.
|
friend |
Definition at line 1165 of file matpackI.h.
|
friend |
Definition at line 1166 of file matpackI.h.
|
friend |
Definition at line 1167 of file matpackI.h.
|
friend |
Definition at line 1168 of file matpackI.h.
|
friend |
Const version of transpose.
Definition at line 1444 of file matpackI.cc.
|
friend |
Returns the transpose.
This creates a special MatrixView for the transpose. The original is not changed!
Definition at line 1450 of file matpackI.cc.
|
friend |
Definition at line 1162 of file matpackI.h.