source: src/LinearAlgebra/Eigenspace.hpp@ 3f9eba

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 3f9eba was 56f73b, checked in by Frederik Heber <heber@…>, 15 years ago

Added config.h also to all header files, code check test ascertain this in the future.

  • as we want to use config.h to pass stuff such as MEMDEBUG, NDEBUG, LOG_OBSERVER, we have to make sure that it is present in each and every file.
  • split up CodeChecks/testsuite.at: each test has its own .at file.
  • Property mode set to 100644
File size: 2.0 KB
Line 
1/*
2 * Eigenspace.hpp
3 *
4 * Created on: Nov 22, 2010
5 * Author: heber
6 */
7
8#ifndef EIGENSPACE_HPP_
9#define EIGENSPACE_HPP_
10
11// include config.h
12#ifdef HAVE_CONFIG_H
13#include <config.h>
14#endif
15
16
17#include <boost/shared_ptr.hpp>
18
19#include <iosfwd>
20#include <set>
21#include <vector>
22#include "MatrixContent.hpp"
23#include "VectorContent.hpp"
24
25/** Container for a matrix and its eigenvectors and -values.
26 *
27 * Note that a set of eigenvectors always forms a subspace within a vector
28 * space, the so-called eigenspace. Hence, the name of this class is
29 * eigenspace.
30 *
31 * Eigenvector are stored internally as a matrix. A returned set of
32 * eigenvectors is just a set of VectorViewContent on row vectors in
33 * this matrix.
34 */
35class Eigenspace
36{
37 friend std::ostream & operator<<(std::ostream &ost, const Eigenspace &_s);
38public:
39 typedef std::vector< boost::shared_ptr<VectorContent> > eigenvectorset;
40 typedef std::vector< double > eigenvalueset;
41 typedef std::set<size_t> indexset;
42
43 Eigenspace(const indexset &_Indices);
44 Eigenspace(const indexset &_Indices, const MatrixContent &_matrix);
45 virtual ~Eigenspace();
46
47 void calculateEigenSubspace();
48
49 // accessors
50 const MatrixContent & getEigenspaceMatrix() const;
51 const VectorContent & getEigenvector(const size_t i) const;
52 const MatrixContent & getEigenvectorMatrix() const;
53 const double getEigenvalue(const size_t i) const;
54 const eigenvalueset & getEigenvalues() const;
55 const indexset & getIndices() const;
56 void setEigenspaceMatrix(const MatrixContent &_content);
57
58 void setEigenpairs(const eigenvectorset &CurrentEigenvectors, const eigenvalueset &CurrentEigenvalues);
59
60 bool contains(const Eigenspace &_indexset);
61
62protected:
63 const indexset Indices;
64 eigenvectorset Eigenvectors;
65 eigenvalueset Eigenvalues;
66 MatrixContent EigenspaceMatrix;
67 MatrixContent EigenvectorMatrix;
68
69private:
70 void createEigenvectorViews();
71
72};
73
74std::ostream & operator<<(std::ostream &ost, const Eigenspace &_s);
75
76#endif /* EIGENSPACE_HPP_ */
Note: See TracBrowser for help on using the repository browser.