source: src/Descriptors/MoleculeDescriptor.cpp@ b8d15ba

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 b8d15ba was 8cce2b, checked in by Tillmann Crueger <crueger@…>, 14 years ago

Made the MoleculeDescriptor and AtomDescriptor use faster STL-Algorithms for finding

  • Property mode set to 100644
File size: 4.8 KB
Line 
1/*
2 * MoleculeDescriptor.cpp
3 *
4 * Created on: Feb 5, 2010
5 * Author: crueger
6 */
7
8#include "Helpers/MemDebug.hpp"
9
10#include "Descriptors/MoleculeDescriptor.hpp"
11#include "Descriptors/MoleculeDescriptor_impl.hpp"
12
13#include "World.hpp"
14#include "Patterns/ObservedContainer_impl.hpp"
15
16#include "molecule.hpp"
17
18#include <boost/bind.hpp>
19#include <iostream>
20
21using namespace std;
22
23typedef World::MoleculeSet::internal_iterator molecules_iter_t;
24
25/************************ Forwarding object **************************************/
26
27
28MoleculeDescriptor::MoleculeDescriptor(impl_ptr _impl) :
29 impl(_impl)
30{}
31
32MoleculeDescriptor::MoleculeDescriptor(const MoleculeDescriptor& src) :
33 impl(src.get_impl())
34{}
35
36MoleculeDescriptor::~MoleculeDescriptor()
37{}
38
39MoleculeDescriptor& MoleculeDescriptor::operator=(MoleculeDescriptor &src){
40 if(&src!=this) {
41 impl=src.get_impl();
42 }
43 return *this;
44}
45
46molecule* MoleculeDescriptor::find(){
47 return impl->find();
48}
49
50std::vector<molecule*> MoleculeDescriptor::findAll(){
51 return impl->findAll();
52}
53
54MoleculeDescriptor::impl_ptr MoleculeDescriptor::get_impl() const{
55 return impl;
56}
57
58
59
60
61/**************************** implementation ********************/
62
63MoleculeDescriptor_impl::MoleculeDescriptor_impl()
64{
65}
66
67MoleculeDescriptor_impl::~MoleculeDescriptor_impl()
68{
69}
70
71World::MoleculeSet& MoleculeDescriptor_impl::getMolecules(){
72 return World::getInstance().molecules;
73}
74
75molecule* MoleculeDescriptor_impl::find() {
76 World::MoleculeSet &molecules = getMolecules();
77 molecules_iter_t res = find_if(molecules.begin_internal(),molecules.end_internal(),boost::bind(&MoleculeDescriptor_impl::predicate,this,_1));
78 return (res!=molecules.end_internal())?((*res).second):0;
79}
80
81vector<molecule*> MoleculeDescriptor_impl::findAll() {
82 vector<molecule*> res;
83 World::MoleculeSet &molecules = getMolecules();
84 for_each(molecules.begin_internal(),
85 molecules.end_internal(),
86 boost::bind(&MoleculeDescriptor_impl::checkAndAdd,
87 this,&res,_1));
88 return res;
89}
90
91void MoleculeDescriptor_impl::checkAndAdd(std::vector<molecule*> *v,std::pair<moleculeId_t,molecule*> p){
92 if(predicate(p)){
93 v->push_back(p.second);
94 }
95}
96
97/************************** Universe and Emptyset *****************/
98
99MoleculeAllDescriptor_impl::MoleculeAllDescriptor_impl()
100{}
101
102MoleculeAllDescriptor_impl::~MoleculeAllDescriptor_impl()
103{}
104
105bool MoleculeAllDescriptor_impl::predicate(std::pair<moleculeId_t,molecule*>){
106 return true;
107}
108
109MoleculeDescriptor AllMolecules(){
110 return MoleculeDescriptor(MoleculeDescriptor::impl_ptr(new MoleculeAllDescriptor_impl));
111}
112
113MoleculeNoneDescriptor_impl::MoleculeNoneDescriptor_impl()
114{}
115
116MoleculeNoneDescriptor_impl::~MoleculeNoneDescriptor_impl()
117{}
118
119bool MoleculeNoneDescriptor_impl::predicate(std::pair<moleculeId_t,molecule*>){
120 return false;
121}
122
123MoleculeDescriptor NoMolecules(){
124 return MoleculeDescriptor(MoleculeDescriptor::impl_ptr(new MoleculeNoneDescriptor_impl));
125}
126
127/************************** Operator stuff ************************/
128
129// AND
130MoleculeAndDescriptor_impl::MoleculeAndDescriptor_impl(MoleculeDescriptor::impl_ptr _lhs, MoleculeDescriptor::impl_ptr _rhs) :
131 lhs(_lhs), rhs(_rhs)
132{}
133
134MoleculeAndDescriptor_impl::~MoleculeAndDescriptor_impl()
135{}
136
137bool MoleculeAndDescriptor_impl::predicate(std::pair<moleculeId_t,molecule*> molecule){
138 return lhs->predicate(molecule) && rhs->predicate(molecule);
139}
140MoleculeDescriptor operator&&(const MoleculeDescriptor &lhs, const MoleculeDescriptor &rhs){
141 MoleculeDescriptor::impl_ptr newImpl = MoleculeDescriptor::impl_ptr(new MoleculeAndDescriptor_impl(lhs.get_impl(),rhs.get_impl()));
142 return MoleculeDescriptor(newImpl);
143}
144
145// OR
146MoleculeOrDescriptor_impl::MoleculeOrDescriptor_impl(MoleculeDescriptor::impl_ptr _lhs ,MoleculeDescriptor::impl_ptr _rhs) :
147 lhs(_lhs), rhs(_rhs)
148{}
149
150MoleculeOrDescriptor_impl::~MoleculeOrDescriptor_impl(){
151}
152
153bool MoleculeOrDescriptor_impl::predicate(std::pair<moleculeId_t,molecule*> molecule){
154 return lhs->predicate(molecule) || rhs->predicate(molecule);
155}
156
157MoleculeDescriptor operator||(const MoleculeDescriptor &lhs, const MoleculeDescriptor &rhs){
158 MoleculeDescriptor::impl_ptr newImpl = MoleculeDescriptor::impl_ptr(new MoleculeOrDescriptor_impl(lhs.get_impl(),rhs.get_impl()));
159 return MoleculeDescriptor(newImpl);
160}
161
162// NOT
163
164MoleculeNotDescriptor_impl::MoleculeNotDescriptor_impl(MoleculeDescriptor::impl_ptr _arg) :
165 arg(_arg)
166{}
167
168
169MoleculeNotDescriptor_impl::~MoleculeNotDescriptor_impl()
170{
171}
172
173bool MoleculeNotDescriptor_impl::predicate(std::pair<moleculeId_t,molecule*> molecule){
174 return !(arg->predicate(molecule));
175}
176
177MoleculeDescriptor operator!(const MoleculeDescriptor &arg){
178 MoleculeDescriptor::impl_ptr newImpl = MoleculeDescriptor::impl_ptr(new MoleculeNotDescriptor_impl(arg.get_impl()));
179 return MoleculeDescriptor(newImpl);
180}
Note: See TracBrowser for help on using the repository browser.