| 1 | /* | 
|---|
| 2 | * Project: MoleCuilder | 
|---|
| 3 | * Description: creates and alters molecular systems | 
|---|
| 4 | * Copyright (C)  2012 University of Bonn. All rights reserved. | 
|---|
| 5 | * | 
|---|
| 6 | * | 
|---|
| 7 | *   This file is part of MoleCuilder. | 
|---|
| 8 | * | 
|---|
| 9 | *    MoleCuilder is free software: you can redistribute it and/or modify | 
|---|
| 10 | *    it under the terms of the GNU General Public License as published by | 
|---|
| 11 | *    the Free Software Foundation, either version 2 of the License, or | 
|---|
| 12 | *    (at your option) any later version. | 
|---|
| 13 | * | 
|---|
| 14 | *    MoleCuilder is distributed in the hope that it will be useful, | 
|---|
| 15 | *    but WITHOUT ANY WARRANTY; without even the implied warranty of | 
|---|
| 16 | *    MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the | 
|---|
| 17 | *    GNU General Public License for more details. | 
|---|
| 18 | * | 
|---|
| 19 | *    You should have received a copy of the GNU General Public License | 
|---|
| 20 | *    along with MoleCuilder.  If not, see <http://www.gnu.org/licenses/>. | 
|---|
| 21 | */ | 
|---|
| 22 |  | 
|---|
| 23 | /* | 
|---|
| 24 | * UndoRedoHelpers.cpp | 
|---|
| 25 | * | 
|---|
| 26 | *  Created on: Apr 5, 2012 | 
|---|
| 27 | *      Author: heber | 
|---|
| 28 | */ | 
|---|
| 29 |  | 
|---|
| 30 |  | 
|---|
| 31 | // include config.h | 
|---|
| 32 | #ifdef HAVE_CONFIG_H | 
|---|
| 33 | #include <config.h> | 
|---|
| 34 | #endif | 
|---|
| 35 |  | 
|---|
| 36 | #include "CodePatterns/MemDebug.hpp" | 
|---|
| 37 |  | 
|---|
| 38 | #include "UndoRedoHelpers.hpp" | 
|---|
| 39 |  | 
|---|
| 40 | #include <boost/foreach.hpp> | 
|---|
| 41 |  | 
|---|
| 42 | #include "Atom/atom.hpp" | 
|---|
| 43 | #include "molecule.hpp" | 
|---|
| 44 | #include "Descriptors/AtomIdDescriptor.hpp" | 
|---|
| 45 | #include "Descriptors/MoleculeIdDescriptor.hpp" | 
|---|
| 46 | #include "CodePatterns/Assert.hpp" | 
|---|
| 47 | #include "CodePatterns/Log.hpp" | 
|---|
| 48 | #include "World.hpp" | 
|---|
| 49 |  | 
|---|
| 50 | bool MoleCuilder::AddAtomsFromAtomicInfo(std::vector<AtomicInfo> &atoms) | 
|---|
| 51 | { | 
|---|
| 52 | size_t i=0; | 
|---|
| 53 | for (; i<atoms.size(); ++i) { | 
|---|
| 54 | // re-create the atom | 
|---|
| 55 | LOG(3, "DEBUG: Re-adding atom " << atoms[i].getId() << "."); | 
|---|
| 56 | atom *Walker = World::getInstance().createAtom(); | 
|---|
| 57 | if (!atoms[i].setAtom(*Walker)) { | 
|---|
| 58 | ELOG(1, "Failed to set id."); | 
|---|
| 59 | World::getInstance().destroyAtom(Walker); | 
|---|
| 60 | break; | 
|---|
| 61 | } | 
|---|
| 62 | } | 
|---|
| 63 | if (i<atoms.size()) { | 
|---|
| 64 | // remove all previous ones, too | 
|---|
| 65 | for (size_t j=0;j<i;++j) | 
|---|
| 66 | World::getInstance().destroyAtom(atoms[j].getId()); | 
|---|
| 67 | // and announce the failure | 
|---|
| 68 | return false; | 
|---|
| 69 | } | 
|---|
| 70 | return true; | 
|---|
| 71 | } | 
|---|
| 72 |  | 
|---|
| 73 | void MoleCuilder::RemoveAtomsFromAtomicInfo(std::vector<AtomicInfo> &atoms) | 
|---|
| 74 | { | 
|---|
| 75 | BOOST_FOREACH(const AtomicInfo &_atom, atoms) { | 
|---|
| 76 | World::getInstance().destroyAtom(_atom.getId()); | 
|---|
| 77 | } | 
|---|
| 78 | } | 
|---|
| 79 |  | 
|---|
| 80 | void MoleCuilder::StoreBondInformationFromAtoms( | 
|---|
| 81 | const std::vector<const atom*> &atoms, | 
|---|
| 82 | std::vector< BondInfo > &bonds) | 
|---|
| 83 | { | 
|---|
| 84 | ASSERT( bonds.empty(), | 
|---|
| 85 | "StoreBondInformationFromAtoms() - give bonds vector is not empty."); | 
|---|
| 86 | bonds.reserve(atoms.size()*4); | 
|---|
| 87 | for (std::vector<const atom*>::const_iterator atomiter = atoms.begin(); | 
|---|
| 88 | atomiter != atoms.end(); ++atomiter) { | 
|---|
| 89 | const BondList & _atom_bonds = (*atomiter)->getListOfBonds(); | 
|---|
| 90 | for(BondList::const_iterator iter = _atom_bonds.begin(); iter != _atom_bonds.end(); ++iter) | 
|---|
| 91 | bonds.push_back( BondInfo(*iter) ); | 
|---|
| 92 | } | 
|---|
| 93 | } | 
|---|
| 94 |  | 
|---|
| 95 | bool MoleCuilder::AddBondsFromBondInfo(const std::vector< BondInfo > &bonds) | 
|---|
| 96 | { | 
|---|
| 97 | bool status = true; | 
|---|
| 98 | for(std::vector< BondInfo >::const_iterator iter = bonds.begin(); | 
|---|
| 99 | iter != bonds.end(); ++iter) | 
|---|
| 100 | if (!(*iter).RecreateBond()) | 
|---|
| 101 | status = false; | 
|---|
| 102 | return status; | 
|---|
| 103 | } | 
|---|
| 104 |  | 
|---|
| 105 | void MoleCuilder::SetAtomsFromAtomicInfo(std::vector<AtomicInfo> &movedatoms) | 
|---|
| 106 | { | 
|---|
| 107 | BOOST_FOREACH( const AtomicInfo &_atominfo, movedatoms) { | 
|---|
| 108 | const atomId_t id = _atominfo.getId(); | 
|---|
| 109 | atom * const _atom = World::getInstance().getAtom(AtomById(id)); | 
|---|
| 110 | ASSERT( _atom != NULL, | 
|---|
| 111 | "MoleCuilder::SetAtomsFromAtomicInfo() - cannot find atom with id " | 
|---|
| 112 | +toString(id)+" in the world."); | 
|---|
| 113 | _atominfo.setAtom( *_atom ); | 
|---|
| 114 | } | 
|---|
| 115 | } | 
|---|
| 116 |  | 
|---|
| 117 | void MoleCuilder::SelectAtomsFromAtomicInfo(std::vector<AtomicInfo> &movedatoms) | 
|---|
| 118 | { | 
|---|
| 119 | BOOST_FOREACH( const AtomicInfo &_atominfo, movedatoms) { | 
|---|
| 120 | const atomId_t id = _atominfo.getId(); | 
|---|
| 121 | World::getInstance().selectAtom(id); | 
|---|
| 122 | } | 
|---|
| 123 | } | 
|---|
| 124 |  | 
|---|
| 125 | void MoleCuilder::ResetAtomPosition(const std::vector<AtomicInfo> &movedatoms, const std::vector<Vector> &MovedToVector) | 
|---|
| 126 | { | 
|---|
| 127 | std::vector<Vector>::const_iterator positer = MovedToVector.begin(); | 
|---|
| 128 | ASSERT(movedatoms.size() == MovedToVector.size(), | 
|---|
| 129 | "MoleCuilder::ResetAtomPosition() -  the number of atoms " | 
|---|
| 130 | +toString(movedatoms.size())+" and the number of positions " | 
|---|
| 131 | +toString(MovedToVector.size())+" is not the same."); | 
|---|
| 132 | BOOST_FOREACH( const AtomicInfo &_atominfo, movedatoms) { | 
|---|
| 133 | const atomId_t id = _atominfo.getId(); | 
|---|
| 134 | atom * const _atom = World::getInstance().getAtom(AtomById(id)); | 
|---|
| 135 | ASSERT( _atom != NULL, | 
|---|
| 136 | "FillSphericalSurfaceAction::performRedo() - cannot find atom with id " | 
|---|
| 137 | +toString(id)+" in the world."); | 
|---|
| 138 | _atom->setPosition( *positer ); | 
|---|
| 139 | ++positer; | 
|---|
| 140 | } | 
|---|
| 141 | } | 
|---|
| 142 |  | 
|---|
| 143 | void MoleCuilder::RemoveMoleculesWithAtomsByIds(const std::vector<moleculeId_t> &ids) | 
|---|
| 144 | { | 
|---|
| 145 | for (std::vector<moleculeId_t>::const_iterator iter = ids.begin(); | 
|---|
| 146 | iter != ids.end(); ++iter) { | 
|---|
| 147 | molecule * const mol = World::getInstance().getMolecule(MoleculeById(*iter)); | 
|---|
| 148 | if (mol != NULL) { | 
|---|
| 149 | mol->removeAtomsinMolecule(); | 
|---|
| 150 | World::getInstance().destroyMolecule(mol); | 
|---|
| 151 | } | 
|---|
| 152 | } | 
|---|
| 153 | } | 
|---|