sine.c


1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52

/* 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 = FALSE;

double sin(double input, int precision){
  double ttl = input;
  int denom = 1;
  double num;
  for (int i=3;i<(2*precision)+2;i+=2) {
    num = pow(input,i);
    denom *= i*(i-1);
    ttl += ((i-1)/2 %2 == 1) ? -num/denom : num/denom;
  }
  return ttl;
}

int main(){
  const int precision = 20;
  const double threshold = 1e-9;

  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 {
      if (DEBUG){
        printf("%f IS ALL GOOD!\n",piValues[i][0]);
      }
    }
  }
  printf("All checks have passed!\n");
  return 0;
}