newmat7.cpp File Reference

#include <ossim/matrix/include.h>
#include <ossim/matrix/newmat.h>
#include <ossim/matrix/newmatrc.h>

Go to the source code of this file.

Macros
#define	REPORT   {}

Functions
bool	operator== (const BaseMatrix &A, const BaseMatrix &B)
bool	operator== (const GeneralMatrix &A, const GeneralMatrix &B)
bool	IsZero (const BaseMatrix &A)
void	CrossProductBody (Real *a, Real *b, Real *c)
Matrix	CrossProduct (const Matrix &A, const Matrix &B)
ReturnMatrix	CrossProductRows (const Matrix &A, const Matrix &B)
ReturnMatrix	CrossProductColumns (const Matrix &A, const Matrix &B)

Macro Definition Documentation

REPORT

#define REPORT   {}

Definition at line 18 of file newmat7.cpp.

Referenced by MultipliedMatrix::Evaluate(), AddedMatrix::Evaluate(), SPMatrix::Evaluate(), KPMatrix::Evaluate(), ConcatenatedMatrix::Evaluate(), StackedMatrix::Evaluate(), SolvedMatrix::Evaluate(), SubtractedMatrix::Evaluate(), InvertedMatrix::Evaluate(), GeneralMatrix::IsEqual(), CroutMatrix::IsEqual(), BandLUMatrix::IsEqual(), GeneralMatrix::IsZero(), IsZero(), GeneralMatrix::MakeSolver(), Matrix::MakeSolver(), GeneralMatrix::MinusEqual(), BaseMatrix::Norm1(), BaseMatrix::NormInfinity(), GeneralMatrix::PlusEqual(), UpperTriangularMatrix::Solver(), LowerTriangularMatrix::Solver(), and CroutMatrix::Solver().

Function Documentation

CrossProduct()

Matrix CrossProduct	(	const Matrix &	A,
		const Matrix &	B
	)

Definition at line 964 of file newmat7.cpp.

{
   REPORT
   int ac = A.Ncols(); int ar = A.Nrows();
   int bc = B.Ncols(); int br = B.Nrows();
   Real* a = A.Store(); Real* b = B.Store();
   if (ac == 3)
   {
      if (bc != 3 || ar != 1 || br != 1)
         { Tracer et("CrossProduct"); IncompatibleDimensionsException(A, B); }
      REPORT
      RowVector C(3);  Real* c = C.Store(); CrossProductBody(a, b, c);
      return C;
   }
   else
   {
      if (ac != 1 || bc != 1 || ar != 3 || br != 3)
         { Tracer et("CrossProduct"); IncompatibleDimensionsException(A, B); }
      REPORT
      ColumnVector C(3);  Real* c = C.Store(); CrossProductBody(a, b, c);
      return C;
   }
}

CrossProductBody()

void CrossProductBody ( Real *  a, Real *  b, Real *  c  )

inline

Definition at line 957 of file newmat7.cpp.

{
   c[0] = a[1] * b[2] - a[2] * b[1];
   c[1] = a[2] * b[0] - a[0] * b[2];
   c[2] = a[0] * b[1] - a[1] * b[0];
}

CrossProductColumns()

ReturnMatrix CrossProductColumns	(	const Matrix &	A,
		const Matrix &	B
	)

Definition at line 1009 of file newmat7.cpp.

{
   REPORT
   int n = A.Ncols();
   if (A.Nrows() != 3 || B.Nrows() != 3 || n != B.Ncols())
      { Tracer et("CrossProductColumns"); IncompatibleDimensionsException(A, B); }
   Matrix C(3, n);
   Real* a = A.Store(); Real* b = B.Store(); Real* c = C.Store();
   Real* an = a+n; Real* an2 = an+n;
   Real* bn = b+n; Real* bn2 = bn+n;
   Real* cn = c+n; Real* cn2 = cn+n;

   int i = n; 
   while (i--)
   {
      *c++   = *an    * *bn2   - *an2   * *bn;
      *cn++  = *an2++ * *b     - *a     * *bn2++;
      *cn2++ = *a++   * *bn++  - *an++  * *b++;
   }

   return C.ForReturn();
}

CrossProductRows()

ReturnMatrix CrossProductRows	(	const Matrix &	A,
		const Matrix &	B
	)

Definition at line 988 of file newmat7.cpp.

{
   REPORT
   int n = A.Nrows();
   if (A.Ncols() != 3 || B.Ncols() != 3 || n != B.Nrows())
      { Tracer et("CrossProductRows"); IncompatibleDimensionsException(A, B); }
   Matrix C(n, 3);
   Real* a = A.Store(); Real* b = B.Store(); Real* c = C.Store();
   if (n--)
   {
      for (;;)
      {
         CrossProductBody(a, b, c);
         if (!(n--)) break;
         a += 3; b += 3; c += 3;
      }
   }

   return C.ForReturn();
}

IsZero()

bool IsZero

(

const BaseMatrix &

A

)

Definition at line 889 of file newmat7.cpp.

References A, GeneralMatrix::Evaluate(), GeneralMatrix::IsZero(), and REPORT.

{
   Tracer tr("BaseMatrix::IsZero");
   REPORT
   GeneralMatrix* gm1 = 0; bool bx;
   Try { gm1=((BaseMatrix&)A).Evaluate(); bx = gm1->IsZero(); }
   CatchAll { if (gm1) gm1->tDelete(); ReThrow; }
   gm1->tDelete();
   return bx;
}

operator==() [1/2]

bool operator==	(	const BaseMatrix &	A,
		const BaseMatrix &	B
	)

Definition at line 812 of file newmat7.cpp.

{
   Tracer tr("BaseMatrix ==");
   REPORT
   GeneralMatrix* gmA = ((BaseMatrix&)A).Evaluate();
   GeneralMatrix* gmB = ((BaseMatrix&)B).Evaluate();

   if (gmA == gmB)                            // same matrix
      { REPORT gmA->tDelete(); return true; }

   if ( gmA->Nrows() != gmB->Nrows() || gmA->Ncols() != gmB->Ncols() )
                                              // different dimensions
      { REPORT gmA->tDelete(); gmB->tDelete(); return false; }

   // check for CroutMatrix or BandLUMatrix
   MatrixType AType = gmA->Type(); MatrixType BType = gmB->Type();
   if (AType.CannotConvert() || BType.CannotConvert() )
   {
      REPORT
      bool bx = gmA->IsEqual(*gmB);
      gmA->tDelete(); gmB->tDelete();
      return bx;
   }

   // is matrix storage the same
   // will need to modify if further matrix structures are introduced
   if (AType == BType && gmA->BandWidth() == gmB->BandWidth())
   {                                          // compare store
      REPORT
      bool bx = RealEqual(gmA->Store(),gmB->Store(),gmA->Storage());
      gmA->tDelete(); gmB->tDelete();
      return bx;
   }

   // matrix storage different - just subtract
   REPORT  return IsZero(*gmA-*gmB);
}

operator==() [2/2]

bool operator==	(	const GeneralMatrix &	A,
		const GeneralMatrix &	B
	)

Definition at line 850 of file newmat7.cpp.

Referenced by MatrixBandWidth::operator!=().

{
   Tracer tr("GeneralMatrix ==");
   // May or may not call tDeletes
   REPORT

   if (&A == &B)                              // same matrix
      { REPORT return true; }

   if ( A.Nrows() != B.Nrows() || A.Ncols() != B.Ncols() )
      { REPORT return false; }                // different dimensions

   // check for CroutMatrix or BandLUMatrix
   MatrixType AType = A.Type(); MatrixType BType = B.Type();
   if (AType.CannotConvert() || BType.CannotConvert() )
      { REPORT  return A.IsEqual(B); }

   // is matrix storage the same
   // will need to modify if further matrix structures are introduced
   if (AType == BType && A.BandWidth() == B.BandWidth())
      { REPORT return RealEqual(A.Store(),B.Store(),A.Storage()); }

   // matrix storage different - just subtract
   REPORT  return IsZero(A-B);
}