getrf.h
Go to the documentation of this file.
1 //=================================================================================================
33 //=================================================================================================
34 
35 #ifndef _BLAZE_MATH_LAPACK_GETRF_H_
36 #define _BLAZE_MATH_LAPACK_GETRF_H_
37 
38 
39 //*************************************************************************************************
40 // Includes
41 //*************************************************************************************************
42 
43 #include <boost/cast.hpp>
50 #include <blaze/util/Assert.h>
51 #include <blaze/util/Complex.h>
53 
54 
55 namespace blaze {
56 
57 //=================================================================================================
58 //
59 // LAPACK FORWARD DECLARATIONS
60 //
61 //=================================================================================================
62 
63 //*************************************************************************************************
65 extern "C" {
66 
67 void sgetrf_( int* m, int* n, float* A, int* lda, int* ipiv, int* info );
68 void dgetrf_( int* m, int* n, double* A, int* lda, int* ipiv, int* info );
69 void cgetrf_( int* m, int* n, float* A, int* lda, int* ipiv, int* info );
70 void zgetrf_( int* m, int* n, double* A, int* lda, int* ipiv, int* info );
71 
72 }
74 //*************************************************************************************************
75 
76 
77 
78 
79 //=================================================================================================
80 //
81 // LAPACK LU DECOMPOSITION FUNCTIONS (GETRF)
82 //
83 //=================================================================================================
84 
85 //*************************************************************************************************
88 inline void getrf( int m, int n, float* A, int lda, int* ipiv, int* info );
89 
90 inline void getrf( int m, int n, double* A, int lda, int* ipiv, int* info );
91 
92 inline void getrf( int m, int n, complex<float>* A, int lda, int* ipiv, int* info );
93 
94 inline void getrf( int m, int n, complex<double>* A, int lda, int* ipiv, int* info );
95 
96 template< typename MT, bool SO >
97 inline void getrf( DenseMatrix<MT,SO>& A, int* ipiv );
99 //*************************************************************************************************
100 
101 
102 //*************************************************************************************************
139 inline void getrf( int m, int n, float* A, int lda, int* ipiv, int* info )
140 {
141  sgetrf_( &m, &n, A, &lda, ipiv, info );
142 }
143 //*************************************************************************************************
144 
145 
146 //*************************************************************************************************
183 inline void getrf( int m, int n, double* A, int lda, int* ipiv, int* info )
184 {
185  dgetrf_( &m, &n, A, &lda, ipiv, info );
186 }
187 //*************************************************************************************************
188 
189 
190 //*************************************************************************************************
227 inline void getrf( int m, int n, complex<float>* A, int lda, int* ipiv, int* info )
228 {
229  BLAZE_STATIC_ASSERT( sizeof( complex<float> ) == 2UL*sizeof( float ) );
230 
231  cgetrf_( &m, &n, reinterpret_cast<float*>( A ), &lda, ipiv, info );
232 }
233 //*************************************************************************************************
234 
235 
236 //*************************************************************************************************
273 inline void getrf( int m, int n, complex<double>* A, int lda, int* ipiv, int* info )
274 {
275  BLAZE_STATIC_ASSERT( sizeof( complex<double> ) == 2UL*sizeof( double ) );
276 
277  zgetrf_( &m, &n, reinterpret_cast<double*>( A ), &lda, ipiv, info );
278 }
279 //*************************************************************************************************
280 
281 
282 //*************************************************************************************************
329 template< typename MT // Type of the dense matrix
330  , bool SO > // Storage order of the dense matrix
331 inline void getrf( DenseMatrix<MT,SO>& A, int* ipiv )
332 {
333  using boost::numeric_cast;
334 
339 
340  int m ( numeric_cast<int>( SO ? (~A).rows() : (~A).columns() ) );
341  int n ( numeric_cast<int>( SO ? (~A).columns() : (~A).rows() ) );
342  int lda ( numeric_cast<int>( (~A).spacing() ) );
343  int info( 0 );
344 
345  if( m == 0 || n == 0 ) {
346  return;
347  }
348 
349  getrf( m, n, (~A).data(), lda, ipiv, &info );
350 
351  BLAZE_INTERNAL_ASSERT( info >= 0, "Invalid argument for LU decomposition" );
352 }
354 //*************************************************************************************************
355 
356 } // namespace blaze
357 
358 #endif
Constraint on the data type.
Log level for high-level information.
Definition: LogLevel.h:80
#define BLAZE_CONSTRAINT_MUST_HAVE_MUTABLE_DATA_ACCESS(T)
Constraint on the data type.In case the given data type T does not provide low-level data access to m...
Definition: MutableDataAccess.h:79
#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
BLAZE_ALWAYS_INLINE size_t rows(const Matrix< MT, SO > &matrix)
Returns the current number of rows of the matrix.
Definition: Matrix.h:308
Constraint on the data type.
#define BLAZE_CONSTRAINT_MUST_NOT_BE_ADAPTOR_TYPE(T)
Constraint on the data type.In case the given data type T is an adaptor type (as for instance LowerMa...
Definition: Adaptor.h:118
Namespace of the Blaze C++ math library.
Definition: Blaze.h:57
Compile time assertion.
Header file for the DenseMatrix base class.
Type ElementType
Type of the sparse matrix elements.
Definition: CompressedMatrix.h:2586
Constraint on the data type.
Constraint on the data type.
const bool spacing
Adding an additional spacing line between two log messages.This setting gives the opportunity to add ...
Definition: Logging.h:70
void getrf(int m, int n, float *A, int lda, int *ipiv, int *info)
LAPACK kernel for the LU decomposition of the given dense general single precision column-major matri...
Definition: getrf.h:139
Header file for run time assertion macros.
#define BLAZE_CONSTRAINT_MUST_BE_BLAS_COMPATIBLE_TYPE(T)
Constraint on the data type.In case the given data type T is not a BLAS compatible data type (i...
Definition: BlasCompatible.h:79
Header file for the IsRowMajorMatrix type trait.
BLAZE_ALWAYS_INLINE size_t columns(const Matrix< MT, SO > &matrix)
Returns the current number of columns of the matrix.
Definition: Matrix.h:324
Header file for the complex data type.
#define BLAZE_STATIC_ASSERT(expr)
Compile time assertion macro.In case of an invalid compile time expression, a compilation error is cr...
Definition: StaticAssert.h:143
#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