angle.cpp

Go to the documentation of this file.

#ifndef ANGLE_IMPLEMENTATION
#define ANGLE_IMPLEMENTATION

/*****************************************************************************\
*                                                                             *
*  Name   : angle                                                             *
*  Author : Chris Koeritz                                                     *
*                                                                             *
*******************************************************************************
* Copyright (c) 1992-$now By Author.  This program is free software; you can  *
* redistribute it and/or modify it under the terms of the GNU General Public  *
* License as published by the Free Software Foundation; either version 2 of   *
* the License or (at your option) any later version.  This is online at:      *
*     http://www.fsf.org/copyleft/gpl.html                                    *
* Please send any updates to: fred@gruntose.com                               *
\*****************************************************************************/

#include "angle.h"

#include <basis/definitions.h>
#include <basis/guards.h>

#include <math.h>

namespace geometric {

template <class contents>
angle<contents>::angle(contents a, angular_units unit) { set(a, unit); }

template <class contents>
angle<contents> angle<contents>::operator - (void) const
{ angle<contents> to_return(*this); to_return *= -1; return to_return; }

template <class contents>
angle<contents> angle<contents>::operator + (const angle<contents> &a) const
{ angle<contents> to_return(*this); to_return += a; return to_return; }

template <class contents>
angle<contents> angle<contents>::operator - (const angle<contents> &a) const
{ angle<contents> to_return(*this); to_return -= a; return to_return; }

template <class contents>
angle<contents> angle<contents>::operator * (contents to_multiply) const
{
  angle<contents> to_return(*this);
  to_return *= to_multiply;
  return to_return;
}

template <class contents>
angle<contents> angle<contents>::operator / (contents to_divide) const
{ angle<contents> to_return(*this); to_return /= to_divide; return to_return; }

template <class contents>
angle<contents> &angle<contents>::operator += (const angle<contents> &a)
{ _theta += a._theta; return *this; }

template <class contents>
angle<contents> &angle<contents>::operator -= (const angle<contents> &a)
{ _theta -= a._theta; return *this; }

template <class contents>
angle<contents> &angle<contents>::operator *= (contents f)
{ _theta *= f; return *this; }

template <class contents>
angle<contents> &angle<contents>::operator /= (contents f)
{ _theta /= f; return *this; }

template <class contents>
contents angle<contents>::sine() const { return sin(_theta); }

template <class contents>
contents angle<contents>::cosine() const { return cos(_theta); }

template <class contents>
contents angle<contents>::tangent() const { return tan(_theta); }

template <class contents>
void angle<contents>::pack(byte_array &packed_form) const
{ basis::attach(packed_form, _theta); }

template <class contents>
bool angle<contents>::unpack(byte_array &packed_form)
{ return basis::detach(packed_form, _theta); }

template <class contents>
contents angle<contents>::to_internal(contents a, angular_units unit) const
{
  switch(unit) {
    case RADIANS: return a;
    case DEGREES: return a * PI_APPROX / 180.0;
    default: return 0;
  }
}

template <class contents>
contents angle<contents>::from_internal(contents a, angular_units unit) const
{
  switch(unit) {
    case RADIANS: return a;
    case DEGREES: return a * 180.0 / PI_APPROX;
    default: return 0;
  }
}

template <class contents>
void angle<contents>::set(contents a, angular_units unit)
{ _theta = to_internal(a, unit); }

template <class contents>
contents angle<contents>::get(angular_units unit) const
{ return from_internal(_theta, unit); }

template <class contents>
angle<contents> angle<contents>::arccosine(contents adjacent,
    contents hypotenuse, outcome &retval)
{
  contents d = adjacent / hypotenuse;
  retval = common::BAD_INPUT;
  bounds_return(d, -1.0, 1.0, angle<contents>());
  retval = common::OKAY;
  return angle<contents>(acos(d), angle<contents>::RADIANS);
}

template <class contents>
angle<contents> angle<contents>::arcsine(contents opposite, contents hypotenuse,
    outcome &retval)
{
  contents d = opposite / hypotenuse;
  retval = common::BAD_INPUT;
  bounds_return(d, -1.0, 1.0, angle<contents>());
  retval = common::OKAY;
  return angle<contents>(asin(d), angle<contents>::RADIANS);
}

template <class contents>
angle<contents> angle<contents>::arctangent(contents adjacent,
    contents opposite, outcome &retval)
{
  retval = common::BAD_INPUT;
  if ( (adjacent == 0.0) && (opposite == 0.0) ) return angle<contents>();
  retval = common::OKAY;
  return angle<contents>(atan2(opposite, adjacent), angle<contents>::RADIANS);
}

} // namespace.

#endif

---