naginterfaces.library.matop.complex_herm_matrix_fun(uplo, a, f, data=None)[source]

complex_herm_matrix_fun computes the matrix function, , of a complex Hermitian matrix . must also be a complex Hermitian matrix.

For full information please refer to the NAG Library document for f01ff

uplostr, length 1

If , the upper triangle of the matrix is stored.

If , the lower triangle of the matrix is stored.

acomplex, array-like, shape

The Hermitian matrix .

fcallable fx = f(x, data=None)

The function evaluates at a number of points .

xfloat, ndarray, shape

The points at which the function is to be evaluated.

dataarbitrary, optional, modifiable in place

User-communication data for callback functions.

fxfloat, array-like, shape

The function values. should return the value , for .

dataarbitrary, optional

User-communication data for callback functions.

acomplex, ndarray, shape

If no exception or warning is raised, the upper or lower triangular part of the matrix function, .

(errno )

An internal error occurred when computing the spectral factorization. Please contact NAG.

(errno )

On entry, .

Constraint: .

(errno )

On entry, .

Constraint: or .

(errno )

The computation of the spectral factorization failed to converge.

(errno )

Termination requested in .


is computed using a spectral factorization of

where is the real diagonal matrix whose diagonal elements, , are the eigenvalues of , is a unitary matrix whose columns are the eigenvectors of and is the conjugate transpose of . is then given by

where is the diagonal matrix whose th diagonal element is . See for example Section 4.5 of Higham (2008). is assumed to be real.


Higham, N J, 2008, Functions of Matrices: Theory and Computation, SIAM, Philadelphia, PA, USA