source: src/LinearAlgebra/Vector.hpp@ b4cf2b

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 b4cf2b was bbf1bd, checked in by Frederik Heber <heber@…>, 15 years ago

Extended VectorContent class.

  • VectorContent formerly has been just a structure to contain the gsl_vector due to forward declaration reasons.
  • now VectorContent is a true wrapper to gsl_vector, i.e. all functionality that is now specific to 3 dimensions has been shifted from GSLVector over to VectorContent.
  • VectorContentView is preserved to allow for VectorContent as a view on a row or column of a matrix.
  • changed and renamed unit test gslvectorunittest -> VectorContentUnitTest
  • GSLVector is not used anymore anywhere
  • one long-sough error was a missing assignment operator filled-in in a wrong manner by the compiler for VectorContent.

Note that:

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