SIFEL/GEFEL/eigsol.h File Reference

Go to the source code of this file.

Functions
void	gen_jacobi (double *a, double *b, double *x, double *w, long n, long ni, double *thresholds, long nthr, double zero)
void	gen_jacobi2 (double *k, double *m, double *x, double *w, long n, long ni, double *gate, long ng, double limit)
void	jacobi (double *a, double *evec, double *eval, long n, long ni, long &ani, double limit, long normalize)

---

Function Documentation

void gen_jacobi	(	double *	a,
		double *	b,
		double *	x,
		double *	w,
		long	n,
		long	ni,
		double *	thresholds,
		long	nthr,
		double	zero
	)

generalized Jacobi iteration method for problems A.x = w.B.x

output

Parameters:

	x	- array containing eigenvectors (stored in columns)
	w	- array containing eigenvalues

input

Parameters:

	a	- matrix a (dense storage)
	b	- matrix b (dense storage)
	n	- number of rows (columns) of matrices a and b
	ni	- maximum number of iterations
	thresholds	- array containing thresholds
	nthr	- number of thresholds
	zero	- computer zero

20.1.2002

Definition at line 179 of file eigsol.cpp.

References limit.

Referenced by subspace_iter_jac().

void gen_jacobi2	(	double *	k,
		double *	m,
		double *	x,
		double *	w,
		long	n,
		long	ni,
		double *	gate,
		long	ng,
		double	limit
	)

Definition at line 306 of file eigsol.cpp.

References ll, nullv(), and ss().

void jacobi	(	double *	a,
		double *	evec,
		double *	eval,
		long	n,
		long	ni,
		long &	ani,
		double	limit,
		long	normalize
	)

function computes all eigenvalues and eigenvectors of the matrix matrix is stored as dense matrix and will be overwritten!

eigenvectors are columns of the matrix evec but matrix is stored in usual way, by rows.

Parameters:

	a	- array containing matrix
	evec	- array containing eigenvectors
	eval	- array containing eigenvalues
	n	- order of matrix a (number of rows or columns)
	ni	- maximum number of iterations
	ani	- number of performed iterations
	limit	- maximum acceptable absolute value of offdiagonal element
	normalize	- switch for normalizing of matrix (on>0, off=0)

!!!! nutno predelat trideni vlastnich vektoru

23.7.2001 updated 17.8.2007 by TKo. - normalizing of matrix

Definition at line 28 of file eigsol.cpp.

Referenced by princ_val().