Semi functional Sine function

1 files changed, 50 insertions, 0 deletions

diff --git a/sine.c b/sine.c
new file mode 100644
index 0000000..871acc0
--- /dev/null
+++ b/sine.c
@@ -0,0 +1,50 @@
+/* sine.h
+ * ----------------------------------------------------------------------------
+ * "THE BEER-WARE LICENSE" (Revision 42):
+ * <benjamin@chausse.xyz> wrote this file. As long as you retain this notice
+ * you can do whatever you want with this stuff. If we meet some day, and you
+ * think this stuff is worth it, you can buy me a beer in return.
+ * Benjamin Chausse
+ * ----------------------------------------------------------------------------
+ */
+
+#include <stdio.h>
+#include "format.h"
+#include "testData.h"
+
+const int DEBUG = TRUE;
+
+float sin(float input, int precision){
+  float ttl = input;
+  int denom = 1;
+  float num;
+  for (int i=3;i<(2*precision)+2;i+=2) {
+    num = mpow(input,i);
+    denom *= i*(i-1);
+    ttl += ((i-1)/2 %2 == 1) ? -1*num/denom : num/denom;
+  }
+  return ttl;
+}
+
+int main(){
+  const int precision = 16;
+  const float threshold = 1e-4;
+
+  for (int i=0; i<COUNT_OF(piValues);i++){
+    if ( abs(piValues[i][1]-sin(piValues[i][0],precision))>threshold ){
+      printf("ERROR: \n"
+             "Deviation between calculated value and expected value\n"
+             "surpasses threshold of %f for number %f\n"
+             "Expected: %f but got: %f\n",
+             threshold,
+             piValues[i][0],
+             piValues[i][1],
+             sin(piValues[i][0],precision));
+      return 1;
+    } else {
+      printf("%f IS ALL GOOD!\n",piValues[i][0]);
+    }
+  }
+
+  return 0;
+}