| 1 | #ifndef VECTOR_HPP_ | 
|---|
| 2 | #define VECTOR_HPP_ | 
|---|
| 3 |  | 
|---|
| 4 | using namespace std; | 
|---|
| 5 |  | 
|---|
| 6 | /*********************************************** includes ***********************************/ | 
|---|
| 7 |  | 
|---|
| 8 | // include config.h | 
|---|
| 9 | #ifdef HAVE_CONFIG_H | 
|---|
| 10 | #include <config.h> | 
|---|
| 11 | #endif | 
|---|
| 12 |  | 
|---|
| 13 | #include <iosfwd> | 
|---|
| 14 |  | 
|---|
| 15 | #include <memory> | 
|---|
| 16 | #include <vector> | 
|---|
| 17 |  | 
|---|
| 18 | #include "defs.hpp" | 
|---|
| 19 | #include "Space.hpp" | 
|---|
| 20 |  | 
|---|
| 21 | /********************************************** declarations *******************************/ | 
|---|
| 22 |  | 
|---|
| 23 | class Vector; | 
|---|
| 24 | class Matrix; | 
|---|
| 25 | struct VectorContent; | 
|---|
| 26 |  | 
|---|
| 27 | typedef std::vector<Vector> pointset; | 
|---|
| 28 |  | 
|---|
| 29 | /** Single vector. | 
|---|
| 30 | * basically, just a x[3] but with helpful functions | 
|---|
| 31 | */ | 
|---|
| 32 | class Vector : public Space{ | 
|---|
| 33 | friend Vector operator*(const Matrix&,const Vector&); | 
|---|
| 34 | friend class Matrix; | 
|---|
| 35 | public: | 
|---|
| 36 | Vector(); | 
|---|
| 37 | Vector(const double x1, const double x2, const double x3); | 
|---|
| 38 | Vector(const Vector& src); | 
|---|
| 39 | virtual ~Vector(); | 
|---|
| 40 |  | 
|---|
| 41 | // Method implemented by forwarding to the Representation | 
|---|
| 42 |  | 
|---|
| 43 | double DistanceSquared(const Vector &y) const; | 
|---|
| 44 | double DistanceToSpace(const Space& space) const; | 
|---|
| 45 | double ScalarProduct(const Vector &y) const; | 
|---|
| 46 | double Angle(const Vector &y) const; | 
|---|
| 47 | bool IsZero() const; | 
|---|
| 48 | bool IsOne() const; | 
|---|
| 49 | bool IsNormalTo(const Vector &normal) const; | 
|---|
| 50 | bool IsEqualTo(const Vector &a) const; | 
|---|
| 51 |  | 
|---|
| 52 | void AddVector(const Vector &y); | 
|---|
| 53 | void SubtractVector(const Vector &y); | 
|---|
| 54 | void VectorProduct(const Vector &y); | 
|---|
| 55 | void ProjectOntoPlane(const Vector &y); | 
|---|
| 56 | void ProjectIt(const Vector &y); | 
|---|
| 57 | Vector Projection(const Vector &y) const; | 
|---|
| 58 | void ScaleAll(const double *factor); | 
|---|
| 59 | void ScaleAll(const Vector &factor); | 
|---|
| 60 | void Scale(const double factor); | 
|---|
| 61 | bool GetOneNormalVector(const Vector &x1); | 
|---|
| 62 | bool MakeNormalTo(const Vector &y1); | 
|---|
| 63 | std::pair<Vector,Vector> partition(const Vector&) const; | 
|---|
| 64 | std::pair<pointset,Vector> partition(const pointset&) const; | 
|---|
| 65 |  | 
|---|
| 66 | // Accessors ussually come in pairs... and sometimes even more than that | 
|---|
| 67 | double& operator[](size_t i); | 
|---|
| 68 | const double& operator[](size_t i) const; | 
|---|
| 69 | double& at(size_t i); | 
|---|
| 70 | const double& at(size_t i) const; | 
|---|
| 71 |  | 
|---|
| 72 | // Assignment operator | 
|---|
| 73 | Vector &operator=(const Vector& src); | 
|---|
| 74 |  | 
|---|
| 75 | // Access to internal structure | 
|---|
| 76 | VectorContent* get(); | 
|---|
| 77 |  | 
|---|
| 78 | // Methods that are derived directly from other methods | 
|---|
| 79 | double Norm() const; | 
|---|
| 80 | double NormSquared() const; | 
|---|
| 81 | void Normalize(); | 
|---|
| 82 | void Zero(); | 
|---|
| 83 | void One(const double one); | 
|---|
| 84 | void LinearCombinationOfVectors(const Vector &x1, const Vector &x2, const Vector &x3, const double * const factors); | 
|---|
| 85 |  | 
|---|
| 86 | // operators for mathematical operations | 
|---|
| 87 | bool operator==(const Vector& b) const; | 
|---|
| 88 | bool operator!=(const Vector& b) const; | 
|---|
| 89 | const Vector& operator+=(const Vector& b); | 
|---|
| 90 | const Vector& operator-=(const Vector& b); | 
|---|
| 91 | Vector const operator+(const Vector& b) const; | 
|---|
| 92 | Vector const operator-(const Vector& b) const; | 
|---|
| 93 |  | 
|---|
| 94 | // Methods inherited from Space | 
|---|
| 95 | virtual double distance(const Vector &point) const; | 
|---|
| 96 | virtual Vector getClosestPoint(const Vector &point) const; | 
|---|
| 97 |  | 
|---|
| 98 | protected: | 
|---|
| 99 |  | 
|---|
| 100 | private: | 
|---|
| 101 | Vector(VectorContent *); | 
|---|
| 102 | VectorContent *content; | 
|---|
| 103 |  | 
|---|
| 104 | }; | 
|---|
| 105 |  | 
|---|
| 106 | // some commonly used and fixed vectors | 
|---|
| 107 | const extern Vector zeroVec; | 
|---|
| 108 | const extern Vector e1; | 
|---|
| 109 | const extern Vector e2; | 
|---|
| 110 | const extern Vector e3; | 
|---|
| 111 |  | 
|---|
| 112 | ostream & operator << (ostream& ost, const Vector &m); | 
|---|
| 113 | const Vector& operator*=(Vector& a, const double m); | 
|---|
| 114 | Vector const operator*(const Vector& a, const double m); | 
|---|
| 115 | Vector const operator*(const double m, const Vector& a); | 
|---|
| 116 |  | 
|---|
| 117 | #endif /*VECTOR_HPP_*/ | 
|---|