6 |
|
* redistribute this software in source and binary code form, provided |
7 |
|
* that the following conditions are met: |
8 |
|
* |
9 |
< |
* 1. Acknowledgement of the program authors must be made in any |
10 |
< |
* publication of scientific results based in part on use of the |
11 |
< |
* program. An acceptable form of acknowledgement is citation of |
12 |
< |
* the article in which the program was described (Matthew |
13 |
< |
* A. Meineke, Charles F. Vardeman II, Teng Lin, Christopher |
14 |
< |
* J. Fennell and J. Daniel Gezelter, "OOPSE: An Object-Oriented |
15 |
< |
* Parallel Simulation Engine for Molecular Dynamics," |
16 |
< |
* J. Comput. Chem. 26, pp. 252-271 (2005)) |
17 |
< |
* |
18 |
< |
* 2. Redistributions of source code must retain the above copyright |
9 |
> |
* 1. Redistributions of source code must retain the above copyright |
10 |
|
* notice, this list of conditions and the following disclaimer. |
11 |
|
* |
12 |
< |
* 3. Redistributions in binary form must reproduce the above copyright |
12 |
> |
* 2. Redistributions in binary form must reproduce the above copyright |
13 |
|
* notice, this list of conditions and the following disclaimer in the |
14 |
|
* documentation and/or other materials provided with the |
15 |
|
* distribution. |
28 |
|
* arising out of the use of or inability to use software, even if the |
29 |
|
* University of Notre Dame has been advised of the possibility of |
30 |
|
* such damages. |
31 |
+ |
* |
32 |
+ |
* SUPPORT OPEN SCIENCE! If you use OpenMD or its source code in your |
33 |
+ |
* research, please cite the appropriate papers when you publish your |
34 |
+ |
* work. Good starting points are: |
35 |
+ |
* |
36 |
+ |
* [1] Meineke, et al., J. Comp. Chem. 26, 252-271 (2005). |
37 |
+ |
* [2] Fennell & Gezelter, J. Chem. Phys. 124, 234104 (2006). |
38 |
+ |
* [3] Sun, Lin & Gezelter, J. Chem. Phys. 128, 234107 (2008). |
39 |
+ |
* [4] Kuang & Gezelter, J. Chem. Phys. 133, 164101 (2010). |
40 |
+ |
* [5] Vardeman, Stocker & Gezelter, J. Chem. Theory Comput. 7, 834 (2011). |
41 |
|
*/ |
42 |
|
|
43 |
|
/** |
55 |
|
#include <list> |
56 |
|
#include <utility> |
57 |
|
|
58 |
< |
namespace oopse { |
58 |
> |
namespace OpenMD { |
59 |
|
|
60 |
|
/** |
61 |
|
* @class IndexListContainer |
65 |
|
class IndexListContainer{ |
66 |
|
public: |
67 |
|
static const int MAX_INTEGER = 2147483647; |
68 |
< |
typedef std::list<std::pair<int, int> >::iterator IndexListContainerIterator; |
68 |
> |
typedef std::list<std::pair<int, int> >::iterator IndexListContainerIterator; |
69 |
|
|
69 |
– |
|
70 |
|
IndexListContainer(int minIndex = 0, int maxIndex = MAX_INTEGER) |
71 |
|
:minIndex_(minIndex), maxIndex_(maxIndex) { |
72 |
< |
|
73 |
< |
indexContainer_.push_back(std::make_pair(minIndex, maxIndex)); |
74 |
< |
} |
72 |
> |
indexContainer_.push_back(std::make_pair(minIndex, maxIndex)); |
73 |
> |
} |
74 |
|
|
75 |
< |
int pop() { |
77 |
< |
|
75 |
> |
int pop() { |
76 |
|
if (indexContainer_.empty()) { |
77 |
|
std::cerr << "" << std::endl; |
78 |
|
} |
79 |
|
|
82 |
– |
IndexListContainerIterator i = indexContainer_.begin(); |
80 |
|
int result; |
81 |
|
|
82 |
|
result = indexContainer_.front().first; |
83 |
|
|
84 |
|
if (indexContainer_.front().first == indexContainer_.front().second) { |
85 |
|
indexContainer_.pop_front(); |
86 |
< |
} else if (indexContainer_.front().first < indexContainer_.front().second) { |
86 |
> |
} else if (indexContainer_.front().first < |
87 |
> |
indexContainer_.front().second) { |
88 |
|
++indexContainer_.front().first; |
89 |
|
} else { |
90 |
|
std::cerr << "" << std::endl; |
91 |
|
} |
92 |
< |
|
92 |
> |
|
93 |
|
return result; |
94 |
|
} |
95 |
|
|
98 |
|
* |
99 |
|
*/ |
100 |
|
void insert(int index) { |
101 |
< |
IndexListContainerIterator insertPos = internalInsert(index, index); |
101 |
> |
IndexListContainerIterator insertPos = internalInsert(index, index); |
102 |
|
merge(insertPos); |
103 |
|
} |
104 |
|
|
108 |
|
* @param endIndex |
109 |
|
*/ |
110 |
|
void insert(int beginIndex, int endIndex) { |
111 |
< |
IndexListContainerIterator insertPos = internalInsert(beginIndex, endIndex); |
111 |
> |
IndexListContainerIterator insertPos = internalInsert(beginIndex, |
112 |
> |
endIndex); |
113 |
|
merge(insertPos); |
114 |
|
} |
115 |
|
|
127 |
|
std::unique(indices.begin(), indices.end()); |
128 |
|
|
129 |
|
std::vector<int>::iterator i; |
131 |
– |
IndexListContainerIterator insertPos; |
130 |
|
int beginIndex; |
131 |
|
|
132 |
|
beginIndex = indices[0]; |
137 |
|
beginIndex = *i; |
138 |
|
} |
139 |
|
} |
142 |
– |
|
143 |
– |
|
140 |
|
} |
141 |
|
|
142 |
|
std::vector<int> getIndicesBefore(int index) { |
160 |
|
indexContainer_.erase(indexContainer_.begin(), i); |
161 |
|
break; |
162 |
|
} else { |
167 |
– |
|
163 |
|
for (int j = (*i).first; j < index; ++j) { |
164 |
|
result.push_back(j); |
165 |
< |
} |
171 |
< |
|
165 |
> |
} |
166 |
|
(*i).first = index; |
167 |
|
indexContainer_.erase(indexContainer_.begin(), i); |
168 |
|
break; |
169 |
|
} |
176 |
– |
|
170 |
|
} |
178 |
– |
|
171 |
|
return result; |
180 |
– |
|
172 |
|
} |
173 |
|
|
174 |
|
int getMaxIndex() { |
187 |
|
if (endIndex > maxIndex_) { |
188 |
|
std::cerr << "" << std::endl; |
189 |
|
} |
190 |
< |
|
200 |
< |
|
201 |
< |
IndexListContainerIterator j; |
202 |
< |
|
190 |
> |
|
191 |
|
IndexListContainerIterator i = indexContainer_.begin(); |
192 |
|
for (; i != indexContainer_.end(); ++i) { |
193 |
|
if ((*i).first > endIndex) { |
282 |
|
void releaseRigidBodyIndex(std::vector<int> indices) { |
283 |
|
rigidBodyIndexContainer_.insert(indices); |
284 |
|
} |
285 |
< |
|
285 |
> |
|
286 |
> |
int getNextCutoffGroupIndex() { |
287 |
> |
return cutoffGroupIndexContainer_.pop(); |
288 |
> |
} |
289 |
> |
|
290 |
> |
std::vector<int> getCutoffGroupIndicesBefore(int index) { |
291 |
> |
return cutoffGroupIndexContainer_.getIndicesBefore(index); |
292 |
> |
} |
293 |
> |
|
294 |
> |
void releaseCutoffGroupIndex(int index) { |
295 |
> |
cutoffGroupIndexContainer_.insert(index); |
296 |
> |
} |
297 |
> |
|
298 |
> |
void releaseCutoffGroupIndex(int beginIndex, int endIndex) { |
299 |
> |
cutoffGroupIndexContainer_.insert(beginIndex, endIndex); |
300 |
> |
} |
301 |
> |
|
302 |
> |
void releaseCutoffGroupIndex(std::vector<int> indices) { |
303 |
> |
cutoffGroupIndexContainer_.insert(indices); |
304 |
> |
} |
305 |
> |
|
306 |
> |
int getNextBondIndex() { |
307 |
> |
return bondIndexContainer_.pop(); |
308 |
> |
} |
309 |
> |
|
310 |
> |
std::vector<int> getBondIndicesBefore(int index) { |
311 |
> |
return bondIndexContainer_.getIndicesBefore(index); |
312 |
> |
} |
313 |
> |
|
314 |
> |
void releaseBondIndex(int index) { |
315 |
> |
bondIndexContainer_.insert(index); |
316 |
> |
} |
317 |
> |
|
318 |
> |
void releaseBondIndex(int beginIndex, int endIndex) { |
319 |
> |
bondIndexContainer_.insert(beginIndex, endIndex); |
320 |
> |
} |
321 |
> |
|
322 |
> |
void releaseBondIndex(std::vector<int> indices) { |
323 |
> |
bondIndexContainer_.insert(indices); |
324 |
> |
} |
325 |
> |
|
326 |
> |
int getNextBendIndex() { |
327 |
> |
return bendIndexContainer_.pop(); |
328 |
> |
} |
329 |
> |
|
330 |
> |
std::vector<int> getBendIndicesBefore(int index) { |
331 |
> |
return bendIndexContainer_.getIndicesBefore(index); |
332 |
> |
} |
333 |
> |
|
334 |
> |
void releaseBendIndex(int index) { |
335 |
> |
bendIndexContainer_.insert(index); |
336 |
> |
} |
337 |
> |
|
338 |
> |
void releaseBendIndex(int beginIndex, int endIndex) { |
339 |
> |
bendIndexContainer_.insert(beginIndex, endIndex); |
340 |
> |
} |
341 |
> |
|
342 |
> |
void releaseBendIndex(std::vector<int> indices) { |
343 |
> |
bendIndexContainer_.insert(indices); |
344 |
> |
} |
345 |
> |
|
346 |
> |
int getNextTorsionIndex() { |
347 |
> |
return torsionIndexContainer_.pop(); |
348 |
> |
} |
349 |
> |
|
350 |
> |
std::vector<int> getTorsionIndicesBefore(int index) { |
351 |
> |
return torsionIndexContainer_.getIndicesBefore(index); |
352 |
> |
} |
353 |
> |
|
354 |
> |
void releaseTorsionIndex(int index) { |
355 |
> |
torsionIndexContainer_.insert(index); |
356 |
> |
} |
357 |
> |
|
358 |
> |
void releaseTorsionIndex(int beginIndex, int endIndex) { |
359 |
> |
torsionIndexContainer_.insert(beginIndex, endIndex); |
360 |
> |
} |
361 |
> |
|
362 |
> |
void releaseTorsionIndex(std::vector<int> indices) { |
363 |
> |
torsionIndexContainer_.insert(indices); |
364 |
> |
} |
365 |
> |
|
366 |
> |
int getNextInversionIndex() { |
367 |
> |
return inversionIndexContainer_.pop(); |
368 |
> |
} |
369 |
> |
|
370 |
> |
std::vector<int> getInversionIndicesBefore(int index) { |
371 |
> |
return inversionIndexContainer_.getIndicesBefore(index); |
372 |
> |
} |
373 |
> |
|
374 |
> |
void releaseInversionIndex(int index) { |
375 |
> |
inversionIndexContainer_.insert(index); |
376 |
> |
} |
377 |
> |
|
378 |
> |
void releaseInversionIndex(int beginIndex, int endIndex) { |
379 |
> |
inversionIndexContainer_.insert(beginIndex, endIndex); |
380 |
> |
} |
381 |
> |
|
382 |
> |
void releaseInversionIndex(std::vector<int> indices) { |
383 |
> |
inversionIndexContainer_.insert(indices); |
384 |
> |
} |
385 |
> |
|
386 |
|
private: |
387 |
|
|
388 |
|
IndexListContainer atomIndexContainer_; |
389 |
|
IndexListContainer rigidBodyIndexContainer_; |
390 |
+ |
IndexListContainer cutoffGroupIndexContainer_; |
391 |
+ |
IndexListContainer bondIndexContainer_; |
392 |
+ |
IndexListContainer bendIndexContainer_; |
393 |
+ |
IndexListContainer torsionIndexContainer_; |
394 |
+ |
IndexListContainer inversionIndexContainer_; |
395 |
|
}; |
396 |
|
|
397 |
< |
} //end namespace oopse |
397 |
> |
} //end namespace OpenMD |
398 |
|
#endif //UTILS_LOCALINDEXMANAGER_HPP |