| 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 | 
 | 
|---|
| 15 | #include "molecule.hpp"
 | 
|---|
| 16 | 
 | 
|---|
| 17 | #include <boost/bind.hpp>
 | 
|---|
| 18 | #include <cassert>
 | 
|---|
| 19 | #include <iostream>
 | 
|---|
| 20 | 
 | 
|---|
| 21 | using namespace std;
 | 
|---|
| 22 | 
 | 
|---|
| 23 | typedef World::MoleculeSet::iterator molecules_iter_t;
 | 
|---|
| 24 | 
 | 
|---|
| 25 | /************************ Forwarding object **************************************/
 | 
|---|
| 26 | 
 | 
|---|
| 27 | 
 | 
|---|
| 28 | MoleculeDescriptor::MoleculeDescriptor(impl_ptr _impl) :
 | 
|---|
| 29 |     impl(_impl)
 | 
|---|
| 30 | {}
 | 
|---|
| 31 | 
 | 
|---|
| 32 | MoleculeDescriptor::MoleculeDescriptor(const MoleculeDescriptor& src) :
 | 
|---|
| 33 |     impl(src.get_impl())
 | 
|---|
| 34 | {}
 | 
|---|
| 35 | 
 | 
|---|
| 36 | MoleculeDescriptor::~MoleculeDescriptor()
 | 
|---|
| 37 | {}
 | 
|---|
| 38 | 
 | 
|---|
| 39 | MoleculeDescriptor& MoleculeDescriptor::operator=(MoleculeDescriptor &src){
 | 
|---|
| 40 |   if(&src!=this) {
 | 
|---|
| 41 |     impl=src.get_impl();
 | 
|---|
| 42 |   }
 | 
|---|
| 43 |   return *this;
 | 
|---|
| 44 | }
 | 
|---|
| 45 | 
 | 
|---|
| 46 | molecule* MoleculeDescriptor::find(){
 | 
|---|
| 47 |   return impl->find();
 | 
|---|
| 48 | }
 | 
|---|
| 49 | 
 | 
|---|
| 50 | std::vector<molecule*> MoleculeDescriptor::findAll(){
 | 
|---|
| 51 |   return impl->findAll();
 | 
|---|
| 52 | }
 | 
|---|
| 53 | 
 | 
|---|
| 54 | MoleculeDescriptor::impl_ptr MoleculeDescriptor::get_impl() const{
 | 
|---|
| 55 |   return impl;
 | 
|---|
| 56 | }
 | 
|---|
| 57 | 
 | 
|---|
| 58 | 
 | 
|---|
| 59 | 
 | 
|---|
| 60 | 
 | 
|---|
| 61 | /**************************** implementation ********************/
 | 
|---|
| 62 | 
 | 
|---|
| 63 | MoleculeDescriptor_impl::MoleculeDescriptor_impl()
 | 
|---|
| 64 | {
 | 
|---|
| 65 | }
 | 
|---|
| 66 | 
 | 
|---|
| 67 | MoleculeDescriptor_impl::~MoleculeDescriptor_impl()
 | 
|---|
| 68 | {
 | 
|---|
| 69 | }
 | 
|---|
| 70 | 
 | 
|---|
| 71 | World::MoleculeSet& MoleculeDescriptor_impl::getMolecules(){
 | 
|---|
| 72 |   return World::getInstance().molecules;
 | 
|---|
| 73 | }
 | 
|---|
| 74 | 
 | 
|---|
| 75 | molecule* MoleculeDescriptor_impl::find() {
 | 
|---|
| 76 |   World::MoleculeSet molecules = getMolecules();
 | 
|---|
| 77 |   molecules_iter_t res = find_if(molecules.begin(),molecules.end(),boost::bind(&MoleculeDescriptor_impl::predicate,this,_1));
 | 
|---|
| 78 |   return (res!=molecules.end())?((*res).second):0;
 | 
|---|
| 79 | }
 | 
|---|
| 80 | 
 | 
|---|
| 81 | vector<molecule*> MoleculeDescriptor_impl::findAll() {
 | 
|---|
| 82 |   vector<molecule*> res;
 | 
|---|
| 83 |   World::MoleculeSet molecules = getMolecules();
 | 
|---|
| 84 |   molecules_iter_t iter;
 | 
|---|
| 85 |   for(iter=molecules.begin();iter!=molecules.end();++iter) {
 | 
|---|
| 86 |     if(predicate(*iter)){
 | 
|---|
| 87 |       res.push_back((*iter).second);
 | 
|---|
| 88 |     }
 | 
|---|
| 89 |   }
 | 
|---|
| 90 |   return res;
 | 
|---|
| 91 | }
 | 
|---|
| 92 | 
 | 
|---|
| 93 | /************************** Universe and Emptyset *****************/
 | 
|---|
| 94 | 
 | 
|---|
| 95 | MoleculeAllDescriptor_impl::MoleculeAllDescriptor_impl()
 | 
|---|
| 96 | {}
 | 
|---|
| 97 | 
 | 
|---|
| 98 | MoleculeAllDescriptor_impl::~MoleculeAllDescriptor_impl()
 | 
|---|
| 99 | {}
 | 
|---|
| 100 | 
 | 
|---|
| 101 | bool MoleculeAllDescriptor_impl::predicate(std::pair<moleculeId_t,molecule*>){
 | 
|---|
| 102 |   return true;
 | 
|---|
| 103 | }
 | 
|---|
| 104 | 
 | 
|---|
| 105 | MoleculeDescriptor AllMolecules(){
 | 
|---|
| 106 |   return MoleculeDescriptor(MoleculeDescriptor::impl_ptr(new MoleculeAllDescriptor_impl));
 | 
|---|
| 107 | }
 | 
|---|
| 108 | 
 | 
|---|
| 109 | MoleculeNoneDescriptor_impl::MoleculeNoneDescriptor_impl()
 | 
|---|
| 110 | {}
 | 
|---|
| 111 | 
 | 
|---|
| 112 | MoleculeNoneDescriptor_impl::~MoleculeNoneDescriptor_impl()
 | 
|---|
| 113 | {}
 | 
|---|
| 114 | 
 | 
|---|
| 115 | bool MoleculeNoneDescriptor_impl::predicate(std::pair<moleculeId_t,molecule*>){
 | 
|---|
| 116 |   return false;
 | 
|---|
| 117 | }
 | 
|---|
| 118 | 
 | 
|---|
| 119 | MoleculeDescriptor NoMolecules(){
 | 
|---|
| 120 |   return MoleculeDescriptor(MoleculeDescriptor::impl_ptr(new MoleculeNoneDescriptor_impl));
 | 
|---|
| 121 | }
 | 
|---|
| 122 | 
 | 
|---|
| 123 | /************************** Operator stuff ************************/
 | 
|---|
| 124 | 
 | 
|---|
| 125 | // AND
 | 
|---|
| 126 | MoleculeAndDescriptor_impl::MoleculeAndDescriptor_impl(MoleculeDescriptor::impl_ptr _lhs, MoleculeDescriptor::impl_ptr _rhs) :
 | 
|---|
| 127 |     lhs(_lhs), rhs(_rhs)
 | 
|---|
| 128 | {}
 | 
|---|
| 129 | 
 | 
|---|
| 130 | MoleculeAndDescriptor_impl::~MoleculeAndDescriptor_impl()
 | 
|---|
| 131 | {}
 | 
|---|
| 132 | 
 | 
|---|
| 133 | bool MoleculeAndDescriptor_impl::predicate(std::pair<moleculeId_t,molecule*> molecule){
 | 
|---|
| 134 |   return lhs->predicate(molecule) && rhs->predicate(molecule);
 | 
|---|
| 135 | }
 | 
|---|
| 136 | MoleculeDescriptor operator&&(const MoleculeDescriptor &lhs, const MoleculeDescriptor &rhs){
 | 
|---|
| 137 |   MoleculeDescriptor::impl_ptr newImpl = MoleculeDescriptor::impl_ptr(new MoleculeAndDescriptor_impl(lhs.get_impl(),rhs.get_impl()));
 | 
|---|
| 138 |   return MoleculeDescriptor(newImpl);
 | 
|---|
| 139 | }
 | 
|---|
| 140 | 
 | 
|---|
| 141 | // OR
 | 
|---|
| 142 | MoleculeOrDescriptor_impl::MoleculeOrDescriptor_impl(MoleculeDescriptor::impl_ptr _lhs ,MoleculeDescriptor::impl_ptr _rhs) :
 | 
|---|
| 143 |     lhs(_lhs), rhs(_rhs)
 | 
|---|
| 144 | {}
 | 
|---|
| 145 | 
 | 
|---|
| 146 | MoleculeOrDescriptor_impl::~MoleculeOrDescriptor_impl(){
 | 
|---|
| 147 | }
 | 
|---|
| 148 | 
 | 
|---|
| 149 | bool MoleculeOrDescriptor_impl::predicate(std::pair<moleculeId_t,molecule*> molecule){
 | 
|---|
| 150 |   return lhs->predicate(molecule) || rhs->predicate(molecule);
 | 
|---|
| 151 | }
 | 
|---|
| 152 | 
 | 
|---|
| 153 | MoleculeDescriptor  operator||(const MoleculeDescriptor &lhs, const MoleculeDescriptor &rhs){
 | 
|---|
| 154 |   MoleculeDescriptor::impl_ptr newImpl = MoleculeDescriptor::impl_ptr(new MoleculeOrDescriptor_impl(lhs.get_impl(),rhs.get_impl()));
 | 
|---|
| 155 |   return MoleculeDescriptor(newImpl);
 | 
|---|
| 156 | }
 | 
|---|
| 157 | 
 | 
|---|
| 158 | // NOT
 | 
|---|
| 159 | 
 | 
|---|
| 160 | MoleculeNotDescriptor_impl::MoleculeNotDescriptor_impl(MoleculeDescriptor::impl_ptr _arg) :
 | 
|---|
| 161 |     arg(_arg)
 | 
|---|
| 162 | {}
 | 
|---|
| 163 | 
 | 
|---|
| 164 | 
 | 
|---|
| 165 | MoleculeNotDescriptor_impl::~MoleculeNotDescriptor_impl()
 | 
|---|
| 166 | {
 | 
|---|
| 167 | }
 | 
|---|
| 168 | 
 | 
|---|
| 169 | bool MoleculeNotDescriptor_impl::predicate(std::pair<moleculeId_t,molecule*> molecule){
 | 
|---|
| 170 |  return !(arg->predicate(molecule));
 | 
|---|
| 171 | }
 | 
|---|
| 172 | 
 | 
|---|
| 173 | MoleculeDescriptor operator!(const MoleculeDescriptor &arg){
 | 
|---|
| 174 |   MoleculeDescriptor::impl_ptr newImpl = MoleculeDescriptor::impl_ptr(new MoleculeNotDescriptor_impl(arg.get_impl()));
 | 
|---|
| 175 |   return MoleculeDescriptor(newImpl);
 | 
|---|
| 176 | }
 | 
|---|