her.hpp

Go to the documentation of this file.

 
//
// Copyright (c) 2017-2021, University of Tennessee. All rights reserved.
// Copyright (c) 2021-2023, University of Colorado Denver. All rights reserved.
//
// This file is part of <T>LAPACK.
// <T>LAPACK is free software: you can redistribute it and/or modify it under
// the terms of the BSD 3-Clause license. See the accompanying LICENSE file.
 
#ifndef TLAPACK_LEGACY_HER_HH
#define TLAPACK_LEGACY_HER_HH
 
#include "tlapack/blas/her.hpp"
#include "tlapack/legacy_api/base/types.hpp"
#include "tlapack/legacy_api/base/utils.hpp"
 
namespace tlapack {
namespace legacy {
 
    template <typename TA, typename TX>
    void her(Layout layout,
             Uplo uplo,
             idx_t n,
             real_type<TA, TX> alpha,  // zher takes double alpha; use real
             TX const* x,
             int_t incx,
             TA* A,
             idx_t lda)
    {
        using internal::create_matrix;
        using real_t = real_type<TA, TX>;
 
        // check arguments
        tlapack_check_false(layout != Layout::ColMajor &&
                            layout != Layout::RowMajor);
        tlapack_check_false(uplo != Uplo::Lower && uplo != Uplo::Upper);
        tlapack_check_false(n < 0);
        tlapack_check_false(incx == 0);
        tlapack_check_false(lda < n);
 
        // quick return
        if (n == 0 || alpha == real_t(0)) return;
 
        // for row major, swap lower <=> upper
        if (layout == Layout::RowMajor) {
            uplo = (uplo == Uplo::Lower ? Uplo::Upper : Uplo::Lower);
        }
 
        // Matrix views
        auto A_ = create_matrix<TA>(A, n, n, lda);
 
        tlapack_expr_with_vector(x_, TX, n, x, incx, her(uplo, alpha, x_, A_));
    }
 
}  // namespace legacy
}  // namespace tlapack
 
#endif  //  #ifndef TLAPACK_LEGACY_HER_HH