unm2r.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_UNM2R_HH
#define TLAPACK_LEGACY_UNM2R_HH
 
#include "tlapack/lapack/unm2r.hpp"
 
namespace tlapack {
namespace legacy {
 
    template <class side_t, class trans_t, typename TA, typename TC>
    int unm2r(side_t side,
              trans_t trans,
              idx_t m,
              idx_t n,
              idx_t k,
              TA* A,
              idx_t lda,
              const real_type<TA, TC>* tau,
              TC* C,
              idx_t ldc)
    {
        using internal::create_matrix;
        using internal::create_vector;
 
        // check arguments
        tlapack_check_false(side != Side::Left && side != Side::Right);
        tlapack_check_false(trans != Op::NoTrans && trans != Op::Trans &&
                            trans != Op::ConjTrans);
        tlapack_check_false(m < 0);
        tlapack_check_false(n < 0);
        const idx_t q = (side == Side::Left) ? m : n;
        tlapack_check_false(k < 0 || k > q);
        tlapack_check_false(lda < q);
        tlapack_check_false(ldc < m);
 
        // quick return
        if ((m == 0) || (n == 0) || (k == 0)) return 0;
 
        // Matrix views
        const auto A_ = create_matrix<TA>((TA*)A, q, k, lda);
        const auto tau_ = create_vector((TA*)tau, k);
        auto C_ = create_matrix<TC>(C, m, n, ldc);
 
        return unm2r(side, trans, A, tau_, C_);
    }
 
}  // namespace legacy
}  // namespace tlapack
 
#endif  // TLAPACK_LEGACY_UNM2R_HH