// ------------------------------------------------------------
// mutil.h
// Matrix utilities
// 5.3.2005, J.Ahrensfeld
//
// ------------------------------------------------------------
#ifndef	MUTIL_H
#define	MUTIL_H

// ------------------------------------------------------------
// ------------------------------------------------------------
void RandMat(double *pA, int N);
void IncMat(double *pA, int nRow, int nCol, char order);
void CopyTrans(double *pA, double *pAT, int nRow, int nCol);
void Eye(double *pI, double value, int N);
void Diag(double *pA, double *pY, int N);
void MakeDiag(double *pD, double *pX, int N);
void PrintMatrix(char *pStr, double *A, int num_row, int num_col);
double smin(double v1, double v2);
double smax(double v1, double v2);
void vsum(double *pX, double *pY, int lda, int N);
void GetPrincipleEigenVector(double *pA, double *pV, int lda);
double dnrm(int N, int M, double *A, char mode);

// ------------------------------------------------------------
#endif //MUTIL_H