Continued with question 3, not finished yet.

simmod/ass1/.gitignore

+speed.*.*
+speed
+fp
+float_double
+report.pdf

simmod/ass1/Makefile

@@ -2,7 +2,7 @@ FLAGS=-Wall -Wextra -std=c99 -pedantic -O0 -lm
 TYPES=float double LD
 OPS=ADD DIV MULT SQRT
 
-all: fp speed report.pdf
+all: fp speed report.pdf float_double
 
 %.pdf: %.tex
 	pdflatex $^
@@ -21,6 +21,9 @@ speed: speed.c
 fp: floating_point.o
 	gcc $(FLAGS) -o $@ $^
 
+float_double: float_vs_double.o
+	gcc $(FLAGS) -o $@ $^
+
 %.o: %.c
 	gcc $(FLAGS) -o $@ -c $^
 

simmod/ass1/float_vs_double.c

+#include <stdlib.h>
+#include <stdio.h>
+
+int main(void) {
+    puts("Using type {TYPE}.");
+
+    puts("Forward summation:");
+    printf("N = 1e8: %e\n", sum_forward(1e8));
+    printf("N = 2e8: %e\n", sum_forward(2e8));
+
+    puts("Backward summation:");
+    printf("N = 1e8: %e\n", sum_backward(1e8));
+    printf("N = 2e8: %e\n", sum_backward(2e8));
+
+    return 0;
+}
+
+{TYPE} sum_forward(int N) {
+    {TYPE} sum = 0;
+
+    for( int i = 1; i <= N; i++ )
+        sum += 1.0 / i;
+
+    return sum;
+}
+
+{TYPE} sum_backward(int N) {
+    {TYPE} sum = 0;
+
+    for( int i = N; i; i-- )
+        sum += 1.0 / i;
+
+    return sum;
+}