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