FMatrix.h

Go to the documentation of this file.

 
#ifndef FMatrix_H
#define FMatrix_H
 
#include <stdio.h>
#include <stdlib.h>
#include <math.h>
#include <assert.h>
#include <string.h>
#include <stdbool.h>
 
#include "ludec.h"
#include "DataTypes.h"
 
 
class FMatrix {
    public:
  int     mxNCols=0;         
  int     mxNRows=0;         
  int     NCols=0;           
  int     NRows=0;           
  int     lsize=0;           
  Real  **Ent=NULL;          
  bool    decomposed=false;  
  int    *indx=NULL;         
  Real  **base1=NULL;        
  void     LUD( );
  void     LUsolve( Real * );
  FMatrix& operator=(FMatrix &);
  template<class T>
  FMatrix( T , T );
  ~FMatrix( );
  void Zero();
  FMatrix *Transpose(bool);
  FMatrix *Invert();
  void     resize( int, int );
  FMatrix *Inv();
};
 
#ifdef __cplusplus
extern "C"
{
#endif
 
  void     Destroy_FMatrix(FMatrix *M );
  void     Zero_FMatrix(FMatrix *M );
  Real     FMatrix_Trace(FMatrix *M);
 
  FMatrix *FM_Rotate(FMatrix *R, double angle, int ax);
  Real     FMatrix_Det( FMatrix * );
  void     FM_Times_FM(FMatrix *A, FMatrix *B, FMatrix *C);
  void     FM_Times_FMinArray(FMatrix *A, FMatrix *B, Real *C);
  void     TransFM_Times_FM(FMatrix *A, FMatrix *B, FMatrix *C);
  void     FM_Times_TransFM(FMatrix *A, FMatrix *B, FMatrix *C);
  void     FM_Times_Vec(Real dx, FMatrix *B, Real*C, Real *result);
  void     B_C_BTrans(Real dx, FMatrix *B, Real *C, FMatrix *result);
  void     FM_Plus_FM(FMatrix *A, FMatrix *B, FMatrix *C);
  void     FM_Times_Scalar(FMatrix *A,Real a);
  Real*    FM_X_Vec( FMatrix *A, Real *b, Real *c );
  void     Print_FMatrix(FMatrix *A);
 
#ifdef __cplusplus
}
 
template<class T>
FMatrix::FMatrix(T NR, T NC)
{
  mxNRows = NRows = NR;
  mxNCols = NCols = NC;
 
  Ent = (Real**) calloc(mxNRows, sizeof(Real*));
  for (int i = 0; i < mxNRows; i++)
  {
    Ent[i] = (Real*) calloc(mxNCols, sizeof(Real));
  }
}
#endif
#endif