RectMatrixRowCol Class Reference

#include <newmatrm.h>

Inheritance diagram for RectMatrixRowCol:

Public Member Functions
Real	operator* (const RectMatrixRowCol &) const
void	AddScaled (const RectMatrixRowCol &, Real)
void	Divide (const RectMatrixRowCol &, Real)
void	Divide (Real)
void	Negate ()
void	Zero ()
Real &	operator[] (int i)
Real	SumSquare () const
Real &	First ()
void	DownDiag ()
void	UpDiag ()

Protected Member Functions
	RectMatrixRowCol (Real *st, int nx, int sp, int sh)
void	Reset (Real *st, int nx, int sp, int sh)

Protected Attributes
Real *	store
int	n
int	spacing
int	shift

Friends
void	ComplexScale (RectMatrixCol &, RectMatrixCol &, Real, Real)
void	Rotate (RectMatrixCol &, RectMatrixCol &, Real, Real)

Detailed Description

Definition at line 16 of file newmatrm.h.

Constructor & Destructor Documentation

RectMatrixRowCol()

RectMatrixRowCol::RectMatrixRowCol ( Real *  st, int  nx, int  sp, int  sh  )

inlineprotected

Definition at line 27 of file newmatrm.h.

      : store(st), n(nx), spacing(sp), shift(sh) {}

Member Function Documentation

AddScaled()

void RectMatrixRowCol::AddScaled	(	const RectMatrixRowCol &	rmrc,
		Real	r
	)

Definition at line 79 of file newmatrm.cpp.

References n, n, REPORT, spacing, and store.

Referenced by SVD().

{
   REPORT
   int i = n; Real* s = store; int d = spacing;
   Real* s1 = rmrc.store; int d1 = rmrc.spacing;
   if (i!=rmrc.n)
   {
      Tracer tr("newmatrm");
      Throw(InternalException("Dimensions differ in AddScaled"));
   }
   // while (i--) { *s += *s1 * r; s += d; s1 += d1; }
   if (i) for (;;)
      { *s += *s1 * r; if (!(--i)) break; s += d; s1 += d1; }
}

Divide() [1/2]

void RectMatrixRowCol::Divide	(	const RectMatrixRowCol &	rmrc,
		Real	r
	)

Definition at line 94 of file newmatrm.cpp.

References n, n, REPORT, spacing, and store.

Referenced by SVD().

{
   REPORT
   int i = n; Real* s = store; int d = spacing;
   Real* s1 = rmrc.store; int d1 = rmrc.spacing;
   if (i!=rmrc.n)
   {
      Tracer tr("newmatrm");
      Throw(InternalException("Dimensions differ in Divide"));
   }
   // while (i--) { *s = *s1 / r; s += d; s1 += d1; }
   if (i) for (;;) { *s = *s1 / r; if (!(--i)) break; s += d; s1 += d1; }
}

Divide() [2/2]

void RectMatrixRowCol::Divide

(

Real

r

)

Definition at line 108 of file newmatrm.cpp.

References n, and REPORT.

{
   REPORT
   int i = n; Real* s = store; int d = spacing;
   // while (i--) { *s /= r; s += d; }
   if (i) for (;;) { *s /= r; if (!(--i)) break; s += d; }
}

DownDiag()

void RectMatrixRowCol::DownDiag ( )

inline

Definition at line 41 of file newmatrm.h.

References n.

Referenced by SVD().

{ store += (shift+spacing); n--; }

First()

Real& RectMatrixRowCol::First ( )

inline

Definition at line 40 of file newmatrm.h.

Referenced by SVD().

{ return *store; }                       // get first element

Negate()

void RectMatrixRowCol::Negate

(

)

Definition at line 116 of file newmatrm.cpp.

References n, and REPORT.

Referenced by SVD().

{
   REPORT
   int i = n; Real* s = store; int d = spacing;
   // while (i--) { *s = - *s; s += d; }
   if (i) for (;;) { *s = - *s; if (!(--i)) break; s += d; }
}

operator*()

Real RectMatrixRowCol::operator*

(

const RectMatrixRowCol &

rmrc

)

const

Definition at line 62 of file newmatrm.cpp.

References n, n, REPORT, spacing, and store.

{
   REPORT
   long_Real sum = 0.0; int i = n;
   Real* s = store; int d = spacing;
   Real* s1 = rmrc.store; int d1 = rmrc.spacing;
   if (i!=rmrc.n)
   {
      Tracer tr("newmatrm");
      Throw(InternalException("Dimensions differ in *"));
   }
   // while (i--) { sum += (long_Real)*s * *s1; s += d; s1 += d1; }
   if (i) for(;;)
      { sum += (long_Real)*s * *s1; if (!(--i)) break; s += d; s1 += d1; }
   return (Real)sum;
}

operator[]()

Real& RectMatrixRowCol::operator[] ( int  i )

inline

Definition at line 38 of file newmatrm.h.

{ return *(store+i*spacing); } // element

Reset()

void RectMatrixRowCol::Reset ( Real *  st, int  nx, int  sp, int  sh  )

inlineprotected

Definition at line 29 of file newmatrm.h.

References n.

Referenced by RectMatrixRow::Reset(), and RectMatrixCol::Reset().

      { store=st; n=nx; spacing=sp; shift=sh; }

SumSquare()

Real RectMatrixRowCol::SumSquare

(

)

const

Definition at line 52 of file newmatrm.cpp.

References n, and REPORT.

Referenced by SVD().

{
   REPORT
   long_Real sum = 0.0; int i = n; Real* s = store; int d = spacing;
   // while (i--) { sum += (long_Real)*s * *s; s += d; }
   if (i) for(;;)
      { sum += (long_Real)*s * *s; if (!(--i)) break; s += d; }
   return (Real)sum;
}

UpDiag()

void RectMatrixRowCol::UpDiag ( )

inline

Definition at line 42 of file newmatrm.h.

References n.

Referenced by SVD().

{ store -= (shift+spacing); n++; }

Zero()

void RectMatrixRowCol::Zero

(

)

Definition at line 124 of file newmatrm.cpp.

References n, and REPORT.

Referenced by SVD().

{
   REPORT
   int i = n; Real* s = store; int d = spacing;
   // while (i--) { *s = 0.0; s += d; }
   if (i) for (;;) { *s = 0.0; if (!(--i)) break; s += d; }
}

Friends And Related Function Documentation

ComplexScale

void ComplexScale ( RectMatrixCol &  U, RectMatrixCol &  V, Real  x, Real  y  )

friend

Definition at line 132 of file newmatrm.cpp.

{
   REPORT
   int n = U.n;
   if (n != V.n)
   {
      Tracer tr("newmatrm");
      Throw(InternalException("Dimensions differ in ComplexScale"));
   }
   Real* u = U.store; Real* v = V.store; 
   int su = U.spacing; int sv = V.spacing;
   //while (n--)
   //{
   //   Real z = *u * x - *v * y;  *v =  *u * y + *v * x;  *u = z;
   //   u += su;  v += sv;
   //}
   if (n) for (;;)
   {
      Real z = *u * x - *v * y;  *v =  *u * y + *v * x;  *u = z;
      if (!(--n)) break;
      u += su;  v += sv;
   }
}

Rotate

void Rotate ( RectMatrixCol &  U, RectMatrixCol &  V, Real  tau, Real  s  )

friend

Definition at line 156 of file newmatrm.cpp.

{
   REPORT
   //  (U, V) = (U, V) * (c, s)  where  tau = s/(1+c), c^2 + s^2 = 1
   int n = U.n;
   if (n != V.n)
   {
      Tracer tr("newmatrm");
      Throw(InternalException("Dimensions differ in Rotate"));
   }
   Real* u = U.store; Real* v = V.store;
   int su = U.spacing; int sv = V.spacing;
   //while (n--)
   //{
   //   Real zu = *u; Real zv = *v;
   //   *u -= s * (zv + zu * tau); *v += s * (zu - zv * tau);
   //   u += su;  v += sv;
   //}
   if (n) for(;;)
   {
      Real zu = *u; Real zv = *v;
      *u -= s * (zv + zu * tau); *v += s * (zu - zv * tau);
      if (!(--n)) break;
      u += su;  v += sv;
   }
}

Member Data Documentation

n

int RectMatrixRowCol::n

protected

Definition at line 24 of file newmatrm.h.

Referenced by AddScaled(), ComplexScale(), Divide(), operator*(), and Rotate().

shift

int RectMatrixRowCol::shift

protected

Definition at line 26 of file newmatrm.h.

spacing

int RectMatrixRowCol::spacing

protected

Definition at line 25 of file newmatrm.h.

Referenced by AddScaled(), ComplexScale(), Divide(), operator*(), and Rotate().

store

Real* RectMatrixRowCol::store

protected

Definition at line 23 of file newmatrm.h.

Referenced by AddScaled(), ComplexScale(), Divide(), operator*(), and Rotate().

---

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

• newmatrm.h
• newmatrm.cpp