source: src/Matrix.hpp@ 51be2a

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 51be2a was 0eb2dc, checked in by Tillmann Crueger <crueger@…>, 14 years ago

Added comparison Operator to matrix class

  • Property mode set to 100644
File size: 2.9 KB
Line 
1/*
2 * Matrix.hpp
3 *
4 * Created on: Jun 25, 2010
5 * Author: crueger
6 */
7
8#ifndef MATRIX_HPP_
9#define MATRIX_HPP_
10
11#include <iosfwd>
12#include "defs.hpp"
13
14/**
15 * Simple class to store matrices with a few extra functions
16 */
17
18class Vector;
19class MatrixContent;
20
21class Matrix
22{
23 friend Vector operator*(const Matrix&,const Vector&);
24public:
25 Matrix();
26
27 /**
28 * construct a matrix from a 3x3 double array that contains all Elements.
29 *
30 * The elements are laid out in the following way:
31 * array -> matrix
32 * 0 -> (0,0)
33 * 1 -> (1,0)
34 * 2 -> (2,0)
35 * 3 -> (0,1)
36 * 4 -> (1,1)
37 * 5 -> (2,1)
38 * 6 -> (0,2)
39 * 7 -> (1,2)
40 * 8 -> (2,2)
41 *
42 */
43 Matrix(const double*);
44 Matrix(const Matrix&);
45 virtual ~Matrix();
46
47 /**
48 * Set the matrix to a unit matrix.
49 */
50 void one();
51
52 /**
53 * Access the matrix at index (i,j)
54 */
55 double &at(size_t i, size_t j);
56 /**
57 * Access the matrix at index (i,j)
58 */
59 const double at(size_t i, size_t j) const;
60
61 /**
62 * Set the matrix at index (i,j).
63 *
64 * Slightly faster than at(i,j)=x
65 */
66 void set(size_t i, size_t j, const double value);
67
68 /**
69 * get the ith row of the matrix as a vector
70 */
71 Vector &row(size_t);
72 const Vector &row(size_t) const;
73
74 /**
75 * get the ith column of the matrix as a vector
76 */
77 Vector &column(size_t);
78 const Vector &column(size_t) const;
79
80 /**
81 * get the diagonal of the matrix as a vector
82 */
83 Vector &diagonal();
84 const Vector &diagonal() const;
85
86 /**
87 * Calculate the determinant of the matrix
88 */
89 double determinant() const;
90
91 /**
92 * Calculate the inverse of the matrix.
93 *
94 * Rather costly, so use precomputation as often as possible.
95 */
96 Matrix invert() const;
97
98 // operators
99 Matrix &operator=(const Matrix&);
100
101 Matrix &operator+=(const Matrix&);
102 Matrix &operator-=(const Matrix&);
103 Matrix &operator*=(const Matrix&);
104
105 Matrix &operator*=(const double);
106
107 Matrix operator+(const Matrix&) const;
108 Matrix operator-(const Matrix&) const;
109 Matrix operator*(const Matrix&) const;
110
111 bool operator==(const Matrix&) const;
112
113private:
114 Matrix(MatrixContent*);
115 void createViews();
116 MatrixContent *content;
117 // we keep around some Vector views of the matrix, to return references
118 Vector* rows_ptr[NDIM];
119 Vector* columns_ptr[NDIM];
120 Vector* diagonal_ptr;
121};
122
123Matrix operator*(const double,const Matrix&);
124Matrix operator*(const Matrix&,const double);
125
126/**
127 * Takes a symmetric matrix that stores the lower diagonal and produces a
128 * full matrix.
129 *
130 * The array is laid out as follows:
131 *
132 * array -> matrix
133 * 0 -> (0,0)
134 * 1 -> (1,0);(0,1)
135 * 2 -> (1,1)
136 * 3 -> (2,0);(0,2)
137 * 4 -> (2,1);(1,2)
138 * 5 -> (2,2)
139 */
140Matrix ReturnFullMatrixforSymmetric(const double * const cell_size);
141
142std::ostream &operator<<(std::ostream&,const Matrix&);
143Vector operator*(const Matrix&,const Vector&);
144Vector& operator*=(Vector&,const Matrix&);
145
146#endif /* MATRIX_HPP_ */
Note: See TracBrowser for help on using the repository browser.