singular_values.plt prints traces instead of singular values

Makefile

@@ -1,13 +1,13 @@
 HEADERS=triangle.h linalg.h queue.h
 
-#SPECIAL_OPTIONS=-O0 -g -D_DEBUG
+SPECIAL_OPTIONS=-O0 -g -D_DEBUG
 #SPECIAL_OPTIONS=-O3 -pg -funroll-loops -fno-inline
-SPECIAL_OPTIONS=-O3 -flto -funroll-loops -Winline
+#SPECIAL_OPTIONS=-O3 -flto -funroll-loops -Winline
 #SPECIAL_OPTIONS=
 
 OPTIONS=-m64 -march=native -mtune=native -std=gnu99 -D_GNU_SOURCE $(SPECIAL_OPTIONS)
 
-all: discreteness singular_values element_path limit_set hyperbolic
+all: discreteness singular_values element_path limit_set hyperbolic ellipticity
 
 singular_values: singular_values.o triangle.o linalg.o
 	gcc $(OPTIONS) -o singular_values triangle.o linalg.o singular_values.o -lm -lgsl -lcblas
@@ -24,6 +24,9 @@ limit_set: limit_set.o linalg.o triangle.o
 hyperbolic: hyperbolic.o triangle.o linalg.o
 	gcc $(OPTIONS) -o hyperbolic hyperbolic.o triangle.o linalg.o -lm -lgsl -lcblas
 
+ellipticity: ellipticity.o triangle.o linalg.o
+	gcc $(OPTIONS) -o ellipticity ellipticity.o triangle.o linalg.o -lm -lgsl -lcblas
+
 singular_values.o: singular_values.c $(HEADERS)
 	gcc $(OPTIONS) -c singular_values.c
 
@@ -45,5 +48,8 @@ limit_set.o: limit_set.c $(HEADERS)
 hyperbolic.o: hyperbolic.c $(HEADERS)
 	gcc $(OPTIONS) -c hyperbolic.c
 
+ellipticity.o: ellipticity.c $(HEADERS)
+	gcc $(OPTIONS) -c ellipticity.c
+
 clean:
-	rm -f singular_values discreteness element_path limit_set hyperbolic triangle.o linalg.o singular_values.o discreteness.o element_path.o limit_set.o hyperbolic.o
+	rm -f singular_values discreteness element_path limit_set hyperbolic ellipticity triangle.o linalg.o singular_values.o discreteness.o element_path.o limit_set.o hyperbolic.o ellipticity.o