Ignore:
File:
1 edited

Legend:

Unmodified
Added
Removed
  • src/Actions/MoleculeAction/RotateAroundSelfByAngleAction.cpp

    r833b15 rb8463f  
    6262  }
    6363
     64  // check whether Axis is valid
     65  if (params.Axis.get().IsZero()) {
     66    STATUS("Rotation Axis must not be zero.");
     67    return Action::failure;
     68  }
     69
     70  // convert from degrees to radian
     71  const double angle_radian = params.angle.get() * M_PI/180.;
     72
    6473  // go through all selected molecules
    6574  BOOST_FOREACH(molecule *mol, selectedMolecules) {
    66     // check whether Axis is valid
    67     if (params.Axis.get().IsZero()) {
    68       STATUS("Rotation Axis must not be zero.");
    69       return Action::failure;
    70     }
    71 
    72     // convert from degrees to radian
    73     params.angle.set(params.angle.get() * M_PI/180.);
    74 
    7575    // Creation Line that is the rotation axis
    7676    const Vector CenterOfGravity = mol->DetermineCenterOfGravity();
    77     LOG(0,  "Center of gravity is " << CenterOfGravity << ".");
     77    LOG(3, "DEBUG: Center of gravity is " << CenterOfGravity << ".");
    7878    Line RotationAxis(CenterOfGravity, params.Axis.get());
    79     LOG(0, "Rotate " << mol->getName() << " around self by " << params.angle.get() << " radian around axis " << RotationAxis << ".");
     79    LOG(2, "DEBUG: Rotate " << mol->getName() << " around self by " << angle_radian << " radian around axis " << RotationAxis << ".");
    8080
    8181    for (molecule::iterator iter = mol->begin(); iter != mol->end(); ++iter) {
    82       (*iter)->setPosition(RotationAxis.rotateVector((*iter)->getPosition(), params.angle.get()));
     82      (*iter)->setPosition(RotationAxis.rotateVector((*iter)->getPosition(), angle_radian));
    8383    }
    84     LOG(0, "done.");
     84    LOG(2, "DEBUG: Rotatating of " << mol->getName() << " done.");
    8585  }
    8686
     
    9191  MoleculeRotateAroundSelfByAngleState *state = assert_cast<MoleculeRotateAroundSelfByAngleState*>(_state.get());
    9292
     93  // convert from degrees to radian
     94  const double angle_radian = -state->params.angle.get() * M_PI/180.;
     95
    9396  BOOST_FOREACH(molecule *mol, state->selectedMolecules) {
    9497    const Vector CenterOfGravity = mol->DetermineCenterOfGravity();
    9598    LOG(0,  "Center of gravity is " <<CenterOfGravity << ".");
    9699    Line RotationAxis(CenterOfGravity, state->params.Axis.get());
    97     LOG(0, "Rotate " << mol->getName() << " around self by " << -state->params.angle.get() << " radian around axis " << RotationAxis << ".");
     100    LOG(0, "Rotate " << mol->getName() << " around self by " << angle_radian << " radian around axis " << RotationAxis << ".");
    98101
    99     for (molecule::iterator iter = mol->begin(); iter != mol->end(); ++iter) {
    100       (*iter)->setPosition(RotationAxis.rotateVector((*iter)->getPosition(), -state->params.angle.get()));
    101     }
     102    for (molecule::iterator iter = mol->begin(); iter != mol->end(); ++iter)
     103      (*iter)->setPosition(RotationAxis.rotateVector((*iter)->getPosition(), angle_radian));
    102104  }
    103105
     
    108110  MoleculeRotateAroundSelfByAngleState *state = assert_cast<MoleculeRotateAroundSelfByAngleState*>(_state.get());
    109111
     112  // convert from degrees to radian
     113  const double angle_radian = state->params.angle.get() * M_PI/180.;
     114
    110115  BOOST_FOREACH(molecule *mol, state->selectedMolecules) {
    111116    const Vector CenterOfGravity = mol->DetermineCenterOfGravity();
    112117    LOG(0,  "Center of gravity is " << CenterOfGravity << ".");
    113118    Line RotationAxis(CenterOfGravity, state->params.Axis.get());
    114     LOG(0, "Rotate " << mol->getName() << " around self by " << state->params.angle.get() << " radian around axis " << RotationAxis << ".");
     119    LOG(0, "Rotate " << mol->getName() << " around self by " << angle_radian << " radian around axis " << RotationAxis << ".");
    115120
    116     for (molecule::iterator iter = mol->begin(); iter != mol->end(); ++iter) {
    117       (*iter)->setPosition(RotationAxis.rotateVector((*iter)->getPosition(), state->params.angle.get()));
    118     }
     121    for (molecule::iterator iter = mol->begin(); iter != mol->end(); ++iter)
     122      (*iter)->setPosition(RotationAxis.rotateVector((*iter)->getPosition(), angle_radian));
    119123  }
    120124
Note: See TracChangeset for help on using the changeset viewer.