eig.cpp

Go to the documentation of this file.

/// @file eigen/eig.cpp
/// @brief Full symmetric eigendecomposition dispatcher.
 
#include "backends/lapack/impl.hpp"
#include "backends/omp/impl.hpp"
#include "backends/seq/impl.hpp"
#include "linalg/eigen/jacobi_eig.hpp"
 
namespace num {
 
EigenResult eig_sym(const Matrix& A, real tol, idx max_sweeps, Backend backend) {
    switch (backend) {
        case Backend::lapack:
            return backends::lapack::eig_sym(A);
        case Backend::omp:
            return backends::omp::eig_sym(A, tol, max_sweeps);
        default:
            return backends::seq::eig_sym(A, tol, max_sweeps);
    }
}
 
} // namespace num