GNU Octave  8.1.0
A high-level interpreted language, primarily intended for numerical computations, mostly compatible with Matlab
Functions
eigs-base.cc File Reference
#include <cmath>
#include <ostream>
#include "Array.h"
#include "CSparse.h"
#include "MatrixType.h"
#include "PermMatrix.h"
#include "chol.h"
#include "dSparse.h"
#include "eigs-base.h"
#include "lo-arpack-proto.h"
#include "lo-blas-proto.h"
#include "lo-error.h"
#include "lo-ieee.h"
#include "lu.h"
#include "mx-ops.h"
#include "oct-locbuf.h"
#include "oct-rand.h"
#include "sparse-chol.h"
#include "sparse-lu.h"
Include dependency graph for eigs-base.cc:

Go to the source code of this file.

Functions

std::string arpack_errno2str (const octave_idx_type &errnum, const std::string &fcn_name)
 
template<typename M >
octave_idx_type EigsComplexNonSymmetricFunc (EigsComplexFunc fcn, octave_idx_type n_arg, const std::string &_typ, Complex sigma, octave_idx_type k_arg, octave_idx_type p_arg, 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)
 
template OCTAVE_API octave_idx_type EigsComplexNonSymmetricFunc< ComplexMatrix > (EigsComplexFunc fcn, 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, const ComplexMatrix &_b, ColumnVector &permB, ComplexColumnVector &cresid, std::ostream &os, double tol, bool rvec, bool cholB, int disp, int maxit)
 
template OCTAVE_API octave_idx_type EigsComplexNonSymmetricFunc< SparseComplexMatrix > (EigsComplexFunc fcn, 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, const SparseComplexMatrix &_b, ColumnVector &permB, ComplexColumnVector &cresid, std::ostream &os, double tol, bool rvec, bool cholB, int disp, int maxit)
 
template<typename M >
octave_idx_type EigsComplexNonSymmetricMatrix (const M &m, const std::string typ, octave_idx_type k_arg, octave_idx_type p_arg, 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)
 
template OCTAVE_API octave_idx_type EigsComplexNonSymmetricMatrix< ComplexMatrix > (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 &cresid, std::ostream &os, double tol, bool rvec, bool cholB, int disp, int maxit)
 
template OCTAVE_API octave_idx_type EigsComplexNonSymmetricMatrix< SparseComplexMatrix > (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 &cresid, std::ostream &os, double tol, bool rvec, bool cholB, int disp, int maxit)
 
template<typename M >
octave_idx_type EigsComplexNonSymmetricMatrixShift (const M &m, Complex sigma, octave_idx_type k_arg, octave_idx_type p_arg, 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)
 
template OCTAVE_API octave_idx_type EigsComplexNonSymmetricMatrixShift< ComplexMatrix > (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 &cresid, std::ostream &os, double tol, bool rvec, bool cholB, int disp, int maxit)
 
template OCTAVE_API octave_idx_type EigsComplexNonSymmetricMatrixShift< SparseComplexMatrix > (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 &cresid, std::ostream &os, double tol, bool rvec, bool cholB, int disp, int maxit)
 
template<typename M >
octave_idx_type EigsRealNonSymmetricFunc (EigsFunc fcn, octave_idx_type n_arg, const std::string &_typ, double sigmar, octave_idx_type k_arg, octave_idx_type p_arg, 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)
 
template OCTAVE_API octave_idx_type EigsRealNonSymmetricFunc< Matrix > (EigsFunc fcn, 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, const Matrix &_b, ColumnVector &permB, ColumnVector &resid, std::ostream &os, double tol, bool rvec, bool cholB, int disp, int maxit)
 
template OCTAVE_API octave_idx_type EigsRealNonSymmetricFunc< SparseMatrix > (EigsFunc fcn, 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, const SparseMatrix &_b, ColumnVector &permB, ColumnVector &resid, std::ostream &os, double tol, bool rvec, bool cholB, int disp, int maxit)
 
template<typename M >
octave_idx_type EigsRealNonSymmetricMatrix (const M &m, const std::string typ, octave_idx_type k_arg, octave_idx_type p_arg, 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)
 
template OCTAVE_API octave_idx_type EigsRealNonSymmetricMatrix< Matrix > (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, bool rvec, bool cholB, int disp, int maxit)
 
template OCTAVE_API octave_idx_type EigsRealNonSymmetricMatrix< SparseMatrix > (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, bool rvec, bool cholB, int disp, int maxit)
 
template<typename M >
octave_idx_type EigsRealNonSymmetricMatrixShift (const M &m, double sigmar, octave_idx_type k_arg, octave_idx_type p_arg, 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)
 
template OCTAVE_API octave_idx_type EigsRealNonSymmetricMatrixShift< Matrix > (const Matrix &m, double sigmar, 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, bool rvec, bool cholB, int disp, int maxit)
 
template OCTAVE_API octave_idx_type EigsRealNonSymmetricMatrixShift< SparseMatrix > (const SparseMatrix &m, double sigmar, 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, bool rvec, bool cholB, int disp, int maxit)
 
template<typename M >
octave_idx_type EigsRealSymmetricFunc (EigsFunc fcn, octave_idx_type n_arg, const std::string &_typ, double sigma, octave_idx_type k_arg, octave_idx_type p_arg, 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)
 
template OCTAVE_API octave_idx_type EigsRealSymmetricFunc< Matrix > (EigsFunc fcn, 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, const Matrix &_b, ColumnVector &permB, ColumnVector &resid, std::ostream &os, double tol, bool rvec, bool cholB, int disp, int maxit)
 
template OCTAVE_API octave_idx_type EigsRealSymmetricFunc< SparseMatrix > (EigsFunc fcn, 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, const SparseMatrix &_b, ColumnVector &permB, ColumnVector &resid, std::ostream &os, double tol, bool rvec, bool cholB, int disp, int maxit)
 
template<typename M >
octave_idx_type EigsRealSymmetricMatrix (const M &m, const std::string typ, octave_idx_type k_arg, octave_idx_type p_arg, 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)
 
template OCTAVE_API octave_idx_type EigsRealSymmetricMatrix< Matrix > (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, bool rvec, bool cholB, int disp, int maxit)
 
template OCTAVE_API octave_idx_type EigsRealSymmetricMatrix< SparseMatrix > (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, bool rvec, bool cholB, int disp, int maxit)
 
template<typename M >
octave_idx_type EigsRealSymmetricMatrixShift (const M &m, double sigma, octave_idx_type k_arg, octave_idx_type p_arg, 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)
 
template OCTAVE_API octave_idx_type EigsRealSymmetricMatrixShift< Matrix > (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, bool rvec, bool cholB, int disp, int maxit)
 
template OCTAVE_API octave_idx_type EigsRealSymmetricMatrixShift< SparseMatrix > (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, bool rvec, bool cholB, int disp, int maxit)
 
template<typename SM , typename M >
static M ltsolve (const SM &L, const ColumnVector &Q, const M &m)
 
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, ColumnVector &r)
 
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, ColumnVector &r)
 
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, ColumnVector &r)
 
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, ColumnVector &r)
 
template<typename M , typename SM >
static octave_idx_type lusolve (const SM &L, const SM &U, M &m)
 
static bool make_cholb (ComplexMatrix &b, ComplexMatrix &bt, ColumnVector &permB)
 
static bool make_cholb (Matrix &b, Matrix &bt, ColumnVector &permB)
 
static bool make_cholb (SparseComplexMatrix &b, SparseComplexMatrix &bt, ColumnVector &permB)
 
static bool make_cholb (SparseMatrix &b, SparseMatrix &bt, ColumnVector &permB)
 
template<typename SM , typename M >
static M utsolve (const SM &U, const ColumnVector &Q, const M &m)
 
static bool vector_product (const ComplexMatrix &m, const Complex *x, Complex *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 SparseMatrix &m, const double *x, double *y)
 
static void warn_convergence (void)
 

Function Documentation

◆ arpack_errno2str()

std::string arpack_errno2str ( const octave_idx_type &  errnum,
const std::string &  fcn_name 
)

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

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

◆ EigsComplexNonSymmetricFunc()

template<typename M >
octave_idx_type EigsComplexNonSymmetricFunc ( EigsComplexFunc  fcn,
octave_idx_type  n_arg,
const std::string &  _typ,
Complex  sigma,
octave_idx_type  k_arg,
octave_idx_type  p_arg,
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 3553 of file eigs-base.cc.

References abs(), arpack_errno2str(), ColumnVector, ComplexColumnVector, current_liboctave_error_handler, d, F77_DBLE_CMPLX_ARG, F77_FUNC(), Array< T, Alloc >::fortran_vec(), Array< T, Alloc >::isempty(), isnan(), ltsolve(), M, make_cholb(), n, NaN, Array< T, Alloc >::numel(), OCTAVE_LOCAL_BUFFER, ComplexColumnVector::resize(), ComplexMatrix::resize(), utsolve(), vector_product(), x, and x_nint().

Referenced by F__eigs__().

◆ EigsComplexNonSymmetricFunc< ComplexMatrix >()

template OCTAVE_API octave_idx_type EigsComplexNonSymmetricFunc< ComplexMatrix > ( EigsComplexFunc  fcn,
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,
const ComplexMatrix &  _b,
ColumnVector &  permB,
ComplexColumnVector &  cresid,
std::ostream &  os,
double  tol,
bool  rvec,
bool  cholB,
int  disp,
int  maxit 
)

◆ EigsComplexNonSymmetricFunc< SparseComplexMatrix >()

template OCTAVE_API octave_idx_type EigsComplexNonSymmetricFunc< SparseComplexMatrix > ( EigsComplexFunc  fcn,
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,
const SparseComplexMatrix &  _b,
ColumnVector &  permB,
ComplexColumnVector &  cresid,
std::ostream &  os,
double  tol,
bool  rvec,
bool  cholB,
int  disp,
int  maxit 
)

◆ EigsComplexNonSymmetricMatrix()

template<typename M >
octave_idx_type EigsComplexNonSymmetricMatrix ( const M &  m,
const std::string  typ,
octave_idx_type  k_arg,
octave_idx_type  p_arg,
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 2927 of file eigs-base.cc.

References arpack_errno2str(), ColumnVector, ComplexColumnVector, current_liboctave_error_handler, d, F77_CONST_DBLE_CMPLX_ARG, F77_DBLE_CMPLX_ARG, F77_FUNC(), Array< T, Alloc >::fortran_vec(), Array< T, Alloc >::isempty(), isnan(), ltsolve(), M, m, make_cholb(), n, NaN, Array< T, Alloc >::numel(), OCTAVE_LOCAL_BUFFER, ComplexColumnVector::resize(), ComplexMatrix::resize(), utsolve(), vector_product(), and x_nint().

Referenced by F__eigs__().

◆ EigsComplexNonSymmetricMatrix< ComplexMatrix >()

template OCTAVE_API octave_idx_type EigsComplexNonSymmetricMatrix< ComplexMatrix > ( 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 &  cresid,
std::ostream &  os,
double  tol,
bool  rvec,
bool  cholB,
int  disp,
int  maxit 
)

◆ EigsComplexNonSymmetricMatrix< SparseComplexMatrix >()

template OCTAVE_API octave_idx_type EigsComplexNonSymmetricMatrix< SparseComplexMatrix > ( 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 &  cresid,
std::ostream &  os,
double  tol,
bool  rvec,
bool  cholB,
int  disp,
int  maxit 
)

◆ EigsComplexNonSymmetricMatrixShift()

template<typename M >
octave_idx_type EigsComplexNonSymmetricMatrixShift ( const M &  m,
Complex  sigma,
octave_idx_type  k_arg,
octave_idx_type  p_arg,
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 3230 of file eigs-base.cc.

References arpack_errno2str(), ComplexColumnVector, current_liboctave_error_handler, d, F77_CONST_DBLE_CMPLX_ARG, F77_DBLE_CMPLX_ARG, F77_FUNC(), Array< T, Alloc >::fortran_vec(), Array< T, Alloc >::isempty(), isnan(), LuAminusSigmaB(), lusolve(), M, m, n, NaN, Array< T, Alloc >::numel(), OCTAVE_LOCAL_BUFFER, Q, r, ComplexColumnVector::resize(), ComplexMatrix::resize(), vector_product(), and x_nint().

Referenced by F__eigs__().

◆ EigsComplexNonSymmetricMatrixShift< ComplexMatrix >()

template OCTAVE_API octave_idx_type EigsComplexNonSymmetricMatrixShift< ComplexMatrix > ( 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 &  cresid,
std::ostream &  os,
double  tol,
bool  rvec,
bool  cholB,
int  disp,
int  maxit 
)

◆ EigsComplexNonSymmetricMatrixShift< SparseComplexMatrix >()

template OCTAVE_API octave_idx_type EigsComplexNonSymmetricMatrixShift< SparseComplexMatrix > ( 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 &  cresid,
std::ostream &  os,
double  tol,
bool  rvec,
bool  cholB,
int  disp,
int  maxit 
)

◆ EigsRealNonSymmetricFunc()

template<typename M >
octave_idx_type EigsRealNonSymmetricFunc ( EigsFunc  fcn,
octave_idx_type  n_arg,
const std::string &  _typ,
double  sigmar,
octave_idx_type  k_arg,
octave_idx_type  p_arg,
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 2485 of file eigs-base.cc.

References abs(), arpack_errno2str(), ColumnVector, current_liboctave_error_handler, d, F77_FUNC(), Array< T, Alloc >::fortran_vec(), imag(), Array< T, Alloc >::isempty(), isnan(), ltsolve(), M, make_cholb(), n, NaN, Array< T, Alloc >::numel(), OCTAVE_LOCAL_BUFFER, ComplexColumnVector::resize(), ComplexMatrix::resize(), utsolve(), vector_product(), x, and x_nint().

Referenced by F__eigs__().

◆ EigsRealNonSymmetricFunc< Matrix >()

template OCTAVE_API octave_idx_type EigsRealNonSymmetricFunc< Matrix > ( EigsFunc  fcn,
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,
const Matrix &  _b,
ColumnVector &  permB,
ColumnVector &  resid,
std::ostream &  os,
double  tol,
bool  rvec,
bool  cholB,
int  disp,
int  maxit 
)

◆ EigsRealNonSymmetricFunc< SparseMatrix >()

template OCTAVE_API octave_idx_type EigsRealNonSymmetricFunc< SparseMatrix > ( EigsFunc  fcn,
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,
const SparseMatrix &  _b,
ColumnVector &  permB,
ColumnVector &  resid,
std::ostream &  os,
double  tol,
bool  rvec,
bool  cholB,
int  disp,
int  maxit 
)

◆ EigsRealNonSymmetricMatrix()

template<typename M >
octave_idx_type EigsRealNonSymmetricMatrix ( const M &  m,
const std::string  typ,
octave_idx_type  k_arg,
octave_idx_type  p_arg,
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 1762 of file eigs-base.cc.

References arpack_errno2str(), ColumnVector, current_liboctave_error_handler, d, F77_FUNC(), Array< T, Alloc >::fortran_vec(), imag(), Array< T, Alloc >::isempty(), isnan(), ltsolve(), M, m, make_cholb(), n, NaN, Array< T, Alloc >::numel(), OCTAVE_LOCAL_BUFFER, ComplexColumnVector::resize(), ComplexMatrix::resize(), utsolve(), vector_product(), and x_nint().

Referenced by F__eigs__().

◆ EigsRealNonSymmetricMatrix< Matrix >()

template OCTAVE_API octave_idx_type EigsRealNonSymmetricMatrix< Matrix > ( 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,
bool  rvec,
bool  cholB,
int  disp,
int  maxit 
)

◆ EigsRealNonSymmetricMatrix< SparseMatrix >()

template OCTAVE_API octave_idx_type EigsRealNonSymmetricMatrix< SparseMatrix > ( 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,
bool  rvec,
bool  cholB,
int  disp,
int  maxit 
)

◆ EigsRealNonSymmetricMatrixShift()

template<typename M >
octave_idx_type EigsRealNonSymmetricMatrixShift ( const M &  m,
double  sigmar,
octave_idx_type  k_arg,
octave_idx_type  p_arg,
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 2113 of file eigs-base.cc.

References arpack_errno2str(), ColumnVector, current_liboctave_error_handler, d, F77_FUNC(), Array< T, Alloc >::fortran_vec(), imag(), Array< T, Alloc >::isempty(), isnan(), LuAminusSigmaB(), lusolve(), M, m, n, NaN, Array< T, Alloc >::numel(), OCTAVE_LOCAL_BUFFER, Q, r, ComplexColumnVector::resize(), ComplexMatrix::resize(), vector_product(), and x_nint().

Referenced by F__eigs__().

◆ EigsRealNonSymmetricMatrixShift< Matrix >()

template OCTAVE_API octave_idx_type EigsRealNonSymmetricMatrixShift< Matrix > ( const Matrix &  m,
double  sigmar,
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,
bool  rvec,
bool  cholB,
int  disp,
int  maxit 
)

◆ EigsRealNonSymmetricMatrixShift< SparseMatrix >()

template OCTAVE_API octave_idx_type EigsRealNonSymmetricMatrixShift< SparseMatrix > ( const SparseMatrix &  m,
double  sigmar,
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,
bool  rvec,
bool  cholB,
int  disp,
int  maxit 
)

◆ EigsRealSymmetricFunc()

template<typename M >
octave_idx_type EigsRealSymmetricFunc ( EigsFunc  fcn,
octave_idx_type  n_arg,
const std::string &  _typ,
double  sigma,
octave_idx_type  k_arg,
octave_idx_type  p_arg,
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 1381 of file eigs-base.cc.

References abs(), arpack_errno2str(), ColumnVector, current_liboctave_error_handler, d, F77_FUNC(), Array< T, Alloc >::fortran_vec(), Array< T, Alloc >::isempty(), isnan(), ltsolve(), M, make_cholb(), n, NaN, Array< T, Alloc >::numel(), OCTAVE_LOCAL_BUFFER, ColumnVector::resize(), Matrix::resize(), utsolve(), vector_product(), x, and x_nint().

Referenced by F__eigs__().

◆ EigsRealSymmetricFunc< Matrix >()

template OCTAVE_API octave_idx_type EigsRealSymmetricFunc< Matrix > ( EigsFunc  fcn,
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,
const Matrix &  _b,
ColumnVector &  permB,
ColumnVector &  resid,
std::ostream &  os,
double  tol,
bool  rvec,
bool  cholB,
int  disp,
int  maxit 
)

◆ EigsRealSymmetricFunc< SparseMatrix >()

template OCTAVE_API octave_idx_type EigsRealSymmetricFunc< SparseMatrix > ( EigsFunc  fcn,
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,
const SparseMatrix &  _b,
ColumnVector &  permB,
ColumnVector &  resid,
std::ostream &  os,
double  tol,
bool  rvec,
bool  cholB,
int  disp,
int  maxit 
)

◆ EigsRealSymmetricMatrix()

template<typename M >
octave_idx_type EigsRealSymmetricMatrix ( const M &  m,
const std::string  typ,
octave_idx_type  k_arg,
octave_idx_type  p_arg,
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 786 of file eigs-base.cc.

References arpack_errno2str(), ColumnVector, current_liboctave_error_handler, d, F77_FUNC(), Array< T, Alloc >::fortran_vec(), Array< T, Alloc >::isempty(), isnan(), ltsolve(), M, m, make_cholb(), n, NaN, Array< T, Alloc >::numel(), OCTAVE_LOCAL_BUFFER, ColumnVector::resize(), Matrix::resize(), utsolve(), vector_product(), and x_nint().

Referenced by F__eigs__().

◆ EigsRealSymmetricMatrix< Matrix >()

template OCTAVE_API octave_idx_type EigsRealSymmetricMatrix< Matrix > ( 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,
bool  rvec,
bool  cholB,
int  disp,
int  maxit 
)

◆ EigsRealSymmetricMatrix< SparseMatrix >()

template OCTAVE_API octave_idx_type EigsRealSymmetricMatrix< SparseMatrix > ( 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,
bool  rvec,
bool  cholB,
int  disp,
int  maxit 
)

◆ EigsRealSymmetricMatrixShift()

template<typename M >
octave_idx_type EigsRealSymmetricMatrixShift ( const M &  m,
double  sigma,
octave_idx_type  k_arg,
octave_idx_type  p_arg,
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 1078 of file eigs-base.cc.

References arpack_errno2str(), ColumnVector, current_liboctave_error_handler, d, F77_FUNC(), Array< T, Alloc >::fortran_vec(), Array< T, Alloc >::isempty(), isnan(), LuAminusSigmaB(), lusolve(), M, m, n, NaN, Array< T, Alloc >::numel(), OCTAVE_LOCAL_BUFFER, Q, r, ColumnVector::resize(), Matrix::resize(), vector_product(), and x_nint().

Referenced by F__eigs__().

◆ EigsRealSymmetricMatrixShift< Matrix >()

template OCTAVE_API octave_idx_type EigsRealSymmetricMatrixShift< Matrix > ( 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,
bool  rvec,
bool  cholB,
int  disp,
int  maxit 
)

◆ EigsRealSymmetricMatrixShift< SparseMatrix >()

template OCTAVE_API octave_idx_type EigsRealSymmetricMatrixShift< SparseMatrix > ( 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,
bool  rvec,
bool  cholB,
int  disp,
int  maxit 
)

◆ ltsolve()

template<typename SM , typename M >
static M ltsolve ( const SM &  L,
const ColumnVector &  Q,
const M &  m 
)
static

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

References MatrixType::Lower, M, m, n, and Q.

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

◆ LuAminusSigmaB() [1/4]

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,
ColumnVector &  r 
)
static

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

References abs(), Array< T, Alloc >::cols(), d, Array< T, Alloc >::data(), Array< T, Alloc >::fortran_vec(), ComplexMatrix::hermitian(), imag(), Array< T, Alloc >::isempty(), isnan(), m, n, NaN, Array< T, Alloc >::numel(), Q, r, real(), Array< T, Alloc >::rows(), warn_convergence(), and Array< T, Alloc >::xelem().

◆ LuAminusSigmaB() [2/4]

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,
ColumnVector &  r 
)
static

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

References abs(), Array< T, Alloc >::cols(), d, Array< T, Alloc >::data(), Array< T, Alloc >::fortran_vec(), Array< T, Alloc >::isempty(), isnan(), m, n, NaN, Array< T, Alloc >::numel(), Q, r, Array< T, Alloc >::rows(), Matrix::transpose(), warn_convergence(), and Array< T, Alloc >::xelem().

◆ LuAminusSigmaB() [3/4]

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,
ColumnVector &  r 
)
static

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

References abs(), d, SparseComplexMatrix::hermitian(), imag(), Sparse< T, Alloc >::isempty(), isnan(), m, Matrix, n, NaN, Array< T, Alloc >::numel(), Q, r, real(), sparse_lu< SparseComplexMatrix >, SparseMatrix::transpose(), warn_convergence(), Sparse< T, Alloc >::xcidx(), Sparse< T, Alloc >::xdata(), and Sparse< T, Alloc >::xridx().

◆ LuAminusSigmaB() [4/4]

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,
ColumnVector &  r 
)
static

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

References abs(), d, Sparse< T, Alloc >::isempty(), isnan(), m, n, NaN, Array< T, Alloc >::numel(), Q, r, sparse_lu< SparseMatrix >, SparseMatrix::transpose(), warn_convergence(), Sparse< T, Alloc >::xcidx(), Sparse< T, Alloc >::xdata(), and Sparse< T, Alloc >::xridx().

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

◆ lusolve()

template<typename M , typename SM >
static octave_idx_type lusolve ( const SM &  L,
const SM &  U,
M &  m 
)
static

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

References MatrixType::Lower, m, and MatrixType::Upper.

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

◆ make_cholb() [1/4]

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

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

References Array< T, Alloc >::cols(), ColumnVector, ComplexMatrix::hermitian(), and n.

◆ make_cholb() [2/4]

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

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

References Array< T, Alloc >::cols(), ColumnVector, n, and Matrix::transpose().

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

◆ make_cholb() [3/4]

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

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

References SparseComplexMatrix::hermitian().

◆ make_cholb() [4/4]

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

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

References sparse_chol< SparseMatrix >, and SparseMatrix::transpose().

◆ utsolve()

template<typename SM , typename M >
static M utsolve ( const SM &  U,
const ColumnVector &  Q,
const M &  m 
)
static

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

References M, m, n, Q, and MatrixType::Upper.

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

◆ vector_product() [1/4]

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

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

References F77_CONST_DBLE_CMPLX_ARG, F77_DBLE_CMPLX_ARG, F77_XFCN, m, and x.

◆ vector_product() [2/4]

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

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

References F77_XFCN, m, and x.

◆ vector_product() [3/4]

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

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

References m, and x.

◆ vector_product() [4/4]

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

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

References m, and x.

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

◆ warn_convergence()

static void warn_convergence ( void  )
static

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

Referenced by LuAminusSigmaB().