MatrixRowCol Class Reference

#include <newmatrc.h>

Inheritance diagram for MatrixRowCol:

Public Member Functions
void	IncrMat ()
void	IncrDiag ()
void	IncrId ()
void	IncrUT ()
void	IncrLT ()
void	Zero ()
void	Add (const MatrixRowCol &)
void	AddScaled (const MatrixRowCol &, Real)
void	Add (const MatrixRowCol &, const MatrixRowCol &)
void	Add (const MatrixRowCol &, Real)
void	NegAdd (const MatrixRowCol &, Real)
void	Sub (const MatrixRowCol &)
void	Sub (const MatrixRowCol &, const MatrixRowCol &)
void	RevSub (const MatrixRowCol &)
void	ConCat (const MatrixRowCol &, const MatrixRowCol &)
void	Multiply (const MatrixRowCol &)
void	Multiply (const MatrixRowCol &, const MatrixRowCol &)
void	KP (const MatrixRowCol &, const MatrixRowCol &)
void	Copy (const MatrixRowCol &)
void	CopyCheck (const MatrixRowCol &)
void	Check (const MatrixRowCol &)
void	Check ()
void	Copy (const Real *&)
void	Copy (const int *&)
void	Copy (Real)
void	Add (Real)
void	Multiply (Real)
Real	SumAbsoluteValue ()
Real	MaximumAbsoluteValue1 (Real r, int &i)
Real	MinimumAbsoluteValue1 (Real r, int &i)
Real	Maximum1 (Real r, int &i)
Real	Minimum1 (Real r, int &i)
Real	Sum ()
void	Inject (const MatrixRowCol &)
void	Negate (const MatrixRowCol &)
void	Multiply (const MatrixRowCol &, Real)
Real *	Data ()
int	Skip ()
int	Storage ()
int	Length ()
void	Skip (int i)
void	Storage (int i)
void	Length (int i)
void	SubRowCol (MatrixRowCol &, int, int) const
	MatrixRowCol ()
	~MatrixRowCol ()

Public Attributes
int	length
int	skip
int	storage
int	rowcol
GeneralMatrix *	gm
Real *	data
LoadAndStoreFlag	cw

Friends
Real	DotProd (const MatrixRowCol &, const MatrixRowCol &)

Detailed Description

Definition at line 33 of file newmatrc.h.

Constructor & Destructor Documentation

MatrixRowCol()

MatrixRowCol::MatrixRowCol ( )

inline

Definition at line 99 of file newmatrc.h.

{}                           // to stop warning messages

~MatrixRowCol()

MatrixRowCol::~MatrixRowCol

(

)

Definition at line 824 of file newmat3.cpp.

References HaveStore, and MONITOR_REAL_DELETE.

{
   if (+(cw*HaveStore))
   {
      MONITOR_REAL_DELETE("Free    (RowCol)",-1,data)  // do not know length
      delete [] data;
   }
}

Member Function Documentation

Add() [1/4]

void MatrixRowCol::Add

(

const MatrixRowCol &

mrc

)

Definition at line 30 of file newmat2.cpp.

References data, REPORT, skip, and storage.

Referenced by ShiftedMatrix::Evaluate(), and GetSubMatrix::operator+=().

{
   // THIS += mrc
   REPORT
   int f = mrc.skip; int l = f + mrc.storage; int lx = skip + storage;
   if (f < skip) f = skip; if (l > lx) l = lx; l -= f;
   if (l<=0) return;
   Real* elx=data+(f-skip); Real* el=mrc.data+(f-mrc.skip);
   while (l--) *elx++ += *el++;
}

Add() [2/4]

void MatrixRowCol::Add	(	const MatrixRowCol &	mrc1,
		const MatrixRowCol &	mrc2
	)

Definition at line 88 of file newmat2.cpp.

References data, REPORT, skip, and storage.

{
   // THIS = mrc1 + mrc2
   int f = skip; int l = skip + storage;
   int f1 = mrc1.skip; int l1 = f1 + mrc1.storage;
   if (f1<f) f1=f; if (l1>l) l1=l;
   int f2 = mrc2.skip; int l2 = f2 + mrc2.storage;
   if (f2<f) f2=f; if (l2>l) l2=l;
   Real* el = data + (f-skip);
   Real* el1 = mrc1.data+(f1-mrc1.skip); Real* el2 = mrc2.data+(f2-mrc2.skip);
   if (f1<f2)
   {
      int i = f1-f; while (i--) *el++ = 0.0;
      if (l1<=f2)                              // disjoint
      {
         REPORT                                // not accessed
         i = l1-f1;     while (i--) *el++ = *el1++;
         i = f2-l1;     while (i--) *el++ = 0.0;
         i = l2-f2;     while (i--) *el++ = *el2++;
         i = l-l2;      while (i--) *el++ = 0.0;
      }
      else
      {
         i = f2-f1;    while (i--) *el++ = *el1++;
         if (l1<=l2)
         {
            REPORT
            i = l1-f2; while (i--) *el++ = *el1++ + *el2++;
            i = l2-l1; while (i--) *el++ = *el2++;
            i = l-l2;  while (i--) *el++ = 0.0;
         }
         else
         {
            REPORT
            i = l2-f2; while (i--) *el++ = *el1++ + *el2++;
            i = l1-l2; while (i--) *el++ = *el1++;
            i = l-l1;  while (i--) *el++ = 0.0;
         }
      }
   }
   else
   {
      int i = f2-f; while (i--) *el++ = 0.0;
      if (l2<=f1)                              // disjoint
      {
         REPORT                                // not accessed
         i = l2-f2;     while (i--) *el++ = *el2++;
         i = f1-l2;     while (i--) *el++ = 0.0;
         i = l1-f1;     while (i--) *el++ = *el1++;
         i = l-l1;      while (i--) *el++ = 0.0;
      }
      else
      {
         i = f1-f2;    while (i--) *el++ = *el2++;
         if (l2<=l1)
         {
            REPORT
            i = l2-f1; while (i--) *el++ = *el1++ + *el2++;
            i = l1-l2; while (i--) *el++ = *el1++;
            i = l-l1;  while (i--) *el++ = 0.0;
         }
         else
         {
            REPORT
            i = l1-f1; while (i--) *el++ = *el1++ + *el2++;
            i = l2-l1; while (i--) *el++ = *el2++;
            i = l-l2;  while (i--) *el++ = 0.0;
         }
      }
   }
}

Add() [3/4]

void MatrixRowCol::Add	(	const MatrixRowCol &	mrc1,
		Real	x
	)

Definition at line 233 of file newmat2.cpp.

References data, REPORT, skip, storage, and x.

{
   // THIS = mrc1 + x
   REPORT
   if (!storage) return;
   int f = mrc1.skip; int l = f + mrc1.storage; int lx = skip + storage;
   if (f < skip) { f = skip; if (l < f) l = f; }
   if (l > lx) { l = lx; if (f > lx) f = lx; }

   Real* elx = data; Real* ely = mrc1.data+(f-mrc1.skip);

   int l1 = f-skip;  while (l1--) *elx++ = x;
       l1 = l-f;     while (l1--) *elx++ = *ely++ + x;
       lx -= l;      while (lx--) *elx++ = x;
}

Add() [4/4]

void MatrixRowCol::Add

(

Real

r

)

Definition at line 564 of file newmat2.cpp.

References REPORT.

{
   // THIS += r
   REPORT
   Real* elx = data; int l = storage; while (l--) *elx++ += r;
}

AddScaled()

void MatrixRowCol::AddScaled	(	const MatrixRowCol &	mrc,
		Real	x
	)

Definition at line 41 of file newmat2.cpp.

References data, REPORT, skip, storage, and x.

{
   REPORT
   // THIS += (mrc * x)
   int f = mrc.skip; int l = f + mrc.storage; int lx = skip + storage;
   if (f < skip) f = skip; if (l > lx) l = lx; l -= f;
   if (l<=0) return;
   Real* elx=data+(f-skip); Real* el=mrc.data+(f-mrc.skip);
   while (l--) *elx++ += *el++ * x;
}

Check() [1/2]

void MatrixRowCol::Check

(

const MatrixRowCol &

mrc1

)

Definition at line 446 of file newmat2.cpp.

References REPORT, skip, and storage.

Referenced by GetSubMatrix::operator+=(), and GetSubMatrix::operator-=().

{
   REPORT
   int f = mrc1.skip; int l = f + mrc1.storage; int lx = skip + storage;
   if (f < skip || l > lx) Throw(ProgramException("Illegal Conversion"));
}

Check() [2/2]

void MatrixRowCol::Check

(

)

Definition at line 454 of file newmat2.cpp.

References REPORT.

{
   REPORT
   if (skip!=0 || storage!=length)
      Throw(ProgramException("Illegal Conversion"));
}

ConCat()

void MatrixRowCol::ConCat	(	const MatrixRowCol &	mrc1,
		const MatrixRowCol &	mrc2
	)

Definition at line 281 of file newmat2.cpp.

References data, length, REPORT, skip, and storage.

Referenced by ConcatenatedMatrix::Evaluate().

{
   // THIS = mrc1 | mrc2
   REPORT
   int f1 = mrc1.skip; int l1 = f1 + mrc1.storage; int lx = skip + storage;
   if (f1 < skip) { f1 = skip; if (l1 < f1) l1 = f1; }
   if (l1 > lx) { l1 = lx; if (f1 > lx) f1 = lx; }

   Real* elx = data;

   int i = f1-skip;  while (i--) *elx++ =0.0;
   i = l1-f1;
   if (i)                       // in case f1 would take ely out of range
      { Real* ely = mrc1.data+(f1-mrc1.skip);  while (i--) *elx++ = *ely++; }

   int f2 = mrc2.skip; int l2 = f2 + mrc2.storage; i = mrc1.length;
   int skipx = l1 - i; lx -= i; // addresses rel to second seg, maybe -ve
   if (f2 < skipx) { f2 = skipx; if (l2 < f2) l2 = f2; }
   if (l2 > lx) { l2 = lx; if (f2 > lx) f2 = lx; }

   i = f2-skipx; while (i--) *elx++ = 0.0;
   i = l2-f2;
   if (i)                       // in case f2 would take ely out of range
      { Real* ely = mrc2.data+(f2-mrc2.skip); while (i--) *elx++ = *ely++; }
   lx -= l2;     while (lx--) *elx++ = 0.0;
}

Copy() [1/4]

void MatrixRowCol::Copy

(

const MatrixRowCol &

mrc1

)

Definition at line 413 of file newmat2.cpp.

References data, REPORT, skip, and storage.

Referenced by GeneralMatrix::Evaluate(), StackedMatrix::Evaluate(), GetSubMatrix::Evaluate(), GetSubMatrix::operator<<(), GetSubMatrix::operator=(), and GeneralMatrix::Transpose().

{
   // THIS = mrc1
   REPORT
   if (!storage) return;
   int f = mrc1.skip; int l = f + mrc1.storage; int lx = skip + storage;
   if (f < skip) { f = skip; if (l < f) l = f; }
   if (l > lx) { l = lx; if (f > lx) f = lx; }

   Real* elx = data; Real* ely = 0;

   if (l-f) ely = mrc1.data+(f-mrc1.skip);

   int l1 = f-skip;  while (l1--) *elx++ = 0.0;
       l1 = l-f;     while (l1--) *elx++ = *ely++;
       lx -= l;      while (lx--) *elx++ = 0.0;
}

Copy() [2/4]

void MatrixRowCol::Copy

(

const Real *&

r

)

Definition at line 530 of file newmat2.cpp.

References REPORT.

{
   // THIS = *r
   REPORT
   Real* elx = data; const Real* ely = r+skip; r += length;
   int l = storage; while (l--) *elx++ = *ely++;
}

Copy() [3/4]

void MatrixRowCol::Copy

(

const int *&

r

)

Definition at line 538 of file newmat2.cpp.

References REPORT.

{
   // THIS = *r
   REPORT
   Real* elx = data; const int* ely = r+skip; r += length;
   int l = storage; while (l--) *elx++ = *ely++;
}

Copy() [4/4]

void MatrixRowCol::Copy

(

Real

r

)

Definition at line 546 of file newmat2.cpp.

References REPORT.

{
   // THIS = r
   REPORT  Real* elx = data; int l = storage; while (l--) *elx++ = r;
}

CopyCheck()

void MatrixRowCol::CopyCheck

(

const MatrixRowCol &

mrc1

)

Definition at line 431 of file newmat2.cpp.

References data, REPORT, skip, and storage.

Referenced by GetSubMatrix::operator=().

{
   REPORT
   if (!storage) return;
   int f = mrc1.skip; int l = f + mrc1.storage; int lx = skip + storage;
   if (f < skip || l > lx) Throw(ProgramException("Illegal Conversion"));

   Real* elx = data; Real* ely = mrc1.data+(f-mrc1.skip);

   int l1 = f-skip;  while (l1--) *elx++ = 0.0;
       l1 = l-f;     while (l1--) *elx++ = *ely++;
       lx -= l;      while (lx--) *elx++ = 0.0;
}

Data()

Real* MatrixRowCol::Data ( )

inline

Definition at line 90 of file newmatrc.h.

Referenced by GetSubMatrix::operator<<(), and GeneralMatrix::Transpose().

{ return data; }

IncrDiag()

void MatrixRowCol::IncrDiag ( )

inline

Definition at line 46 of file newmatrc.h.

Referenced by DiagonalMatrix::NextCol(), IdentityMatrix::NextCol(), and DiagonalMatrix::NextRow().

{ rowcol++; skip++; data++; }

IncrId()

void MatrixRowCol::IncrId ( )

inline

Definition at line 47 of file newmatrc.h.

Referenced by IdentityMatrix::NextCol(), and IdentityMatrix::NextRow().

{ rowcol++; skip++; }

IncrLT()

void MatrixRowCol::IncrLT ( )

inline

Definition at line 49 of file newmatrc.h.

Referenced by LowerTriangularMatrix::NextRow().

{ rowcol++; data += storage; storage++; }

IncrMat()

void MatrixRowCol::IncrMat ( )

inline

Definition at line 45 of file newmatrc.h.

Referenced by Matrix::NextRow().

{ rowcol++; data += storage; }   // used by NextRow

IncrUT()

void MatrixRowCol::IncrUT ( )

inline

Definition at line 48 of file newmatrc.h.

Referenced by UpperTriangularMatrix::NextRow().

{ rowcol++; data += storage; storage--; skip++; }

Inject()

void MatrixRowCol::Inject

(

const MatrixRowCol &

mrc

)

Definition at line 63 of file newmat2.cpp.

References data, REPORT, skip, and storage.

Referenced by GeneralMatrix::Inject(), and GetSubMatrix::Inject().

{
   REPORT
   int f = mrc.skip; int l = f + mrc.storage; int lx = skip + storage;
   if (f < skip) f = skip; if (l > lx) l = lx; l -= f;
   if (l<=0) return;
   Real* elx=data+(f-skip); Real* ely=mrc.data+(f-mrc.skip);
   while (l--) *elx++ = *ely++;
}

KP()

void MatrixRowCol::KP	(	const MatrixRowCol &	mrc1,
		const MatrixRowCol &	mrc2
	)

Definition at line 346 of file newmat2.cpp.

References data, if(), length, skip, and storage.

{
   int f = skip; int s = storage; Real* el = data; int i;

   i = mrc1.skip * mrc2.length;
   if (i > f)
   {
      i -= f; f = 0; if (i > s) { i = s; s = 0; }  else s -= i;
      while (i--) *el++ = 0.0;
      if (s == 0) return;
   }
   else f -= i;

   i = mrc1.storage; Real* el1 = mrc1.data;
   int mrc2_skip = mrc2.skip; int mrc2_storage = mrc2.storage;
   int mrc2_length = mrc2.length;
   int mrc2_remain = mrc2_length - mrc2_skip - mrc2_storage;
   while (i--)
   {
      int j; Real* el2 = mrc2.data; Real vel1 = *el1;
      if (f == 0 && mrc2_length <= s)
      {
         j = mrc2_skip; s -= j;    while (j--) *el++ = 0.0;
         j = mrc2_storage; s -= j; while (j--) *el++ = vel1 * *el2++;
         j = mrc2_remain; s -= j;  while (j--) *el++ = 0.0;
      }
      else if (f >= mrc2_length) f -= mrc2_length;
      else
      {
         j = mrc2_skip;
         if (j > f)
         {
            j -= f; f = 0; if (j > s) { j = s; s = 0; } else s -= j;
            while (j--) *el++ = 0.0;
         }
         else f -= j;

         j = mrc2_storage;
         if (j > f)
         {
            j -= f; el2 += f; f = 0; if (j > s) { j = s; s = 0; } else s -= j;
            while (j--) *el++ = vel1 * *el2++;
         }
         else f -= j;

         j = mrc2_remain;
         if (j > f)
         {
            j -= f; f = 0; if (j > s) { j = s; s = 0; } else s -= j;
            while (j--) *el++ = 0.0;
         }
         else f -= j;
      }
      if (s == 0) return;
      ++el1;
   }

   i = (mrc1.length - mrc1.skip - mrc1.storage) * mrc2.length;
   if (i > f)
   {
      i -= f; if (i > s) i = s;
      while (i--) *el++ = 0.0;
   }
}

Length() [1/2]

int MatrixRowCol::Length ( )

inline

Definition at line 93 of file newmatrc.h.

{ return length; }             // length of row or column

Length() [2/2]

void MatrixRowCol::Length ( int  i )

inline

Definition at line 96 of file newmatrc.h.

{ length=i; }

Maximum1()

Real MatrixRowCol::Maximum1	(	Real	r,
		int &	i
	)

Definition at line 602 of file newmat2.cpp.

References REPORT.

{
   REPORT
   Real* elx = data; int l = storage; int li = -1;
   while (l--) { Real f = *elx++; if (r <= f) { r = f; li = l; } }
   i = (li >= 0) ? storage - li + skip : 0;
   return r;
}

MaximumAbsoluteValue1()

Real MatrixRowCol::MaximumAbsoluteValue1	(	Real	r,
		int &	i
	)

Definition at line 582 of file newmat2.cpp.

References REPORT.

{
   REPORT
   Real* elx = data; int l = storage; int li = -1;
   while (l--) { Real f = std::fabs(*elx++); if (r <= f) { r = f; li = l; } }
   i = (li >= 0) ? storage - li + skip : 0;
   return r;
}

Minimum1()

Real MatrixRowCol::Minimum1	(	Real	r,
		int &	i
	)

Definition at line 612 of file newmat2.cpp.

References REPORT.

{
   REPORT
   Real* elx = data; int l = storage; int li = -1;
   while (l--) { Real f = *elx++; if (r >= f) { r = f; li = l; } }
   i = (li >= 0) ? storage - li + skip : 0;
   return r;
}

MinimumAbsoluteValue1()

Real MatrixRowCol::MinimumAbsoluteValue1	(	Real	r,
		int &	i
	)

Definition at line 592 of file newmat2.cpp.

References REPORT.

{
   REPORT
   Real* elx = data; int l = storage; int li = -1;
   while (l--) { Real f = std::fabs(*elx++); if (r >= f) { r = f; li = l; } }
   i = (li >= 0) ? storage - li + skip : 0;
   return r;
}

Multiply() [1/4]

void MatrixRowCol::Multiply

(

const MatrixRowCol &

mrc1

)

Definition at line 308 of file newmat2.cpp.

References data, REPORT, skip, and storage.

Referenced by ScaledMatrix::Evaluate(), and GetSubMatrix::operator*=().

{
   REPORT
   if (!storage) return;
   int f = mrc1.skip; int l = f + mrc1.storage; int lx = skip + storage;
   if (f < skip) { f = skip; if (l < f) l = f; }
   if (l > lx) { l = lx; if (f > lx) f = lx; }

   Real* elx = data; Real* ely = mrc1.data+(f-mrc1.skip);

   int l1 = f-skip;  while (l1--) *elx++ = 0;
       l1 = l-f;     while (l1--) *elx++ *= *ely++;
       lx -= l;      while (lx--) *elx++ = 0;
}

Multiply() [2/4]

void MatrixRowCol::Multiply	(	const MatrixRowCol &	mrc1,
		const MatrixRowCol &	mrc2
	)

Definition at line 324 of file newmat2.cpp.

References data, REPORT, skip, and storage.

{
   int f = skip; int l = skip + storage;
   int f1 = mrc1.skip; int l1 = f1 + mrc1.storage;
   if (f1<f) f1=f; if (l1>l) l1=l;
   int f2 = mrc2.skip; int l2 = f2 + mrc2.storage;
   if (f2<f) f2=f; if (l2>l) l2=l;
   Real* el = data + (f-skip); int i;
   if (f1<f2) f1 = f2; if (l1>l2) l1 = l2;
   if (l1<=f1) { REPORT i = l-f; while (i--) *el++ = 0.0; }  // disjoint
   else
   {
      REPORT
      Real* el1 = mrc1.data+(f1-mrc1.skip);
      Real* el2 = mrc2.data+(f1-mrc2.skip);
      i = f1-f ;    while (i--) *el++ = 0.0;
      i = l1-f1;    while (i--) *el++ = *el1++ * *el2++;
      i = l-l1;     while (i--) *el++ = 0.0;
   }
}

Multiply() [3/4]

void MatrixRowCol::Multiply

(

Real

r

)

Definition at line 558 of file newmat2.cpp.

References REPORT.

{
   // THIS *= r
   REPORT  Real* elx = data; int l = storage; while (l--) *elx++ *= r;
}

Multiply() [4/4]

void MatrixRowCol::Multiply	(	const MatrixRowCol &	mrc1,
		Real	s
	)

Definition at line 480 of file newmat2.cpp.

References data, REPORT, skip, and storage.

{
   // THIS = mrc1 * s
   REPORT
   if (!storage) return;
   int f = mrc1.skip; int l = f + mrc1.storage; int lx = skip + storage;
   if (f < skip) { f = skip; if (l < f) l = f; }
   if (l > lx) { l = lx; if (f > lx) f = lx; }

   Real* elx = data; Real* ely = mrc1.data+(f-mrc1.skip);

   int l1 = f-skip;  while (l1--) *elx++ = 0.0;
       l1 = l-f;     while (l1--) *elx++ = *ely++ * s;
       lx -= l;      while (lx--) *elx++ = 0.0;
}

NegAdd()

void MatrixRowCol::NegAdd	(	const MatrixRowCol &	mrc1,
		Real	x
	)

Definition at line 249 of file newmat2.cpp.

References data, REPORT, skip, storage, and x.

Referenced by NegShiftedMatrix::Evaluate().

{
   // THIS = x - mrc1
   REPORT
   if (!storage) return;
   int f = mrc1.skip; int l = f + mrc1.storage; int lx = skip + storage;
   if (f < skip) { f = skip; if (l < f) l = f; }
   if (l > lx) { l = lx; if (f > lx) f = lx; }

   Real* elx = data; Real* ely = mrc1.data+(f-mrc1.skip);

   int l1 = f-skip;  while (l1--) *elx++ = x;
       l1 = l-f;     while (l1--) *elx++ = x - *ely++;
       lx -= l;      while (lx--) *elx++ = x;
}

Negate()

void MatrixRowCol::Negate

(

const MatrixRowCol &

mrc1

)

Definition at line 464 of file newmat2.cpp.

References data, REPORT, skip, and storage.

Referenced by NegatedMatrix::Evaluate().

{
   // THIS = -mrc1
   REPORT
   if (!storage) return;
   int f = mrc1.skip; int l = f + mrc1.storage; int lx = skip + storage;
   if (f < skip) { f = skip; if (l < f) l = f; }
   if (l > lx) { l = lx; if (f > lx) f = lx; }

   Real* elx = data; Real* ely = mrc1.data+(f-mrc1.skip);

   int l1 = f-skip;  while (l1--) *elx++ = 0.0;
       l1 = l-f;     while (l1--) *elx++ = - *ely++;
       lx -= l;      while (lx--) *elx++ = 0.0;
}

RevSub()

void MatrixRowCol::RevSub

(

const MatrixRowCol &

mrc1

)

Definition at line 265 of file newmat2.cpp.

References data, REPORT, skip, and storage.

{
   // THIS = mrc - THIS
   REPORT
   if (!storage) return;
   int f = mrc1.skip; int l = f + mrc1.storage; int lx = skip + storage;
   if (f < skip) { f = skip; if (l < f) l = f; }
   if (l > lx) { l = lx; if (f > lx) f = lx; }

   Real* elx = data; Real* ely = mrc1.data+(f-mrc1.skip);

   int l1 = f-skip;  while (l1--) { *elx = - *elx; elx++; }
       l1 = l-f;     while (l1--) { *elx = *ely++ - *elx; elx++; }
       lx -= l;      while (lx--) { *elx = - *elx; elx++; }
}

Skip() [1/2]

int MatrixRowCol::Skip ( )

inline

Definition at line 91 of file newmatrc.h.

{ return skip; }                 // number of elements skipped

Skip() [2/2]

void MatrixRowCol::Skip ( int  i )

inline

Definition at line 94 of file newmatrc.h.

{ skip=i; }

Storage() [1/2]

int MatrixRowCol::Storage ( )

inline

Definition at line 92 of file newmatrc.h.

Referenced by GetSubMatrix::operator<<().

{ return storage; }           // number of elements stored

Storage() [2/2]

void MatrixRowCol::Storage ( int  i )

inline

Definition at line 95 of file newmatrc.h.

{ storage=i; }

Sub() [1/2]

void MatrixRowCol::Sub

(

const MatrixRowCol &

mrc

)

Definition at line 52 of file newmat2.cpp.

References data, REPORT, skip, and storage.

Referenced by GetSubMatrix::operator-=().

{
   REPORT
   // THIS -= mrc
   int f = mrc.skip; int l = f + mrc.storage; int lx = skip + storage;
   if (f < skip) f = skip; if (l > lx) l = lx; l -= f;
   if (l<=0) return;
   Real* elx=data+(f-skip); Real* el=mrc.data+(f-mrc.skip);
   while (l--) *elx++ -= *el++;
}

Sub() [2/2]

void MatrixRowCol::Sub	(	const MatrixRowCol &	mrc1,
		const MatrixRowCol &	mrc2
	)

Definition at line 160 of file newmat2.cpp.

References data, REPORT, skip, and storage.

{
   // THIS = mrc1 - mrc2
   int f = skip; int l = skip + storage;
   int f1 = mrc1.skip; int l1 = f1 + mrc1.storage;
   if (f1<f) f1=f; if (l1>l) l1=l;
   int f2 = mrc2.skip; int l2 = f2 + mrc2.storage;
   if (f2<f) f2=f; if (l2>l) l2=l;
   Real* el = data + (f-skip);
   Real* el1 = mrc1.data+(f1-mrc1.skip); Real* el2 = mrc2.data+(f2-mrc2.skip);
   if (f1<f2)
   {
      int i = f1-f; while (i--) *el++ = 0.0;
      if (l1<=f2)                              // disjoint
      {
         REPORT                                // not accessed
         i = l1-f1;     while (i--) *el++ = *el1++;
         i = f2-l1;     while (i--) *el++ = 0.0;
         i = l2-f2;     while (i--) *el++ = - *el2++;
         i = l-l2;      while (i--) *el++ = 0.0;
      }
      else
      {
         i = f2-f1;    while (i--) *el++ = *el1++;
         if (l1<=l2)
         {
            REPORT
            i = l1-f2; while (i--) *el++ = *el1++ - *el2++;
            i = l2-l1; while (i--) *el++ = - *el2++;
            i = l-l2;  while (i--) *el++ = 0.0;
         }
         else
         {
            REPORT
            i = l2-f2; while (i--) *el++ = *el1++ - *el2++;
            i = l1-l2; while (i--) *el++ = *el1++;
            i = l-l1;  while (i--) *el++ = 0.0;
         }
      }
   }
   else
   {
      int i = f2-f; while (i--) *el++ = 0.0;
      if (l2<=f1)                              // disjoint
      {
         REPORT                                // not accessed
         i = l2-f2;     while (i--) *el++ = - *el2++;
         i = f1-l2;     while (i--) *el++ = 0.0;
         i = l1-f1;     while (i--) *el++ = *el1++;
         i = l-l1;      while (i--) *el++ = 0.0;
      }
      else
      {
         i = f1-f2;    while (i--) *el++ = - *el2++;
         if (l2<=l1)
         {
            REPORT
            i = l2-f1; while (i--) *el++ = *el1++ - *el2++;
            i = l1-l2; while (i--) *el++ = *el1++;
            i = l-l1;  while (i--) *el++ = 0.0;
         }
         else
         {
            REPORT
            i = l1-f1; while (i--) *el++ = *el1++ - *el2++;
            i = l2-l1; while (i--) *el++ = - *el2++;
            i = l-l2;  while (i--) *el++ = 0.0;
         }
      }
   }
}

SubRowCol()

void MatrixRowCol::SubRowCol	(	MatrixRowCol &	mrc,
		int	skip1,
		int	l1
	)		const

Definition at line 629 of file newmat2.cpp.

References cw, data, length, skip, and storage.

Referenced by GetSubMatrix::Evaluate(), GetSubMatrix::Inject(), GetSubMatrix::operator*=(), GetSubMatrix::operator+=(), GetSubMatrix::operator-=(), GetSubMatrix::operator<<(), and GetSubMatrix::operator=().

{
   mrc.length = l1;  int d = skip - skip1;
   if (d<0) { mrc.skip = 0; mrc.data = data - d; }
   else  { mrc.skip = d; mrc.data = data; }
   d = skip + storage - skip1;
   d = ((l1 < d) ? l1 : d) - mrc.skip;  mrc.storage = (d < 0) ? 0 : d;
   mrc.cw = 0;
}

Sum()

Real MatrixRowCol::Sum

(

)

Definition at line 621 of file newmat2.cpp.

References REPORT.

{
   REPORT
   Real sum = 0.0; Real* elx = data; int l = storage;
   while (l--) sum += *elx++;
   return sum;
}

SumAbsoluteValue()

Real MatrixRowCol::SumAbsoluteValue

(

)

Definition at line 571 of file newmat2.cpp.

References REPORT.

Referenced by BaseMatrix::Norm1(), and BaseMatrix::NormInfinity().

{
   REPORT
   Real sum = 0.0; Real* elx = data; int l = storage;
   while (l--) sum += std::fabs(*elx++);
   return sum;
}

Zero()

void MatrixRowCol::Zero

(

)

Definition at line 552 of file newmat2.cpp.

References REPORT.

{
   // THIS = 0
   REPORT  Real* elx = data; int l = storage; while (l--) *elx++ = 0;
}

Friends And Related Function Documentation

DotProd

Real DotProd ( const MatrixRowCol &  mrc1, const MatrixRowCol &  mrc2  )

friend

Definition at line 74 of file newmat2.cpp.

{
   REPORT                                         // not accessed
   int f = mrc1.skip; int f2 = mrc2.skip;
   int l = f + mrc1.storage; int l2 = f2 + mrc2.storage;
   if (f < f2) f = f2; if (l > l2) l = l2; l -= f;
   if (l<=0) return 0.0;

   Real* el1=mrc1.data+(f-mrc1.skip); Real* el2=mrc2.data+(f-mrc2.skip);
   Real sum = 0.0;
   while (l--) sum += *el1++ * *el2++;
   return sum;
}

Member Data Documentation

cw

LoadAndStoreFlag MatrixRowCol::cw

Definition at line 44 of file newmatrc.h.

Referenced by Matrix::GetCol(), SymmetricMatrix::GetCol(), UpperTriangularMatrix::GetCol(), LowerTriangularMatrix::GetCol(), DiagonalMatrix::GetCol(), RowVector::GetCol(), BandMatrix::GetCol(), SymmetricBandMatrix::GetCol(), IdentityMatrix::GetCol(), SymmetricMatrix::GetRow(), SymmetricBandMatrix::GetRow(), GeneralMatrix::NextCol(), Matrix::NextCol(), DiagonalMatrix::NextCol(), RowVector::NextCol(), IdentityMatrix::NextCol(), GeneralMatrix::NextRow(), GetSubMatrix::operator<<(), Matrix::RestoreCol(), and SubRowCol().

data

Real* MatrixRowCol::data

Definition at line 43 of file newmatrc.h.

Referenced by Add(), AddScaled(), ConCat(), Copy(), CopyCheck(), DotProd(), Matrix::GetCol(), SymmetricMatrix::GetCol(), UpperTriangularMatrix::GetCol(), LowerTriangularMatrix::GetCol(), DiagonalMatrix::GetCol(), RowVector::GetCol(), BandMatrix::GetCol(), SymmetricBandMatrix::GetCol(), IdentityMatrix::GetCol(), Matrix::GetRow(), SymmetricMatrix::GetRow(), UpperTriangularMatrix::GetRow(), LowerTriangularMatrix::GetRow(), DiagonalMatrix::GetRow(), BandMatrix::GetRow(), SymmetricBandMatrix::GetRow(), IdentityMatrix::GetRow(), Inject(), KP(), Multiply(), NegAdd(), Negate(), Matrix::NextCol(), DiagonalMatrix::NextCol(), RowVector::NextCol(), IdentityMatrix::NextCol(), BandMatrix::NextRow(), operator<<(), Matrix::RestoreCol(), SymmetricMatrix::RestoreCol(), UpperTriangularMatrix::RestoreCol(), LowerTriangularMatrix::RestoreCol(), RowVector::RestoreCol(), BandMatrix::RestoreCol(), SymmetricBandMatrix::RestoreCol(), RevSub(), UpperTriangularMatrix::Solver(), LowerTriangularMatrix::Solver(), DiagonalMatrix::Solver(), CroutMatrix::Solver(), UpperBandMatrix::Solver(), LowerBandMatrix::Solver(), BandLUMatrix::Solver(), IdentityMatrix::Solver(), Sub(), and SubRowCol().

gm

GeneralMatrix* MatrixRowCol::gm

Definition at line 42 of file newmatrc.h.

length

int MatrixRowCol::length

Definition at line 38 of file newmatrc.h.

Referenced by ConCat(), Matrix::GetCol(), SymmetricMatrix::GetCol(), UpperTriangularMatrix::GetCol(), LowerTriangularMatrix::GetCol(), DiagonalMatrix::GetCol(), RowVector::GetCol(), BandMatrix::GetCol(), SymmetricBandMatrix::GetCol(), IdentityMatrix::GetCol(), Matrix::GetRow(), SymmetricMatrix::GetRow(), UpperTriangularMatrix::GetRow(), LowerTriangularMatrix::GetRow(), DiagonalMatrix::GetRow(), BandMatrix::GetRow(), SymmetricBandMatrix::GetRow(), IdentityMatrix::GetRow(), KP(), and SubRowCol().

rowcol

int MatrixRowCol::rowcol

Definition at line 41 of file newmatrc.h.

Referenced by Matrix::GetCol(), SymmetricMatrix::GetCol(), UpperTriangularMatrix::GetCol(), LowerTriangularMatrix::GetCol(), DiagonalMatrix::GetCol(), RowVector::GetCol(), BandMatrix::GetCol(), SymmetricBandMatrix::GetCol(), IdentityMatrix::GetCol(), Matrix::GetRow(), SymmetricMatrix::GetRow(), UpperTriangularMatrix::GetRow(), LowerTriangularMatrix::GetRow(), DiagonalMatrix::GetRow(), BandMatrix::GetRow(), SymmetricBandMatrix::GetRow(), IdentityMatrix::GetRow(), GeneralMatrix::NextCol(), Matrix::NextCol(), DiagonalMatrix::NextCol(), RowVector::NextCol(), IdentityMatrix::NextCol(), GeneralMatrix::NextRow(), BandMatrix::NextRow(), Matrix::RestoreCol(), SymmetricMatrix::RestoreCol(), UpperTriangularMatrix::RestoreCol(), LowerTriangularMatrix::RestoreCol(), RowVector::RestoreCol(), BandMatrix::RestoreCol(), and SymmetricBandMatrix::RestoreCol().

skip

int MatrixRowCol::skip

Definition at line 39 of file newmatrc.h.

Referenced by Add(), AddScaled(), Check(), ConCat(), Copy(), CopyCheck(), DotProd(), Matrix::GetCol(), SymmetricMatrix::GetCol(), UpperTriangularMatrix::GetCol(), LowerTriangularMatrix::GetCol(), DiagonalMatrix::GetCol(), RowVector::GetCol(), BandMatrix::GetCol(), SymmetricBandMatrix::GetCol(), IdentityMatrix::GetCol(), Matrix::GetRow(), SymmetricMatrix::GetRow(), UpperTriangularMatrix::GetRow(), LowerTriangularMatrix::GetRow(), DiagonalMatrix::GetRow(), BandMatrix::GetRow(), SymmetricBandMatrix::GetRow(), IdentityMatrix::GetRow(), Inject(), KP(), Multiply(), NegAdd(), Negate(), BandMatrix::NextRow(), operator<<(), RevSub(), UpperTriangularMatrix::Solver(), LowerTriangularMatrix::Solver(), DiagonalMatrix::Solver(), CroutMatrix::Solver(), UpperBandMatrix::Solver(), LowerBandMatrix::Solver(), BandLUMatrix::Solver(), IdentityMatrix::Solver(), Sub(), and SubRowCol().

storage

int MatrixRowCol::storage

Definition at line 40 of file newmatrc.h.

Referenced by Add(), AddScaled(), Check(), ConCat(), Copy(), CopyCheck(), DotProd(), Matrix::GetCol(), SymmetricMatrix::GetCol(), UpperTriangularMatrix::GetCol(), LowerTriangularMatrix::GetCol(), DiagonalMatrix::GetCol(), RowVector::GetCol(), BandMatrix::GetCol(), SymmetricBandMatrix::GetCol(), IdentityMatrix::GetCol(), Matrix::GetRow(), SymmetricMatrix::GetRow(), UpperTriangularMatrix::GetRow(), LowerTriangularMatrix::GetRow(), DiagonalMatrix::GetRow(), BandMatrix::GetRow(), SymmetricBandMatrix::GetRow(), IdentityMatrix::GetRow(), Inject(), KP(), Multiply(), NegAdd(), Negate(), BandMatrix::NextRow(), operator<<(), BandMatrix::RestoreCol(), SymmetricBandMatrix::RestoreCol(), RevSub(), UpperTriangularMatrix::Solver(), LowerTriangularMatrix::Solver(), DiagonalMatrix::Solver(), CroutMatrix::Solver(), UpperBandMatrix::Solver(), LowerBandMatrix::Solver(), BandLUMatrix::Solver(), IdentityMatrix::Solver(), Sub(), and SubRowCol().

---

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

• newmatrc.h
• newmat2.cpp
• newmat3.cpp