Changeset c8165c for src/Actions

Timestamp:

May 8, 2017, 2:01:48 PM (8 years ago)

Author:

Frederik Heber <frederik.heber@…>

Branches:

ForceAnnealing_goodresults, ForceAnnealing_tocheck

Children:

db0833

Parents:

73faf4 (diff), c1446cd (diff)
Note: this is a merge changeset, the changes displayed below correspond to the merge itself.
Use the (diff) links above to see all the changes relative to each parent.

Message:

Merge branch 'Docu_Python_wait' into Candidate_v1.6.1

Conflicts:

tests/Python/AllActions/options.dat

• two options introduced at same place, both get in.

Location:

src/Actions

Files:

• Action.hpp (modified) (1 diff)
• GlobalListOfActions.hpp (modified) (1 diff)
• Makefile.am (modified) (3 diffs)
• MoleculeAction/RotateAroundBondAction.cpp (added)
• MoleculeAction/RotateAroundBondAction.def (added)
• MoleculeAction/RotateAroundBondAction.hpp (added)
• MoleculeAction/RotateToPrincipalAxisSystemAction.cpp (modified) (2 diffs)
• MoleculeAction/RotateToPrincipalAxisSystemAction.def (modified) (2 diffs)
• PotentialAction/ParsePotentialsAction.cpp (modified) (1 diff)
• WorldAction/StepWorldTimeAction.cpp (modified) (3 diffs)
• WorldAction/StepWorldTimeAction.def (modified) (1 diff)

src/Actions/Action.hpp

@@ -29 +29 @@
  */
 #ifndef STATUS
-#define STATUS(msg) pushStatus(msg)
+#define STATUS(msg) \
+  pushStatus(msg); \
+  LOG(0, "STATUS: " << msg)
 #endif
 

src/Actions/GlobalListOfActions.hpp

@@ -88 +88 @@
   (MoleculeLoad) \
   (MoleculeRemove) \
+  (MoleculeRotateAroundBond) \
   (MoleculeRotateAroundSelfByAngle) \
   (MoleculeRotateToPrincipalAxisSystem) \

src/Actions/Makefile.am

@@ -339 +339 @@
   Actions/MoleculeAction/LoadAction.cpp \
   Actions/MoleculeAction/RemoveAction.cpp \
+  Actions/MoleculeAction/RotateAroundBondAction.cpp \
   Actions/MoleculeAction/RotateAroundSelfByAngleAction.cpp \
   Actions/MoleculeAction/RotateToPrincipalAxisSystemAction.cpp \
@@ -357 +358 @@
   Actions/MoleculeAction/LoadAction.hpp \
   Actions/MoleculeAction/RemoveAction.hpp \
+  Actions/MoleculeAction/RotateAroundBondAction.hpp \
   Actions/MoleculeAction/RotateAroundSelfByAngleAction.hpp \
   Actions/MoleculeAction/RotateToPrincipalAxisSystemAction.hpp \
@@ -375 +377 @@
   Actions/MoleculeAction/LoadAction.def \
   Actions/MoleculeAction/RemoveAction.def \
+  Actions/MoleculeAction/RotateAroundBondAction.def \
   Actions/MoleculeAction/RotateAroundSelfByAngleAction.def \
   Actions/MoleculeAction/RotateToPrincipalAxisSystemAction.def \

src/Actions/MoleculeAction/RotateToPrincipalAxisSystemAction.cpp

@@ -35 +35 @@
 //#include "CodePatterns/MemDebug.hpp"
 
+#include "CodePatterns/Assert.hpp"
 #include "CodePatterns/Log.hpp"
-#include "CodePatterns/Verbose.hpp"
+
+#include "Actions/UndoRedoHelpers.hpp"
+#include "Atom/AtomicInfo.hpp"
 #include "LinearAlgebra/Line.hpp"
 #include "LinearAlgebra/RealSpaceMatrix.hpp"
@@ -58 +61 @@
   molecule *mol = NULL;
 
+  std::vector<AtomicInfo> UndoInfo;
   for (World::MoleculeSelectionIterator iter = World::getInstance().beginMoleculeSelection(); iter != World::getInstance().endMoleculeSelection(); ++iter) {
     mol = iter->second;
     LOG(0, "Converting to prinicipal axis system.");
 
-    RealSpaceMatrix InertiaTensor = mol->getInertiaTensor();
+    // gather undo information: store position of all atoms of molecule
+    UndoInfo.reserve(UndoInfo.size()+mol->size());
+    {
+      for (molecule::const_iterator iter = const_cast<molecule const *>(mol)->begin();
+          iter != const_cast<molecule const *>(mol)->end();
+          ++iter) {
+        const atom * const Walker = *iter;
+        UndoInfo.push_back(AtomicInfo(*Walker));
+      }
+    }
 
+    // rotate
+//    RealSpaceMatrix InertiaTensor = mol->getInertiaTensor();
     mol->RotateToPrincipalAxisSystem(params.Axis.get());
 
     // summing anew for debugging (resulting matrix has to be diagonal!)
-    InertiaTensor = mol->getInertiaTensor();
+    RealSpaceMatrix InertiaTensor = mol->getInertiaTensor();
   }
-  return Action::success;
+
+  MoleculeRotateToPrincipalAxisSystemState *UndoState =
+      new MoleculeRotateToPrincipalAxisSystemState(UndoInfo, params);
+  return ActionState::ptr(UndoState);
 }
 
 ActionState::ptr MoleculeRotateToPrincipalAxisSystemAction::performUndo(ActionState::ptr _state) {
-//  MoleculeRotateToPrincipalAxisSystemState *state = assert_cast<MoleculeRotateToPrincipalAxisSystemState*>(_state.get());
+  MoleculeRotateToPrincipalAxisSystemState *state =
+      assert_cast<MoleculeRotateToPrincipalAxisSystemState*>(_state.get());
 
-//  string newName = state->mol->getName();
-//  state->mol->setName(state->lastName);
+  // set stored old state
+  SetAtomsFromAtomicInfo(state->undoinfo);
 
-  STATUS("Undo of MoleculeRotateToPrincipalAxisSystemAction not implemented.");
-  return Action::failure;
+  return ActionState::ptr(_state);
 }
 
 ActionState::ptr MoleculeRotateToPrincipalAxisSystemAction::performRedo(ActionState::ptr _state){
-  STATUS("Redo of MoleculeRotateToPrincipalAxisSystemAction not implemented.");
-  return Action::failure;
+  MoleculeRotateToPrincipalAxisSystemState *state =
+      assert_cast<MoleculeRotateToPrincipalAxisSystemState*>(_state.get());
+
+  for (World::MoleculeSelectionIterator iter = World::getInstance().beginMoleculeSelection();
+      iter != World::getInstance().endMoleculeSelection(); ++iter) {
+    molecule * const mol = iter->second;
+    mol->RotateToPrincipalAxisSystem(state->params.Axis.get());
+  }
+
+  return ActionState::ptr(_state);
 }
 

src/Actions/MoleculeAction/RotateToPrincipalAxisSystemAction.def

@@ -7 +7 @@
 
 // all includes and forward declarations necessary for non-integral types below
+#include <vector>
+
 #include "Actions/Values.hpp"
+#include "Atom/AtomicInfo.hpp"
 #include "LinearAlgebra/Vector.hpp"
 
@@ -23 +26 @@
 (VectorNotZeroValidator())
 
-#undef statetypes
-#undef statereferences
+#define statetypes (std::vector<AtomicInfo>)
+#define statereferences (undoinfo)
 
 // some defines for all the names, you may use ACTION, STATE and PARAMS

src/Actions/PotentialAction/ParsePotentialsAction.cpp

@@ -73 +73 @@
         deserialize();
       } catch (SerializerMissingValueException &e) {
-        if (const std::string *key = boost::get_error_info<SerializerKey>(e))
+        if (const std::string *key = boost::get_error_info<SerializerKey>(e)) {
           STATUS("Missing value when parsing information for potential "+*key+".");
-        else
+        } else
           STATUS("Missing value parsing information for potential with unknown key.");
         return Action::failure;
       } catch (SerializerIllegalKeyException &e) {
-        if (const std::string *key = boost::get_error_info<SerializerKey>(e))
+        if (const std::string *key = boost::get_error_info<SerializerKey>(e)) {
           STATUS("Illegal key parsing information for potential "+*key+".");
-        else
+        } else {
           STATUS("Illegal key parsing information for potential with unknown key.");
+        }
         return Action::failure;
       }

src/Actions/WorldAction/StepWorldTimeAction.cpp

@@ -47 +47 @@
 using namespace MoleCuilder;
 
-// if both are given, we use backwards
-static int getSteps(const unsigned int _forward, const unsigned int _backward)
-{
-  int steps = 0;
-  if (_backward > 0)
-    steps = -_backward;
-  else
-    steps = _forward;
-
-  return steps;
-}
-
 // and construct the stuff
 #include "StepWorldTimeAction.def"
@@ -69 +57 @@
   WorldStepWorldTimeState *UndoState = new WorldStepWorldTimeState(oldtime, params);
 
-  const int steps = getSteps(params.steps_forward.get(), params.steps_backward.get());
-
-  if ((oldtime + steps) < 0) {
+  if ((oldtime + params.steps.get()) < 0) {
     ELOG(1, "Cannot step back before time step #0.");
     delete UndoState;
     return Action::failure;
   }
-  World::getInstance().setTime(oldtime+steps);
+  World::getInstance().setTime(oldtime+params.steps.get());
   LOG(0, "Current time step is now: " << WorldTime::getTime() << ".");
   return ActionState::ptr(UndoState);
@@ -93 +79 @@
   WorldStepWorldTimeState *state = assert_cast<WorldStepWorldTimeState*>(_state.get());
 
-  const int steps = getSteps(
-      state->params.steps_forward.get(),
-      state->params.steps_backward.get());
-
-  World::getInstance().setTime(state->oldtime+steps);
+  World::getInstance().setTime(state->oldtime+state->params.steps.get());
   LOG(0, "Current time step is now: " << WorldTime::getTime() << ".");
 

src/Actions/WorldAction/StepWorldTimeAction.def

@@ -14 +14 @@
 // ValueStorage by the token "Z" -> first column: int, Z, "Z"
 // "undefine" if no parameters are required, use (NOPARAM_DEFAULT) for each (undefined) default value
-#define paramtypes (unsigned int)(unsigned int)
-#define paramtokens ("steps-forward")("steps-backward")
-#define paramdescriptions ("how many steps to take forward or backward")("how many steps to take forward or backward")
-#define paramdefaults (PARAM_DEFAULT(1))(PARAM_DEFAULT(0))
-#define paramreferences (steps_forward)(steps_backward)
+#define paramtypes (int)
+#define paramtokens ("step-world-time")
+#define paramdescriptions ("how many steps to take forward (positive) or backward (negative)")
+#define paramdefaults (PARAM_DEFAULT(1))
+#define paramreferences (steps)
 #define paramvalids \
-(DummyValidator< unsigned int >()) \
-(DummyValidator< unsigned int >())
+(DummyValidator< int >())
 
 #define statetypes (unsigned int)