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