Matrix3.h

Go to the documentation of this file.

//                                                         //
// Copyright (c) 2003-2014 by The University of Queensland //
// Centre for Geoscience Computing                         //
// http://earth.uq.edu.au/centre-geoscience-computing      //
//                                                         //
// Primary Business: Brisbane, Queensland, Australia       //
// Licensed under the Open Software License version 3.0    //
// http://www.opensource.org/licenses/osl-3.0.php          //
//                                                         //

#ifndef __MATRIX3_H
#define __MATRIX3_H

#define DO_INLINE_MATRIX3 1

#if DO_INLINE_MATRIX3 >= 1
#define MATRIX3_INLINE inline
#else
#define MATRIX3_INLINE
#endif

// --- IO includes ---
#include <iostream>
using std::ostream;

class MatSingularError
{
public:
  MATRIX3_INLINE MatSingularError()
  {
  };
};


class Vec3 ;
class Matrix3
{
private:
  double m[3][3];

public:
  MATRIX3_INLINE Matrix3();
  MATRIX3_INLINE Matrix3(const Vec3&, const Vec3&,const Vec3&);
  MATRIX3_INLINE Matrix3(const double[3][3]);
  MATRIX3_INLINE Matrix3(const Matrix3&);
  MATRIX3_INLINE virtual ~Matrix3();

  MATRIX3_INLINE double det();
  Vec3 solve(const Vec3&) const;    
  Vec3 solve_homogeneous() const;
  void invert(); 
  MATRIX3_INLINE Matrix3 inv(); 
  MATRIX3_INLINE void transpose(); 
  MATRIX3_INLINE Matrix3 trans() const; 
  MATRIX3_INLINE Vec3 operator *(const Vec3&) const;
  MATRIX3_INLINE Matrix3 operator *(double) const;
  MATRIX3_INLINE Matrix3 operator /(double) const;
  MATRIX3_INLINE Matrix3 operator *(const Matrix3&) const;
  MATRIX3_INLINE Matrix3& operator +=(const Matrix3&);
  MATRIX3_INLINE Matrix3 operator +(const Matrix3&)const;
  MATRIX3_INLINE Matrix3 operator -(const Matrix3&)const;
  MATRIX3_INLINE bool operator==(const Matrix3&)const;
  MATRIX3_INLINE Matrix3 &operator=(const Matrix3&);
  MATRIX3_INLINE double trace() const;
  MATRIX3_INLINE double norm() const;

  MATRIX3_INLINE double operator()(int i, int j) const {return m[i][j];}

  MATRIX3_INLINE double& operator()(int i, int j){return m[i][j];}

  MATRIX3_INLINE friend Matrix3 operator*(double,const Matrix3&);

  MATRIX3_INLINE friend Matrix3 star(const Vec3&);

  MATRIX3_INLINE static Matrix3 Unit();

  void eigen(Vec3&,Vec3&,Vec3&,double&,double&,double&);

  // output
  MATRIX3_INLINE friend ostream& operator<<(ostream&,const Matrix3&);
}; 

#if DO_INLINE_MATRIX3 >= 1
#include "Foundation/Matrix3.hpp"
#endif

#endif //__MATRIX3_H