| [342f33f] | 1 | #ifndef VECTOR_HPP_
 | 
|---|
 | 2 | #define VECTOR_HPP_
 | 
|---|
 | 3 | 
 | 
|---|
| [e9b8bb] | 4 | class Vector;
 | 
|---|
| [7d96c6] | 5 | 
 | 
|---|
| [342f33f] | 6 | /** Single vector.
 | 
|---|
 | 7 |  * basically, just a x[3] but with helpful functions
 | 
|---|
 | 8 |  */
 | 
|---|
| [e9b8bb] | 9 | class Vector {   
 | 
|---|
| [342f33f] | 10 |   public:
 | 
|---|
 | 11 |     double x[NDIM];
 | 
|---|
 | 12 | 
 | 
|---|
| [e9b8bb] | 13 |   Vector();
 | 
|---|
 | 14 |   Vector(double x1, double x2, double x3);
 | 
|---|
 | 15 |   ~Vector();
 | 
|---|
| [342f33f] | 16 | 
 | 
|---|
| [e9b8bb] | 17 |   double Distance(const Vector *y) const;
 | 
|---|
 | 18 |   double PeriodicDistance(const Vector *y, const double *cell_size) const;
 | 
|---|
 | 19 |   double ScalarProduct(const Vector *y) const;
 | 
|---|
 | 20 |   double Projection(const Vector *y) const;
 | 
|---|
| [342f33f] | 21 |   double Norm() const ;
 | 
|---|
| [e9b8bb] | 22 |   double Angle(Vector *y) const;
 | 
|---|
| [342f33f] | 23 | 
 | 
|---|
| [e9b8bb] | 24 |   void AddVector(const Vector *y);
 | 
|---|
 | 25 |   void SubtractVector(const Vector *y);
 | 
|---|
 | 26 |   void CopyVector(const Vector *y);
 | 
|---|
 | 27 |   void RotateVector(const Vector *y, const double alpha);
 | 
|---|
 | 28 |   void ProjectOntoPlane(const Vector *y);
 | 
|---|
| [3dfad8] | 29 |   void Zero(); 
 | 
|---|
| [498a9f] | 30 |   void One(double one);
 | 
|---|
 | 31 |   void Init(double x1, double x2, double x3);
 | 
|---|
| [342f33f] | 32 |   void Normalize();
 | 
|---|
| [e9b8bb] | 33 |   void Translate(const Vector *x);
 | 
|---|
 | 34 |   void Mirror(const Vector *x);
 | 
|---|
| [342f33f] | 35 |   void Scale(double **factor);
 | 
|---|
 | 36 |   void Scale(double *factor);
 | 
|---|
 | 37 |   void Scale(double factor);
 | 
|---|
 | 38 |   void MatrixMultiplication(double *M);
 | 
|---|
 | 39 |   void InverseMatrixMultiplication(double *M);
 | 
|---|
 | 40 |   void KeepPeriodic(ofstream *out, double *matrix);
 | 
|---|
| [e9b8bb] | 41 |   void LinearCombinationOfVectors(const Vector *x1, const Vector *x2, const Vector *x3, double *factors);
 | 
|---|
| [342f33f] | 42 |   
 | 
|---|
| [e9b8bb] | 43 |   double CutsPlaneAt(Vector *A, Vector *B, Vector *C);
 | 
|---|
 | 44 |   bool GetOneNormalVector(const Vector *x1);
 | 
|---|
 | 45 |   bool MakeNormalVector(const Vector *y1);
 | 
|---|
 | 46 |   bool MakeNormalVector(const Vector *y1, const Vector *y2);
 | 
|---|
 | 47 |   bool MakeNormalVector(const Vector *x1, const Vector *x2, const Vector *x3);
 | 
|---|
 | 48 |   bool SolveSystem(Vector *x1, Vector *x2, Vector *y, double alpha, double beta, double c);
 | 
|---|
 | 49 |   bool LSQdistance(Vector **vectors, int dim);
 | 
|---|
| [342f33f] | 50 | 
 | 
|---|
 | 51 |   void AskPosition(double *cell_size, bool check);
 | 
|---|
 | 52 |   bool Output(ofstream *out) const;
 | 
|---|
 | 53 | };
 | 
|---|
 | 54 | 
 | 
|---|
| [e9b8bb] | 55 | ofstream& operator<<(ofstream& ost, Vector& m);
 | 
|---|
 | 56 | Vector& operator+=(Vector& a, const Vector& b);
 | 
|---|
 | 57 | Vector& operator*=(Vector& a, const double m);
 | 
|---|
 | 58 | Vector& operator*(const Vector& a, const double m);
 | 
|---|
 | 59 | Vector& operator+(const Vector& a, const Vector& b);
 | 
|---|
| [342f33f] | 60 | 
 | 
|---|
 | 61 | #endif /*VECTOR_HPP_*/
 | 
|---|