larft.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_LARFT_HH
#define TLAPACK_LEGACY_LARFT_HH
 
#include "tlapack/lapack/larft.hpp"
 
namespace tlapack {
namespace legacy {
 
    template <class direction_t, class storev_t, typename scalar_t>
    int larft(direction_t direction,
              storev_t storev,
              idx_t n,
              idx_t k,
              const scalar_t* V,
              idx_t ldV,
              const scalar_t* tau,
              scalar_t* T,
              idx_t ldT)
    {
        using internal::create_matrix;
        using internal::create_vector;
 
        // check arguments
        tlapack_check_false(direction != Direction::Forward &&
                            direction != Direction::Backward);
        tlapack_check_false(storev != StoreV::Columnwise &&
                            storev != StoreV::Rowwise);
        tlapack_check_false(n < 0);
        tlapack_check_false(k < 1);
        tlapack_check_false(ldV < ((storev == StoreV::Columnwise) ? n : k));
        tlapack_check_false(ldT < k);
 
        // Quick return
        if (n == 0 || k == 0) return 0;
 
        // Matrix views
        auto V_ = (storev == StoreV::Columnwise)
                      ? create_matrix<scalar_t>((scalar_t*)V, n, k, ldV)
                      : create_matrix<scalar_t>((scalar_t*)V, k, n, ldV);
        auto tau_ = create_vector((scalar_t*)tau, k);
        auto T_ = create_matrix<scalar_t>(T, k, k, ldT);
 
        return larft(direction, storev, V_, tau_, T_);
    }
 
}  // namespace legacy
}  // namespace tlapack
 
#endif  // TLAPACK_LEGACY_LARFT_HH