Merged simmod/modsim directories.

.gitignore

 *.log
 *.aux
+*.o
+*.pdf
 *.swp
 *.swo
 *.toc
-*.o
-ass*.tar.gz
+*.gz
 compilerbouw/
 robotica/

simmod/ass1/src/Makefile → modsim/ass1/Makefile

@@ -4,31 +4,29 @@ SPEED_TYPES=float double LD
 SUM_TYPES=float double
 OPS=ADD DIV MULT SQRT
 
-all: fp speed sum kahan pr
-
-#all: fp speed highlight report.pdf sum kahan pr
-#
-#highlight: floating_point.tex extra_precision.tex \
-#	speed.tex sum.tex kahan_sum.tex Makefile.tex
-#
-#s%.tex: s%.c
-#	pygmentize -O style=colorful -o $@ $^
-#
-#kahan_sum.tex: kahan_sum.c
-#	pygmentize -O style=colorful -o $@ $^
-#
-#Makefile.tex: Makefile
-#	pygmentize -O style=colorful -o $@ $^
-#
-#extra_precision.tex: extra_precision.c
-#	pygmentize -O style=colorful -o $@ $^
-#
-#floating_point.tex: floating_point.c
-#	pygmentize -O style=colorful -o $@ $^
-#
-#%.pdf: %.tex
-#	pdflatex $^
-#	pdflatex $^
+all: fp speed highlight report.pdf sum kahan pr
+
+highlight: floating_point.tex extra_precision.tex \
+	speed.tex sum.tex kahan_sum.tex Makefile.tex
+
+s%.tex: s%.c
+	pygmentize -O style=colorful -o $@ $^
+
+kahan_sum.tex: kahan_sum.c
+	pygmentize -O style=colorful -o $@ $^
+
+Makefile.tex: Makefile
+	pygmentize -O style=colorful -o $@ $^
+
+extra_precision.tex: extra_precision.c
+	pygmentize -O style=colorful -o $@ $^
+
+floating_point.tex: floating_point.c
+	pygmentize -O style=colorful -o $@ $^
+
+%.pdf: %.tex
+	pdflatex $^
+	pdflatex $^
 
 speed: speed.c
 	for t in $(SPEED_TYPES); do \

modsim/ass2/.gitignore

+q1
+q2
+q3
+q4
+q5
+*.o

modsim/ass2/Makefile

+CC=clang
+CFLAGS=-Wall -Wextra -pedantic -std=c99 -D_GNU_SOURCE
+LFLAGS=-lm
+
+all: q1 q2 q3
+
+q1: q1.o
+	$(CC) $(CFLAGS) $(LFLAGS) -o $@ $^
+
+q2: q2.o
+	$(CC) $(CFLAGS) $(LFLAGS) -o $@ $^
+
+%.o: %.c
+	$(CC) $(CFLAGS) $(LFLAGS) -o $@ -c $^
+
+clean:
+	for i in `seq 5`; do \
+		rm -vf q$$i; \
+	done;
+	rm -vf *.o

modsim/ass2/q1.c

+#include <stdlib.h>
+#include <stdio.h>
+#include <math.h>
+
+#define H 1e-3
+#define TABLE_LINE(func, x) (printf("%-24s%.12f\t%.12f\n", #x, \
+            slope_right(func, x, H), slope_central(func, x, H)))
+
+typedef double (*func_ptr)(double x);
+
+double slope_right(func_ptr func, double x, double h) {
+    return (func(x + h) - func(x)) / h;
+}
+
+double slope_central(func_ptr func, double x, double h) {
+    return (func(x + h) - func(x - h)) / (2 * h);
+}
+
+int main(void) {
+    puts("x\t\t\tright\t\tcentral");
+    TABLE_LINE(&sin, M_PI / 3);
+    TABLE_LINE(&sin, 100 * M_PI + M_PI / 3);
+    TABLE_LINE(&sin, 1e12 * M_PI + M_PI / 3);
+
+    return 0;
+}

modsim/ass2/q2.c

+#include <stdlib.h>
+#include <stdio.h>
+#include <math.h>
+
+#define EPSILON 1e-11
+
+typedef double (*func_ptr)(double x);
+
+double bisec(func_ptr f, double left, double right, int *steps) {
+    int i;
+    double mid;
+
+    for( i = 1; fabs(right - left) > 2 * EPSILON; i++ ) {
+        mid = (right + left) / 2;
+
+        if( f(left) * f(mid) < 0 )
+            right = mid;
+        else if( f(right) * f(mid) < 0 )
+            left = mid;
+        else
+            break;
+    }
+
+    *steps = i;
+
+    return mid;
+}
+
+double func(double x) {
+    return x * sin(x) - 1;
+}
+
+int main(void) {
+    int steps;
+
+    printf("zero point: %.20f\n", bisec(&func, 0, 2, &steps));
+    printf("Steps: %d\n", steps);
+
+    return 0;
+}

modsim/ass2/q3.c

+#include <stdlib.h>
+#include <stdio.h>
+#include <math.h>
+
+#define EPSILON 1e-11
+
+typedef double (*func_ptr)(double x);
+
+double bisec(func_ptr f, double left, double right, int *steps) {
+    int i;
+    double mid;
+
+    for( i = 1; fabs(right - left) > 2 * EPSILON; i++ ) {
+        mid = (right + left) / 2;
+
+        if( f(left) * f(mid) < 0 )
+            right = mid;
+        else if( f(right) * f(mid) < 0 )
+            left = mid;
+        else
+            break;
+    }
+
+    *steps = i;
+
+    return mid;
+}
+
+double func(double x) {
+    return x * sin(x) - 1;
+}
+
+int main(void) {
+    int steps;
+
+    printf("zero point: %.20f\n", bisec(&func, 0, 2, &steps));
+    printf("Steps: %d\n", steps);
+
+    return 0;
+}

os/extra/src/report/.gitignore

+*.log
+*.toc
+*.aux
+report.pdf
+report-dot2tex*
+dot2tex.cache

os/extra/src/report/Makefile

+all:
+	pdflatex report.tex