lacpy.hpp

Go to the documentation of this file.

 
//
// 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_LACPY_HH
#define TLAPACK_LEGACY_LACPY_HH
 
#include "tlapack/lapack/lacpy.hpp"
 
namespace tlapack {
namespace legacy {
 
    template <class uplo_t, typename TA, typename TB>
    void lacpy(
        uplo_t uplo, idx_t m, idx_t n, const TA* A, idx_t lda, TB* B, idx_t ldb)
    {
        using internal::create_matrix;
 
        // check arguments
        tlapack_check_false(uplo != Uplo::Lower && uplo != Uplo::Upper &&
                            uplo != Uplo::General);
 
        // Matrix views
        const auto A_ = create_matrix<TA>((TA*)A, m, n, lda);
        auto B_ = create_matrix<TB>(B, m, n, ldb);
 
        lacpy(uplo, A_, B_);
    }
 
    template <typename TA, typename TB>
    void lacpy(MatrixType matrixtype,
               idx_t m,
               idx_t n,
               const TA* A,
               idx_t lda,
               TB* B,
               idx_t ldb)
    {
        if (matrixtype == MatrixType::Upper) {
            lacpy(UPPER_TRIANGLE, m, n, A, lda, B, ldb);
        }
        else if (matrixtype == MatrixType::Lower) {
            lacpy(LOWER_TRIANGLE, m, n, A, lda, B, ldb);
        }
        else {
            lacpy(GENERAL, m, n, A, lda, B, ldb);
        }
    }
 
}  // namespace legacy
}  // namespace tlapack
 
#endif  // TLAPACK_LEGACY_LACPY_HH