newmat.h File Reference

#include <ossim/matrix/include.h>
#include <ossim/matrix/myexcept.h>

Go to the source code of this file.

Classes
class	LogAndSign
class	MatrixType
class	MatrixBandWidth
class	ArrayLengthSpecifier
class	BaseMatrix
class	GeneralMatrix
class	Matrix
class	SquareMatrix
class	nricMatrix
class	SymmetricMatrix
class	UpperTriangularMatrix
class	LowerTriangularMatrix
class	DiagonalMatrix
class	RowVector
class	ColumnVector
class	CroutMatrix
class	BandMatrix
class	UpperBandMatrix
class	LowerBandMatrix
class	SymmetricBandMatrix
class	BandLUMatrix
class	IdentityMatrix
class	GenericMatrix
class	MultipliedMatrix
class	AddedMatrix
class	SPMatrix
class	KPMatrix
class	ConcatenatedMatrix
class	StackedMatrix
class	SolvedMatrix
class	SubtractedMatrix
class	ShiftedMatrix
class	NegShiftedMatrix
class	ScaledMatrix
class	NegatedMatrix
class	TransposedMatrix
class	ReversedMatrix
class	InvertedMatrix
class	RowedMatrix
class	ColedMatrix
class	DiagedMatrix
class	MatedMatrix
class	ReturnMatrix
class	GetSubMatrix
class	LinearEquationSolver
class	MatrixInput
class	SimpleIntArray
class	NPDException
class	ConvergenceException
class	SingularException
class	OverflowException
class	ProgramException
class	IndexException
class	VectorException
class	NotSquareException
class	SubMatrixDimensionException
class	IncompatibleDimensionsException
class	NotDefinedException
class	CannotBuildException
class	InternalException

Macros
#define	NEWMAT_LIB   0
#define	MatrixTypeUnSp   0

Functions
void	MatrixErrorNoSpace (const void *)
bool	operator== (const GeneralMatrix &A, const GeneralMatrix &B)
bool	operator== (const BaseMatrix &A, const BaseMatrix &B)
bool	operator!= (const GeneralMatrix &A, const GeneralMatrix &B)
bool	operator!= (const BaseMatrix &A, const BaseMatrix &B)
bool	operator<= (const BaseMatrix &A, const BaseMatrix &)
bool	operator>= (const BaseMatrix &A, const BaseMatrix &)
bool	operator< (const BaseMatrix &A, const BaseMatrix &)
bool	operator> (const BaseMatrix &A, const BaseMatrix &)
bool	IsZero (const BaseMatrix &A)
OSSIM_DLL Matrix	CrossProduct (const Matrix &A, const Matrix &B)
OSSIM_DLL ReturnMatrix	CrossProductRows (const Matrix &A, const Matrix &B)
OSSIM_DLL ReturnMatrix	CrossProductColumns (const Matrix &A, const Matrix &B)
OSSIM_DLL Real	DotProduct (const Matrix &A, const Matrix &B)
LogAndSign	LogDeterminant (const BaseMatrix &B)
Real	Determinant (const BaseMatrix &B)
Real	SumSquare (const BaseMatrix &B)
Real	NormFrobenius (const BaseMatrix &B)
Real	Trace (const BaseMatrix &B)
Real	SumAbsoluteValue (const BaseMatrix &B)
Real	Sum (const BaseMatrix &B)
Real	MaximumAbsoluteValue (const BaseMatrix &B)
Real	MinimumAbsoluteValue (const BaseMatrix &B)
Real	Maximum (const BaseMatrix &B)
Real	Minimum (const BaseMatrix &B)
Real	Norm1 (const BaseMatrix &B)
Real	Norm1 (RowVector &RV)
Real	NormInfinity (const BaseMatrix &B)
Real	NormInfinity (ColumnVector &CV)
bool	IsZero (const GeneralMatrix &A)

Macro Definition Documentation

MatrixTypeUnSp

#define MatrixTypeUnSp   0

Definition at line 271 of file newmat.h.

NEWMAT_LIB

#define NEWMAT_LIB   0

Definition at line 6 of file newmat.h.

Function Documentation

CrossProduct()

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

CrossProductColumns()

OSSIM_DLL 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()

OSSIM_DLL 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();
}

Determinant()

Real Determinant ( const BaseMatrix &  B )

inline

Definition at line 1753 of file newmat.h.

   { return B.Determinant(); }

DotProduct()

OSSIM_DLL Real DotProduct	(	const Matrix &	A,
		const Matrix &	B
	)

Definition at line 526 of file newmat8.cpp.

{
   REPORT
   int n = A.storage;
   if (n != B.storage) Throw(IncompatibleDimensionsException(A,B));
   Real sum = 0.0; Real* a = A.store; Real* b = B.store;
   while (n--) sum += *a++ * *b++;
   return sum;
}

IsZero() [1/2]

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

IsZero() [2/2]

bool IsZero ( const GeneralMatrix &  A )

inline

Definition at line 1773 of file newmat.h.

References A.

{ return A.IsZero(); }

LogDeterminant()

LogAndSign LogDeterminant ( const BaseMatrix &  B )

inline

Definition at line 1751 of file newmat.h.

   { return B.LogDeterminant(); }

MatrixErrorNoSpace()

void MatrixErrorNoSpace

(

const void *

)

Definition at line 292 of file newmatex.cpp.

Referenced by BandLUMatrix::BandLUMatrix(), CroutMatrix::CroutMatrix(), RowedMatrix::Evaluate(), ColedMatrix::Evaluate(), DiagedMatrix::Evaluate(), MatedMatrix::Evaluate(), GeneralMatrix::GeneralMatrix(), Matrix::GetCol(), SymmetricMatrix::GetCol(), UpperTriangularMatrix::GetCol(), LowerTriangularMatrix::GetCol(), BandMatrix::GetCol(), SymmetricBandMatrix::GetCol(), SymmetricMatrix::GetRow(), SymmetricBandMatrix::GetRow(), GeneralMatrix::GetStore(), Matrix::Image(), SquareMatrix::Image(), nricMatrix::Image(), SymmetricMatrix::Image(), UpperTriangularMatrix::Image(), LowerTriangularMatrix::Image(), DiagonalMatrix::Image(), RowVector::Image(), ColumnVector::Image(), BandMatrix::Image(), UpperBandMatrix::Image(), LowerBandMatrix::Image(), SymmetricBandMatrix::Image(), IdentityMatrix::Image(), nricMatrix::MakeRowPointer(), GeneralMatrix::MakeSolver(), Matrix::MakeSolver(), BandMatrix::MakeSolver(), SymmetricBandMatrix::MakeSolver(), MatrixType::New(), GeneralMatrix::ReSize(), GeneralMatrix::reuse(), RowVector::Transpose(), and ColumnVector::Transpose().

{ if (!v) Throw(Bad_alloc()); }

Maximum()

Real Maximum ( const BaseMatrix &  B )

inline

Definition at line 1766 of file newmat.h.

Referenced by SVD().

{ return B.Maximum(); }

MaximumAbsoluteValue()

Real MaximumAbsoluteValue ( const BaseMatrix &  B )

inline

Definition at line 1762 of file newmat.h.

   { return B.MaximumAbsoluteValue(); }

Minimum()

Real Minimum ( const BaseMatrix &  B )

inline

Definition at line 1767 of file newmat.h.

Referenced by SVD().

{ return B.Minimum(); }

MinimumAbsoluteValue()

Real MinimumAbsoluteValue ( const BaseMatrix &  B )

inline

Definition at line 1764 of file newmat.h.

   { return B.MinimumAbsoluteValue(); }

Norm1() [1/2]

Real Norm1 ( const BaseMatrix &  B )

inline

Definition at line 1768 of file newmat.h.

{ return B.Norm1(); }

Norm1() [2/2]

Real Norm1 ( RowVector &  RV )

inline

Definition at line 1769 of file newmat.h.

References GeneralMatrix::MaximumAbsoluteValue().

{ return RV.MaximumAbsoluteValue(); }

NormFrobenius()

Real NormFrobenius ( const BaseMatrix &  B )

inline

Definition at line 1756 of file newmat.h.

{ return B.NormFrobenius(); }

NormInfinity() [1/2]

Real NormInfinity ( const BaseMatrix &  B )

inline

Definition at line 1770 of file newmat.h.

{ return B.NormInfinity(); }

NormInfinity() [2/2]

Real NormInfinity ( ColumnVector &  CV )

inline

Definition at line 1771 of file newmat.h.

References GeneralMatrix::MaximumAbsoluteValue().

   { return CV.MaximumAbsoluteValue(); }

operator!=() [1/2]

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

inline

Definition at line 1723 of file newmat.h.

   { return ! (A==B); }

operator!=() [2/2]

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

inline

Definition at line 1725 of file newmat.h.

   { return ! (A==B); }

operator<()

bool operator< ( const BaseMatrix &  A, const BaseMatrix &    )

inline

Definition at line 1734 of file newmat.h.

References A.

   { A.IEQND(); return true; }

operator<=()

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

inline

Definition at line 1730 of file newmat.h.

References A.

   { A.IEQND(); return true; }

operator==() [1/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);
}

operator==() [2/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>()

bool operator> ( const BaseMatrix &  A, const BaseMatrix &    )

inline

Definition at line 1736 of file newmat.h.

References A.

   { A.IEQND(); return true; }

operator>=()

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

inline

Definition at line 1732 of file newmat.h.

References A.

   { A.IEQND(); return true; }

Sum()

Real Sum ( const BaseMatrix &  B )

inline

Definition at line 1760 of file newmat.h.

   { return B.Sum(); }

SumAbsoluteValue()

Real SumAbsoluteValue ( const BaseMatrix &  B )

inline

Definition at line 1758 of file newmat.h.

   { return B.SumAbsoluteValue(); }

SumSquare()

Real SumSquare ( const BaseMatrix &  B )

inline

Definition at line 1755 of file newmat.h.

Referenced by MLE_D_FI::NextPoint(), and BaseMatrix::NormFrobenius().

{ return B.SumSquare(); }

Trace()

Real Trace ( const BaseMatrix &  B )

inline

Definition at line 1757 of file newmat.h.

Referenced by IdentityMatrix::Sum(), and IdentityMatrix::SumAbsoluteValue().

{ return B.Trace(); }