Subversion Repositories gelsvn

Rev

Rev 12 | Go to most recent revision | Details | Last modification | View Log | RSS feed

Rev Author Line No. Line
2 bj 1 
#ifndef __CGLA_ARITHSQMAT3X3FLOAT_H__
2 
#define __CGLA_ARITHSQMAT3X3FLOAT_H__
3 
 
4 
#include "ExceptionStandard.h"
5 
#include "CGLA.h"
6 
#include "Vec3f.h"
7 
#include "ArithSqMatFloat.h"
8 
 
9 
namespace CGLA {
10 
 
11 
	CGLA_DERIVEEXCEPTION(Mat3x3fException)
12 
	CGLA_DERIVEEXCEPTION(Mat3x3fSingular)
13 
 
14 
	/** 3 by 3 float matrix.
15 
			This class will typically be used for rotation or
16 
			scaling matrices for 3D vectors. */
17 
	template<class V, class M>
18 
	class ArithSqMat3x3Float: public ArithSqMatFloat<V, M, 3>
19 
	{
20 
	public:
21 
 
22 
		/// Vector type
23 
		typedef V VectorType;
24 
 
25 
		/// The type of a matrix element
26 
		typedef typename V::ScalarType ScalarType;
27 
 
28 
	public:
29 
 
30 
		/// Construct matrix from 3 Vec3f vectors.
31 
		ArithSqMat3x3Float(V _a, V _b, V _c):
32 
			ArithSqMatFloat<V, M, 3> (_a,_b,_c) {}
33 
 
34 
		/// Construct the 0 matrix
35 
		ArithSqMat3x3Float() {}
36 
 
37 
		/// Construct a matrix from a single scalar value.
38 
		explicit ArithSqMat3x3Float(ScalarType a):
39 
			ArithSqMatFloat<V, M, 3>(a) {}
40 
 
41 
	};
42 
 
43 
	/// Invert 3x3 matrix
44 
	template<class V, class M>
45 
	M invert(const ArithSqMat3x3Float<V,M>&);
46 
 
47 
	/** Compute determinant. There is a more generic function for
48 
			computing determinants of square matrices (ArithSqMat). This one
49 
			is faster but works only on Mat3x3f */
50 
	template<class V, class M>
51 
	inline
52 
	typename ArithSqMat3x3Float<V,M>::ScalarType
53 
 	determinant(const ArithSqMat3x3Float<V,M>& m)
54 
	{
55 
		return
56 
			m[0][0]*m[1][1]*m[2][2] +
57 
			m[0][1]*m[1][2]*m[2][0] +
58 
			m[0][2]*m[1][0]*m[2][1] -
59 
			m[0][2]*m[1][1]*m[2][0] -
60 
			m[0][0]*m[1][2]*m[2][1] -
61 
			m[0][1]*m[1][0]*m[2][2] ;
62 
	}
63 
 
64 
 
65 
}
66 
#endif
67 
 
68 
 
69 
 
70 
 
71 
 
72 
 
73