| 372 | 
  | 
      return fortranInitialized_; | 
| 373 | 
  | 
    } | 
| 374 | 
  | 
         | 
| 375 | 
+ | 
    bool getCalcBoxDipole() { | 
| 376 | 
+ | 
      return calcBoxDipole_; | 
| 377 | 
+ | 
    } | 
| 378 | 
+ | 
 | 
| 379 | 
  | 
    //below functions are just forward functions | 
| 380 | 
  | 
    //To compose or to inherit is always a hot debate. In general, is-a relation need subclassing, in the | 
| 381 | 
  | 
    //the other hand, has-a relation need composing. | 
| 455 | 
  | 
    /** Figure out which polynomial type to use for the switching function */ | 
| 456 | 
  | 
    void setupSwitchingFunction(); | 
| 457 | 
  | 
 | 
| 458 | 
+ | 
    /** Determine if we need to accumulate the simulation box dipole */ | 
| 459 | 
+ | 
    void setupAccumulateBoxDipole(); | 
| 460 | 
+ | 
 | 
| 461 | 
  | 
    /** Calculates the number of degress of freedom in the whole system */ | 
| 462 | 
  | 
    void calcNdf(); | 
| 463 | 
  | 
    void calcNdfRaw(); | 
| 541 | 
  | 
 | 
| 542 | 
  | 
    bool fortranInitialized_; /**< flag indicate whether fortran side is initialized */ | 
| 543 | 
  | 
 | 
| 544 | 
+ | 
    bool calcBoxDipole_; /**< flag to indicate whether or not we calculate the simulation box dipole moment */ | 
| 545 | 
+ | 
 | 
| 546 | 
  | 
#ifdef IS_MPI | 
| 547 | 
  | 
    //in Parallel version, we need MolToProc | 
| 548 | 
  | 
  public: |