source: src/vector.hpp@ 7ac4af

Action_Thermostats Add_AtomRandomPerturbation Add_FitFragmentPartialChargesAction Add_RotateAroundBondAction Add_SelectAtomByNameAction Added_ParseSaveFragmentResults AddingActions_SaveParseParticleParameters Adding_Graph_to_ChangeBondActions Adding_MD_integration_tests Adding_ParticleName_to_Atom Adding_StructOpt_integration_tests AtomFragments Automaking_mpqc_open AutomationFragmentation_failures Candidate_v1.5.4 Candidate_v1.6.0 Candidate_v1.6.1 Candidate_v1.7.0 ChangeBugEmailaddress ChangingTestPorts ChemicalSpaceEvaluator CombiningParticlePotentialParsing Combining_Subpackages Debian_Package_split Debian_package_split_molecuildergui_only Disabling_MemDebug Docu_Python_wait EmpiricalPotential_contain_HomologyGraph EmpiricalPotential_contain_HomologyGraph_documentation Enable_parallel_make_install Enhance_userguide Enhanced_StructuralOptimization Enhanced_StructuralOptimization_continued Example_ManyWaysToTranslateAtom Exclude_Hydrogens_annealWithBondGraph FitPartialCharges_GlobalError Fix_BoundInBox_CenterInBox_MoleculeActions Fix_ChargeSampling_PBC Fix_ChronosMutex Fix_FitPartialCharges Fix_FitPotential_needs_atomicnumbers Fix_ForceAnnealing Fix_IndependentFragmentGrids Fix_ParseParticles Fix_ParseParticles_split_forward_backward_Actions Fix_PopActions Fix_QtFragmentList_sorted_selection Fix_Restrictedkeyset_FragmentMolecule Fix_StatusMsg Fix_StepWorldTime_single_argument Fix_Verbose_Codepatterns Fix_fitting_potentials Fixes ForceAnnealing_goodresults ForceAnnealing_oldresults ForceAnnealing_tocheck ForceAnnealing_with_BondGraph ForceAnnealing_with_BondGraph_continued ForceAnnealing_with_BondGraph_continued_betteresults ForceAnnealing_with_BondGraph_contraction-expansion FragmentAction_writes_AtomFragments FragmentMolecule_checks_bonddegrees GeometryObjects Gui_Fixes Gui_displays_atomic_force_velocity ImplicitCharges IndependentFragmentGrids IndependentFragmentGrids_IndividualZeroInstances IndependentFragmentGrids_IntegrationTest IndependentFragmentGrids_Sole_NN_Calculation JobMarket_RobustOnKillsSegFaults JobMarket_StableWorkerPool JobMarket_unresolvable_hostname_fix MoreRobust_FragmentAutomation ODR_violation_mpqc_open PartialCharges_OrthogonalSummation PdbParser_setsAtomName PythonUI_with_named_parameters QtGui_reactivate_TimeChanged_changes Recreated_GuiChecks Rewrite_FitPartialCharges RotateToPrincipalAxisSystem_UndoRedo SaturateAtoms_findBestMatching SaturateAtoms_singleDegree StoppableMakroAction Subpackage_CodePatterns Subpackage_JobMarket Subpackage_LinearAlgebra Subpackage_levmar Subpackage_mpqc_open Subpackage_vmg Switchable_LogView ThirdParty_MPQC_rebuilt_buildsystem TrajectoryDependenant_MaxOrder TremoloParser_IncreasedPrecision TremoloParser_MultipleTimesteps TremoloParser_setsAtomName Ubuntu_1604_changes stable
Last change on this file since 7ac4af was 4b94bb, checked in by Tillmann Crueger <crueger@…>, 15 years ago

Moved Vector-Matrix operations to matrix class
  • Property mode set to 100755
File size: 3.2 KB
Line 
1#ifndef VECTOR_HPP_
2#define VECTOR_HPP_
3
4using namespace std;
5
6/*********************************************** includes ***********************************/
7
8// include config.h
9#ifdef HAVE_CONFIG_H
10#include <config.h>
11#endif
12
13#include <iostream>
14#include <gsl/gsl_vector.h>
15
16#include <memory>
17#include <vector>
18
19#include "defs.hpp"
20#include "Space.hpp"
21
22/********************************************** declarations *******************************/
23
24class Vector;
25class Matrix;
26
27typedef std::vector<Vector> pointset;
28
29/** Single vector.
30 * basically, just a x[3] but with helpful functions
31 */
32class Vector : public Space{
33 friend Vector operator*(const Matrix&,const Vector&);
34public:
35 Vector();
36 Vector(const double x1, const double x2, const double x3);
37 Vector(const Vector& src);
38 virtual ~Vector();
39
40 // Method implemented by forwarding to the Representation
41
42 double DistanceSquared(const Vector &y) const;
43 double DistanceToSpace(const Space& space) const;
44 double ScalarProduct(const Vector &y) const;
45 double Angle(const Vector &y) const;
46 bool IsZero() const;
47 bool IsOne() const;
48 bool IsNormalTo(const Vector &normal) const;
49 bool IsEqualTo(const Vector &a) const;
50
51 void AddVector(const Vector &y);
52 void SubtractVector(const Vector &y);
53 void VectorProduct(const Vector &y);
54 void ProjectOntoPlane(const Vector &y);
55 void ProjectIt(const Vector &y);
56 Vector Projection(const Vector &y) const;
57 void ScaleAll(const double *factor);
58 void ScaleAll(const Vector &factor);
59 void Scale(const double factor);
60 bool GetOneNormalVector(const Vector &x1);
61 bool MakeNormalTo(const Vector &y1);
62 std::pair<Vector,Vector> partition(const Vector&) const;
63 std::pair<pointset,Vector> partition(const pointset&) const;
64
65 // Accessors ussually come in pairs... and sometimes even more than that
66 double& operator[](size_t i);
67 const double& operator[](size_t i) const;
68 double& at(size_t i);
69 const double& at(size_t i) const;
70
71 // Assignment operator
72 Vector &operator=(const Vector& src);
73
74 // Access to internal structure
75 gsl_vector* get();
76
77 // Methods that are derived directly from other methods
78 double Norm() const;
79 double NormSquared() const;
80 void Normalize();
81 void Zero();
82 void One(const double one);
83 void LinearCombinationOfVectors(const Vector &x1, const Vector &x2, const Vector &x3, const double * const factors);
84
85 // operators for mathematical operations
86 bool operator==(const Vector& b) const;
87 bool operator!=(const Vector& b) const;
88 const Vector& operator+=(const Vector& b);
89 const Vector& operator-=(const Vector& b);
90 Vector const operator+(const Vector& b) const;
91 Vector const operator-(const Vector& b) const;
92
93 // Methods inherited from Space
94 virtual double distance(const Vector &point) const;
95 virtual Vector getClosestPoint(const Vector &point) const;
96
97protected:
98
99private:
100 Vector(gsl_vector *);
101 gsl_vector *content;
102
103};
104
105// some commonly used and fixed vectors
106const extern Vector zeroVec;
107const extern Vector e1;
108const extern Vector e2;
109const extern Vector e3;
110
111ostream & operator << (ostream& ost, const Vector &m);
112const Vector& operator*=(Vector& a, const double m);
113Vector const operator*(const Vector& a, const double m);
114Vector const operator*(const double m, const Vector& a);
115
116#endif /*VECTOR_HPP_*/
Note: See TracBrowser for help on using the repository browser.