dbandedsymmatrix_functions.h

Go to the documentation of this file.

#ifndef DBANDEDSYMMATRIX_FUNCTIONS
#define DBANDEDSYMMATRIX_FUNCTIONS

#include "lapack_functions.h"

void DoubleBandedSymmetricMatrixVectorProduct(const tLVectorIndex& nX,
                                            const tLVectorIncrement& incX,
                                            const double* x,
                                            const tLVectorIndex& nRows,//number of viewed rows of A
                                            const tLVectorIndex& nBands,//number of sub diagonals of A
                                            const double* A,
                                            const lapack_real& alpha,
                                            const lapack_real& beta,
                                            const tLVectorIndex& nY,
                                            const tLVectorIncrement& incY,
                                            double* y);



//C=alpha*S.B+beta.C or C:=alpha.B.S+beta C where
// - S is band symetric and
// - B is full matrix and
// - C is full matrix
void DoubleBandedSymmetricMatrixMatrixProduct(const tBoolean& isLeftSide,
                                            const tLVectorIndex& nSRows,
                                            const tLVectorIndex& nSBands,
                                            const double* S,
                                            const tLVectorIndex& nBRows,
                                            const tLVectorIndex& nBCols,
                                            const tLVectorIncrement& incB,
                                            const tLVectorIndex& ldB,
                                            const double* B,
                                            const lapack_real& alpha,
                                            const lapack_real& beta,
                                            const tLVectorIndex& nCRows,
                                            const tLVectorIndex& nCCols,
                                            const tLVectorIncrement& incC,
                                            const tLVectorIndex& ldC,
                                            double *C);

tBoolean DoubleBandedSymmetricMatrixEigenValues(const tLVectorIndex& nRowsA,
                                              const tLVectorIndex& nBandsA,
                                              double* A,
                                              const tLVectorIndex& nU,
                                              double* U,
                                              const tLVectorIndex& ldW,
                                              double* W);

tBoolean DoubleBandedSymmetricMatrixBoundedEigenValues(const tLVectorIndex& nRowsA,
                                                     const tLVectorIndex& nBandsA,
                                                     double* A,
                                                     const tLVectorIndex& nU,//of size M=N  or M=iu-il+1
                                                     double* U,
                                                     const tLVectorIndex& minIndex,
                                                     const tLVectorIndex& maxIndex,
                                                     const double& minBound,
                                                     const double& maxBound,
                                                     const tLVectorIndex& ldQ,
                                                     double *Q,//of size nRowsA x nRowsA
                                                     const tLVectorIndex& ldW,//Z of size nRows  x M
                                                     double *W);
    
  
    
tBoolean DoubleBandedSymmetricMatrixCholeskyFactorization(const tLVectorIndex& nRows,
                                                        const tLVectorIndex& nBands,
                                                        double* A);

#endif