| 111 | 
  | 
       | 
| 112 | 
  | 
      if( !strcmp( info.name, name ) ){ | 
| 113 | 
  | 
        sprintf( painCave.errMsg, | 
| 114 | 
< | 
                 "Duplicate TraPPE_Ex atom type \"%s\" found in " | 
| 115 | 
< | 
                 "the TraPPE_ExFF param file./n", | 
| 114 | 
> | 
                 "Duplicate DUFF atom type \"%s\" found in " | 
| 115 | 
> | 
                 "the DUFF param file./n", | 
| 116 | 
  | 
                 name ); | 
| 117 | 
  | 
        painCave.isFatal = 1; | 
| 118 | 
  | 
        simError(); | 
| 194 | 
  | 
       | 
| 195 | 
  | 
      if(dup){ | 
| 196 | 
  | 
        sprintf( painCave.errMsg, | 
| 197 | 
< | 
                 "Duplicate TraPPE_Ex bond type \"%s - %s\" found in " | 
| 198 | 
< | 
                 "the TraPPE_ExFF param file./n", | 
| 197 | 
> | 
                 "Duplicate DUFF bond type \"%s - %s\" found in " | 
| 198 | 
> | 
                 "the DUFF param file./n", | 
| 199 | 
  | 
                 nameA, nameB ); | 
| 200 | 
  | 
        painCave.isFatal = 1; | 
| 201 | 
  | 
        simError(); | 
| 264 | 
  | 
 | 
| 265 | 
  | 
      if(dup){ | 
| 266 | 
  | 
        sprintf( painCave.errMsg, | 
| 267 | 
< | 
                 "Duplicate TraPPE_Ex bend type \"%s - %s - %s\" found in " | 
| 268 | 
< | 
                 "the TraPPE_ExFF param file./n", | 
| 267 | 
> | 
                 "Duplicate DUFF bend type \"%s - %s - %s\" found in " | 
| 268 | 
> | 
                 "the DUFF param file./n", | 
| 269 | 
  | 
                 nameA, nameB, nameC ); | 
| 270 | 
  | 
        painCave.isFatal = 1; | 
| 271 | 
  | 
        simError(); | 
| 349 | 
  | 
       | 
| 350 | 
  | 
      if(dup){ | 
| 351 | 
  | 
        sprintf( painCave.errMsg, | 
| 352 | 
< | 
                 "Duplicate TraPPE_Ex torsion type \"%s - %s - %s - %s\" found in " | 
| 353 | 
< | 
                 "the TraPPE_ExFF param file./n", nameA, nameB, nameC, nameD ); | 
| 352 | 
> | 
                 "Duplicate DUFF torsion type \"%s - %s - %s - %s\" found in " | 
| 353 | 
> | 
                 "the DUFF param file./n", nameA, nameB, nameC, nameD ); | 
| 354 | 
  | 
        painCave.isFatal = 1; | 
| 355 | 
  | 
        simError(); | 
| 356 | 
  | 
      } | 
| 418 | 
  | 
//**************************************************************** | 
| 419 | 
  | 
 | 
| 420 | 
  | 
 | 
| 421 | 
< | 
TraPPE_ExFF::TraPPE_ExFF(){ | 
| 421 | 
> | 
DUFF::DUFF(){ | 
| 422 | 
  | 
 | 
| 423 | 
  | 
  char fileName[200]; | 
| 424 | 
  | 
  char* ffPath_env = "FORCE_PARAM_PATH"; | 
| 537 | 
  | 
     | 
| 538 | 
  | 
    // generate the force file name | 
| 539 | 
  | 
     | 
| 540 | 
< | 
    strcpy( fileName, "TraPPE_Ex.frc" );  | 
| 540 | 
> | 
    strcpy( fileName, "DUFF.frc" );  | 
| 541 | 
  | 
    //    fprintf( stderr,"Trying to open %s\n", fileName ); | 
| 542 | 
  | 
     | 
| 543 | 
  | 
    // attempt to open the file in the current directory first. | 
| 576 | 
  | 
#ifdef IS_MPI | 
| 577 | 
  | 
  } | 
| 578 | 
  | 
   | 
| 579 | 
< | 
  sprintf( checkPointMsg, "TraPPE_ExFF file opened sucessfully." ); | 
| 579 | 
> | 
  sprintf( checkPointMsg, "DUFF file opened sucessfully." ); | 
| 580 | 
  | 
  MPIcheckPoint(); | 
| 581 | 
  | 
   | 
| 582 | 
  | 
#endif // is_mpi | 
| 583 | 
  | 
} | 
| 584 | 
  | 
 | 
| 585 | 
  | 
 | 
| 586 | 
< | 
TraPPE_ExFF::~TraPPE_ExFF(){ | 
| 586 | 
> | 
DUFF::~DUFF(){ | 
| 587 | 
  | 
 | 
| 588 | 
  | 
  if( headAtomType != NULL ) delete headAtomType; | 
| 589 | 
  | 
  if( headBondType != NULL ) delete headBondType; | 
| 601 | 
  | 
#endif // is_mpi | 
| 602 | 
  | 
} | 
| 603 | 
  | 
 | 
| 604 | 
< | 
void TraPPE_ExFF::cleanMe( void ){ | 
| 604 | 
> | 
void DUFF::cleanMe( void ){ | 
| 605 | 
  | 
 | 
| 606 | 
  | 
#ifdef IS_MPI | 
| 607 | 
  | 
   | 
| 620 | 
  | 
} | 
| 621 | 
  | 
 | 
| 622 | 
  | 
 | 
| 623 | 
< | 
void TraPPE_ExFF::initForceField( int ljMixRule ){ | 
| 623 | 
> | 
void DUFF::initForceField( int ljMixRule ){ | 
| 624 | 
  | 
   | 
| 625 | 
  | 
  initFortran( ljMixRule, entry_plug->useReactionField ); | 
| 626 | 
  | 
} | 
| 627 | 
  | 
 | 
| 628 | 
  | 
 | 
| 629 | 
< | 
void TraPPE_ExFF::readParams( void ){ | 
| 629 | 
> | 
void DUFF::readParams( void ){ | 
| 630 | 
  | 
 | 
| 631 | 
  | 
  int i, a, b, c, d; | 
| 632 | 
  | 
  int identNum; | 
| 698 | 
  | 
    // send out the linked list to all the other processes | 
| 699 | 
  | 
 | 
| 700 | 
  | 
    sprintf( checkPointMsg, | 
| 701 | 
< | 
             "TraPPE_ExFF atom structures read successfully." ); | 
| 701 | 
> | 
             "DUFF atom structures read successfully." ); | 
| 702 | 
  | 
    MPIcheckPoint(); | 
| 703 | 
  | 
 | 
| 704 | 
  | 
    currentAtomType = headAtomType->next; //skip the first element who is a place holder. | 
| 710 | 
  | 
      sendFrcStruct( &atomInfo, mpiAtomStructType ); | 
| 711 | 
  | 
 | 
| 712 | 
  | 
      sprintf( checkPointMsg, | 
| 713 | 
< | 
               "successfully sent TraPPE_Ex force type: \"%s\"\n", | 
| 713 | 
> | 
               "successfully sent DUFF force type: \"%s\"\n", | 
| 714 | 
  | 
               atomInfo.name ); | 
| 715 | 
  | 
      MPIcheckPoint(); | 
| 716 | 
  | 
 | 
| 790 | 
  | 
       | 
| 791 | 
  | 
#ifdef IS_MPI | 
| 792 | 
  | 
  sprintf( checkPointMsg, | 
| 793 | 
< | 
           "TraPPE_ExFF atom structures successfully sent to fortran\n" ); | 
| 793 | 
> | 
           "DUFF atom structures successfully sent to fortran\n" ); | 
| 794 | 
  | 
  MPIcheckPoint(); | 
| 795 | 
  | 
#endif // is_mpi | 
| 796 | 
  | 
 | 
| 844 | 
  | 
    // send out the linked list to all the other processes | 
| 845 | 
  | 
     | 
| 846 | 
  | 
    sprintf( checkPointMsg, | 
| 847 | 
< | 
             "TraPPE_Ex bond structures read successfully." ); | 
| 847 | 
> | 
             "DUFF bond structures read successfully." ); | 
| 848 | 
  | 
    MPIcheckPoint(); | 
| 849 | 
  | 
     | 
| 850 | 
  | 
    currentBondType = headBondType->next; | 
| 874 | 
  | 
  } | 
| 875 | 
  | 
 | 
| 876 | 
  | 
  sprintf( checkPointMsg, | 
| 877 | 
< | 
           "TraPPE_ExFF bond structures broadcast successfully." ); | 
| 877 | 
> | 
           "DUFF bond structures broadcast successfully." ); | 
| 878 | 
  | 
  MPIcheckPoint(); | 
| 879 | 
  | 
 | 
| 880 | 
  | 
#endif // is_mpi | 
| 927 | 
  | 
    // send out the linked list to all the other processes | 
| 928 | 
  | 
 | 
| 929 | 
  | 
    sprintf( checkPointMsg, | 
| 930 | 
< | 
             "TraPPE_Ex bend structures read successfully." ); | 
| 930 | 
> | 
             "DUFF bend structures read successfully." ); | 
| 931 | 
  | 
    MPIcheckPoint(); | 
| 932 | 
  | 
 | 
| 933 | 
  | 
    currentBendType = headBendType->next; | 
| 957 | 
  | 
  } | 
| 958 | 
  | 
 | 
| 959 | 
  | 
  sprintf( checkPointMsg, | 
| 960 | 
< | 
           "TraPPE_ExFF bend structures broadcast successfully." ); | 
| 960 | 
> | 
           "DUFF bend structures broadcast successfully." ); | 
| 961 | 
  | 
  MPIcheckPoint(); | 
| 962 | 
  | 
 | 
| 963 | 
  | 
#endif // is_mpi | 
| 1012 | 
  | 
    // send out the linked list to all the other processes | 
| 1013 | 
  | 
     | 
| 1014 | 
  | 
    sprintf( checkPointMsg, | 
| 1015 | 
< | 
             "TraPPE_Ex torsion structures read successfully." ); | 
| 1015 | 
> | 
             "DUFF torsion structures read successfully." ); | 
| 1016 | 
  | 
    MPIcheckPoint(); | 
| 1017 | 
  | 
     | 
| 1018 | 
  | 
    currentTorsionType = headTorsionType->next; | 
| 1042 | 
  | 
  } | 
| 1043 | 
  | 
 | 
| 1044 | 
  | 
  sprintf( checkPointMsg, | 
| 1045 | 
< | 
           "TraPPE_ExFF torsion structures broadcast successfully." ); | 
| 1045 | 
> | 
           "DUFF torsion structures broadcast successfully." ); | 
| 1046 | 
  | 
  MPIcheckPoint(); | 
| 1047 | 
  | 
 | 
| 1048 | 
  | 
#endif // is_mpi | 
| 1052 | 
  | 
 | 
| 1053 | 
  | 
 | 
| 1054 | 
  | 
 | 
| 1055 | 
< | 
void TraPPE_ExFF::initializeAtoms( int nAtoms, Atom** the_atoms ){ | 
| 1055 | 
> | 
void DUFF::initializeAtoms( int nAtoms, Atom** the_atoms ){ | 
| 1056 | 
  | 
   | 
| 1057 | 
  | 
   | 
| 1058 | 
  | 
  ////////////////////////////////////////////////// | 
| 1141 | 
  | 
      else{ | 
| 1142 | 
  | 
         | 
| 1143 | 
  | 
        sprintf( painCave.errMsg, | 
| 1144 | 
< | 
                "TraPPE_ExFF error: Atom \"%s\" is a dipole, yet no standard" | 
| 1144 | 
> | 
                "DUFF error: Atom \"%s\" is a dipole, yet no standard" | 
| 1145 | 
  | 
                 " orientation was specifed in the BASS file.\n", | 
| 1146 | 
  | 
                 currentAtomType->name ); | 
| 1147 | 
  | 
        painCave.isFatal = 1; | 
| 1151 | 
  | 
    else{ | 
| 1152 | 
  | 
      if( the_atoms[i]->isDirectional() ){ | 
| 1153 | 
  | 
        sprintf( painCave.errMsg, | 
| 1154 | 
< | 
                 "TraPPE_ExFF error: Atom \"%s\" was given a standard" | 
| 1154 | 
> | 
                 "DUFF error: Atom \"%s\" was given a standard" | 
| 1155 | 
  | 
                 "orientation in the BASS file, yet it is not a dipole.\n", | 
| 1156 | 
  | 
                 currentAtomType->name); | 
| 1157 | 
  | 
        painCave.isFatal = 1; | 
| 1161 | 
  | 
  } | 
| 1162 | 
  | 
} | 
| 1163 | 
  | 
 | 
| 1164 | 
< | 
void TraPPE_ExFF::initializeBonds( int nBonds, Bond** bondArray, | 
| 1164 | 
> | 
void DUFF::initializeBonds( int nBonds, Bond** bondArray, | 
| 1165 | 
  | 
                                   bond_pair* the_bonds ){ | 
| 1166 | 
  | 
  int i,a,b; | 
| 1167 | 
  | 
  char* atomA; | 
| 1199 | 
  | 
  } | 
| 1200 | 
  | 
} | 
| 1201 | 
  | 
 | 
| 1202 | 
< | 
void TraPPE_ExFF::initializeBends( int nBends, Bend** bendArray, | 
| 1202 | 
> | 
void DUFF::initializeBends( int nBends, Bend** bendArray, | 
| 1203 | 
  | 
                                   bend_set* the_bends ){ | 
| 1204 | 
  | 
   | 
| 1205 | 
  | 
  QuadraticBend* qBend; | 
| 1280 | 
  | 
  } | 
| 1281 | 
  | 
} | 
| 1282 | 
  | 
 | 
| 1283 | 
< | 
void TraPPE_ExFF::initializeTorsions( int nTorsions, Torsion** torsionArray, | 
| 1283 | 
> | 
void DUFF::initializeTorsions( int nTorsions, Torsion** torsionArray, | 
| 1284 | 
  | 
                                      torsion_set* the_torsions ){ | 
| 1285 | 
  | 
 | 
| 1286 | 
  | 
  int i, a, b, c, d; | 
| 1327 | 
  | 
  } | 
| 1328 | 
  | 
} | 
| 1329 | 
  | 
 | 
| 1330 | 
< | 
void TraPPE_ExFF::fastForward( char* stopText, char* searchOwner ){ | 
| 1330 | 
> | 
void DUFF::fastForward( char* stopText, char* searchOwner ){ | 
| 1331 | 
  | 
 | 
| 1332 | 
  | 
  int foundText = 0; | 
| 1333 | 
  | 
  char* the_token; | 
| 1515 | 
  | 
    } | 
| 1516 | 
  | 
    else{ | 
| 1517 | 
  | 
      sprintf( painCave.errMsg,  | 
| 1518 | 
< | 
               "Unknown TraPPE_Ex bond type \"%s\" at line %d\n", | 
| 1518 | 
> | 
               "Unknown DUFF bond type \"%s\" at line %d\n", | 
| 1519 | 
  | 
               info.type, | 
| 1520 | 
  | 
               lineNum ); | 
| 1521 | 
  | 
      painCave.isFatal = 1; | 
| 1604 | 
  | 
     | 
| 1605 | 
  | 
    else{ | 
| 1606 | 
  | 
      sprintf( painCave.errMsg,  | 
| 1607 | 
< | 
               "Unknown TraPPE_Ex bend type \"%s\" at line %d\n", | 
| 1607 | 
> | 
               "Unknown DUFF bend type \"%s\" at line %d\n", | 
| 1608 | 
  | 
               info.type, | 
| 1609 | 
  | 
               lineNum ); | 
| 1610 | 
  | 
      painCave.isFatal = 1; | 
| 1702 | 
  | 
     | 
| 1703 | 
  | 
    else{ | 
| 1704 | 
  | 
      sprintf( painCave.errMsg,  | 
| 1705 | 
< | 
               "Unknown TraPPE_Ex torsion type \"%s\" at line %d\n", | 
| 1705 | 
> | 
               "Unknown DUFF torsion type \"%s\" at line %d\n", | 
| 1706 | 
  | 
               info.type, | 
| 1707 | 
  | 
               lineNum ); | 
| 1708 | 
  | 
      painCave.isFatal = 1; |