Loading [MathJax]/extensions/tex2jax.js
<T>LAPACK 0.1.2
C++ Template Linear Algebra PACKage
All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Modules Pages Concepts
her.hpp
Go to the documentation of this file.
1
3//
4// Copyright (c) 2017-2021, University of Tennessee. All rights reserved.
5// Copyright (c) 2025, University of Colorado Denver. All rights reserved.
6//
7// This file is part of <T>LAPACK.
8// <T>LAPACK is free software: you can redistribute it and/or modify it under
9// the terms of the BSD 3-Clause license. See the accompanying LICENSE file.
10
11#ifndef TLAPACK_BLAS_HER_HH
12#define TLAPACK_BLAS_HER_HH
13
15
16namespace tlapack {
17
45 enable_if_t<(
46 /* Requires: */
48 int> = 0,
49 class T = type_t<matrixA_t>,
53void her(Uplo uplo, const alpha_t& alpha, const vectorX_t& x, matrixA_t& A)
54{
55 // data traits
56 using idx_t = size_type<matrixA_t>;
58
59 // constants
60 const idx_t n = nrows(A);
61
62 // check arguments
64 tlapack_check_false(size(x) != n);
65 tlapack_check_false(ncols(A) != n);
66
67 if (uplo == Uplo::Upper) {
68 for (idx_t j = 0; j < n; ++j) {
69 const scalar_t tmp = alpha * conj(x[j]);
70 for (idx_t i = 0; i < j; ++i)
71 A(i, j) += x[i] * tmp;
72 A(j, j) =
73 real(A(j, j)) + real(x[j]) * real(tmp) - imag(x[j]) * imag(tmp);
74 }
75 }
76 else {
77 for (idx_t j = 0; j < n; ++j) {
78 const scalar_t tmp = alpha * conj(x[j]);
79 A(j, j) =
80 real(A(j, j)) + real(x[j]) * real(tmp) - imag(x[j]) * imag(tmp);
81 for (idx_t i = j + 1; i < n; ++i)
82 A(i, j) += x[i] * tmp;
83 }
84 }
85}
86
87#ifdef TLAPACK_USE_LAPACKPP
88
92 class T = type_t<matrixA_t>,
96void her(Uplo uplo, const alpha_t alpha, const vectorX_t& x, matrixA_t& A)
97{
98 // Legacy objects
99 auto A_ = legacy_matrix(A);
100 auto x_ = legacy_vector(x);
101
102 // Constants to forward
103 constexpr Layout L = layout<matrixA_t>;
104 const auto& n = A_.n;
105
106 return ::blas::her((::blas::Layout)L, (::blas::Uplo)uplo, n, alpha, x_.ptr,
107 x_.inc, A_.ptr, A_.ldim);
108}
109
110#endif
111
112} // namespace tlapack
113
114#endif // #ifndef TLAPACK_BLAS_HER_HH
#define TLAPACK_LEGACY_VECTOR
Macro for tlapack::concepts::LegacyVector compatible with C++17.
Definition concepts.hpp:954
#define TLAPACK_LEGACY_MATRIX
Macro for tlapack::concepts::LegacyMatrix compatible with C++17.
Definition concepts.hpp:951
#define TLAPACK_VECTOR
Macro for tlapack::concepts::Vector compatible with C++17.
Definition concepts.hpp:906
#define TLAPACK_MATRIX
Macro for tlapack::concepts::Matrix compatible with C++17.
Definition concepts.hpp:896
#define TLAPACK_REAL
Macro for tlapack::concepts::Real compatible with C++17.
Definition concepts.hpp:918
void her(Uplo uplo, const alpha_t &alpha, const vectorX_t &x, matrixA_t &A)
Hermitian matrix rank-1 update:
Definition her.hpp:53
#define tlapack_check_false(cond)
Throw an error if cond is true.
Definition exceptionHandling.hpp:113
Concept for types that represent tlapack::Uplo.
Sort the numbers in D in increasing order (if ID = 'I') or in decreasing order (if ID = 'D' ).
Definition arrayTraits.hpp:15
enable_if_t<(!allow_optblas< T1, Ts... >), int > disable_if_allow_optblas_t
Disable if the list of types allows optimized BLAS library.
Definition utils.hpp:303
typename traits::real_type_traits< Types..., int >::type real_type
The common real type of the list of types.
Definition scalar_type_traits.hpp:113
constexpr real_type< T > real(const T &x) noexcept
Extends std::real() to real datatypes.
Definition utils.hpp:71
constexpr T conj(const T &x) noexcept
Extends std::conj() to real datatypes.
Definition utils.hpp:100
constexpr real_type< T > imag(const T &x) noexcept
Extends std::imag() to real datatypes.
Definition utils.hpp:86
Uplo
Definition types.hpp:50
@ Upper
0 <= i <= j, 0 <= j <= n.
@ Lower
0 <= i <= m, 0 <= j <= i.
Layout
Definition types.hpp:29