yodaUtility/sgp4/cVector.cpp

//
// cVector.cpp
//
// Copyright (c) 2001-2003 Michael F. Henry
//
#include "stdafx.h"

#include "math.h"
#include "cVector.h"

//*****************************************************************************
// Multiply each component in the vector by 'factor'.
//*****************************************************************************
void cVector::Mul(double factor)
{
   m_x *= factor;
   m_y *= factor;
   m_z *= factor;
   m_w *= fabs(factor);
}

//*****************************************************************************
// Subtract a vector from this one.
//*****************************************************************************
void cVector::Sub(const cVector& vec)
{
   m_x -= vec.m_x;
   m_y -= vec.m_y;
   m_z -= vec.m_z;
   m_w -= vec.m_w;
}

//*****************************************************************************
// Calculate the angle between this vector and another
//*****************************************************************************
double cVector::Angle(const cVector& vec) const
{
  return acos(Dot(vec) / (Magnitude() * vec.Magnitude()));
}

//*****************************************************************************
//
//*****************************************************************************
double cVector::Magnitude() const
{
  return sqrt((m_x * m_x) + 
              (m_y * m_y) + 
              (m_z * m_z));
}

//*****************************************************************************
// Return the dot product
//*****************************************************************************
double cVector::Dot(const cVector& vec) const
{
   return (m_x * vec.m_x) +
          (m_y * vec.m_y) +
          (m_z * vec.m_z);
}
1	//
2	// cVector.cpp
3	//
4	// Copyright (c) 2001-2003 Michael F. Henry
5	//
6	#include "stdafx.h"
7
8	#include "math.h"
9	#include "cVector.h"
10
11	//*****************************************************************************
12	// Multiply each component in the vector by 'factor'.
13	//*****************************************************************************
14	void cVector::Mul(double factor)
15	{
16	m_x *= factor;
17	m_y *= factor;
18	m_z *= factor;
19	m_w *= fabs(factor);
20	}
21
22	//*****************************************************************************
23	// Subtract a vector from this one.
24	//*****************************************************************************
25	void cVector::Sub(const cVector& vec)
26	{
27	m_x -= vec.m_x;
28	m_y -= vec.m_y;
29	m_z -= vec.m_z;
30	m_w -= vec.m_w;
31	}
32
33	//*****************************************************************************
34	// Calculate the angle between this vector and another
35	//*****************************************************************************
36	double cVector::Angle(const cVector& vec) const
37	{
38	return acos(Dot(vec) / (Magnitude() * vec.Magnitude()));
39	}
40
41	//*****************************************************************************
42	//
43	//*****************************************************************************
44	double cVector::Magnitude() const
45	{
46	return sqrt((m_x * m_x) +
47	(m_y * m_y) +
48	(m_z * m_z));
49	}
50
51	//*****************************************************************************
52	// Return the dot product
53	//*****************************************************************************
54	double cVector::Dot(const cVector& vec) const
55	{
56	return (m_x * vec.m_x) +
57	(m_y * vec.m_y) +
58	(m_z * vec.m_z);
59	}