newmatrm.cpp File Reference

#include <cmath>
#include <ossim/matrix/newmat.h>
#include <ossim/matrix/newmatrm.h>

Go to the source code of this file.

Macros
#define	WANT_MATH
#define	REPORT   {}

Functions
void	ComplexScale (RectMatrixCol &U, RectMatrixCol &V, Real x, Real y)
void	Rotate (RectMatrixCol &U, RectMatrixCol &V, Real tau, Real s)
Real	pythag (Real f, Real g, Real &c, Real &s)

Macro Definition Documentation

REPORT

#define REPORT   {}

Definition at line 18 of file newmatrm.cpp.

Referenced by RectMatrixRowCol::AddScaled(), ComplexScale(), RectMatrixRowCol::Divide(), RectMatrixRowCol::Negate(), RectMatrixRowCol::operator*(), pythag(), RectMatrixRow::Reset(), RectMatrixCol::Reset(), Rotate(), RectMatrixRowCol::SumSquare(), and RectMatrixRowCol::Zero().

WANT_MATH

#define WANT_MATH

Definition at line 5 of file newmatrm.cpp.

Function Documentation

ComplexScale()

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

Definition at line 132 of file newmatrm.cpp.

References RectMatrixRowCol::n, n, REPORT, RectMatrixRowCol::spacing, RectMatrixRowCol::store, x, and y.

Referenced by SVD().

{
   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;
   }
}

pythag()

Real pythag	(	Real	f,
		Real	g,
		Real &	c,
		Real &	s
	)

Definition at line 186 of file newmatrm.cpp.

References REPORT, and sq.

Referenced by DowndateCholesky(), LeftCircularUpdateCholesky(), RightCircularUpdateCholesky(), SVD(), and UpdateCholesky().

{
   if (f==0 && g==0) { c=1.0; s=0.0; return 0.0; }
   Real af = f>=0 ? f : -f;
   Real ag = g>=0 ? g : -g;
   if (ag<af)
   {
      REPORT
      Real h = g/f; Real sq = std::sqrt(1.0+h*h);
      if (f<0) sq = -sq;           // make return value non-negative
      c = 1.0/sq; s = h/sq; return sq*f;
   }
   else
   {
      REPORT
      Real h = f/g; Real sq = std::sqrt(1.0+h*h);
      if (g<0) sq = -sq;
      s = 1.0/sq; c = h/sq; return sq*g;
   }
}

Rotate()

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

Definition at line 156 of file newmatrm.cpp.

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

{
   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;
   }
}