| [bcf653] | 1 | /*
 | 
|---|
 | 2 |  * Project: MoleCuilder
 | 
|---|
 | 3 |  * Description: creates and alters molecular systems
 | 
|---|
| [0aa122] | 4 |  * Copyright (C)  2010-2012 University of Bonn. All rights reserved.
 | 
|---|
| [bcf653] | 5 |  * Please see the LICENSE file or "Copyright notice" in builder.cpp for details.
 | 
|---|
 | 6 |  */
 | 
|---|
 | 7 | 
 | 
|---|
| [22c44bf] | 8 | /*
 | 
|---|
 | 9 |  * RotateAroundOriginByAngleAction.cpp
 | 
|---|
 | 10 |  *
 | 
|---|
 | 11 |  *  Created on: Aug 06, 2010
 | 
|---|
 | 12 |  *      Author: heber
 | 
|---|
 | 13 |  */
 | 
|---|
 | 14 | 
 | 
|---|
| [bf3817] | 15 | // include config.h
 | 
|---|
 | 16 | #ifdef HAVE_CONFIG_H
 | 
|---|
 | 17 | #include <config.h>
 | 
|---|
 | 18 | #endif
 | 
|---|
 | 19 | 
 | 
|---|
| [ad011c] | 20 | #include "CodePatterns/MemDebug.hpp"
 | 
|---|
| [22c44bf] | 21 | 
 | 
|---|
| [ad011c] | 22 | #include "CodePatterns/Log.hpp"
 | 
|---|
 | 23 | #include "CodePatterns/Verbose.hpp"
 | 
|---|
| [22c44bf] | 24 | #include "LinearAlgebra/Line.hpp"
 | 
|---|
 | 25 | #include "LinearAlgebra/Vector.hpp"
 | 
|---|
 | 26 | #include "molecule.hpp"
 | 
|---|
 | 27 | 
 | 
|---|
 | 28 | 
 | 
|---|
| [146cff2] | 29 | #include <cmath>
 | 
|---|
| [22c44bf] | 30 | #include <iostream>
 | 
|---|
 | 31 | #include <fstream>
 | 
|---|
 | 32 | #include <string>
 | 
|---|
 | 33 | 
 | 
|---|
| [1fd675] | 34 | #include "Actions/AtomAction/RotateAroundOriginByAngleAction.hpp"
 | 
|---|
| [22c44bf] | 35 | 
 | 
|---|
| [ce7fdc] | 36 | using namespace MoleCuilder;
 | 
|---|
 | 37 | 
 | 
|---|
| [1fd675] | 38 | // and construct the stuff
 | 
|---|
 | 39 | #include "RotateAroundOriginByAngleAction.def"
 | 
|---|
 | 40 | #include "Action_impl_pre.hpp"
 | 
|---|
 | 41 | /** =========== define the function ====================== */
 | 
|---|
| [22c44bf] | 42 | Action::state_ptr AtomRotateAroundOriginByAngleAction::performCall() {
 | 
|---|
 | 43 |   std::vector<atom *> selectedAtoms = World::getInstance().getSelectedAtoms();
 | 
|---|
 | 44 | 
 | 
|---|
 | 45 |   // check whether Axis is valid
 | 
|---|
| [1fd675] | 46 |   if (params.Axis.IsZero())
 | 
|---|
| [22c44bf] | 47 |     return Action::failure;
 | 
|---|
 | 48 | 
 | 
|---|
 | 49 |   // convert from degrees to radian
 | 
|---|
| [1fd675] | 50 |   params.angle *= M_PI/180.;
 | 
|---|
| [22c44bf] | 51 | 
 | 
|---|
 | 52 |   // Creation Line that is the rotation axis
 | 
|---|
| [1fd675] | 53 |   Line RotationAxis(Vector(0.,0.,0.), params.Axis);
 | 
|---|
| [22c44bf] | 54 | 
 | 
|---|
| [47d041] | 55 |   LOG(0, "Rotate around origin by " << params.angle << " radian, axis from origin to " << params.Axis << ".");
 | 
|---|
| [22c44bf] | 56 |   // TODO: use AtomSet::rotate?
 | 
|---|
 | 57 |   for (std::vector<atom *>::iterator iter = selectedAtoms.begin(); iter != selectedAtoms.end(); ++iter) {
 | 
|---|
| [1fd675] | 58 |     (*iter)->setPosition(RotationAxis.rotateVector((*iter)->getPosition(), params.angle));
 | 
|---|
| [22c44bf] | 59 |   }
 | 
|---|
| [47d041] | 60 |   LOG(0, "done.");
 | 
|---|
| [1fd675] | 61 |   return Action::state_ptr(new AtomRotateAroundOriginByAngleState(World::getInstance().getSelectedAtoms(), params));
 | 
|---|
| [22c44bf] | 62 | }
 | 
|---|
 | 63 | 
 | 
|---|
 | 64 | Action::state_ptr AtomRotateAroundOriginByAngleAction::performUndo(Action::state_ptr _state) {
 | 
|---|
 | 65 |   AtomRotateAroundOriginByAngleState *state = assert_cast<AtomRotateAroundOriginByAngleState*>(_state.get());
 | 
|---|
 | 66 | 
 | 
|---|
 | 67 |   // Creation Line that is the rotation axis
 | 
|---|
| [1fd675] | 68 |   Line RotationAxis(Vector(0.,0.,0.), state->params.Axis);
 | 
|---|
| [22c44bf] | 69 | 
 | 
|---|
 | 70 |   for (std::vector<atom *>::iterator iter = state->selectedAtoms.begin(); iter != state->selectedAtoms.end(); ++iter) {
 | 
|---|
| [1fd675] | 71 |     (*iter)->setPosition(RotationAxis.rotateVector((*iter)->getPosition(), -state->params.angle));
 | 
|---|
| [22c44bf] | 72 |   }
 | 
|---|
 | 73 | 
 | 
|---|
 | 74 |   return Action::state_ptr(_state);
 | 
|---|
 | 75 | }
 | 
|---|
 | 76 | 
 | 
|---|
 | 77 | Action::state_ptr AtomRotateAroundOriginByAngleAction::performRedo(Action::state_ptr _state){
 | 
|---|
 | 78 |   AtomRotateAroundOriginByAngleState *state = assert_cast<AtomRotateAroundOriginByAngleState*>(_state.get());
 | 
|---|
 | 79 | 
 | 
|---|
 | 80 |   // Creation Line that is the rotation axis
 | 
|---|
| [1fd675] | 81 |   Line RotationAxis(Vector(0.,0.,0.), state->params.Axis);
 | 
|---|
| [22c44bf] | 82 | 
 | 
|---|
 | 83 |   for (std::vector<atom *>::iterator iter = state->selectedAtoms.begin(); iter != state->selectedAtoms.end(); ++iter) {
 | 
|---|
| [1fd675] | 84 |     (*iter)->setPosition(RotationAxis.rotateVector((*iter)->getPosition(), state->params.angle));
 | 
|---|
| [22c44bf] | 85 |   }
 | 
|---|
 | 86 | 
 | 
|---|
 | 87 |   return Action::state_ptr(_state);
 | 
|---|
 | 88 | }
 | 
|---|
 | 89 | 
 | 
|---|
 | 90 | bool AtomRotateAroundOriginByAngleAction::canUndo() {
 | 
|---|
 | 91 |   return true;
 | 
|---|
 | 92 | }
 | 
|---|
 | 93 | 
 | 
|---|
 | 94 | bool AtomRotateAroundOriginByAngleAction::shouldUndo() {
 | 
|---|
 | 95 |   return true;
 | 
|---|
 | 96 | }
 | 
|---|
| [1fd675] | 97 | /** =========== end of function ====================== */
 | 
|---|