lq.hpp

Go to the documentation of this file.

// SPDX-License-Identifier: BSD-2-Clause
 
#ifndef pRC_ALGORITHMS_LQ_H
#define pRC_ALGORITHMS_LQ_H
 
#include <prc/algorithms/qr.hpp>
#include <prc/core/tensor/declarations.hpp>
 
namespace pRC
{
    template<Size B = 32, class X, IsTensorish R = RemoveReference<X>>
        requires IsFloat<Value<R>> && (R::Dimension == 2)
    static inline constexpr auto lq(X &&input)
    {
        auto const [q, r] = qr<B>(adjoint(forward<X>(input)));
        return Tuple(eval(adjoint(r)), eval(adjoint(q)));
    }
}
#endif // pRC_ALGORITHMS_LQ_H