35 #ifndef _BLAZE_MATH_EXPRESSIONS_SMATTRANSPOSER_H_
36 #define _BLAZE_MATH_EXPRESSIONS_SMATTRANSPOSER_H_
71 class SMatTransposer :
public SparseMatrix< SMatTransposer<MT,SO>, SO >
93 enum { smpAssignable = MT::smpAssignable };
148 return sm_.cbegin(i);
164 return sm_.cbegin(i);
223 template<
typename Other >
240 template<
typename Other >
256 return sm_.columns();
294 return sm_.insert( j, i, value );
309 sm_.reserve( nonzeros );
327 inline void reserve(
size_t i,
size_t nonzeros ) {
328 sm_.reserve( i, nonzeros );
360 sm_.append( j, i, value, check );
388 template<
typename Other >
391 return sm_.canAlias( alias );
401 template<
typename Other >
404 return sm_.isAliased( alias );
415 return sm_.canSMPAssign();
430 template<
typename MT2 >
441 const size_t m(
rows() );
443 for(
size_t i=0UL; i<m; ++i ) {
444 for( RhsIterator element=(~rhs).
begin(i); element!=(~rhs).
end(i); ++element )
445 sm_.append( element->index(), i, element->value() );
462 template<
typename MT2 >
473 const size_t m(
rows() );
477 std::vector<size_t> rowLengths( m, 0UL );
478 for(
size_t j=0UL; j<n; ++j ) {
479 for( RhsIterator element=(~rhs).
begin(j); element!=(~rhs).
end(j); ++element )
480 ++rowLengths[element->index()];
484 for(
size_t i=0UL; i<m; ++i ) {
485 sm_.reserve( i, rowLengths[i] );
489 for(
size_t j=0UL; j<n; ++j ) {
490 for( RhsIterator element=(~rhs).
begin(j); element!=(~rhs).
end(j); ++element ) {
491 sm_.append( j, element->index(), element->value() );
529 template<
typename MT >
552 enum { smpAssignable = MT::smpAssignable };
597 return sm_.cbegin(j);
608 return sm_.cbegin(j);
652 template<
typename Other >
669 template<
typename Other >
684 inline size_t rows()
const {
685 return sm_.columns();
694 inline size_t columns()
const {
704 inline void reset() {
723 return sm_.insert( j, i, value );
737 inline void reserve(
size_t nonzeros ) {
738 sm_.reserve( nonzeros );
753 inline void reserve(
size_t i,
size_t nonzeros ) {
754 sm_.reserve( i, nonzeros );
786 sm_.append( j, i, value, check );
814 template<
typename Other >
815 inline bool canAlias(
const Other* alias )
const
817 return sm_.canAlias( alias );
827 template<
typename Other >
828 inline bool isAliased(
const Other* alias )
const
830 return sm_.isAliased( alias );
841 return sm_.canSMPAssign();
856 template<
typename MT2 >
857 inline void assign(
const SparseMatrix<MT2,false>& rhs )
867 const size_t m(
rows() );
871 std::vector<size_t> columnLengths( n, 0UL );
872 for(
size_t i=0UL; i<m; ++i ) {
873 for( RhsIterator element=(~rhs).
begin(i); element!=(~rhs).
end(i); ++element )
874 ++columnLengths[element->index()];
878 for(
size_t j=0UL; j<n; ++j ) {
879 sm_.reserve( j, columnLengths[j] );
883 for(
size_t i=0UL; i<m; ++i ) {
884 for( RhsIterator element=(~rhs).
begin(i); element!=(~rhs).
end(i); ++element ) {
885 sm_.append( element->index(), i, element->value() );
902 template<
typename MT2 >
903 inline void assign(
const SparseMatrix<MT2,true>& rhs )
915 for(
size_t j=0UL; j<n; ++j ) {
916 for( RhsIterator element=(~rhs).
begin(j); element!=(~rhs).
end(j); ++element )
917 sm_.append( j, element->index(), element->value() );
956 template<
typename MT
958 inline void reset( SMatTransposer<MT,SO>& m )
976 template<
typename MT,
bool SO >
977 struct SubmatrixTrait< SMatTransposer<MT,SO> >
MT OppositeType
Result type with opposite storage order for expression template evaluations.
Definition: SMatTransposer.h:77
MT::ConstReference ConstReference
Reference to a constant matrix value.
Definition: SMatTransposer.h:83
MT::ElementType ElementType
Resulting element type.
Definition: SMatTransposer.h:79
Constraint on the data type.
ConstIterator begin(size_t i) const
Returns an iterator to the first non-zero element of row/column i.
Definition: SMatTransposer.h:147
#define BLAZE_USER_ASSERT(expr, msg)
Run time assertion macro for user checks.In case of an invalid run time expression, the program execution is terminated. The BLAZE_USER_ASSERT macro can be disabled by setting the BLAZE_USER_ASSERT flag to zero or by defining NDEBUG during the compilation.
Definition: Assert.h:117
Iterator begin(size_t i)
Returns an iterator to the first non-zero element of row/column i.
Definition: SMatTransposer.h:131
Efficient implementation of a compressed matrix.The CompressedMatrix class template is the represent...
Definition: CompressedMatrix.h:205
#define BLAZE_CONSTRAINT_MUST_NOT_BE_COMPUTATION_TYPE(T)
Constraint on the data type.In case the given data type T is a computational expression (i...
Definition: Computation.h:118
MT::ReturnType ReturnType
Return type for expression template evaluations.
Definition: SMatTransposer.h:80
void reset()
Resets the matrix elements.
Definition: SMatTransposer.h:275
SMatTransposer< MT, SO > This
Type of this SMatTransposer instance.
Definition: SMatTransposer.h:75
bool canAlias(const Other *alias) const
Returns whether the matrix can alias with the given address alias.
Definition: SMatTransposer.h:389
BLAZE_ALWAYS_INLINE size_t nonZeros(const Matrix< MT, SO > &matrix)
Returns the total number of non-zero elements in the matrix.
Definition: Matrix.h:386
MT::ResultType TransposeType
Transpose type for expression template evaluations.
Definition: SMatTransposer.h:78
Iterator end(size_t i)
Returns an iterator just past the last non-zero element of row/column i.
Definition: SMatTransposer.h:179
EnableIf< IsNumeric< Other >, SMatTransposer >::Type & operator*=(Other rhs)
Multiplication assignment operator for the multiplication between a matrix and a scalar value ( )...
Definition: SMatTransposer.h:224
Iterator insert(size_t i, size_t j, const ElementType &value)
Inserting an element into the sparse matrix.
Definition: SMatTransposer.h:293
ConstIterator cend(size_t i) const
Returns an iterator just past the last non-zero element of row/column i.
Definition: SMatTransposer.h:211
Base class for sparse matrices.The SparseMatrix class is a base class for all sparse matrix classes...
Definition: Forward.h:107
Header file for the SparseMatrix base class.
Constraint on the data type.
MT::Iterator Iterator
Iterator over non-constant elements.
Definition: SMatTransposer.h:84
bool isAliased(const Other *alias) const
Returns whether the matrix is aliased with the given address alias.
Definition: SMatTransposer.h:402
void reserve(size_t i, size_t nonzeros)
Setting the minimum capacity of a specific row/column of the sparse matrix.
Definition: SMatTransposer.h:327
#define BLAZE_CONSTRAINT_MUST_BE_COLUMN_MAJOR_MATRIX_TYPE(T)
Constraint on the data type.In case the given data type T is not a column-major dense or sparse matri...
Definition: StorageOrder.h:161
const Element * ConstIterator
Iterator over constant elements.
Definition: CompressedMatrix.h:2482
void assign(const SparseMatrix< MT2, false > &rhs)
Implementation of the transpose assignment of a row-major sparse matrix.
Definition: SMatTransposer.h:431
MT::TransposeType ResultType
Result type for expression template evaluations.
Definition: SMatTransposer.h:76
ConstIterator end(size_t i) const
Returns an iterator just past the last non-zero element of row/column i.
Definition: SMatTransposer.h:195
CompressedMatrix< Type, false > TransposeType
Transpose type for expression template evaluations.
Definition: CompressedMatrix.h:2475
Constraints on the storage order of matrix types.
Type ElementType
Type of the sparse matrix elements.
Definition: CompressedMatrix.h:2476
void append(size_t i, size_t j, const ElementType &value, bool check=false)
Appending an element to the specified row/column of the sparse matrix.
Definition: SMatTransposer.h:359
const Type & ConstReference
Reference to a constant matrix value.
Definition: CompressedMatrix.h:2480
size_t rows() const
Returns the current number of rows of the matrix.
Definition: SMatTransposer.h:255
Header file for the EnableIf class template.
Header file for the IsNumeric type trait.
MT::ConstIterator ConstIterator
Iterator over constant elements.
Definition: SMatTransposer.h:85
void reserve(size_t nonzeros)
Setting the minimum capacity of the sparse matrix.
Definition: SMatTransposer.h:308
Expression object for the transposition of a sparse matrix.The SMatTransposer class is a wrapper obje...
Definition: Forward.h:103
#define BLAZE_CONSTRAINT_MUST_BE_ROW_MAJOR_MATRIX_TYPE(T)
Constraint on the data type.In case the given data type T is not a row-major dense or sparse matrix t...
Definition: StorageOrder.h:81
const Type & ReturnType
Return type for expression template evaluations.
Definition: CompressedMatrix.h:2477
Header file for run time assertion macros.
Header file for the submatrix trait.
SMatTransposer(MT &sm)
Constructor for the SMatTransposer class.
Definition: SMatTransposer.h:101
Substitution Failure Is Not An Error (SFINAE) class.The EnableIf class template is an auxiliary tool ...
Definition: EnableIf.h:184
MT & sm_
The sparse matrix operand.
Definition: SMatTransposer.h:499
Element * Iterator
Iterator over non-constant elements.
Definition: CompressedMatrix.h:2481
ConstReference operator()(size_t i, size_t j) const
2D-access to the matrix elements.
Definition: SMatTransposer.h:113
BLAZE_ALWAYS_INLINE void reset(const NonNumericProxy< MT > &proxy)
Resetting the represented element to the default initial values.
Definition: NonNumericProxy.h:833
ConstIterator cbegin(size_t i) const
Returns an iterator to the first non-zero element of row/column i.
Definition: SMatTransposer.h:163
const This & CompositeType
Data type for composite expression templates.
Definition: SMatTransposer.h:81
void assign(const SparseMatrix< MT2, true > &rhs)
Implementation of the transpose assignment of a column-major sparse matrix.
Definition: SMatTransposer.h:463
bool canSMPAssign() const
Returns whether the matrix can be used in SMP assignments.
Definition: SMatTransposer.h:413
void finalize(size_t i)
Finalizing the element insertion of a row/column.
Definition: SMatTransposer.h:377
EnableIf< IsNumeric< Other >, SMatTransposer >::Type & operator/=(Other rhs)
Division assignment operator for the division of a matrix by a scalar value ( ).
Definition: SMatTransposer.h:241
This ResultType
Result type for expression template evaluations.
Definition: CompressedMatrix.h:2473
Header file for basic type definitions.
MatrixAccessProxy< This > Reference
Reference to a non-constant matrix value.
Definition: CompressedMatrix.h:2479
MT::Reference Reference
Reference to a non-constant matrix value.
Definition: SMatTransposer.h:82
#define BLAZE_INTERNAL_ASSERT(expr, msg)
Run time assertion macro for internal checks.In case of an invalid run time expression, the program execution is terminated. The BLAZE_INTERNAL_ASSERT macro can be disabled by setting the BLAZE_USER_ASSERTION flag to zero or by defining NDEBUG during the compilation.
Definition: Assert.h:101
#define BLAZE_CONSTRAINT_MUST_BE_SPARSE_MATRIX_TYPE(T)
Constraint on the data type.In case the given data type T is not a sparse, N-dimensional matrix type...
Definition: SparseMatrix.h:79
size_t columns() const
Returns the current number of columns of the matrix.
Definition: SMatTransposer.h:265