SymmetricBandMatrix Class Reference

#include <newmat.h>

Inheritance diagram for SymmetricBandMatrix:

Public Member Functions
	SymmetricBandMatrix ()
	~SymmetricBandMatrix ()
	SymmetricBandMatrix (int n, int lb)
	SymmetricBandMatrix (const BaseMatrix &)
void	operator= (const BaseMatrix &)
void	operator= (Real f)
void	operator= (const SymmetricBandMatrix &m)
Real &	operator() (int, int)
Real &	element (int, int)
Real	operator() (int, int) const
Real	element (int, int) const
MatrixType	Type () const
	SymmetricBandMatrix (const SymmetricBandMatrix &gm)
GeneralMatrix *	MakeSolver ()
Real	SumSquare () const
Real	SumAbsoluteValue () const
Real	Sum () const
Real	MaximumAbsoluteValue () const
Real	MinimumAbsoluteValue () const
Real	Maximum () const
Real	Minimum () const
Real	Trace () const
LogAndSign	LogDeterminant () const
void	GetRow (MatrixRowCol &)
void	GetCol (MatrixRowCol &)
void	GetCol (MatrixColX &)
void	RestoreCol (MatrixRowCol &)
void	RestoreCol (MatrixColX &)
GeneralMatrix *	Transpose (TransposedMatrix *, MatrixType)
void	ReSize (int, int)
void	ReSize (const GeneralMatrix &A)
bool	SameStorageType (const GeneralMatrix &A) const
void	ReSizeForAdd (const GeneralMatrix &A, const GeneralMatrix &B)
void	ReSizeForSP (const GeneralMatrix &A, const GeneralMatrix &B)
MatrixBandWidth	BandWidth () const
void	SetParameters (const GeneralMatrix *)
void	operator<< (const Real *r)
void	operator<< (const int *r)
void	operator<< (const BaseMatrix &X)
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	Solver (MatrixColX &, const MatrixColX &)
virtual void	RestoreRow (MatrixRowCol &)
virtual void	NextRow (MatrixRowCol &)
virtual void	NextCol (MatrixRowCol &)
virtual void	NextCol (MatrixColX &)
Real	MaximumAbsoluteValue1 (int &i) const
Real	MinimumAbsoluteValue1 (int &i) const
Real	Maximum1 (int &i) const
Real	Minimum1 (int &i) const
Real	MaximumAbsoluteValue2 (int &i, int &j) const
Real	MinimumAbsoluteValue2 (int &i, int &j) const
Real	Maximum2 (int &i, int &j) const
Real	Minimum2 (int &i, int &j) const
virtual bool	IsEqual (const GeneralMatrix &) const
void	CheckStore () const
	operator ReturnMatrix () const
ReturnMatrix	ForReturn () const
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

Public Attributes
int	lower

Private Member Functions
GeneralMatrix *	Image () const
void	CornerClear () const
short	SimpleAddOK (const GeneralMatrix *gm)

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
void	CheckConversion (const BaseMatrix &)
void	ReSize (int, int, int)
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 1056 of file newmat.h.

Constructor & Destructor Documentation

SymmetricBandMatrix() [1/4]

SymmetricBandMatrix::SymmetricBandMatrix ( )

inline

Definition at line 1063 of file newmat.h.

Referenced by Image().

:GeneralMatrix() { lower=0; CornerClear(); }

~SymmetricBandMatrix()

SymmetricBandMatrix::~SymmetricBandMatrix ( )

inline

Definition at line 1064 of file newmat.h.

{}

SymmetricBandMatrix() [2/4]

SymmetricBandMatrix::SymmetricBandMatrix ( int  n, int  lb  )

inline

Definition at line 1065 of file newmat.h.

References n.

{ ReSize(n,lb); CornerClear(); }

SymmetricBandMatrix() [3/4]

SymmetricBandMatrix::SymmetricBandMatrix

(

const BaseMatrix &

M

)

Definition at line 429 of file bandmat.cpp.

References GeneralMatrix::Evaluate(), REPORT, and MatrixType::SB.

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

SymmetricBandMatrix() [4/4]

SymmetricBandMatrix::SymmetricBandMatrix ( const SymmetricBandMatrix &  gm )

inline

Definition at line 1079 of file newmat.h.

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

Member Function Documentation

BandWidth()

MatrixBandWidth SymmetricBandMatrix::BandWidth ( ) const

virtual

Reimplemented from BaseMatrix.

Definition at line 530 of file bandmat.cpp.

References REPORT.

   { REPORT return MatrixBandWidth(lower,lower); }

CornerClear()

void SymmetricBandMatrix::CornerClear ( ) const

private

Definition at line 515 of file bandmat.cpp.

References REPORT.

{
   // set unused parts of BandMatrix to zero
   REPORT
   int i = lower; Real* s = store; int bw = lower + 1;
   if (i) for(;;)
   {
      int j = i;
      Real* sj = s;
      while (j--) *sj++ = 0.0;
      if (!(--i)) break;
      s += bw;
   }
}

element() [1/2]

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

Definition at line 837 of file newmat6.cpp.

References n, and REPORT.

{
   REPORT
   int w = lower+1;
   if (m>=n)
   {
      REPORT
      int i = lower+n-m;
      if ( m>=nrows || n<0 || i<0 )
         Throw(IndexException(m,n,*this,true));
      return store[w*m+i];
   }
   else
   {
      REPORT
      int i = lower+m-n;
      if ( n>=nrows || m<0 || i<0 )
         Throw(IndexException(m,n,*this,true));
      return store[w*n+i];
   }
}

element() [2/2]

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

Definition at line 859 of file newmat6.cpp.

References n, and REPORT.

{
   REPORT
   int w = lower+1;
   if (m>=n)
   {
      REPORT
      int i = lower+n-m;
      if ( m>=nrows || n<0 || i<0 )
         Throw(IndexException(m,n,*this,true));
      return store[w*m+i];
   }
   else
   {
      REPORT
      int i = lower+m-n;
      if ( n>=nrows || m<0 || i<0 )
         Throw(IndexException(m,n,*this,true));
      return store[w*n+i];
   }
}

GetCol() [1/2]

void SymmetricBandMatrix::GetCol ( MatrixRowCol &  mrc )

virtual

Implements GeneralMatrix.

Definition at line 685 of file newmat3.cpp.

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

{
   // do not allow StoreHere
   if (+(mrc.cw*StoreHere))
      Throw(InternalException("SymmetricBandMatrix::GetCol(MatrixRowCol&)"));

   int c=mrc.rowcol; int w1 = lower+1; mrc.length=nrows;
   REPORT
   int s = c-lower; int o = c*w1;
   if (s<0) { w1 += s; o -= s; s = 0; }
   mrc.skip = s;

   if (+(mrc.cw*DirectPart))
   { REPORT  mrc.data = store+o; mrc.storage = w1; }
   else
   {
      // do not allow StoreOnExit and !DirectPart
      if (+(mrc.cw*StoreOnExit))
         Throw(InternalException("SymmetricBandMatrix::GetCol(MatrixRowCol&)"));
      int w = w1+lower; s += w-ncols; Real* ColCopy;
      if (s>0) w -= s; mrc.storage = w; int w2 = w-w1;

      if ( +(mrc.cw*HaveStore) ) { REPORT ColCopy = mrc.data; }
      else
      {
         REPORT ColCopy = new Real [2*lower+1]; MatrixErrorNoSpace(ColCopy);
         MONITOR_REAL_NEW("Make (SmBGetCol)",2*lower+1,ColCopy)
         mrc.cw += HaveStore; mrc.data = ColCopy;
      }

      if (+(mrc.cw*LoadOnEntry))
      {
         REPORT
         Real* Mstore = store+o;
         while (w1--) *ColCopy++ = *Mstore++;
         Mstore--;
         while (w2--) { Mstore += lower; *ColCopy++ = *Mstore; }
      }
   }
}

GetCol() [2/2]

void SymmetricBandMatrix::GetCol ( MatrixColX &  mrc )

virtual

Implements GeneralMatrix.

Definition at line 726 of file newmat3.cpp.

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

{
   int c=mrc.rowcol; int w1 = lower+1; mrc.length=nrows;
   if (+(mrc.cw*DirectPart))
   {
      REPORT
      int b = c*w1+lower;
      mrc.skip = c; c += w1-nrows; w1 -= c; mrc.storage = w1;
      Real* ColCopy = mrc.data = mrc.store+mrc.skip;
      if (+(mrc.cw*LoadOnEntry))
      {
         REPORT
         Real* Mstore = store+b;
         // while (w1--) { *ColCopy++ = *Mstore; Mstore += lower; }
         if (w1) for (;;)
            { *ColCopy++ = *Mstore; if (!(--w1)) break; Mstore += lower; }
      }
   }
   else
   {
      REPORT
      // do not allow StoreOnExit and !DirectPart
      if (+(mrc.cw*StoreOnExit))
         Throw(InternalException("SymmetricBandMatrix::GetCol(MatrixColX&)"));
      int s = c-lower; int o = c*w1;
      if (s<0) { w1 += s; o -= s; s = 0; }
      mrc.skip = s;

      int w = w1+lower; s += w-ncols;
      if (s>0) w -= s; mrc.storage = w; int w2 = w-w1;

      Real* ColCopy = mrc.data = mrc.store+mrc.skip;

      if (+(mrc.cw*LoadOnEntry))
      {
         REPORT
         Real* Mstore = store+o;
         while (w1--) *ColCopy++ = *Mstore++;
         Mstore--;
         while (w2--) { Mstore += lower; *ColCopy++ = *Mstore; }
      }

   }
}

GetRow()

void SymmetricBandMatrix::GetRow ( MatrixRowCol &  mrc )

virtual

Implements GeneralMatrix.

Definition at line 648 of file newmat3.cpp.

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

{
   REPORT
   int r=mrc.rowcol; int s = r-lower; int w1 = lower+1; int o = r*w1;
   mrc.length = ncols;
   if (s<0) { w1 += s; o -= s; s = 0; }
   mrc.skip = s;

   if (+(mrc.cw*DirectPart))
      { REPORT  mrc.data = store+o; mrc.storage = w1; }
   else
   {
      // do not allow StoreOnExit and !DirectPart
      if (+(mrc.cw*StoreOnExit))
         Throw(InternalException("SymmetricBandMatrix::GetRow(MatrixRowCol&)"));
      int w = w1+lower; s += w-ncols; Real* RowCopy;
      if (s>0) w -= s; mrc.storage = w; int w2 = w-w1;
      if (!(mrc.cw*HaveStore))
      {
         REPORT
         RowCopy = new Real [2*lower+1]; MatrixErrorNoSpace(RowCopy);
         MONITOR_REAL_NEW("Make (SmBGetRow)",2*lower+1,RowCopy)
         mrc.cw += HaveStore; mrc.data = RowCopy;
      }
      else { REPORT  RowCopy = mrc.data; }

      if (+(mrc.cw*LoadOnEntry))
      {
          REPORT
         Real* Mstore = store+o;
         while (w1--) *RowCopy++ = *Mstore++;
         Mstore--;
         while (w2--) { Mstore += lower; *RowCopy++ = *Mstore; }
      }
   }
}

Image()

GeneralMatrix * SymmetricBandMatrix::Image ( ) const

privatevirtual

Reimplemented from GeneralMatrix.

Definition at line 817 of file newmat4.cpp.

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

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

LogDeterminant()

LogAndSign SymmetricBandMatrix::LogDeterminant ( ) const

virtual

Reimplemented from GeneralMatrix.

Definition at line 440 of file bandmat.cpp.

References REPORT.

{
   REPORT
   BandLUMatrix C(*this); return C.LogDeterminant();
}

MakeSolver()

GeneralMatrix * SymmetricBandMatrix::MakeSolver ( )

virtual

Reimplemented from GeneralMatrix.

Definition at line 422 of file bandmat.cpp.

References MatrixErrorNoSpace(), and REPORT.

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

Maximum()

Real SymmetricBandMatrix::Maximum ( ) const

inlinevirtual

Reimplemented from GeneralMatrix.

Definition at line 1088 of file newmat.h.

References GeneralMatrix::Maximum2().

{ int i, j; return GeneralMatrix::Maximum2(i, j); }

MaximumAbsoluteValue()

Real SymmetricBandMatrix::MaximumAbsoluteValue ( ) const

inlinevirtual

Reimplemented from GeneralMatrix.

Definition at line 1084 of file newmat.h.

References GeneralMatrix::MaximumAbsoluteValue().

      { CornerClear(); return GeneralMatrix::MaximumAbsoluteValue(); }

Minimum()

Real SymmetricBandMatrix::Minimum ( ) const

inlinevirtual

Reimplemented from GeneralMatrix.

Definition at line 1089 of file newmat.h.

References GeneralMatrix::Minimum2().

{ int i, j; return GeneralMatrix::Minimum2(i, j); }

MinimumAbsoluteValue()

Real SymmetricBandMatrix::MinimumAbsoluteValue ( ) const

inlinevirtual

Reimplemented from GeneralMatrix.

Definition at line 1086 of file newmat.h.

References GeneralMatrix::MinimumAbsoluteValue2().

      { int i, j; return GeneralMatrix::MinimumAbsoluteValue2(i, j); }

operator()() [1/2]

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

Definition at line 119 of file newmat6.cpp.

References n, and REPORT.

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

operator()() [2/2]

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

Definition at line 231 of file newmat6.cpp.

References n, and REPORT.

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

operator<<() [1/3]

void SymmetricBandMatrix::operator<<

(

const Real *

r

)

Definition at line 456 of file newmat5.cpp.

{ Throw(ProgramException("Cannot use array read with a BandMatrix")); }

operator<<() [2/3]

void SymmetricBandMatrix::operator<<

(

const int *

r

)

Definition at line 459 of file newmat5.cpp.

{ Throw(ProgramException("Cannot use array read with a BandMatrix")); }

operator<<() [3/3]

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

inline

Definition at line 1107 of file newmat.h.

References GeneralMatrix::operator<<().

{ GeneralMatrix::operator<<(X); }

operator=() [1/3]

void SymmetricBandMatrix::operator=

(

const BaseMatrix &

X

)

Definition at line 508 of file bandmat.cpp.

References REPORT, and MatrixType::SB.

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

operator=() [2/3]

void SymmetricBandMatrix::operator= ( Real  f )

inline

Definition at line 1068 of file newmat.h.

References GeneralMatrix::operator=().

{ GeneralMatrix::operator=(f); }

operator=() [3/3]

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

inline

Definition at line 1069 of file newmat.h.

{ Eq(m); }

ReSize() [1/2]

void SymmetricBandMatrix::ReSize	(	int	n,
		int	lb
	)

Definition at line 449 of file bandmat.cpp.

References n, REPORT, and GeneralMatrix::ReSize().

{
   REPORT
   Tracer tr("SymmetricBandMatrix::ReSize");
   if (lb<0) Throw(ProgramException("Undefined bandwidth"));
   lower = (lb<=n) ? lb : n-1;
   GeneralMatrix::ReSize(n,n,n*(lower+1));
}

ReSize() [2/2]

void SymmetricBandMatrix::ReSize ( const GeneralMatrix &  A )

virtual

Reimplemented from GeneralMatrix.

Definition at line 458 of file bandmat.cpp.

References A, MatrixBandWidth::Lower(), n, REPORT, and MatrixBandWidth::Upper().

{
   REPORT
   int n = A.Nrows();
   if (n != A.Ncols())
   {
      Tracer tr("SymmetricBandMatrix::ReSize(GM)");
      Throw(NotSquareException(*this));
   }
   MatrixBandWidth mbw = A.BandWidth(); int b = mbw.Lower();
   if (b != mbw.Upper())
   {
      Tracer tr("SymmetricBandMatrix::ReSize(GM)");
      Throw(ProgramException("Upper and lower band-widths not equal"));
   }
   ReSize(n, b);
}

ReSizeForAdd()

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

virtual

Reimplemented from GeneralMatrix.

Definition at line 483 of file bandmat.cpp.

{
   REPORT
   Tracer tr("SymmetricBandMatrix::ReSizeForAdd");
   MatrixBandWidth A_BW = A.BandWidth(); MatrixBandWidth B_BW = B.BandWidth();
   if ((A_BW.Lower() < 0) | (B_BW.Lower() < 0))
         Throw(ProgramException("Can't ReSize to SymmetricBandMatrix" ));
   // already know A and B are square
   ReSize(A.Nrows(), my_max(A_BW.Lower(), B_BW.Lower()));
}

ReSizeForSP()

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

virtual

Reimplemented from GeneralMatrix.

Definition at line 495 of file bandmat.cpp.

{
   REPORT
   Tracer tr("SymmetricBandMatrix::ReSizeForSP");
   MatrixBandWidth A_BW = A.BandWidth(); MatrixBandWidth B_BW = B.BandWidth();
   if ((A_BW.Lower() < 0) | (B_BW.Lower() < 0))
         Throw(ProgramException("Can't ReSize to SymmetricBandMatrix" ));
   // already know A and B are square
   ReSize(A.Nrows(), my_min(A_BW.Lower(), B_BW.Lower()));
}

RestoreCol() [1/2]

void SymmetricBandMatrix::RestoreCol ( MatrixRowCol &  )

inlinevirtual

Reimplemented from GeneralMatrix.

Definition at line 1095 of file newmat.h.

{}

RestoreCol() [2/2]

void SymmetricBandMatrix::RestoreCol ( MatrixColX &  mrc )

virtual

Reimplemented from GeneralMatrix.

Definition at line 771 of file newmat3.cpp.

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

{
   REPORT
   int c = mrc.rowcol;
   Real* Mstore = store + c*lower+c+lower;
   Real* Cstore = mrc.data; int w = mrc.storage;
   // while (w--) { *Mstore = *Cstore++; Mstore += lower; }
   if (w) for (;;)
      { *Mstore = *Cstore++; if (!(--w)) break; Mstore += lower; }
}

SameStorageType()

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

virtual

Reimplemented from GeneralMatrix.

Definition at line 476 of file bandmat.cpp.

References A, and REPORT.

{
   if (Type() != A.Type()) { REPORT return false; }
   REPORT
   return BandWidth() == A.BandWidth();
}

SetParameters()

void SymmetricBandMatrix::SetParameters ( const GeneralMatrix *  gmx )

virtual

Reimplemented from GeneralMatrix.

Definition at line 446 of file bandmat.cpp.

References BaseMatrix::BandWidth(), MatrixBandWidth::lower, and REPORT.

{ REPORT lower = gmx->BandWidth().lower; }

SimpleAddOK()

short SymmetricBandMatrix::SimpleAddOK ( const GeneralMatrix *  gm )

privatevirtual

Reimplemented from GeneralMatrix.

Definition at line 73 of file bandmat.cpp.

References lower, and REPORT.

{
   const SymmetricBandMatrix* bm = (const SymmetricBandMatrix*)gm;
   if (bm->lower == lower) { REPORT return 0; }
   else if (bm->lower > lower) { REPORT return 1; }
   else { REPORT return 2; }
}

Sum()

Real SymmetricBandMatrix::Sum ( ) const

virtual

Reimplemented from GeneralMatrix.

Definition at line 556 of file bandmat.cpp.

References REPORT.

{
   REPORT
   CornerClear();
   Real sum1=0.0; Real sum2=0.0; Real* s=store; int i=nrows; int l=lower;
   while (i--)
      { int j = l; while (j--) sum2 += *s++; sum1 += *s++; }
   ((GeneralMatrix&)*this).tDelete(); return sum1 + 2.0 * sum2;
}

SumAbsoluteValue()

Real SymmetricBandMatrix::SumAbsoluteValue ( ) const

virtual

Reimplemented from GeneralMatrix.

Definition at line 546 of file bandmat.cpp.

References REPORT.

{
   REPORT
   CornerClear();
   Real sum1=0.0; Real sum2=0.0; Real* s=store; int i=nrows; int l=lower;
   while (i--)
   { int j = l; while (j--) sum2 += std::fabs(*s++); sum1 += std::fabs(*s++); }
   ((GeneralMatrix&)*this).tDelete(); return sum1 + 2.0 * sum2;
}

SumSquare()

Real SymmetricBandMatrix::SumSquare ( ) const

virtual

Reimplemented from GeneralMatrix.

Definition at line 536 of file bandmat.cpp.

References REPORT, and square().

{
   REPORT
   CornerClear();
   Real sum1=0.0; Real sum2=0.0; Real* s=store; int i=nrows; int l=lower;
   while (i--)
      { int j = l; while (j--) sum2 += square(*s++); sum1 += square(*s++); }
   ((GeneralMatrix&)*this).tDelete(); return sum1 + 2.0 * sum2;
}

Trace()

Real SymmetricBandMatrix::Trace ( ) const

virtual

Reimplemented from BaseMatrix.

Definition at line 592 of file newmat8.cpp.

References REPORT.

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

Transpose()

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

virtual

Reimplemented from GeneralMatrix.

Definition at line 437 of file bandmat.cpp.

References REPORT.

{ REPORT  return Evaluate(mt); }

Type()

MatrixType SymmetricBandMatrix::Type ( ) const

virtual

Implements GeneralMatrix.

Definition at line 425 of file newmat4.cpp.

References MatrixType::SB.

{ return MatrixType::SB; }

Member Data Documentation

lower

int SymmetricBandMatrix::lower

Definition at line 1062 of file newmat.h.

Referenced by Cholesky(), and SimpleAddOK().

---

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

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