ungl2.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_UNGL2_HH
#define TLAPACK_UNGL2_HH
 
#include "tlapack/base/utils.hpp"
#include "tlapack/lapack/ungq_level2.hpp"
 
namespace tlapack {
 
template <class T, TLAPACK_SMATRIX matrix_t, TLAPACK_VECTOR vector_t>
constexpr WorkInfo ungl2_worksize(const matrix_t& Q, const vector_t& tauw)
{
    using idx_t = size_type<matrix_t>;
    using range = pair<idx_t, idx_t>;
 
    // constants
    const idx_t k = nrows(Q);
 
    if (k > 1) {
        auto&& C = rows(Q, range{1, k});
        return larf_worksize<T>(RIGHT_SIDE, FORWARD, ROWWISE_STORAGE, row(Q, 0),
                                tauw[0], C);
    }
    return WorkInfo(0);
}
 
template <TLAPACK_SMATRIX matrix_t,
          TLAPACK_VECTOR vector_t,
          TLAPACK_WORKSPACE work_t>
int ungl2_work(matrix_t& Q, const vector_t& tauw, work_t& work)
{
    return ungq_level2_work(FORWARD, ROWWISE_STORAGE, Q, tauw, work);
}
 
template <TLAPACK_SMATRIX matrix_t, TLAPACK_VECTOR vector_t>
int ungl2(matrix_t& Q, const vector_t& tauw)
{
    using T = type_t<matrix_t>;
 
    // functor
    Create<matrix_t> new_matrix;
 
    // Allocates workspace
    WorkInfo workinfo = ungl2_worksize<T>(Q, tauw);
    std::vector<T> work_;
    auto work = new_matrix(work_, workinfo.m, workinfo.n);
 
    return ungl2_work(Q, tauw, work);
}
}  // namespace tlapack
#endif  // TLAPACK_UNGL2_HH