| 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 |  * AddEmptyBoundaryAction.cpp
 | 
|---|
| 10 |  *
 | 
|---|
| 11 |  *  Created on: May 8, 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 "Helpers/Log.hpp"
 | 
|---|
| 24 | #include "LinearAlgebra/Vector.hpp"
 | 
|---|
| 25 | #include "World.hpp"
 | 
|---|
| 26 | 
 | 
|---|
| 27 | #include <iostream>
 | 
|---|
| 28 | #include <string>
 | 
|---|
| 29 | #include <vector>
 | 
|---|
| 30 | 
 | 
|---|
| 31 | using namespace std;
 | 
|---|
| 32 | 
 | 
|---|
| 33 | #include "Actions/WorldAction/AddEmptyBoundaryAction.hpp"
 | 
|---|
| 34 | 
 | 
|---|
| 35 | // and construct the stuff
 | 
|---|
| 36 | #include "AddEmptyBoundaryAction.def"
 | 
|---|
| 37 | #include "Action_impl_pre.hpp"
 | 
|---|
| 38 | /** =========== define the function ====================== */
 | 
|---|
| 39 | Action::state_ptr WorldAddEmptyBoundaryAction::performCall() {
 | 
|---|
| 40 |   Vector Min;
 | 
|---|
| 41 |   Vector Max;
 | 
|---|
| 42 |   int j=0;
 | 
|---|
| 43 | 
 | 
|---|
| 44 |   // obtain information
 | 
|---|
| 45 |   getParametersfromValueStorage();
 | 
|---|
| 46 | 
 | 
|---|
| 47 |   // get maximum and minimum
 | 
|---|
| 48 |   vector<atom *> AllAtoms = World::getInstance().getAllAtoms();
 | 
|---|
| 49 |   ASSERT(AllAtoms.size() > 0, "There must be atoms present for AddingEmptyBoundary.");
 | 
|---|
| 50 |   vector<atom *>::iterator AtomRunner = AllAtoms.begin();
 | 
|---|
| 51 |   Min = (*AtomRunner)->getPosition();
 | 
|---|
| 52 |   Max = (*AtomRunner)->getPosition();
 | 
|---|
| 53 |   for (; AtomRunner != AllAtoms.end(); ++AtomRunner) {
 | 
|---|
| 54 |     for (int i=0;i<NDIM;i++) {
 | 
|---|
| 55 |       if ((*AtomRunner)->at(i) > Max[i])
 | 
|---|
| 56 |         Max[i] = (*AtomRunner)->at(i);
 | 
|---|
| 57 |       if ((*AtomRunner)->at(i) < Min[i])
 | 
|---|
| 58 |         Min[i] = (*AtomRunner)->at(i);
 | 
|---|
| 59 |     }
 | 
|---|
| 60 |   }
 | 
|---|
| 61 |   // set new box size
 | 
|---|
| 62 |   double * const cell_size = new double[6];
 | 
|---|
| 63 |   for (j=0;j<6;j++)
 | 
|---|
| 64 |     cell_size[j] = 0.;
 | 
|---|
| 65 |   j=-1;
 | 
|---|
| 66 |   for (int i=0;i<NDIM;i++) {
 | 
|---|
| 67 |     j += i+1;
 | 
|---|
| 68 |     cell_size[j] = (Max[i]-Min[i]+2.*params.boundary[i]);
 | 
|---|
| 69 |   }
 | 
|---|
| 70 |   World::getInstance().setDomain(cell_size);
 | 
|---|
| 71 |   delete[] cell_size;
 | 
|---|
| 72 |   // translate all atoms, such that Min is aty (0,0,0)
 | 
|---|
| 73 |   AtomRunner = AllAtoms.begin();
 | 
|---|
| 74 |   for (; AtomRunner != AllAtoms.end(); ++AtomRunner)
 | 
|---|
| 75 |     *(*AtomRunner) -= Min - params.boundary;
 | 
|---|
| 76 |   return Action::success;
 | 
|---|
| 77 | }
 | 
|---|
| 78 | 
 | 
|---|
| 79 | Action::state_ptr WorldAddEmptyBoundaryAction::performUndo(Action::state_ptr _state) {
 | 
|---|
| 80 | //  ParserLoadXyzState *state = assert_cast<ParserLoadXyzState*>(_state.get());
 | 
|---|
| 81 | 
 | 
|---|
| 82 |   return Action::failure;
 | 
|---|
| 83 | //  string newName = state->mol->getName();
 | 
|---|
| 84 | //  state->mol->setName(state->lastName);
 | 
|---|
| 85 | //
 | 
|---|
| 86 | //  return Action::state_ptr(new ParserLoadXyzState(state->mol,newName));
 | 
|---|
| 87 | }
 | 
|---|
| 88 | 
 | 
|---|
| 89 | Action::state_ptr WorldAddEmptyBoundaryAction::performRedo(Action::state_ptr _state){
 | 
|---|
| 90 |   return Action::failure;
 | 
|---|
| 91 | }
 | 
|---|
| 92 | 
 | 
|---|
| 93 | bool WorldAddEmptyBoundaryAction::canUndo() {
 | 
|---|
| 94 |   return false;
 | 
|---|
| 95 | }
 | 
|---|
| 96 | 
 | 
|---|
| 97 | bool WorldAddEmptyBoundaryAction::shouldUndo() {
 | 
|---|
| 98 |   return false;
 | 
|---|
| 99 | }
 | 
|---|
| 100 | /** =========== end of function ====================== */
 | 
|---|