 -            return T(table.factorial[n]);
 +            return T(const_tables::factorial[n]);
      }
- 
+
      template <>
 @@ -163,7 +703,7 @@ namespace tables {
          } else if (n > MAXDOUBLEFACT)

-            return T(table.invsqrt_int[n]);
+            return T(const_tables::invsqrt_int[n]);
     }
 
     template <>
-    Multi invsqrt_int(int n) {
+    inline Multi invsqrt_int(int n) {

+                return T(const_tables::half_factorial_pos[n]);
         }
     }

diff --git a/starry/tables.py b/starry/tables.py
new file mode 100644
index 0000000..0ad0dac
