Initial commit

diff --git a/Makefile b/Makefile
new file mode 100644 (file)
index 0000000..f66c367
--- /dev/null
+++ b/Makefile
@@ -0,0 +1,11 @@
+CC=gcc
+
+all: verify.c
+       $(CC) verify.c -o verify
+
+debug: verify.c
+       $(CC) verify.c -o verify -DDEBUG
+
+clean:
+       rm -f verify
+

diff --git a/README.md b/README.md
new file mode 100644 (file)
index 0000000..ac9fdb0
--- /dev/null
+++ b/README.md
@@ -0,0 +1,4 @@
+# Pearson verify
+This program is meant to be used to verify the ouput of the pearson section of https://github.com/sm13294/dv1674_a2.\
+\
+This program compares 64 bit floating point numbers stored in a separate line each from \<file1\> and \<file2\>. If all the pairs of numbers are the same up to the 15th decimal place, 0 will be returned. If not, this program will check if all numbers are the same up to the 11th decimal place, and if so 1 will be returned. If [stop] is 1, this program will stop as soon as this happens. If any numbers are not the same up to the 11th decimal place, 2 will be returned. If [stop] is 2, this program will stop as soon as this happens. If there is an error while executing this program, -1 will be returned.

diff --git a/verify.c b/verify.c
new file mode 100644 (file)
index 0000000..b2bdd91
--- /dev/null
+++ b/verify.c
@@ -0,0 +1,91 @@
+#include <stdio.h>
+#include <stdlib.h>
+#include <math.h>
+
+#define ERROR_15 (0.000000000000001)
+#define ERROR_11 (0.00000000001)
+
+
+int main(int argc, char *argv[]){
+
+       int ret = 0;
+       int stop = 0;
+       FILE *f1 = NULL;
+       FILE *f2 = NULL;
+       long unsigned int line = 1;
+
+       if(argc < 3){
+               fprintf(stderr, "ERROR: Wrong usage.\n");
+               fprintf(stderr, "Usage: %s <file1> <file2> [stop]\n", argv[0]);
+               fprintf(stderr,
+                       "\tThis program compares 64 bit floating point numbers stored in a separate\n"
+                       "\tline each from <file1> and <file2>. If all the pairs of numbers are the\n"
+                       "\tsame up to the 15th decimal place, 0 will be returned. If not, this program\n"
+                       "\twill check if all numbers are the same up to the 11th decimal place, and\n"
+                       "\tif so 1 will be returned. If [stop] is 1, this program will stop as soon\n"
+                       "\tas this happens. If any numbers are not the same up to the 11th decimal\n"
+                       "\tplace, 2 will be returned. If [stop] is 2, this program will stop as soon\n"
+                       "\tas this happens. If there is an error while executing this program, -1 will\n"
+                       "\tbe returned. Compile with -DDEBUG (make debug) or define DEBUG in the program\n"
+                       "\tto print debug information.\n"
+               );
+               return -1;
+       }
+
+       if(argc > 3){
+               stop = atoi(argv[3]);
+       }
+
+       f1 = fopen(argv[1], "rb");
+       if(!f1){
+               fprintf(stderr, "ERROR: Cannot open file '%s'.\n", argv[1]);
+               return -1;
+       }
+       f2 = fopen(argv[2], "rb");
+       if(!f1){
+               fprintf(stderr, "ERROR: Cannot open file '%s'.\n", argv[2]);
+               return -1;
+       }
+
+       while(!feof(f1) && !feof(f2)){
+               double d1, d2;
+       
+               if(fscanf(f1, "%lg\n", &d1) != 1){
+                       fprintf(stderr, "ERROR: Cannot read number from file '%s' at line %lu.\n", argv[1], line);
+                       return -1;
+               }
+               if(fscanf(f2, "%lg\n", &d2) != 1){
+                       fprintf(stderr, "ERROR: Cannot read number from file '%s' at line %lu.\n", argv[1], line);
+                       return -1;
+               }
+               
+               double error = fabs(d1-d2);
+               if(error >= ERROR_15){
+                       if(error < ERROR_11){
+                               ret = 1;
+                       }else{
+                               ret = 2;
+                       }
+                       if(stop > 0 && stop == ret){
+                               goto end;
+                       }
+               }
+               line++;
+       }
+
+       if((feof(f1) && !feof(f2)) || (!feof(f1) && feof(f2))){
+               fprintf(stderr, "ERROR: Different number of lines in files '%s' and '%s'.\n", argv[1], argv[2]);
+               return -1;
+       }
+
+end:
+       fclose(f1);
+       fclose(f2);
+
+#ifdef DEBUG
+       printf("%lu lines read\n", line);
+       printf("Return value: %d\n", ret);
+#endif
+
+       return ret;
+}