source: src/Descriptors/MoleculeDescriptor.cpp@ e215c1

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 e215c1 was 02ce36, checked in by Frederik Heber <heber@…>, 13 years ago

We now require CodePatterns 1.2.0.

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