Blaze 3.9
ungrq.h
Go to the documentation of this file.
1//=================================================================================================
33//=================================================================================================
34
35#ifndef _BLAZE_MATH_LAPACK_UNGRQ_H_
36#define _BLAZE_MATH_LAPACK_UNGRQ_H_
37
38
39//*************************************************************************************************
40// Includes
41//*************************************************************************************************
42
43#include <memory>
44#include <blaze/math/Aliases.h>
53#include <blaze/util/Assert.h>
56
57
58namespace blaze {
59
60//=================================================================================================
61//
62// LAPACK FUNCTIONS TO RECONSTRUCT Q FROM A RQ DECOMPOSITION (UNGRQ)
63//
64//=================================================================================================
65
66//*************************************************************************************************
69template< typename MT, bool SO >
70void ungrq( DenseMatrix<MT,SO>& A, const ElementType_t<MT>* tau );
72//*************************************************************************************************
73
74
75//*************************************************************************************************
120template< typename MT, bool SO >
121inline void ungrq( DenseMatrix<MT,SO>& A, const ElementType_t<MT>* tau )
122{
129
130 using ET = ElementType_t<MT>;
131
132 blas_int_t n ( numeric_cast<blas_int_t>( SO ? (*A).columns() : (*A).rows() ) );
133 blas_int_t m ( numeric_cast<blas_int_t>( SO ? (*A).rows() : (*A).columns() ) );
134 blas_int_t k ( min( m, n ) );
135 blas_int_t lda ( numeric_cast<blas_int_t>( (*A).spacing() ) );
136 blas_int_t info( 0 );
137
138 if( k == 0 ) {
139 return;
140 }
141
142 blas_int_t lwork( k*lda );
143 const std::unique_ptr<ET[]> work( new ET[lwork] );
144
145 if( SO ) {
146 const size_t offset( ( m > n )?( m - n ):( 0UL ) );
147 ungrq( k, n, k, (*A).data()+offset, lda, tau, work.get(), lwork, &info );
148 }
149 else {
150 const size_t offset( ( m < n )?( n - m ):( 0UL ) );
151 ungql( m, k, k, (*A).data(offset), lda, tau, work.get(), lwork, &info );
152 }
153
154 BLAZE_INTERNAL_ASSERT( info == 0, "Invalid argument for Q reconstruction" );
155}
156//*************************************************************************************************
157
158} // namespace blaze
159
160#endif
Constraint on the data type.
Header file for auxiliary alias declarations.
typename T::ElementType ElementType_t
Alias declaration for nested ElementType type definitions.
Definition: Aliases.h:190
Header file for run time assertion macros.
Constraint on the data type.
Constraint on the data type.
Constraint on the data type.
Cast operators for numeric types.
Header file for the CLAPACK ungql wrapper functions.
Header file for the CLAPACK ungrq wrapper functions.
Base class for dense matrices.
Definition: DenseMatrix.h:82
Constraint on the data type.
Constraint on the data type.
Header file for the DenseMatrix base class.
#define BLAZE_CONSTRAINT_MUST_BE_COMPLEX_TYPE(T)
Constraint on the data type.
Definition: Complex.h:62
decltype(auto) min(const DenseMatrix< MT1, SO1 > &lhs, const DenseMatrix< MT2, SO2 > &rhs)
Computes the componentwise minimum of the dense matrices lhs and rhs.
Definition: DMatDMatMapExpr.h:1339
void ungrq(DenseMatrix< MT, SO > &A, const ElementType_t< MT > *tau)
LAPACK kernel for the reconstruction of the orthogonal matrix Q from a RQ decomposition.
Definition: ungrq.h:121
void ungql(blas_int_t m, blas_int_t n, blas_int_t k, complex< float > *A, blas_int_t lda, const complex< float > *tau, complex< float > *work, blas_int_t lwork, blas_int_t *info)
LAPACK kernel for the reconstruction of the orthogonal matrix Q from a QL decomposition.
Definition: ungql.h:127
#define BLAZE_CONSTRAINT_MUST_BE_BLAS_COMPATIBLE_TYPE(T)
Constraint on the data type.
Definition: BLASCompatible.h:61
#define BLAZE_CONSTRAINT_MUST_BE_CONTIGUOUS_TYPE(T)
Constraint on the data type.
Definition: Contiguous.h:61
#define BLAZE_CONSTRAINT_MUST_NOT_BE_COMPUTATION_TYPE(T)
Constraint on the data type.
Definition: Computation.h:81
#define BLAZE_CONSTRAINT_MUST_NOT_BE_ADAPTOR_TYPE(T)
Constraint on the data type.
Definition: Adaptor.h:81
#define BLAZE_CONSTRAINT_MUST_HAVE_MUTABLE_DATA_ACCESS(T)
Constraint on the data type.
Definition: MutableDataAccess.h:61
int32_t blas_int_t
Signed integer type used in the BLAS/LAPACK wrapper functions.
Definition: Types.h:64
#define BLAZE_INTERNAL_ASSERT(expr, msg)
Run time assertion macro for internal checks.
Definition: Assert.h:101