| 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 |  * VerletIntegrationAction.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 "CodePatterns/MemDebug.hpp"
 | 
|---|
| 21 | 
 | 
|---|
| 22 | #include "CodePatterns/Log.hpp"
 | 
|---|
| 23 | #include "molecule.hpp"
 | 
|---|
| 24 | #include "CodePatterns/Verbose.hpp"
 | 
|---|
| 25 | #include "World.hpp"
 | 
|---|
| 26 | 
 | 
|---|
| 27 | #include <iostream>
 | 
|---|
| 28 | #include <fstream>
 | 
|---|
| 29 | #include <string>
 | 
|---|
| 30 | 
 | 
|---|
| 31 | using namespace std;
 | 
|---|
| 32 | 
 | 
|---|
| 33 | #include "Actions/MoleculeAction/VerletIntegrationAction.hpp"
 | 
|---|
| 34 | 
 | 
|---|
| 35 | // and construct the stuff
 | 
|---|
| 36 | #include "VerletIntegrationAction.def"
 | 
|---|
| 37 | #include "Action_impl_pre.hpp"
 | 
|---|
| 38 | /** =========== define the function ====================== */
 | 
|---|
| 39 | Action::state_ptr MoleculeVerletIntegrationAction::performCall() {
 | 
|---|
| 40 |   molecule *mol = NULL;
 | 
|---|
| 41 | 
 | 
|---|
| 42 |   // obtain information
 | 
|---|
| 43 |   getParametersfromValueStorage();
 | 
|---|
| 44 | 
 | 
|---|
| 45 |   for (World::MoleculeSelectionIterator iter = World::getInstance().beginMoleculeSelection(); iter != World::getInstance().endMoleculeSelection(); ++iter) {
 | 
|---|
| 46 |     mol = iter->second;
 | 
|---|
| 47 |     DoLog(1) && (Log() << Verbose(1) << "Parsing forces file and Verlet integrating." << endl);
 | 
|---|
| 48 |     // TODO: sollte besser stream nutzen, nicht filename direkt (es sei denn, ist prefix), besser fuer unit test
 | 
|---|
| 49 |     char outputname[MAXSTRINGSIZE];
 | 
|---|
| 50 |     strcpy(outputname, params.forcesfile.string().c_str());
 | 
|---|
| 51 |     if (!mol->VerletForceIntegration(outputname, *(World::getInstance().getConfig()), 0))
 | 
|---|
| 52 |       DoLog(2) && (Log() << Verbose(2) << "File " << params.forcesfile << " not found." << endl);
 | 
|---|
| 53 |     else
 | 
|---|
| 54 |       DoLog(2) && (Log() << Verbose(2) << "File " << params.forcesfile << " found and parsed." << endl);
 | 
|---|
| 55 |   }
 | 
|---|
| 56 |   return Action::success;
 | 
|---|
| 57 | }
 | 
|---|
| 58 | 
 | 
|---|
| 59 | Action::state_ptr MoleculeVerletIntegrationAction::performUndo(Action::state_ptr _state) {
 | 
|---|
| 60 | //  MoleculeVerletIntegrationState *state = assert_cast<MoleculeVerletIntegrationState*>(_state.get());
 | 
|---|
| 61 | 
 | 
|---|
| 62 | //  string newName = state->mol->getName();
 | 
|---|
| 63 | //  state->mol->setName(state->lastName);
 | 
|---|
| 64 | 
 | 
|---|
| 65 |   return Action::failure;
 | 
|---|
| 66 | }
 | 
|---|
| 67 | 
 | 
|---|
| 68 | Action::state_ptr MoleculeVerletIntegrationAction::performRedo(Action::state_ptr _state){
 | 
|---|
| 69 |   // Undo and redo have to do the same for this action
 | 
|---|
| 70 |   return performUndo(_state);
 | 
|---|
| 71 | }
 | 
|---|
| 72 | 
 | 
|---|
| 73 | bool MoleculeVerletIntegrationAction::canUndo() {
 | 
|---|
| 74 |   return false;
 | 
|---|
| 75 | }
 | 
|---|
| 76 | 
 | 
|---|
| 77 | bool MoleculeVerletIntegrationAction::shouldUndo() {
 | 
|---|
| 78 |   return false;
 | 
|---|
| 79 | }
 | 
|---|
| 80 | /** =========== end of function ====================== */
 | 
|---|