dbandeduppermatrix_functions.h

Go to the documentation of this file.

#ifndef DBANDEDUPPERMATRIX_FUNCTIONS
#define DBANDEDUPPERMATRIX_FUNCTIONS

#include "lapack_functions.h"


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

//B := alpha*op( A )*B, or B := alpha*B*op( A ),
void DoubleBandedUpperMatrixMatrixProduct( const tBoolean& isLeft,
                                         const tBoolean& isTransA,
                                         const tLVectorIndex& nABands,
                                         const double *A,//upper packed matrix
                                         const double& alpha,
                                         const tLVectorIndex& nBRows,
                                         const tLVectorIndex &nBCols,
                                         const tLVectorIncrement &incB,
                                         const tLVectorIndex &ldB,
                                         double* B);

//solve op(A).X=B and B:=X where X & B are  vectors 
void DoubleBandedUpperMatrixVectorSolve( const tBoolean& isTransA,
                                       const tLVectorIndex& nARows,
                                       const tLVectorIndex& nABands,
                                       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 DoubleBandedUpperMatrixSolve( const tBoolean& isTransA,
                                     const tLVectorIndex& nARows,
                                     const tLVectorIndex& nABands,
                                     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);

                                        
#endif