LowerTriangularMatrix Class Reference

#include <newmat.h>

Inheritance diagram for LowerTriangularMatrix:

Public Member Functions
	LowerTriangularMatrix ()
	~LowerTriangularMatrix ()
	LowerTriangularMatrix (ArrayLengthSpecifier)
	LowerTriangularMatrix (const LowerTriangularMatrix &gm)
	LowerTriangularMatrix (const BaseMatrix &M)
void	operator= (const BaseMatrix &)
void	operator= (Real f)
void	operator= (const LowerTriangularMatrix &m)
Real &	operator() (int, int)
Real &	element (int, int)
Real	operator() (int, int) const
Real	element (int, int) const
MatrixType	Type () const
GeneralMatrix *	MakeSolver ()
void	Solver (MatrixColX &, const MatrixColX &)
LogAndSign	LogDeterminant () const
Real	Trace () const
void	GetRow (MatrixRowCol &)
void	GetCol (MatrixRowCol &)
void	GetCol (MatrixColX &)
void	RestoreCol (MatrixRowCol &)
void	RestoreCol (MatrixColX &c)
void	NextRow (MatrixRowCol &)
void	ReSize (int)
void	ReSize (const GeneralMatrix &A)
MatrixBandWidth	BandWidth () const
void	operator+= (const LowerTriangularMatrix &M)
void	operator-= (const LowerTriangularMatrix &M)
void	operator+= (Real f)
void	operator-= (Real f)
Public Member Functions inherited from GeneralMatrix
GeneralMatrix *	Evaluate (MatrixType mt=MatrixTypeUnSp)
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 void	RestoreRow (MatrixRowCol &)
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
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)
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
Real	Norm1 () const
Real	NormInfinity () const
void	IEQND () const

Private Member Functions
GeneralMatrix *	Image () const

Additional Inherited Members
Protected Member Functions inherited from GeneralMatrix
	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 inherited from GeneralMatrix
int	tag
int	nrows
int	ncols
int	storage
Real *	store

Detailed Description

Definition at line 731 of file newmat.h.

Constructor & Destructor Documentation

LowerTriangularMatrix() [1/4]

LowerTriangularMatrix::LowerTriangularMatrix ( )

inline

Definition at line 735 of file newmat.h.

Referenced by Image().

:GeneralMatrix() {}

~LowerTriangularMatrix()

LowerTriangularMatrix::~LowerTriangularMatrix ( )

inline

Definition at line 736 of file newmat.h.

{}

LowerTriangularMatrix() [2/4]

LowerTriangularMatrix::LowerTriangularMatrix

(

ArrayLengthSpecifier

n

)

Definition at line 63 of file newmat4.cpp.

   : GeneralMatrix(tristore(n.Value()))
{ REPORT nrows=n.Value(); ncols=n.Value(); }

LowerTriangularMatrix() [3/4]

LowerTriangularMatrix::LowerTriangularMatrix ( const LowerTriangularMatrix &  gm )

inline

Definition at line 738 of file newmat.h.

:GeneralMatrix(gm) { GetMatrix(&gm); }

LowerTriangularMatrix() [4/4]

LowerTriangularMatrix::LowerTriangularMatrix

(

const BaseMatrix &

M

)

Definition at line 137 of file newmat4.cpp.

References GeneralMatrix::Evaluate(), GeneralMatrix::GetMatrix(), MatrixType::LT, and REPORT.

{
   REPORT // CheckConversion(M);
   // MatrixConversionCheck mcc;
   GeneralMatrix* gmx=((BaseMatrix&)M).Evaluate(MatrixType::LT);
   GetMatrix(gmx);
}

Member Function Documentation

BandWidth()

MatrixBandWidth LowerTriangularMatrix::BandWidth ( ) const

virtual

Reimplemented from BaseMatrix.

Definition at line 438 of file newmat4.cpp.

References REPORT.

   { REPORT return MatrixBandWidth(-1,0); }

element() [1/2]

Real & LowerTriangularMatrix::element	(	int	m,
		int	n
	)

Definition at line 709 of file newmat6.cpp.

References n, and REPORT.

Referenced by QRZT(), and UpdateQRZT().

{
   REPORT
   if (n<0 || m<n || m>=nrows)
      Throw(IndexException(m,n,*this,true));
   return store[tristore(m)+n];
}

element() [2/2]

Real LowerTriangularMatrix::element	(	int	m,
		int	n
	)		const

Definition at line 717 of file newmat6.cpp.

References n, and REPORT.

{
   REPORT
   if (n<0 || m<n || m>=nrows)
      Throw(IndexException(m,n,*this,true));
   return store[tristore(m)+n];
}

GetCol() [1/2]

void LowerTriangularMatrix::GetCol ( MatrixRowCol &  mrc )

virtual

Implements GeneralMatrix.

Definition at line 341 of file newmat3.cpp.

References MatrixRowCol::cw, MatrixRowCol::data, HaveStore, MatrixRowCol::length, LoadOnEntry, MatrixErrorNoSpace(), MONITOR_REAL_NEW, REPORT, MatrixRowCol::rowcol, MatrixRowCol::skip, MatrixRowCol::storage, and StoreHere.

{
   REPORT
   int col=mrc.rowcol; mrc.skip=col; mrc.length=nrows;
   int i=nrows-col; mrc.storage=i; Real* ColCopy;
   if ( +(mrc.cw*(StoreHere+HaveStore)) )
      { REPORT  ColCopy = mrc.data; }
   else
   {
      REPORT                                            // not accessed
      ColCopy = new Real [nrows]; MatrixErrorNoSpace(ColCopy);
      MONITOR_REAL_NEW("Make (LT GetCol)",nrows,ColCopy)
      mrc.cw += HaveStore; mrc.data = ColCopy;
   }

   if (+(mrc.cw*LoadOnEntry))
   {
      REPORT
      Real* Mstore = store+(col*(col+3))/2;
      // while (i--) { *ColCopy++ = *Mstore; Mstore += ++col; }
      if (i) for (;;)
         { *ColCopy++ = *Mstore; if (!(--i)) break; Mstore += ++col; }
   }
}

GetCol() [2/2]

void LowerTriangularMatrix::GetCol ( MatrixColX &  mrc )

virtual

Implements GeneralMatrix.

Definition at line 366 of file newmat3.cpp.

References MatrixRowCol::cw, MatrixRowCol::data, MatrixRowCol::length, LoadOnEntry, REPORT, MatrixRowCol::rowcol, MatrixRowCol::skip, MatrixRowCol::storage, and MatrixColX::store.

{
   REPORT
   int col=mrc.rowcol; mrc.skip=col; mrc.length=nrows;
   int i=nrows-col; mrc.storage=i; mrc.data = mrc.store + col;

   if (+(mrc.cw*LoadOnEntry))
   {
      REPORT  Real* ColCopy = mrc.data;
      Real* Mstore = store+(col*(col+3))/2;
      // while (i--) { *ColCopy++ = *Mstore; Mstore += ++col; }
      if (i) for (;;)
         { *ColCopy++ = *Mstore; if (!(--i)) break; Mstore += ++col; }
   }
}

GetRow()

void LowerTriangularMatrix::GetRow ( MatrixRowCol &  mrc )

virtual

Implements GeneralMatrix.

Definition at line 334 of file newmat3.cpp.

References MatrixRowCol::data, MatrixRowCol::length, REPORT, MatrixRowCol::rowcol, MatrixRowCol::skip, and MatrixRowCol::storage.

{
   REPORT
   int row=mrc.rowcol; mrc.skip=0; mrc.storage=row+1; mrc.length=ncols;
   mrc.data=store+(row*(row+1))/2;
}

Image()

GeneralMatrix * LowerTriangularMatrix::Image ( ) const

privatevirtual

Reimplemented from GeneralMatrix.

Definition at line 768 of file newmat4.cpp.

References LowerTriangularMatrix(), MatrixErrorNoSpace(), and REPORT.

{
   REPORT
   GeneralMatrix* gm = new LowerTriangularMatrix(*this);
   MatrixErrorNoSpace(gm); return gm;
}

LogDeterminant()

LogAndSign LowerTriangularMatrix::LogDeterminant ( ) const

virtual

Reimplemented from GeneralMatrix.

Definition at line 657 of file newmat8.cpp.

References REPORT.

{
   REPORT
   int i = nrows; LogAndSign sum; Real* s = store; int j = 2;
   // while (i--) { sum *= *s; s += j++; }
   if (i) for(;;) { sum *= *s; if (!(--i)) break; s += j++; }
   ((GeneralMatrix&)*this).tDelete(); return sum;
}

MakeSolver()

GeneralMatrix* LowerTriangularMatrix::MakeSolver ( )

inlinevirtual

Reimplemented from GeneralMatrix.

Definition at line 752 of file newmat.h.

{ return this; } // for solving

NextRow()

void LowerTriangularMatrix::NextRow ( MatrixRowCol &  mrc )

virtual

Reimplemented from GeneralMatrix.

Definition at line 392 of file newmat3.cpp.

References MatrixRowCol::IncrLT(), and REPORT.

{ REPORT mrc.IncrLT(); }

operator()() [1/2]

Real & LowerTriangularMatrix::operator()	(	int	m,
		int	n
	)

Definition at line 55 of file newmat6.cpp.

References n, and REPORT.

{
   REPORT
   if (n<=0 || m<n || m>nrows)
      Throw(IndexException(m,n,*this));
   return store[tristore(m-1)+n-1];
}

operator()() [2/2]

Real LowerTriangularMatrix::operator()	(	int	m,
		int	n
	)		const

Definition at line 167 of file newmat6.cpp.

References n, and REPORT.

{
   REPORT
   if (n<=0 || m<n || m>nrows)
      Throw(IndexException(m,n,*this));
   return store[tristore(m-1)+n-1];
}

operator+=() [1/2]

void LowerTriangularMatrix::operator+= ( const LowerTriangularMatrix &  M )

inline

Definition at line 765 of file newmat.h.

{ PlusEqual(M); }

operator+=() [2/2]

void LowerTriangularMatrix::operator+= ( Real  f )

inline

Definition at line 767 of file newmat.h.

References GeneralMatrix::operator+=().

{ GeneralMatrix::operator+=(f); }

operator-=() [1/2]

void LowerTriangularMatrix::operator-= ( const LowerTriangularMatrix &  M )

inline

Definition at line 766 of file newmat.h.

{ MinusEqual(M); }

operator-=() [2/2]

void LowerTriangularMatrix::operator-= ( Real  f )

inline

Definition at line 768 of file newmat.h.

References GeneralMatrix::operator-=().

{ GeneralMatrix::operator-=(f); }

operator=() [1/3]

void LowerTriangularMatrix::operator=

(

const BaseMatrix &

X

)

Definition at line 401 of file newmat6.cpp.

References MatrixType::LT, and REPORT.

{
   REPORT //CheckConversion(X);
   // MatrixConversionCheck mcc;
   Eq(X,MatrixType::LT);
}

operator=() [2/3]

void LowerTriangularMatrix::operator= ( Real  f )

inline

Definition at line 741 of file newmat.h.

References GeneralMatrix::operator=().

{ GeneralMatrix::operator=(f); }

operator=() [3/3]

void LowerTriangularMatrix::operator= ( const LowerTriangularMatrix &  m )

inline

Definition at line 742 of file newmat.h.

{ Eq(m); }

ReSize() [1/2]

void LowerTriangularMatrix::ReSize

(

int

nr

)

Definition at line 253 of file newmat4.cpp.

References REPORT, and GeneralMatrix::ReSize().

Referenced by QRZT(), and ReSize().

{ REPORT GeneralMatrix::ReSize(nr,nr,tristore(nr)); }

ReSize() [2/2]

void LowerTriangularMatrix::ReSize ( const GeneralMatrix &  A )

virtual

Reimplemented from GeneralMatrix.

Definition at line 343 of file newmat4.cpp.

References A, n, REPORT, and ReSize().

{
   REPORT
   int n = A.Nrows();
   if (n != A.Ncols())
   {
      Tracer tr("LowerTriangularMatrix::ReSize(GM)");
      Throw(NotSquareException(*this));
   }
   ReSize(n);
}

RestoreCol() [1/2]

void LowerTriangularMatrix::RestoreCol ( MatrixRowCol &  mrc )

virtual

Reimplemented from GeneralMatrix.

Definition at line 382 of file newmat3.cpp.

References MatrixRowCol::data, REPORT, and MatrixRowCol::rowcol.

{
   REPORT
   int col=mrc.rowcol; Real* Cstore = mrc.data;
   Real* Mstore = store+(col*(col+3))/2; int i=nrows-col;
   //while (i--) { *Mstore = *Cstore++; Mstore += ++col; }
   if (i) for (;;)
      { *Mstore = *Cstore++; if (!(--i)) break; Mstore += ++col; }
}

RestoreCol() [2/2]

void LowerTriangularMatrix::RestoreCol ( MatrixColX &  c )

inlinevirtual

Reimplemented from GeneralMatrix.

Definition at line 760 of file newmat.h.

References RestoreCol().

Referenced by RestoreCol().

{ RestoreCol((MatrixRowCol&)c); }

Solver()

void LowerTriangularMatrix::Solver ( MatrixColX &  mcout, const MatrixColX &  mcin  )

virtual

Reimplemented from GeneralMatrix.

Definition at line 70 of file newmat7.cpp.

References MatrixRowCol::data, REPORT, MatrixRowCol::skip, and MatrixRowCol::storage.

{
   REPORT
   int i = mcin.skip-mcout.skip; Real* elx = mcin.data-i;
   while (i-- > 0) *elx++ = 0.0;
   int nc = mcin.skip; i = nc+mcin.storage; elx = mcin.data+mcin.storage;
   int nr = mcout.skip+mcout.storage; int j = nr-i; i = nr-nc;
   while (j-- > 0) *elx++ = 0.0;
   Real* el = mcin.data; Real* Ael = store + (nc*(nc+1))/2; j = 0;
   while (i-- > 0)
   {
      elx = el; Real sum = 0.0; int jx = j++; Ael += nc;
      while (jx--) sum += *Ael++ * *elx++;
      *elx = (*elx - sum) / *Ael++;
   }
}

Trace()

Real LowerTriangularMatrix::Trace ( ) const

virtual

Reimplemented from BaseMatrix.

Definition at line 565 of file newmat8.cpp.

References REPORT.

{
   REPORT
   int i = nrows; Real sum = 0.0; Real* s = store; int j = 2;
   // while (i--) { sum += *s; s += j++; }
   if (i) for (;;) { sum += *s; if (!(--i)) break; s += j++; }
   ((GeneralMatrix&)*this).tDelete(); return sum;
}

Type()

MatrixType LowerTriangularMatrix::Type ( ) const

virtual

Implements GeneralMatrix.

Definition at line 417 of file newmat4.cpp.

References MatrixType::LT.

{ return MatrixType::LT; }

---

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

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