dpackeduppermatrix_functions.h

Go to the documentation of this file.

#ifndef DPACKEDUPPERMATRIX_FUNCTIONS
#define DPACKEDUPPERMATRIX_FUNCTIONS

#include "lapack_functions.h"


//X:=A.X
//or
//X:=tA.X
void DoublePackedUpperMatrixVectorProduct( const tBoolean& isTransA,
                                           const tLVectorIndex& nRowsA,//number of viewed rows of A
                                           const double* A,
                                           const tLVectorIndex& nX,
                                           const tLVectorIncrement& incX,
                                           double* x);

//solve op(A).X=B and B:=X where X & B are  vectors 
void DoublePackedUpperMatrixVectorSolve( const tBoolean& isTransA,
                                         const tLVectorIndex& nARows,
                                         const double *A,//upper packed matrix
                                         const tLVectorIndex& nB,
                                         const tLVectorIncrement &incB,
                                         double* B);
//solve op(A).X=B and B:=X where X & B are  matrices of size nARows x nBcols  
tBoolean DoublePackedUpperMatrixSolve( const tBoolean& isTransA,
                                       const tLVectorIndex& nARows,
                                       const double *A,//upper matrix
                                       const tLVectorIndex& nBCols,//number of columns  of B
                                       const tLVectorIndex& ldB,//memory distance between B[i][j] and B[i][j+1]
                                       double* B);

//B:=alpha.op(A).B
void DoublePackedUpperMatrixMatrixProduct(const tBoolean& isLeft,
                                          const tBoolean& isTransA,
                                          const double* A,
                                          const double& alpha,
                                          const tLVectorIndex& nRowsB,const tLVectorIndex& nColsB,
                                          const tLVectorIncrement& incB,const tLVectorIndex& ldB,
                                          double* B);

//A:=A^-1 inverse the matrix of nRows number of rows 
tBoolean DoublePackedUpperMatrixInverse(const tLVectorIndex& nRows,
                                        double* A);
                                        
#endif