src/gromacs/topology/ifunc.cpp

   1 /*
   2  * This file is part of the GROMACS molecular simulation package.
   3  *
   4  * Copyright (c) 1991-2000, University of Groningen, The Netherlands.
   5  * Copyright (c) 2001-2004, The GROMACS development team.
   6  * Copyright (c) 2013,2014,2015,2017,2018,2019, by the GROMACS development team, led by
   7  * Mark Abraham, David van der Spoel, Berk Hess, and Erik Lindahl,
   8  * and including many others, as listed in the AUTHORS file in the
   9  * top-level source directory and at http://www.gromacs.org.
  10  *
  11  * GROMACS is free software; you can redistribute it and/or
  12  * modify it under the terms of the GNU Lesser General Public License
  13  * as published by the Free Software Foundation; either version 2.1
  14  * of the License, or (at your option) any later version.
  15  *
  16  * GROMACS is distributed in the hope that it will be useful,
  17  * but WITHOUT ANY WARRANTY; without even the implied warranty of
  18  * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
  19  * Lesser General Public License for more details.
  20  *
  21  * You should have received a copy of the GNU Lesser General Public
  22  * License along with GROMACS; if not, see
  23  * http://www.gnu.org/licenses, or write to the Free Software Foundation,
  24  * Inc., 51 Franklin Street, Fifth Floor, Boston, MA  02110-1301  USA.
  25  *
  26  * If you want to redistribute modifications to GROMACS, please
  27  * consider that scientific software is very special. Version
  28  * control is crucial - bugs must be traceable. We will be happy to
  29  * consider code for inclusion in the official distribution, but
  30  * derived work must not be called official GROMACS. Details are found
  31  * in the README & COPYING files - if they are missing, get the
  32  * official version at http://www.gromacs.org.
  33  *
  34  * To help us fund GROMACS development, we humbly ask that you cite
  35  * the research papers on the package. Check out http://www.gromacs.org.
  36  */
  37 /* This file is completely threadsafe - keep it that way! */
  38 #include "gmxpre.h"
  39
  40 #include "ifunc.h"
  41
  42 #define  def_bonded(str, lstr, nra, nrpa, nrpb) \
  43     {str, lstr, (nra), (nrpa), (nrpb), IF_BOND}
  44
  45 #define  def_bondedz(str, lstr, nra, nrpa, nrpb) \
  46     {str, lstr, (nra), (nrpa), (nrpb), IF_BOND | IF_LIMZERO}
  47
  48 #define  def_bondedt(str, lstr, nra, nrpa, nrpb) \
  49     {str, lstr, (nra), (nrpa), (nrpb), IF_BOND | IF_TABULATED}
  50
  51 #define  def_bondedtz(str, lstr, nra, nrpa, nrpb) \
  52     {str, lstr, (nra), (nrpa), (nrpb), IF_BOND | IF_TABULATED | IF_LIMZERO}
  53
  54 #define   def_angle(str, lstr, nra, nrpa, nrpb) \
  55     {str, lstr, (nra), (nrpa), (nrpb), IF_BOND | IF_ATYPE}
  56
  57 #define    def_bond(str, lstr, nra, nrpa, nrpb) \
  58     {str, lstr, (nra), (nrpa), (nrpb), IF_BOND | IF_CHEMBOND | IF_BTYPE}
  59
  60 #define    def_bondt(str, lstr, nra, nrpa, nrpb) \
  61     {str, lstr, (nra), (nrpa), (nrpb), IF_BOND | IF_CHEMBOND | IF_TABULATED}
  62
  63 #define  def_bondnb(str, lstr, nra, nrpa, nrpb) \
  64     {str, lstr, (nra), (nrpa), (nrpb), IF_BOND | IF_CHEMBOND}
  65
  66 #define   def_vsite(str, lstr, nra, nrpa) \
  67     {str, lstr, (nra), (nrpa),     0, IF_VSITE}
  68
  69 #define     def_shk(str, lstr, nra, nrpa, nrpb) \
  70     {str, lstr, (nra), (nrpa), (nrpb), IF_CONSTRAINT}
  71
  72 #define   def_shkcb(str, lstr, nra, nrpa, nrpb) \
  73     {str, lstr, (nra), (nrpa), (nrpb), IF_CONSTRAINT | IF_CHEMBOND}
  74
  75 #define      def_nb(str, lstr, nra, nrp) \
  76     {str, lstr, (nra), (nrp),     0, IF_NULL}
  77
  78 #define    def_nofc(str, lstr) \
  79     {str, lstr,    0,     0,     0, IF_NULL}
  80
  81 /* this MUST correspond to the enum in src/gromacs/topology/ifunc.h */
  82 const t_interaction_function interaction_function[F_NRE] =
  83 {
  84     def_bond    ("BONDS",    "Bond",            2, 2, 2),
  85     def_bond    ("G96BONDS", "G96Bond",         2, 2, 2),
  86     def_bond    ("MORSE",    "Morse",           2, 3, 3),
  87     def_bond    ("CUBICBONDS", "Cubic Bonds",    2, 3, 0),
  88     def_bondnb  ("CONNBONDS", "Connect Bonds",   2, 0, 0),
  89     def_bonded  ("HARMONIC", "Harmonic Pot.",   2, 2, 2),
  90     def_bondnb  ("FENEBONDS", "FENE Bonds",     2, 2, 0),
  91     def_bondt   ("TABBONDS", "Tab. Bonds",      2, 2, 2),
  92     def_bondedtz("TABBONDSNC", "Tab. Bonds NC", 2, 2, 2),
  93     def_bonded  ("RESTRAINTPOT", "Restraint Pot.", 2, 4, 4),
  94     def_angle   ("ANGLES",   "Angle",           3, 2, 2),
  95     def_angle   ("G96ANGLES", "G96Angle",        3, 2, 2),
  96     def_angle   ("RESTRANGLES", "Restricted Angles", 3, 2, 2),
  97     def_angle   ("LINEAR_ANGLES", "Lin. Angle", 3, 2, 2),
  98     def_bonded  ("CROSS_BOND_BOND", "Bond-Cross", 3, 3, 0),
  99     def_bonded  ("CROSS_BOND_ANGLE", "BA-Cross",   3, 4, 0),
 100     def_angle   ("UREY_BRADLEY", "U-B",          3, 4, 4),
 101     def_angle   ("QANGLES", "Quartic Angles",    3, 6, 0),
 102     def_bondedt ("TABANGLES", "Tab. Angles",    3, 2, 2),
 103     def_bonded  ("PDIHS",    "Proper Dih.",     4, 3, 3),
 104     def_bonded  ("RBDIHS",   "Ryckaert-Bell.",  4, 6, 6),
 105     def_bonded  ("RESTRDIHS",  "Restricted Dih.",     4, 2, 2),
 106     def_bonded  ("CBTDIHS",   "CBT Dih.",  4, 6, 6),
 107     def_bonded  ("FOURDIHS", "Fourier Dih.",    4, 4, 4),
 108     def_bonded  ("IDIHS",    "Improper Dih.",   4, 2, 2),
 109     def_bonded  ("PIDIHS",   "Improper Dih.",   4, 3, 3),
 110     def_bondedt ("TABDIHS", "Tab. Dih.",        4, 2, 2),
 111     def_bonded  ("CMAP",  "CMAP Dih.",          5, -1, -1),
 112     def_nofc    ("GB12",     "GB 1-2 Pol. (unused)" ),
 113     def_nofc    ("GB13",     "GB 1-3 Pol. (unused)" ),
 114     def_nofc    ("GB14",     "GB 1-4 Pol. (unused)" ),
 115     def_nofc    ("GBPOL",    "GB Polarization (unused)" ),
 116     def_nofc    ("NPSOLVATION", "Nonpolar Sol. (unused)" ),
 117     def_bondedz ("LJ14",     "LJ-14",           2, 2, 2),
 118     def_nofc    ("COUL14",   "Coulomb-14"                                           ),
 119     def_bondedz ("LJC14_Q",  "LJC-14 q",        2, 5, 0),
 120     def_bondedz ("LJC_NB",   "LJC Pairs NB",    2, 4, 0),
 121     def_nb      ("LJ_SR",    "LJ (SR)",         2, 2                                ),
 122     def_nb      ("BHAM",     "Buck.ham (SR)",   2, 3                                ),
 123     def_nofc    ("LJ_LR",    "LJ (unused)"                                          ),
 124     def_nofc    ("BHAM_LR",  "B.ham (unused)"                                       ),
 125     def_nofc    ("DISPCORR", "Disper. corr."                                        ),
 126     def_nofc    ("COUL_SR",  "Coulomb (SR)"                                         ),
 127     def_nofc    ("COUL_LR",  "Coul (unused)"                                        ),
 128     def_nofc    ("RF_EXCL",  "RF excl."                                             ),
 129     def_nofc    ("COUL_RECIP", "Coul. recip."                                       ),
 130     def_nofc    ("LJ_RECIP", "LJ recip."                                            ),
 131     def_nofc    ("DPD",      "DPD"                                                  ),
 132     def_bondnb  ("POLARIZATION", "Polarization", 2, 1, 0),
 133     def_bonded  ("WATERPOL", "Water Pol.",      5, 6, 0),
 134     def_bonded  ("THOLE",    "Thole Pol.",      4, 3, 0),
 135     def_bondnb  ("ANHARM_POL", "Anharm. Pol.", 2, 3, 0),
 136     def_bonded  ("POSRES",   "Position Rest.",  1, 3, 3),
 137     def_bonded  ("FBPOSRES", "Flat-bottom posres", 1, 3, 0),
 138     def_bonded  ("DISRES",   "Dis. Rest.",      2, 6, 0),
 139     def_nofc    ("DISRESVIOL",   "D.R.Viol. (nm)"                                       ),
 140     def_bonded  ("ORIRES",   "Orient. Rest.",   2, 6, 0),
 141     def_nofc    ("ORDEV",    "Ori. R. RMSD"                                         ),
 142     def_bonded  ("ANGRES",   "Angle Rest.",     4, 3, 3),
 143     def_bonded  ("ANGRESZ",  "Angle Rest. Z",   2, 3, 3),
 144     def_bonded  ("DIHRES",   "Dih. Rest.",      4, 3, 3),
 145     def_nofc    ("DIHRESVIOL",  "Dih. Rest. Viol."                                     ), /* obsolete */
 146     def_shkcb   ("CONSTR",   "Constraint",      2, 1, 1                             ),
 147     def_shk     ("CONSTRNC", "Constr. No Conn.", 2, 1, 1                             ),
 148     def_shkcb   ("SETTLE",   "Settle",          3, 2, 0                             ),
 149     def_vsite   ("VSITE2",   "Virtual site 2",  3, 1                                ),
 150     def_vsite   ("VSITE3",   "Virtual site 3",  4, 2                                ),
 151     def_vsite   ("VSITE3FD", "Virtual site 3fd", 4, 2                                ),
 152     def_vsite   ("VSITE3FAD", "Virtual site 3fad", 4, 2                               ),
 153     def_vsite   ("VSITE3OUT", "Virtual site 3out", 4, 3                               ),
 154     def_vsite   ("VSITE4FD", "Virtual site 4fd", 5, 3                               ),
 155     def_vsite   ("VSITE4FDN", "Virtual site 4fdn", 5, 3                               ),
 156     def_vsite   ("VSITEN",   "Virtual site N",   2, 2                               ),
 157     def_nofc    ("COM_PULL", "COM Pull En."     ),
 158     def_nofc    ("EQM",      "Quantum En."      ),
 159     def_nofc    ("EPOT",     "Potential"        ),
 160     def_nofc    ("EKIN",     "Kinetic En."      ),
 161     def_nofc    ("ETOT",     "Total Energy"     ),
 162     def_nofc    ("ECONS",    "Conserved En."    ),
 163     def_nofc    ("TEMP",     "Temperature"      ),
 164     def_nofc    ("VTEMP",    "Vir. Temp. (not used)"      ),
 165     /* Note that pressure names can not be more than 8 char's,
 166      * because " (bar)" is appended to them.
 167      */
 168     def_nofc    ("PDISPCORR", "Pres. DC"         ),
 169     def_nofc    ("PRES",     "Pressure"         ),
 170     def_nofc    ("DH/DL_CON", "dH/dl constr."    ), /* obsolete */
 171     def_nofc    ("DV/DL",    "dVremain/dl"      ),
 172     def_nofc    ("DK/DL",    "dEkin/dl"         ),
 173     def_nofc    ("DVC/DL",   "dVcoul/dl"        ),
 174     def_nofc    ("DVV/DL",   "dVvdw/dl"         ),
 175     def_nofc    ("DVB/DL",   "dVbonded/dl"      ),
 176     def_nofc    ("DVR/DL",   "dVrestraint/dl"   ),
 177     def_nofc    ("DVT/DL",   "dVtemperature/dl" )
 178 };