35 #ifndef _BLAZE_MATH_CUSTOMMATRIX_H_ 36 #define _BLAZE_MATH_CUSTOMMATRIX_H_ 72 template<
typename Type
76 class Rand< CustomMatrix<Type,AF,PF,SO> >
82 inline void randomize( CustomMatrix<Type,AF,PF,SO>& matrix )
const;
84 template<
typename Arg >
85 inline void randomize( CustomMatrix<Type,AF,PF,SO>& matrix,
const Arg&
min,
const Arg&
max )
const;
100 template<
typename Type
104 inline void Rand< CustomMatrix<Type,AF,PF,SO> >
::randomize( CustomMatrix<Type,AF,PF,SO>& matrix )
const 108 const size_t m( matrix.rows() );
109 const size_t n( matrix.columns() );
111 for(
size_t i=0UL; i<m; ++i ) {
112 for(
size_t j=0UL; j<n; ++j ) {
130 template<
typename Type
134 template<
typename Arg >
135 inline void Rand< CustomMatrix<Type,AF,PF,SO> >
::randomize( CustomMatrix<Type,AF,PF,SO>& matrix,
136 const Arg&
min,
const Arg&
max )
const 140 const size_t m( matrix.rows() );
141 const size_t n( matrix.columns() );
143 for(
size_t i=0UL; i<m; ++i ) {
144 for(
size_t j=0UL; j<n; ++j ) {
169 template<
typename Type
173 void makeSymmetric( CustomMatrix<Type,AF,PF,SO>& matrix )
181 const size_t n( matrix.rows() );
183 for(
size_t i=0UL; i<n; ++i ) {
184 for(
size_t j=0UL; j<i; ++j ) {
186 matrix(j,i) = matrix(i,j);
207 template<
typename Type
212 void makeSymmetric( CustomMatrix<Type,AF,PF,SO>& matrix,
const Arg& min,
const Arg& max )
220 const size_t n( matrix.rows() );
222 for(
size_t i=0UL; i<n; ++i ) {
223 for(
size_t j=0UL; j<i; ++j ) {
225 matrix(j,i) = matrix(i,j);
244 template<
typename Type
248 void makeHermitian( CustomMatrix<Type,AF,PF,SO>& matrix )
254 using BT = UnderlyingBuiltin_<Type>;
260 const size_t n( matrix.rows() );
262 for(
size_t i=0UL; i<n; ++i ) {
263 for(
size_t j=0UL; j<i; ++j ) {
265 matrix(j,i) =
conj( matrix(i,j) );
267 matrix(i,i) = rand<BT>();
286 template<
typename Type
291 void makeHermitian( CustomMatrix<Type,AF,PF,SO>& matrix,
const Arg& min,
const Arg& max )
297 using BT = UnderlyingBuiltin_<Type>;
303 const size_t n( matrix.rows() );
305 for(
size_t i=0UL; i<n; ++i ) {
306 for(
size_t j=0UL; j<i; ++j ) {
308 matrix(j,i) =
conj( matrix(i,j) );
310 matrix(i,i) = rand<BT>(
real( min ),
real( max ) );
327 template<
typename Type
331 void makePositiveDefinite( CustomMatrix<Type,AF,PF,SO>& matrix )
341 const size_t n( matrix.rows() );
344 matrix *=
ctrans( matrix );
346 for(
size_t i=0UL; i<n; ++i ) {
347 matrix(i,i) += Type(n);
#define BLAZE_THROW_INVALID_ARGUMENT(MESSAGE)
Macro for the emission of a std::invalid_argument exception.This macro encapsulates the default way o...
Definition: Exception.h:235
Constraint on the data type.
decltype(auto) real(const DenseMatrix< MT, SO > &dm)
Returns a matrix containing the real part of each single element of dm.
Definition: DMatMapExpr.h:1387
void randomize(T &&value)
Randomization of a given variable.
Definition: Random.h:929
const ElementType_< MT > min(const DenseMatrix< MT, SO > &dm)
Returns the smallest element of the dense matrix.
Definition: DenseMatrix.h:1903
Implementation of a random number generator.
const ElementType_< MT > max(const DenseMatrix< MT, SO > &dm)
Returns the largest element of the dense matrix.
Definition: DenseMatrix.h:1950
Header file for the implementation of a customizable matrix.
void randomize(T &value) const
Randomization of the given variable with a new value in the range .
Definition: Random.h:292
Namespace of the Blaze C++ math library.
Definition: Blaze.h:58
decltype(auto) ctrans(const DenseMatrix< MT, SO > &dm)
Returns the conjugate transpose matrix of dm.
Definition: DMatMapExpr.h:1359
Header file for the UnderlyingBuiltin type trait.
Header file for the exception macros of the math module.
Header file for the implementation of a dynamic MxN matrix.
Header file for the complete DynamicVector implementation.
Header file for the conjugate shim.
Header file for run time assertion macros.
#define BLAZE_CONSTRAINT_MUST_BE_NUMERIC_TYPE(T)
Constraint on the data type.In case the given data type T is not a numeric (integral or floating poin...
Definition: Numeric.h:61
bool isHermitian(const DenseMatrix< MT, SO > &dm)
Checks if the given dense matrix is Hermitian.
Definition: DenseMatrix.h:919
bool isSymmetric(const DenseMatrix< MT, SO > &dm)
Checks if the given dense matrix is symmetric.
Definition: DenseMatrix.h:841
Header file for all basic DenseMatrix functionality.
Header file for the complete IdentityMatrix implementation.
Header file for the real shim.
decltype(auto) conj(const DenseMatrix< MT, SO > &dm)
Returns a matrix containing the complex conjugate of each single element of dm.
Definition: DMatMapExpr.h:1321
BLAZE_ALWAYS_INLINE bool isSquare(const Matrix< MT, SO > &matrix) noexcept
Checks if the given matrix is a square matrix.
Definition: Matrix.h:908
#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