SIFEL/MEFEL/SRC/edsolver.h File Reference

#include <stdio.h>

Go to the source code of this file.

Functions
void	inverse_iteration (double *x, double *w)
void	solve_eigen_dynamics (double *x, double *w)
void	subspace_iter_jac (double *x, double *w)
void	subspace_iter_ortho (double *x, double *w)
void	subspace_shift_iter_ortho (double *x, double *w)

---

Function Documentation

void inverse_iteration	(	double *	x,
		double *	w
	)

function computes the smallest eigenvalue and appropriate eigenvector of the generalized problem Ax=wBx the eigenvector of the smallest eigenvalue is stored in Lsrs in lhs

Parameters:

w

- the smallest eigenvalue

3.11.2001

Definition at line 73 of file edsolver.cpp.

References eigvalsol::aerres, eigvalsol::anies, cmulv(), copyv(), probdesc::eigsol, eigvalsol::erres, gmatrix::gmxv(), Gtm, Mespr, Mmat, Mp, Ndofm, eigvalsol::nies, Out, p, Smat, slesolv::solve_system(), ss(), probdesc::ssle, and probdesc::zero.

Referenced by parallel_solution_eigen_dynamics(), and solve_eigen_dynamics().

void solve_eigen_dynamics	(	double *	x,
		double *	w
	)

Definition at line 13 of file edsolver.cpp.

References compute_req_val(), probdesc::eigsol, inv_iteration, inverse_iteration(), mass_matrix(), Mp, eigvalsol::neigv, print_close(), print_eigenvalues(), print_eigenvectors(), print_err(), print_init(), print_step(), shifted_subspace_it_gsortho, stiffness_matrix(), subspace_it_gsortho, subspace_it_jacobi, subspace_iter_jac(), subspace_iter_ortho(), subspace_shift_iter_ortho(), and eigvalsol::teigsol.

Referenced by response_spectrum_method(), and solve_mefel_deterministic_problem().

void subspace_iter_jac	(	double *	x,
		double *	w
	)

function is not ready for computation initial vectors must be assembled

function computes eigenvalues and eigenvectors subspace iteration with Jacobi method of rotations is used

Parameters:

w

- array containing eigenvalues

4.11.2001

Definition at line 316 of file edsolver.cpp.

References eigvalsol::aerres, eigvalsol::anies, cmulv(), copyv(), probdesc::eigsol, eigvalsol::erres, gen_jacobi(), gmatrix::gmxv(), Gtm, eigvalsol::jacthr, Mespr, Mmat, Mp, mtxmccr(), Ndofm, eigvalsol::neigv, eigvalsol::nev, eigvalsol::nies, eigvalsol::nijmr, eigvalsol::njacthr, Out, p, Smat, slesolv::solve_system(), ss(), probdesc::ssle, and probdesc::zero.

Referenced by parallel_solution_eigen_dynamics(), and solve_eigen_dynamics().

void subspace_iter_ortho	(	double *	x,
		double *	w
	)

function computes eigenvalues and eigenvectors subspace iteration with Gram-Schmidt orthonormalization is used

Parameters:

w

- array containing eigenvalues

4.11.2001

Definition at line 151 of file edsolver.cpp.

References eigvalsol::aerres, eigvalsol::anies, cmulv(), copyv(), probdesc::eigsol, eigvalsol::erres, gmatrix::gmxv(), Gtm, Mespr, Mmat, Mp, Ndofm, eigvalsol::neigv, eigvalsol::nev, eigvalsol::nies, Out, p, Smat, slesolv::solve_system(), ss(), probdesc::ssle, and probdesc::zero.

Referenced by parallel_solution_eigen_dynamics(), and solve_eigen_dynamics().

void subspace_shift_iter_ortho	(	double *	x,
		double *	w
	)

function computes eigenvalues and eigenvectors subspace iteration with Gram-Schmidt orthonormalization is used shift is applied in order to solve systems with singular stiffness matrix

Parameters:

w

- array containing eigenvalues

6. 2. 2014

Definition at line 527 of file edsolver.cpp.

References gmatrix::addgm(), eigvalsol::aerres, eigvalsol::anies, cmulv(), copyv(), probdesc::eigsol, eigvalsol::erres, gmatrix::gmxv(), Gtm, Mespr, Mmat, Mp, Ndofm, eigvalsol::neigv, eigvalsol::nev, eigvalsol::nies, Out, p, eigvalsol::shift, Smat, slesolv::solve_system(), ss(), probdesc::ssle, and probdesc::zero.

Referenced by solve_eigen_dynamics().