| 1 | #ifndef VECTOR_HPP_
 | 
|---|
| 2 | #define VECTOR_HPP_
 | 
|---|
| 3 | 
 | 
|---|
| 4 | class Vector;
 | 
|---|
| 5 | 
 | 
|---|
| 6 | /** Single vector.
 | 
|---|
| 7 |  * basically, just a x[3] but with helpful functions
 | 
|---|
| 8 |  */
 | 
|---|
| 9 | class Vector {   
 | 
|---|
| 10 |   public:
 | 
|---|
| 11 |     double x[NDIM];
 | 
|---|
| 12 | 
 | 
|---|
| 13 |   Vector();
 | 
|---|
| 14 |   Vector(double x1, double x2, double x3);
 | 
|---|
| 15 |   ~Vector();
 | 
|---|
| 16 | 
 | 
|---|
| 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;
 | 
|---|
| 21 |   double Norm() const ;
 | 
|---|
| 22 |   double Angle(Vector *y) const;
 | 
|---|
| 23 | 
 | 
|---|
| 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);
 | 
|---|
| 29 |   void Zero(); 
 | 
|---|
| 30 |   void One(double one);
 | 
|---|
| 31 |   void Init(double x1, double x2, double x3);
 | 
|---|
| 32 |   void Normalize();
 | 
|---|
| 33 |   void Translate(const Vector *x);
 | 
|---|
| 34 |   void Mirror(const Vector *x);
 | 
|---|
| 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);
 | 
|---|
| 41 |   void LinearCombinationOfVectors(const Vector *x1, const Vector *x2, const Vector *x3, double *factors);
 | 
|---|
| 42 |   
 | 
|---|
| 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);
 | 
|---|
| 50 | 
 | 
|---|
| 51 |   void AskPosition(double *cell_size, bool check);
 | 
|---|
| 52 |   bool Output(ofstream *out) const;
 | 
|---|
| 53 | };
 | 
|---|
| 54 | 
 | 
|---|
| 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);
 | 
|---|
| 60 | 
 | 
|---|
| 61 | #endif /*VECTOR_HPP_*/
 | 
|---|