Blaze 3.9
sytrs.h
Go to the documentation of this file.
1//=================================================================================================
33//=================================================================================================
34
35#ifndef _BLAZE_MATH_LAPACK_SYTRS_H_
36#define _BLAZE_MATH_LAPACK_SYTRS_H_
37
38
39//*************************************************************************************************
40// Includes
41//*************************************************************************************************
42
43#include <blaze/math/Aliases.h>
55#include <blaze/util/Assert.h>
58
59
60namespace blaze {
61
62//=================================================================================================
63//
64// LAPACK LDLT-BASED SUBSTITUTION FUNCTIONS (SYTRS)
65//
66//=================================================================================================
67
68//*************************************************************************************************
71template< typename MT, bool SO, typename VT, bool TF >
72void sytrs( const DenseMatrix<MT,SO>& A, DenseVector<VT,TF>& b, char uplo,
73 const blas_int_t* ipiv );
74
75template< typename MT1, bool SO1, typename MT2, bool SO2 >
76void sytrs( const DenseMatrix<MT1,SO1>& A, DenseMatrix<MT2,SO2>& B, char uplo,
77 const blas_int_t* ipiv );
79//*************************************************************************************************
80
81
82//*************************************************************************************************
169template< typename MT // Type of the system matrix
170 , bool SO // Storage order of the system matrix
171 , typename VT // Type of the right-hand side vector
172 , bool TF > // Transpose flag of the right-hand side vector
173inline void sytrs( const DenseMatrix<MT,SO>& A, DenseVector<VT,TF>& b, char uplo,
174 const blas_int_t* ipiv )
175{
181
186
187 if( !isSquare( *A ) ) {
188 BLAZE_THROW_INVALID_ARGUMENT( "Invalid non-square matrix provided" );
189 }
190
191 if( (*b).size() != (*A).rows() ) {
192 BLAZE_THROW_INVALID_ARGUMENT( "Invalid right-hand side vector provided" );
193 }
194
195 if( uplo != 'L' && uplo != 'U' ) {
196 BLAZE_THROW_INVALID_ARGUMENT( "Invalid uplo argument provided" );
197 }
198
199 blas_int_t n ( numeric_cast<blas_int_t>( (*A).rows() ) );
200 blas_int_t nrhs( 1 );
201 blas_int_t lda ( numeric_cast<blas_int_t>( (*A).spacing() ) );
202 blas_int_t ldb ( numeric_cast<blas_int_t>( (*b).size() ) );
203 blas_int_t info( 0 );
204
205 if( n == 0 ) {
206 return;
207 }
208
209 if( IsRowMajorMatrix_v<MT> ) {
210 ( uplo == 'L' )?( uplo = 'U' ):( uplo = 'L' );
211 }
212
213 sytrs( uplo, n, nrhs, (*A).data(), lda, ipiv, (*b).data(), ldb, &info );
214
215 BLAZE_INTERNAL_ASSERT( info == 0, "Invalid function argument" );
216}
217//*************************************************************************************************
218
219
220//*************************************************************************************************
307template< typename MT1 // Type of the system matrix
308 , bool SO1 // Storage order of the system matrix
309 , typename MT2 // Type of the right-hand side matrix
310 , bool SO2 > // Storage order of the right-hand side matrix
311inline void sytrs( const DenseMatrix<MT1,SO1>& A, DenseMatrix<MT2,SO2>& B, char uplo,
312 const blas_int_t* ipiv )
313{
319
325
326 if( !isSquare( *A ) ) {
327 BLAZE_THROW_INVALID_ARGUMENT( "Invalid non-square matrix provided" );
328 }
329
330 if( uplo != 'L' && uplo != 'U' ) {
331 BLAZE_THROW_INVALID_ARGUMENT( "Invalid uplo argument provided" );
332 }
333
334 blas_int_t n ( numeric_cast<blas_int_t>( (*A).rows() ) );
335 blas_int_t mrhs( numeric_cast<blas_int_t>( SO2 ? (*B).rows() : (*B).columns() ) );
336 blas_int_t nrhs( numeric_cast<blas_int_t>( SO2 ? (*B).columns() : (*B).rows() ) );
337 blas_int_t lda ( numeric_cast<blas_int_t>( (*A).spacing() ) );
338 blas_int_t ldb ( numeric_cast<blas_int_t>( (*B).spacing() ) );
339 blas_int_t info( 0 );
340
341 if( n != mrhs ) {
342 BLAZE_THROW_INVALID_ARGUMENT( "Invalid right-hand side matrix provided" );
343 }
344
345 if( n == 0 ) {
346 return;
347 }
348
349 if( IsRowMajorMatrix_v<MT1> ) {
350 ( uplo == 'L' )?( uplo = 'U' ):( uplo = 'L' );
351 }
352
353 sytrs( uplo, n, nrhs, (*A).data(), lda, ipiv, (*B).data(), ldb, &info );
354
355 BLAZE_INTERNAL_ASSERT( info == 0, "Invalid function argument" );
356}
357//*************************************************************************************************
358
359} // namespace blaze
360
361#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.
Header file for the IsRowMajorMatrix type trait.
Constraint on the data type.
Cast operators for numeric types.
Data type constraint.
Header file for the CLAPACK sytrs wrapper functions.
Base class for dense matrices.
Definition: DenseMatrix.h:82
Base class for N-dimensional dense vectors.
Definition: DenseVector.h:77
Constraint on the data type.
Header file for the DenseMatrix base class.
Header file for the DenseVector base class.
#define BLAZE_CONSTRAINT_MUST_BE_SAME_TYPE(A, B)
Data type constraint.
Definition: SameType.h:71
void sytrs(const DenseMatrix< MT1, SO1 > &A, DenseMatrix< MT2, SO2 > &B, char uplo, const blas_int_t *ipiv)
LAPACK kernel for the substitution step of solving a symmetric indefinite linear system of equations ...
Definition: sytrs.h:311
#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
#define BLAZE_CONSTRAINT_MUST_HAVE_CONST_DATA_ACCESS(T)
Constraint on the data type.
Definition: ConstDataAccess.h:61
int32_t blas_int_t
Signed integer type used in the BLAS/LAPACK wrapper functions.
Definition: Types.h:64
bool isSquare(const Matrix< MT, SO > &matrix) noexcept
Checks if the given matrix is a square matrix.
Definition: Matrix.h:1383
#define BLAZE_INTERNAL_ASSERT(expr, msg)
Run time assertion macro for internal checks.
Definition: Assert.h:101
#define BLAZE_THROW_INVALID_ARGUMENT(MESSAGE)
Macro for the emission of a std::invalid_argument exception.
Definition: Exception.h:235
Header file for the exception macros of the math module.