Move complex numbers to their own header

2 files changed, 68 insertions, 61 deletions

diff --git a/complex.h b/complex.h
new file mode 100644
index 0000000..0dcde9c
--- /dev/null
+++ b/complex.h
@@ -0,0 +1,67 @@
+#ifndef COMPLEX_H
+
+#include <math.h>
+
+typedef struct
+{
+    double re;
+    double im;
+} complex_t;
+
+static inline complex_t
+complex_mul (complex_t a, complex_t b)
+{
+    complex_t r;
+
+    r.re = a.re * b.re - a.im * b.im;
+    r.im = a.re * b.im + a.im * b.re;
+
+    return r;
+}
+
+static inline complex_t
+complex_add (complex_t a, complex_t b)
+{
+    complex_t r;
+
+    r.re = a.re + b.re;
+    r.im = a.im + b.im;
+
+    return r;
+}
+
+static inline complex_t
+complex_sub (complex_t a, complex_t b)
+{
+    complex_t r;
+
+    r.re = a.re - b.re;
+    r.im = a.im - b.im;
+
+    return r;
+}
+
+static inline double
+complex_mag (complex_t a)
+{
+    return sqrt (a.re * a.re + a.im * a.im);
+}
+
+static inline double
+complex_arg (complex_t a)
+{
+    return fmod (atan2 (a.im, a.re) + 2 * M_PI, 2 * M_PI);
+}
+
+static inline complex_t
+complex_from_mag_arg (double mag, double arg)
+{
+    complex_t r;
+
+    r.re = mag * cos (arg);
+    r.im = mag * sin (arg);
+
+    return r;
+}
+
+#endif
diff --git a/fft.h b/fft.h
index 61d4ba9..2900542 100644
--- a/fft.h
+++ b/fft.h
@@ -1,4 +1,5 @@
 #include <math.h>
+#include "complex.h"
 
 #ifndef FALSE
 #define FALSE (0)
@@ -8,67 +9,6 @@
 #define TRUE (1)
 #endif
 
-typedef struct
-{
-    double re;
-    double im;
-} complex_t;
-
-static inline complex_t
-complex_mul (complex_t a, complex_t b)
-{
-    complex_t r;
-
-    r.re = a.re * b.re - a.im * b.im;
-    r.im = a.re * b.im + a.im * b.re;
-
-    return r;
-}
-
-static inline complex_t
-complex_add (complex_t a, complex_t b)
-{
-    complex_t r;
-
-    r.re = a.re + b.re;
-    r.im = a.im + b.im;
-
-    return r;
-}
-
-static inline complex_t
-complex_sub (complex_t a, complex_t b)
-{
-    complex_t r;
-
-    r.re = a.re - b.re;
-    r.im = a.im - b.im;
-
-    return r;
-}
-
-static inline double
-complex_mag (complex_t a)
-{
-    return sqrt (a.re * a.re + a.im * a.im);
-}
-static inline double
-complex_arg (complex_t a)
-{
-    return fmod (atan2 (a.im, a.re) + 2 * M_PI, 2 * M_PI);
-}
-
-static inline complex_t
-complex_from_mag_arg (double mag, double arg)
-{
-    complex_t r;
-
-    r.re = mag * cos (arg);
-    r.im = mag * sin (arg);
-
-    return r;
-}
-
 static inline double
 rad_to_degree (double theta)
 {