ergo
mat Namespace Reference

Namespaces

namespace  arn
 

Classes

class  Acceptable
 
class  AcceptableMaxIter
 
struct  AccessMap
 
class  Allocator
 
class  AllocatorManager
 
struct  ATAMatrix
 
struct  CongrTransErrorMatrix
 
class  DebugLevelHigh
 
class  DebugLevelLow
 
class  DebugLevelMedium
 
struct  DiffMatrix
 
class  EuclTruncationBase
 
class  EuclTruncationCongrTransMeasure
 Truncation of general matrices with impact on matrix triple multiply as error measure. More...
 
class  EuclTruncationGeneral
 Truncation of general matrices. More...
 
class  EuclTruncationSymm
 Truncation of symmetric matrices. More...
 
class  EuclTruncationSymmElementLevel
 Truncation of symmetric matrices at the element level (used for mixed norm truncation) More...
 
class  EuclTruncationSymmWithZ
 Truncation of symmetric matrices with Z. More...
 
class  Failure
 
class  FileWritable
 Write and read objects to/from file. More...
 
struct  Gblas
 
class  Interval
 
class  Matrix
 Matrix class and heart of the matrix library. More...
 
class  Matrix< Treal >
 
class  MatrixBase
 Base class for matrix API. More...
 
class  MatrixGeneral
 Normal matrix. More...
 
class  MatrixHierarchicBase
 Base class for Matrix and Matrix specialization. More...
 
class  MatrixSymmetric
 Symmetric matrix. More...
 
class  MatrixTriangular
 Upper non-unit triangular matrix. More...
 
class  MemUsage
 
class  Params
 
class  SingletonForTimings
 
class  SizesAndBlocks
 Describes dimensions of matrix and its blocks on all levels. More...
 
class  TC2
 Trace correcting purification. More...
 
class  Time
 
struct  TripleMatrix
 
class  ValidPtr
 Smart pointer class to control access to object. More...
 
class  Vector
 Vector class. More...
 
class  Vector< Treal >
 
class  VectorGeneral
 
class  VectorHierarchicBase
 Base class for Vector and Vector specialization. More...
 
struct  XmY
 This proxy expresses the result of substraction of two objects, of possibly different types, TX and TY. More...
 
struct  XpY
 This proxy expresses the result of addition of two objects, of possibly different types, TX and TY. More...
 
struct  Xtrans
 This proxy expresses the result of transposition of an object of type TX. More...
 
struct  XY
 This proxy expresses the result of multiplication of two objects, of possibly different types, TX and TY. More...
 
struct  XYZ
 This proxy expresses the result of multiplication of three objects, of possibly different types, TX, TY, and TZ. More...
 
struct  XYZpUV
 This proxy expresses the result of multiplication of three objects added to two other multiplied objects. More...
 

Enumerations

enum  property { zero , ful }
 
enum  normType { frobNorm , euclNorm , mixedNorm }
 
enum  side { left , right }
 
enum  inchversion { unstable , stable }
 
enum  matrix_type { matrix_matr , matrix_symm , matrix_triang }
 

Functions

template<>
floatallocateElements< float > (int n)
 
template<>
void freeElements (float *ptr)
 
template<>
doubleallocateElements< double > (int n)
 
template<>
void freeElements (double *ptr)
 
template<class T >
TallocateElements (int n)
 
template<class T >
void freeElements (T *ptr)
 
template<>
floatallocateElements< float > (int n)
 
template<>
void freeElements (float *ptr)
 
template<>
doubleallocateElements< double > (int n)
 
template<>
void freeElements (double *ptr)
 
template<typename Treal >
int sign (Treal value)
 Sign function returns the sign of the input.
 
template<typename Treal , typename Tfun >
Treal bisection (Tfun const &fun, Treal min, Treal max, Treal const tol)
 Bisection algorithm for root finding.
 
static long int get_file_size (const char *fileName)
 
static void copy_file (const char *sourceFileName, const char *destFileName)
 
template<class Treal >
static Treal maxdiff (const Treal *f1, const Treal *f2, int size)
 
template<class Treal >
static Treal maxdiff_tri (const Treal *f1, const Treal *f2, int size)
 
template<class Treal >
static Treal frobdiff (const Treal *f1, const Treal *f2, int size)
 
template<typename Treal , typename Trealonfile >
static void fileread (Treal *ptr, int size, FILE *file)
 
template<typename Treal , typename Tmatrix >
static void read_matrix (Tmatrix &A, char const *const matrixPath, int const size)
 
template<typename Treal , typename Trealonfile , typename Tmatrix >
static void read_sparse_matrix (Tmatrix &A, char const *const rowPath, char const *const colPath, char const *const valPath, int const nval)
 
template<typename Treal >
static void read_xyz (Treal *x, Treal *y, Treal *z, char *atomsPath, int const natoms, int const size)
 
template<typename Treal >
Interval< TrealsqrtInt (Interval< Treal > const &other)
 
template<typename Treal >
std::ostream & operator<< (std::ostream &s, Interval< Treal > const &in)
 
template<typename Tmatrix , typename Treal >
Interval< TrealeuclIfSmall (Tmatrix const &M, Treal const requestedAbsAccuracy, Treal const requestedRelAccuracy, Treal const maxAbsVal, typename Tmatrix::VectorType *const eVecPtr=0, int maxIter=-1)
 Returns interval containing the Euclidean norm of the matrix M.
 
template<class T >
static void gemm (const char *ta, const char *tb, const int *n, const int *k, const int *l, const T *alpha, const T *A, const int *lda, const T *B, const int *ldb, const T *beta, T *C, const int *ldc)
 
template<class T >
static void pptrf (const char *uplo, const int *n, T *ap, int *info)
 
template<class T >
static void spgst (const int *itype, const char *uplo, const int *n, T *ap, const T *bp, int *info)
 
template<class T >
static void tptri (const char *uplo, const char *diag, const int *n, T *ap, int *info)
 
template<class T >
static void trmm (const char *side, const char *uplo, const char *transa, const char *diag, const int *m, const int *n, const T *alpha, const T *A, const int *lda, T *B, const int *ldb)
 
template<class T >
static void sygv (const int *itype, const char *jobz, const char *uplo, const int *n, T *A, const int *lda, T *B, const int *ldb, T *w, T *work, const int *lwork, int *info)
 
template<class T >
static void ggev (const char *jobbl, const char *jobvr, const int *n, T *A, const int *lda, T *B, const int *ldb, T *alphar, T *alphai, T *beta, T *vl, const int *ldvl, T *vr, const int *ldvr, T *work, const int *lwork, int *info)
 
template<class T >
static void potrf (const char *uplo, const int *n, T *A, const int *lda, int *info)
 
template<class T >
static void trtri (const char *uplo, const char *diag, const int *n, T *A, const int *lda, int *info)
 
template<class T >
static void syrk (const char *uplo, const char *trans, const int *n, const int *k, const T *alpha, const T *A, const int *lda, const T *beta, T *C, const int *ldc)
 
template<class T >
static void symm (const char *side, const char *uplo, const int *m, const int *n, const T *alpha, const T *A, const int *lda, const T *B, const int *ldb, const T *beta, T *C, const int *ldc)
 
template<class T >
static void pocon (const char *uplo, const int *n, const T *A, const int *lda, const T *anorm, T *rcond, T *work, int *iwork, int *info)
 
template<class T >
static void stevx (const char *jobz, const char *range, const int *n, T *d, T *e, const T *vl, const T *vu, const int *il, const int *iu, const T *abstol, int *m, T *w, T *z, const int *ldz, T *work, int *iwork, int *ifail, int *info)
 
template<class T >
static void stevr (const char *jobz, const char *range, const int *n, T *d, T *e, const T *vl, const T *vu, const int *il, const int *iu, const T *abstol, int *m, T *w, T *z, const int *ldz, int *isuppz, T *work, int *lwork, int *iwork, int *liwork, int *info)
 
template<class T >
static void syev (const char *jobz, const char *uplo, const int *n, T *a, const int *lda, T *w, T *work, const int *lwork, int *info)
 
template<class T >
static void gemv (const char *ta, const int *m, const int *n, const T *alpha, const T *A, const int *lda, const T *x, const int *incx, const T *beta, T *y, const int *incy)
 
template<class T >
static void symv (const char *uplo, const int *n, const T *alpha, const T *A, const int *lda, const T *x, const int *incx, const T *beta, T *y, const int *incy)
 
template<class T >
static void trmv (const char *uplo, const char *trans, const char *diag, const int *n, const T *A, const int *lda, T *x, const int *incx)
 
template<class T >
static void scal (const int *n, const T *da, T *dx, const int *incx)
 
template<class T >
static T dot (const int *n, const T *dx, const int *incx, const T *dy, const int *incy)
 
template<class T >
static void axpy (const int *n, const T *da, const T *dx, const int *incx, T *dy, const int *incy)
 
template<>
void gemm< double > (const char *ta, const char *tb, const int *n, const int *k, const int *l, const double *alpha, const double *A, const int *lda, const double *B, const int *ldb, const double *beta, double *C, const int *ldc)
 
template<>
void pptrf< double > (const char *uplo, const int *n, double *ap, int *info)
 
template<>
void spgst< double > (const int *itype, const char *uplo, const int *n, double *ap, const double *bp, int *info)
 
template<>
void tptri< double > (const char *uplo, const char *diag, const int *n, double *ap, int *info)
 
template<>
void trmm< double > (const char *side, const char *uplo, const char *transa, const char *diag, const int *m, const int *n, const double *alpha, const double *A, const int *lda, double *B, const int *ldb)
 
template<>
void sygv< double > (const int *itype, const char *jobz, const char *uplo, const int *n, double *A, const int *lda, double *B, const int *ldb, double *w, double *work, const int *lwork, int *info)
 
template<>
void ggev< double > (const char *jobbl, const char *jobvr, const int *n, double *A, const int *lda, double *B, const int *ldb, double *alphar, double *alphai, double *beta, double *vl, const int *ldvl, double *vr, const int *ldvr, double *work, const int *lwork, int *info)
 
template<>
void potrf< double > (const char *uplo, const int *n, double *A, const int *lda, int *info)
 
template<>
void trtri< double > (const char *uplo, const char *diag, const int *n, double *A, const int *lda, int *info)
 
template<>
void syrk< double > (const char *uplo, const char *trans, const int *n, const int *k, const double *alpha, const double *A, const int *lda, const double *beta, double *C, const int *ldc)
 
template<>
void symm< double > (const char *side, const char *uplo, const int *m, const int *n, const double *alpha, const double *A, const int *lda, const double *B, const int *ldb, const double *beta, double *C, const int *ldc)
 
template<>
void pocon< double > (const char *uplo, const int *n, const double *A, const int *lda, const double *anorm, double *rcond, double *work, int *iwork, int *info)
 
template<>
void stevx< double > (const char *jobz, const char *range, const int *n, double *d, double *e, const double *vl, const double *vu, const int *il, const int *iu, const double *abstol, int *m, double *w, double *z, const int *ldz, double *work, int *iwork, int *ifail, int *info)
 
template<>
void stevr< double > (const char *jobz, const char *range, const int *n, double *d, double *e, const double *vl, const double *vu, const int *il, const int *iu, const double *abstol, int *m, double *w, double *z, const int *ldz, int *isuppz, double *work, int *lwork, int *iwork, int *liwork, int *info)
 
template<>
void syev< double > (const char *jobz, const char *uplo, const int *n, double *a, const int *lda, double *w, double *work, const int *lwork, int *info)
 
template<>
void gemv< double > (const char *ta, const int *m, const int *n, const double *alpha, const double *A, const int *lda, const double *x, const int *incx, const double *beta, double *y, const int *incy)
 
template<>
void symv< double > (const char *uplo, const int *n, const double *alpha, const double *A, const int *lda, const double *x, const int *incx, const double *beta, double *y, const int *incy)
 
template<>
void trmv< double > (const char *uplo, const char *trans, const char *diag, const int *n, const double *A, const int *lda, double *x, const int *incx)
 
template<>
void scal< double > (const int *n, const double *da, double *dx, const int *incx)
 
template<>
double dot< double > (const int *n, const double *dx, const int *incx, const double *dy, const int *incy)
 
template<>
void axpy< double > (const int *n, const double *da, const double *dx, const int *incx, double *dy, const int *incy)
 
template<>
void gemm< float > (const char *ta, const char *tb, const int *n, const int *k, const int *l, const float *alpha, const float *A, const int *lda, const float *B, const int *ldb, const float *beta, float *C, const int *ldc)
 
template<>
void pptrf< float > (const char *uplo, const int *n, float *ap, int *info)
 
template<>
void spgst< float > (const int *itype, const char *uplo, const int *n, float *ap, const float *bp, int *info)
 
template<>
void tptri< float > (const char *uplo, const char *diag, const int *n, float *ap, int *info)
 
template<>
void trmm< float > (const char *side, const char *uplo, const char *transa, const char *diag, const int *m, const int *n, const float *alpha, const float *A, const int *lda, float *B, const int *ldb)
 
template<>
void sygv< float > (const int *itype, const char *jobz, const char *uplo, const int *n, float *A, const int *lda, float *B, const int *ldb, float *w, float *work, const int *lwork, int *info)
 
template<>
void ggev< float > (const char *jobbl, const char *jobvr, const int *n, float *A, const int *lda, float *B, const int *ldb, float *alphar, float *alphai, float *beta, float *vl, const int *ldvl, float *vr, const int *ldvr, float *work, const int *lwork, int *info)
 
template<>
void potrf< float > (const char *uplo, const int *n, float *A, const int *lda, int *info)
 
template<>
void trtri< float > (const char *uplo, const char *diag, const int *n, float *A, const int *lda, int *info)
 
template<>
void syrk< float > (const char *uplo, const char *trans, const int *n, const int *k, const float *alpha, const float *A, const int *lda, const float *beta, float *C, const int *ldc)
 
template<>
void symm< float > (const char *side, const char *uplo, const int *m, const int *n, const float *alpha, const float *A, const int *lda, const float *B, const int *ldb, const float *beta, float *C, const int *ldc)
 
template<>
void pocon< float > (const char *uplo, const int *n, const float *A, const int *lda, const float *anorm, float *rcond, float *work, int *iwork, int *info)
 
template<>
void stevx< float > (const char *jobz, const char *range, const int *n, float *d, float *e, const float *vl, const float *vu, const int *il, const int *iu, const float *abstol, int *m, float *w, float *z, const int *ldz, float *work, int *iwork, int *ifail, int *info)
 
template<>
void stevr< float > (const char *jobz, const char *range, const int *n, float *d, float *e, const float *vl, const float *vu, const int *il, const int *iu, const float *abstol, int *m, float *w, float *z, const int *ldz, int *isuppz, float *work, int *lwork, int *iwork, int *liwork, int *info)
 
template<>
void syev< float > (const char *jobz, const char *uplo, const int *n, float *a, const int *lda, float *w, float *work, const int *lwork, int *info)
 
template<>
void gemv< float > (const char *ta, const int *m, const int *n, const float *alpha, const float *A, const int *lda, const float *x, const int *incx, const float *beta, float *y, const int *incy)
 
template<>
void symv< float > (const char *uplo, const int *n, const float *alpha, const float *A, const int *lda, const float *x, const int *incx, const float *beta, float *y, const int *incy)
 
template<>
void trmv< float > (const char *uplo, const char *trans, const char *diag, const int *n, const float *A, const int *lda, float *x, const int *incx)
 
template<>
void scal< float > (const int *n, const float *da, float *dx, const int *incx)
 
template<>
void axpy< float > (const int *n, const float *da, const float *dx, const int *incx, float *dy, const int *incy)
 
template<class Treal >
static void fulltopacked (const Treal *full, Treal *packed, const int size)
 
template<class Treal >
static void packedtofull (const Treal *packed, Treal *full, const int size)
 
template<class Treal >
static void tripackedtofull (const Treal *packed, Treal *full, const int size)
 
template<class Treal >
static void trifulltofull (Treal *trifull, const int size)
 
template<class RandomIt >
static void random_shuffle (RandomIt first, RandomIt last)
 
normType getNormType (const char *normStr)
 
std::string getNormTypeString (normType nType)
 
template<typename Treal >
static Treal getMachineEpsilon ()
 
template<typename Treal >
static Treal getRelPrecision ()
 
template<typename TX >
Xtrans< TXtranspose (TX const &A)
 Transposition.
 
template<typename TX >
Xtrans< TXtranspose (const Xtrans< TX > &xtrans)
 Transposition.
 
template<typename TX , typename TY >
XY< TX, TYoperator* (Xtrans< TX > const &trAA, Xtrans< TY > const &trBB)
 Multiplication of two transposition proxys holding objects of type TX and TY respectively.
 
template<typename TX , typename TY >
XY< TX, TYoperator* (TX const &AA, Xtrans< TY > const &trBB)
 Multiplication of an object of type TX with a tranposition proxy
holding an object of type TY.
 
template<typename TX , typename TY >
XY< TX, TYoperator* (Xtrans< TX > const &trAA, TY const &BB)
 Multiplication of a tranposition proxy holding an object of type TX with an object of type TY.
 
template<typename TX , typename TY >
XY< TX, TYoperator* (TX const &AA, TY const &BB)
 Multiplication of an object of type TX with an object of type TY.
 
template<typename TX , typename TY , typename TZ >
XYZ< TX, TY, TZoperator* (XY< TX, TY > const &AB, Xtrans< TZ > const &trCC)
 Multiplication of a multiplication proxy XY with a transposition proxy Xtrans.
 
template<typename TX , typename TY , typename TZ >
XYZ< TX, TY, TZoperator* (XY< TX, TY > const &AB, TZ const &CC)
 Multiplication of a multiplication proxy XY with an object of type TZ.
 
template<typename TX , typename TY , typename TZ , typename TU , typename TV >
XYZpUV< TX, TY, TZ, TU, TVoperator+ (XYZ< TX, TY, TZ > const &ABC, XY< TU, TV > const &DE)
 Addition of two multiplication proxys XYZ and XY.
 
template<typename TX , typename TY >
XpY< TX, TYoperator+ (TX const &AA, TY const &BB)
 Addition of two objects of type TX and TY.
 
template<typename TX , typename TY >
XmY< TX, TYoperator- (TX const &AA, TY const &BB)
 Substraction of two objects of type TX and TY.
 
template<typename Treal , typename Tmatrix >
Treal trace (const XYZ< Treal, MatrixGeneral< Treal, Tmatrix >, MatrixGeneral< Treal, Tmatrix > > &smm)
 
template<typename Treal , typename Tmatrix , typename Top >
Treal accumulate (MatrixSymmetric< Treal, Tmatrix > &A, Top &op)
 Performs operation specified in 'op' on all nonzero matrix elements and sums up the result and returns it.
 
template<class Treal >
void quicksort (const Treal *value, int *index, int low, int high)
 
template<typename Treal , typename Tvector >
Treal operator* (Xtrans< VectorGeneral< Treal, Tvector > > const &xT, VectorGeneral< Treal, Tvector > const &y)
 transpose(x) * y Scalar (dot) product of two vectors
 

Detailed Description

design principles

Heart of matrix library: MatrixHierarchicBase is the base class both for higher levels in the hierarchic Matrix (first part of this file) and for the specialization at the lowest level (second part of this file). Hence, functionality common to these class templates can be implemented in MatrixHierarchicBase.

The interface: The interface is in MatrixBase.h MatrixSymmetric.h MatrixGeneral.h MatrixTriangular.h.

Functions added to this file will not be visible. They have to be called from the interface classes.

Example usage: Check out API_test.cc

Enumeration Type Documentation

◆ inchversion

Enumerator
unstable 
stable 

◆ matrix_type

Enumerator
matrix_matr 
matrix_symm 
matrix_triang 

◆ normType

Enumerator
frobNorm 
euclNorm 
mixedNorm 

◆ property

Enumerator
zero 
ful 

◆ side

Enumerator
left 
right 

Function Documentation

◆ accumulate()

Treal mat::accumulate ( MatrixSymmetric< Treal, Tmatrix > &  A,
Top op 
)

Performs operation specified in 'op' on all nonzero matrix elements and sums up the result and returns it.

References A, and getMachineEpsilon().

Referenced by testAccumulation().

◆ allocateElements()

template<class T >
T * mat::allocateElements ( int  n)

References getMachineEpsilon().

◆ allocateElements< double >() [1/2]

◆ allocateElements< double >() [2/2]

◆ allocateElements< float >() [1/2]

◆ allocateElements< float >() [2/2]

◆ axpy()

template<class T >
static void mat::axpy ( const int n,
const T da,
const T dx,
const int incx,
T dy,
const int incy 
)
inlinestatic

◆ axpy< double >()

template<>
void mat::axpy< double > ( const int n,
const double da,
const double dx,
const int incx,
double dy,
const int incy 
)
inline

◆ axpy< float >()

template<>
void mat::axpy< float > ( const int n,
const float da,
const float dx,
const int incx,
float dy,
const int incy 
)
inline

◆ bisection()

template<typename Treal , typename Tfun >
Treal mat::bisection ( Tfun const &  fun,
Treal  min,
Treal  max,
Treal const  tol 
)

Bisection algorithm for root finding.

The bisection method finds the root of a function in the interval [min, max], or more precisely the place where the function changes sign. It is assumed that the function only changes sign once in the given interval.

The function is given by a class that has a member function named eval that evaluates the function in the given point.

References getMachineEpsilon(), max, min(), sign(), and template_blas_fabs().

Referenced by mat::TC2< Treal, Tmatrix >::fermi_level(), mat::TC2< Treal, Tmatrix >::homo(), and mat::TC2< Treal, Tmatrix >::lumo().

◆ copy_file()

static void mat::copy_file ( const char sourceFileName,
const char destFileName 
)
static

◆ dot()

template<class T >
static T mat::dot ( const int n,
const T dx,
const int incx,
const T dy,
const int incy 
)
inlinestatic

◆ dot< double >()

template<>
double mat::dot< double > ( const int n,
const double dx,
const int incx,
const double dy,
const int incy 
)
inline

◆ euclIfSmall()

Interval< Treal > mat::euclIfSmall ( Tmatrix const &  M,
Treal const  requestedAbsAccuracy,
Treal const  requestedRelAccuracy,
Treal const  maxAbsVal,
typename Tmatrix::VectorType *const  eVecPtr = 0,
int  maxIter = -1 
)

◆ fileread()

static void mat::fileread ( Treal ptr,
int  size,
FILE file 
)
static

References getMachineEpsilon().

◆ freeElements() [1/5]

template<>
void mat::freeElements ( double ptr)

◆ freeElements() [2/5]

template<>
void mat::freeElements ( double ptr)

◆ freeElements() [3/5]

◆ freeElements() [4/5]

◆ freeElements() [5/5]

template<class T >
void mat::freeElements ( T ptr)

◆ frobdiff()

template<class Treal >
static Treal mat::frobdiff ( const Treal f1,
const Treal f2,
int  size 
)
static

◆ fulltopacked()

template<class Treal >
static void mat::fulltopacked ( const Treal full,
Treal packed,
const int  size 
)
static

References getMachineEpsilon().

◆ gemm()

template<class T >
static void mat::gemm ( const char ta,
const char tb,
const int n,
const int k,
const int l,
const T alpha,
const T A,
const int lda,
const T B,
const int ldb,
const T beta,
T C,
const int ldc 
)
inlinestatic

◆ gemm< double >()

template<>
void mat::gemm< double > ( const char ta,
const char tb,
const int n,
const int k,
const int l,
const double alpha,
const double A,
const int lda,
const double B,
const int ldb,
const double beta,
double C,
const int ldc 
)
inline

◆ gemm< float >()

template<>
void mat::gemm< float > ( const char ta,
const char tb,
const int n,
const int k,
const int l,
const float alpha,
const float A,
const int lda,
const float B,
const int ldb,
const float beta,
float C,
const int ldc 
)
inline

◆ gemv()

template<class T >
static void mat::gemv ( const char ta,
const int m,
const int n,
const T alpha,
const T A,
const int lda,
const T x,
const int incx,
const T beta,
T y,
const int incy 
)
inlinestatic

◆ gemv< double >()

template<>
void mat::gemv< double > ( const char ta,
const int m,
const int n,
const double alpha,
const double A,
const int lda,
const double x,
const int incx,
const double beta,
double y,
const int incy 
)
inline

◆ gemv< float >()

template<>
void mat::gemv< float > ( const char ta,
const int m,
const int n,
const float alpha,
const float A,
const int lda,
const float x,
const int incx,
const float beta,
float y,
const int incy 
)
inline

◆ get_file_size()

static long int mat::get_file_size ( const char fileName)
static

◆ getMachineEpsilon()

template<typename Treal >
static Treal mat::getMachineEpsilon ( )
inlinestatic

References getMachineEpsilon().

Referenced by accumulate(), mat::Matrix< Treal >::add(), mat::Matrix< Treal, Telement >::add(), mat::Matrix< Treal, Telement >::add_abs_col_sums(), mat::Matrix< Treal >::add_abs_col_sums(), mat::MatrixSymmetric< Treal, Tmatrix >::add_values(), mat::Vector< Treal, Telement >::addFromFull(), mat::Matrix< Treal, Telement >::addIdentity(), mat::Matrix< Treal >::addIdentity(), mat::SingletonForTimings::addTime(), mat::Matrix< Treal, Telement >::addValues(), mat::Matrix< Treal >::addValues(), mat::Matrix< Treal, Telement >::addValues(), mat::Matrix< Treal >::addValues(), mat::Allocator< Treal >::alloc(), mat::AllocatorManager< Treal >::alloc(), mat::Matrix< Treal, Telement >::allocate(), mat::Matrix< Treal >::allocate(), mat::Vector< Treal, Telement >::allocate(), mat::Vector< Treal >::allocate(), allocateElements(), mat::Allocator< Treal >::Allocator(), mat::AllocatorManager< Treal >::AllocatorManager(), mat::DebugLevelHigh::assertAlways(), mat::DebugLevelHigh::assertDebug(), mat::VectorGeneral< Treal, Tvector >::assign_from_full(), assign_from_full_matrix(), mat::MatrixGeneral< Treal, Tmatrix >::assign_from_sparse(), mat::MatrixSymmetric< Treal, Tmatrix >::assign_from_sparse(), mat::MatrixTriangular< Treal, Tmatrix >::assign_from_sparse(), mat::MatrixGeneral< Treal, Tmatrix >::assign_from_sparse(), mat::MatrixSymmetric< Treal, Tmatrix >::assign_from_sparse(), mat::MatrixGeneral< Treal, Tmatrix >::assign_from_sparse(), mat::MatrixSymmetric< Treal, Tmatrix >::assign_from_sparse(), mat::Matrix< Treal >::assignDiffFrobNormsLowestLevel(), mat::Matrix< Treal, Telement >::assignDiffFrobNormsLowestLevel(), mat::Matrix< Treal >::assignFrobNormsLowestLevel(), mat::Matrix< Treal, Telement >::assignFrobNormsLowestLevel(), mat::Matrix< Treal, Telement >::assignFromFull(), mat::Matrix< Treal >::assignFromFull(), mat::MatrixGeneral< Treal, Tmatrix >::assignFromFull(), mat::MatrixSymmetric< Treal, Tmatrix >::assignFromFull(), mat::MatrixSymmetric< Treal, Tmatrix >::assignFromFull(), mat::Matrix< Treal, Telement >::assignFromSparse(), mat::Matrix< Treal >::assignFromSparse(), mat::Matrix< Treal, Telement >::assignFromSparse(), mat::Matrix< Treal >::assignFromSparse(), axpy(), mat::Vector< Treal >::axpy(), mat::Vector< Treal, Telement >::axpy(), axpy< double >(), axpy< float >(), bisection(), BoysFuncAccurate(), PurificationGeneral< MatrixType >::compute_spectrum_bounds(), mat::arn::LanczosLargestMagnitudeEig< Treal, Tmatrix, Tvector >::computeEigenPairTri(), mat::arn::LanczosSeveralLargestEig< Treal, Tmatrix, Tvector >::computeEigenPairTri(), mat::arn::LanczosLargestMagnitudeEig< Treal, Tmatrix, Tvector >::converged(), mat::arn::LanczosLargestMagnitudeEigIfSmall< Treal, Tmatrix, Tvector >::converged(), mat::arn::LanczosSeveralLargestEig< Treal, Tmatrix, Tvector >::converged(), mat::arn::LanczosSeveralLargestEig< Treal, Tmatrix, Tvector >::converged_ith(), copy_file(), mat::FileWritable::copyFromFile(), mat::FileWritable::copyToFile(), mat::arn::Lanczos< Treal, Tmatrix, Tvector >::copyTridiag(), mat::arn::LanczosSeveralLargestEig< Treal, Tmatrix, Tvector >::copyTridiag(), mat::Interval< Treal >::cover(), mat::MatrixSymmetric< Treal, Tmatrix >::diff(), mat::MatrixSymmetric< Treal, Tmatrix >::diffIfSmall(), dot(), mat::Vector< Treal, Telement >::dot(), dot< double >(), dotIsBroken(), per::Perturbation< Treal, Tmatrix, Tvector >::dryRun(), mat::MatrixGeneral< Treal, Tmatrix >::eucl(), mat::MatrixSymmetric< Treal, Tmatrix >::eucl(), mat::MatrixTriangular< Treal, Tmatrix >::eucl(), mat::MatrixSymmetric< Treal, Tmatrix >::eucl_diff(), mat::MatrixGeneral< Treal, Tmatrix >::eucl_thresh(), mat::MatrixTriangular< Treal, Tmatrix >::eucl_thresh(), mat::MatrixSymmetric< Treal, Tmatrix >::eucl_thresh(), mat::MatrixTriangular< Treal, Tmatrix >::eucl_thresh_congr_trans_measure(), mat::MatrixSymmetric< Treal, Tmatrix >::euclDiffIfSmall(), euclIfSmall(), mat::EuclTruncationSymm< Tmatrix, Treal >::euclIfSmall(), mat::EuclTruncationSymmWithZ< Tmatrix, TmatrixZ, Treal >::euclIfSmall(), mat::EuclTruncationGeneral< Tmatrix, Treal >::euclIfSmall(), mat::EuclTruncationCongrTransMeasure< Tmatrix, TmatrixB, Treal >::euclIfSmall(), mat::TC2< Treal, Tmatrix >::Fun::eval(), fileread(), mat::FileWritable::fileSize(), mat::FileWritable::FileWritable(), mat::FileWritable::FileWritable(), mat::Allocator< Treal >::free(), mat::AllocatorManager< Treal >::free(), mat::Matrix< Treal >::frob_squared_thresh(), mat::Matrix< Treal, Telement >::frob_squared_thresh(), frobdiff(), frobdiff(), mat::Matrix< Treal, Telement >::frobSquared(), mat::Matrix< Treal >::frobSquared(), mat::Matrix< Treal >::frobSquaredDiff(), mat::Matrix< Treal, Telement >::frobSquaredDiff(), mat::Matrix< Treal >::frobThreshElementLevel(), mat::Matrix< Treal, Telement >::frobThreshElementLevel(), mat::Matrix< Treal, Telement >::frobThreshLowestLevel(), mat::Matrix< Treal, Telement >::fullMatrix(), mat::Matrix< Treal >::fullMatrix(), mat::MatrixGeneral< Treal, Tmatrix >::fullMatrix(), mat::MatrixSymmetric< Treal, Tmatrix >::fullMatrix(), fulltopacked(), mat::Vector< Treal, Telement >::fullVector(), mat::Vector< Treal >::fullVector(), mat::TC2< Treal, Tmatrix >::Fun::Fun(), mat::Matrix< Treal, Telement >::geAccumulateWith(), mat::Matrix< Treal >::geAccumulateWith(), mat::Matrix< Treal >::gemm(), mat::Matrix< Treal, Telement >::gemm(), gemm(), gemm< double >(), gemm< float >(), mat::Matrix< Treal, Telement >::gemm_upper_tr_only(), mat::Vector< Treal >::gemv(), mat::Vector< Treal, Telement >::gemv(), gemv(), gemv< double >(), gemv< float >(), mat::Matrix< Treal, Telement >::gershgorin(), mat::Matrix< Treal >::gershgorin(), mat::MatrixGeneral< Treal, Tmatrix >::get_all_values(), mat::MatrixSymmetric< Treal, Tmatrix >::get_all_values(), get_dens_from_cmo_FermiDiracDistr(), mat::Matrix< Treal, Telement >::get_diagonal(), mat::Matrix< Treal >::get_diagonal(), PurificationGeneral< MatrixType >::get_epsilon(), get_eucl_norm_try_different_acc(), get_file_size(), mat::arn::LanczosSeveralLargestEig< Treal, Tmatrix, Tvector >::get_ith_eigenpair(), mat::MatrixGeneral< Treal, Tmatrix >::get_values(), mat::MatrixSymmetric< Treal, Tmatrix >::get_values(), get_wall_seconds(), mat::Time::get_wall_seconds(), mat::Matrix< Treal, Telement >::getAllValues(), mat::SizesAndBlocks::getBlockSizeVector(), mat::DiffMatrix< Tmatrix, Treal >::getCols(), mat::ATAMatrix< Tmatrix, Treal >::getCols(), mat::TripleMatrix< Tmatrix, Tmatrix2, Treal >::getCols(), mat::CongrTransErrorMatrix< Tmatrix, Tmatrix2, Treal >::getCols(), mat::MatrixBase< Treal, Tmatrix >::getCols(), mat::MatrixHierarchicBase< Treal, Telement >::getCols(), mat::arn::MatrixTridiagSymmetric< Treal >::getEigsByIndex(), mat::arn::MatrixTridiagSymmetric< Treal >::getEigsByInterval(), mat::arn::Lanczos< Treal, Tmatrix, Tvector >::getEigVector(), mat::arn::LanczosSeveralLargestEig< Treal, Tmatrix, Tvector >::getEigVector(), SCF_general::GetEuclideanNormOfMatrix(), mat::Matrix< Treal, Telement >::getFrobSqElementLevel(), mat::Matrix< Treal >::getFrobSqElementLevel(), mat::Matrix< Treal, Telement >::getFrobSqLowestLevel(), mat::Matrix< Treal >::getFrobSqLowestLevel(), mat::EuclTruncationSymm< Tmatrix, Treal >::getFrobTruncBounds(), mat::EuclTruncationSymmWithZ< Tmatrix, TmatrixZ, Treal >::getFrobTruncBounds(), mat::EuclTruncationGeneral< Tmatrix, Treal >::getFrobTruncBounds(), mat::EuclTruncationCongrTransMeasure< Tmatrix, TmatrixB, Treal >::getFrobTruncBounds(), mat::arn::LanczosLargestMagnitudeEig< Treal, Tmatrix, Tvector >::getLargestMagnitudeEig(), mat::arn::LanczosLargestMagnitudeEig< Treal, Tmatrix, Tvector >::getLargestMagnitudeEigPair(), getMachineEpsilon(), mat::Params::getMatrixParallelLevel(), mat::MemUsage::getMemUsage(), getNormType(), getNormTypeString(), mat::Params::getNProcs(), mat::MemUsage::getNumberFromBuffer(), mat::MatrixSymmetric< Treal, Tmatrix >::getPermutedAndSymmetrized(), mat::MatrixBase< Treal, Tmatrix >::getPermutedIndexes(), getRelPrecision(), mat::MatrixBase< Treal, Tmatrix >::getRows(), mat::MatrixHierarchicBase< Treal, Telement >::getRows(), mat::MatrixSymmetric< Treal, Tmatrix >::getSizesAndBlocksForFrobNormMat(), mat::SizesAndBlocks::getSizesAndBlocksForLowerLevel(), mat::AllocatorManager< Treal >::getStatistics(), mat::FileWritable::getStatsCount(), mat::FileWritable::getStatsFileSizes(), mat::FileWritable::getStatsTime(), mat::Matrix< Treal, Telement >::getValues(), mat::Matrix< Treal, Telement >::getValues(), mat::Matrix< Treal >::getValues(), mat::Matrix< Treal >::getValues(), ggev(), ggev< double >(), ggev< float >(), mat::ValidPtr< Tobj >::haveDataStructureSet(), mat::MatrixTriangular< Treal, Tmatrix >::inch(), mat::MatrixTriangular< Treal, Tmatrix >::inch(), mat::arn::Lanczos< Treal, Tmatrix, Tvector >::increaseCapacity(), mat::arn::LanczosSeveralLargestEig< Treal, Tmatrix, Tvector >::increaseCapacity(), mat::arn::MatrixTridiagSymmetric< Treal >::increaseCapacity(), mat::AllocatorManager< Treal >::init(), mat::MatrixBase< Treal, Tmatrix >::inMemorySet(), mat::VectorGeneral< Treal, Tvector >::inMemorySet(), mat::ValidPtr< Tobj >::inMemorySet(), mat::AllocatorManager< Treal >::instance(), mat::FileWritable::Stats::instance(), mat::SingletonForTimings::instance(), mat::FileWritable::Manager::instance_prot(), mat::Interval< Treal >::intersect(), mat::Interval< Treal >::intersect_always_non_empty(), mat::Interval< Treal >::invPuriStep(), mat::Interval< Treal >::invPuriStep(), mat::arn::Lanczos< Treal, Tmatrix, Tvector >::Lanczos(), mat::arn::LanczosSeveralLargestEig< Treal, Tmatrix, Tvector >::LanczosSeveralLargestEig(), main(), main(), mainFun(), mainFun(), mat::MatrixBase< Treal, Tmatrix >::MatrixBase(), mat::MatrixHierarchicBase< Treal, Telement >::MatrixHierarchicBase(), mat::MatrixSymmetric< Treal, Tmatrix >::MatrixSymmetric(), mat::DiffMatrix< Tmatrix, Treal >::matVecProd(), mat::MatrixSymmetric< Treal, Tmatrix >::matVecProd(), mat::TripleMatrix< Tmatrix, Tmatrix2, Treal >::matVecProd(), mat::CongrTransErrorMatrix< Tmatrix, Tmatrix2, Treal >::matVecProd(), mat::Matrix< Treal, Telement >::maxAbsValue(), mat::Matrix< Treal >::maxAbsValue(), maxdiff(), maxdiff(), maxdiff_tri(), maxdiff_tri(), mat::Matrix< Treal, Telement >::memory_usage(), mat::Matrix< Treal >::memory_usage(), mat::MatrixSymmetric< Treal, Tmatrix >::mixed_diff(), mat::MatrixSymmetric< Treal, Tmatrix >::mixed_norm(), mat::MatrixSymmetric< Treal, Tmatrix >::mixed_norm_thresh(), mat::Matrix< Treal, Telement >::nnz(), mat::Matrix< Treal, Telement >::nosymToSym(), mat::VectorHierarchicBase< Treal, Telement >::operator(), mat::VectorHierarchicBase< Treal, Telement >::operator(), mat::Interval< Treal >::operator*(), operator*(), operator*(), operator*(), operator*(), operator*(), operator*(), operator*(), mat::Matrix< Treal, Telement >::operator*=(), mat::Matrix< Treal >::operator*=(), mat::Vector< Treal, Telement >::operator*=(), mat::Interval< Treal >::operator+(), mat::Interval< Treal >::operator+(), operator+(), operator+(), mat::VectorGeneral< Treal, Tvector >::operator+=(), mat::MatrixGeneral< Treal, Tmatrix >::operator+=(), mat::MatrixSymmetric< Treal, Tmatrix >::operator+=(), mat::MatrixGeneral< Treal, Tmatrix >::operator+=(), mat::VectorGeneral< Treal, Tvector >::operator+=(), mat::MatrixGeneral< Treal, Tmatrix >::operator+=(), mat::MatrixSymmetric< Treal, Tmatrix >::operator+=(), mat::MatrixGeneral< Treal, Tmatrix >::operator+=(), mat::VectorGeneral< Treal, Tvector >::operator+=(), mat::MatrixGeneral< Treal, Tmatrix >::operator+=(), mat::MatrixSymmetric< Treal, Tmatrix >::operator+=(), mat::MatrixGeneral< Treal, Tmatrix >::operator+=(), mat::MatrixSymmetric< Treal, Tmatrix >::operator+=(), mat::MatrixTriangular< Treal, Tmatrix >::operator+=(), mat::Interval< Treal >::operator-(), mat::Interval< Treal >::operator-(), operator-(), mat::MatrixGeneral< Treal, Tmatrix >::operator-=(), mat::MatrixSymmetric< Treal, Tmatrix >::operator-=(), mat::MatrixSymmetric< Treal, Tmatrix >::operator-=(), mat::Interval< Treal >::operator/(), operator<<(), mat::MatrixBase< Treal, Tmatrix >::operator=(), mat::MatrixSymmetric< Treal, Tmatrix >::operator=(), mat::MatrixHierarchicBase< Treal, Telement >::operator=(), mat::MatrixTriangular< Treal, Tmatrix >::operator=(), mat::MatrixGeneral< Treal, Tmatrix >::operator=(), mat::VectorGeneral< Treal, Tvector >::operator=(), mat::VectorHierarchicBase< Treal, Telement >::operator=(), mat::MatrixBase< Treal, Tmatrix >::operator=(), mat::MatrixGeneral< Treal, Tmatrix >::operator=(), mat::MatrixGeneral< Treal, Tmatrix >::operator=(), mat::MatrixGeneral< Treal, Tmatrix >::operator=(), mat::VectorGeneral< Treal, Tvector >::operator=(), mat::MatrixGeneral< Treal, Tmatrix >::operator=(), mat::MatrixGeneral< Treal, Tmatrix >::operator=(), mat::VectorGeneral< Treal, Tvector >::operator=(), mat::MatrixSymmetric< Treal, Tmatrix >::operator=(), mat::MatrixGeneral< Treal, Tmatrix >::operator=(), mat::MatrixSymmetric< Treal, Tmatrix >::operator=(), mat::MatrixGeneral< Treal, Tmatrix >::operator=(), mat::MatrixGeneral< Treal, Tmatrix >::operator=(), mat::VectorGeneral< Treal, Tvector >::operator=(), mat::MatrixGeneral< Treal, Tmatrix >::operator=(), mat::MatrixSymmetric< Treal, Tmatrix >::operator=(), mat::MatrixGeneral< Treal, Tmatrix >::operator=(), mat::VectorGeneral< Treal, Tvector >::operator=(), mat::MatrixGeneral< Treal, Tmatrix >::operator=(), mat::MatrixGeneral< Treal, Tmatrix >::operator=(), mat::MatrixSymmetric< Treal, Tmatrix >::operator=(), mat::MatrixGeneral< Treal, Tmatrix >::operator=(), mat::VectorGeneral< Treal, Tvector >::operator=(), mat::MatrixGeneral< Treal, Tmatrix >::operator=(), mat::MatrixGeneral< Treal, Tmatrix >::operator=(), mat::MatrixSymmetric< Treal, Tmatrix >::operator=(), mat::VectorGeneral< Treal, Tvector >::operator=(), mat::FileWritable::operator=(), mat::Matrix< Treal, Telement >::operator=(), mat::Matrix< Treal >::operator=(), mat::MatrixGeneral< Treal, Tmatrix >::operator=(), mat::MatrixSymmetric< Treal, Tmatrix >::operator=(), mat::MatrixTriangular< Treal, Tmatrix >::operator=(), mat::Vector< Treal, Telement >::operator=(), mat::Vector< Treal >::operator=(), mat::VectorGeneral< Treal, Tvector >::operator=(), mat::SizesAndBlocks::operator=(), mat::MatrixSymmetric< Treal, Tmatrix >::operator=(), mat::MatrixGeneral< Treal, Tmatrix >::operator=(), mat::MatrixSymmetric< Treal, Tmatrix >::operator=(), mat::MatrixSymmetric< Treal, Tmatrix >::operator=(), mat::MatrixBase< Treal, Tmatrix >::operator==(), mat::MatrixHierarchicBase< Treal, Telement >::operator==(), mat::SizesAndBlocks::operator==(), mat::Interval< Treal >::overlap(), mat::Allocator< Treal >::ownsPtr(), packedtofull(), pocon(), pocon< double >(), pocon< float >(), potrf(), potrf< double >(), potrf< float >(), pptrf(), pptrf< double >(), pptrf< float >(), mat::TC2< Treal, Tmatrix >::purify(), mat::Interval< Treal >::puriStep(), mat::Interval< Treal >::puriStep(), mat::DiffMatrix< Tmatrix, Treal >::quickEuclBounds(), mat::ATAMatrix< Tmatrix, Treal >::quickEuclBounds(), mat::TripleMatrix< Tmatrix, Tmatrix2, Treal >::quickEuclBounds(), mat::CongrTransErrorMatrix< Tmatrix, Tmatrix2, Treal >::quickEuclBounds(), mat::MatrixSymmetric< Treal, Tmatrix >::quickEuclBounds(), quicksort(), mat::Matrix< Treal, Telement >::random(), mat::Matrix< Treal >::random(), mat::Vector< Treal, Telement >::random(), mat::Vector< Treal >::random(), random_shuffle(), mat::Matrix< Treal, Telement >::randomZeroStructure(), mat::Matrix< Treal >::randomZeroStructure(), mat::MatrixGeneral< Treal, Tmatrix >::read_from_buffer(), mat::MatrixSymmetric< Treal, Tmatrix >::read_from_buffer(), mat::MatrixTriangular< Treal, Tmatrix >::read_from_buffer(), read_matrix(), read_sparse_matrix(), read_xyz(), mat::FileWritable::readFromFile(), mat::Matrix< Treal, Telement >::readFromFile(), mat::Matrix< Treal >::readFromFile(), mat::Vector< Treal, Telement >::readFromFile(), mat::Vector< Treal >::readFromFile(), mat::MatrixBase< Treal, Tmatrix >::readFromFileBase(), mat::VectorGeneral< Treal, Tvector >::readFromFileProt(), mat::FileWritable::Manager::registerObj(), mat::arn::LanczosLargestMagnitudeEig< Treal, Tmatrix, Tvector >::rerun(), mat::MatrixHierarchicBase< Treal, Telement >::resetCols(), mat::MatrixHierarchicBase< Treal, Telement >::resetRows(), mat::VectorHierarchicBase< Treal, Telement >::resetRows(), mat::VectorGeneral< Treal, Tvector >::resetSizesAndBlocks(), mat::MatrixBase< Treal, Tmatrix >::resetSizesAndBlocks(), mat::arn::Lanczos< Treal, Tmatrix, Tvector >::restart(), mat::arn::LanczosSeveralLargestEig< Treal, Tmatrix, Tvector >::run(), mat::EuclTruncationBase< Tmatrix, Treal >::run(), scal(), scal< double >(), scal< float >(), SCF_general::SCF_general(), mat::arn::LanczosSeveralLargestEig< Treal, Tmatrix, Tvector >::selective_orth(), expRule< Treal >::set(), mat::arn::LanczosLargestMagnitudeEig< Treal, Tmatrix, Tvector >::setAbsTol(), mat::arn::LanczosSeveralLargestEig< Treal, Tmatrix, Tvector >::setAbsTol(), mat::Matrix< Treal, Telement >::setElementsByRule(), mat::Matrix< Treal >::setElementsByRule(), mat::Params::setMatrixParallelLevel(), mat::Params::setNProcs(), mat::FileWritable::setPath(), mat::arn::LanczosLargestMagnitudeEig< Treal, Tmatrix, Tvector >::setRelTol(), mat::arn::LanczosSeveralLargestEig< Treal, Tmatrix, Tvector >::setRelTol(), mat::SizesAndBlocks::setup(), mat::SizesAndBlocks::SizesAndBlocks(), mat::SizesAndBlocks::SizesAndBlocks(), small_calculation_core(), spgst(), spgst< double >(), spgst< float >(), sqrtInt(), mat::Matrix< Treal >::ssmm(), mat::Matrix< Treal, Telement >::ssmm(), mat::Matrix< Treal, Telement >::ssmm_upper_tr_only(), mat::arn::Lanczos< Treal, Tmatrix, Tvector >::step(), mat::arn::LanczosSeveralLargestEig< Treal, Tmatrix, Tvector >::step(), stevr(), stevr< double >(), stevr< float >(), stevx(), stevx< double >(), stevx< float >(), mat::MatrixHierarchicBase< Treal, Telement >::swap(), mat::ValidPtr< Tobj >::swap(), mat::Matrix< Treal, Telement >::sy_gershgorin(), mat::Matrix< Treal >::sy_gershgorin(), mat::Matrix< Treal, Telement >::sy_nnz(), mat::Matrix< Treal, Telement >::sy_nvalues(), mat::Matrix< Treal >::sy_trace_ab(), mat::Matrix< Treal, Telement >::sy_trace_ab(), mat::Matrix< Treal, Telement >::syAccumulateWith(), mat::Matrix< Treal >::syAccumulateWith(), mat::Matrix< Treal, Telement >::syAddValues(), mat::Matrix< Treal >::syAddValues(), mat::Matrix< Treal >::syAssignDiffFrobNormsLowestLevel(), mat::Matrix< Treal, Telement >::syAssignDiffFrobNormsLowestLevel(), mat::Matrix< Treal >::syAssignFrobNormsLowestLevel(), mat::Matrix< Treal, Telement >::syAssignFrobNormsLowestLevel(), mat::Matrix< Treal, Telement >::syAssignFromSparse(), mat::Matrix< Treal >::syAssignFromSparse(), syev(), syev< double >(), syev< float >(), mat::Matrix< Treal, Telement >::syFrobSquared(), mat::Matrix< Treal >::syFrobSquared(), mat::Matrix< Treal >::syFrobSquaredDiff(), mat::Matrix< Treal, Telement >::syFrobSquaredDiff(), mat::Matrix< Treal, Telement >::syFullMatrix(), mat::Matrix< Treal >::syFullMatrix(), mat::Matrix< Treal, Telement >::syGetValues(), mat::Matrix< Treal >::syGetValues(), sygv(), sygv< double >(), sygv< float >(), mat::Matrix< Treal >::syInch(), mat::Matrix< Treal, Telement >::syInch(), symm(), mat::Matrix< Treal >::symm(), mat::Matrix< Treal, Telement >::symm(), symm< double >(), symm< float >(), mat::Matrix< Treal, Telement >::symToNosym(), mat::Vector< Treal >::symv(), mat::Vector< Treal, Telement >::symv(), symv(), symv< double >(), symv< float >(), mat::Matrix< Treal, Telement >::syRandom(), mat::Matrix< Treal >::syRandom(), mat::Matrix< Treal, Telement >::syRandomZeroStructure(), mat::Matrix< Treal >::syRandomZeroStructure(), syrk(), mat::Matrix< Treal >::syrk(), mat::Matrix< Treal, Telement >::syrk(), syrk< double >(), syrk< float >(), mat::Matrix< Treal, Telement >::sySetElementsByRule(), mat::Matrix< Treal >::sySetElementsByRule(), mat::Matrix< Treal >::sysq(), mat::Matrix< Treal, Telement >::sysq(), mat::Matrix< Treal, Telement >::sytr_upper_tr_only(), mat::Matrix< Treal, Telement >::syUpTriFullMatrix(), mat::Matrix< Treal >::syUpTriFullMatrix(), mat::TC2< Treal, Tmatrix >::TC2(), test_small(), test_small(), testAccumulation(), mat::Time::toc(), tptri(), tptri< double >(), tptri< float >(), mat::Matrix< Treal, Telement >::trace(), mat::Matrix< Treal >::trace(), trace(), mat::Matrix< Treal >::trace_ab(), mat::Matrix< Treal, Telement >::trace_ab(), mat::Matrix< Treal >::trace_aTb(), mat::Matrix< Treal, Telement >::trace_aTb(), transpose(), mat::Matrix< Treal >::transpose(), mat::Matrix< Treal, Telement >::transpose(), transpose(), trifulltofull(), tripackedtofull(), trmm(), mat::Matrix< Treal >::trmm(), mat::Matrix< Treal, Telement >::trmm(), trmm< double >(), trmm< float >(), mat::Matrix< Treal >::trmm_upper_tr_only(), mat::Matrix< Treal, Telement >::trmm_upper_tr_only(), mat::Vector< Treal >::trmv(), mat::Vector< Treal, Telement >::trmv(), trmv(), trmv< double >(), trmv< float >(), mat::Matrix< Treal, Telement >::trSetElementsByRule(), trtri(), trtri< double >(), trtri< float >(), mat::Matrix< Treal >::truncateAccordingToSparsityPattern(), mat::Matrix< Treal, Telement >::truncateAccordingToSparsityPattern(), try_integral_diffs(), mat::FileWritable::Manager::unRegisterObj(), mat::VectorGeneral< Treal, Tvector >::VectorGeneral(), mat::VectorGeneral< Treal, Tvector >::VectorGeneral(), mat::VectorHierarchicBase< Treal, Telement >::VectorHierarchicBase(), verify_idempotency(), mat::SizesAndBlocks::whichBlock(), mat::MatrixGeneral< Treal, Tmatrix >::write_to_buffer(), mat::MatrixSymmetric< Treal, Tmatrix >::write_to_buffer(), mat::MatrixTriangular< Treal, Tmatrix >::write_to_buffer(), mat::MatrixBase< Treal, Tmatrix >::write_to_buffer_base(), mat::MatrixBase< Treal, Tmatrix >::write_to_buffer_count(), mat::FileWritable::writeAndReadAll(), mat::FileWritable::writeToFile(), mat::Matrix< Treal, Telement >::writeToFile(), mat::Matrix< Treal >::writeToFile(), mat::Vector< Treal, Telement >::writeToFile(), mat::Vector< Treal >::writeToFile(), mat::MatrixBase< Treal, Tmatrix >::writeToFileBase(), mat::VectorGeneral< Treal, Tvector >::writeToFileProt(), x_times_ln_x(), mat::AllocatorManager< Treal >::~AllocatorManager(), mat::FileWritable::~FileWritable(), and mat::arn::LanczosSeveralLargestEig< Treal, Tmatrix, Tvector >::~LanczosSeveralLargestEig().

◆ getNormType()

normType mat::getNormType ( const char normStr)

References euclNorm, frobNorm, getMachineEpsilon(), and mixedNorm.

Referenced by es_run().

◆ getNormTypeString()

std::string mat::getNormTypeString ( normType  nType)

◆ getRelPrecision()

template<typename Treal >
static Treal mat::getRelPrecision ( )
inlinestatic

References getMachineEpsilon().

◆ ggev()

template<class T >
static void mat::ggev ( const char jobbl,
const char jobvr,
const int n,
T A,
const int lda,
T B,
const int ldb,
T alphar,
T alphai,
T beta,
T vl,
const int ldvl,
T vr,
const int ldvr,
T work,
const int lwork,
int info 
)
inlinestatic

◆ ggev< double >()

template<>
void mat::ggev< double > ( const char jobbl,
const char jobvr,
const int n,
double A,
const int lda,
double B,
const int ldb,
double alphar,
double alphai,
double beta,
double vl,
const int ldvl,
double vr,
const int ldvr,
double work,
const int lwork,
int info 
)
inline

◆ ggev< float >()

template<>
void mat::ggev< float > ( const char jobbl,
const char jobvr,
const int n,
float A,
const int lda,
float B,
const int ldb,
float alphar,
float alphai,
float beta,
float vl,
const int ldvl,
float vr,
const int ldvr,
float work,
const int lwork,
int info 
)
inline

◆ maxdiff()

template<class Treal >
static Treal mat::maxdiff ( const Treal f1,
const Treal f2,
int  size 
)
static

◆ maxdiff_tri()

template<class Treal >
static Treal mat::maxdiff_tri ( const Treal f1,
const Treal f2,
int  size 
)
static

◆ operator*() [1/7]

template<typename TX , typename TY >
XY< TX, TY > mat::operator* ( TX const &  AA,
TY const &  BB 
)
inline

Multiplication of an object of type TX with an object of type TY.

Returns multiplication proxy XY.

See also
XY
operator*(Xtrans<TX> const &, Xtrans<TY> const &)
operator*(TX const &, Xtrans<TY> const &)
operator*(Xtrans<TX> const &, TY const &)

References getMachineEpsilon().

◆ operator*() [2/7]

template<typename TX , typename TY >
XY< TX, TY > mat::operator* ( TX const &  AA,
Xtrans< TY > const &  trBB 
)
inline

Multiplication of an object of type TX with a tranposition proxy
holding an object of type TY.

Returns multiplication proxy XY.

See also
XY
Xtrans
operator*(Xtrans<TX> const &, Xtrans<TY> const &)
operator*(Xtrans<TX> const &, TY const &)
operator*(TX const &, TY const &)

References getMachineEpsilon().

◆ operator*() [3/7]

template<typename TX , typename TY >
XY< TX, TY > mat::operator* ( Xtrans< TX > const &  trAA,
TY const &  BB 
)
inline

Multiplication of a tranposition proxy holding an object of type TX with an object of type TY.

Returns multiplication proxy XY.

See also
XY
Xtrans
operator*(Xtrans<TX> const &, Xtrans<TY> const &)
operator*(TX const &, Xtrans<TY> const &)
operator*(TX const &, TY const &)

References getMachineEpsilon().

◆ operator*() [4/7]

template<typename TX , typename TY >
XY< TX, TY > mat::operator* ( Xtrans< TX > const &  trAA,
Xtrans< TY > const &  trBB 
)
inline

Multiplication of two transposition proxys holding objects of type TX and TY respectively.

Returns multiplication proxy XY.

See also
XY
Xtrans
operator*(TX const &, Xtrans<TY> const &)
operator*(Xtrans<TX> const &, TY const &)
operator*(TX const &, TY const &)

References getMachineEpsilon().

◆ operator*() [5/7]

Treal mat::operator* ( Xtrans< VectorGeneral< Treal, Tvector > > const &  xT,
VectorGeneral< Treal, Tvector > const &  y 
)

transpose(x) * y Scalar (dot) product of two vectors

References getMachineEpsilon(), and mat::VectorGeneral< Treal, Tvector >::getVector().

◆ operator*() [6/7]

template<typename TX , typename TY , typename TZ >
XYZ< TX, TY, TZ > mat::operator* ( XY< TX, TY > const &  AB,
TZ const &  CC 
)
inline

Multiplication of a multiplication proxy XY with an object of type TZ.

Returns multiplication proxy XYZ.

See also
XY
XYZ

References getMachineEpsilon().

◆ operator*() [7/7]

template<typename TX , typename TY , typename TZ >
XYZ< TX, TY, TZ > mat::operator* ( XY< TX, TY > const &  AB,
Xtrans< TZ > const &  trCC 
)
inline

Multiplication of a multiplication proxy XY with a transposition proxy Xtrans.

Returns multiplication proxy XYZ.

See also
XY
XYZ
Xtrans

References getMachineEpsilon().

◆ operator+() [1/2]

template<typename TX , typename TY >
XpY< TX, TY > mat::operator+ ( TX const &  AA,
TY const &  BB 
)
inline

Addition of two objects of type TX and TY.

See also
XpY

References getMachineEpsilon().

◆ operator+() [2/2]

XYZpUV< TX, TY, TZ, TU, TV > mat::operator+ ( XYZ< TX, TY, TZ > const &  ABC,
XY< TU, TV > const &  DE 
)
inline

Addition of two multiplication proxys XYZ and XY.

Returns multiplication and addition proxy XYZpUV.

See also
XY
XYZ
XYZpUV

References getMachineEpsilon().

◆ operator-()

template<typename TX , typename TY >
XmY< TX, TY > mat::operator- ( TX const &  AA,
TY const &  BB 
)
inline

Substraction of two objects of type TX and TY.

See also
XmY

References getMachineEpsilon().

◆ operator<<()

template<typename Treal >
std::ostream & mat::operator<< ( std::ostream &  s,
Interval< Treal > const &  in 
)

References getMachineEpsilon().

◆ packedtofull()

template<class Treal >
static void mat::packedtofull ( const Treal packed,
Treal full,
const int  size 
)
static

References getMachineEpsilon().

Referenced by mainFun().

◆ pocon()

template<class T >
static void mat::pocon ( const char uplo,
const int n,
const T A,
const int lda,
const T anorm,
T rcond,
T work,
int iwork,
int info 
)
inlinestatic

◆ pocon< double >()

template<>
void mat::pocon< double > ( const char uplo,
const int n,
const double A,
const int lda,
const double anorm,
double rcond,
double work,
int iwork,
int info 
)
inline

◆ pocon< float >()

template<>
void mat::pocon< float > ( const char uplo,
const int n,
const float A,
const int lda,
const float anorm,
float rcond,
float work,
int iwork,
int info 
)
inline

◆ potrf()

template<class T >
static void mat::potrf ( const char uplo,
const int n,
T A,
const int lda,
int info 
)
inlinestatic

◆ potrf< double >()

template<>
void mat::potrf< double > ( const char uplo,
const int n,
double A,
const int lda,
int info 
)
inline

◆ potrf< float >()

template<>
void mat::potrf< float > ( const char uplo,
const int n,
float A,
const int lda,
int info 
)
inline

◆ pptrf()

template<class T >
static void mat::pptrf ( const char uplo,
const int n,
T ap,
int info 
)
inlinestatic

References getMachineEpsilon(), and template_lapack_pptrf().

Referenced by mainFun().

◆ pptrf< double >()

template<>
void mat::pptrf< double > ( const char uplo,
const int n,
double ap,
int info 
)
inline

◆ pptrf< float >()

template<>
void mat::pptrf< float > ( const char uplo,
const int n,
float ap,
int info 
)
inline

◆ quicksort()

template<class Treal >
void mat::quicksort ( const Treal value,
int index,
int  low,
int  high 
)

References getMachineEpsilon(), and quicksort().

Referenced by quicksort().

◆ random_shuffle()

template<class RandomIt >
static void mat::random_shuffle ( RandomIt  first,
RandomIt  last 
)
static

References getMachineEpsilon().

Referenced by mainFun().

◆ read_matrix()

static void mat::read_matrix ( Tmatrix A,
char const *const  matrixPath,
int const  size 
)
static

References A, and getMachineEpsilon().

◆ read_sparse_matrix()

static void mat::read_sparse_matrix ( Tmatrix A,
char const *const  rowPath,
char const *const  colPath,
char const *const  valPath,
int const  nval 
)
static

◆ read_xyz()

template<typename Treal >
static void mat::read_xyz ( Treal x,
Treal y,
Treal z,
char atomsPath,
int const  natoms,
int const  size 
)
static

References getMachineEpsilon().

◆ scal()

template<class T >
static void mat::scal ( const int n,
const T da,
T dx,
const int incx 
)
inlinestatic

◆ scal< double >()

template<>
void mat::scal< double > ( const int n,
const double da,
double dx,
const int incx 
)
inline

◆ scal< float >()

template<>
void mat::scal< float > ( const int n,
const float da,
float dx,
const int incx 
)
inline

◆ sign()

template<typename Treal >
int mat::sign ( Treal  value)
inline

Sign function returns the sign of the input.

1 for positive, -1 for negative and 0 for zero.

Referenced by bisection().

◆ spgst()

template<class T >
static void mat::spgst ( const int itype,
const char uplo,
const int n,
T ap,
const T bp,
int info 
)
inlinestatic

◆ spgst< double >()

template<>
void mat::spgst< double > ( const int itype,
const char uplo,
const int n,
double ap,
const double bp,
int info 
)
inline

◆ spgst< float >()

template<>
void mat::spgst< float > ( const int itype,
const char uplo,
const int n,
float ap,
const float bp,
int info 
)
inline

◆ sqrtInt()

template<typename Treal >
Interval< Treal > mat::sqrtInt ( Interval< Treal > const &  other)

◆ stevr()

template<class T >
static void mat::stevr ( const char jobz,
const char range,
const int n,
T d,
T e,
const T vl,
const T vu,
const int il,
const int iu,
const T abstol,
int m,
T w,
T z,
const int ldz,
int isuppz,
T work,
int lwork,
int iwork,
int liwork,
int info 
)
inlinestatic

◆ stevr< double >()

template<>
void mat::stevr< double > ( const char jobz,
const char range,
const int n,
double d,
double e,
const double vl,
const double vu,
const int il,
const int iu,
const double abstol,
int m,
double w,
double z,
const int ldz,
int isuppz,
double work,
int lwork,
int iwork,
int liwork,
int info 
)
inline

◆ stevr< float >()

template<>
void mat::stevr< float > ( const char jobz,
const char range,
const int n,
float d,
float e,
const float vl,
const float vu,
const int il,
const int iu,
const float abstol,
int m,
float w,
float z,
const int ldz,
int isuppz,
float work,
int lwork,
int iwork,
int liwork,
int info 
)
inline

◆ stevx()

template<class T >
static void mat::stevx ( const char jobz,
const char range,
const int n,
T d,
T e,
const T vl,
const T vu,
const int il,
const int iu,
const T abstol,
int m,
T w,
T z,
const int ldz,
T work,
int iwork,
int ifail,
int info 
)
inlinestatic

◆ stevx< double >()

template<>
void mat::stevx< double > ( const char jobz,
const char range,
const int n,
double d,
double e,
const double vl,
const double vu,
const int il,
const int iu,
const double abstol,
int m,
double w,
double z,
const int ldz,
double work,
int iwork,
int ifail,
int info 
)
inline

◆ stevx< float >()

template<>
void mat::stevx< float > ( const char jobz,
const char range,
const int n,
float d,
float e,
const float vl,
const float vu,
const int il,
const int iu,
const float abstol,
int m,
float w,
float z,
const int ldz,
float work,
int iwork,
int ifail,
int info 
)
inline

◆ syev()

template<class T >
static void mat::syev ( const char jobz,
const char uplo,
const int n,
T a,
const int lda,
T w,
T work,
const int lwork,
int info 
)
inlinestatic

◆ syev< double >()

template<>
void mat::syev< double > ( const char jobz,
const char uplo,
const int n,
double a,
const int lda,
double w,
double work,
const int lwork,
int info 
)
inline

◆ syev< float >()

template<>
void mat::syev< float > ( const char jobz,
const char uplo,
const int n,
float a,
const int lda,
float w,
float work,
const int lwork,
int info 
)
inline

◆ sygv()

template<class T >
static void mat::sygv ( const int itype,
const char jobz,
const char uplo,
const int n,
T A,
const int lda,
T B,
const int ldb,
T w,
T work,
const int lwork,
int info 
)
inlinestatic

◆ sygv< double >()

template<>
void mat::sygv< double > ( const int itype,
const char jobz,
const char uplo,
const int n,
double A,
const int lda,
double B,
const int ldb,
double w,
double work,
const int lwork,
int info 
)
inline

◆ sygv< float >()

template<>
void mat::sygv< float > ( const int itype,
const char jobz,
const char uplo,
const int n,
float A,
const int lda,
float B,
const int ldb,
float w,
float work,
const int lwork,
int info 
)
inline

◆ symm()

template<class T >
static void mat::symm ( const char side,
const char uplo,
const int m,
const int n,
const T alpha,
const T A,
const int lda,
const T B,
const int ldb,
const T beta,
T C,
const int ldc 
)
inlinestatic

◆ symm< double >()

template<>
void mat::symm< double > ( const char side,
const char uplo,
const int m,
const int n,
const double alpha,
const double A,
const int lda,
const double B,
const int ldb,
const double beta,
double C,
const int ldc 
)
inline

◆ symm< float >()

template<>
void mat::symm< float > ( const char side,
const char uplo,
const int m,
const int n,
const float alpha,
const float A,
const int lda,
const float B,
const int ldb,
const float beta,
float C,
const int ldc 
)
inline

◆ symv()

template<class T >
static void mat::symv ( const char uplo,
const int n,
const T alpha,
const T A,
const int lda,
const T x,
const int incx,
const T beta,
T y,
const int incy 
)
inlinestatic

◆ symv< double >()

template<>
void mat::symv< double > ( const char uplo,
const int n,
const double alpha,
const double A,
const int lda,
const double x,
const int incx,
const double beta,
double y,
const int incy 
)
inline

◆ symv< float >()

template<>
void mat::symv< float > ( const char uplo,
const int n,
const float alpha,
const float A,
const int lda,
const float x,
const int incx,
const float beta,
float y,
const int incy 
)
inline

◆ syrk()

template<class T >
static void mat::syrk ( const char uplo,
const char trans,
const int n,
const int k,
const T alpha,
const T A,
const int lda,
const T beta,
T C,
const int ldc 
)
inlinestatic

◆ syrk< double >()

template<>
void mat::syrk< double > ( const char uplo,
const char trans,
const int n,
const int k,
const double alpha,
const double A,
const int lda,
const double beta,
double C,
const int ldc 
)
inline

◆ syrk< float >()

template<>
void mat::syrk< float > ( const char uplo,
const char trans,
const int n,
const int k,
const float alpha,
const float A,
const int lda,
const float beta,
float C,
const int ldc 
)
inline

◆ tptri()

template<class T >
static void mat::tptri ( const char uplo,
const char diag,
const int n,
T ap,
int info 
)
inlinestatic

References getMachineEpsilon(), and template_lapack_tptri().

Referenced by mainFun().

◆ tptri< double >()

template<>
void mat::tptri< double > ( const char uplo,
const char diag,
const int n,
double ap,
int info 
)
inline

◆ tptri< float >()

template<>
void mat::tptri< float > ( const char uplo,
const char diag,
const int n,
float ap,
int info 
)
inline

◆ trace()

◆ transpose() [1/2]

template<typename TX >
Xtrans< TX > mat::transpose ( const Xtrans< TX > &  xtrans)
inline

Transposition.

Returns a transposition proxy of an object of type Xtrans<TX>. Only for correct treatment of repeated transposition, e.g. transpose(transpose(A))

See also
Xtrans
transpose(TX const &)

References getMachineEpsilon().

◆ transpose() [2/2]

◆ trifulltofull()

template<class Treal >
static void mat::trifulltofull ( Treal trifull,
const int  size 
)
static

◆ tripackedtofull()

template<class Treal >
static void mat::tripackedtofull ( const Treal packed,
Treal full,
const int  size 
)
static

References getMachineEpsilon().

Referenced by mainFun().

◆ trmm()

template<class T >
static void mat::trmm ( const char side,
const char uplo,
const char transa,
const char diag,
const int m,
const int n,
const T alpha,
const T A,
const int lda,
T B,
const int ldb 
)
inlinestatic

◆ trmm< double >()

template<>
void mat::trmm< double > ( const char side,
const char uplo,
const char transa,
const char diag,
const int m,
const int n,
const double alpha,
const double A,
const int lda,
double B,
const int ldb 
)
inline

◆ trmm< float >()

template<>
void mat::trmm< float > ( const char side,
const char uplo,
const char transa,
const char diag,
const int m,
const int n,
const float alpha,
const float A,
const int lda,
float B,
const int ldb 
)
inline

◆ trmv()

template<class T >
static void mat::trmv ( const char uplo,
const char trans,
const char diag,
const int n,
const T A,
const int lda,
T x,
const int incx 
)
inlinestatic

◆ trmv< double >()

template<>
void mat::trmv< double > ( const char uplo,
const char trans,
const char diag,
const int n,
const double A,
const int lda,
double x,
const int incx 
)
inline

◆ trmv< float >()

template<>
void mat::trmv< float > ( const char uplo,
const char trans,
const char diag,
const int n,
const float A,
const int lda,
float x,
const int incx 
)
inline

◆ trtri()

template<class T >
static void mat::trtri ( const char uplo,
const char diag,
const int n,
T A,
const int lda,
int info 
)
inlinestatic

◆ trtri< double >()

template<>
void mat::trtri< double > ( const char uplo,
const char diag,
const int n,
double A,
const int lda,
int info 
)
inline

◆ trtri< float >()

template<>
void mat::trtri< float > ( const char uplo,
const char diag,
const int n,
float A,
const int lda,
int info 
)
inline