jacobi_eig.cpp File Reference

Sequential cyclic Jacobi eigensolver. More...

#include "impl.hpp"
#include <algorithm>
#include <cmath>
#include <stdexcept>

Go to the source code of this file.

Namespaces
namespace	num
namespace	num::backends
namespace	num::backends::seq

Functions
EigenResult	num::backends::seq::eig_sym (const Matrix &A, real tol, idx max_sweeps)

Detailed Description

Sequential cyclic Jacobi eigensolver.

Uses \(t=\operatorname{sign}(\tau)/(|\tau|+\sqrt{1+\tau^2})\) with \(\tau=(a_{qq}-a_{pp})/(2a_{pq})\).

Definition in file jacobi_eig.cpp.