| 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 | * RotateAroundOriginByAngleAction.cpp | 
|---|
| 10 | * | 
|---|
| 11 | *  Created on: Aug 06, 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 "CodePatterns/Verbose.hpp" | 
|---|
| 24 | #include "LinearAlgebra/Line.hpp" | 
|---|
| 25 | #include "LinearAlgebra/Vector.hpp" | 
|---|
| 26 | #include "molecule.hpp" | 
|---|
| 27 |  | 
|---|
| 28 |  | 
|---|
| 29 | #include <cmath> | 
|---|
| 30 | #include <iostream> | 
|---|
| 31 | #include <fstream> | 
|---|
| 32 | #include <string> | 
|---|
| 33 |  | 
|---|
| 34 | using namespace std; | 
|---|
| 35 |  | 
|---|
| 36 | #include "Actions/AtomAction/RotateAroundOriginByAngleAction.hpp" | 
|---|
| 37 |  | 
|---|
| 38 | // and construct the stuff | 
|---|
| 39 | #include "RotateAroundOriginByAngleAction.def" | 
|---|
| 40 | #include "Action_impl_pre.hpp" | 
|---|
| 41 | /** =========== define the function ====================== */ | 
|---|
| 42 | Action::state_ptr AtomRotateAroundOriginByAngleAction::performCall() { | 
|---|
| 43 | std::vector<atom *> selectedAtoms = World::getInstance().getSelectedAtoms(); | 
|---|
| 44 |  | 
|---|
| 45 | // obtain information | 
|---|
| 46 | getParametersfromValueStorage(); | 
|---|
| 47 |  | 
|---|
| 48 | // check whether Axis is valid | 
|---|
| 49 | if (params.Axis.IsZero()) | 
|---|
| 50 | return Action::failure; | 
|---|
| 51 |  | 
|---|
| 52 | // convert from degrees to radian | 
|---|
| 53 | params.angle *= M_PI/180.; | 
|---|
| 54 |  | 
|---|
| 55 | // Creation Line that is the rotation axis | 
|---|
| 56 | Line RotationAxis(Vector(0.,0.,0.), params.Axis); | 
|---|
| 57 |  | 
|---|
| 58 | DoLog(0) && (Log() << Verbose(0) << "Rotate around origin by " << params.angle << " radian, axis from origin to " << params.Axis << "." << endl); | 
|---|
| 59 | // TODO: use AtomSet::rotate? | 
|---|
| 60 | for (std::vector<atom *>::iterator iter = selectedAtoms.begin(); iter != selectedAtoms.end(); ++iter) { | 
|---|
| 61 | (*iter)->setPosition(RotationAxis.rotateVector((*iter)->getPosition(), params.angle)); | 
|---|
| 62 | } | 
|---|
| 63 | DoLog(0) && (Log() << Verbose(0) << "done." << endl); | 
|---|
| 64 | return Action::state_ptr(new AtomRotateAroundOriginByAngleState(World::getInstance().getSelectedAtoms(), params)); | 
|---|
| 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 | 
|---|
| 71 | Line RotationAxis(Vector(0.,0.,0.), state->params.Axis); | 
|---|
| 72 |  | 
|---|
| 73 | for (std::vector<atom *>::iterator iter = state->selectedAtoms.begin(); iter != state->selectedAtoms.end(); ++iter) { | 
|---|
| 74 | (*iter)->setPosition(RotationAxis.rotateVector((*iter)->getPosition(), -state->params.angle)); | 
|---|
| 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 | 
|---|
| 84 | Line RotationAxis(Vector(0.,0.,0.), state->params.Axis); | 
|---|
| 85 |  | 
|---|
| 86 | for (std::vector<atom *>::iterator iter = state->selectedAtoms.begin(); iter != state->selectedAtoms.end(); ++iter) { | 
|---|
| 87 | (*iter)->setPosition(RotationAxis.rotateVector((*iter)->getPosition(), state->params.angle)); | 
|---|
| 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 | } | 
|---|
| 100 | /** =========== end of function ====================== */ | 
|---|