From ba0ddea5e1ff2d27fd77e3aa7263466705c64c2e Mon Sep 17 00:00:00 2001 From: Anders Hafreager Date: Sun, 28 May 2017 15:44:12 -0700 Subject: [PATCH] Using correct ndegree instead of nnn --- src/compute_hexorder_atom.cpp | 8 ++++---- 1 file changed, 4 insertions(+), 4 deletions(-) diff --git a/src/compute_hexorder_atom.cpp b/src/compute_hexorder_atom.cpp index 93b84080bc..013036f364 100644 --- a/src/compute_hexorder_atom.cpp +++ b/src/compute_hexorder_atom.cpp @@ -248,7 +248,7 @@ inline void ComputeHexOrderAtom::calc_qn_complex(double delx, double dely, doubl double x = delx*rinv; double y = dely*rinv; std::complex z(x, y); - std::complex zn = pow(z, nnn); + std::complex zn = pow(z, ndegree); u = real(zn); v = imag(zn); } @@ -259,9 +259,9 @@ inline void ComputeHexOrderAtom::calc_qn_complex(double delx, double dely, doubl inline void ComputeHexOrderAtom::calc_qn_trig(double delx, double dely, double &u, double &v) { double ntheta; if(fabs(delx) <= MY_EPSILON) { - if(dely > 0.0) ntheta = nnn * MY_PI / 2.0; - else ntheta = nnn * 3.0 * MY_PI / 2.0; - } else ntheta = nnn * atan(dely / delx); + if(dely > 0.0) ntheta = ndegree * MY_PI / 2.0; + else ntheta = ndegree * 3.0 * MY_PI / 2.0; + } else ntheta = ndegree * atan(dely / delx); u = cos(ntheta); v = sin(ntheta); }