Subversion Repositories gelsvn

Rev

Rev 569 | Only display areas with differences | Ignore whitespace | Details | Blame | Last modification | View Log | RSS feed

Rev 569 Rev 595
- 
 
 1 
/* ----------------------------------------------------------------------- *
- 
 
 2 
 * This file is part of GEL, http://www.imm.dtu.dk/GEL
- 
 
 3 
 * Copyright (C) the authors and DTU Informatics
- 
 
 4 
 * For license and list of authors, see ../../doc/intro.pdf
- 
 
 5 
 * @brief Abstract 3D floating point vector class
- 
 
 6 
 * ----------------------------------------------------------------------- */
- 
 
 7 
 
- 
 
 8 
/** @file ArithVec3Float.h
- 
 
 9 
 * @brief Abstract 3D floating point vector class
- 
 
 10 
 */
- 
 
 11 
 
1 
#ifndef __CGLA__ARITHVEC3FLOAT_H__
 12 
#ifndef __CGLA__ARITHVEC3FLOAT_H__
2 
#define __CGLA__ARITHVEC3FLOAT_H__
 13 
#define __CGLA__ARITHVEC3FLOAT_H__
3 
 
 14 
 
4 
#include "ArithVecFloat.h"
 15 
#include "ArithVecFloat.h"
5 
 
 16 
 
6 
namespace CGLA {
 17 
namespace CGLA {
7 
 
 18 
 
8 
	template<class T, class V>
 19 
	template<class T, class V>
9 
	class ArithVec3Float: public ArithVecFloat<T,V,3>
 20 
	class ArithVec3Float: public ArithVecFloat<T,V,3>
10 
	{
 21 
	{
11 
	public:
 22 
	public:
12 
 
 23 
 
13 
		/// Construct a 3D float vector.
 24 
		/// Construct a 3D float vector.
14 
		ArithVec3Float(T a, T b, T c): ArithVecFloat<T,V,3>(a,b,c) {}
 25 
		ArithVec3Float(T a, T b, T c): ArithVecFloat<T,V,3>(a,b,c) {}
15 
 
 26 
 
16 
		/// Construct a 3D float vector.
 27 
		/// Construct a 3D float vector.
17 
		ArithVec3Float() {}
 28 
		ArithVec3Float() {}
18 
 
 29 
 
19 
		/** Get the vector in spherical coordinates.
 30 
		/** Get the vector in spherical coordinates.
20 
				The first argument (theta) is inclination from the vertical axis.
 31 
				The first argument (theta) is inclination from the vertical axis.
21 
				The second argument (phi) is the angle of rotation about the vertical
32 
				The second argument (phi) is the angle of rotation about the vertical
22 
				axis. The third argument (r) is the length of the vector. */
 33 
				axis. The third argument (r) is the length of the vector. */
23 
		void get_spherical( T&, T&, T& ) const;
 34 
		void get_spherical( T&, T&, T& ) const;
24 
 
 35 
 
25 
		/** Assign the vector in spherical coordinates.
 36 
		/** Assign the vector in spherical coordinates.
26 
				The first argument (theta) is inclination from the vertical axis.
 37 
				The first argument (theta) is inclination from the vertical axis.
27 
				The second argument (phi) is the angle of rotation about the vertical
38 
				The second argument (phi) is the angle of rotation about the vertical
28 
				axis. The third argument (r) is the length of the vector. */
 39 
				axis. The third argument (r) is the length of the vector. */
29 
		void set_spherical( T, T, T);
 40 
		void set_spherical( T, T, T);
30
41
31 
	};
 42 
	};
32 
 
 43 
 
33 
	/// Returns cross product of arguments
 44 
	/// Returns cross product of arguments
34 
	template<class T, class V>
 45 
	template<class T, class V>
35 
	inline V cross( const ArithVec3Float<T,V>& x,
46 
	inline V cross( const ArithVec3Float<T,V>& x,
36 
									const ArithVec3Float<T,V>& y )
47 
									const ArithVec3Float<T,V>& y )
37 
	{
 48 
	{
38 
		return V( x[1] * y[2] - x[2] * y[1],
49 
		return V( x[1] * y[2] - x[2] * y[1],
39 
							x[2] * y[0] - x[0] * y[2],
50 
							x[2] * y[0] - x[0] * y[2],
40 
							x[0] * y[1] - x[1] * y[0] );
 51 
							x[0] * y[1] - x[1] * y[0] );
41 
	}
 52 
	}
42 
 
 53 
 
43 
	/** Compute basis of orthogonal plane.
 54 
	/** Compute basis of orthogonal plane.
44 
			Given a vector compute two vectors that are orthogonal to it and
55 
			Given a vector compute two vectors that are orthogonal to it and
45 
			to each other. */
 56 
			to each other. */
46 
	template<class T, class V>
 57 
	template<class T, class V>
47 
	void orthogonal(const ArithVec3Float<T,V>&,
 58 
	void orthogonal(const ArithVec3Float<T,V>&,
48 
									ArithVec3Float<T,V>&,
 59 
									ArithVec3Float<T,V>&,
49 
									ArithVec3Float<T,V>&);
 60 
									ArithVec3Float<T,V>&);
50 
 
 61 
 
51 
  /** Build an orthonormal basis from a 3d unit vector [Frisvad 2012].
 62 
  /** Build an orthonormal basis from a 3d unit vector [Frisvad 2012].
52 
      Given a unit vector compute two unit vectors that are orthogonal to
 63 
      Given a unit vector compute two unit vectors that are orthogonal to
53 
      it and to each other. */
 64 
      it and to each other. */
54 
  template<class T, class V>
 65 
  template<class T, class V>
55 
  void onb(const ArithVec3Float<T,V>&,
 66 
  void onb(const ArithVec3Float<T,V>&,
56 
           ArithVec3Float<T,V>&,
 67 
           ArithVec3Float<T,V>&,
57 
           ArithVec3Float<T,V>&);
 68 
           ArithVec3Float<T,V>&);
58 
}
 69 
}
59 
 
 70 
 
60 
#endif
 71 
#endif
61 
 
 72 
 
62 
 
 73