source: src/LinearAlgebra/Vector.hpp@ 111387

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 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 111387 was 3da2fb, checked in by Frederik Heber <heber@…>, 14 years ago

New module MatrixVector_ops that contains matrix vector products.

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