| 279 |
|
zPos.push_back((*parameters)[searchResult].zPos); |
| 280 |
|
// cout << "index: "<< (*parameters)[searchResult].zconsIndex |
| 281 |
|
// <<"\tzPos = " << (*parameters)[searchResult].zPos << endl; |
| 282 |
– |
kz.push_back((*parameters)[searchResult]. kRatio * zForceConst); |
| 282 |
|
|
| 283 |
+ |
kz.push_back((*parameters)[searchResult]. kRatio * zForceConst); |
| 284 |
|
molecules[i].getCOM(COM); |
| 285 |
|
} |
| 286 |
|
else |
| 394 |
|
zconsMols.push_back(&molecules[i]); |
| 395 |
|
zPos.push_back((*parameters)[index].zPos); |
| 396 |
|
kz.push_back((*parameters)[index].kRatio * zForceConst); |
| 397 |
– |
|
| 397 |
|
massOfZConsMols.push_back(molecules[i].getTotalMass()); |
| 398 |
|
|
| 399 |
|
molecules[i].getCOM(COM); |
| 931 |
|
MPI_Allreduce(&totalFZ_local, &totalFZ, 1, MPI_DOUBLE,MPI_SUM, MPI_COMM_WORLD); |
| 932 |
|
#endif |
| 933 |
|
|
| 934 |
< |
cout << "before substracting harmonic force from moving molecuels " |
| 935 |
< |
<< "total force is " << calcTotalForce() << endl; |
| 934 |
> |
//cout << "before substracting harmonic force from moving molecuels " |
| 935 |
> |
// << "total force is " << calcTotalForce() << endl; |
| 936 |
|
|
| 937 |
|
force[0]= 0; |
| 938 |
|
force[1]= 0; |
| 950 |
|
} |
| 951 |
|
} |
| 952 |
|
|
| 953 |
< |
cout << "after substracting harmonic force from moving molecuels " |
| 954 |
< |
<< "total force is " << calcTotalForce() << endl; |
| 953 |
> |
//cout << "after substracting harmonic force from moving molecuels " |
| 954 |
> |
// << "total force is " << calcTotalForce() << endl; |
| 955 |
|
|
| 956 |
|
} |
| 957 |
|
|