| 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 |  |  |