NmMatrix.h

Go to the documentation of this file.

#pragma once
 
//
// (C) Copyright 2005-2024 by Graebert GmbH.
//
// Permission to use, copy, modify, and distribute this software in
// object code form for any purpose and without fee is hereby granted,
// provided that the above copyright notice appears in all copies and
// that both that copyright notice and the limited warranty and
// restricted rights notice below appear in all supporting
// documentation.
//
// GRAEBERT PROVIDES THIS PROGRAM "AS IS" AND WITH ALL FAULTS.
// GRAEBERT SPECIFICALLY DISCLAIMS ANY IMPLIED WARRANTY OF
// MERCHANTABILITY OR FITNESS FOR A PARTICULAR USE.  GRAEBERT GMBH
// DOES NOT WARRANT THAT THE OPERATION OF THE PROGRAM WILL BE
// UNINTERRUPTED OR ERROR FREE.
 
class NmIntVec;
class NmVector;
 
class ARX_API NmMatrix
{
public:
          NmMatrix( NmMatrix const &);
          NmMatrix( NmVector const &);
          NmMatrix(int);
          NmMatrix(int,double);
          NmMatrix(int,int);
          NmMatrix(void);
          ~NmMatrix(void);
          NmMatrix &  invert(double);
          NmMatrix &  invertLUs(double,bool &);
          NmMatrix &  operator*=( NmMatrix const &);
          NmMatrix &  operator*=(double);
          NmMatrix &  operator+=( NmMatrix const &);
          NmMatrix &  operator-=( NmMatrix const &);
          NmMatrix &  operator/=(double);
          NmMatrix &  operator=( NmMatrix const &);
          NmMatrix &  operator=( NmVector const &);
          NmMatrix &  operator=(double);
          NmMatrix &  reset(int);
          NmMatrix &  reset(int,int);
          NmMatrix &  reset(int,int,double,double);
          NmMatrix &  setColumn(int, NmVector const &);
          NmMatrix &  setDiagonal(double);
          NmMatrix &  setRow(int, NmVector const &);
          NmMatrix &  transpose(void);
          NmMatrix  operator*( NmMatrix const &)const ;
          NmMatrix  operator+( NmMatrix const &)const ;
          NmMatrix  operator-( NmMatrix const &)const ;
          NmVector  column(int)const ;
          NmVector  diagonal(void)const ;
          NmVector  operator*( NmVector const &)const ;
          NmVector  row(int)const ;
         double &  operator()(int,int);
         double  determinant(double);
         double  operator()(int,int)const ;
         int  LUdecompose( NmIntVec &);
         int  colDim(void)const ;
         int  decomposeLU( NmIntVec &,double);
         int  decomposeLUScaling( NmIntVec &,double &,double);
         int  eigensystem( NmVector &);
         int  exequal( NmVector &, NmIntVec &, NmIntVec &, NmVector &, NmIntVec &, NmIntVec &,double &,double);
         int  feasible( NmVector &, NmIntVec &, NmIntVec &, NmVector &, NmIntVec &, NmIntVec &,double &,double);
         int  feasible( NmVector &, NmIntVec &, NmIntVec &,double);
         int  optimal( NmVector &, NmIntVec &, NmIntVec &, NmVector &, NmIntVec &, NmIntVec &,double &,double);
         int  rowDim(void)const ;
         int  solve( NmVector &,double);
         int  symmetricLU( NmVector &,double);
         void  dim(int &,int &)const ;
         void  force( NmVector &,double);
         void  solveAfterLUdecompose( NmVector &, NmIntVec const &)const ;
         void  solveLU( NmIntVec const &, NmVector &)const ;
         void  solveLUs( NmIntVec const &, NmVector &)const ;
 
private:
         int  eigen( NmVector &, NmVector &);
         int  exequal( NmVector &, NmIntVec &, NmIntVec &,double);
         int  isUpperTriangular(void);
         void  band( NmVector &, NmVector &);
         void  pivot( NmVector &, NmIntVec &, NmIntVec &, NmVector &, NmIntVec &, NmIntVec &,double &,int,int,double);
         void  pivot( NmVector &, NmIntVec &, NmIntVec &,int,int,double);
 
};