close issue #5, fix CPP algorithm

madwyn · Nov 30, 2015 · ed286cf · ed286cf
1 parent 505b3c1
commit ed286cf
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/src/Cube2Cyl.hpp b/src/Cube2Cyl.hpp
@@ -367,7 +367,7 @@ inline void Cube2Cyl::calCubeFace(const double& theta, const double& phi) {
     }
 
     // find out which segment the line strikes to
-    phiThreshold = atan2(S_RADIUS / cos(normTheta), S_RADIUS);
+    phiThreshold = atan2(S_RADIUS, S_RADIUS / cos(normTheta));
 
     // in the top segment
     if (phi > phiThreshold) {