Special cases working

softwareQinc · Sep 5, 2023 · 0a466aa · 0a466aa
1 parent cb0cd02
commit 0a466aa
Showing 1 changed file with 8 additions and 2 deletions.
diff --git a/include/grid_synth/exact_synthesis.hpp b/include/grid_synth/exact_synthesis.hpp
@@ -9,8 +9,14 @@ namespace grid_synth {
 
 // Returns known common cases for multiples of theta = pi/4. w =
 // 1/(sqrt(omega)).
-inline str_t check_common_cases(real_t theta, real_t eps) {
-    if(abs(theta-real_t("0.25")) < eps) {
+inline str_t check_common_cases(real_t theta, const real_t& eps) {
+
+  while(theta > real_t("2"))
+      theta = theta - real_t("2");
+  while(theta < 0)
+      theta = theta + real_t("2");
+
+  if(abs(theta-real_t("0.25")) < eps) {
         return "T w";
     }
     else if(abs(theta-real_t("0.5") < eps)) {