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 
)
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   
)
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
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
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
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
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
template octave_idx_type lusolve ( const SparseMatrix ,
const SparseMatrix ,
Matrix  
)
static
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
static bool make_cholb ( SparseMatrix b,
SparseMatrix bt,
ColumnVector permB 
)
static
static bool make_cholb ( ComplexMatrix b,
ComplexMatrix bt,
ColumnVector permB 
)
static
static bool make_cholb ( SparseComplexMatrix b,
SparseComplexMatrix bt,
ColumnVector permB 
)
static
template ComplexMatrix utsolve ( const SparseComplexMatrix ,
const ColumnVector ,
const ComplexMatrix  
)
static
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
static bool vector_product ( const Matrix m,
const double x,
double y 
)
static
static bool vector_product ( const SparseComplexMatrix m,
const Complex x,
Complex y 
)
static
static bool vector_product ( const ComplexMatrix m,
const Complex x,
Complex y 
)
static
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.

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