UpperTriangularMatrix Class Reference

#include <newmat.h>

Inheritance diagram for UpperTriangularMatrix:

Public Member Functions
	UpperTriangularMatrix ()
	~UpperTriangularMatrix ()
	UpperTriangularMatrix (ArrayLengthSpecifier)
void	operator= (const BaseMatrix &)
void	operator= (const UpperTriangularMatrix &m)
	UpperTriangularMatrix (const BaseMatrix &)
	UpperTriangularMatrix (const UpperTriangularMatrix &gm)
void	operator= (Real f)
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 UpperTriangularMatrix &M)
void	operator-= (const UpperTriangularMatrix &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 690 of file newmat.h.

Constructor & Destructor Documentation

UpperTriangularMatrix() [1/4]

UpperTriangularMatrix::UpperTriangularMatrix ( )

inline

Definition at line 694 of file newmat.h.

Referenced by Image().

:GeneralMatrix() {}

~UpperTriangularMatrix()

UpperTriangularMatrix::~UpperTriangularMatrix ( )

inline

Definition at line 695 of file newmat.h.

{}

UpperTriangularMatrix() [2/4]

UpperTriangularMatrix::UpperTriangularMatrix

(

ArrayLengthSpecifier

n

)

Definition at line 59 of file newmat4.cpp.

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

UpperTriangularMatrix() [3/4]

UpperTriangularMatrix::UpperTriangularMatrix

(

const BaseMatrix &

M

)

Definition at line 129 of file newmat4.cpp.

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

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

UpperTriangularMatrix() [4/4]

UpperTriangularMatrix::UpperTriangularMatrix ( const UpperTriangularMatrix &  gm )

inline

Definition at line 700 of file newmat.h.

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

Member Function Documentation

BandWidth()

MatrixBandWidth UpperTriangularMatrix::BandWidth ( ) const

virtual

Reimplemented from BaseMatrix.

Definition at line 435 of file newmat4.cpp.

References REPORT.

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

element() [1/2]

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

Definition at line 693 of file newmat6.cpp.

References n, and REPORT.

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

element() [2/2]

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

Definition at line 701 of file newmat6.cpp.

References n, and REPORT.

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

GetCol() [1/2]

void UpperTriangularMatrix::GetCol ( MatrixRowCol &  mrc )

virtual

Implements GeneralMatrix.

Definition at line 280 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
   mrc.skip=0; int i=mrc.rowcol+1; mrc.storage=i;
   mrc.length=nrows; Real* ColCopy;
   if ( !(mrc.cw*(StoreHere+HaveStore)) )
   {
      REPORT                                              // not accessed
      ColCopy = new Real [nrows]; MatrixErrorNoSpace(ColCopy);
      MONITOR_REAL_NEW("Make (UT GetCol)",nrows,ColCopy)
      mrc.data = ColCopy; mrc.cw += HaveStore;
   }
   else { REPORT ColCopy = mrc.data; }
   if (+(mrc.cw*LoadOnEntry))
   {
      REPORT
      Real* Mstore = store+mrc.rowcol; int j = ncols;
      // while (i--) { *ColCopy++ = *Mstore; Mstore += --j; }
      if (i) for (;;)
         { *ColCopy++ = *Mstore; if (!(--i)) break; Mstore += --j; }
   }
}

GetCol() [2/2]

void UpperTriangularMatrix::GetCol ( MatrixColX &  mrc )

virtual

Implements GeneralMatrix.

Definition at line 303 of file newmat3.cpp.

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

{
   REPORT
   mrc.skip=0; int i=mrc.rowcol+1; mrc.storage=i;
   mrc.length=nrows;
   if (+(mrc.cw*LoadOnEntry))
   {
      REPORT
      Real* ColCopy = mrc.data;
      Real* Mstore = store+mrc.rowcol; int j = ncols;
      // while (i--) { *ColCopy++ = *Mstore; Mstore += --j; }
      if (i) for (;;)
         { *ColCopy++ = *Mstore; if (!(--i)) break; Mstore += --j; }
   }
}

GetRow()

void UpperTriangularMatrix::GetRow ( MatrixRowCol &  mrc )

virtual

Implements GeneralMatrix.

Definition at line 272 of file newmat3.cpp.

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

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

Image()

GeneralMatrix * UpperTriangularMatrix::Image ( ) const

privatevirtual

Reimplemented from GeneralMatrix.

Definition at line 761 of file newmat4.cpp.

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

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

LogDeterminant()

LogAndSign UpperTriangularMatrix::LogDeterminant ( ) const

virtual

Reimplemented from GeneralMatrix.

Definition at line 666 of file newmat8.cpp.

References REPORT.

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

MakeSolver()

GeneralMatrix* UpperTriangularMatrix::MakeSolver ( )

inlinevirtual

Reimplemented from GeneralMatrix.

Definition at line 711 of file newmat.h.

{ return this; } // for solving

NextRow()

void UpperTriangularMatrix::NextRow ( MatrixRowCol &  mrc )

virtual

Reimplemented from GeneralMatrix.

Definition at line 329 of file newmat3.cpp.

References MatrixRowCol::IncrUT(), and REPORT.

{ REPORT mrc.IncrUT(); }

operator()() [1/2]

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

Definition at line 47 of file newmat6.cpp.

References n, and REPORT.

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

operator()() [2/2]

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

Definition at line 159 of file newmat6.cpp.

References n, and REPORT.

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

operator+=() [1/2]

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

inline

Definition at line 724 of file newmat.h.

{ PlusEqual(M); }

operator+=() [2/2]

void UpperTriangularMatrix::operator+= ( Real  f )

inline

Definition at line 726 of file newmat.h.

References GeneralMatrix::operator+=().

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

operator-=() [1/2]

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

inline

Definition at line 725 of file newmat.h.

{ MinusEqual(M); }

operator-=() [2/2]

void UpperTriangularMatrix::operator-= ( Real  f )

inline

Definition at line 727 of file newmat.h.

References GeneralMatrix::operator-=().

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

operator=() [1/3]

void UpperTriangularMatrix::operator=

(

const BaseMatrix &

X

)

Definition at line 394 of file newmat6.cpp.

References REPORT, and MatrixType::UT.

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

operator=() [2/3]

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

inline

Definition at line 698 of file newmat.h.

{ Eq(m); }

operator=() [3/3]

void UpperTriangularMatrix::operator= ( Real  f )

inline

Definition at line 701 of file newmat.h.

References GeneralMatrix::operator=().

{ GeneralMatrix::operator=(f); }

ReSize() [1/2]

void UpperTriangularMatrix::ReSize

(

int

nr

)

Definition at line 250 of file newmat4.cpp.

References REPORT, and GeneralMatrix::ReSize().

Referenced by QRZ(), and ReSize().

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

ReSize() [2/2]

void UpperTriangularMatrix::ReSize ( const GeneralMatrix &  A )

virtual

Reimplemented from GeneralMatrix.

Definition at line 331 of file newmat4.cpp.

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

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

RestoreCol() [1/2]

void UpperTriangularMatrix::RestoreCol ( MatrixRowCol &  mrc )

virtual

Reimplemented from GeneralMatrix.

Definition at line 319 of file newmat3.cpp.

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

{
  REPORT
  Real* Mstore = store+mrc.rowcol; int i=mrc.rowcol+1; int j = ncols;
  Real* Cstore = mrc.data;
  // while (i--) { *Mstore = *Cstore++; Mstore += --j; }
  if (i) for (;;)
     { *Mstore = *Cstore++; if (!(--i)) break; Mstore += --j; }
}

RestoreCol() [2/2]

void UpperTriangularMatrix::RestoreCol ( MatrixColX &  c )

inlinevirtual

Reimplemented from GeneralMatrix.

Definition at line 719 of file newmat.h.

References RestoreCol().

Referenced by RestoreCol().

{ RestoreCol((MatrixRowCol&)c); }

Solver()

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

virtual

Reimplemented from GeneralMatrix.

Definition at line 51 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 nr = mcin.skip+mcin.storage;
   elx = mcin.data+mcin.storage; Real* el = elx;
   int j = mcout.skip+mcout.storage-nr; int nc = ncols-nr; i = nr-mcout.skip;
   while (j-- > 0) *elx++ = 0.0;
   Real* Ael = store + (nr*(2*ncols-nr+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 = (*elx - sum) / *(--Ael);
   }
}

Trace()

Real UpperTriangularMatrix::Trace ( ) const

virtual

Reimplemented from BaseMatrix.

Definition at line 574 of file newmat8.cpp.

References REPORT.

{
   REPORT
   int i = nrows; Real sum = 0.0; Real* s = store;
   while (i) { sum += *s; s += i--; }             // won t cause a problem
   ((GeneralMatrix&)*this).tDelete(); return sum;
}

Type()

MatrixType UpperTriangularMatrix::Type ( ) const

virtual

Implements GeneralMatrix.

Definition at line 416 of file newmat4.cpp.

References MatrixType::UT.

{ return MatrixType::UT; }

---

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

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