/Users/jab/Documents/Teaching/02585/GEL2_and_demos/GEL/src/CGLA/ArithSqMat3x3Float.h

#ifndef __CGLA_ARITHSQMAT3X3FLOAT_H__
#define __CGLA_ARITHSQMAT3X3FLOAT_H__

#include "ExceptionStandard.h"
#include "CGLA.h"
#include "Vec3f.h"
#include "ArithSqMatFloat.h"

namespace CGLA 
{

  CGLA_DERIVEEXCEPTION(Mat3x3fException)
    CGLA_DERIVEEXCEPTION(Mat3x3fSingular)

    template<class V, class M>
    class ArithSqMat3x3Float: public ArithSqMatFloat<V, M, 3>
    {
    public:

      typedef V VectorType;

      typedef typename V::ScalarType ScalarType;

    public:

      ArithSqMat3x3Float(V _a, V _b, V _c): 
        ArithSqMatFloat<V, M, 3> (_a,_b,_c) {}
  
      ArithSqMat3x3Float() {}

      explicit ArithSqMat3x3Float(ScalarType a): 
        ArithSqMatFloat<V, M, 3>(a) {}

    };

  template<class V, class M>
    M invert(const ArithSqMat3x3Float<V,M>&);

  template<class V, class M>
    inline 
    typename ArithSqMat3x3Float<V,M>::ScalarType
    determinant(const ArithSqMat3x3Float<V,M>& m)
    {
      return 
        m[0][0]*m[1][1]*m[2][2] +
        m[0][1]*m[1][2]*m[2][0] +
        m[0][2]*m[1][0]*m[2][1] -
        m[0][2]*m[1][1]*m[2][0] -
        m[0][0]*m[1][2]*m[2][1] -
        m[0][1]*m[1][0]*m[2][2] ;
    }


}
#endif