[342f33f] | 1 | #ifndef VECTOR_HPP_
|
---|
| 2 | #define VECTOR_HPP_
|
---|
| 3 |
|
---|
[357fba] | 4 | using namespace std;
|
---|
| 5 |
|
---|
[f66195] | 6 | /*********************************************** includes ***********************************/
|
---|
| 7 |
|
---|
| 8 | // include config.h
|
---|
| 9 | #ifdef HAVE_CONFIG_H
|
---|
| 10 | #include <config.h>
|
---|
| 11 | #endif
|
---|
[357fba] | 12 |
|
---|
[c75de1] | 13 | #include <iostream>
|
---|
[54a746] | 14 | #include <gsl/gsl_vector.h>
|
---|
| 15 | #include <gsl/gsl_multimin.h>
|
---|
| 16 |
|
---|
[f66195] | 17 | #include "defs.hpp"
|
---|
| 18 |
|
---|
| 19 | /********************************************** declarations *******************************/
|
---|
[7d96c6] | 20 |
|
---|
[342f33f] | 21 | /** Single vector.
|
---|
| 22 | * basically, just a x[3] but with helpful functions
|
---|
| 23 | */
|
---|
[e4ea46] | 24 | class Vector {
|
---|
[042f82] | 25 | public:
|
---|
| 26 | double x[NDIM];
|
---|
[6ac7ee] | 27 |
|
---|
[042f82] | 28 | Vector();
|
---|
[fb73b8] | 29 | Vector(const double x1, const double x2, const double x3);
|
---|
[042f82] | 30 | ~Vector();
|
---|
[6ac7ee] | 31 |
|
---|
[776b64] | 32 | double Distance(const Vector * const y) const;
|
---|
| 33 | double DistanceSquared(const Vector * const y) const;
|
---|
[e138de] | 34 | double DistanceToPlane(const Vector * const PlaneNormal, const Vector * const PlaneOffset) const;
|
---|
[776b64] | 35 | double PeriodicDistance(const Vector * const y, const double * const cell_size) const;
|
---|
| 36 | double PeriodicDistanceSquared(const Vector * const y, const double * const cell_size) const;
|
---|
| 37 | double ScalarProduct(const Vector * const y) const;
|
---|
[d4d0dd] | 38 | double Norm() const;
|
---|
| 39 | double NormSquared() const;
|
---|
[776b64] | 40 | double Angle(const Vector * const y) const;
|
---|
[54a746] | 41 | bool IsZero() const;
|
---|
| 42 | bool IsOne() const;
|
---|
[776b64] | 43 | bool IsNormalTo(const Vector * const normal) const;
|
---|
[b998c3] | 44 | bool IsEqualTo(const Vector * const a) const;
|
---|
[6ac7ee] | 45 |
|
---|
[776b64] | 46 | void AddVector(const Vector * const y);
|
---|
| 47 | void SubtractVector(const Vector * const y);
|
---|
| 48 | void CopyVector(const Vector * const y);
|
---|
| 49 | void CopyVector(const Vector &y);
|
---|
| 50 | void RotateVector(const Vector * const y, const double alpha);
|
---|
| 51 | void VectorProduct(const Vector * const y);
|
---|
| 52 | void ProjectOntoPlane(const Vector * const y);
|
---|
| 53 | void ProjectIt(const Vector * const y);
|
---|
| 54 | Vector Projection(const Vector * const y) const;
|
---|
[042f82] | 55 | void Zero();
|
---|
[776b64] | 56 | void One(const double one);
|
---|
| 57 | void Init(const double x1, const double x2, const double x3);
|
---|
[042f82] | 58 | void Normalize();
|
---|
[776b64] | 59 | void Translate(const Vector * const x);
|
---|
| 60 | void Mirror(const Vector * const x);
|
---|
| 61 | void Scale(const double ** const factor);
|
---|
| 62 | void Scale(const double * const factor);
|
---|
| 63 | void Scale(const double factor);
|
---|
| 64 | void MatrixMultiplication(const double * const M);
|
---|
| 65 | void InverseMatrixMultiplication(const double * const M);
|
---|
[e138de] | 66 | void KeepPeriodic(const double * const matrix);
|
---|
[776b64] | 67 | void LinearCombinationOfVectors(const Vector * const x1, const Vector * const x2, const Vector * const x3, const double * const factors);
|
---|
| 68 | double CutsPlaneAt(const Vector * const A, const Vector * const B, const Vector * const C) const;
|
---|
[e138de] | 69 | bool GetIntersectionWithPlane(const Vector * const PlaneNormal, const Vector * const PlaneOffset, const Vector * const Origin, const Vector * const LineVector);
|
---|
| 70 | bool GetIntersectionOfTwoLinesOnPlane(const Vector * const Line1a, const Vector * const Line1b, const Vector * const Line2a, const Vector * const Line2b, const Vector *Normal = NULL);
|
---|
[776b64] | 71 | bool GetOneNormalVector(const Vector * const x1);
|
---|
| 72 | bool MakeNormalVector(const Vector * const y1);
|
---|
| 73 | bool MakeNormalVector(const Vector * const y1, const Vector * const y2);
|
---|
| 74 | bool MakeNormalVector(const Vector * const x1, const Vector * const x2, const Vector * const x3);
|
---|
| 75 | bool SolveSystem(Vector * x1, Vector * x2, Vector * y, const double alpha, const double beta, const double c);
|
---|
| 76 | bool LSQdistance(const Vector ** vectors, int dim);
|
---|
| 77 | void AskPosition(const double * const cell_size, const bool check);
|
---|
[e138de] | 78 | void Output() const;
|
---|
[776b64] | 79 | bool IsInParallelepiped(const Vector &offset, const double * const parallelepiped) const;
|
---|
| 80 | void WrapPeriodically(const double * const M, const double * const Minv);
|
---|
[2ededc2] | 81 |
|
---|
| 82 | Vector& operator=(const Vector &src);
|
---|
| 83 |
|
---|
[342f33f] | 84 | };
|
---|
| 85 |
|
---|
[9c20aa] | 86 | ostream & operator << (ostream& ost, const Vector &m);
|
---|
[ef9df36] | 87 | bool operator==(const Vector& a, const Vector& b);
|
---|
[54a746] | 88 | Vector& operator+=(Vector& a, const Vector& b);
|
---|
| 89 | Vector& operator-=(Vector& a, const Vector& b);
|
---|
| 90 | Vector& operator*=(Vector& a, const double m);
|
---|
| 91 | Vector& operator*(const Vector& a, const double m);
|
---|
| 92 | Vector& operator*(const double m, const Vector& a);
|
---|
| 93 | Vector& operator+(const Vector& a, const Vector& b);
|
---|
| 94 | Vector& operator-(const Vector& a, const Vector& b);
|
---|
| 95 |
|
---|
[342f33f] | 96 |
|
---|
| 97 | #endif /*VECTOR_HPP_*/
|
---|