#include <cfloat>
#include <cmath>
#include <vector>
#include <iostream>
#include "f77-fcn.h"
#include "quit.h"
#include "SparsedbleLU.h"
#include "SparseCmplxLU.h"
#include "dSparse.h"
#include "CSparse.h"
#include "MatrixType.h"
#include "SparsedbleCHOL.h"
#include "SparseCmplxCHOL.h"
#include "oct-rand.h"
#include "dbleCHOL.h"
#include "CmplxCHOL.h"
#include "dbleLU.h"
#include "CmplxLU.h"

Include dependency graph for eigs-base.cc:

This graph shows which files directly or indirectly include this file:

Go to the source code of this file.

Typedefs
typedef ComplexColumnVector(*	EigsComplexFunc )(const ComplexColumnVector &x, int &eigs_error)
typedef ColumnVector(*	EigsFunc )(const ColumnVector &x, int &eigs_error)

Functions
octave_idx_type	EigsComplexNonSymmetricFunc (EigsComplexFunc fun, octave_idx_type n, const std::string &_typ, Complex sigma, octave_idx_type k, octave_idx_type p, octave_idx_type &info, ComplexMatrix &eig_vec, ComplexColumnVector &eig_val, ComplexColumnVector &cresid, std::ostream &os, double tol, bool rvec, bool, int disp, int maxit)
template<class M >
octave_idx_type	EigsComplexNonSymmetricMatrix (const M &m, const std::string typ, octave_idx_type k, octave_idx_type p, octave_idx_type &info, ComplexMatrix &eig_vec, ComplexColumnVector &eig_val, const M &_b, ColumnVector &permB, ComplexColumnVector &cresid, std::ostream &os, double tol, bool rvec, bool cholB, int disp, int maxit)
octave_idx_type	EigsComplexNonSymmetricMatrix (const ComplexMatrix &m, const std::string typ, octave_idx_type k, octave_idx_type p, octave_idx_type &info, ComplexMatrix &eig_vec, ComplexColumnVector &eig_val, const ComplexMatrix &b, ColumnVector &permB, ComplexColumnVector &resid, std::ostream &os, double tol=std::numeric_limits< double >::epsilon(), bool rvec=false, bool cholB=0, int disp=0, int maxit=300)
octave_idx_type	EigsComplexNonSymmetricMatrix (const SparseComplexMatrix &m, const std::string typ, octave_idx_type k, octave_idx_type p, octave_idx_type &info, ComplexMatrix &eig_vec, ComplexColumnVector &eig_val, const SparseComplexMatrix &b, ColumnVector &permB, ComplexColumnVector &resid, std::ostream &os, double tol=std::numeric_limits< double >::epsilon(), bool rvec=false, bool cholB=0, int disp=0, int maxit=300)
template<class M >
octave_idx_type	EigsComplexNonSymmetricMatrixShift (const M &m, Complex sigma, octave_idx_type k, octave_idx_type p, octave_idx_type &info, ComplexMatrix &eig_vec, ComplexColumnVector &eig_val, const M &_b, ColumnVector &permB, ComplexColumnVector &cresid, std::ostream &os, double tol, bool rvec, bool cholB, int disp, int maxit)
octave_idx_type	EigsComplexNonSymmetricMatrixShift (const ComplexMatrix &m, Complex sigma, octave_idx_type k, octave_idx_type p, octave_idx_type &info, ComplexMatrix &eig_vec, ComplexColumnVector &eig_val, const ComplexMatrix &b, ColumnVector &permB, ComplexColumnVector &resid, std::ostream &os, double tol=std::numeric_limits< double >::epsilon(), bool rvec=false, bool cholB=0, int disp=0, int maxit=300)
octave_idx_type	EigsComplexNonSymmetricMatrixShift (const SparseComplexMatrix &m, Complex sigma, octave_idx_type k, octave_idx_type p, octave_idx_type &info, ComplexMatrix &eig_vec, ComplexColumnVector &eig_val, const SparseComplexMatrix &b, ColumnVector &permB, ComplexColumnVector &resid, std::ostream &os, double tol=std::numeric_limits< double >::epsilon(), bool rvec=false, bool cholB=0, int disp=0, int maxit=300)
octave_idx_type	EigsRealNonSymmetricFunc (EigsFunc fun, octave_idx_type n, const std::string &_typ, double sigmar, octave_idx_type k, octave_idx_type p, octave_idx_type &info, ComplexMatrix &eig_vec, ComplexColumnVector &eig_val, ColumnVector &resid, std::ostream &os, double tol, bool rvec, bool, int disp, int maxit)
template<class M >
octave_idx_type	EigsRealNonSymmetricMatrix (const M &m, const std::string typ, octave_idx_type k, octave_idx_type p, octave_idx_type &info, ComplexMatrix &eig_vec, ComplexColumnVector &eig_val, const M &_b, ColumnVector &permB, ColumnVector &resid, std::ostream &os, double tol, bool rvec, bool cholB, int disp, int maxit)
octave_idx_type	EigsRealNonSymmetricMatrix (const Matrix &m, const std::string typ, octave_idx_type k, octave_idx_type p, octave_idx_type &info, ComplexMatrix &eig_vec, ComplexColumnVector &eig_val, const Matrix &b, ColumnVector &permB, ColumnVector &resid, std::ostream &os, double tol=std::numeric_limits< double >::epsilon(), bool rvec=false, bool cholB=0, int disp=0, int maxit=300)
octave_idx_type	EigsRealNonSymmetricMatrix (const SparseMatrix &m, const std::string typ, octave_idx_type k, octave_idx_type p, octave_idx_type &info, ComplexMatrix &eig_vec, ComplexColumnVector &eig_val, const SparseMatrix &b, ColumnVector &permB, ColumnVector &resid, std::ostream &os, double tol=std::numeric_limits< double >::epsilon(), bool rvec=false, bool cholB=0, int disp=0, int maxit=300)
template<class M >
octave_idx_type	EigsRealNonSymmetricMatrixShift (const M &m, double sigmar, octave_idx_type k, octave_idx_type p, octave_idx_type &info, ComplexMatrix &eig_vec, ComplexColumnVector &eig_val, const M &_b, ColumnVector &permB, ColumnVector &resid, std::ostream &os, double tol, bool rvec, bool cholB, int disp, int maxit)
octave_idx_type	EigsRealNonSymmetricMatrixShift (const Matrix &m, double sigma, octave_idx_type k, octave_idx_type p, octave_idx_type &info, ComplexMatrix &eig_vec, ComplexColumnVector &eig_val, const Matrix &b, ColumnVector &permB, ColumnVector &resid, std::ostream &os, double tol=std::numeric_limits< double >::epsilon(), bool rvec=false, bool cholB=0, int disp=0, int maxit=300)
octave_idx_type	EigsRealNonSymmetricMatrixShift (const SparseMatrix &m, double sigma, octave_idx_type k, octave_idx_type p, octave_idx_type &info, ComplexMatrix &eig_vec, ComplexColumnVector &eig_val, const SparseMatrix &b, ColumnVector &permB, ColumnVector &resid, std::ostream &os, double tol=std::numeric_limits< double >::epsilon(), bool rvec=false, bool cholB=0, int disp=0, int maxit=300)
octave_idx_type	EigsRealSymmetricFunc (EigsFunc fun, octave_idx_type n, const std::string &_typ, double sigma, octave_idx_type k, octave_idx_type p, octave_idx_type &info, Matrix &eig_vec, ColumnVector &eig_val, ColumnVector &resid, std::ostream &os, double tol, bool rvec, bool, int disp, int maxit)
template<class M >
octave_idx_type	EigsRealSymmetricMatrix (const M &m, const std::string typ, octave_idx_type k, octave_idx_type p, octave_idx_type &info, Matrix &eig_vec, ColumnVector &eig_val, const M &_b, ColumnVector &permB, ColumnVector &resid, std::ostream &os, double tol, bool rvec, bool cholB, int disp, int maxit)
octave_idx_type	EigsRealSymmetricMatrix (const Matrix &m, const std::string typ, octave_idx_type k, octave_idx_type p, octave_idx_type &info, Matrix &eig_vec, ColumnVector &eig_val, const Matrix &b, ColumnVector &permB, ColumnVector &resid, std::ostream &os, double tol=std::numeric_limits< double >::epsilon(), bool rvec=false, bool cholB=0, int disp=0, int maxit=300)
octave_idx_type	EigsRealSymmetricMatrix (const SparseMatrix &m, const std::string typ, octave_idx_type k, octave_idx_type p, octave_idx_type &info, Matrix &eig_vec, ColumnVector &eig_val, const SparseMatrix &b, ColumnVector &permB, ColumnVector &resid, std::ostream &os, double tol=std::numeric_limits< double >::epsilon(), bool rvec=false, bool cholB=0, int disp=0, int maxit=300)
template<class M >
octave_idx_type	EigsRealSymmetricMatrixShift (const M &m, double sigma, octave_idx_type k, octave_idx_type p, octave_idx_type &info, Matrix &eig_vec, ColumnVector &eig_val, const M &_b, ColumnVector &permB, ColumnVector &resid, std::ostream &os, double tol, bool rvec, bool cholB, int disp, int maxit)
octave_idx_type	EigsRealSymmetricMatrixShift (const Matrix &m, double sigma, octave_idx_type k, octave_idx_type p, octave_idx_type &info, Matrix &eig_vec, ColumnVector &eig_val, const Matrix &b, ColumnVector &permB, ColumnVector &resid, std::ostream &os, double tol=std::numeric_limits< double >::epsilon(), bool rvec=false, bool cholB=0, int disp=0, int maxit=300)
octave_idx_type	EigsRealSymmetricMatrixShift (const SparseMatrix &m, double sigma, octave_idx_type k, octave_idx_type p, octave_idx_type &info, Matrix &eig_vec, ColumnVector &eig_val, const SparseMatrix &b, ColumnVector &permB, ColumnVector &resid, std::ostream &os, double tol=std::numeric_limits< double >::epsilon(), bool rvec=false, bool cholB=0, int disp=0, int maxit=300)
F77_RET_T	F77_FUNC (dsaupd, DSAUPD)(octave_idx_type &
F77_RET_T	F77_FUNC (dseupd, DSEUPD)(const octave_idx_type &
F77_RET_T	F77_FUNC (dnaupd, DNAUPD)(octave_idx_type &
F77_RET_T	F77_FUNC (dneupd, DNEUPD)(const octave_idx_type &
F77_RET_T	F77_FUNC (znaupd, ZNAUPD)(octave_idx_type &
F77_RET_T	F77_FUNC (zneupd, ZNEUPD)(const octave_idx_type &
F77_RET_T	F77_FUNC (dgemv, DGEMV)(F77_CONST_CHAR_ARG_DECL
F77_RET_T	F77_FUNC (zgemv, ZGEMV)(F77_CONST_CHAR_ARG_DECL
static ComplexMatrix	ltsolve (const SparseComplexMatrix &, const ColumnVector &, const ComplexMatrix &)
static Matrix	ltsolve (const SparseMatrix &, const ColumnVector &, const Matrix &,)
static ComplexMatrix	ltsolve (const ComplexMatrix &, const ColumnVector &, const ComplexMatrix &)
static Matrix	ltsolve (const Matrix &, const ColumnVector &, const Matrix &,)
template<class SM , class M >
static M	ltsolve (const SM &L, const ColumnVector &Q, const M &m)
static bool	LuAminusSigmaB (const SparseMatrix &m, const SparseMatrix &b, bool cholB, const ColumnVector &permB, double sigma, SparseMatrix &L, SparseMatrix &U, octave_idx_type P, octave_idx_type Q)
static bool	LuAminusSigmaB (const Matrix &m, const Matrix &b, bool cholB, const ColumnVector &permB, double sigma, Matrix &L, Matrix &U, octave_idx_type P, octave_idx_type Q)
static bool	LuAminusSigmaB (const SparseComplexMatrix &m, const SparseComplexMatrix &b, bool cholB, const ColumnVector &permB, Complex sigma, SparseComplexMatrix &L, SparseComplexMatrix &U, octave_idx_type P, octave_idx_type Q)
static bool	LuAminusSigmaB (const ComplexMatrix &m, const ComplexMatrix &b, bool cholB, const ColumnVector &permB, Complex sigma, ComplexMatrix &L, ComplexMatrix &U, octave_idx_type P, octave_idx_type Q)
static octave_idx_type	lusolve (const SparseMatrix &, const SparseMatrix &, Matrix &)
static octave_idx_type	lusolve (const SparseComplexMatrix &, const SparseComplexMatrix &, ComplexMatrix &)
static octave_idx_type	lusolve (const Matrix &, const Matrix &, Matrix &)
static octave_idx_type	lusolve (const ComplexMatrix &, const ComplexMatrix &, ComplexMatrix &)
template<class M , class SM >
static octave_idx_type	lusolve (const SM &L, const SM &U, M &m)
static bool	make_cholb (Matrix &b, Matrix &bt, ColumnVector &permB)
static bool	make_cholb (SparseMatrix &b, SparseMatrix &bt, ColumnVector &permB)
static bool	make_cholb (ComplexMatrix &b, ComplexMatrix &bt, ColumnVector &permB)
static bool	make_cholb (SparseComplexMatrix &b, SparseComplexMatrix &bt, ColumnVector &permB)
static ComplexMatrix	utsolve (const SparseComplexMatrix &, const ColumnVector &, const ComplexMatrix &)
static Matrix	utsolve (const SparseMatrix &, const ColumnVector &, const Matrix &)
static ComplexMatrix	utsolve (const ComplexMatrix &, const ColumnVector &, const ComplexMatrix &)
static Matrix	utsolve (const Matrix &, const ColumnVector &, const Matrix &)
template<class SM , class M >
static M	utsolve (const SM &U, const ColumnVector &Q, const M &m)
static bool	vector_product (const SparseMatrix &m, const double x, double y)
static bool	vector_product (const Matrix &m, const double x, double y)
static bool	vector_product (const SparseComplexMatrix &m, const Complex x, Complex y)
static bool	vector_product (const ComplexMatrix &m, const Complex x, Complex y)

Variables
F77_RET_T const octave_idx_type const octave_idx_type const double double const octave_idx_type double const octave_idx_type octave_idx_type octave_idx_type double double const octave_idx_type octave_idx_type &F77_CHAR_ARG_LEN_DECL	F77_CHAR_ARG_LEN_DECL
F77_RET_T	F77_CONST_CHAR_ARG_DECL

Typedef Documentation

typedef ComplexColumnVector(* EigsComplexFunc)(const ComplexColumnVector &x, int &eigs_error)

Definition at line 50 of file eigs-base.cc.

typedef ColumnVector(* EigsFunc)(const ColumnVector &x, int &eigs_error)

Definition at line 48 of file eigs-base.cc.

Function Documentation

octave_idx_type EigsComplexNonSymmetricFunc	(	EigsComplexFunc	fun,
		octave_idx_type	n,
		const std::string &	_typ,
		Complex	sigma,
		octave_idx_type	k,
		octave_idx_type	p,
		octave_idx_type &	info,
		ComplexMatrix &	eig_vec,
		ComplexColumnVector &	eig_val,
		ComplexColumnVector &	cresid,
		std::ostream &	os,
		double	tol,
		bool	rvec,
		bool	,
		int	disp,
		int	maxit
	)

Definition at line 3407 of file eigs-base.cc.

References abs(), d, octave_rand::distribution(), F77_CHAR_ARG_LEN, F77_CONST_CHAR_ARG2, f77_exception_encountered, F77_FUNC(), Array< T >::fortran_vec(), Array< T >::is_empty(), OCTAVE_LOCAL_BUFFER, octave_NaN, ComplexColumnVector::resize(), ComplexMatrix::resize(), octave_rand::vector(), x, and xisnan().

Referenced by DEFUN_DLD().

template<class M >

octave_idx_type EigsComplexNonSymmetricMatrix	(	const M &	m,
		const std::string	typ,
		octave_idx_type	k,
		octave_idx_type	p,
		octave_idx_type &	info,
		ComplexMatrix &	eig_vec,
		ComplexColumnVector &	eig_val,
		const M &	_b,
		ColumnVector &	permB,
		ComplexColumnVector &	cresid,
		std::ostream &	os,
		double	tol,
		bool	rvec,
		bool	cholB,
		int	disp,
		int	maxit
	)

Definition at line 2750 of file eigs-base.cc.

References d, D_NINT(), octave_rand::distribution(), F77_CHAR_ARG_LEN, F77_CONST_CHAR_ARG2, f77_exception_encountered, F77_FUNC(), Array< T >::fortran_vec(), Array< T >::is_empty(), Array< T >::length(), ltsolve(), M, make_cholb(), OCTAVE_LOCAL_BUFFER, octave_NaN, ComplexColumnVector::resize(), ComplexMatrix::resize(), utsolve(), octave_rand::vector(), vector_product(), and xisnan().

Referenced by DEFUN_DLD().

octave_idx_type EigsComplexNonSymmetricMatrix	(	const ComplexMatrix &	m,
		const std::string	typ,
		octave_idx_type	k,
		octave_idx_type	p,
		octave_idx_type &	info,
		ComplexMatrix &	eig_vec,
		ComplexColumnVector &	eig_val,
		const ComplexMatrix &	b,
		ColumnVector &	permB,
		ComplexColumnVector &	resid,
		std::ostream &	os,
		double	tol = `std::numeric_limits< double >::epsilon()`,
		bool	rvec = `false`,
		bool	cholB = `0`,
		int	disp = `0`,
		int	maxit = `300`
	)

octave_idx_type EigsComplexNonSymmetricMatrix	(	const SparseComplexMatrix &	m,
		const std::string	typ,
		octave_idx_type	k,
		octave_idx_type	p,
		octave_idx_type &	info,
		ComplexMatrix &	eig_vec,
		ComplexColumnVector &	eig_val,
		const SparseComplexMatrix &	b,
		ColumnVector &	permB,
		ComplexColumnVector &	resid,
		std::ostream &	os,
		double	tol = `std::numeric_limits< double >::epsilon()`,
		bool	rvec = `false`,
		bool	cholB = `0`,
		int	disp = `0`,
		int	maxit = `300`
	)

template<class M >

octave_idx_type EigsComplexNonSymmetricMatrixShift	(	const M &	m,
		Complex	sigma,
		octave_idx_type	k,
		octave_idx_type	p,
		octave_idx_type &	info,
		ComplexMatrix &	eig_vec,
		ComplexColumnVector &	eig_val,
		const M &	_b,
		ColumnVector &	permB,
		ComplexColumnVector &	cresid,
		std::ostream &	os,
		double	tol,
		bool	rvec,
		bool	cholB,
		int	disp,
		int	maxit
	)

Definition at line 3071 of file eigs-base.cc.

References d, D_NINT(), octave_rand::distribution(), F77_CHAR_ARG_LEN, F77_CONST_CHAR_ARG2, f77_exception_encountered, F77_FUNC(), Array< T >::fortran_vec(), Array< T >::is_empty(), Array< T >::length(), LuAminusSigmaB(), lusolve(), M, OCTAVE_LOCAL_BUFFER, octave_NaN, Q, ComplexColumnVector::resize(), ComplexMatrix::resize(), octave_rand::vector(), vector_product(), and xisnan().

Referenced by DEFUN_DLD().

octave_idx_type EigsComplexNonSymmetricMatrixShift	(	const ComplexMatrix &	m,
		Complex	sigma,
		octave_idx_type	k,
		octave_idx_type	p,
		octave_idx_type &	info,
		ComplexMatrix &	eig_vec,
		ComplexColumnVector &	eig_val,
		const ComplexMatrix &	b,
		ColumnVector &	permB,
		ComplexColumnVector &	resid,
		std::ostream &	os,
		double	tol = `std::numeric_limits< double >::epsilon()`,
		bool	rvec = `false`,
		bool	cholB = `0`,
		int	disp = `0`,
		int	maxit = `300`
	)

octave_idx_type EigsComplexNonSymmetricMatrixShift	(	const SparseComplexMatrix &	m,
		Complex	sigma,
		octave_idx_type	k,
		octave_idx_type	p,
		octave_idx_type &	info,
		ComplexMatrix &	eig_vec,
		ComplexColumnVector &	eig_val,
		const SparseComplexMatrix &	b,
		ColumnVector &	permB,
		ComplexColumnVector &	resid,
		std::ostream &	os,
		double	tol = `std::numeric_limits< double >::epsilon()`,
		bool	rvec = `false`,
		bool	cholB = `0`,
		int	disp = `0`,
		int	maxit = `300`
	)

octave_idx_type EigsRealNonSymmetricFunc	(	EigsFunc	fun,
		octave_idx_type	n,
		const std::string &	_typ,
		double	sigmar,
		octave_idx_type	k,
		octave_idx_type	p,
		octave_idx_type &	info,
		ComplexMatrix &	eig_vec,
		ComplexColumnVector &	eig_val,
		ColumnVector &	resid,
		std::ostream &	os,
		double	tol,
		bool	rvec,
		bool	,
		int	disp,
		int	maxit
	)

Definition at line 2437 of file eigs-base.cc.

References abs(), d, octave_rand::distribution(), F77_CHAR_ARG_LEN, F77_CONST_CHAR_ARG2, f77_exception_encountered, F77_FUNC(), Array< T >::fortran_vec(), imag(), Array< T >::is_empty(), OCTAVE_LOCAL_BUFFER, octave_NaN, ComplexColumnVector::resize(), ComplexMatrix::resize(), octave_rand::vector(), x, and xisnan().

Referenced by DEFUN_DLD().

template<class M >

octave_idx_type EigsRealNonSymmetricMatrix	(	const M &	m,
		const std::string	typ,
		octave_idx_type	k,
		octave_idx_type	p,
		octave_idx_type &	info,
		ComplexMatrix &	eig_vec,
		ComplexColumnVector &	eig_val,
		const M &	_b,
		ColumnVector &	permB,
		ColumnVector &	resid,
		std::ostream &	os,
		double	tol,
		bool	rvec,
		bool	cholB,
		int	disp,
		int	maxit
	)

Definition at line 1686 of file eigs-base.cc.

References d, D_NINT(), octave_rand::distribution(), F77_CHAR_ARG_LEN, F77_CONST_CHAR_ARG2, f77_exception_encountered, F77_FUNC(), Array< T >::fortran_vec(), imag(), Array< T >::is_empty(), Array< T >::length(), ltsolve(), M, make_cholb(), OCTAVE_LOCAL_BUFFER, octave_NaN, ComplexColumnVector::resize(), ComplexMatrix::resize(), utsolve(), octave_rand::vector(), vector_product(), and xisnan().

Referenced by DEFUN_DLD().

octave_idx_type EigsRealNonSymmetricMatrix	(	const Matrix &	m,
		const std::string	typ,
		octave_idx_type	k,
		octave_idx_type	p,
		octave_idx_type &	info,
		ComplexMatrix &	eig_vec,
		ComplexColumnVector &	eig_val,
		const Matrix &	b,
		ColumnVector &	permB,
		ColumnVector &	resid,
		std::ostream &	os,
		double	tol = `std::numeric_limits< double >::epsilon()`,
		bool	rvec = `false`,
		bool	cholB = `0`,
		int	disp = `0`,
		int	maxit = `300`
	)

octave_idx_type EigsRealNonSymmetricMatrix	(	const SparseMatrix &	m,
		const std::string	typ,
		octave_idx_type	k,
		octave_idx_type	p,
		octave_idx_type &	info,
		ComplexMatrix &	eig_vec,
		ComplexColumnVector &	eig_val,
		const SparseMatrix &	b,
		ColumnVector &	permB,
		ColumnVector &	resid,
		std::ostream &	os,
		double	tol = `std::numeric_limits< double >::epsilon()`,
		bool	rvec = `false`,
		bool	cholB = `0`,
		int	disp = `0`,
		int	maxit = `300`
	)

template<class M >

octave_idx_type EigsRealNonSymmetricMatrixShift	(	const M &	m,
		double	sigmar,
		octave_idx_type	k,
		octave_idx_type	p,
		octave_idx_type &	info,
		ComplexMatrix &	eig_vec,
		ComplexColumnVector &	eig_val,
		const M &	_b,
		ColumnVector &	permB,
		ColumnVector &	resid,
		std::ostream &	os,
		double	tol,
		bool	rvec,
		bool	cholB,
		int	disp,
		int	maxit
	)

Definition at line 2055 of file eigs-base.cc.

References d, D_NINT(), octave_rand::distribution(), F77_CHAR_ARG_LEN, F77_CONST_CHAR_ARG2, f77_exception_encountered, F77_FUNC(), Array< T >::fortran_vec(), imag(), Array< T >::is_empty(), Array< T >::length(), LuAminusSigmaB(), lusolve(), M, OCTAVE_LOCAL_BUFFER, octave_NaN, Q, ComplexColumnVector::resize(), ComplexMatrix::resize(), octave_rand::vector(), vector_product(), and xisnan().

Referenced by DEFUN_DLD().

octave_idx_type EigsRealNonSymmetricMatrixShift	(	const Matrix &	m,
		double	sigma,
		octave_idx_type	k,
		octave_idx_type	p,
		octave_idx_type &	info,
		ComplexMatrix &	eig_vec,
		ComplexColumnVector &	eig_val,
		const Matrix &	b,
		ColumnVector &	permB,
		ColumnVector &	resid,
		std::ostream &	os,
		double	tol = `std::numeric_limits< double >::epsilon()`,
		bool	rvec = `false`,
		bool	cholB = `0`,
		int	disp = `0`,
		int	maxit = `300`
	)

octave_idx_type EigsRealNonSymmetricMatrixShift	(	const SparseMatrix &	m,
		double	sigma,
		octave_idx_type	k,
		octave_idx_type	p,
		octave_idx_type &	info,
		ComplexMatrix &	eig_vec,
		ComplexColumnVector &	eig_val,
		const SparseMatrix &	b,
		ColumnVector &	permB,
		ColumnVector &	resid,
		std::ostream &	os,
		double	tol = `std::numeric_limits< double >::epsilon()`,
		bool	rvec = `false`,
		bool	cholB = `0`,
		int	disp = `0`,
		int	maxit = `300`
	)

octave_idx_type EigsRealSymmetricFunc	(	EigsFunc	fun,
		octave_idx_type	n,
		const std::string &	_typ,
		double	sigma,
		octave_idx_type	k,
		octave_idx_type	p,
		octave_idx_type &	info,
		Matrix &	eig_vec,
		ColumnVector &	eig_val,
		ColumnVector &	resid,
		std::ostream &	os,
		double	tol,
		bool	rvec,
		bool	,
		int	disp,
		int	maxit
	)

Definition at line 1421 of file eigs-base.cc.

References abs(), d, octave_rand::distribution(), F77_CHAR_ARG_LEN, F77_CONST_CHAR_ARG2, f77_exception_encountered, F77_FUNC(), Array< T >::fortran_vec(), Array< T >::is_empty(), OCTAVE_LOCAL_BUFFER, octave_NaN, ColumnVector::resize(), Matrix::resize(), octave_rand::vector(), x, and xisnan().

Referenced by DEFUN_DLD().

template<class M >

octave_idx_type EigsRealSymmetricMatrix	(	const M &	m,
		const std::string	typ,
		octave_idx_type	k,
		octave_idx_type	p,
		octave_idx_type &	info,
		Matrix &	eig_vec,
		ColumnVector &	eig_val,
		const M &	_b,
		ColumnVector &	permB,
		ColumnVector &	resid,
		std::ostream &	os,
		double	tol,
		bool	rvec,
		bool	cholB,
		int	disp,
		int	maxit
	)

Definition at line 773 of file eigs-base.cc.

References d, D_NINT(), octave_rand::distribution(), F77_CHAR_ARG_LEN, F77_CONST_CHAR_ARG2, f77_exception_encountered, F77_FUNC(), Array< T >::fortran_vec(), Array< T >::is_empty(), Array< T >::length(), ltsolve(), M, make_cholb(), OCTAVE_LOCAL_BUFFER, octave_NaN, ColumnVector::resize(), Matrix::resize(), utsolve(), octave_rand::vector(), vector_product(), and xisnan().

Referenced by DEFUN_DLD().

octave_idx_type EigsRealSymmetricMatrix	(	const Matrix &	m,
		const std::string	typ,
		octave_idx_type	k,
		octave_idx_type	p,
		octave_idx_type &	info,
		Matrix &	eig_vec,
		ColumnVector &	eig_val,
		const Matrix &	b,
		ColumnVector &	permB,
		ColumnVector &	resid,
		std::ostream &	os,
		double	tol = `std::numeric_limits< double >::epsilon()`,
		bool	rvec = `false`,
		bool	cholB = `0`,
		int	disp = `0`,
		int	maxit = `300`
	)

octave_idx_type EigsRealSymmetricMatrix	(	const SparseMatrix &	m,
		const std::string	typ,
		octave_idx_type	k,
		octave_idx_type	p,
		octave_idx_type &	info,
		Matrix &	eig_vec,
		ColumnVector &	eig_val,
		const SparseMatrix &	b,
		ColumnVector &	permB,
		ColumnVector &	resid,
		std::ostream &	os,
		double	tol = `std::numeric_limits< double >::epsilon()`,
		bool	rvec = `false`,
		bool	cholB = `0`,
		int	disp = `0`,
		int	maxit = `300`
	)

template<class M >

octave_idx_type EigsRealSymmetricMatrixShift	(	const M &	m,
		double	sigma,
		octave_idx_type	k,
		octave_idx_type	p,
		octave_idx_type &	info,
		Matrix &	eig_vec,
		ColumnVector &	eig_val,
		const M &	_b,
		ColumnVector &	permB,
		ColumnVector &	resid,
		std::ostream &	os,
		double	tol,
		bool	rvec,
		bool	cholB,
		int	disp,
		int	maxit
	)

Definition at line 1094 of file eigs-base.cc.

References d, D_NINT(), octave_rand::distribution(), F77_CHAR_ARG_LEN, F77_CONST_CHAR_ARG2, f77_exception_encountered, F77_FUNC(), Array< T >::fortran_vec(), Array< T >::is_empty(), Array< T >::length(), LuAminusSigmaB(), lusolve(), M, OCTAVE_LOCAL_BUFFER, octave_NaN, Q, ColumnVector::resize(), Matrix::resize(), octave_rand::vector(), vector_product(), and xisnan().

Referenced by DEFUN_DLD().

octave_idx_type EigsRealSymmetricMatrixShift	(	const Matrix &	m,
		double	sigma,
		octave_idx_type	k,
		octave_idx_type	p,
		octave_idx_type &	info,
		Matrix &	eig_vec,
		ColumnVector &	eig_val,
		const Matrix &	b,
		ColumnVector &	permB,
		ColumnVector &	resid,
		std::ostream &	os,
		double	tol = `std::numeric_limits< double >::epsilon()`,
		bool	rvec = `false`,
		bool	cholB = `0`,
		int	disp = `0`,
		int	maxit = `300`
	)

octave_idx_type EigsRealSymmetricMatrixShift	(	const SparseMatrix &	m,
		double	sigma,
		octave_idx_type	k,
		octave_idx_type	p,
		octave_idx_type &	info,
		Matrix &	eig_vec,
		ColumnVector &	eig_val,
		const SparseMatrix &	b,
		ColumnVector &	permB,
		ColumnVector &	resid,
		std::ostream &	os,
		double	tol = `std::numeric_limits< double >::epsilon()`,
		bool	rvec = `false`,
		bool	cholB = `0`,
		int	disp = `0`,
		int	maxit = `300`
	)

F77_RET_T F77_FUNC	(	dsaupd	,
		DSAUPD
	)

F77_RET_T F77_FUNC	(	dseupd	,
		DSEUPD
	)		const

F77_RET_T F77_FUNC	(	dnaupd	,
		DNAUPD
	)

F77_RET_T F77_FUNC	(	dneupd	,
		DNEUPD
	)		const

F77_RET_T F77_FUNC	(	znaupd	,
		ZNAUPD
	)

F77_RET_T F77_FUNC	(	zneupd	,
		ZNEUPD
	)		const

F77_RET_T F77_FUNC	(	dgemv	,
		DGEMV
	)

F77_RET_T F77_FUNC	(	zgemv	,
		ZGEMV
	)

template ComplexMatrix ltsolve	(	const SparseComplexMatrix &	,
		const ColumnVector &	,
		const ComplexMatrix &
	)

static

Referenced by EigsComplexNonSymmetricMatrix(), EigsRealNonSymmetricMatrix(), and EigsRealSymmetricMatrix().

template Matrix ltsolve	(	const SparseMatrix &	,
		const ColumnVector &	,
		const Matrix &
	)

static

template ComplexMatrix ltsolve	(	const ComplexMatrix &	,
		const ColumnVector &	,
		const ComplexMatrix &
	)

static

template Matrix ltsolve	(	const Matrix &	,
		const ColumnVector &	,
		const Matrix &
	)

static

template<class SM , class M >

static M ltsolve	(	const SM &	L,
		const ColumnVector &	Q,
		const M &	m
	)

static

Definition at line 229 of file eigs-base.cc.

References Array< T >::fortran_vec(), MatrixType::Lower, M, and ComplexMatrix::rcond().

static bool LuAminusSigmaB	(	const SparseMatrix &	m,
		const SparseMatrix &	b,
		bool	cholB,
		const ColumnVector &	permB,
		double	sigma,
		SparseMatrix &	L,
		SparseMatrix &	U,
		octave_idx_type *	P,
		octave_idx_type *	Q
	)

static

static bool LuAminusSigmaB	(	const Matrix &	m,
		const Matrix &	b,
		bool	cholB,
		const ColumnVector &	permB,
		double	sigma,
		Matrix &	L,
		Matrix &	U,
		octave_idx_type *	P,
		octave_idx_type *	Q
	)

static

Definition at line 521 of file eigs-base.cc.

References abs(), Array< T >::cols(), d, Array< T >::fortran_vec(), Array< T >::is_empty(), base_lu< lu_type >::L(), Array< T >::length(), octave_NaN, base_lu< lu_type >::P(), ComplexMatrix::rcond(), Array< T >::rows(), PermMatrix::transpose(), Matrix::transpose(), base_lu< lu_type >::U(), Array< T >::xelem(), and xisnan().

static bool LuAminusSigmaB	(	const SparseComplexMatrix &	m,
		const SparseComplexMatrix &	b,
		bool	cholB,
		const ColumnVector &	permB,
		Complex	sigma,
		SparseComplexMatrix &	L,
		SparseComplexMatrix &	U,
		octave_idx_type *	P,
		octave_idx_type *	Q
	)

static

Definition at line 598 of file eigs-base.cc.

References abs(), sparse_base_lu< lu_type, lu_elt_type, p_type, p_elt_type >::col_perm(), d, SparseComplexMatrix::hermitian(), Sparse< T >::is_empty(), sparse_base_lu< lu_type, lu_elt_type, p_type, p_elt_type >::L(), Array< T >::length(), octave_NaN, ComplexMatrix::rcond(), sparse_base_lu< lu_type, lu_elt_type, p_type, p_elt_type >::row_perm(), Sparse< T >::rows(), SparseMatrix::transpose(), sparse_base_lu< lu_type, lu_elt_type, p_type, p_elt_type >::U(), Sparse< T >::xcidx(), Sparse< T >::xdata(), xisnan(), and Sparse< T >::xridx().

static bool LuAminusSigmaB	(	const ComplexMatrix &	m,
		const ComplexMatrix &	b,
		bool	cholB,
		const ColumnVector &	permB,
		Complex	sigma,
		ComplexMatrix &	L,
		ComplexMatrix &	U,
		octave_idx_type *	P,
		octave_idx_type *	Q
	)

static

Definition at line 695 of file eigs-base.cc.

References abs(), Array< T >::cols(), d, Array< T >::fortran_vec(), ComplexMatrix::hermitian(), Array< T >::is_empty(), base_lu< lu_type >::L(), Array< T >::length(), octave_NaN, base_lu< lu_type >::P(), ComplexMatrix::rcond(), Array< T >::rows(), PermMatrix::transpose(), base_lu< lu_type >::U(), Array< T >::xelem(), and xisnan().

template octave_idx_type lusolve	(	const SparseMatrix &	,
		const SparseMatrix &	,
		Matrix &
	)

static

Referenced by EigsComplexNonSymmetricMatrixShift(), EigsRealNonSymmetricMatrixShift(), and EigsRealSymmetricMatrixShift().

template octave_idx_type lusolve	(	const SparseComplexMatrix &	,
		const SparseComplexMatrix &	,
		ComplexMatrix &
	)

static

template octave_idx_type lusolve	(	const Matrix &	,
		const Matrix &	,
		Matrix &
	)

static

template octave_idx_type lusolve	(	const ComplexMatrix &	,
		const ComplexMatrix &	,
		ComplexMatrix &
	)

static

template<class M , class SM >

static octave_idx_type lusolve	(	const SM &	L,
		const SM &	U,
		M &	m
	)

static

Definition at line 211 of file eigs-base.cc.

References ComplexMatrix::rcond(), and MatrixType::Upper.

static bool make_cholb	(	Matrix &	b,
		Matrix &	bt,
		ColumnVector &	permB
	)

static

Definition at line 348 of file eigs-base.cc.

References CHOL::chol_matrix(), Array< T >::cols(), and Matrix::transpose().

Referenced by EigsComplexNonSymmetricMatrix(), EigsRealNonSymmetricMatrix(), and EigsRealSymmetricMatrix().

static bool make_cholb	(	SparseMatrix &	b,
		SparseMatrix &	bt,
		ColumnVector &	permB
	)

static

Definition at line 368 of file eigs-base.cc.

References SparseCHOL::L(), SparseCHOL::P(), SparseCHOL::perm(), and SparseMatrix::transpose().

static bool make_cholb	(	ComplexMatrix &	b,
		ComplexMatrix &	bt,
		ColumnVector &	permB
	)

static

Definition at line 385 of file eigs-base.cc.

References ComplexCHOL::chol_matrix(), Array< T >::cols(), and ComplexMatrix::hermitian().

static bool make_cholb	(	SparseComplexMatrix &	b,
		SparseComplexMatrix &	bt,
		ColumnVector &	permB
	)

static

Definition at line 405 of file eigs-base.cc.

References SparseComplexMatrix::hermitian(), SparseComplexCHOL::L(), SparseComplexCHOL::P(), and SparseComplexCHOL::perm().

template ComplexMatrix utsolve	(	const SparseComplexMatrix &	,
		const ColumnVector &	,
		const ComplexMatrix &
	)

static

Referenced by EigsComplexNonSymmetricMatrix(), EigsRealNonSymmetricMatrix(), and EigsRealSymmetricMatrix().

template Matrix utsolve	(	const SparseMatrix &	,
		const ColumnVector &	,
		const Matrix &
	)

static

template ComplexMatrix utsolve	(	const ComplexMatrix &	,
		const ColumnVector &	,
		const ComplexMatrix &
	)

static

template Matrix utsolve	(	const Matrix &	,
		const ColumnVector &	,
		const Matrix &
	)

static

template<class SM , class M >

static M utsolve	(	const SM &	U,
		const ColumnVector &	Q,
		const M &	m
	)

static

Definition at line 256 of file eigs-base.cc.

References Array< T >::fortran_vec(), M, ComplexMatrix::rcond(), and MatrixType::Upper.

static bool vector_product	(	const SparseMatrix &	m,
		const double *	x,
		double *	y
	)

static

Definition at line 275 of file eigs-base.cc.

References Sparse< T >::cidx(), Sparse< T >::cols(), Sparse< T >::data(), and Sparse< T >::ridx().

Referenced by EigsComplexNonSymmetricMatrix(), EigsComplexNonSymmetricMatrixShift(), EigsRealNonSymmetricMatrix(), EigsRealNonSymmetricMatrixShift(), EigsRealSymmetricMatrix(), and EigsRealSymmetricMatrixShift().

static bool vector_product	(	const Matrix &	m,
		const double *	x,
		double *	y
	)

static

Definition at line 290 of file eigs-base.cc.

References Array< T >::cols(), Array< T >::data(), F77_CHAR_ARG_LEN, F77_CONST_CHAR_ARG2, f77_exception_encountered, F77_XFCN, Array< T >::rows(), and x.

static bool vector_product	(	const SparseComplexMatrix &	m,
		const Complex *	x,
		Complex *	y
	)

static

Definition at line 311 of file eigs-base.cc.

References Sparse< T >::cidx(), Sparse< T >::cols(), Sparse< T >::data(), and Sparse< T >::ridx().

static bool vector_product	(	const ComplexMatrix &	m,
		const Complex *	x,
		Complex *	y
	)

static

Definition at line 327 of file eigs-base.cc.

References Array< T >::cols(), Array< T >::data(), F77_CHAR_ARG_LEN, F77_CONST_CHAR_ARG2, f77_exception_encountered, F77_XFCN, Array< T >::rows(), and x.

Variable Documentation

F77_RET_T const octave_idx_type const octave_idx_type const Complex const Complex const octave_idx_type const Complex const octave_idx_type const Complex Complex const octave_idx_type& F77_CHAR_ARG_LEN_DECL

Definition at line 57 of file eigs-base.cc.

F77_RET_T octave_idx_type Complex Complex const octave_idx_type const Complex Complex const octave_idx_type F77_CONST_CHAR_ARG_DECL

Definition at line 57 of file eigs-base.cc.

Typedefs

Functions

Variables

Typedef Documentation

Function Documentation

Variable Documentation