Contains routines to diagonalize mass matrices. More...
#include <limits>#include <cctype>#include <complex>#include <algorithm>#include <tuple>#include <Eigen/Core>#include <Eigen/SVD>#include <Eigen/Eigenvalues>#include <unsupported/Eigen/MatrixFunctions>
Include dependency graph for Linalg.hpp:
This graph shows which files directly or indirectly include this file:
Go to the source code of this file.
Classes | |
| struct | himalaya::FlipSignOp< Real > |
Namespaces | |
| himalaya | |
Macros | |
| #define | MAX_(i, j) (((i) > (j)) ? (i) : (j)) |
| #define | MIN_(i, j) (((i) < (j)) ? (i) : (j)) |
Typedefs | |
| using | himalaya::V2 = Eigen::Vector2d |
| real 2-vector More... | |
| using | himalaya::RM22 = Eigen::Matrix2d |
| real 2x2 matrix More... | |
Functions | |
| template<class Real , class Scalar , int M, int N> | |
| void | himalaya::svd_eigen (const Eigen::Matrix< Scalar, M, N > &m, Eigen::Array< Real, MIN_(M, N), 1 > &s, Eigen::Matrix< Scalar, M, M > *u, Eigen::Matrix< Scalar, N, N > *vh) |
| template<class Real , class Scalar , int N> | |
| void | himalaya::hermitian_eigen (const Eigen::Matrix< Scalar, N, N > &m, Eigen::Array< Real, N, 1 > &w, Eigen::Matrix< Scalar, N, N > *z) |
| template<int M, int N, class Real > | |
| void | himalaya::disna (const char &JOB, const Eigen::Array< Real, MIN_(M, N), 1 > &D, Eigen::Array< Real, MIN_(M, N), 1 > &SEP, int &INFO) |
| template<class Real , class Scalar , int M, int N> | |
| void | himalaya::svd_internal (const Eigen::Matrix< Scalar, M, N > &m, Eigen::Array< Real, MIN_(M, N), 1 > &s, Eigen::Matrix< Scalar, M, M > *u, Eigen::Matrix< Scalar, N, N > *vh) |
| template<class Real , class Scalar , int M, int N> | |
| void | himalaya::svd_errbd (const Eigen::Matrix< Scalar, M, N > &m, Eigen::Array< Real, MIN_(M, N), 1 > &s, Eigen::Matrix< Scalar, M, M > *u=0, Eigen::Matrix< Scalar, N, N > *vh=0, Real *s_errbd=0, Eigen::Array< Real, MIN_(M, N), 1 > *u_errbd=0, Eigen::Array< Real, MIN_(M, N), 1 > *v_errbd=0) |
| template<class Real , class Scalar , int M, int N> | |
| void | himalaya::svd (const Eigen::Matrix< Scalar, M, N > &m, Eigen::Array< Real, MIN_(M, N), 1 > &s, Eigen::Matrix< Scalar, M, M > &u, Eigen::Matrix< Scalar, N, N > &vh) |
| template<class Real , class Scalar , int M, int N> | |
| void | himalaya::svd (const Eigen::Matrix< Scalar, M, N > &m, Eigen::Array< Real, MIN_(M, N), 1 > &s, Eigen::Matrix< Scalar, M, M > &u, Eigen::Matrix< Scalar, N, N > &vh, Real &s_errbd) |
| template<class Real , class Scalar , int M, int N> | |
| void | himalaya::svd (const Eigen::Matrix< Scalar, M, N > &m, Eigen::Array< Real, MIN_(M, N), 1 > &s, Eigen::Matrix< Scalar, M, M > &u, Eigen::Matrix< Scalar, N, N > &vh, Real &s_errbd, Eigen::Array< Real, MIN_(M, N), 1 > &u_errbd, Eigen::Array< Real, MIN_(M, N), 1 > &v_errbd) |
| template<class Real , class Scalar , int M, int N> | |
| void | himalaya::svd (const Eigen::Matrix< Scalar, M, N > &m, Eigen::Array< Real, MIN_(M, N), 1 > &s) |
| template<class Real , class Scalar , int M, int N> | |
| void | himalaya::svd (const Eigen::Matrix< Scalar, M, N > &m, Eigen::Array< Real, MIN_(M, N), 1 > &s, Real &s_errbd) |
| template<class Real , class Scalar , int N> | |
| void | himalaya::diagonalize_hermitian_internal (const Eigen::Matrix< Scalar, N, N > &m, Eigen::Array< Real, N, 1 > &w, Eigen::Matrix< Scalar, N, N > *z) |
| template<class Real , class Scalar , int N> | |
| void | himalaya::diagonalize_hermitian_errbd (const Eigen::Matrix< Scalar, N, N > &m, Eigen::Array< Real, N, 1 > &w, Eigen::Matrix< Scalar, N, N > *z=0, Real *w_errbd=0, Eigen::Array< Real, N, 1 > *z_errbd=0) |
| template<class Real , class Scalar , int N> | |
| void | himalaya::diagonalize_hermitian (const Eigen::Matrix< Scalar, N, N > &m, Eigen::Array< Real, N, 1 > &w, Eigen::Matrix< Scalar, N, N > &z) |
| template<class Real , class Scalar , int N> | |
| void | himalaya::diagonalize_hermitian (const Eigen::Matrix< Scalar, N, N > &m, Eigen::Array< Real, N, 1 > &w, Eigen::Matrix< Scalar, N, N > &z, Real &w_errbd) |
| template<class Real , class Scalar , int N> | |
| void | himalaya::diagonalize_hermitian (const Eigen::Matrix< Scalar, N, N > &m, Eigen::Array< Real, N, 1 > &w, Eigen::Matrix< Scalar, N, N > &z, Real &w_errbd, Eigen::Array< Real, N, 1 > &z_errbd) |
| template<class Real , class Scalar , int N> | |
| void | himalaya::diagonalize_hermitian (const Eigen::Matrix< Scalar, N, N > &m, Eigen::Array< Real, N, 1 > &w) |
| template<class Real , class Scalar , int N> | |
| void | himalaya::diagonalize_hermitian (const Eigen::Matrix< Scalar, N, N > &m, Eigen::Array< Real, N, 1 > &w, Real &w_errbd) |
| template<class Real , int N> | |
| void | himalaya::diagonalize_symmetric_errbd (const Eigen::Matrix< std::complex< Real >, N, N > &m, Eigen::Array< Real, N, 1 > &s, Eigen::Matrix< std::complex< Real >, N, N > *u=0, Real *s_errbd=0, Eigen::Array< Real, N, 1 > *u_errbd=0) |
| template<class Real , int N> | |
| void | himalaya::diagonalize_symmetric (const Eigen::Matrix< std::complex< Real >, N, N > &m, Eigen::Array< Real, N, 1 > &s, Eigen::Matrix< std::complex< Real >, N, N > &u) |
| template<class Real , int N> | |
| void | himalaya::diagonalize_symmetric (const Eigen::Matrix< std::complex< Real >, N, N > &m, Eigen::Array< Real, N, 1 > &s, Eigen::Matrix< std::complex< Real >, N, N > &u, Real &s_errbd) |
| template<class Real , int N> | |
| void | himalaya::diagonalize_symmetric (const Eigen::Matrix< std::complex< Real >, N, N > &m, Eigen::Array< Real, N, 1 > &s, Eigen::Matrix< std::complex< Real >, N, N > &u, Real &s_errbd, Eigen::Array< Real, N, 1 > &u_errbd) |
| template<class Real , int N> | |
| void | himalaya::diagonalize_symmetric (const Eigen::Matrix< std::complex< Real >, N, N > &m, Eigen::Array< Real, N, 1 > &s) |
| template<class Real , int N> | |
| void | himalaya::diagonalize_symmetric (const Eigen::Matrix< std::complex< Real >, N, N > &m, Eigen::Array< Real, N, 1 > &s, Real &s_errbd) |
| template<class Real , int N> | |
| void | himalaya::diagonalize_symmetric_errbd (const Eigen::Matrix< Real, N, N > &m, Eigen::Array< Real, N, 1 > &s, Eigen::Matrix< std::complex< Real >, N, N > *u=0, Real *s_errbd=0, Eigen::Array< Real, N, 1 > *u_errbd=0) |
| template<class Real , int N> | |
| void | himalaya::diagonalize_symmetric (const Eigen::Matrix< Real, N, N > &m, Eigen::Array< Real, N, 1 > &s, Eigen::Matrix< std::complex< Real >, N, N > &u) |
| template<class Real , int N> | |
| void | himalaya::diagonalize_symmetric (const Eigen::Matrix< Real, N, N > &m, Eigen::Array< Real, N, 1 > &s, Eigen::Matrix< std::complex< Real >, N, N > &u, Real &s_errbd) |
| template<class Real , int N> | |
| void | himalaya::diagonalize_symmetric (const Eigen::Matrix< Real, N, N > &m, Eigen::Array< Real, N, 1 > &s, Eigen::Matrix< std::complex< Real >, N, N > &u, Real &s_errbd, Eigen::Array< Real, N, 1 > &u_errbd) |
| template<class Real , int N> | |
| void | himalaya::diagonalize_symmetric (const Eigen::Matrix< Real, N, N > &m, Eigen::Array< Real, N, 1 > &s) |
| template<class Real , int N> | |
| void | himalaya::diagonalize_symmetric (const Eigen::Matrix< Real, N, N > &m, Eigen::Array< Real, N, 1 > &s, Real &s_errbd) |
| template<class Real , class Scalar , int M, int N> | |
| void | himalaya::reorder_svd_errbd (const Eigen::Matrix< Scalar, M, N > &m, Eigen::Array< Real, MIN_(M, N), 1 > &s, Eigen::Matrix< Scalar, M, M > *u=0, Eigen::Matrix< Scalar, N, N > *vh=0, Real *s_errbd=0, Eigen::Array< Real, MIN_(M, N), 1 > *u_errbd=0, Eigen::Array< Real, MIN_(M, N), 1 > *v_errbd=0) |
| template<class Real , class Scalar , int M, int N> | |
| void | himalaya::reorder_svd (const Eigen::Matrix< Scalar, M, N > &m, Eigen::Array< Real, MIN_(M, N), 1 > &s, Eigen::Matrix< Scalar, M, M > &u, Eigen::Matrix< Scalar, N, N > &vh) |
| template<class Real , class Scalar , int M, int N> | |
| void | himalaya::reorder_svd (const Eigen::Matrix< Scalar, M, N > &m, Eigen::Array< Real, MIN_(M, N), 1 > &s, Eigen::Matrix< Scalar, M, M > &u, Eigen::Matrix< Scalar, N, N > &vh, Real &s_errbd) |
| template<class Real , class Scalar , int M, int N> | |
| void | himalaya::reorder_svd (const Eigen::Matrix< Scalar, M, N > &m, Eigen::Array< Real, MIN_(M, N), 1 > &s, Eigen::Matrix< Scalar, M, M > &u, Eigen::Matrix< Scalar, N, N > &vh, Real &s_errbd, Eigen::Array< Real, MIN_(M, N), 1 > &u_errbd, Eigen::Array< Real, MIN_(M, N), 1 > &v_errbd) |
| template<class Real , class Scalar , int M, int N> | |
| void | himalaya::reorder_svd (const Eigen::Matrix< Scalar, M, N > &m, Eigen::Array< Real, MIN_(M, N), 1 > &s) |
| template<class Real , class Scalar , int M, int N> | |
| void | himalaya::reorder_svd (const Eigen::Matrix< Scalar, M, N > &m, Eigen::Array< Real, MIN_(M, N), 1 > &s, Real &s_errbd) |
| template<class Real , int N> | |
| void | himalaya::reorder_diagonalize_symmetric_errbd (const Eigen::Matrix< std::complex< Real >, N, N > &m, Eigen::Array< Real, N, 1 > &s, Eigen::Matrix< std::complex< Real >, N, N > *u=0, Real *s_errbd=0, Eigen::Array< Real, N, 1 > *u_errbd=0) |
| template<class Real , int N> | |
| void | himalaya::reorder_diagonalize_symmetric_errbd (const Eigen::Matrix< Real, N, N > &m, Eigen::Array< Real, N, 1 > &s, Eigen::Matrix< std::complex< Real >, N, N > *u=0, Real *s_errbd=0, Eigen::Array< Real, N, 1 > *u_errbd=0) |
| template<class Real , class Scalar , int N> | |
| void | himalaya::reorder_diagonalize_symmetric (const Eigen::Matrix< Scalar, N, N > &m, Eigen::Array< Real, N, 1 > &s, Eigen::Matrix< std::complex< Real >, N, N > &u) |
| template<class Real , class Scalar , int N> | |
| void | himalaya::reorder_diagonalize_symmetric (const Eigen::Matrix< Scalar, N, N > &m, Eigen::Array< Real, N, 1 > &s, Eigen::Matrix< std::complex< Real >, N, N > &u, Real &s_errbd) |
| template<class Real , class Scalar , int N> | |
| void | himalaya::reorder_diagonalize_symmetric (const Eigen::Matrix< Scalar, N, N > &m, Eigen::Array< Real, N, 1 > &s, Eigen::Matrix< std::complex< Real >, N, N > &u, Real &s_errbd, Eigen::Array< Real, N, 1 > &u_errbd) |
| template<class Real , class Scalar , int N> | |
| void | himalaya::reorder_diagonalize_symmetric (const Eigen::Matrix< Scalar, N, N > &m, Eigen::Array< Real, N, 1 > &s) |
| template<class Real , class Scalar , int N> | |
| void | himalaya::reorder_diagonalize_symmetric (const Eigen::Matrix< Scalar, N, N > &m, Eigen::Array< Real, N, 1 > &s, Real &s_errbd) |
| template<class Real , class Scalar , int M, int N> | |
| void | himalaya::fs_svd_errbd (const Eigen::Matrix< Scalar, M, N > &m, Eigen::Array< Real, MIN_(M, N), 1 > &s, Eigen::Matrix< Scalar, M, M > *u=0, Eigen::Matrix< Scalar, N, N > *v=0, Real *s_errbd=0, Eigen::Array< Real, MIN_(M, N), 1 > *u_errbd=0, Eigen::Array< Real, MIN_(M, N), 1 > *v_errbd=0) |
| template<class Real , class Scalar , int M, int N> | |
| void | himalaya::fs_svd (const Eigen::Matrix< Scalar, M, N > &m, Eigen::Array< Real, MIN_(M, N), 1 > &s, Eigen::Matrix< Scalar, M, M > &u, Eigen::Matrix< Scalar, N, N > &v) |
| template<class Real , class Scalar , int M, int N> | |
| void | himalaya::fs_svd (const Eigen::Matrix< Scalar, M, N > &m, Eigen::Array< Real, MIN_(M, N), 1 > &s, Eigen::Matrix< Scalar, M, M > &u, Eigen::Matrix< Scalar, N, N > &v, Real &s_errbd) |
| template<class Real , class Scalar , int M, int N> | |
| void | himalaya::fs_svd (const Eigen::Matrix< Scalar, M, N > &m, Eigen::Array< Real, MIN_(M, N), 1 > &s, Eigen::Matrix< Scalar, M, M > &u, Eigen::Matrix< Scalar, N, N > &v, Real &s_errbd, Eigen::Array< Real, MIN_(M, N), 1 > &u_errbd, Eigen::Array< Real, MIN_(M, N), 1 > &v_errbd) |
| template<class Real , class Scalar , int M, int N> | |
| void | himalaya::fs_svd (const Eigen::Matrix< Scalar, M, N > &m, Eigen::Array< Real, MIN_(M, N), 1 > &s) |
| template<class Real , class Scalar , int M, int N> | |
| void | himalaya::fs_svd (const Eigen::Matrix< Scalar, M, N > &m, Eigen::Array< Real, MIN_(M, N), 1 > &s, Real &s_errbd) |
| template<class Real , int M, int N> | |
| void | himalaya::fs_svd (const Eigen::Matrix< Real, M, N > &m, Eigen::Array< Real, MIN_(M, N), 1 > &s, Eigen::Matrix< std::complex< Real >, M, M > &u, Eigen::Matrix< std::complex< Real >, N, N > &v) |
| template<class Real , int M, int N> | |
| void | himalaya::fs_svd (const Eigen::Matrix< Real, M, N > &m, Eigen::Array< Real, MIN_(M, N), 1 > &s, Eigen::Matrix< std::complex< Real >, M, M > &u, Eigen::Matrix< std::complex< Real >, N, N > &v, Real &s_errbd) |
| template<class Real , int M, int N> | |
| void | himalaya::fs_svd (const Eigen::Matrix< Real, M, N > &m, Eigen::Array< Real, MIN_(M, N), 1 > &s, Eigen::Matrix< std::complex< Real >, M, M > &u, Eigen::Matrix< std::complex< Real >, N, N > &v, Real &s_errbd, Eigen::Array< Real, MIN_(M, N), 1 > &u_errbd, Eigen::Array< Real, MIN_(M, N), 1 > &v_errbd) |
| template<class Real , class Scalar , int N> | |
| void | himalaya::fs_diagonalize_symmetric_errbd (const Eigen::Matrix< Scalar, N, N > &m, Eigen::Array< Real, N, 1 > &s, Eigen::Matrix< std::complex< Real >, N, N > *u=0, Real *s_errbd=0, Eigen::Array< Real, N, 1 > *u_errbd=0) |
| template<class Real , class Scalar , int N> | |
| void | himalaya::fs_diagonalize_symmetric (const Eigen::Matrix< Scalar, N, N > &m, Eigen::Array< Real, N, 1 > &s, Eigen::Matrix< std::complex< Real >, N, N > &u) |
| template<class Real , class Scalar , int N> | |
| void | himalaya::fs_diagonalize_symmetric (const Eigen::Matrix< Scalar, N, N > &m, Eigen::Array< Real, N, 1 > &s, Eigen::Matrix< std::complex< Real >, N, N > &u, Real &s_errbd) |
| template<class Real , class Scalar , int N> | |
| void | himalaya::fs_diagonalize_symmetric (const Eigen::Matrix< Scalar, N, N > &m, Eigen::Array< Real, N, 1 > &s, Eigen::Matrix< std::complex< Real >, N, N > &u, Real &s_errbd, Eigen::Array< Real, N, 1 > &u_errbd) |
| template<class Real , class Scalar , int N> | |
| void | himalaya::fs_diagonalize_symmetric (const Eigen::Matrix< Scalar, N, N > &m, Eigen::Array< Real, N, 1 > &s) |
| template<class Real , class Scalar , int N> | |
| void | himalaya::fs_diagonalize_symmetric (const Eigen::Matrix< Scalar, N, N > &m, Eigen::Array< Real, N, 1 > &s, Real &s_errbd) |
| template<class Real , class Scalar , int N> | |
| void | himalaya::fs_diagonalize_hermitian_errbd (const Eigen::Matrix< Scalar, N, N > &m, Eigen::Array< Real, N, 1 > &w, Eigen::Matrix< Scalar, N, N > *z=0, Real *w_errbd=0, Eigen::Array< Real, N, 1 > *z_errbd=0) |
| template<class Real , class Scalar , int N> | |
| void | himalaya::fs_diagonalize_hermitian (const Eigen::Matrix< Scalar, N, N > &m, Eigen::Array< Real, N, 1 > &w, Eigen::Matrix< Scalar, N, N > &z) |
| template<class Real , class Scalar , int N> | |
| void | himalaya::fs_diagonalize_hermitian (const Eigen::Matrix< Scalar, N, N > &m, Eigen::Array< Real, N, 1 > &w, Eigen::Matrix< Scalar, N, N > &z, Real &w_errbd) |
| template<class Real , class Scalar , int N> | |
| void | himalaya::fs_diagonalize_hermitian (const Eigen::Matrix< Scalar, N, N > &m, Eigen::Array< Real, N, 1 > &w, Eigen::Matrix< Scalar, N, N > &z, Real &w_errbd, Eigen::Array< Real, N, 1 > &z_errbd) |
| template<class Real , class Scalar , int N> | |
| void | himalaya::fs_diagonalize_hermitian (const Eigen::Matrix< Scalar, N, N > &m, Eigen::Array< Real, N, 1 > &w) |
| template<class Real , class Scalar , int N> | |
| void | himalaya::fs_diagonalize_hermitian (const Eigen::Matrix< Scalar, N, N > &m, Eigen::Array< Real, N, 1 > &w, Real &w_errbd) |
| std::tuple< V2, V2, V2, V2 > | himalaya::fs_diagonalize_hermitian_perturbatively (const RM22 &m0, const RM22 &m1, const RM22 &m2, const RM22 &m3) |
Detailed Description
Contains routines to diagonalize mass matrices.
- Note
- This file has been taken from the FlexibleSUSY 2.2.0 package.
Definition in file Linalg.hpp.
Macro Definition Documentation
◆ MAX_
| #define MAX_ | ( | i, | |
| j | |||
| ) | (((i) > (j)) ? (i) : (j)) |
Definition at line 41 of file Linalg.hpp.
◆ MIN_
| #define MIN_ | ( | i, | |
| j | |||
| ) | (((i) < (j)) ? (i) : (j)) |
Definition at line 42 of file Linalg.hpp.
