#include "eigsol.h"
#include "vector.h"
#include <math.h>
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) |
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
x | - array containing eigenvectors (stored in columns) | |
w | - array containing eigenvalues |
input
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.
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.
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().