MatrixTools.hh

Go to the documentation of this file.

#ifndef _MATRIX_TOOLS_H_
#define _MATRIX_TOOLS_H_

#include <fstream.h>
#include "ocsTools.hh"

// Written by Ron Holzloehner. All Matrices are assumed to be square.



void CopyMatrix(double ** LHS_Matrix, double ** RHS_Matrix, int dim);


void DenseMatrixProduct(double **A, double **B, double **R,
                        bool Transpose_B, bool R_is_symmetric, int dim);



void DenseMatrixSimilarityProduct(double **A, double **B, double **R,
                                  double ** Work, bool Transpose_A, 
                                  bool R_is_symmetric,int dim);


void DenseMatrixScale(double **A, int dim, double scale);

 
void CopyLowerToUpperHalf(double **A, int dim);


void SymmetrizeMatrix(double **A, int dim);


void TransposeMatrix(double **A, int dim);


void ZeroMatrix(double **A, int dim);


void DenseMatrixFilter(double **A, int dim, double *filter);


void WriteMatrix(double **A, int dim, ofstream *F, double scaler);


void ReadMatrix(double **A, int dim, ifstream *F, double scaling);


void DenseMatrixVectorProduct(double * Out, double **A, double * In, int dim);



void DenseVectorMatrixProduct(double *Out, double *In, double **A, int dim);


void AddToPrincipalDiag(double **A, int dim, double add);


void AddVectorToPrincipalDiag(double **A, double *V, 
                              double ScalarMultiple,int dim);


double ScalarProductNormalized(cplx *aa, cplx *bb, int len);


double VectorAngle(cplx *aa, cplx *bb, int len);



void SchmidtOrthogonalizer(cplx **vec, int num, int len,
                           double *alpha);

void SchmidtOrthogonalizer(cplx **vec, int num, int len, 
                           double *alpha, int k1);




double InfNormMatrix(double **A,int dim);


double OneNormMatrix(double **A,int dim);


double InnerProduct(double * V, double *W, int Length);

#endif // _MATRIX_TOOLS_H_