201 |
|
short int have_pos_iter_tolerance; |
202 |
|
double etaTolerance; |
203 |
|
short int have_eta_tolerance; |
204 |
– |
|
205 |
– |
double volume; |
204 |
|
|
205 |
|
}; |
206 |
|
|
318 |
|
public: |
319 |
|
|
320 |
|
NPTf ( SimInfo *theInfo, ForceFields* the_ff); |
321 |
< |
virtual ~NPTf() {}; |
321 |
> |
virtual ~NPTf(); |
322 |
|
|
323 |
|
virtual void integrateStep( int calcPot, int calcStress ){ |
324 |
|
calcStress = 1; |
325 |
|
T::integrateStep( calcPot, calcStress ); |
328 |
– |
accIntegralOfChidt(); |
326 |
|
} |
327 |
|
|
328 |
|
virtual double getConservedQuantity(void); |
343 |
|
|
344 |
|
virtual int readyCheck(); |
345 |
|
|
349 |
– |
void accIntegralOfChidt(void) { integralOfChidt += dt * chi;} |
346 |
|
|
347 |
|
// chi and eta are the propagated degrees of freedom |
348 |
|
|
349 |
|
double chi; |
350 |
|
double eta[3][3]; |
351 |
|
double NkBT; |
352 |
+ |
double fkBT; |
353 |
+ |
|
354 |
+ |
int Nparticles; |
355 |
|
|
356 |
+ |
double *oldPos; |
357 |
+ |
double *oldVel; |
358 |
+ |
double *oldJi; |
359 |
+ |
|
360 |
|
double integralOfChidt; |
361 |
|
|
362 |
|
// targetTemp, targetPressure, and tauBarostat must be set. |
373 |
|
short int have_chi_tolerance; |
374 |
|
double posIterTolerance; |
375 |
|
short int have_pos_iter_tolerance; |
376 |
< |
|
376 |
> |
double etaTolerance; |
377 |
> |
short int have_eta_tolerance; |
378 |
|
}; |
379 |
|
|
380 |
|
template<typename T> class NPTxym : public T{ |