larfb.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_LARFB_HH
#define TLAPACK_LEGACY_LARFB_HH
 
#include "tlapack/lapack/larfb.hpp"
 
namespace tlapack {
namespace legacy {
 
    template <class side_t,
              class trans_t,
              class direction_t,
              class storev_t,
              typename TV,
              typename TC>
    int larfb(side_t side,
              trans_t trans,
              direction_t direction,
              storev_t storev,
              idx_t m,
              idx_t n,
              idx_t k,
              TV const* V,
              idx_t ldv,
              TV const* T,
              idx_t ldt,
              TC* C,
              idx_t ldc)
    {
        using internal::create_matrix;
 
        // check arguments
        tlapack_check_false(side != Side::Left && side != Side::Right);
        tlapack_check_false(trans != Op::NoTrans && trans != Op::ConjTrans &&
                            ((trans != Op::Trans) || is_complex<TV>));
        tlapack_check_false(direction != Direction::Backward &&
                            direction != Direction::Forward);
        tlapack_check_false(storev != StoreV::Columnwise &&
                            storev != StoreV::Rowwise);
 
        // Quick return
        if (m <= 0 || n <= 0) return 0;
 
        // Views
        const auto V_ = (storev == StoreV::Columnwise)
                            ? create_matrix<TV>(
                                  (TV*)V, (side == Side::Left) ? m : n, k, ldv)
                            : create_matrix<TV>(
                                  (TV*)V, k, (side == Side::Left) ? m : n, ldv);
        const auto T_ = create_matrix<TV>((TV*)T, k, k, ldt);
        auto C_ = create_matrix<TC>(C, m, n, ldc);
 
        return larfb(side, trans, direction, storev, V_, T_, C_);
    }
 
}  // namespace legacy
}  // namespace tlapack
 
#endif  // TLAPACK_LEGACY_LARFB_HH