Subversion Repositories gelsvn

#ifndef __CGLA_UNITVECTOR_H_

#ifndef __CGLA_UNITVECTOR_H__

#define __CGLA_UNITVECTOR_H_

#define __CGLA_UNITVECTOR_H__

	namespace TT=TableTrigonometry;

  namespace TT=TableTrigonometry;

	/** The UnitVector stores a unit length vector as two angles.

  /** The UnitVector stores a unit length vector as two angles.

	A vector stored as two (fix point) angles is much smaller than

  A vector stored as two (fix point) angles is much smaller than

	vector stored in the usual way. On a 32 bit architecture this

  vector stored in the usual way. On a 32 bit architecture this

	class should take up four bytes. not too bad. */

  class should take up four bytes. not too bad. */

	class UnitVector

  class UnitVector

	{

    {

		TT::Angle theta, phi;

      TT::Angle theta, phi;

		void encode(const Vec3f& v)

      void encode(const Vec3f& v)

		{

	{

			theta = TT::t_atan(std::sqrt(CGLA::sqr(v[0])+CGLA::sqr(v[1])), v[2]);

	  theta = TT::t_atan(std::sqrt(CGLA::sqr(v[0])+CGLA::sqr(v[1])), v[2]);

			phi   = TT::t_atan(v[0],v[1]);

	  phi   = TT::t_atan(v[0],v[1]);

		}

	}

	public:

    public:

		/// Construct unitvector from normal vector

      /// Construct unitvector from normal vector

		explicit UnitVector(const Vec3f& v) {encode(v);}

      explicit UnitVector(const Vec3f& v) {encode(v);}

		/// Construct default unit vector

      /// Construct default unit vector

		explicit UnitVector(): theta(0), phi(0) {}

      explicit UnitVector(): theta(0), phi(0) {}

		/// Get theta angle

      /// Get theta angle

		float t() const {return TT::angle2float(theta);}

      float t() const {return TT::angle2float(theta);}

	/// Inline output operator.

      /// Reconstruct Vec3f from unit vector

	inline std::ostream& operator<<(std::ostream& os, const UnitVector& u)

      operator Vec3f() const

		os << "<v=" << u.t() << " h=" << u.f() << ">";

	  float costf = TT::t_cos(theta);

		return os;

		       TT::t_sin(theta));