From ef385531fefa48663bfcf4bab62016e4f97b9736 Mon Sep 17 00:00:00 2001 From: Axel Kohlmeyer Date: Mon, 7 Nov 2011 20:53:37 -0500 Subject: [PATCH] fix bug with coefficient parsing for sdk angles --- src/USER-MISC/angle_sdk.cpp | 5 +++-- 1 file changed, 3 insertions(+), 2 deletions(-) diff --git a/src/USER-MISC/angle_sdk.cpp b/src/USER-MISC/angle_sdk.cpp index ba5c86fe92..93b8c19a26 100644 --- a/src/USER-MISC/angle_sdk.cpp +++ b/src/USER-MISC/angle_sdk.cpp @@ -244,7 +244,7 @@ void AngleSDK::coeff(int narg, char **arg) double k_one = force->numeric(arg[1]); double theta0_one = force->numeric(arg[2]); - double repscale_one = 1.0; + double repscale_one; // backward compatibility with old cg/cmm style input: // this had @@ -254,7 +254,8 @@ void AngleSDK::coeff(int narg, char **arg) if (narg == 6) { repscale_one = force->numeric(arg[4]); if (repscale_one > 0.0) repscale_one = 1.0; - } else if (narg == 3) repscale_one = force->numeric(arg[3]); + } else if (narg == 4) repscale_one = force->numeric(arg[3]); + else if (narg == 3) repscale_one = 1.0; else error->all(FLERR,"Incorrect args for angle coefficients"); // convert theta0 from degrees to radians and store coefficients