Subversion Repositories gelsvn

Rev

Rev 89 | Go to most recent revision | Details | Compare with Previous | Last modification | View Log | RSS feed

Rev Author Line No. Line
2 bj 1 
#ifndef __CGLA_MAT4X4_H__
2 
#define __CGLA_MAT4X4_H__
3 
 
4 
#include "ExceptionStandard.h"
5 
#include "CGLA.h"
6 
#include "Vec3f.h"
7 
#include "Vec3Hf.h"
8 
#include "Vec4f.h"
9 
#include "ArithSqMat4x4Float.h"
10 
 
11 
 
12 jab 12 
namespace CGLA
13 
{
2 bj 14 
 
89 jab 15 
  /** \brief 4x4 float matrix.
12 jab 16 
      This class is useful for transformations such as perspective projections
17 
      or translation where 3x3 matrices do not suffice. */
18 
  class Mat4x4f: public ArithSqMat4x4Float<Vec4f, Mat4x4f>
19 
    {
20 
    public:
2 bj 21 
 
12 jab 22 
      /// Construct a Mat4x4f from four Vec4f vectors
23 
      Mat4x4f(Vec4f _a, Vec4f _b, Vec4f _c, Vec4f _d):
24 
	ArithSqMat4x4Float<Vec4f, Mat4x4f> (_a,_b,_c,_d) {}
2 bj 25 
 
12 jab 26 
      /// Construct the NaN matrix
27 
      Mat4x4f() {}
2 bj 28 
 
12 jab 29 
      /// Construct a matrix with identical elements.
43 jrf 30 
      explicit Mat4x4f(float a) : ArithSqMat4x4Float<Vec4f, Mat4x4f> (a) {}
12 jab 31 
    };
10 jab 32 
 
12 jab 33 
  /// Create a rotation _matrix. Rotates about one of the major axes.
34 
  Mat4x4f rotation_Mat4x4f(CGLA::Axis axis, float angle);
10 jab 35 
 
12 jab 36 
  /// Create a translation matrix
37 
  Mat4x4f translation_Mat4x4f(const Vec3f&);
2 bj 38 
 
12 jab 39 
  /// Create a scaling matrix.
40 
  Mat4x4f scaling_Mat4x4f(const Vec3f&);
2 bj 41 
 
12 jab 42 
  /// Create an identity matrix.
43 
  inline Mat4x4f identity_Mat4x4f()
44 
    {
501 jrf 45 
      return Mat4x4f(Vec4f(1.0f,0.0f,0.0f,0.0f),
46 
		     Vec4f(0.0f,1.0f,0.0f,0.0f),
47 
		     Vec4f(0.0f,0.0f,1.0f,0.0f),
48 
		     Vec4f(0.0f,0.0f,0.0f,1.0f));
12 jab 49 
    }
50 jrf 50 
 
51 
  /** Compute inverse assuming that the upper-left 3x3 sub-matrix is
52 
      orthonormal (which is the case if the transformation is only
53 
      a concatenation of rotations and translations).
54 
  */
55 
  inline Mat4x4f invert_ortho(const Mat4x4f& m)
56 
  {
57 
    Vec3f rx(m[0][0], m[1][0], m[2][0]);
58 
    Vec3f ry(m[0][1], m[1][1], m[2][1]);
59 
    Vec3f rz(m[0][2], m[1][2], m[2][2]);
60 
    Vec3f t(m[0][3], m[1][3], m[2][3]);
61 
 
62 
    return Mat4x4f(Vec4f(rx, -dot(t, rx)),
63 
		   Vec4f(ry, -dot(t, ry)),
64 
		   Vec4f(rz, -dot(t, rz)),
501 jrf 65 
		   Vec4f(0.0f, 0.0f, 0.0f, 1.0f));
50 jrf 66 
  }
2 bj 67 
}
68 
#endif
69 
 
70 
 
71 
 
72 
 
73 
 
74 
 
75