#include <math.h>
#include "edsolver.h"
#include "global.h"
#include "globmat.h"
#include "mechprint.h"
#include "elemswitch.h"
#include "gmatrix.h"
#include "eigsol.h"
#include "vector.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) |
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
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
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
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
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().