WebSVN – gelsvn – Diff – /branches/imesh-and-hmesh/src/CGLA/ArithSqMatFloat.h


#ifndef __CGLA_ARITHSQMATFLOAT_H__
#define __CGLA_ARITHSQMATFLOAT_H__
 
#include "ArithMatFloat.h"
 
namespace CGLA 
{
 
  /** Template for square matrices.
      Some functions like trace and determinant work only on
      square matrices. To express this in the class hierarchy,
      ArithSqMatFloat was created. ArithSqMatFloat is derived from ArithMat
      and contains a few extra facilities applicable only to
      square matrices.
  */
  template <class VT, class MT, int ROWS>
    class ArithSqMatFloat: public ArithMatFloat<VT,VT,MT,ROWS> 
    { 
    public:
 
      /// Vector type
      typedef VT VectorType;
 
      /// The type of a matrix element
      typedef typename VT::ScalarType ScalarType;
 
    protected:
 
      /// Construct 0 matrix
      ArithSqMatFloat() {}
 
      /// Construct matrix where all values are equal to constructor argument.
      explicit ArithSqMatFloat(ScalarType _a):
	ArithMatFloat<VT,VT,MT,ROWS>(_a) {}
 
      /// Construct 2x2 Matrix from two vectors
      ArithSqMatFloat(VT _a, VT _b): 
	ArithMatFloat<VT,VT,MT,ROWS>(_a,_b) {}
 
      /// Construct 3x3 Matrix from three vectors
      ArithSqMatFloat(VT _a, VT _b, VT _c): 
	ArithMatFloat<VT,VT,MT,ROWS>(_a,_b,_c) {}
 
      /// Construct 4x4 Matrix from four vectors
      ArithSqMatFloat(VT _a, VT _b, VT _c, VT _d): 
	ArithMatFloat<VT,VT,MT,ROWS>(_a,_b,_c,_d) {}
		
    public:
 
      /** Assignment multiplication of matrices. 
	  This function is not very efficient. This because we need a temporary
	  matrix anyway, so it can't really be made efficient. */
      const MT& operator*=(const MT& m2)
	{
	  (*this) = (*this) * m2;
	  return static_cast<const MT&>(*this);
	}
		
      void identity()
	{
	  for(int i=0;i<ROWS;++i)
	    {
	      for(int j=0;j<ROWS;++j)
		(*this)[i][j] = ScalarType(0);
	      (*this)[i][i] = ScalarType(1);
	    }
	}
 
		
    };
 
  /** Multiply two matrices derived from same type, 
      producing a new of same type */
  template <class VT, class MT, int ROWS>
    inline MT operator*(const ArithSqMatFloat<VT,MT,ROWS>& m1,
			const ArithSqMatFloat<VT,MT,ROWS>& m2) 
    {
      MT n;
      for(int i=0;i<ROWS;i++)
	for(int j=0;j<ROWS;j++)
	  {
	    n[i][j] = 0;
	    for(int k=0;k<ROWS;k++)
	      n[i][j] += m1[i][k] * m2[k][j]; 
	  }
      return n;
    }
 
  /** Compute the transpose of a square matrix. This function returns
      the transpose of its argument. */
  template <class VT, class MT, int ROWS>
    inline MT transpose(const ArithSqMatFloat<VT,MT,ROWS>& m) 
    {
      MT m_new;
      for(int i=0;i<MT::get_v_dim();i++)
	for(int j=0;j<MT::get_h_dim();j++)
	  m_new[i][j] = m[j][i];
      return m_new;
    }
 
  /// Compute trace. Works only for sq. matrices.
  template <class VT, class MT, int ROWS>
    inline typename MT::ScalarType trace(const ArithSqMatFloat<VT,MT,ROWS>& M)
    {
      typename ArithSqMatFloat<VT,MT,ROWS>::ScalarType s=0;
      for(int i=0;i<ROWS;i++)
	s += M[i][i];
      return s;
    }
 
}
#endif
 

Subversion Repositories gelsvn

(root)/branches/imesh-and-hmesh/src/CGLA/ArithSqMatFloat.h – Rev 12 → 22

Rev 12		Rev 22
1	`#ifndef __CGLA_ARITHSQMATFLOAT_H__`	1	`#ifndef __CGLA_ARITHSQMATFLOAT_H__`
2	`#define __CGLA_ARITHSQMATFLOAT_H__`	2	`#define __CGLA_ARITHSQMATFLOAT_H__`
3		3
4	`#include "ArithMatFloat.h"`	4	`#include "ArithMatFloat.h"`
5		5
6	`namespace CGLA`	6	`namespace CGLA`
7	`{`	7	`{`
8		8
9	`/** Template for square matrices.`	9	`/** Template for square matrices.`
10	`Some functions like trace and determinant work only on`	10	`Some functions like trace and determinant work only on`
11	`square matrices. To express this in the class hierarchy,`	11	`square matrices. To express this in the class hierarchy,`
12	`ArithSqMatFloat was created. ArithSqMatFloat is derived from ArithMat`	12	`ArithSqMatFloat was created. ArithSqMatFloat is derived from ArithMat`
13	`and contains a few extra facilities applicable only to`	13	`and contains a few extra facilities applicable only to`
14	`square matrices.`	14	`square matrices.`
15	`*/`	15	`*/`
16	`template <class VT, class MT, int ROWS>`	16	`template <class VT, class MT, int ROWS>`
17	`class ArithSqMatFloat: public ArithMatFloat<VT,VT,MT,ROWS>`	17	`class ArithSqMatFloat: public ArithMatFloat<VT,VT,MT,ROWS>`
18	`{`	18	`{`
19	`public:`	19	`public:`
20		20
21	`/// Vector type`	21	`/// Vector type`
22	`typedef VT VectorType;`	22	`typedef VT VectorType;`
23		23
24	`/// The type of a matrix element`	24	`/// The type of a matrix element`
25	`typedef typename VT::ScalarType ScalarType;`	25	`typedef typename VT::ScalarType ScalarType;`
26		26
27	`protected:`	27	`protected:`
28		28
29	`/// Construct 0 matrix`	29	`/// Construct 0 matrix`
30	`ArithSqMatFloat() {}`	30	`ArithSqMatFloat() {}`
31		31
32	`/// Construct matrix where all values are equal to constructor argument.`	32	`/// Construct matrix where all values are equal to constructor argument.`
33	`explicit ArithSqMatFloat(ScalarType _a):`	33	`explicit ArithSqMatFloat(ScalarType _a):`
34	`ArithMatFloat<VT,VT,MT,ROWS>(_a) {}`	34	`ArithMatFloat<VT,VT,MT,ROWS>(_a) {}`
35		35
36	`/// Construct 2x2 Matrix from two vectors`	36	`/// Construct 2x2 Matrix from two vectors`
37	`ArithSqMatFloat(VT _a, VT _b):`	37	`ArithSqMatFloat(VT _a, VT _b):`
38	`ArithMatFloat<VT,VT,MT,ROWS>(_a,_b) {}`	38	`ArithMatFloat<VT,VT,MT,ROWS>(_a,_b) {}`
39		39
40	`/// Construct 3x3 Matrix from three vectors`	40	`/// Construct 3x3 Matrix from three vectors`
41	`ArithSqMatFloat(VT _a, VT _b, VT _c):`	41	`ArithSqMatFloat(VT _a, VT _b, VT _c):`
42	`ArithMatFloat<VT,VT,MT,ROWS>(_a,_b,_c) {}`	42	`ArithMatFloat<VT,VT,MT,ROWS>(_a,_b,_c) {}`
43		43
44	`/// Construct 4x4 Matrix from four vectors`	44	`/// Construct 4x4 Matrix from four vectors`
45	`ArithSqMatFloat(VT _a, VT _b, VT _c, VT _d):`	45	`ArithSqMatFloat(VT _a, VT _b, VT _c, VT _d):`
46	`ArithMatFloat<VT,VT,MT,ROWS>(_a,_b,_c,_d) {}`	46	`ArithMatFloat<VT,VT,MT,ROWS>(_a,_b,_c,_d) {}`
47		47
48	`public:`	48	`public:`
49		49
50	`/** Assignment multiplication of matrices.`	50	`/** Assignment multiplication of matrices.`
51	`This function is not very efficient. This because we need a temporary`	51	`This function is not very efficient. This because we need a temporary`
52	`matrix anyway, so it can't really be made efficient. */`	52	`matrix anyway, so it can't really be made efficient. */`
53	`const MT& operator*=(const MT& m2)`	53	`const MT& operator*=(const MT& m2)`
54	`{`	54	`{`
55	`(this) = (this) * m2;`	55	`(this) = (this) * m2;`
56	`return static_cast<const MT&>(*this);`	56	`return static_cast<const MT&>(*this);`
57	`}`	57	`}`
58		58
59	`void identity()`	59	`void identity()`
60	`{`	60	`{`
61	`for(int i=0;i<ROWS;++i)`	61	`for(int i=0;i<ROWS;++i)`
62	`{`	62	`{`
63	`for(int j=0;j<ROWS;++j)`	63	`for(int j=0;j<ROWS;++j)`
64	`(*this)[i][j] = ScalarType(0);`	64	`(*this)[i][j] = ScalarType(0);`
65	`(*this)[i][i] = ScalarType(1);`	65	`(*this)[i][i] = ScalarType(1);`
66	`}`	66	`}`
67	`}`	67	`}`
68		68
69		69
70	`};`	70	`};`
71		71
72	`/** Multiply two matrices derived from same type,`	72	`/** Multiply two matrices derived from same type,`
73	`producing a new of same type */`	73	`producing a new of same type */`
74	`template <class VT, class MT, int ROWS>`	74	`template <class VT, class MT, int ROWS>`
75	`inline MT operator*(const ArithSqMatFloat<VT,MT,ROWS>& m1,`	75	`inline MT operator*(const ArithSqMatFloat<VT,MT,ROWS>& m1,`
76	`const ArithSqMatFloat<VT,MT,ROWS>& m2)`	76	`const ArithSqMatFloat<VT,MT,ROWS>& m2)`
77	`{`	77	`{`
78	`MT n;`	78	`MT n;`
79	`for(int i=0;i<ROWS;i++)`	79	`for(int i=0;i<ROWS;i++)`
80	`for(int j=0;j<ROWS;j++)`	80	`for(int j=0;j<ROWS;j++)`
81	`{`	81	`{`
82	`n[i][j] = 0;`	82	`n[i][j] = 0;`
83	`for(int k=0;k<ROWS;k++)`	83	`for(int k=0;k<ROWS;k++)`
84	`n[i][j] += m1[i][k] * m2[k][j];`	84	`n[i][j] += m1[i][k] * m2[k][j];`
85	`}`	85	`}`
86	`return n;`	86	`return n;`
87	`}`	87	`}`
88		88
89	`/** Compute the transpose of a square matrix. This function returns`	89	`/** Compute the transpose of a square matrix. This function returns`
90	`the transpose of its argument. */`	90	`the transpose of its argument. */`
91	`template <class VT, class MT, int ROWS>`	91	`template <class VT, class MT, int ROWS>`
92	`inline MT transpose(const ArithSqMatFloat<VT,MT,ROWS>& m)`	92	`inline MT transpose(const ArithSqMatFloat<VT,MT,ROWS>& m)`
93	`{`	93	`{`
94	`MT m_new;`	94	`MT m_new;`
95	`for(int i=0;i<MT::get_v_dim();i++)`	95	`for(int i=0;i<MT::get_v_dim();i++)`
96	`for(int j=0;j<MT::get_h_dim();j++)`	96	`for(int j=0;j<MT::get_h_dim();j++)`
97	`m_new[i][j] = m[j][i];`	97	`m_new[i][j] = m[j][i];`
98	`return m_new;`	98	`return m_new;`
99	`}`	99	`}`
100		100
101	`/// Compute trace. Works only for sq. matrices.`	101	`/// Compute trace. Works only for sq. matrices.`
102	`template <class VT, class MT, int ROWS>`	102	`template <class VT, class MT, int ROWS>`
103	`inline typename MT::ScalarType trace(const ArithSqMatFloat<VT,MT,ROWS>& M)`	103	`inline typename MT::ScalarType trace(const ArithSqMatFloat<VT,MT,ROWS>& M)`
104	`{`	104	`{`
105	`typename ArithSqMatFloat<VT,MT,ROWS>::ScalarType s=0;`	105	`typename ArithSqMatFloat<VT,MT,ROWS>::ScalarType s=0;`
106	`for(int i=0;i<ROWS;i++)`	106	`for(int i=0;i<ROWS;i++)`
107	`s += M[i][i];`	107	`s += M[i][i];`
108	`return s;`	108	`return s;`
109	`}`	109	`}`
110		110
111	`}`	111	`}`
112	`#endif`	112	`#endif`
113		113