commit!

Committer:
narshu
Date:
Fri Jun 15 20:40:17 2012 +0000
Revision:
0:42026f893a2d

        

Who changed what in which revision?

UserRevisionLine numberNew contents of line
narshu 0:42026f893a2d 1 #ifndef GEOMETRYFUNCS_H
narshu 0:42026f893a2d 2 #define GEOMETRYFUNCS_H
narshu 0:42026f893a2d 3
narshu 0:42026f893a2d 4 #include <tvmet/Matrix.h>
narshu 0:42026f893a2d 5
narshu 0:42026f893a2d 6 template <typename T>
narshu 0:42026f893a2d 7 Matrix <T, 2, 2> Rotmatrix(T theta) {
narshu 0:42026f893a2d 8 Matrix <T, 2, 2> outmatrix;
narshu 0:42026f893a2d 9 outmatrix = cos(theta), -sin(theta),
narshu 0:42026f893a2d 10 sin(theta), cos(theta);
narshu 0:42026f893a2d 11 return outmatrix;
narshu 0:42026f893a2d 12 }
narshu 0:42026f893a2d 13
narshu 0:42026f893a2d 14 // rectifies angle to range -PI to PI
narshu 0:42026f893a2d 15 template <typename T>
narshu 0:42026f893a2d 16 T rectifyAng (T ang_in) {
narshu 0:42026f893a2d 17 ang_in -= (floor(ang_in/(2*PI)))*2*PI;
narshu 0:42026f893a2d 18 if (ang_in < -PI) {
narshu 0:42026f893a2d 19 ang_in += 2*PI;
narshu 0:42026f893a2d 20 }
narshu 0:42026f893a2d 21 if (ang_in > PI) {
narshu 0:42026f893a2d 22 ang_in -= 2*PI;
narshu 0:42026f893a2d 23 }
narshu 0:42026f893a2d 24
narshu 0:42026f893a2d 25 return ang_in;
narshu 0:42026f893a2d 26 }
narshu 0:42026f893a2d 27
narshu 0:42026f893a2d 28 #endif //GEOMETRYFUNCS_H