29 Numerics library [numerics]

29.9 Basic linear algebra algorithms [linalg]

29.9.15 BLAS 3 algorithms [linalg.algs.blas3]

29.9.15.5 Rank-2k update of a symmetric or Hermitian matrix [linalg.algs.blas3.rank2k]

1
#
[Note 1: 
These functions correspond to the BLAS functions xSYR2K and xHER2K[bib].
— end note]
2
#
The following elements apply to all functions in [linalg.algs.blas3.rank2k].
3
#
For any function F in this subclause with a parameter named t, an InMat3 template parameter, and a function parameter InMat3 E, t applies to accesses done through the parameter E.
F only accesses the triangle of E specified by t.
For accesses of diagonal elements E[i, i], F only uses the value real-if-needed(E[i, i]) if the name of F starts with hermitian.
For accesses E[i, j] outside the triangle specified by t, F only uses the value
  • (3.1)
    conj-if-needed(E[j, i]) if the name of F starts with hermitian, or
  • (3.2)
    E[j, i] if the name of F starts with symmetric.
4
#
Mandates:
  • (4.1)
    If OutMat has layout_blas_packed layout, then the layout's Triangle template argument has the same type as the function's Triangle template argument;
  • (4.2)
    If the function has an InMat3 template parameter and if InMat3 has layout_blas_packed layout, then the layout's Triangle template argument has the same type as the function's Triangle template argument;
  • (4.3)
    possibly-multipliable<decltype(A), decltype(transposed(B)), decltype(C)>() is true;
  • (4.4)
    possibly-multipliable<decltype(B), decltype(transposed(A)), decltype(C)>() is true; and
  • (4.5)
    possibly-addable<decltype(C), decltype(E), decltype(C)>() is true for those overloads with an E parameter.
5
#
Preconditions:
  • (5.1)
    multipliable(A, transposed(B), C) is true,
  • (5.2)
    multipliable(B, transposed(A), C) is true, and
    [Note 2: 
    This and the previous imply that C is square.
    — end note]
  • (5.3)
    addable(C, E, C) is true for those overloads with an E parameter.
6
#
Complexity: .
7
#
Remarks: C may alias E for those overloads with an E parameter.
🔗
template<in-matrix InMat1, in-matrix InMat2, possibly-packed-out-matrix OutMat, class Triangle> void symmetric_matrix_rank_2k_update(InMat1 A, InMat2 B, OutMat C, Triangle t); template<class ExecutionPolicy, in-matrix InMat1, in-matrix InMat2, possibly-packed-out-matrix OutMat, class Triangle> void symmetric_matrix_rank_2k_update(ExecutionPolicy&& exec, InMat1 A, InMat2 B, OutMat C, Triangle t);
8
#
Effects: Computes .
🔗
template<in-matrix InMat1, in-matrix InMat2, possibly-packed-out-matrix OutMat, class Triangle> void hermitian_matrix_rank_2k_update(InMat1 A, InMat2 B, OutMat C, Triangle t); template<class ExecutionPolicy, in-matrix InMat1, in-matrix InMat2, possibly-packed-out-matrix OutMat, class Triangle> void hermitian_matrix_rank_2k_update(ExecutionPolicy&& exec, InMat1 A, InMat2 B, OutMat C, Triangle t);
9
#
Effects: Computes .
🔗
template<in-matrix InMat1, in-matrix InMat2, in-matrix InMat3, possibly-packed-out-matrix OutMat, class Triangle> void symmetric_matrix_rank_2k_update(InMat1 A, InMat2 B, InMat3 E, OutMat C, Triangle t); template<class ExecutionPolicy, in-matrix InMat1, in-matrix InMat2, in-matrix InMat3, possibly-packed-out-matrix OutMat, class Triangle> void symmetric_matrix_rank_2k_update(ExecutionPolicy&& exec, InMat1 A, InMat2 B, InMat3 E, OutMat C, Triangle t);
Computes .
🔗
template<in-matrix InMat1, in-matrix InMat2, in-matrix InMat3, possibly-packed-out-matrix OutMat, class Triangle> void hermitian_matrix_rank_2k_update(InMat1 A, InMat2 B, InMat3 E, OutMat C, Triangle t); template<class ExecutionPolicy, in-matrix InMat1, in-matrix InMat2, in-matrix InMat3, possibly-packed-out-matrix OutMat, class Triangle> void hermitian_matrix_rank_2k_update(ExecutionPolicy&& exec, InMat1 A, InMat2 B, InMat3 E, OutMat C, Triangle t);
10
#
Effects: Computes .