Blaze 3.9
geqrf.h
Go to the documentation of this file.
1//=================================================================================================
33//=================================================================================================
34
35#ifndef _BLAZE_MATH_LAPACK_GEQRF_H_
36#define _BLAZE_MATH_LAPACK_GEQRF_H_
37
38
39//*************************************************************************************************
40// Includes
41//*************************************************************************************************
42
43#include <memory>
44#include <blaze/math/Aliases.h>
53#include <blaze/util/Assert.h>
55
56
57namespace blaze {
58
59//=================================================================================================
60//
61// LAPACK QR DECOMPOSITION FUNCTIONS (GEQRF)
62//
63//=================================================================================================
64
65//*************************************************************************************************
68template< typename MT, bool SO >
69void geqrf( DenseMatrix<MT,SO>& A, ElementType_t<MT>* tau );
71//*************************************************************************************************
72
73
74//*************************************************************************************************
116template< typename MT // Type of the dense matrix
117 , bool SO > // Storage order of the dense matrix
119{
125
126 using ET = ElementType_t<MT>;
127
128 blas_int_t m ( numeric_cast<blas_int_t>( SO ? (*A).rows() : (*A).columns() ) );
129 blas_int_t n ( numeric_cast<blas_int_t>( SO ? (*A).columns() : (*A).rows() ) );
130 blas_int_t lda ( numeric_cast<blas_int_t>( (*A).spacing() ) );
131 blas_int_t info( 0 );
132
133 if( m == 0 || n == 0 ) {
134 return;
135 }
136
137 blas_int_t lwork( ( SO ? n : m ) * lda );
138 const std::unique_ptr<ET[]> work( new ET[lwork] );
139
140 if( SO ) {
141 geqrf( m, n, (*A).data(), lda, tau, work.get(), lwork, &info );
142 }
143 else {
144 gelqf( m, n, (*A).data(), lda, tau, work.get(), lwork, &info );
145 }
146
147 BLAZE_INTERNAL_ASSERT( info == 0, "Invalid argument for QR decomposition" );
148}
149//*************************************************************************************************
150
151} // namespace blaze
152
153#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 gelqf wrapper functions.
Header file for the CLAPACK geqrf wrapper functions.
Base class for dense matrices.
Definition: DenseMatrix.h:82
Constraint on the data type.
Header file for the DenseMatrix base class.
void geqrf(DenseMatrix< MT, SO > &A, ElementType_t< MT > *tau)
LAPACK kernel for the QR decomposition of the given dense matrix.
Definition: geqrf.h:118
void gelqf(blas_int_t m, blas_int_t n, float *A, blas_int_t lda, float *tau, float *work, blas_int_t lwork, blas_int_t *info)
LAPACK kernel for the LQ decomposition of the given dense single precision column-major matrix.
Definition: gelqf.h:153
#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