CORDIC.rar_cordic_cordic C_cordic sin cos

/* A testbench for sine and cosine using CORDIC. It generates angles between 1 and NUM_DEGREE and provides those to the cordic sin and cos function. Then it compares the result to corresponding sin/cos result computed from the same functions in math.h. It computes the error of the cordic sin/cos result vs. the math.h results. The values and errors are logged in the file out.dat for debugging purposes. The cumulative errors are printed. VARIABLES: double theta: Input angle double &s: Reference to the sine part double &c: Reference to the cos part error_sin= [abs(s-zs)/zs]*100; error_cos= [abs(c-zc)/zc]*100; Total_Error_Sin = sum(error_sin) Total_error_Cos = sum(error_cos) */ #include <math.h> #include "cordic.h" #include <stdio.h> #include <stdlib.h> double abs_double(double var){ if ( var < 0) var = -var; return var; } int main(int argc, char **argv) { FILE *fp; FILE *fs; FILE *fc; cos_sin_TYPE s; //sin output cos_sin_TYPE c; //cos output theta_TYPE radian; //radian input double zs, zc; // sin and cos values calculated from math.h. //Error checking double Total_Error_Sin=0.0; double Total_error_Cos=0.0; double error_sin=0.0, error_cos=0.0; fp=fopen("out.dat","w"); fs=fopen("sin.dat","w"); fc=fopen("cos.dat","w"); for(int i=1;i<NUM_DEGREE;i++) { radian = i*M_PI/180; cordic(radian, s, c); zs = sin((double)radian); zc = cos((double)radian); error_sin=(abs_double((double)s-zs)/zs)*100.0; error_cos=(abs_double((double)c-zc)/zc)*100.0; Total_Error_Sin=Total_Error_Sin+error_sin; Total_error_Cos=Total_error_Cos+error_cos; fprintf(fp, "degree=%d, radian=%f, cos=%f, sin=%f\n", i, (double)radian, (double)c, (double)s); fprintf(fs, "%f\n",(double)error_sin); fprintf(fc, "%f\n",(double)error_cos); } fclose(fp); fclose(fs); fclose(fc); printf ("Total_Error_Sin=%f, Total_error_Cos=%f, \n", Total_Error_Sin, Total_error_Cos); return 0; }