recexp_sprandsym.cc File Reference

Test serial recursive expansion and computation of homo and lumo eigenvectors of a sparse matrix with a given eigenspectrum. More...

#include "purification_sp2.h"
#include "purification_sp2acc.h"
#include "matrix_typedefs.h"
#include "realtype.h"
#include "matrix_utilities.h"
#include "integral_matrix_wrappers.h"
#include "SizesAndBlocks.h"
#include "Matrix.h"
#include "Vector.h"
#include "MatrixSymmetric.h"
#include "MatrixTriangular.h"
#include "MatrixGeneral.h"
#include "VectorGeneral.h"
#include "output.h"
#include <iostream>
#include <fstream>
#include <sstream>
#include <string.h>
#include "random_matrices.h"
#include "get_eigenvectors.h"

Macros
#define	SQRT_EPSILON_REAL   template_blas_sqrt(mat::getMachineEpsilon<real>())
#define	MAX_DOUBLE   std::numeric_limits<real>::max()
#define	MIN_DOUBLE   std::numeric_limits<real>::min()

Typedefs
typedef ergo_real	real
typedef symmMatrix	MatrixType
typedef normalMatrix	MatrixTypeGeneral
typedef MatrixType::VectorType	VectorType

Functions
static bool	abs_compare (real a, real b)
int	main (int argc, char *argv[])

Variables
real	TOL_ERR_SUBS_DEFAULT = SQRT_EPSILON_REAL
real	TOL_TRACE_ERROR_DEFAULT = SQRT_EPSILON_REAL
real	SCALAR_TOL = SQRT_EPSILON_REAL

Detailed Description

Test serial recursive expansion and computation of homo and lumo eigenvectors of a sparse matrix with a given eigenspectrum.

Matrix is generated using Givens rotations starting from a diagonal matrix with a required eigenvalues, so eigenvectors of a result matrix are explicitly known and saved into the matrix Q.

Author

Anastasia Kruchinina responsible

Macro Definition Documentation

#define MAX_DOUBLE   std::numeric_limits<real>::max()

Referenced by main().

#define MIN_DOUBLE   std::numeric_limits<real>::min()

Referenced by main().

#define SQRT_EPSILON_REAL   template_blas_sqrt(mat::getMachineEpsilon<real>())

Typedef Documentation

typedef symmMatrix MatrixType

typedef normalMatrix MatrixTypeGeneral

typedef ergo_real real

typedef MatrixType::VectorType VectorType

Function Documentation

static bool abs_compare ( real  a, real  b  )

static

References template_blas_fabs().

Referenced by main().

int main	(	int	argc,
		char *	argv[]
	)

References abs_compare(), PurificationGeneral< MatrixType >::clear(), PuriInfo::converged, distance(), mat::Interval< Treal >::empty(), mat::euclNorm, mat::VectorGeneral< Treal, Tvector >::fullvector(), mat::MatrixGeneral< Treal, Tmatrix >::get_values(), mat::Params::getMatrixParallelLevel(), mat::Params::getNProcs(), PurificationGeneral< MatrixType >::info, PurificationGeneral< MatrixType >::initialize(), MAX_DOUBLE, MIN_DOUBLE, mat::MatrixSymmetric< Treal, Tmatrix >::nnz(), PuriInfo::print_collected_info_printf(), PurificationGeneral< MatrixType >::PurificationStart(), read_vector(), SCALAR_TOL, PurificationGeneral< MatrixType >::set_eigenvectors_params(), mat::Params::setMatrixParallelLevel(), mat::Params::setNProcs(), sprandsym(), template_blas_fabs(), TOL_ERR_SUBS_DEFAULT, TOL_TRACE_ERROR_DEFAULT, mat::MatrixBase< Treal, Tmatrix >::trace(), write_vector(), and PurificationGeneral< MatrixType >::X.

Variable Documentation

real SCALAR_TOL = SQRT_EPSILON_REAL

real TOL_ERR_SUBS_DEFAULT = SQRT_EPSILON_REAL

real TOL_TRACE_ERROR_DEFAULT = SQRT_EPSILON_REAL