GNU Octave  4.0.0
A high-level interpreted language, primarily intended for numerical computations, mostly compatible with Matlab
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
Typedefs | Functions | Variables
eigs-base.cc File Reference
#include <cfloat>
#include <cmath>
#include <vector>
#include <iostream>
#include "f77-fcn.h"
#include "oct-locbuf.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)
 
static void warn_convergence (void)
 

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 51 of file eigs-base.cc.

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

Definition at line 49 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 3397 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 G__eigs__().

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 2740 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 G__eigs__().

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 3061 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 G__eigs__().

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 2427 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 G__eigs__().

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 1676 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 G__eigs__().

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 2045 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 G__eigs__().

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 1411 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 G__eigs__().

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 763 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 G__eigs__().

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 1084 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 G__eigs__().

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   
)

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

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 239 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

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

References abs(), sparse_base_lu< lu_type, lu_elt_type, p_type, p_elt_type >::col_perm(), d, 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(), warn_convergence(), Sparse< T >::xcidx(), Sparse< T >::xdata(), xisnan(), and Sparse< T >::xridx().

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

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 526 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(), warn_convergence(), 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(), warn_convergence(), 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 690 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(), warn_convergence(), 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 221 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 358 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 378 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 395 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 415 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 266 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 285 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 300 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 321 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 337 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 void warn_convergence ( void  )
static

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

Referenced by LuAminusSigmaB().

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 58 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 58 of file eigs-base.cc.

Generated on Thu Jun 4 2015 23:30:39 for GNU Octave by   doxygen 1.8.8