Changeset 71910a for molecuilder/src/vector.hpp

Timestamp:

Apr 7, 2010, 3:45:38 PM (16 years ago)

Author:

Tillmann Crueger <crueger@…>

Children:

0f55b2

Parents:

770138

Message:

Made data internal data-structure of vector class private

• Replaced occurences of access to internals with operator
• moved Vector-class into LinAlg-Module
• Reworked Vector to allow clean modularization
• Added Plane class to describe arbitrary planes in 3d space

File:

• molecuilder/src/vector.hpp (modified) (3 diffs)

molecuilder/src/vector.hpp

@@ -24 +24 @@
 class Vector {
   public:
-    double x[NDIM];
 
   Vector();
   Vector(const double x1, const double x2, const double x3);
+  Vector(const Vector& src);
   ~Vector();
 
@@ -48 +48 @@
   void CopyVector(const Vector * const y);
   void CopyVector(const Vector &y);
-  void RotateVector(const Vector * const y, const double alpha);
   void VectorProduct(const Vector * const y);
   void ProjectOntoPlane(const Vector * const y);
@@ -63 +62 @@
   void Scale(const double factor);
   void MatrixMultiplication(const double * const M);
-  void InverseMatrixMultiplication(const double * const M);
+  bool InverseMatrixMultiplication(const double * const M);
   void KeepPeriodic(const double * const matrix);
   void LinearCombinationOfVectors(const Vector * const x1, const Vector * const x2, const Vector * const x3, const double * const factors);
   double CutsPlaneAt(const Vector * const A, const Vector * const B, const Vector * const C) const;
-  bool GetIntersectionWithPlane(const Vector * const PlaneNormal, const Vector * const PlaneOffset, const Vector * const Origin, const Vector * const LineVector);
-  bool GetIntersectionOfTwoLinesOnPlane(const Vector * const Line1a, const Vector * const Line1b, const Vector * const Line2a, const Vector * const Line2b, const Vector *Normal = NULL);
   bool GetOneNormalVector(const Vector * const x1);
-  bool MakeNormalVector(const Vector * const y1);
-  bool MakeNormalVector(const Vector * const y1, const Vector * const y2);
-  bool MakeNormalVector(const Vector * const x1, const Vector * const x2, const Vector * const x3);
-  bool SolveSystem(Vector * x1, Vector * x2, Vector * y, const double alpha, const double beta, const double c);
-  bool LSQdistance(const Vector ** vectors, int dim);
-  void AskPosition(const double * const cell_size, const bool check);
-  void Output() const;
+  bool MakeNormalTo(const Vector &y1);
+  //bool SolveSystem(Vector * x1, Vector * x2, Vector * y, const double alpha, const double beta, const double c);
   bool IsInParallelepiped(const Vector &offset, const double * const parallelepiped) const;
   void WrapPeriodically(const double * const M, const double * const Minv);
+
+  // Accessors ussually come in pairs... and sometimes even more than that
+  double& operator[](size_t i);
+  const double& operator[](size_t i) const;
+  double& at(size_t i);
+  const double& at(size_t i) const;
+
+  // Assignment operator
+  Vector &operator=(const Vector& src);
+
+  // Access to internal structure
+  double* get();
+
+private:
+  double x[NDIM];
 
 };