From 2da64ab8e3f31e2f107d95ae4a99748bf827d3df Mon Sep 17 00:00:00 2001
From: Liu Yuting L16101157 <101034703+L16101157@users.noreply.github.com>
Date: Sun, 6 Mar 2022 14:13:09 +0800
Subject: [PATCH 1/2] Update 5_practice.c

---
 5_practice.c | 3 +++
 1 file changed, 3 insertions(+)

diff --git a/5_practice.c b/5_practice.c
index 70d6651..092a658 100644
--- a/5_practice.c
+++ b/5_practice.c
@@ -2,6 +2,9 @@
 #include <stdlib.h>	// for memory allocation
 #include <time.h>	// for time calculation
 #include <math.h>	// for sine and cosine functions
+
+//For pull request
+
 int main() {
 	// Declare all the variables
 	int k, n, N;

From 317b8aaf6d5b4fbcaa1626ff19edd331aae6b3eb Mon Sep 17 00:00:00 2001
From: Liu Yuting L16101157 <101034703+L16101157@users.noreply.github.com>
Date: Wed, 9 Mar 2022 21:09:02 +0800
Subject: [PATCH 2/2] Raw DFT

---
 5_practice.c | 120 ++++++++++++++++++++++++++++++++++-----------------
 1 file changed, 81 insertions(+), 39 deletions(-)

diff --git a/5_practice.c b/5_practice.c
index 092a658..8d5262c 100644
--- a/5_practice.c
+++ b/5_practice.c
@@ -1,43 +1,85 @@
-#include <stdio.h>	// for printf function
-#include <stdlib.h>	// for memory allocation
-#include <time.h>	// for time calculation
-#include <math.h>	// for sine and cosine functions
+//
+//  DFT
+//
+//  Created by Liu Yuting on 2022/3/9.
+//
 
-//For pull request
+#include <stdio.h>
+#include <stdlib.h>
+#include <math.h>
+#include <time.h>
 
-int main() {
-	// Declare all the variables
-	int k, n, N;
-	double *x, *yr, *yi;
-	time_t t;
+int main(int argc, const char * argv[]) {
+    // Declare variables.
+    int N=7;
+    int n, k;
+    int index;
+   
+    double *x, *yr, *yi, *cos_array, *sin_array;
+    x=(double *) malloc(N*sizeof(double));
+    yr=(double *) malloc(N*sizeof(double));
+    yi=(double *) malloc(N*sizeof(double));
+    cos_array=(double *) malloc(N*sizeof(double));
+    sin_array=(double *) malloc(N*sizeof(double));
+    for(n=0; n<N; n++) {
+        //Initialize x.
+        x[n]=n;
+    }
+    time_t t;
+   
+    t=clock();
+   
+    cos_array[0]=1.0;
+    sin_array[0]=0.0;
+    cos_array[1]=cos(2*M_PI/N);
+    sin_array[1]=sin(2*M_PI/N);
+   
+    for(n=2; n<N; n++) {
+        //Compute table of cos and sin.
+        cos_array[n]=cos_array[1]*cos_array[n-1]-sin_array[1]*sin_array[n-1];
+        sin_array[n]=cos_array[1]*sin_array[n-1]+sin_array[1]*cos_array[n-1];
+    }
+   
+    for(k=0; k<N; k++) {
+        //Compute DFT.
+        yr[k]=0.0;
+        yi[k]=0.0;
+        for(n=0; n<N; n++) {
+            //Compute yr[k], yi[k].
+            index=(n*k)%N;
+            yr[k]+=x[n]*cos_array[index];
+            yi[k]+=x[n]*sin_array[index];
+        }
+    }
+   
+    t=clock()-t;
+   
+    printf("x:");
+    for(n=0; n<N; n++) {
+        //Print values of x.
+        printf(" %f", x[n]);
+    }
+    printf("\n");
+   
+    printf("yr:");
+    for(k=0; k<N; k++) {
+        //Print values of yr.
+        printf(" %f", yr[k]);
+    }
+    printf("\n");
+   
+    printf("yi:");
+    for(k=0; k<N; k++) {
+        //Print values of yi.
+        printf(" %f", yi[k]);
+    }
+    printf("\n");
+   
+    printf("Time: %f s\n", 1.0*t/CLOCKS_PER_SEC);
+   
+    free(yi);
+    free(yr);
+    free(x);
 
-	// Input the number N
-	
-	
-	// Locate the memory for x, yr, yi;
-
-	// Initial setting for x, for example, x[k] = k
-	
-	
-	
-	
-	t = clock();	
-	// yr[n]+i*yi[n] = sum(exp(-i 2 Pi k n / N)*x[k], k=0..N-1), n=0..N-1 
-
-
-
-
-
-
-
-	// output the results
-	t = clock() - t;
-	printf("%d ms for discrete Fourier Transform of %d elements\n", t, N);
-
-	// free the memory located by x, yr, yi
-	
-	
-	
-	return 100;
+    return 0;
 }
-