| [342f33f] | 1 | #ifndef VECTOR_HPP_ | 
|---|
|  | 2 | #define VECTOR_HPP_ | 
|---|
|  | 3 |  | 
|---|
| [7d96c6] | 4 | class vector; | 
|---|
|  | 5 |  | 
|---|
| [342f33f] | 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(); | 
|---|
| [498a9f] | 14 | vector(double x1, double x2, double x3); | 
|---|
| [342f33f] | 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); | 
|---|
| [498a9f] | 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(); | 
|---|
|  | 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 |  | 
|---|
| [110ceb] | 43 | double CutsPlaneAt(vector *A, vector *B, vector *C); | 
|---|
| [342f33f] | 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_*/ | 
|---|