GeneralMatrix Class Reference

#include <newmat.h>

Inheritance diagram for GeneralMatrix:

Public Member Functions
GeneralMatrix *	Evaluate (MatrixType mt=MatrixTypeUnSp)
virtual MatrixType	Type () const =0
int	Nrows () const
int	Ncols () const
int	Storage () const
Real *	Store () const
virtual	~GeneralMatrix ()
void	tDelete ()
bool	reuse ()
void	Protect ()
int	Tag () const
bool	IsZero () const
void	Release ()
void	Release (int t)
void	ReleaseAndDelete ()
void	operator<< (const Real *)
void	operator<< (const int *)
void	operator<< (const BaseMatrix &X)
void	Inject (const GeneralMatrix &)
void	operator+= (const BaseMatrix &)
void	operator-= (const BaseMatrix &)
void	operator*= (const BaseMatrix &)
void	operator|= (const BaseMatrix &)
void	operator &= (const BaseMatrix &)
void	operator+= (Real)
void	operator-= (Real r)
void	operator*= (Real)
void	operator/= (Real r)
virtual GeneralMatrix *	MakeSolver ()
virtual void	Solver (MatrixColX &, const MatrixColX &)
virtual void	GetRow (MatrixRowCol &)=0
virtual void	RestoreRow (MatrixRowCol &)
virtual void	NextRow (MatrixRowCol &)
virtual void	GetCol (MatrixRowCol &)=0
virtual void	GetCol (MatrixColX &)=0
virtual void	RestoreCol (MatrixRowCol &)
virtual void	RestoreCol (MatrixColX &)
virtual void	NextCol (MatrixRowCol &)
virtual void	NextCol (MatrixColX &)
Real	SumSquare () const
Real	SumAbsoluteValue () const
Real	Sum () const
Real	MaximumAbsoluteValue1 (int &i) const
Real	MinimumAbsoluteValue1 (int &i) const
Real	Maximum1 (int &i) const
Real	Minimum1 (int &i) const
Real	MaximumAbsoluteValue () const
Real	MaximumAbsoluteValue2 (int &i, int &j) const
Real	MinimumAbsoluteValue () const
Real	MinimumAbsoluteValue2 (int &i, int &j) const
Real	Maximum () const
Real	Maximum2 (int &i, int &j) const
Real	Minimum () const
Real	Minimum2 (int &i, int &j) const
LogAndSign	LogDeterminant () const
virtual bool	IsEqual (const GeneralMatrix &) const
void	CheckStore () const
virtual void	SetParameters (const GeneralMatrix *)
	operator ReturnMatrix () const
ReturnMatrix	ForReturn () const
virtual bool	SameStorageType (const GeneralMatrix &A) const
virtual void	ReSizeForAdd (const GeneralMatrix &A, const GeneralMatrix &B)
virtual void	ReSizeForSP (const GeneralMatrix &A, const GeneralMatrix &B)
virtual void	ReSize (const GeneralMatrix &A)
MatrixInput	operator<< (Real)
MatrixInput	operator<< (int f)
void	CleanUp ()
Public Member Functions inherited from BaseMatrix
AddedMatrix	operator+ (const BaseMatrix &) const
MultipliedMatrix	operator* (const BaseMatrix &) const
SubtractedMatrix	operator- (const BaseMatrix &) const
ConcatenatedMatrix	operator| (const BaseMatrix &) const
StackedMatrix	operator & (const BaseMatrix &) const
ShiftedMatrix	operator+ (Real) const
ScaledMatrix	operator* (Real) const
ScaledMatrix	operator/ (Real) const
ShiftedMatrix	operator- (Real) const
TransposedMatrix	t () const
NegatedMatrix	operator- () const
ReversedMatrix	Reverse () const
InvertedMatrix	i () const
RowedMatrix	AsRow () const
ColedMatrix	AsColumn () const
DiagedMatrix	AsDiagonal () const
MatedMatrix	AsMatrix (int, int) const
GetSubMatrix	SubMatrix (int, int, int, int) const
GetSubMatrix	SymSubMatrix (int, int) const
GetSubMatrix	Row (int) const
GetSubMatrix	Rows (int, int) const
GetSubMatrix	Column (int) const
GetSubMatrix	Columns (int, int) const
Real	AsScalar () const
Real	Determinant () const
Real	NormFrobenius () const
virtual Real	Trace () const
Real	Norm1 () const
Real	NormInfinity () const
virtual MatrixBandWidth	BandWidth () const
void	IEQND () const

Protected Member Functions
	GeneralMatrix ()
	GeneralMatrix (ArrayLengthSpecifier)
void	Add (GeneralMatrix *, Real)
void	Add (Real)
void	NegAdd (GeneralMatrix *, Real)
void	NegAdd (Real)
void	Multiply (GeneralMatrix *, Real)
void	Multiply (Real)
void	Negate (GeneralMatrix *)
void	Negate ()
void	ReverseElements ()
void	ReverseElements (GeneralMatrix *)
void	operator= (Real)
Real *	GetStore ()
GeneralMatrix *	BorrowStore (GeneralMatrix *, MatrixType)
void	GetMatrix (const GeneralMatrix *)
void	Eq (const BaseMatrix &, MatrixType)
void	Eq (const GeneralMatrix &)
void	Eq (const BaseMatrix &, MatrixType, bool)
void	Eq2 (const BaseMatrix &, MatrixType)
int	search (const BaseMatrix *) const
virtual GeneralMatrix *	Transpose (TransposedMatrix *, MatrixType)
void	CheckConversion (const BaseMatrix &)
void	ReSize (int, int, int)
virtual short	SimpleAddOK (const GeneralMatrix *)
virtual void	MiniCleanUp ()
void	PlusEqual (const GeneralMatrix &gm)
void	MinusEqual (const GeneralMatrix &gm)
void	PlusEqual (Real f)
void	MinusEqual (Real f)

Protected Attributes
int	tag
int	nrows
int	ncols
int	storage
Real *	store

Private Member Functions
virtual GeneralMatrix *	Image () const

Friends
class	Matrix
class	SquareMatrix
class	nricMatrix
class	SymmetricMatrix
class	UpperTriangularMatrix
class	LowerTriangularMatrix
class	DiagonalMatrix
class	CroutMatrix
class	RowVector
class	ColumnVector
class	BandMatrix
class	LowerBandMatrix
class	UpperBandMatrix
class	SymmetricBandMatrix
class	BaseMatrix
class	AddedMatrix
class	MultipliedMatrix
class	SubtractedMatrix
class	SPMatrix
class	KPMatrix
class	ConcatenatedMatrix
class	StackedMatrix
class	SolvedMatrix
class	ShiftedMatrix
class	NegShiftedMatrix
class	ScaledMatrix
class	TransposedMatrix
class	ReversedMatrix
class	NegatedMatrix
class	InvertedMatrix
class	RowedMatrix
class	ColedMatrix
class	DiagedMatrix
class	MatedMatrix
class	GetSubMatrix
class	ReturnMatrix
class	LinearEquationSolver
class	GenericMatrix

Detailed Description

Definition at line 386 of file newmat.h.

Constructor & Destructor Documentation

GeneralMatrix() [1/2]

GeneralMatrix::GeneralMatrix ( )

protected

Definition at line 33 of file newmat4.cpp.

{ store=0; storage=0; nrows=0; ncols=0; tag=-1; }

GeneralMatrix() [2/2]

GeneralMatrix::GeneralMatrix ( ArrayLengthSpecifier  s )

protected

Definition at line 36 of file newmat4.cpp.

References MatrixErrorNoSpace(), MONITOR_REAL_NEW, REPORT, and ArrayLengthSpecifier::Value().

{
   REPORT
   storage=s.Value(); tag=-1;
   if (storage)
   {
      store = new Real [storage]; MatrixErrorNoSpace(store);
      MONITOR_REAL_NEW("Make (GenMatrix)",storage,store)
   }
   else store = 0;
}

~GeneralMatrix()

GeneralMatrix::~GeneralMatrix ( )

virtual

Definition at line 161 of file newmat4.cpp.

References MONITOR_REAL_DELETE, storage, and store.

{
   if (store)
   {
      MONITOR_REAL_DELETE("Free (GenMatrix)",storage,store)
      delete [] store;
   }
}

Member Function Documentation

Add() [1/2]

void GeneralMatrix::Add ( GeneralMatrix *  gm1, Real  f  )

protected

Definition at line 311 of file newmat5.cpp.

References REPORT, and store.

Referenced by Matrix::operator+=(), SquareMatrix::operator+=(), nricMatrix::operator+=(), SymmetricMatrix::operator+=(), RowVector::operator+=(), ColumnVector::operator+=(), Matrix::operator-=(), SquareMatrix::operator-=(), nricMatrix::operator-=(), SymmetricMatrix::operator-=(), RowVector::operator-=(), and ColumnVector::operator-=().

{
   REPORT
   Real* s1=gm1->store; Real* s=store; int i=(storage >> 2);
   while (i--)
   { *s++ = *s1++ + f; *s++ = *s1++ + f; *s++ = *s1++ + f; *s++ = *s1++ + f; }
   i = storage & 3; while (i--) *s++ = *s1++ + f;
}

Add() [2/2]

void GeneralMatrix::Add ( Real  f )

protected

Definition at line 320 of file newmat5.cpp.

References REPORT.

{
   REPORT
   Real* s=store; int i=(storage >> 2);
   while (i--) { *s++ += f; *s++ += f; *s++ += f; *s++ += f; }
   i = storage & 3; while (i--) *s++ += f;
}

BorrowStore()

GeneralMatrix * GeneralMatrix::BorrowStore ( GeneralMatrix *  gmx, MatrixType  mt  )

protected

Definition at line 617 of file newmat4.cpp.

References Compare(), Evaluate(), GetStore(), REPORT, store, tag, tDelete(), and Type().

{
   if (!mt)
   {
      if (tag == -1) { REPORT gmx->tag = -2; gmx->store = store; }
      else { REPORT gmx->tag = 0; gmx->store = GetStore(); }
   }
   else if (Compare(gmx->Type(),mt))
   { REPORT  gmx->tag = 0; gmx->store = GetStore(); }
   else
   {
      REPORT gmx->tag = -2; gmx->store = store;
      gmx = gmx->Evaluate(mt); gmx->tag = 0; tDelete();
   }

   return gmx;
}

CheckConversion()

void GeneralMatrix::CheckConversion ( const BaseMatrix &  )

protected

CheckStore()

void GeneralMatrix::CheckStore

(

)

const

Definition at line 864 of file newmat4.cpp.

References REPORT, and store.

{
   REPORT
   if (!store)
      Throw(ProgramException("NRIC accessing matrix with unset dimensions"));
}

CleanUp()

void GeneralMatrix::CleanUp ( )

virtual

Reimplemented from BaseMatrix.

Reimplemented in BandLUMatrix, CroutMatrix, ColumnVector, RowVector, and nricMatrix.

Definition at line 874 of file newmat4.cpp.

References MONITOR_REAL_DELETE, ncols, nrows, REPORT, storage, store, and tag.

Referenced by nricMatrix::CleanUp(), RowVector::CleanUp(), ColumnVector::CleanUp(), CroutMatrix::CleanUp(), and BandLUMatrix::CleanUp().

{
   // set matrix dimensions to zero, delete storage
   REPORT
   if (store && storage)
   {
      MONITOR_REAL_DELETE("Free (CleanUp)    ",storage,store)
      REPORT delete [] store;
   }
   store=0; storage=0; nrows=0; ncols=0; tag = -1;
}

Eq() [1/3]

void GeneralMatrix::Eq ( const BaseMatrix &  X, MatrixType  mt  )

protected

Definition at line 637 of file newmat4.cpp.

References Evaluate(), GetMatrix(), MONITOR_REAL_DELETE, Protect(), Release(), REPORT, BaseMatrix::search(), storage, and store.

Referenced by Eq().

{
#ifdef DO_SEARCH
   int counter=X.search(this);
   if (counter==0)
   {
      REPORT
      if (store)
      {
         MONITOR_REAL_DELETE("Free (operator=)",storage,store)
         REPORT delete [] store; storage = 0; store = 0;
      }
   }
   else { REPORT Release(counter); }
   GeneralMatrix* gmx = ((BaseMatrix&)X).Evaluate(mt);
   if (gmx!=this) { REPORT GetMatrix(gmx); }
   else { REPORT }
   Protect();
#else
   GeneralMatrix* gmx = ((BaseMatrix&)X).Evaluate(mt);
   if (gmx!=this)
   {
      REPORT
      if (store)
      {
         MONITOR_REAL_DELETE("Free (operator=)",storage,store)
         REPORT delete [] store; storage = 0; store = 0;
      }
      GetMatrix(gmx);
   }
   else { REPORT }
   Protect();
#endif
}

Eq() [2/3]

void GeneralMatrix::Eq ( const GeneralMatrix &  X )

protected

Definition at line 677 of file newmat4.cpp.

References GetMatrix(), MONITOR_REAL_DELETE, Protect(), REPORT, storage, and store.

{
   GeneralMatrix* gmx = (GeneralMatrix*)&X;
   if (gmx!=this)
   {
      REPORT
      if (store)
      {
         MONITOR_REAL_DELETE("Free (operator=)",storage,store)
         REPORT delete [] store; storage = 0; store = 0;
      }
      GetMatrix(gmx);
   }
   else { REPORT }
   Protect();
}

Eq() [3/3]

void GeneralMatrix::Eq ( const BaseMatrix &  X, MatrixType  mt, bool  ldok  )

protected

Definition at line 695 of file newmat4.cpp.

References Eq(), REPORT, and MatrixType::SetDataLossOK().

{
   REPORT
   if (ldok) mt.SetDataLossOK();
   Eq(X, mt);
}

Eq2()

void GeneralMatrix::Eq2 ( const BaseMatrix &  X, MatrixType  mt  )

protected

Definition at line 702 of file newmat4.cpp.

References Evaluate(), GetMatrix(), Protect(), and REPORT.

{
   GeneralMatrix* gmx = ((BaseMatrix&)X).Evaluate(mt);
   if (gmx!=this) { REPORT GetMatrix(gmx); }  // simplify GetMatrix ?
   else { REPORT }
   Protect();
}

Evaluate()

GeneralMatrix * GeneralMatrix::Evaluate ( MatrixType  mt = MatrixTypeUnSp )

virtual

Implements BaseMatrix.

Definition at line 79 of file newmat5.cpp.

References Compare(), MatrixRowCol::Copy(), DirectPart, LoadOnEntry, MatrixType::New(), MatrixRow::Next(), REPORT, and StoreOnExit.

Referenced by BaseMatrix::AsScalar(), BandLUMatrix::BandLUMatrix(), BandMatrix::BandMatrix(), BorrowStore(), CroutMatrix::CroutMatrix(), BaseMatrix::Determinant(), DiagonalMatrix::DiagonalMatrix(), Eq(), Eq2(), AddedMatrix::Evaluate(), SPMatrix::Evaluate(), ConcatenatedMatrix::Evaluate(), StackedMatrix::Evaluate(), SubtractedMatrix::Evaluate(), ShiftedMatrix::Evaluate(), NegShiftedMatrix::Evaluate(), ScaledMatrix::Evaluate(), NegatedMatrix::Evaluate(), ReversedMatrix::Evaluate(), InvertedMatrix::Evaluate(), RowedMatrix::Evaluate(), ColedMatrix::Evaluate(), DiagedMatrix::Evaluate(), ReturnMatrix::Evaluate(), GenericMatrix::GenericMatrix(), IdentityMatrix::IdentityMatrix(), IsZero(), LinearEquationSolver::LinearEquationSolver(), BaseMatrix::LogDeterminant(), LowerBandMatrix::LowerBandMatrix(), LowerTriangularMatrix::LowerTriangularMatrix(), Matrix::Matrix(), BaseMatrix::Maximum(), BaseMatrix::Maximum1(), BaseMatrix::Maximum2(), BaseMatrix::MaximumAbsoluteValue(), BaseMatrix::MaximumAbsoluteValue1(), BaseMatrix::MaximumAbsoluteValue2(), BaseMatrix::Minimum(), BaseMatrix::Minimum1(), BaseMatrix::Minimum2(), BaseMatrix::MinimumAbsoluteValue(), BaseMatrix::MinimumAbsoluteValue1(), BaseMatrix::MinimumAbsoluteValue2(), BaseMatrix::Norm1(), BaseMatrix::NormInfinity(), operator*=(), GenericMatrix::operator*=(), operator+=(), GenericMatrix::operator+=(), GetSubMatrix::operator+=(), operator-=(), GenericMatrix::operator-=(), GetSubMatrix::operator-=(), operator<<(), GetSubMatrix::operator<<(), GenericMatrix::operator=(), GetSubMatrix::operator=(), operator|=(), GenericMatrix::operator|=(), GetSubMatrix::SetUpLHS(), BaseMatrix::Sum(), BaseMatrix::SumAbsoluteValue(), BaseMatrix::SumSquare(), SymmetricBandMatrix::SymmetricBandMatrix(), SymmetricMatrix::SymmetricMatrix(), BaseMatrix::Trace(), UpperBandMatrix::UpperBandMatrix(), and UpperTriangularMatrix::UpperTriangularMatrix().

{
   if (Compare(this->Type(),mt)) { REPORT return this; }
   REPORT
   GeneralMatrix* gmx = mt.New(nrows,ncols,this);
   MatrixRow mr(this, LoadOnEntry);
   MatrixRow mrx(gmx, StoreOnExit+DirectPart);
   int i=nrows;
   while (i--) { mrx.Copy(mr); mrx.Next(); mr.Next(); }
   tDelete(); gmx->ReleaseAndDelete(); return gmx;
}

ForReturn()

ReturnMatrix GeneralMatrix::ForReturn

(

)

const

Definition at line 206 of file newmat4.cpp.

References Image(), REPORT, and ReturnMatrix.

Referenced by Cholesky().

{
   REPORT
   GeneralMatrix* gm = Image(); gm->ReleaseAndDelete();
   return ReturnMatrix(gm);
}

GetCol() [1/2]

virtual void GeneralMatrix::GetCol ( MatrixRowCol &  )

pure virtual

Implemented in IdentityMatrix, BandLUMatrix, SymmetricBandMatrix, BandMatrix, CroutMatrix, RowVector, DiagonalMatrix, LowerTriangularMatrix, UpperTriangularMatrix, SymmetricMatrix, and Matrix.

GetCol() [2/2]

virtual void GeneralMatrix::GetCol ( MatrixColX &  )

pure virtual

Implemented in IdentityMatrix, BandLUMatrix, SymmetricBandMatrix, BandMatrix, CroutMatrix, RowVector, DiagonalMatrix, LowerTriangularMatrix, UpperTriangularMatrix, SymmetricMatrix, and Matrix.

GetMatrix()

void GeneralMatrix::GetMatrix ( const GeneralMatrix *  gmx )

protected

Definition at line 610 of file newmat4.cpp.

References ncols, Ncols(), nrows, Nrows(), REPORT, SetParameters(), storage, store, and tag.

Referenced by CroutMatrix::CroutMatrix(), DiagonalMatrix::DiagonalMatrix(), Eq(), Eq2(), IdentityMatrix::IdentityMatrix(), LowerTriangularMatrix::LowerTriangularMatrix(), Matrix::Matrix(), SquareMatrix::SquareMatrix(), SymmetricMatrix::SymmetricMatrix(), and UpperTriangularMatrix::UpperTriangularMatrix().

{
   REPORT  tag=-1; nrows=gmx->Nrows(); ncols=gmx->Ncols();
   storage=gmx->storage; SetParameters(gmx);
   store=((GeneralMatrix*)gmx)->GetStore();
}

GetRow()

virtual void GeneralMatrix::GetRow ( MatrixRowCol &  )

pure virtual

Implemented in IdentityMatrix, BandLUMatrix, SymmetricBandMatrix, BandMatrix, CroutMatrix, DiagonalMatrix, LowerTriangularMatrix, UpperTriangularMatrix, SymmetricMatrix, and Matrix.

GetStore()

Real * GeneralMatrix::GetStore ( )

protected

Definition at line 587 of file newmat4.cpp.

References MatrixErrorNoSpace(), MONITOR_REAL_NEW, and storage.

Referenced by BorrowStore().

{
   if (tag<0 || tag>1)
   {
      Real* s;
      if (storage)
      {
         s = new Real [storage]; MatrixErrorNoSpace(s);
         MONITOR_REAL_NEW("Make  (GetStore)",storage,s)
         BlockCopy(storage, store, s);
      }
      else s = 0;
      if (tag > 1) { REPORT tag--; }
      else if (tag < -1) { REPORT store = 0; delete this; } // borrowed store
      else { REPORT }
      return s;
   }
   Real* s = store;                             // CleanUp - done later
   if (tag==0) { REPORT store = 0; delete this; }
   else { REPORT  MiniCleanUp(); }
   return s;
}

Image()

GeneralMatrix * GeneralMatrix::Image ( ) const

privatevirtual

Reimplemented in IdentityMatrix, SymmetricBandMatrix, LowerBandMatrix, UpperBandMatrix, BandMatrix, ColumnVector, RowVector, DiagonalMatrix, LowerTriangularMatrix, UpperTriangularMatrix, SymmetricMatrix, nricMatrix, SquareMatrix, and Matrix.

Definition at line 838 of file newmat4.cpp.

Referenced by ForReturn(), GenericMatrix::GenericMatrix(), LinearEquationSolver::LinearEquationSolver(), GenericMatrix::operator*=(), GenericMatrix::operator+=(), GenericMatrix::operator-=(), GenericMatrix::operator=(), and GenericMatrix::operator|=().

{
  //   bool dummy = true;
  //   if (dummy)                                   // get rid of warning message
  Throw(InternalException("Cannot apply Image to this matrix type"));
  return 0;
}

Inject()

void GeneralMatrix::Inject

(

const GeneralMatrix &

X

)

Definition at line 714 of file newmat4.cpp.

References DirectPart, BaseMatrix::i(), MatrixRowCol::Inject(), LoadOnEntry, ncols, MatrixRow::Next(), nrows, REPORT, and StoreOnExit.

{
   REPORT
   Tracer tr("Inject");
   if (nrows != X.nrows || ncols != X.ncols)
      Throw(IncompatibleDimensionsException());
   MatrixRow mr((GeneralMatrix*)&X, LoadOnEntry);
   MatrixRow mrx(this, LoadOnEntry+StoreOnExit+DirectPart);
   int i=nrows;
   while (i--) { mrx.Inject(mr); mrx.Next(); mr.Next(); }
}

IsEqual()

bool GeneralMatrix::IsEqual ( const GeneralMatrix &  A ) const

virtual

Reimplemented in BandLUMatrix, and CroutMatrix.

Definition at line 903 of file newmat7.cpp.

References A, and REPORT.

{
   Tracer tr("GeneralMatrix IsEqual");
   if (A.Type() != Type())                       // not same types
      { REPORT return false; }
   if (&A == this)                               // same matrix
      { REPORT  return true; }
   if (A.nrows != nrows || A.ncols != ncols)
                                                 // different dimensions
   { REPORT return false; }
   // is matrix storage the same - compare store
   REPORT
   return RealEqual(A.store,store,storage);
}

IsZero()

bool GeneralMatrix::IsZero

(

)

const

Definition at line 876 of file newmat7.cpp.

References REPORT.

Referenced by IsZero().

{
   REPORT
   Real* s=store; int i = storage >> 2;
   while (i--)
   {
      if (*s++) return false; if (*s++) return false;
      if (*s++) return false; if (*s++) return false;
   }
   i = storage & 3; while (i--) if (*s++) return false;
   return true;
}

LogDeterminant()

LogAndSign GeneralMatrix::LogDeterminant ( ) const

virtual

Reimplemented from BaseMatrix.

Reimplemented in IdentityMatrix, BandLUMatrix, SymmetricBandMatrix, LowerBandMatrix, UpperBandMatrix, BandMatrix, CroutMatrix, DiagonalMatrix, LowerTriangularMatrix, and UpperTriangularMatrix.

Definition at line 688 of file newmat8.cpp.

References REPORT.

{
   REPORT
   Tracer tr("LogDeterminant");
   if (nrows != ncols) Throw(NotSquareException(*this));
   CroutMatrix C(*this); return C.LogDeterminant();
}

MakeSolver()

GeneralMatrix * GeneralMatrix::MakeSolver ( )

virtual

Reimplemented in IdentityMatrix, BandLUMatrix, SymmetricBandMatrix, LowerBandMatrix, UpperBandMatrix, BandMatrix, CroutMatrix, DiagonalMatrix, LowerTriangularMatrix, UpperTriangularMatrix, and Matrix.

Definition at line 24 of file newmat7.cpp.

References MatrixErrorNoSpace(), and REPORT.

Referenced by LinearEquationSolver::LinearEquationSolver().

{
   REPORT
   GeneralMatrix* gm = new CroutMatrix(*this);
   MatrixErrorNoSpace(gm); gm->ReleaseAndDelete(); return gm;
}

Maximum()

Real GeneralMatrix::Maximum ( ) const

virtual

Reimplemented from BaseMatrix.

Reimplemented in SymmetricBandMatrix, and BandMatrix.

Definition at line 254 of file newmat8.cpp.

References REPORT.

{
   REPORT
   if (storage == 0) NullMatrixError(this);
   int l = storage - 1; Real* s = store; Real maxval = *s++;
   while (l--) { Real a = *s++; if (maxval < a) maxval = a; }
   ((GeneralMatrix&)*this).tDelete(); return maxval;
}

Maximum1()

Real GeneralMatrix::Maximum1 ( int &  i ) const

virtual

Reimplemented from BaseMatrix.

Definition at line 263 of file newmat8.cpp.

References REPORT.

Referenced by Matrix::Maximum2().

{
   REPORT
   if (storage == 0) NullMatrixError(this);
   int l = storage - 1; Real* s = store; Real maxval = *s++; int li = l;
   while (l--) { Real a = *s++; if (maxval <= a) { maxval = a; li = l; } }
   i = storage - li;
   ((GeneralMatrix&)*this).tDelete(); return maxval;
}

Maximum2()

Real GeneralMatrix::Maximum2 ( int &  i, int &  j  ) const

virtual

Reimplemented from BaseMatrix.

Reimplemented in Matrix.

Definition at line 322 of file newmat8.cpp.

References REPORT.

Referenced by BandMatrix::Maximum(), and SymmetricBandMatrix::Maximum().

{
   REPORT
   if (storage == 0) NullMatrixError(this);
   Real maxval = -FloatingPointPrecision::Maximum(); int nr = Nrows();
   MatrixRow mr((GeneralMatrix*)this, LoadOnEntry+DirectPart);
   for (int r = 1; r <= nr; r++)
   {
      int c; maxval = mr.Maximum1(maxval, c);
      if (c > 0) { i = r; j = c; }
      mr.Next();
   }
   ((GeneralMatrix&)*this).tDelete(); return maxval;
}

MaximumAbsoluteValue()

Real GeneralMatrix::MaximumAbsoluteValue ( ) const

virtual

Reimplemented from BaseMatrix.

Reimplemented in SymmetricBandMatrix, and BandMatrix.

Definition at line 214 of file newmat8.cpp.

References REPORT.

Referenced by BandMatrix::MaximumAbsoluteValue(), SymmetricBandMatrix::MaximumAbsoluteValue(), Norm1(), and NormInfinity().

{
   REPORT
   if (storage == 0) NullMatrixError(this);
   Real maxval = 0.0; int l = storage; Real* s = store;
   while (l--) { Real a = fabs(*s++); if (maxval < a) maxval = a; }
   ((GeneralMatrix&)*this).tDelete(); return maxval;
}

MaximumAbsoluteValue1()

Real GeneralMatrix::MaximumAbsoluteValue1 ( int &  i ) const

virtual

Reimplemented from BaseMatrix.

Definition at line 223 of file newmat8.cpp.

References REPORT.

Referenced by Matrix::MaximumAbsoluteValue2().

{
   REPORT
   if (storage == 0) NullMatrixError(this);
   Real maxval = 0.0; int l = storage; Real* s = store; int li = storage;
   while (l--)
      { Real a = fabs(*s++); if (maxval <= a) { maxval = a; li = l; }  }
   i = storage - li;
   ((GeneralMatrix&)*this).tDelete(); return maxval;
}

MaximumAbsoluteValue2()

Real GeneralMatrix::MaximumAbsoluteValue2 ( int &  i, int &  j  ) const

virtual

Reimplemented from BaseMatrix.

Reimplemented in Matrix.

Definition at line 292 of file newmat8.cpp.

References REPORT.

{
   REPORT
   if (storage == 0) NullMatrixError(this);
   Real maxval = 0.0; int nr = Nrows();
   MatrixRow mr((GeneralMatrix*)this, LoadOnEntry+DirectPart);
   for (int r = 1; r <= nr; r++)
   {
      int c; maxval = mr.MaximumAbsoluteValue1(maxval, c);
      if (c > 0) { i = r; j = c; }
      mr.Next();
   }
   ((GeneralMatrix&)*this).tDelete(); return maxval;
}

MiniCleanUp()

virtual void GeneralMatrix::MiniCleanUp ( )

inlineprotectedvirtual

Reimplemented in BandLUMatrix, CroutMatrix, ColumnVector, RowVector, and nricMatrix.

Definition at line 421 of file newmat.h.

Referenced by nricMatrix::MiniCleanUp(), CroutMatrix::MiniCleanUp(), BandLUMatrix::MiniCleanUp(), and tDelete().

{ store = 0; storage = 0; nrows = 0; ncols = 0; tag = -1;}

Minimum()

Real GeneralMatrix::Minimum ( ) const

virtual

Reimplemented from BaseMatrix.

Reimplemented in SymmetricBandMatrix, and BandMatrix.

Definition at line 273 of file newmat8.cpp.

References REPORT.

{
   REPORT
   if (storage == 0) NullMatrixError(this);
   int l = storage - 1; Real* s = store; Real minval = *s++;
   while (l--) { Real a = *s++; if (minval > a) minval = a; }
   ((GeneralMatrix&)*this).tDelete(); return minval;
}

Minimum1()

Real GeneralMatrix::Minimum1 ( int &  i ) const

virtual

Reimplemented from BaseMatrix.

Definition at line 282 of file newmat8.cpp.

References REPORT.

Referenced by Matrix::Minimum2().

{
   REPORT
   if (storage == 0) NullMatrixError(this);
   int l = storage - 1; Real* s = store; Real minval = *s++; int li = l;
   while (l--) { Real a = *s++; if (minval >= a) { minval = a; li = l; } }
   i = storage - li;
   ((GeneralMatrix&)*this).tDelete(); return minval;
}

Minimum2()

Real GeneralMatrix::Minimum2 ( int &  i, int &  j  ) const

virtual

Reimplemented from BaseMatrix.

Reimplemented in Matrix.

Definition at line 337 of file newmat8.cpp.

References REPORT.

Referenced by BandMatrix::Minimum(), and SymmetricBandMatrix::Minimum().

{
   REPORT
   if (storage == 0) NullMatrixError(this);
   Real minval = FloatingPointPrecision::Maximum(); int nr = Nrows();
   MatrixRow mr((GeneralMatrix*)this, LoadOnEntry+DirectPart);
   for (int r = 1; r <= nr; r++)
   {
      int c; minval = mr.Minimum1(minval, c);
      if (c > 0) { i = r; j = c; }
      mr.Next();
   }
   ((GeneralMatrix&)*this).tDelete(); return minval;
}

MinimumAbsoluteValue()

Real GeneralMatrix::MinimumAbsoluteValue ( ) const

virtual

Reimplemented from BaseMatrix.

Reimplemented in SymmetricBandMatrix, and BandMatrix.

Definition at line 234 of file newmat8.cpp.

References REPORT.

{
   REPORT
   if (storage == 0) NullMatrixError(this);
   int l = storage - 1; Real* s = store; Real minval = fabs(*s++);
   while (l--) { Real a = fabs(*s++); if (minval > a) minval = a; }
   ((GeneralMatrix&)*this).tDelete(); return minval;
}

MinimumAbsoluteValue1()

Real GeneralMatrix::MinimumAbsoluteValue1 ( int &  i ) const

virtual

Reimplemented from BaseMatrix.

Definition at line 243 of file newmat8.cpp.

References REPORT.

Referenced by Matrix::MinimumAbsoluteValue2().

{
   REPORT
   if (storage == 0) NullMatrixError(this);
   int l = storage - 1; Real* s = store; Real minval = fabs(*s++); int li = l;
   while (l--)
      { Real a = fabs(*s++); if (minval >= a) { minval = a; li = l; }  }
   i = storage - li;
   ((GeneralMatrix&)*this).tDelete(); return minval;
}

MinimumAbsoluteValue2()

Real GeneralMatrix::MinimumAbsoluteValue2 ( int &  i, int &  j  ) const

virtual

Reimplemented from BaseMatrix.

Reimplemented in Matrix.

Definition at line 307 of file newmat8.cpp.

References REPORT.

Referenced by BandMatrix::MinimumAbsoluteValue(), and SymmetricBandMatrix::MinimumAbsoluteValue().

{
   REPORT
   if (storage == 0)  NullMatrixError(this);
   Real minval = FloatingPointPrecision::Maximum(); int nr = Nrows();
   MatrixRow mr((GeneralMatrix*)this, LoadOnEntry+DirectPart);
   for (int r = 1; r <= nr; r++)
   {
      int c; minval = mr.MinimumAbsoluteValue1(minval, c);
      if (c > 0) { i = r; j = c; }
      mr.Next();
   }
   ((GeneralMatrix&)*this).tDelete(); return minval;
}

MinusEqual() [1/2]

void GeneralMatrix::MinusEqual ( const GeneralMatrix &  gm )

protected

Definition at line 178 of file newmat7.cpp.

References ncols, nrows, and REPORT.

{
   REPORT
   if (nrows != gm.nrows || ncols != gm.ncols)
      Throw(IncompatibleDimensionsException(*this, gm));
   SubtractFrom(this, &gm);
}

MinusEqual() [2/2]

void GeneralMatrix::MinusEqual ( Real  f )

protected

Multiply() [1/2]

void GeneralMatrix::Multiply ( GeneralMatrix *  gm1, Real  f  )

protected

Definition at line 368 of file newmat5.cpp.

References REPORT, and store.

{
   REPORT
   Real* s1=gm1->store; Real* s=store;  int i=(storage >> 2);
   while (i--)
   { *s++ = *s1++ * f; *s++ = *s1++ * f; *s++ = *s1++ * f; *s++ = *s1++ * f; }
   i = storage & 3; while (i--) *s++ = *s1++ * f;
}

Multiply() [2/2]

void GeneralMatrix::Multiply ( Real  f )

protected

Definition at line 377 of file newmat5.cpp.

References REPORT.

{
   REPORT
   Real* s=store; int i=(storage >> 2);
   while (i--) { *s++ *= f; *s++ *= f; *s++ *= f; *s++ *= f; }
   i = storage & 3; while (i--) *s++ *= f;
}

Ncols()

int GeneralMatrix::Ncols ( ) const

inline

Definition at line 431 of file newmat.h.

Referenced by AddedMatrix::Evaluate(), SPMatrix::Evaluate(), ConcatenatedMatrix::Evaluate(), StackedMatrix::Evaluate(), SubtractedMatrix::Evaluate(), FFT2(), FFT2I(), GetMatrix(), GetSubMatrix::Inject(), BaseMatrix::Norm1(), GetSubMatrix::operator+=(), GetSubMatrix::operator-=(), GetSubMatrix::operator<<(), GetSubMatrix::operator=(), QRZ(), QRZT(), RectMatrixRow::Reset(), RectMatrixCol::Reset(), GetSubMatrix::SetUpLHS(), UpdateQRZ(), and UpdateQRZT().

{ return ncols; }

NegAdd() [1/2]

void GeneralMatrix::NegAdd ( GeneralMatrix *  gm1, Real  f  )

protected

Definition at line 328 of file newmat5.cpp.

References REPORT, and store.

{
   REPORT
   Real* s1=gm1->store; Real* s=store; int i=(storage >> 2);
   while (i--)
   { *s++ = f - *s1++; *s++ = f - *s1++; *s++ = f - *s1++; *s++ = f - *s1++; }
   i = storage & 3; while (i--) *s++ = f - *s1++;
}

NegAdd() [2/2]

void GeneralMatrix::NegAdd ( Real  f )

protected

Definition at line 337 of file newmat5.cpp.

References REPORT.

{
   REPORT
   Real* s=store; int i=(storage >> 2);
   while (i--)
   {
      *s = f - *s; s++; *s = f - *s; s++;
      *s = f - *s; s++; *s = f - *s; s++;
   }
   i = storage & 3; while (i--)  { *s = f - *s; s++; }
}

Negate() [1/2]

void GeneralMatrix::Negate ( GeneralMatrix *  gm1 )

protected

Definition at line 349 of file newmat5.cpp.

References REPORT, and store.

{
   // change sign of elements
   REPORT
   Real* s1=gm1->store; Real* s=store; int i=(storage >> 2);
   while (i--)
   { *s++ = -(*s1++); *s++ = -(*s1++); *s++ = -(*s1++); *s++ = -(*s1++); }
   i = storage & 3; while(i--) *s++ = -(*s1++);
}

Negate() [2/2]

void GeneralMatrix::Negate ( )

protected

Definition at line 359 of file newmat5.cpp.

References REPORT.

{
   REPORT
   Real* s=store; int i=(storage >> 2);
   while (i--)
   { *s = -(*s); s++; *s = -(*s); s++; *s = -(*s); s++; *s = -(*s); s++; }
   i = storage & 3; while(i--) { *s = -(*s); s++; }
}

NextCol() [1/2]

void GeneralMatrix::NextCol ( MatrixRowCol &  mrc )

virtual

Reimplemented in IdentityMatrix, RowVector, DiagonalMatrix, and Matrix.

Definition at line 90 of file newmat3.cpp.

References MatrixRowCol::cw, REPORT, MatrixRowCol::rowcol, and StoreOnExit.

{
   REPORT                                                // 423
   if (+(mrc.cw*StoreOnExit)) { REPORT this->RestoreCol(mrc); }
   mrc.rowcol++;
   if (mrc.rowcol<ncols) { REPORT this->GetCol(mrc); }
   else { REPORT mrc.cw -= StoreOnExit; }
}

NextCol() [2/2]

void GeneralMatrix::NextCol ( MatrixColX &  mrc )

virtual

Reimplemented in IdentityMatrix, RowVector, DiagonalMatrix, and Matrix.

Definition at line 99 of file newmat3.cpp.

References MatrixRowCol::cw, REPORT, MatrixRowCol::rowcol, and StoreOnExit.

{
   REPORT                                                // 423
   if (+(mrc.cw*StoreOnExit)) { REPORT this->RestoreCol(mrc); }
   mrc.rowcol++;
   if (mrc.rowcol<ncols) { REPORT this->GetCol(mrc); }
   else { REPORT mrc.cw -= StoreOnExit; }
}

NextRow()

void GeneralMatrix::NextRow ( MatrixRowCol &  mrc )

virtual

Reimplemented in IdentityMatrix, BandMatrix, DiagonalMatrix, LowerTriangularMatrix, UpperTriangularMatrix, and Matrix.

Definition at line 81 of file newmat3.cpp.

References MatrixRowCol::cw, REPORT, MatrixRowCol::rowcol, and StoreOnExit.

{
   REPORT
   if (+(mrc.cw*StoreOnExit)) { REPORT this->RestoreRow(mrc); }
   mrc.rowcol++;
   if (mrc.rowcol<nrows) { REPORT this->GetRow(mrc); }
   else { REPORT mrc.cw -= StoreOnExit; }
}

Nrows()

int GeneralMatrix::Nrows ( ) const

inline

Definition at line 430 of file newmat.h.

Referenced by KPMatrix::BandWidth(), Cholesky(), DCT(), DCT_II(), DCT_II_inverse(), DCT_inverse(), DowndateCholesky(), DST(), DST_II(), DST_II_inverse(), DST_inverse(), AddedMatrix::Evaluate(), SPMatrix::Evaluate(), ConcatenatedMatrix::Evaluate(), StackedMatrix::Evaluate(), SubtractedMatrix::Evaluate(), ShiftedMatrix::Evaluate(), NegShiftedMatrix::Evaluate(), ScaledMatrix::Evaluate(), NegatedMatrix::Evaluate(), FFT(), FFT2(), FFT2I(), FFTI(), FindMaximum2::Fit(), NonLinearLeastSquares::Fit(), GetMatrix(), GetSubMatrix::Inject(), LeftCircularUpdateCholesky(), BaseMatrix::NormInfinity(), GetSubMatrix::operator+=(), GetSubMatrix::operator-=(), operator<<(), GetSubMatrix::operator<<(), GetSubMatrix::operator=(), QRZ(), QRZT(), RealFFT(), RectMatrixCol::Reset(), RightCircularUpdateCholesky(), GetSubMatrix::SetUpLHS(), UpdateCholesky(), UpdateQRZ(), and UpdateQRZT().

{ return nrows; }          // get dimensions

operator &=()

void GeneralMatrix::operator&=

(

const BaseMatrix &

)

operator ReturnMatrix()

GeneralMatrix::operator ReturnMatrix

(

)

const

Definition at line 197 of file newmat4.cpp.

References REPORT.

{
   REPORT
   GeneralMatrix* gm = Image(); gm->ReleaseAndDelete();
   return ReturnMatrix(gm);
}

operator*=() [1/2]

void GeneralMatrix::operator*=

(

const BaseMatrix &

X

)

Definition at line 495 of file newmat6.cpp.

References Evaluate(), and REPORT.

{
   REPORT
   Tracer tr("GeneralMatrix::operator*=");
   // MatrixConversionCheck mcc;
   Protect();                                   // so it cannot get deleted
                        // during Evaluate
   GeneralMatrix* gm = ((BaseMatrix&)X).Evaluate();
   MultipliedMatrix am(this,gm);
   if (gm==this) Release(2); else Release();
   Eq2(am,Type());
}

operator*=() [2/2]

void GeneralMatrix::operator*=

(

Real

r

)

Definition at line 543 of file newmat6.cpp.

References REPORT.

{
   REPORT
   Tracer tr("GeneralMatrix::operator*=(Real)");
   // MatrixConversionCheck mcc;
   ScaledMatrix am(this,r);
   Release(); Eq2(am,Type());
}

operator+=() [1/2]

void GeneralMatrix::operator+=

(

const BaseMatrix &

X

)

Definition at line 469 of file newmat6.cpp.

References Evaluate(), and REPORT.

Referenced by UpperTriangularMatrix::operator+=(), LowerTriangularMatrix::operator+=(), and DiagonalMatrix::operator+=().

{
   REPORT
   Tracer tr("GeneralMatrix::operator+=");
   // MatrixConversionCheck mcc;
   Protect();                                   // so it cannot get deleted
                        // during Evaluate
   GeneralMatrix* gm = ((BaseMatrix&)X).Evaluate();
   AddedMatrix am(this,gm);
   if (gm==this) Release(2); else Release();
   Eq2(am,Type());
}

operator+=() [2/2]

void GeneralMatrix::operator+=

(

Real

r

)

Definition at line 534 of file newmat6.cpp.

References REPORT.

{
   REPORT
   Tracer tr("GeneralMatrix::operator+=(Real)");
   // MatrixConversionCheck mcc;
   ShiftedMatrix am(this,r);
   Release(); Eq2(am,Type());
}

operator-=() [1/2]

void GeneralMatrix::operator-=

(

const BaseMatrix &

X

)

Definition at line 482 of file newmat6.cpp.

References Evaluate(), and REPORT.

Referenced by UpperTriangularMatrix::operator-=(), LowerTriangularMatrix::operator-=(), and DiagonalMatrix::operator-=().

{
   REPORT
   Tracer tr("GeneralMatrix::operator-=");
   // MatrixConversionCheck mcc;
   Protect();                                   // so it cannot get deleted
                        // during Evaluate
   GeneralMatrix* gm = ((BaseMatrix&)X).Evaluate();
   SubtractedMatrix am(this,gm);
   if (gm==this) Release(2); else Release();
   Eq2(am,Type());
}

operator-=() [2/2]

void GeneralMatrix::operator-= ( Real  r )

inline

Definition at line 454 of file newmat.h.

{ operator+=(-r); }

operator/=()

void GeneralMatrix::operator/= ( Real  r )

inline

Definition at line 456 of file newmat.h.

{ operator*=(1.0/r); }

operator<<() [1/5]

void GeneralMatrix::operator<<

(

const Real *

r

)

Definition at line 422 of file newmat6.cpp.

References REPORT.

Referenced by BandMatrix::operator<<(), and SymmetricBandMatrix::operator<<().

{
   REPORT
   int i = storage; Real* s=store;
   while(i--) *s++ = *r++;
}

operator<<() [2/5]

void GeneralMatrix::operator<<

(

const int *

r

)

Definition at line 430 of file newmat6.cpp.

References REPORT.

{
   REPORT
   int i = storage; Real* s=store;
   while(i--) *s++ = *r++;
}

operator<<() [3/5]

void GeneralMatrix::operator<< ( const BaseMatrix &  X )

inline

Definition at line 445 of file newmat.h.

      { Eq(X,this->Type(),true); }              // = without checking type

operator<<() [4/5]

MatrixInput GeneralMatrix::operator<<

(

Real

f

)

Definition at line 401 of file newmat5.cpp.

References n, and REPORT.

{
   REPORT
   Tracer et("MatrixInput");
   int n = Storage();
   if (n<=0) Throw(ProgramException("Loading data to zero length matrix"));
   Real* r; r = Store(); *r = f; n--;
   return MatrixInput(n, r+1);
}

operator<<() [5/5]

MatrixInput GeneralMatrix::operator<< ( int  f )

inline

Definition at line 1777 of file newmat.h.

{ return *this << (Real)f; }

operator=()

void GeneralMatrix::operator= ( Real  f )

protected

Definition at line 342 of file newmat6.cpp.

References REPORT.

Referenced by Matrix::operator=(), SquareMatrix::operator=(), nricMatrix::operator=(), SymmetricMatrix::operator=(), UpperTriangularMatrix::operator=(), LowerTriangularMatrix::operator=(), DiagonalMatrix::operator=(), RowVector::operator=(), ColumnVector::operator=(), BandMatrix::operator=(), UpperBandMatrix::operator=(), LowerBandMatrix::operator=(), SymmetricBandMatrix::operator=(), and IdentityMatrix::operator=().

{ REPORT int i=storage; Real* s=store; while (i--) { *s++ = f; } }

operator|=()

void GeneralMatrix::operator|=

(

const BaseMatrix &

X

)

Definition at line 508 of file newmat6.cpp.

References Evaluate(), and REPORT.

{
   REPORT
   Tracer tr("GeneralMatrix::operator|=");
   // MatrixConversionCheck mcc;
   Protect();                                   // so it cannot get deleted
                        // during Evaluate
   GeneralMatrix* gm = ((BaseMatrix&)X).Evaluate();
   ConcatenatedMatrix am(this,gm);
   if (gm==this) Release(2); else Release();
   Eq2(am,Type());
}

PlusEqual() [1/2]

void GeneralMatrix::PlusEqual ( const GeneralMatrix &  gm )

protected

Definition at line 148 of file newmat7.cpp.

References ncols, nrows, and REPORT.

{
   REPORT
   if (nrows != gm.nrows || ncols != gm.ncols)
      Throw(IncompatibleDimensionsException(*this, gm));
   AddTo(this, &gm);
}

PlusEqual() [2/2]

void GeneralMatrix::PlusEqual ( Real  f )

protected

Protect()

void GeneralMatrix::Protect ( )

inline

Definition at line 437 of file newmat.h.

Referenced by Eq(), Eq2(), LinearEquationSolver::LinearEquationSolver(), MatrixType::New(), GenericMatrix::operator*=(), GenericMatrix::operator+=(), GenericMatrix::operator-=(), and GenericMatrix::operator|=().

{ tag=-1; }                   // cannot delete or reuse

Release() [1/2]

void GeneralMatrix::Release ( )

inline

Definition at line 440 of file newmat.h.

Referenced by Cholesky(), Eq(), GenericMatrix::operator*=(), GenericMatrix::operator+=(), GenericMatrix::operator-=(), GenericMatrix::operator|=(), and GenericMatrix::Release().

{ tag=1; }                    // del store after next use

Release() [2/2]

void GeneralMatrix::Release ( int  t )

inline

Definition at line 441 of file newmat.h.

{ tag=t; }               // del store after t accesses

ReleaseAndDelete()

void GeneralMatrix::ReleaseAndDelete ( )

inline

Definition at line 442 of file newmat.h.

Referenced by SPMatrix::Evaluate(), ConcatenatedMatrix::Evaluate(), StackedMatrix::Evaluate(), ShiftedMatrix::Evaluate(), NegShiftedMatrix::Evaluate(), ScaledMatrix::Evaluate(), NegatedMatrix::Evaluate(), and Transpose().

{ tag=0; }           // delete matrix after use

ReSize() [1/2]

void GeneralMatrix::ReSize ( int  nr, int  nc, int  s  )

protected

Definition at line 216 of file newmat4.cpp.

References MatrixErrorNoSpace(), MONITOR_REAL_DELETE, MONITOR_REAL_NEW, ncols, nrows, REPORT, storage, store, and tag.

Referenced by Matrix::ReSize(), SquareMatrix::ReSize(), SymmetricMatrix::ReSize(), UpperTriangularMatrix::ReSize(), LowerTriangularMatrix::ReSize(), DiagonalMatrix::ReSize(), RowVector::ReSize(), ColumnVector::ReSize(), BandMatrix::ReSize(), SymmetricBandMatrix::ReSize(), IdentityMatrix::ReSize(), ReSizeForAdd(), and ReSizeForSP().

{
   REPORT
   if (store)
   {
      MONITOR_REAL_DELETE("Free (ReDimensi)",storage,store)
      delete [] store;
   }
   storage=s; nrows=nr; ncols=nc; tag=-1;
   if (s)
   {
      store = new Real [storage]; MatrixErrorNoSpace(store);
      MONITOR_REAL_NEW("Make (ReDimensi)",storage,store)
   }
   else store = 0;
}

ReSize() [2/2]

void GeneralMatrix::ReSize ( const GeneralMatrix &  A )

virtual

Reimplemented in IdentityMatrix, SymmetricBandMatrix, LowerBandMatrix, UpperBandMatrix, BandMatrix, ColumnVector, RowVector, DiagonalMatrix, LowerTriangularMatrix, UpperTriangularMatrix, SymmetricMatrix, nricMatrix, SquareMatrix, and Matrix.

Definition at line 367 of file newmat4.cpp.

References REPORT.

{
   REPORT
   Tracer tr("GeneralMatrix::ReSize(GM)");
   Throw(NotDefinedException("ReSize", "this type of matrix"));
}

ReSizeForAdd()

void GeneralMatrix::ReSizeForAdd ( const GeneralMatrix &  A, const GeneralMatrix &  B  )

virtual

Reimplemented in SymmetricBandMatrix, and BandMatrix.

Definition at line 374 of file newmat4.cpp.

References A, REPORT, and ReSize().

{ REPORT ReSize(A); }

ReSizeForSP()

void GeneralMatrix::ReSizeForSP ( const GeneralMatrix &  A, const GeneralMatrix &  B  )

virtual

Reimplemented in SymmetricBandMatrix, and BandMatrix.

Definition at line 377 of file newmat4.cpp.

References A, REPORT, and ReSize().

{ REPORT ReSize(A); }

RestoreCol() [1/2]

virtual void GeneralMatrix::RestoreCol ( MatrixRowCol &  )

inlinevirtual

Reimplemented in SymmetricBandMatrix, BandMatrix, RowVector, LowerTriangularMatrix, UpperTriangularMatrix, SymmetricMatrix, and Matrix.

Definition at line 464 of file newmat.h.

{}    // Restore matrix col

RestoreCol() [2/2]

virtual void GeneralMatrix::RestoreCol ( MatrixColX &  )

inlinevirtual

Reimplemented in SymmetricBandMatrix, BandMatrix, RowVector, LowerTriangularMatrix, UpperTriangularMatrix, SymmetricMatrix, and Matrix.

Definition at line 465 of file newmat.h.

{}      // Restore matrix col

RestoreRow()

virtual void GeneralMatrix::RestoreRow ( MatrixRowCol &  )

inlinevirtual

Definition at line 460 of file newmat.h.

{}    // Restore matrix row

reuse()

bool GeneralMatrix::reuse

(

)

Definition at line 568 of file newmat4.cpp.

References MatrixErrorNoSpace(), MONITOR_REAL_NEW, REPORT, storage, and tag.

Referenced by AddedMatrix::Evaluate(), SPMatrix::Evaluate(), and SubtractedMatrix::Evaluate().

{
   if (tag < -1)                 // borrowed storage
   {
      if (storage)
      {
         REPORT
         Real* s = new Real [storage]; MatrixErrorNoSpace(s);
         MONITOR_REAL_NEW("Make     (reuse)",storage,s)
         BlockCopy(storage, store, s); store = s;
      }
      else { REPORT MiniCleanUp(); }                // CleanUp
      tag = 0; return true;
   }
   if (tag < 0 ) { REPORT return false; }
   if (tag <= 1 )  { REPORT return true; }
   REPORT tag--; return false;
}

ReverseElements() [1/2]

void GeneralMatrix::ReverseElements ( )

protected

Definition at line 472 of file newmat5.cpp.

References n, REPORT, and x.

Referenced by ReversedMatrix::Evaluate().

{
   // reversing in place
   REPORT
   int n = Storage(); Real* x = Store(); Real* rx = x + n;
   n /= 2;
   while (n--) { Real t = *(--rx); *rx = *x; *(x++) = t; }
}

ReverseElements() [2/2]

void GeneralMatrix::ReverseElements ( GeneralMatrix *  gm )

protected

Definition at line 464 of file newmat5.cpp.

References n, REPORT, Store(), and x.

{
   // reversing into a new matrix
   REPORT
   int n = Storage(); Real* rx = Store() + n; Real* x = gm->Store();
   while (n--) *(--rx) = *(x++);
}

SameStorageType()

bool GeneralMatrix::SameStorageType ( const GeneralMatrix &  A ) const

virtual

Reimplemented in SymmetricBandMatrix, and BandMatrix.

Definition at line 386 of file newmat4.cpp.

References A, REPORT, and Type().

{
   REPORT
   return Type() == A.Type();
}

search()

int GeneralMatrix::search ( const BaseMatrix *  s ) const

protectedvirtual

Implements BaseMatrix.

Definition at line 395 of file newmat4.cpp.

References REPORT.

Referenced by GenericMatrix::search().

{ REPORT return (s==this) ? 1 : 0; }

SetParameters()

virtual void GeneralMatrix::SetParameters ( const GeneralMatrix *  )

inlinevirtual

Reimplemented in SymmetricBandMatrix, and BandMatrix.

Definition at line 487 of file newmat.h.

Referenced by GetMatrix().

{}

SimpleAddOK()

virtual short GeneralMatrix::SimpleAddOK ( const GeneralMatrix *  )

inlineprotectedvirtual

Reimplemented in SymmetricBandMatrix, and BandMatrix.

Definition at line 419 of file newmat.h.

Referenced by AddedMatrix::Evaluate(), SPMatrix::Evaluate(), and SubtractedMatrix::Evaluate().

{ return 0; }

Solver()

virtual void GeneralMatrix::Solver ( MatrixColX &  , const MatrixColX &    )

inlinevirtual

Reimplemented in IdentityMatrix, BandLUMatrix, LowerBandMatrix, UpperBandMatrix, CroutMatrix, DiagonalMatrix, LowerTriangularMatrix, and UpperTriangularMatrix.

Definition at line 458 of file newmat.h.

{}

Storage()

int GeneralMatrix::Storage ( ) const

inline

Definition at line 432 of file newmat.h.

Referenced by SortAscending(), and SortDescending().

{ return storage; }

Store()

Real* GeneralMatrix::Store ( ) const

inline

Definition at line 433 of file newmat.h.

Referenced by Cholesky(), DCT_II(), DCT_II_inverse(), DCT_inverse(), DST_II(), DST_II_inverse(), DST_inverse(), FFT(), nricMatrix::MakeRowPointer(), QRZ(), QRZT(), RectMatrixRow::Reset(), RectMatrixCol::Reset(), ReverseElements(), SortAscending(), SortDescending(), UpdateQRZ(), and UpdateQRZT().

{ return store; }

Sum()

Real GeneralMatrix::Sum ( ) const

virtual

Reimplemented from BaseMatrix.

Reimplemented in IdentityMatrix, SymmetricBandMatrix, BandMatrix, and SymmetricMatrix.

Definition at line 168 of file newmat8.cpp.

References REPORT.

Referenced by BandMatrix::Sum().

{
   REPORT
   Real sum = 0.0; int i = storage; Real* s = store;
   while (i--) sum += *s++;
   ((GeneralMatrix&)*this).tDelete(); return sum;
}

SumAbsoluteValue()

Real GeneralMatrix::SumAbsoluteValue ( ) const

virtual

Reimplemented from BaseMatrix.

Reimplemented in IdentityMatrix, SymmetricBandMatrix, BandMatrix, and SymmetricMatrix.

Definition at line 160 of file newmat8.cpp.

References REPORT.

Referenced by BandMatrix::SumAbsoluteValue().

{
   REPORT
   Real sum = 0.0; int i = storage; Real* s = store;
   while (i--) sum += fabs(*s++);
   ((GeneralMatrix&)*this).tDelete(); return sum;
}

SumSquare()

Real GeneralMatrix::SumSquare ( ) const

virtual

Reimplemented from BaseMatrix.

Reimplemented in IdentityMatrix, SymmetricBandMatrix, BandMatrix, and SymmetricMatrix.

Definition at line 152 of file newmat8.cpp.

References REPORT, and square().

Referenced by DowndateCholesky(), and BandMatrix::SumSquare().

{
   REPORT
   Real sum = 0.0; int i = storage; Real* s = store;
   while (i--) sum += square(*s++);
   ((GeneralMatrix&)*this).tDelete(); return sum;
}

Tag()

int GeneralMatrix::Tag ( ) const

inline

Definition at line 438 of file newmat.h.

{ return tag; }

tDelete()

void GeneralMatrix::tDelete

(

)

Definition at line 535 of file newmat4.cpp.

References MiniCleanUp(), MONITOR_REAL_DELETE, REPORT, storage, store, and tag.

Referenced by BorrowStore(), AddedMatrix::Evaluate(), SPMatrix::Evaluate(), ConcatenatedMatrix::Evaluate(), StackedMatrix::Evaluate(), SubtractedMatrix::Evaluate(), ReversedMatrix::Evaluate(), GetSubMatrix::operator+=(), GetSubMatrix::operator-=(), operator<<(), GetSubMatrix::operator<<(), and GetSubMatrix::operator=().

{
   if (tag<0)
   {
      if (tag<-1) { REPORT store = 0; delete this; return; }  // borrowed
      else { REPORT return; }   // not a temporary matrix - leave alone
   }
   if (tag==1)
   {
      if (store)
      {
         REPORT  MONITOR_REAL_DELETE("Free   (tDelete)",storage,store)
         delete [] store;
      }
      MiniCleanUp(); return;                           // CleanUp
   }
   if (tag==0) { REPORT delete this; return; }

   REPORT tag--; return;
}

Transpose()

GeneralMatrix * GeneralMatrix::Transpose ( TransposedMatrix *  tm, MatrixType  mt  )

protectedvirtual

Reimplemented in IdentityMatrix, SymmetricBandMatrix, ColumnVector, RowVector, DiagonalMatrix, and SymmetricMatrix.

Definition at line 27 of file newmat5.cpp.

References Compare(), MatrixRowCol::Copy(), MatrixRowCol::Data(), DirectPart, LoadOnEntry, MatrixType::New(), MatrixRow::Next(), MatrixCol::Next(), ReleaseAndDelete(), REPORT, and StoreOnExit.

Referenced by TransposedMatrix::Evaluate().

{
   GeneralMatrix* gm1;

   if (Compare(Type().t(),mt))
   {
      REPORT
      gm1 = mt.New(ncols,nrows,tm);
      for (int i=0; i<ncols; i++)
      {
         MatrixRow mr(gm1, StoreOnExit+DirectPart, i);
         MatrixCol mc(this, mr.Data(), LoadOnEntry, i);
      }
   }
   else
   {
      REPORT
      gm1 = mt.New(ncols,nrows,tm);
      MatrixRow mr(this, LoadOnEntry);
      MatrixCol mc(gm1, StoreOnExit+DirectPart);
      int i = nrows;
      while (i--) { mc.Copy(mr); mr.Next(); mc.Next(); }
   }
   tDelete(); gm1->ReleaseAndDelete(); return gm1;
}

Type()

virtual MatrixType GeneralMatrix::Type ( ) const

pure virtual

Implemented in IdentityMatrix, BandLUMatrix, SymmetricBandMatrix, LowerBandMatrix, UpperBandMatrix, BandMatrix, CroutMatrix, ColumnVector, RowVector, DiagonalMatrix, LowerTriangularMatrix, UpperTriangularMatrix, SymmetricMatrix, SquareMatrix, and Matrix.

Referenced by SolvedMatrix::BandWidth(), InvertedMatrix::BandWidth(), BorrowStore(), AddedMatrix::Evaluate(), SPMatrix::Evaluate(), ConcatenatedMatrix::Evaluate(), StackedMatrix::Evaluate(), SubtractedMatrix::Evaluate(), ShiftedMatrix::Evaluate(), NegShiftedMatrix::Evaluate(), ScaledMatrix::Evaluate(), NegatedMatrix::Evaluate(), ReversedMatrix::Evaluate(), GetSubMatrix::operator=(), and SameStorageType().

Friends And Related Function Documentation

AddedMatrix

friend class AddedMatrix

friend

Definition at line 515 of file newmat.h.

BandMatrix

friend class BandMatrix

friend

Definition at line 510 of file newmat.h.

BaseMatrix

friend class BaseMatrix

friend

Definition at line 514 of file newmat.h.

ColedMatrix

friend class ColedMatrix

friend

Definition at line 531 of file newmat.h.

ColumnVector

friend class ColumnVector

friend

Definition at line 509 of file newmat.h.

ConcatenatedMatrix

friend class ConcatenatedMatrix

friend

Definition at line 520 of file newmat.h.

CroutMatrix

friend class CroutMatrix

friend

Definition at line 507 of file newmat.h.

DiagedMatrix

friend class DiagedMatrix

friend

Definition at line 532 of file newmat.h.

DiagonalMatrix

friend class DiagonalMatrix

friend

Definition at line 506 of file newmat.h.

GenericMatrix

friend class GenericMatrix

friend

Definition at line 537 of file newmat.h.

GetSubMatrix

friend class GetSubMatrix

friend

Definition at line 534 of file newmat.h.

InvertedMatrix

friend class InvertedMatrix

friend

Definition at line 529 of file newmat.h.

KPMatrix

friend class KPMatrix

friend

Definition at line 519 of file newmat.h.

LinearEquationSolver

friend class LinearEquationSolver

friend

Definition at line 536 of file newmat.h.

LowerBandMatrix

friend class LowerBandMatrix

friend

Definition at line 511 of file newmat.h.

LowerTriangularMatrix

friend class LowerTriangularMatrix

friend

Definition at line 505 of file newmat.h.

MatedMatrix

friend class MatedMatrix

friend

Definition at line 533 of file newmat.h.

Matrix

friend class Matrix

friend

Definition at line 500 of file newmat.h.

MultipliedMatrix

friend class MultipliedMatrix

friend

Definition at line 516 of file newmat.h.

NegatedMatrix

friend class NegatedMatrix

friend

Definition at line 528 of file newmat.h.

NegShiftedMatrix

friend class NegShiftedMatrix

friend

Definition at line 524 of file newmat.h.

nricMatrix

friend class nricMatrix

friend

Definition at line 502 of file newmat.h.

ReturnMatrix

friend class ReturnMatrix

friend

Definition at line 535 of file newmat.h.

Referenced by ForReturn().

ReversedMatrix

friend class ReversedMatrix

friend

Definition at line 527 of file newmat.h.

RowedMatrix

friend class RowedMatrix

friend

Definition at line 530 of file newmat.h.

RowVector

friend class RowVector

friend

Definition at line 508 of file newmat.h.

ScaledMatrix

friend class ScaledMatrix

friend

Definition at line 525 of file newmat.h.

ShiftedMatrix

friend class ShiftedMatrix

friend

Definition at line 523 of file newmat.h.

SolvedMatrix

friend class SolvedMatrix

friend

Definition at line 522 of file newmat.h.

SPMatrix

friend class SPMatrix

friend

Definition at line 518 of file newmat.h.

SquareMatrix

friend class SquareMatrix

friend

Definition at line 501 of file newmat.h.

StackedMatrix

friend class StackedMatrix

friend

Definition at line 521 of file newmat.h.

SubtractedMatrix

friend class SubtractedMatrix

friend

Definition at line 517 of file newmat.h.

SymmetricBandMatrix

friend class SymmetricBandMatrix

friend

Definition at line 513 of file newmat.h.

SymmetricMatrix

friend class SymmetricMatrix

friend

Definition at line 503 of file newmat.h.

TransposedMatrix

friend class TransposedMatrix

friend

Definition at line 526 of file newmat.h.

UpperBandMatrix

friend class UpperBandMatrix

friend

Definition at line 512 of file newmat.h.

UpperTriangularMatrix

friend class UpperTriangularMatrix

friend

Definition at line 504 of file newmat.h.

Member Data Documentation

ncols

int GeneralMatrix::ncols

protected

Definition at line 391 of file newmat.h.

Referenced by CleanUp(), ColumnVector::CleanUp(), ColumnVector::ColumnVector(), CroutMatrix::CroutMatrix(), DiagonalMatrix::DiagonalMatrix(), RowedMatrix::Evaluate(), ColedMatrix::Evaluate(), DiagedMatrix::Evaluate(), MatedMatrix::Evaluate(), GetMatrix(), Inject(), nricMatrix::MakeRowPointer(), Matrix::Matrix(), MinusEqual(), SquareMatrix::operator=(), PlusEqual(), ReSize(), and SquareMatrix::SquareMatrix().

nrows

int GeneralMatrix::nrows

protected

Definition at line 391 of file newmat.h.

Referenced by CleanUp(), RowVector::CleanUp(), CroutMatrix::CleanUp(), BandLUMatrix::CleanUp(), CroutMatrix::CroutMatrix(), nricMatrix::DeleteRowPointer(), DiagonalMatrix::DiagonalMatrix(), RowedMatrix::Evaluate(), ColedMatrix::Evaluate(), DiagedMatrix::Evaluate(), MatedMatrix::Evaluate(), GetMatrix(), Inject(), nricMatrix::MakeRowPointer(), Matrix::Matrix(), CroutMatrix::MiniCleanUp(), BandLUMatrix::MiniCleanUp(), MinusEqual(), SquareMatrix::operator=(), PlusEqual(), ReSize(), RowVector::RowVector(), SquareMatrix::SquareMatrix(), and CroutMatrix::~CroutMatrix().

storage

int GeneralMatrix::storage

protected

Definition at line 392 of file newmat.h.

Referenced by CleanUp(), Eq(), RowedMatrix::Evaluate(), ColedMatrix::Evaluate(), DiagedMatrix::Evaluate(), MatedMatrix::Evaluate(), GetMatrix(), GetStore(), ReSize(), reuse(), tDelete(), and ~GeneralMatrix().

store

Real* GeneralMatrix::store

protected

Definition at line 393 of file newmat.h.

Referenced by Add(), BorrowStore(), CheckStore(), CleanUp(), Eq(), GetMatrix(), Multiply(), NegAdd(), Negate(), ReSize(), IdentityMatrix::ReSize(), tDelete(), and ~GeneralMatrix().

tag

int GeneralMatrix::tag

protected

Definition at line 390 of file newmat.h.

Referenced by BorrowStore(), CleanUp(), GetMatrix(), ReSize(), reuse(), and tDelete().

---

The documentation for this class was generated from the following files:

• newmat.h
• newmat3.cpp
• newmat4.cpp
• newmat5.cpp
• newmat6.cpp
• newmat7.cpp
• newmat8.cpp