| 1 | /*
 | 
|---|
| 2 |  * Project: MoleCuilder
 | 
|---|
| 3 |  * Description: creates and alters molecular systems
 | 
|---|
| 4 |  * Copyright (C)  2010 University of Bonn. All rights reserved.
 | 
|---|
| 5 |  * Please see the LICENSE file or "Copyright notice" in builder.cpp for details.
 | 
|---|
| 6 |  */
 | 
|---|
| 7 | 
 | 
|---|
| 8 | /*
 | 
|---|
| 9 |  * AtomDescriptor.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/AtomDescriptor.hpp"
 | 
|---|
| 23 | #include "Descriptors/AtomDescriptor_impl.hpp"
 | 
|---|
| 24 | 
 | 
|---|
| 25 | #include "World.hpp"
 | 
|---|
| 26 | #include "atom.hpp"
 | 
|---|
| 27 | #include "CodePatterns/ObservedContainer_impl.hpp"
 | 
|---|
| 28 | 
 | 
|---|
| 29 | #include <boost/bind.hpp>
 | 
|---|
| 30 | 
 | 
|---|
| 31 | #include <iostream>
 | 
|---|
| 32 | 
 | 
|---|
| 33 | using namespace std;
 | 
|---|
| 34 | 
 | 
|---|
| 35 | typedef World::AtomSet::internal_iterator atoms_iter_t;
 | 
|---|
| 36 | 
 | 
|---|
| 37 | /************************ Forwarding object **************************************/
 | 
|---|
| 38 | 
 | 
|---|
| 39 | 
 | 
|---|
| 40 | AtomDescriptor::AtomDescriptor(impl_ptr _impl) :
 | 
|---|
| 41 |     impl(_impl)
 | 
|---|
| 42 | {}
 | 
|---|
| 43 | 
 | 
|---|
| 44 | AtomDescriptor::AtomDescriptor(const AtomDescriptor& src) :
 | 
|---|
| 45 |     impl(src.get_impl())
 | 
|---|
| 46 | {}
 | 
|---|
| 47 | 
 | 
|---|
| 48 | AtomDescriptor::~AtomDescriptor()
 | 
|---|
| 49 | {}
 | 
|---|
| 50 | 
 | 
|---|
| 51 | AtomDescriptor& AtomDescriptor::operator=(AtomDescriptor &src){
 | 
|---|
| 52 |   if(&src!=this) {
 | 
|---|
| 53 |     impl=src.get_impl();
 | 
|---|
| 54 |   }
 | 
|---|
| 55 |   return *this;
 | 
|---|
| 56 | }
 | 
|---|
| 57 | 
 | 
|---|
| 58 | atom* AtomDescriptor::find(){
 | 
|---|
| 59 |   return impl->find();
 | 
|---|
| 60 | }
 | 
|---|
| 61 | 
 | 
|---|
| 62 | std::vector<atom*> AtomDescriptor::findAll(){
 | 
|---|
| 63 |   return impl->findAll();
 | 
|---|
| 64 | }
 | 
|---|
| 65 | 
 | 
|---|
| 66 | AtomDescriptor::impl_ptr AtomDescriptor::get_impl() const{
 | 
|---|
| 67 |   return impl;
 | 
|---|
| 68 | }
 | 
|---|
| 69 | 
 | 
|---|
| 70 | 
 | 
|---|
| 71 | 
 | 
|---|
| 72 | 
 | 
|---|
| 73 | /**************************** implementation ********************/
 | 
|---|
| 74 | 
 | 
|---|
| 75 | AtomDescriptor_impl::AtomDescriptor_impl()
 | 
|---|
| 76 | {
 | 
|---|
| 77 | }
 | 
|---|
| 78 | 
 | 
|---|
| 79 | AtomDescriptor_impl::~AtomDescriptor_impl()
 | 
|---|
| 80 | {
 | 
|---|
| 81 | }
 | 
|---|
| 82 | 
 | 
|---|
| 83 | World::AtomSet& AtomDescriptor_impl::getAtoms(){
 | 
|---|
| 84 |   return World::getInstance().atoms;
 | 
|---|
| 85 | }
 | 
|---|
| 86 | 
 | 
|---|
| 87 | atom* AtomDescriptor_impl::find() {
 | 
|---|
| 88 |   World::AtomSet &atoms = getAtoms();
 | 
|---|
| 89 |   atoms_iter_t res = find_if(atoms.begin_internal(),atoms.end_internal(),boost::bind(&AtomDescriptor_impl::predicate,this,_1));
 | 
|---|
| 90 |   return (res!=atoms.end_internal())?((*res).second):0;
 | 
|---|
| 91 | }
 | 
|---|
| 92 | 
 | 
|---|
| 93 | vector<atom*> AtomDescriptor_impl::findAll() {
 | 
|---|
| 94 |   vector<atom*> res;
 | 
|---|
| 95 |   World::AtomSet atoms = getAtoms();
 | 
|---|
| 96 |   for_each(atoms.begin_internal(),
 | 
|---|
| 97 |            atoms.end_internal(),
 | 
|---|
| 98 |            boost::bind(&AtomDescriptor_impl::checkAndAdd,
 | 
|---|
| 99 |                        this,&res,_1));
 | 
|---|
| 100 |   return res;
 | 
|---|
| 101 | }
 | 
|---|
| 102 | 
 | 
|---|
| 103 | void AtomDescriptor_impl::checkAndAdd(std::vector<atom*> *v,std::pair<atomId_t,atom*> p){
 | 
|---|
| 104 |   if(predicate(p)){
 | 
|---|
| 105 |     v->push_back(p.second);
 | 
|---|
| 106 |   }
 | 
|---|
| 107 | }
 | 
|---|
| 108 | 
 | 
|---|
| 109 | /************************** Universe and Emptyset *****************/
 | 
|---|
| 110 | 
 | 
|---|
| 111 | AtomAllDescriptor_impl::AtomAllDescriptor_impl()
 | 
|---|
| 112 | {}
 | 
|---|
| 113 | 
 | 
|---|
| 114 | AtomAllDescriptor_impl::~AtomAllDescriptor_impl()
 | 
|---|
| 115 | {}
 | 
|---|
| 116 | 
 | 
|---|
| 117 | bool AtomAllDescriptor_impl::predicate(std::pair<atomId_t,atom*>){
 | 
|---|
| 118 |   return true;
 | 
|---|
| 119 | }
 | 
|---|
| 120 | 
 | 
|---|
| 121 | AtomDescriptor AllAtoms(){
 | 
|---|
| 122 |   return AtomDescriptor(AtomDescriptor::impl_ptr(new AtomAllDescriptor_impl));
 | 
|---|
| 123 | }
 | 
|---|
| 124 | 
 | 
|---|
| 125 | AtomNoneDescriptor_impl::AtomNoneDescriptor_impl()
 | 
|---|
| 126 | {}
 | 
|---|
| 127 | 
 | 
|---|
| 128 | AtomNoneDescriptor_impl::~AtomNoneDescriptor_impl()
 | 
|---|
| 129 | {}
 | 
|---|
| 130 | 
 | 
|---|
| 131 | bool AtomNoneDescriptor_impl::predicate(std::pair<atomId_t,atom*>){
 | 
|---|
| 132 |   return false;
 | 
|---|
| 133 | }
 | 
|---|
| 134 | 
 | 
|---|
| 135 | AtomDescriptor NoAtoms(){
 | 
|---|
| 136 |   return AtomDescriptor(AtomDescriptor::impl_ptr(new AtomNoneDescriptor_impl));
 | 
|---|
| 137 | }
 | 
|---|
| 138 | 
 | 
|---|
| 139 | /************************** Operator stuff ************************/
 | 
|---|
| 140 | 
 | 
|---|
| 141 | // AND
 | 
|---|
| 142 | AtomAndDescriptor_impl::AtomAndDescriptor_impl(AtomDescriptor::impl_ptr _lhs, AtomDescriptor::impl_ptr _rhs) :
 | 
|---|
| 143 |     lhs(_lhs), rhs(_rhs)
 | 
|---|
| 144 | {}
 | 
|---|
| 145 | 
 | 
|---|
| 146 | AtomAndDescriptor_impl::~AtomAndDescriptor_impl()
 | 
|---|
| 147 | {}
 | 
|---|
| 148 | 
 | 
|---|
| 149 | bool AtomAndDescriptor_impl::predicate(std::pair<atomId_t,atom*> atom){
 | 
|---|
| 150 |   return lhs->predicate(atom) && rhs->predicate(atom);
 | 
|---|
| 151 | }
 | 
|---|
| 152 | AtomDescriptor operator&&(const AtomDescriptor &lhs, const AtomDescriptor &rhs){
 | 
|---|
| 153 |   AtomDescriptor::impl_ptr newImpl = AtomDescriptor::impl_ptr(new AtomAndDescriptor_impl(lhs.get_impl(),rhs.get_impl()));
 | 
|---|
| 154 |   return AtomDescriptor(newImpl);
 | 
|---|
| 155 | }
 | 
|---|
| 156 | 
 | 
|---|
| 157 | // OR
 | 
|---|
| 158 | AtomOrDescriptor_impl::AtomOrDescriptor_impl(AtomDescriptor::impl_ptr _lhs ,AtomDescriptor::impl_ptr _rhs) :
 | 
|---|
| 159 |     lhs(_lhs), rhs(_rhs)
 | 
|---|
| 160 | {}
 | 
|---|
| 161 | 
 | 
|---|
| 162 | AtomOrDescriptor_impl::~AtomOrDescriptor_impl(){
 | 
|---|
| 163 | }
 | 
|---|
| 164 | 
 | 
|---|
| 165 | bool AtomOrDescriptor_impl::predicate(std::pair<atomId_t,atom*> atom){
 | 
|---|
| 166 |   return lhs->predicate(atom) || rhs->predicate(atom);
 | 
|---|
| 167 | }
 | 
|---|
| 168 | 
 | 
|---|
| 169 | AtomDescriptor  operator||(const AtomDescriptor &lhs, const AtomDescriptor &rhs){
 | 
|---|
| 170 |   AtomDescriptor::impl_ptr newImpl = AtomDescriptor::impl_ptr(new AtomOrDescriptor_impl(lhs.get_impl(),rhs.get_impl()));
 | 
|---|
| 171 |   return AtomDescriptor(newImpl);
 | 
|---|
| 172 | }
 | 
|---|
| 173 | 
 | 
|---|
| 174 | // NOT
 | 
|---|
| 175 | 
 | 
|---|
| 176 | AtomNotDescriptor_impl::AtomNotDescriptor_impl(AtomDescriptor::impl_ptr _arg) :
 | 
|---|
| 177 |     arg(_arg)
 | 
|---|
| 178 | {}
 | 
|---|
| 179 | 
 | 
|---|
| 180 | 
 | 
|---|
| 181 | AtomNotDescriptor_impl::~AtomNotDescriptor_impl()
 | 
|---|
| 182 | {
 | 
|---|
| 183 | }
 | 
|---|
| 184 | 
 | 
|---|
| 185 | bool AtomNotDescriptor_impl::predicate(std::pair<atomId_t,atom*> atom){
 | 
|---|
| 186 |  return !(arg->predicate(atom));
 | 
|---|
| 187 | }
 | 
|---|
| 188 | 
 | 
|---|
| 189 | AtomDescriptor operator!(const AtomDescriptor &arg){
 | 
|---|
| 190 |   AtomDescriptor::impl_ptr newImpl = AtomDescriptor::impl_ptr(new AtomNotDescriptor_impl(arg.get_impl()));
 | 
|---|
| 191 |   return AtomDescriptor(newImpl);
 | 
|---|
| 192 | }
 | 
|---|