| 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 | * FillVoidWithMoleculeAction.cpp | 
|---|
| 10 | * | 
|---|
| 11 | *  Created on: May 10, 2010 | 
|---|
| 12 | *      Author: heber | 
|---|
| 13 | */ | 
|---|
| 14 |  | 
|---|
| 15 | // include config.h | 
|---|
| 16 | #ifdef HAVE_CONFIG_H | 
|---|
| 17 | #include <config.h> | 
|---|
| 18 | #endif | 
|---|
| 19 |  | 
|---|
| 20 | #include "Helpers/MemDebug.hpp" | 
|---|
| 21 |  | 
|---|
| 22 | #include "atom.hpp" | 
|---|
| 23 | #include "bondgraph.hpp" | 
|---|
| 24 | #include "boundary.hpp" | 
|---|
| 25 | #include "config.hpp" | 
|---|
| 26 | #include "molecule.hpp" | 
|---|
| 27 | #include "Helpers/Verbose.hpp" | 
|---|
| 28 | #include "World.hpp" | 
|---|
| 29 |  | 
|---|
| 30 | #include "Parser/MpqcParser.hpp" | 
|---|
| 31 | #include "Parser/PcpParser.hpp" | 
|---|
| 32 | #include "Parser/PdbParser.hpp" | 
|---|
| 33 | #include "Parser/TremoloParser.hpp" | 
|---|
| 34 | #include "Parser/XyzParser.hpp" | 
|---|
| 35 | #include "Parser/FormatParserStorage.hpp" | 
|---|
| 36 |  | 
|---|
| 37 | #include <iostream> | 
|---|
| 38 | #include <string> | 
|---|
| 39 |  | 
|---|
| 40 | using namespace std; | 
|---|
| 41 |  | 
|---|
| 42 | #include "Actions/MoleculeAction/FillVoidWithMoleculeAction.hpp" | 
|---|
| 43 |  | 
|---|
| 44 | // and construct the stuff | 
|---|
| 45 | #include "FillVoidWithMoleculeAction.def" | 
|---|
| 46 | #include "Action_impl_pre.hpp" | 
|---|
| 47 | /** =========== define the function ====================== */ | 
|---|
| 48 | Action::state_ptr MoleculeFillVoidWithMoleculeAction::performCall() { | 
|---|
| 49 |  | 
|---|
| 50 | // obtain information | 
|---|
| 51 | getParametersfromValueStorage(); | 
|---|
| 52 |  | 
|---|
| 53 | DoLog(1) && (Log() << Verbose(1) << "Filling Box with water molecules, lengths(" << params.lengths[0] << "," << params.lengths[1] << "," << params.lengths[2] << "), distances (" << params.distances[0] << "," << params.distances[1] << "," << params.distances[2] << "), DoRotate " << params.DoRotate << "." << endl); | 
|---|
| 54 | // construct water molecule | 
|---|
| 55 | molecule *filler = World::getInstance().createMolecule(); | 
|---|
| 56 | std::string FilenameSuffix = params.fillername.string().substr(params.fillername.string().find_last_of('.')+1, params.fillername.string().length()); | 
|---|
| 57 | ifstream input; | 
|---|
| 58 | input.open(params.fillername.string().c_str()); | 
|---|
| 59 | switch (FormatParserStorage::getInstance().getTypeFromSuffix(FilenameSuffix)) { | 
|---|
| 60 | case mpqc: | 
|---|
| 61 | { | 
|---|
| 62 | MpqcParser mpqcparser; | 
|---|
| 63 | mpqcparser.load(&input); | 
|---|
| 64 | break; | 
|---|
| 65 | } | 
|---|
| 66 | case pcp: | 
|---|
| 67 | { | 
|---|
| 68 | PcpParser pcpparser; | 
|---|
| 69 | pcpparser.load(&input); | 
|---|
| 70 | break; | 
|---|
| 71 | } | 
|---|
| 72 | case pdb: | 
|---|
| 73 | { | 
|---|
| 74 | PdbParser pdbparser; | 
|---|
| 75 | pdbparser.load(&input); | 
|---|
| 76 | break; | 
|---|
| 77 | } | 
|---|
| 78 | case tremolo: | 
|---|
| 79 | { | 
|---|
| 80 | TremoloParser tremoloparser; | 
|---|
| 81 | tremoloparser.load(&input); | 
|---|
| 82 | break; | 
|---|
| 83 | } | 
|---|
| 84 | case xyz: | 
|---|
| 85 | { | 
|---|
| 86 | XyzParser xyzparser; | 
|---|
| 87 | xyzparser.load(&input); | 
|---|
| 88 | break; | 
|---|
| 89 | } | 
|---|
| 90 | default: | 
|---|
| 91 | DoeLog(0) && (eLog()<< Verbose(0) << "Could not parse filler molecule from " << params.fillername << "." << endl); | 
|---|
| 92 | break; | 
|---|
| 93 | } | 
|---|
| 94 | World::MoleculeIterator iter = World::getInstance().getMoleculeIter(); | 
|---|
| 95 | for (; iter != World::getInstance().moleculeEnd(); ++iter) | 
|---|
| 96 | filler = *iter; // get last molecule | 
|---|
| 97 | filler->SetNameFromFilename(params.fillername.string().c_str()); | 
|---|
| 98 | World::getInstance().getConfig()->BG->ConstructBondGraph(filler); | 
|---|
| 99 |  | 
|---|
| 100 | // call routine | 
|---|
| 101 | double distance[NDIM]; | 
|---|
| 102 | for (int i=0;i<NDIM;i++) | 
|---|
| 103 | distance[i] = params.distances[i]; | 
|---|
| 104 | FillVoidWithMolecule(filler, *(World::getInstance().getConfig()), distance, params.lengths[0], params.lengths[1], params.lengths[2], params.DoRotate); | 
|---|
| 105 |  | 
|---|
| 106 | return Action::success; | 
|---|
| 107 | } | 
|---|
| 108 |  | 
|---|
| 109 | Action::state_ptr MoleculeFillVoidWithMoleculeAction::performUndo(Action::state_ptr _state) { | 
|---|
| 110 | //  MoleculeFillVoidWithMoleculeState *state = assert_cast<MoleculeFillVoidWithMoleculeState*>(_state.get()); | 
|---|
| 111 |  | 
|---|
| 112 | //  string newName = state->mol->getName(); | 
|---|
| 113 | //  state->mol->setName(state->lastName); | 
|---|
| 114 |  | 
|---|
| 115 | return Action::failure; | 
|---|
| 116 | } | 
|---|
| 117 |  | 
|---|
| 118 | Action::state_ptr MoleculeFillVoidWithMoleculeAction::performRedo(Action::state_ptr _state){ | 
|---|
| 119 | // Undo and redo have to do the same for this action | 
|---|
| 120 | return performUndo(_state); | 
|---|
| 121 | } | 
|---|
| 122 |  | 
|---|
| 123 | bool MoleculeFillVoidWithMoleculeAction::canUndo() { | 
|---|
| 124 | return false; | 
|---|
| 125 | } | 
|---|
| 126 |  | 
|---|
| 127 | bool MoleculeFillVoidWithMoleculeAction::shouldUndo() { | 
|---|
| 128 | return false; | 
|---|
| 129 | } | 
|---|
| 130 | /** =========== end of function ====================== */ | 
|---|