lantr.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_LANTR_HH
#define TLAPACK_LEGACY_LANTR_HH
 
#include "tlapack/base/types.hpp"
#include "tlapack/lapack/lantr.hpp"
 
namespace tlapack {
namespace legacy {
 
    template <typename TA>
    real_type<TA> lantr(Norm normType,
                        Uplo uplo,
                        Diag diag,
                        idx_t m,
                        idx_t n,
                        const TA* A,
                        idx_t lda)
    {
        using internal::create_matrix;
 
        // check arguments
        tlapack_check_false(normType != Norm::Fro && normType != Norm::Inf &&
                            normType != Norm::Max && normType != Norm::One);
        tlapack_check_false(uplo != Uplo::Lower && uplo != Uplo::Upper);
        tlapack_check_false(diag != Diag::NonUnit && diag != Diag::Unit);
 
        // quick return
        if (m == 0 || n == 0) return 0;
 
        // Matrix views
        auto A_ = create_matrix<TA>((TA*)A, m, n, lda);
 
        return lantr(normType, uplo, diag, A_);
    }
 
}  // namespace legacy
}  // namespace tlapack
 
#endif  // TLAPACK_LEGACY_LANTR_HH