| [97ebf8] | 1 | /*
 | 
|---|
 | 2 |  * FragmentationAction.cpp
 | 
|---|
 | 3 |  *
 | 
|---|
 | 4 |  *  Created on: May 9, 2010
 | 
|---|
 | 5 |  *      Author: heber
 | 
|---|
 | 6 |  */
 | 
|---|
 | 7 | 
 | 
|---|
| [112b09] | 8 | #include "Helpers/MemDebug.hpp"
 | 
|---|
 | 9 | 
 | 
|---|
| [97ebf8] | 10 | #include "Actions/FragmentationAction/FragmentationAction.hpp"
 | 
|---|
| [0430e3] | 11 | #include "Actions/ActionRegistry.hpp"
 | 
|---|
| [97ebf8] | 12 | #include "atom.hpp"
 | 
|---|
| [a3fded] | 13 | #include "bondgraph.hpp"
 | 
|---|
| [97ebf8] | 14 | #include "config.hpp"
 | 
|---|
| [952f38] | 15 | #include "Helpers/Log.hpp"
 | 
|---|
| [97ebf8] | 16 | #include "molecule.hpp"
 | 
|---|
 | 17 | #include "Descriptors/MoleculeDescriptor.hpp"
 | 
|---|
 | 18 | #include "stackclass.hpp"
 | 
|---|
 | 19 | #include "World.hpp"
 | 
|---|
 | 20 | 
 | 
|---|
 | 21 | #include <iostream>
 | 
|---|
 | 22 | #include <string>
 | 
|---|
 | 23 | 
 | 
|---|
 | 24 | using namespace std;
 | 
|---|
 | 25 | 
 | 
|---|
 | 26 | #include "UIElements/UIFactory.hpp"
 | 
|---|
 | 27 | #include "UIElements/Dialog.hpp"
 | 
|---|
| [861874] | 28 | #include "Actions/ValueStorage.hpp"
 | 
|---|
| [97ebf8] | 29 | 
 | 
|---|
| [e4b5de] | 30 | const char FragmentationFragmentationAction::NAME[] = "fragment-mol";
 | 
|---|
| [97ebf8] | 31 | 
 | 
|---|
 | 32 | FragmentationFragmentationAction::FragmentationFragmentationAction() :
 | 
|---|
 | 33 |   Action(NAME)
 | 
|---|
 | 34 | {}
 | 
|---|
 | 35 | 
 | 
|---|
 | 36 | FragmentationFragmentationAction::~FragmentationFragmentationAction()
 | 
|---|
 | 37 | {}
 | 
|---|
 | 38 | 
 | 
|---|
| [4ff081] | 39 | void FragmentationFragmentation(std::string &path, double distance, int order) {
 | 
|---|
 | 40 |   ValueStorage::getInstance().setCurrentValue(FragmentationFragmentationAction::NAME, path);
 | 
|---|
 | 41 |   ValueStorage::getInstance().setCurrentValue("distance", distance);
 | 
|---|
 | 42 |   ValueStorage::getInstance().setCurrentValue("order", order);
 | 
|---|
 | 43 |   ActionRegistry::getInstance().getActionByName(FragmentationFragmentationAction::NAME)->call(Action::NonInteractive);
 | 
|---|
 | 44 | };
 | 
|---|
 | 45 | 
 | 
|---|
| [047878] | 46 | Dialog* FragmentationFragmentationAction::fillDialog(Dialog *dialog) {
 | 
|---|
 | 47 |   ASSERT(dialog,"No Dialog given when filling action dialog");
 | 
|---|
| [70d9b9] | 48 | 
 | 
|---|
 | 49 |   dialog->queryString(NAME, ValueStorage::getInstance().getDescription(NAME));
 | 
|---|
 | 50 |   dialog->queryDouble("distance", ValueStorage::getInstance().getDescription("distance"));
 | 
|---|
 | 51 |   dialog->queryInt("order", ValueStorage::getInstance().getDescription("order"));
 | 
|---|
 | 52 | 
 | 
|---|
 | 53 |   return dialog;
 | 
|---|
 | 54 | }
 | 
|---|
 | 55 | 
 | 
|---|
 | 56 | Action::state_ptr FragmentationFragmentationAction::performCall() {
 | 
|---|
| [e4b5de] | 57 |   clock_t start,end;
 | 
|---|
 | 58 |   molecule *mol = NULL;
 | 
|---|
 | 59 |   double distance = -1.;
 | 
|---|
 | 60 |   int order = 0;
 | 
|---|
| [35b698] | 61 |   std::string path;
 | 
|---|
| [e4b5de] | 62 |   config *configuration = World::getInstance().getConfig();
 | 
|---|
 | 63 |   int ExitFlag = 0;
 | 
|---|
 | 64 | 
 | 
|---|
| [70d9b9] | 65 |   ValueStorage::getInstance().queryCurrentValue(NAME, path);
 | 
|---|
 | 66 |   ValueStorage::getInstance().queryCurrentValue("distance", distance);
 | 
|---|
 | 67 |   ValueStorage::getInstance().queryCurrentValue("order", order);
 | 
|---|
| [97ebf8] | 68 | 
 | 
|---|
| [70d9b9] | 69 |   for (World::MoleculeSelectionIterator iter = World::getInstance().beginMoleculeSelection(); iter != World::getInstance().endMoleculeSelection(); ++iter) {
 | 
|---|
 | 70 |     mol = iter->second;
 | 
|---|
| [e4b5de] | 71 |     ASSERT(mol != NULL, "No molecule has been picked for fragmentation.");
 | 
|---|
 | 72 |     DoLog(0) && (Log() << Verbose(0) << "Fragmenting molecule with bond distance " << distance << " angstroem, order of " << order << "." << endl);
 | 
|---|
 | 73 |     DoLog(0) && (Log() << Verbose(0) << "Creating connection matrix..." << endl);
 | 
|---|
 | 74 |     start = clock();
 | 
|---|
 | 75 |     mol->CreateAdjacencyList(distance, configuration->GetIsAngstroem(), &BondGraph::CovalentMinMaxDistance, NULL);
 | 
|---|
 | 76 |     DoLog(0) && (Log() << Verbose(0) << "Fragmenting molecule with current connection matrix ..." << endl);
 | 
|---|
 | 77 |     if (mol->hasBondStructure()) {
 | 
|---|
| [35b698] | 78 |       ExitFlag = mol->FragmentMolecule(order, path);
 | 
|---|
| [e4b5de] | 79 |     }
 | 
|---|
 | 80 |     World::getInstance().setExitFlag(ExitFlag);
 | 
|---|
 | 81 |     end = clock();
 | 
|---|
 | 82 |     DoLog(0) && (Log() << Verbose(0) << "Clocks for this operation: " << (end-start) << ", time: " << ((double)(end-start)/CLOCKS_PER_SEC) << "s." << endl);
 | 
|---|
| [97ebf8] | 83 |   }
 | 
|---|
| [70d9b9] | 84 |   return Action::success;
 | 
|---|
| [97ebf8] | 85 | }
 | 
|---|
 | 86 | 
 | 
|---|
 | 87 | Action::state_ptr FragmentationFragmentationAction::performUndo(Action::state_ptr _state) {
 | 
|---|
| [70d9b9] | 88 |   return Action::success;
 | 
|---|
| [97ebf8] | 89 | }
 | 
|---|
 | 90 | 
 | 
|---|
 | 91 | Action::state_ptr FragmentationFragmentationAction::performRedo(Action::state_ptr _state){
 | 
|---|
| [70d9b9] | 92 |   return Action::success;
 | 
|---|
| [97ebf8] | 93 | }
 | 
|---|
 | 94 | 
 | 
|---|
 | 95 | bool FragmentationFragmentationAction::canUndo() {
 | 
|---|
| [70d9b9] | 96 |   return true;
 | 
|---|
| [97ebf8] | 97 | }
 | 
|---|
 | 98 | 
 | 
|---|
 | 99 | bool FragmentationFragmentationAction::shouldUndo() {
 | 
|---|
| [70d9b9] | 100 |   return true;
 | 
|---|
| [97ebf8] | 101 | }
 | 
|---|
 | 102 | 
 | 
|---|
 | 103 | const string FragmentationFragmentationAction::getName() {
 | 
|---|
 | 104 |   return NAME;
 | 
|---|
 | 105 | }
 | 
|---|