Changeset 4eee8f for molecuilder/src/molecule_dynamics.cpp
- Timestamp:
- Apr 30, 2010, 10:07:23 AM (16 years ago)
- Children:
- 59e7832
- Parents:
- f70c2a
- File:
-
- 1 edited
-
molecuilder/src/molecule_dynamics.cpp (modified) (4 diffs)
Legend:
- Unmodified
- Added
- Removed
-
molecuilder/src/molecule_dynamics.cpp
rf70c2a r4eee8f 48 48 // check whether either is zero() 49 49 if ((Norm1 < MYEPSILON) && (Norm2 < MYEPSILON)) { 50 tmp = Walker->Trajectory.R.at(Params.startstep). Distance(Runner->Trajectory.R.at(Params.startstep));50 tmp = Walker->Trajectory.R.at(Params.startstep).distance(Runner->Trajectory.R.at(Params.startstep)); 51 51 } else if (Norm1 < MYEPSILON) { 52 52 Sprinter = Params.PermutationMap[Walker->nr]; // find first target point … … 66 66 // Log() << Verbose(0) << " and "; 67 67 // Log() << Verbose(0) << trajectory2; 68 tmp = Walker->Trajectory.R.at(Params.startstep). Distance(Runner->Trajectory.R.at(Params.startstep));68 tmp = Walker->Trajectory.R.at(Params.startstep).distance(Runner->Trajectory.R.at(Params.startstep)); 69 69 // Log() << Verbose(0) << " with distance " << tmp << "." << endl; 70 70 } else { // determine distance by finding minimum distance … … 166 166 // first term: distance to target 167 167 Runner = Params.PermutationMap[Walker->nr]; // find target point 168 tmp = (Walker->Trajectory.R.at(Params.startstep). Distance(Runner->Trajectory.R.at(Params.endstep)));168 tmp = (Walker->Trajectory.R.at(Params.startstep).distance(Runner->Trajectory.R.at(Params.endstep))); 169 169 tmp *= Params.IsAngstroem ? 1. : 1./AtomicLengthToAngstroem; 170 170 result += Params.PenaltyConstants[0] * tmp; … … 225 225 while(Runner->next != mol->end) { 226 226 Runner = Runner->next; 227 Params.DistanceList[Walker->nr]->insert( DistancePair(Walker->Trajectory.R.at(Params.startstep). Distance(Runner->Trajectory.R.at(Params.endstep)), Runner) );227 Params.DistanceList[Walker->nr]->insert( DistancePair(Walker->Trajectory.R.at(Params.startstep).distance(Runner->Trajectory.R.at(Params.endstep)), Runner) ); 228 228 } 229 229 }
Note:
See TracChangeset
for help on using the changeset viewer.
