1 files changed, 50 insertions, 0 deletions
diff --git a/src/z_prime.c b/src/z_prime.c
new file mode 100644
index 0000000..2af18b2
--- /dev/null
+++ b/src/z_prime.c
@@ -0,0 +1,50 @@
+#include "zexy.h"
+#include <math.h>
+
+
+static t_class *prime_class;
+
+typedef struct _prime {
+  t_object  x_obj;
+} t_prime;
+
+
+void prime_float(t_prime *x, t_float f)
+{
+
+  unsigned int i=f;
+  unsigned int max_divisor;
+  int divisor=1;
+
+  if (f<2)return;
+
+  if (!(i%2)){
+    if (i==2)outlet_bang(x->x_obj.ob_outlet);
+    return;
+  }
+
+  max_divisor = sqrt(f)+1;
+
+  while ((divisor+=2)<max_divisor)
+    if (!(i%divisor)) return;
+
+  outlet_bang(x->x_obj.ob_outlet);
+}
+
+void *prime_new(void)
+{
+  t_prime *x = (t_prime *)pd_new(prime_class);
+
+  outlet_new(&x->x_obj, &s_float);
+
+  return (void *)x;
+}
+
+void z_prime_setup(void) {
+  prime_class = class_new(gensym("prime"),
+			  (t_newmethod)prime_new,
+			  0, sizeof(t_prime),
+			  CLASS_DEFAULT, 0);
+
+  class_addfloat(prime_class, prime_float);
+}