Changeset 5f4960 for molecuilder
- Timestamp:
- Aug 21, 2009, 12:17:25 PM (16 years ago)
- Children:
- 8b7e0a
- Parents:
- 84b811
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
molecuilder/src/molecules.cpp
r84b811 r5f4960 1168 1168 /** Evaluates the potential energy used for constrained molecular dynamics. 1169 1169 * \f$V_i^{con} = c^{bond} \cdot | r_{P(i)} - R_i | + sum_{i \neq j} C^{min} \cdot \frac{1}{C_{ij}} + C^{inj} \Bigl (1 - \theta \bigl (\prod_{i \neq j} (P(i) - P(j)) \bigr ) \Bigr )\f$ 1170 * where the first term points to the target in minimum distance, the second is a penalty for trajectories lying too close to each other (\f$C_{ij} $ is minimum distance between1170 * where the first term points to the target in minimum distance, the second is a penalty for trajectories lying too close to each other (\f$C_{ij}\f$ is minimum distance between 1171 1171 * trajectories i and j) and the third term is a penalty for two atoms trying to each the same target point. 1172 1172 * Note that for the second term we have to solve the following linear system: … … 3722 3722 * \param *&SortIndex Mapping array of size molecule::AtomCount 3723 3723 * \return true - success, false - failure of SortIndex alloc 3724 * \todo do we really need this still as the IonType may appear in any order due to recent changes 3724 3725 */ 3725 3726 bool molecule::CreateMappingLabelsToConfigSequence(ofstream *out, int *&SortIndex)
Note:
See TracChangeset
for help on using the changeset viewer.