+#if 0
+
//
// VECTOR.H - vector class definition
//
};
#endif // __VECTOR_H__
+
+#else
+
+//
+// vector.h (Last modified: 6/28/2001)
+//
+// Various structures used for 3 dimensional imaging
+//
+// by James L. Hammons
+// (C) 2001 Underground Software
+//
+
+#ifndef __VECTOR_H__
+#define __VECTOR_H__
+
+// What we'll do here is create the vector type and use typedef to alias Point to it. Yeah, that's it.
+
+class Vector
+{
+ public:
+// Vector(double xx = 0, double yy = 0, double zz = 0);
+ Vector(double x1 = 0, double y1 = 0, double z1 = 0, double x2 = 0, double y2 = 0, double z2 = 0);
+ Vector(Vector head, Vector tail); // Create vector from two points
+ Vector operator=(Vector const v);
+ Vector operator+(Vector const v);
+ Vector operator-(Vector const v);
+ Vector operator-(void); // Unary negation
+ Vector operator*(double const v); // Vector times constant (double)
+ Vector operator*(float const v); // Vector times constant (float)
+ Vector operator/(double const v); // Vector divided by constant (double)
+ Vector operator/(float const v); // Vector divided by constant (float)
+ Vector operator*(Vector const v); // Vector product
+ double Dot(Vector const v); // Dot product
+
+ Vector& operator*=(double const v); // Vector times constant self-assignment
+ Vector& operator/=(double const v); // Vector divided by constant self-assignment
+ Vector& operator+=(Vector const v); // Vector plus Vector self-assignment
+ Vector& operator+=(double const v); // Vector plus constant self-assignment
+ Vector& operator-=(Vector const v); // Vector minus Vector self-assignment
+ Vector& operator-=(double const v); // Vector minus constant self-assignment
+
+ bool operator==(Vector const v); // Check for equality
+ bool operator!=(Vector const v); // Check for inequality
+
+ Vector Unit(void);
+ double Magnitude(void);
+ double Angle(void); // Angle of this vector WRT positive X-axis
+ double Angle(Vector v); // Angle between these two vectors
+ bool isZero(double epsilon = 1e-6);
+
+ // Class methods
+
+ static double Dot(Vector v1, Vector v2);
+ static double Magnitude(Vector v1, Vector v2);
+
+ public:
+ double x, y, z;
+};
+
+typedef Vector Point;
+
+#endif // __VECTOR_H__
+
+#endif
\ No newline at end of file