| [bcf653] | 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 |  | 
|---|
| [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 |  | 
|---|
|  | 34 | using namespace std; | 
|---|
|  | 35 |  | 
|---|
| [1fd675] | 36 | #include "Actions/AtomAction/RotateAroundOriginByAngleAction.hpp" | 
|---|
| [22c44bf] | 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 |  | 
|---|
| [1fd675] | 45 | // obtain information | 
|---|
|  | 46 | getParametersfromValueStorage(); | 
|---|
| [22c44bf] | 47 |  | 
|---|
|  | 48 | // check whether Axis is valid | 
|---|
| [1fd675] | 49 | if (params.Axis.IsZero()) | 
|---|
| [22c44bf] | 50 | return Action::failure; | 
|---|
|  | 51 |  | 
|---|
|  | 52 | // convert from degrees to radian | 
|---|
| [1fd675] | 53 | params.angle *= M_PI/180.; | 
|---|
| [22c44bf] | 54 |  | 
|---|
|  | 55 | // Creation Line that is the rotation axis | 
|---|
| [1fd675] | 56 | Line RotationAxis(Vector(0.,0.,0.), params.Axis); | 
|---|
| [22c44bf] | 57 |  | 
|---|
| [1fd675] | 58 | DoLog(0) && (Log() << Verbose(0) << "Rotate around origin by " << params.angle << " radian, axis from origin to " << params.Axis << "." << endl); | 
|---|
| [22c44bf] | 59 | // TODO: use AtomSet::rotate? | 
|---|
|  | 60 | for (std::vector<atom *>::iterator iter = selectedAtoms.begin(); iter != selectedAtoms.end(); ++iter) { | 
|---|
| [1fd675] | 61 | (*iter)->setPosition(RotationAxis.rotateVector((*iter)->getPosition(), params.angle)); | 
|---|
| [22c44bf] | 62 | } | 
|---|
|  | 63 | DoLog(0) && (Log() << Verbose(0) << "done." << endl); | 
|---|
| [1fd675] | 64 | return Action::state_ptr(new AtomRotateAroundOriginByAngleState(World::getInstance().getSelectedAtoms(), params)); | 
|---|
| [22c44bf] | 65 | } | 
|---|
|  | 66 |  | 
|---|
|  | 67 | Action::state_ptr AtomRotateAroundOriginByAngleAction::performUndo(Action::state_ptr _state) { | 
|---|
|  | 68 | AtomRotateAroundOriginByAngleState *state = assert_cast<AtomRotateAroundOriginByAngleState*>(_state.get()); | 
|---|
|  | 69 |  | 
|---|
|  | 70 | // Creation Line that is the rotation axis | 
|---|
| [1fd675] | 71 | Line RotationAxis(Vector(0.,0.,0.), state->params.Axis); | 
|---|
| [22c44bf] | 72 |  | 
|---|
|  | 73 | for (std::vector<atom *>::iterator iter = state->selectedAtoms.begin(); iter != state->selectedAtoms.end(); ++iter) { | 
|---|
| [1fd675] | 74 | (*iter)->setPosition(RotationAxis.rotateVector((*iter)->getPosition(), -state->params.angle)); | 
|---|
| [22c44bf] | 75 | } | 
|---|
|  | 76 |  | 
|---|
|  | 77 | return Action::state_ptr(_state); | 
|---|
|  | 78 | } | 
|---|
|  | 79 |  | 
|---|
|  | 80 | Action::state_ptr AtomRotateAroundOriginByAngleAction::performRedo(Action::state_ptr _state){ | 
|---|
|  | 81 | AtomRotateAroundOriginByAngleState *state = assert_cast<AtomRotateAroundOriginByAngleState*>(_state.get()); | 
|---|
|  | 82 |  | 
|---|
|  | 83 | // Creation Line that is the rotation axis | 
|---|
| [1fd675] | 84 | Line RotationAxis(Vector(0.,0.,0.), state->params.Axis); | 
|---|
| [22c44bf] | 85 |  | 
|---|
|  | 86 | for (std::vector<atom *>::iterator iter = state->selectedAtoms.begin(); iter != state->selectedAtoms.end(); ++iter) { | 
|---|
| [1fd675] | 87 | (*iter)->setPosition(RotationAxis.rotateVector((*iter)->getPosition(), state->params.angle)); | 
|---|
| [22c44bf] | 88 | } | 
|---|
|  | 89 |  | 
|---|
|  | 90 | return Action::state_ptr(_state); | 
|---|
|  | 91 | } | 
|---|
|  | 92 |  | 
|---|
|  | 93 | bool AtomRotateAroundOriginByAngleAction::canUndo() { | 
|---|
|  | 94 | return true; | 
|---|
|  | 95 | } | 
|---|
|  | 96 |  | 
|---|
|  | 97 | bool AtomRotateAroundOriginByAngleAction::shouldUndo() { | 
|---|
|  | 98 | return true; | 
|---|
|  | 99 | } | 
|---|
| [1fd675] | 100 | /** =========== end of function ====================== */ | 
|---|