u

src/yapsut/periodic_stats.py

@@ -24,6 +24,17 @@ def periodic_stats(alpha,deg=False) :
    return np.arctan2(rr[1],rr[0]), np.nanvar(np.arctan2(vp,up))**0.5
 
 if __name__=='__main__' :
-   alpha_deg=np.arange(350,371)
+   print("Test1")
+   alpha=np.arange(350,371)
+   print(alpha)
    print(np.deg2rad(alpha-360).mean(),np.deg2rad(alpha-360).std())
    print(periodic_stats(np.deg2rad(alpha)))
+   print()
+
+   print("Test2")
+   alpha=np.array([-170,-171,-172,-173,-174,-175,-176,-177,-178,-179,-180,180,179,178,177,176,175,174,173,172,171,170])
+   print(alpha)
+   print(np.mod(alpha,360))
+   print(np.deg2rad(np.mod(alpha,360)).mean(),np.deg2rad(np.mod(alpha,360)).std())
+   print(periodic_stats(np.deg2rad(alpha)))
+   print()