Use constexpr for physical constants and move them into gmx namespace

[alexxy/gromacs.git] / src / gromacs / listed_forces / restcbt.cpp

diff --git a/src/gromacs/listed_forces/restcbt.cpp b/src/gromacs/listed_forces/restcbt.cpp
index 2da9dcf8fec015aec0d135b40ab8c274beb77d79..68326ebca3929eff65dd0a1d0766fc87e38c0719 100644 (file)
--- a/src/gromacs/listed_forces/restcbt.cpp
+++ b/src/gromacs/listed_forces/restcbt.cpp
@@ -1,7 +1,7 @@
 /*
  * This file is part of the GROMACS molecular simulation package.
  *
- * Copyright (c) 2014,2015,2017,2019, by the GROMACS development team, led by
+ * Copyright (c) 2014,2015,2017,2019,2021, by the GROMACS development team, led by
  * Mark Abraham, David van der Spoel, Berk Hess, and Erik Lindahl,
  * and including many others, as listed in the AUTHORS file in the
  * top-level source directory and at http://www.gromacs.org.
@@ -78,7 +78,7 @@ void compute_factors_restangles(int             type,
     double term_theta_theta_equil;
 
     k_bending          = forceparams[type].harmonic.krA;
-    theta_equil        = forceparams[type].harmonic.rA * DEG2RAD;
+    theta_equil        = forceparams[type].harmonic.rA * gmx::c_deg2Rad;
     theta_equil        = M_PI - theta_equil;
     cosine_theta_equil = cos(theta_equil);
 
@@ -135,7 +135,7 @@ void compute_factors_restrdihs(int             type,
     real norm_phi;
 
     /* Read parameters phi0 and k_torsion */
-    phi0        = forceparams[type].pdihs.phiA * DEG2RAD;
+    phi0        = forceparams[type].pdihs.phiA * gmx::c_deg2Rad;
     cosine_phi0 = cos(phi0);
     k_torsion   = forceparams[type].pdihs.cpA;