| [eee966] | 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(); | 
|---|
| [e4afb4] | 56 | std::string FilenameSuffix = params.fillername.string().substr(params.fillername.string().find_last_of('.')+1, params.fillername.string().length()); | 
|---|
| [eee966] | 57 | ifstream input; | 
|---|
| [e4afb4] | 58 | input.open(params.fillername.string().c_str()); | 
|---|
| [eee966] | 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 | 
|---|
| [e4afb4] | 97 | filler->SetNameFromFilename(params.fillername.string().c_str()); | 
|---|
| [eee966] | 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 ====================== */ | 
|---|