getrf.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_GETRF_HH
#define TLAPACK_GETRF_HH
 
#include "tlapack/base/utils.hpp"
#include "tlapack/lapack/getrf_level0.hpp"
#include "tlapack/lapack/getrf_recursive.hpp"
 
namespace tlapack {
 
enum class GetrfVariant : char { Level0 = '0', Recursive = 'R' };
 
struct GetrfOpts {
    GetrfVariant variant = GetrfVariant::Recursive;
};
 
template <TLAPACK_MATRIX matrix_t, TLAPACK_VECTOR piv_t>
int getrf(matrix_t& A, piv_t& piv, const GetrfOpts& opts = {})
{
    // Call variant
    if (opts.variant == GetrfVariant::Recursive)
        return getrf_recursive(A, piv);
    else
        return getrf_level0(A, piv);
}
 
}  // namespace tlapack
 
#endif  // TLAPACK_GETRF_HH