GNU Octave 7.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 fun, 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 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, 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 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, 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 fun, 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 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, 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 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, 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 fun, 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 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, 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 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, 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  fun,
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, octave::rand::distribution(), F77_DBLE_CMPLX_ARG, F77_FUNC(), Array< T, Alloc >::fortran_vec(), Array< T, Alloc >::isempty(), octave::math::isnan(), ltsolve(), M, make_cholb(), NaN, Array< T, Alloc >::numel(), OCTAVE_LOCAL_BUFFER, ComplexColumnVector::resize(), ComplexMatrix::resize(), utsolve(), octave::rand::vector(), vector_product(), x, and octave::math::x_nint().

Referenced by F__eigs__().

◆ EigsComplexNonSymmetricFunc< ComplexMatrix >()

template OCTAVE_API octave_idx_type EigsComplexNonSymmetricFunc< ComplexMatrix > ( 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,
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  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,
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, octave::rand::distribution(), F77_CONST_DBLE_CMPLX_ARG, F77_DBLE_CMPLX_ARG, F77_FUNC(), Array< T, Alloc >::fortran_vec(), Array< T, Alloc >::isempty(), octave::math::isnan(), ltsolve(), M, make_cholb(), NaN, Array< T, Alloc >::numel(), OCTAVE_LOCAL_BUFFER, ComplexColumnVector::resize(), ComplexMatrix::resize(), utsolve(), octave::rand::vector(), vector_product(), and octave::math::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, octave::rand::distribution(), F77_CONST_DBLE_CMPLX_ARG, F77_DBLE_CMPLX_ARG, F77_FUNC(), Array< T, Alloc >::fortran_vec(), Array< T, Alloc >::isempty(), octave::math::isnan(), LuAminusSigmaB(), lusolve(), M, NaN, Array< T, Alloc >::numel(), OCTAVE_LOCAL_BUFFER, Q, ComplexColumnVector::resize(), ComplexMatrix::resize(), octave::rand::vector(), vector_product(), and octave::math::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  fun,
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, octave::rand::distribution(), F77_FUNC(), Array< T, Alloc >::fortran_vec(), imag(), Array< T, Alloc >::isempty(), octave::math::isnan(), ltsolve(), M, make_cholb(), NaN, Array< T, Alloc >::numel(), OCTAVE_LOCAL_BUFFER, ComplexColumnVector::resize(), ComplexMatrix::resize(), utsolve(), octave::rand::vector(), vector_product(), x, and octave::math::x_nint().

Referenced by F__eigs__().

◆ EigsRealNonSymmetricFunc< Matrix >()

template OCTAVE_API octave_idx_type EigsRealNonSymmetricFunc< Matrix > ( 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,
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  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,
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, octave::rand::distribution(), F77_FUNC(), Array< T, Alloc >::fortran_vec(), imag(), Array< T, Alloc >::isempty(), octave::math::isnan(), ltsolve(), M, make_cholb(), NaN, Array< T, Alloc >::numel(), OCTAVE_LOCAL_BUFFER, ComplexColumnVector::resize(), ComplexMatrix::resize(), utsolve(), octave::rand::vector(), vector_product(), and octave::math::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, octave::rand::distribution(), F77_FUNC(), Array< T, Alloc >::fortran_vec(), imag(), Array< T, Alloc >::isempty(), octave::math::isnan(), LuAminusSigmaB(), lusolve(), M, NaN, Array< T, Alloc >::numel(), OCTAVE_LOCAL_BUFFER, Q, ComplexColumnVector::resize(), ComplexMatrix::resize(), octave::rand::vector(), vector_product(), and octave::math::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  fun,
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, octave::rand::distribution(), F77_FUNC(), Array< T, Alloc >::fortran_vec(), Array< T, Alloc >::isempty(), octave::math::isnan(), ltsolve(), M, make_cholb(), NaN, Array< T, Alloc >::numel(), OCTAVE_LOCAL_BUFFER, ColumnVector::resize(), Matrix::resize(), utsolve(), octave::rand::vector(), vector_product(), x, and octave::math::x_nint().

Referenced by F__eigs__().

◆ EigsRealSymmetricFunc< Matrix >()

template OCTAVE_API octave_idx_type EigsRealSymmetricFunc< Matrix > ( 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,
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  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,
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, octave::rand::distribution(), F77_FUNC(), Array< T, Alloc >::fortran_vec(), Array< T, Alloc >::isempty(), octave::math::isnan(), ltsolve(), M, make_cholb(), NaN, Array< T, Alloc >::numel(), OCTAVE_LOCAL_BUFFER, ColumnVector::resize(), Matrix::resize(), utsolve(), octave::rand::vector(), vector_product(), and octave::math::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, octave::rand::distribution(), F77_FUNC(), Array< T, Alloc >::fortran_vec(), Array< T, Alloc >::isempty(), octave::math::isnan(), LuAminusSigmaB(), lusolve(), M, NaN, Array< T, Alloc >::numel(), OCTAVE_LOCAL_BUFFER, Q, ColumnVector::resize(), Matrix::resize(), octave::rand::vector(), vector_product(), and octave::math::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, 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(), octave::math::isnan(), octave::math::lu< T >::L(), NaN, Array< T, Alloc >::numel(), octave::math::lu< T >::P_vec(), Q, real(), Array< T, Alloc >::rows(), octave::math::lu< T >::U(), 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(), octave::math::isnan(), octave::math::lu< T >::L(), NaN, Array< T, Alloc >::numel(), octave::math::lu< T >::P_vec(), Q, Array< T, Alloc >::rows(), Matrix::transpose(), octave::math::lu< T >::U(), 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(), octave::math::sparse_lu< lu_type >::col_perm(), d, SparseComplexMatrix::hermitian(), imag(), Sparse< T, Alloc >::isempty(), octave::math::isnan(), octave::math::sparse_lu< lu_type >::L(), Matrix, NaN, Array< T, Alloc >::numel(), Q, octave::math::sparse_lu< lu_type >::R(), real(), octave::math::sparse_lu< lu_type >::row_perm(), Sparse< T, Alloc >::rows(), SparseMatrix::transpose(), octave::math::sparse_lu< lu_type >::U(), 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(), octave::math::sparse_lu< lu_type >::col_perm(), d, Sparse< T, Alloc >::isempty(), octave::math::isnan(), octave::math::sparse_lu< lu_type >::L(), NaN, Array< T, Alloc >::numel(), Q, octave::math::sparse_lu< lu_type >::R(), octave::math::sparse_lu< lu_type >::row_perm(), Sparse< T, Alloc >::rows(), SparseMatrix::transpose(), octave::math::sparse_lu< lu_type >::U(), 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, 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 octave::math::chol< T >::chol_matrix(), Array< T, Alloc >::cols(), ColumnVector, and ComplexMatrix::hermitian().

◆ 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 octave::math::chol< T >::chol_matrix(), Array< T, Alloc >::cols(), ColumnVector, 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(), octave::math::sparse_chol< chol_type >::L(), and octave::math::sparse_chol< chol_type >::perm().

◆ 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 octave::math::sparse_chol< chol_type >::L(), octave::math::sparse_chol< chol_type >::perm(), 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, 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 Array< T, Alloc >::cols(), Array< T, Alloc >::data(), F77_CONST_DBLE_CMPLX_ARG, F77_DBLE_CMPLX_ARG, F77_XFCN, Array< T, Alloc >::rows(), 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 Array< T, Alloc >::cols(), Array< T, Alloc >::data(), F77_XFCN, Array< T, Alloc >::rows(), 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 Sparse< T, Alloc >::cidx(), Sparse< T, Alloc >::cols(), Sparse< T, Alloc >::data(), Sparse< T, Alloc >::ridx(), 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 Sparse< T, Alloc >::cidx(), Sparse< T, Alloc >::cols(), Sparse< T, Alloc >::data(), Sparse< T, Alloc >::ridx(), 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().