From 1b52f921966d463c816b13bccdd1adc6991b1d30 Mon Sep 17 00:00:00 2001
From: Frederik Beimgraben <frederik@beimgraben.net>
Date: Sat, 11 Jan 2025 18:15:53 +0100
Subject: [PATCH] Add a few tests

---
 .gitignore          |   2 +-
 MMS24-25.c          |  54 +++++++++++++++++++++++--------
 Makefile            |  12 +++++--
 testInterpolation   | Bin 16040 -> 0 bytes
 testInterpolation.c |  31 ------------------
 tests               | Bin 0 -> 16328 bytes
 tests.c             |  77 ++++++++++++++++++++++++++++++++++++++++++++
 7 files changed, 128 insertions(+), 48 deletions(-)
 delete mode 100755 testInterpolation
 delete mode 100644 testInterpolation.c
 create mode 100755 tests
 create mode 100644 tests.c

diff --git a/.gitignore b/.gitignore
index 254bdcf..cb50d2d 100644
--- a/.gitignore
+++ b/.gitignore
@@ -54,4 +54,4 @@ dkms.conf
 
 main
 
-result.txt
\ No newline at end of file
+*result.txt
\ No newline at end of file
diff --git a/MMS24-25.c b/MMS24-25.c
index e210532..1f08442 100644
--- a/MMS24-25.c
+++ b/MMS24-25.c
@@ -17,12 +17,39 @@
 #include <stdlib.h>
 #include "MMS24-25.h"
 #include <stdio.h>
+#include <stdarg.h>
 
-// NULL Pointer
+// Define NULL if not already provided
 #ifndef NULL
 #define NULL 0
 #endif
 
+void _error(char* fmt,...) {
+    printf("<\x1B[31;1m ERROR \x1B[0m: \x1B[34m'");
+
+    va_list args;
+    va_start(args, fmt);
+    printf(fmt, args);
+    va_end(args);
+
+    printf("'\x1B[0m >\n");
+
+    exit(-1);
+}
+
+void _debug(char* fmt,...) {
+    printf("[\x1B[31;1m DEBUG \x1B[0m] \x1B[34m'");
+
+    va_list args;
+    va_start(args, fmt);
+    printf(fmt, args);
+    va_end(args);
+
+    printf("'\x1B[0m\n");
+
+    exit(-1);
+}
+
 double interpolateLine(
         double dX1,
         double dY1,
@@ -30,13 +57,18 @@ double interpolateLine(
         double dY2,
         double dXq
     ) {
-    double y = (
-        (dY2 - dY1) / (dX2 - dX1)
-    ) * (dXq + dX1);
+    if (dX2 == dX1)
+        _error("interpolateLine: dX1 and dX2 must differ!");
+
+    double y = dY1 * (
+        (dX2 - dXq) / (dX2 - dX1)
+    ) + dY2 * (
+        (dXq - dX1) / (dX2 - dX1)
+    );
 
     #ifdef _DEBUG_
         // DEBUG: print the parameters and result
-        printf("interpolateLine(%lf, %lf, %lf, %lf, %lf) -> %lf\n",
+        _debug("interpolateLine(%lf, %lf, %lf, %lf, %lf) -> %lf\n",
            dX1,dY1,dX2,dY2,dXq,y);
     #endif
 
@@ -77,10 +109,8 @@ double *generateSineValues(
     double* sineArray = NULL;
     sineArray = (double*) malloc(sizeof(double) * iNumValues);
 
-    if (sineArray == NULL) {
-        printf("Failed to allocate memory for SineArray\n");
-        exit(-1);
-    }
+    if (sineArray == NULL) 
+        _error("*generateSineValues: Failed to allocate memory for Sine Array\n");
 
     for (int i=0; i<iNumValues;i++){
         sineArray[i] = dAmp*sin(((2*M_PI)/iNumSamplesPerPeriod) * i);
@@ -94,10 +124,8 @@ void writeDoubleArray(double *pdOut, int iLen, char *pcOutName) {
      */
 
     FILE* datei = fopen(pcOutName, "w");
-    if (datei == NULL) {
-        printf("Error while opening file '%s'.\n", pcOutName);
-        exit(-1);
-    }
+    if (datei == NULL)
+        _error("writeDoubleArray: Error while opening file '%s'.\n", pcOutName);
 
     for (int i = 0; i < iLen; i++) {
         fprintf(datei, "%f\n", pdOut[i]);
diff --git a/Makefile b/Makefile
index 5da918e..b36eebb 100644
--- a/Makefile
+++ b/Makefile
@@ -2,7 +2,13 @@ CC=gcc
 SOURCE_FILE=MMS24-25.c
 OPTS=-lm
 
-all: testInterpolation
+all: clean tests run
 
-testInterpolation:
-	$(CC) testInterpolation.c $(SOURCE_FILE) -o testInterpolation $(OPTS)
\ No newline at end of file
+clean:
+	rm ./tests || true
+
+tests:
+	$(CC) tests.c $(SOURCE_FILE) -o tests $(OPTS)
+
+run:
+	./tests
\ No newline at end of file
diff --git a/testInterpolation b/testInterpolation
deleted file mode 100755
index af0685cd68ddc296e39227d44bb43d10245ded64..0000000000000000000000000000000000000000
GIT binary patch
literal 0
HcmV?d00001

literal 16040
zcmeHOeQX@n5r22ihxxG2PLtyBl?^4}p!i}NCk+jD&WYDqgPpjEO${PkF1|hcs`H(D
zw?}N0q><cGu9nqMRe)OlFqKMMP(@S_1)8>JOtC|0<$zj66bXxl>Z}9^!L1C0TxZ_S
zocA_ogs9Yiy3g7-Z+<iLX5Q@E^V{9G-&R|j-7c5lQY0D#a_KS?@sXhNtJDDE6PrW{
z)Ypn@#1h~aNz9b{OafA~oOI{YV#14nqFn_x3ec-fSTN-Xk)qu~rHmXf6*LaCTZ~4;
z+v#IcC<IfM=lm&X3`LQ3jyBQwm|E>HH|Z)<KV6?sA0fL@qHKpO+D(z&6xlI7ivdAl
z%K5~c;IWDFX{SDrVe!dY2)o??+1Y6m=`%$iSa!c-1iy&CcCr(z&3^2Z{V-*JUk5wP
z%jYe8ns<@Eh0cSI^3RmgE7nKC`!?3Cj|8hC;b?N8YGBjGs*QEtM9f<&L%_JI;lOFL
zeMhHo9~NoE*xn^vlD`BnEN<dij=C5l@@S6*F+8v!T@DAH7L=FVxbpG4kN#rp_U7AP
zXqfnA)r;R-Dvi*d4`_bw5(3R-1+=j*r2u{#V1@MWS^)pS0{Gwp_`m}AHGnH{*oz;4
zLh*DgfHyCI?_2<12e<-<y&MMMIy;EGzonziAJpS|PdH)d@s74_kyupk2<(gK!td|t
zi$(njBM>+Iep9ueVj>(B{YfJs^ntJ;La}~5Dgs6<B>Do8NUR$=f@#@J9rXo3JJcT!
zM~x6jx(5ROP&gWhgzwWus5=r%=mNC+aL#VTR=3&O-rBNto4?ju>)jx@eRCZS0p=*o
zH~ur*i-5-#%xCcs=-`}H;7DE>hN-<sJOOPiF0dX2%GaGKFOqV^o_X*gxP0B>S!iS7
z-$VW#I$Cbxlk|JX_+G->eWsjo=F?kD96E;>=F6l7r%PN4rz|-4g@Th7ykx!`F>S$1
zEja&>E(MH!I;FJOXDT7i*D@kyKCFd6oc`5FzQQMo#BmH<6rb5beLT2S`pgz${M&*^
zwGV4k_foSd>rFnBLH*0jrmTHFlR^A)v+2}M2Am8y8E`V-WWdS5|49Zu^j!O~Hgu*~
z8(uuNUI^{GX~UhJ(uQ6t9+w%&Zg>UY?3!1h_N?+j3;AA@%}i&Zta%Z6Tzbt+N&W@o
zamh84mi(`g$0gFtQOQ4nJT8r9MkW7K<Z($fb42n#Mjn?!GlwPrL*#LpGSe^l?<0>(
zx0$`bWBfJmwxmwnqov-~hEAPr?@(*fHOIA)`Zqu)yW;O&A!b8f&#LbL3j-Xf4?qXn
z$l|+@+%jjZfT8~$lyK;mWT!o=4x`>UQP8iE`WtRSvi5?OI;)Mnw?!M9E7DvqYbP!k
z%b|l2>Yz9~9rCPdnzzUJ4%fd6Vv)SLQyZ%PCN{OyJ4T5%T>lqfvMaiv&UC|IUS9kX
zG+lef=i?m@&9^cC*xyVK=ypOMaN^oi2ajo~<k1QBRm{%N!DCB|N*LX`QuH#BX}yrm
z0$WdPJ+Q;-s~~t91;gr5Z72;VU+Pu^6w*A0(ppNLB3T*9mRn>eNmfC!xliZqrb(ud
z>`jX-L$XSer7g0v6Y3=9{d73wc|U#5bCWt5@>HHY=jqyW_ME4zGXquIG*l{jp90xG
z;CKU$6L3sI*wJ)m19ZRxm;<#r3B&Lcu%w4<JcrXX69%_K>ZEJv;AGLE2O%5B$0?b<
zbDnY)Je8xT@=owxzGn(5eG(l*>uwc8Xgi5+f$Pkm?{l7eUO_oDy3qR0F)&f*0MfDA
zgT|d;40cedm?5+~1INW<fM8F+dK`nrYany%d^U?#a};FiU?!D3g$^*fwka*uIXR)8
z2UT&pVOTw{J)(^)|0A;KMsuNG0WqmlcHyKHcFBIkkb2$)3Y*na2VOXGWDh4qNgX_s
zO2U}RA?yit4*MNe&w!J$=^~t7^$;G+mN{+c;G9UVgy?TKqtCqlNjCcjoRXOq(3wzY
z(e`uz&hBe)cIEk8ydQKqgM8X&Ay_<}hv=Ms0{zHyI@~#{b-$*KT>e+=NgM8*(_F8}
z6O6;d6XqULXJHhxo7FQ19>&w3JF(T8nK#H5PM?$7$cjI~-H?|08xG<Da4|oKHBuq<
zGK}RZsWb5$I?UAO1~dRI9?)j+cTvu@VYqC}+4my2dLZY@oPcv5;{>EoPaXIj>S7#e
zz>gdsm)E%bP0sxaV*pNZGT>yu$$*mqCj(9foD4V_a5C`UX8_;ZVA+xE+o#8sSV(D%
z_ayuDsF6@s*DD*@LuVq;qbr7<Fj}JUdaFMc2^is6RJk?SRnq|VowDw#m2FvfHHbB#
z+^Yz9NkvwAYdES`u8D-!DHm63m8u4Wb9T*va6}I(Mof`!zygM@^yz)E_@EMs#g*OA
zRbxCJ7!>;{P&FP0e19)kDfs#<9PLp;$lSOlaijOrTu2pFO&y=CUg>jPRkRsy*X27F
z<^60H-$SN>?gje%*=%+c=#q2UY#Qj3K&OFT@nJTLFPq*1s?bX(d{>Lbb>A-G8YpvJ
zb@|fbBOpT@_b_xppEsaF(^C10?0RV8OSjSYOdq1Othuz}HqT}Iiw}z}S8l%P=IgIR
zb+o-3j&g{l4ShA1mfi2(R<Z=1K7b1;KLBvU2ife;fac`xQ4~x0qX3VCJd!8x$;poa
z?t!t%8*uL5J)V=l3G&~7{Fyxar=%Pqr#KmKGT>yu$$*mqCj(9foD4V_a5C_}$pG)0
z<9%{;J0kb$VFb*NQaJF&#VoW@kN3L$XoU%{Br)&LyNY;x;fm!-lL+3YcQx@W|JQ|V
z452h_isQX*8Z~&2+ce4fMsJMz;XP^<#8X(r&F{6#(7rtL?iKd$kp|y~;R|CdSC~X>
zB7pG$3WVcMlOFG>WB>4~M~(-?QWai<O1_fp;RToEIS#l@ll-M5|A=@RE_wOi6W%{=
z_q>tjK{L@ih=z#XOZ0xC4-tKw=(9w*U*4CwecQILDV4Bud96}Y2YXFxl-lary6R2U
zu<5k0=$4vl?xX<4dAvw`khP0&eiXt{uMn;P?h-46Js)e)K)`0x)v}SRZiSW$W)aUf
z7VLt?GG;6kPrhQ|x%sxlD}?4_#H7PO3ib2o0{E{0FB27_%Lffh!F%&U{Y)V4_6WYd
z+2p+i^#3gJQgH=cNRwnR4HW&_<9}D`lahttCnbu7=Vlq;<?ag6?qkHH*C39&hZU1w
zMB;UT7aG?Vz{|u+vDZh6ME(hMnJCSDzj6cNO{8znTf6MPR6I}njO8;Hu<0cHEy@Gt
z7fXonC%<I!@|g-W2MJ$7xO^4^cuL}?uC@Nj0{D{);Ku;RyxHTQ1YEIpOzpP;zpeoO
zDd5Z8@L#}o-HMAI?RmRJTifp1x~kd@(DXOmxudbIWgBcXhl|>;YX0V3jcux5YuY9J
z+go>TZEW@LY;NAIcKACQx3;RV%iX?L-DlZt{#m=_aRYo#7k8=S4){F4uLlE0fHt}3
zOJE0lSoG*oJr0iw@$ujtfk;wMKs0D#?wsdM^ZAyxJ^KD&BIfT6M1#10zhx&_2E$Q*
zGNA`)Gk(6c9u0zzeTf9Qkel=O$HRu+6idQ40OnK5xSr?^MDifpgHdzSe_r44Q9vG#
z-v;CXem{I|;8%Avedf0X`Cka+cK}H<gzzQ?`;5Rophn!Jz1%XR6y9je(7ip;q_;l~
zPgmo{AVG~=TdIsek40f$G7KNjgo8q=_XZNZ!W$fnf(w%xaZ_@y9#6o7<UGa?vbY`z
zpaHe|BZlzGFujI80Bw1qz44fw30}R|ABqS1bboIUY?(7H`~mrdSXyu!Iu7)OyTNM=
z`hv<U9hjZLSrp#xSRXzGFK`KcAqdw|aOu(0_KUxJ<8>5wj=I+J6<~1f!t0EMzVBii
zUe_>>5Av<{_V2tDRwT;5zYOy;a`yiSEohkH^K5Ie+kXRu|Aw^?+w=O6DUY9)O}X#W
zxX%DIaBax;ynbXVcOZfs$B4P(@$)(^3>b#Y_PictItDUpvFHC=AjGvE+w=O4>9u6c
z@v|P&Lm<bsBlEm|WZF*le2zGNg2f0lFsE$K>sO{}Qe^*Z&*OiX?6;C0uZNlP`WfwU
z&1|=S3^2?I_n-U!LH|eiJ^@r@v7dij&sgmhvST`Gn&oQ#_r&z)Hv2RQn7(dv$adTz
zuA8j>c|F9mgcXU}x)V>6J?CF2q{tLIw-nAFL%)Rvs^siflZ0p!D-yN)KMuInUL!lE
z=WX_O|F44__JGQ309E<^j=zhv>od>vEwJ%e?D_t0vI={fpD!ji&$In$5a8V5`ygI_
z^LMWNJuA+CENsW;aTYKfKil(nU+MLzEkvbBay8pC{RkvhdtT3H$bN6m(rmCDQ#Tk}
z?Rg!nRh#-8CwIekIUbA=<M??U*G2Yree`85CFGyfb0W6Sy{IJ#+goXF?&Yk7I6@ur
t{{#FT#&f{&@b3rS7jx_ub5UGpaz(bz*)Qwj^16`yWfb_AZ3Z@0{0H~DM1TMQ

diff --git a/testInterpolation.c b/testInterpolation.c
deleted file mode 100644
index e50bdb9..0000000
--- a/testInterpolation.c
+++ /dev/null
@@ -1,31 +0,0 @@
-//
-// Created by frederik on 1/11/25.
-//
-
-#include "MMS24-25.h"
-#include <stdio.h>
-#include <stdlib.h>
-
-
-int main (int iArgC, char** ppcArgV){
-    int iRes = -1;
-
-    if (iArgC != 6){
-        printf("Number of Arguments != 6\n");
-        printf("Usage testInterpolation <dX1> <dY1> <dX2> <dY2> <dX>\n");
-        exit(iRes);
-    };
-
-    double dX1 = atof(ppcArgV[1]);
-    double dY1 = atof(ppcArgV[2]);
-    double dX2 = atof(ppcArgV[3]);
-    double dY2 = atof(ppcArgV[4]);
-    double dX = atof(ppcArgV[5]);
-
-    for(int i=0; i<4; i++){
-        double dRes = interpolateLine(dX1, dY1, dX2, dY2, dX+(double)i);
-        printf ("%f\n",dRes);
-    }
-
-    return iRes;
-}
diff --git a/tests b/tests
new file mode 100755
index 0000000000000000000000000000000000000000..0b0af51e5982bc4400102625d7c7a9fc66513516
GIT binary patch
literal 16328
zcmeHOe{dAneSdcXTgJdi4m8*#_Cok22M!4Yk-^x<j}wbW0s<0&!T30Id+5f}oqRt)
zP$0B|Fc(!ba^2XT+S4JKNlMa5!z7N&I3YbuZDda}Bd4{8W|Br@<K{G$Z3WA<h$+|4
z_r3RBR*TMc>x|n!?%nKue}CWiefPb$Z+G8)HP}$^^>_rAGI5_EYM?@K0$O0!X#^l3
z>ckYB?-ZXGH-VlYc}gBo08&*hddI?Y%`XEayBUOx13#zOf+0spknD<;eq8}Wq1!=r
z6Uj)^Zg}=q<&h!F?fMKQ7cMwi2gA^5N4d!@Qg*pjW8o>Cm(v=u9g$?$ukHG^9m8P?
z2!$c%lfs~*PUq7N1CUV($g>b`yH0KAhU>IGL#l(y);H6U|4n_HwcWrh<;M-viY*v&
zeP4zh<>liRzRX*6eZ{UrK<A&KPOn(p7TdGDc5z!|QClpY>0H!Vw|vp^T7N3xUm`=G
zxT<kbpIX1MS$L0%9P!xRBRrD61w1NVO|zVIDMq5no(l2kph9pO4!SHPziD36iWe?C
zdgtFA9=#Nq>7Bjo>Xkbnf?pvjagZO%&ruxer1<a7YVi@PSUCd0ZQ}^*i}0x)V;ud>
zMf87AgikFVuii(C@Oz8!KM%fv!(H}+F`oY?i|`*Q!fz<Tj~3xC1mD2nE+@g5?Ph|0
z7{_g5s`&9|l<?nJC$yodF%+?q*4|huZ6%u;*S00%R#SLSn<YY_z3qv3D3uN;)1i=%
zs^cmq6Pb8Kq+)T=kx8du2a}eELoKZjhoa$Fo3J`#X<;RkiKK`oI;^;854W`?TA+yT
z-=ZsM55qp%k&MODQL#_MP^hId9E!%`;kMX8OGI1R5-Cf-s-4Eja>7~^w!UG*nzf-N
z{w4lpg7I1#TUtvE_M)GA(B~&$JkeR<6TnqCMiYL|LX*affg^KA3`46-ybPR*2Y3#U
z<W#1`VpC{zz6P9lN`HdXC;FL3AIEU<iUHtMcz))2v7ZJNN}ry8nSWICcL$`{o=2I^
z-LLp?j+m~OJ_lb<wNlvc;PYHTf{PBma{v!G_}o_MC?L<`W&oNKxReJJMtq)Ih*J@u
zvyl4ozNknsBcKq{r#5gg0tzEO9b6U!6h=JG5pk*mbXK~P6e`bk0Y#9$zU4|rz`>ua
zROH!o@Nq%yvf06(YI93SDFdYplrm7tKq&*I41Amnyyv^~4`$a;x!FDO++rcj<GHlA
z&~JAAQ~7ymURd@r*oC=2$JuvR02tA&BpVzkpv--qXj<wF_DlLXqG<^-n3MDmiKZpb
z;Au&Jk7!zo44#zqw}_@C$Y8Ie|C(r8dJG<w^j{K9%Z0%XNq?PaT51gLmh_(!O-rZ2
z&5}MyG%Yy>1E8sWHE(aozVeWneckN3G`zVfSd*)H$?RFxgh~q2f8rNnB<lCwbqpp*
zK+mc#!-3f|@c{z&jHag}HGcyoO3f{W0pDFmN$<Rd(4Qssm)$|2>YAAyHqX6tzj<!7
z%=ElyzH%)+4G!A1gYv>a)OXjq>-H4i(N({JSY(zqn_a6OAk@shnVw>HulgB?!t^?v
z^DSu1ixZy#>e+Svdc5bO*KxW&<Zo~<berMB{Ij6}>Fw+SaLnwFw@dI1N}QSPX$5G=
z9;dW6WP3URj^)yZuPS)dSJl|-t7<;!tJ-ziS7qgVRqg%0s?30|>hN&H?4B@pfpmM^
zgiTEj$vTU*nVG$C25ga^ajj6eaJ(8ax=a_2*8;01gPvf8nLS<y((E}-R|^zfB(vum
zq-}N|S8=`uS9veI!1EF1Ise#y77D%4peiHg7gY%9iTS6%f9d?E+uC8QQ0>;<jQ~D*
zYp1~rTHC8U=f8b*{MLRS{7<{JhO@QzB7jfc+IjGT*1nD>0J_xw@L{3wkGQUd{3;0O
zT5iRhv>h6KN<;OlD&*+0uVP|30=n$04E7OGrip$cu0fq77eSD?GI(-u4dx=7JT*AY
zNuD13cg$S2^vvKHCz%`kx+YIXeg5F7gB7HyaNn_qoU)ulmh&B3qh+VaPRkVTJ2p?r
zcIDu#`3z2tr^xMTvzuHmcmv0;aJ+(}FKYVYx%{J;J%Js`nc3g~ba#<KGhvP7<fBZt
zIZX3pK&Oen#q@aDw}adt(RIxEL!m(Ac~YXNCO&tgQsIhJihqb^jUK#}qI(bMaC#dI
z+R3#g<D;m?Y1;c}(h!atn}{p@9;R>_Q#8}(XW~k~<>`f+Abf47>sd)yMX*2WTO7>!
zPDC*~5RFi^n1WW~GrkjbT723T1JmNuz7vzQ7>U%@z_j?3@5HZFM+=^`#lW=qr0>Mv
zLyU+Q>WC4r$p0z2Fs7$nqk4394OL_hUp2FtVbyo%!$xAC_UQc(bFe*E9U1*2j9%5C
zrErBLoKf!cBOmY(J_y05J00KaY~oW0<P%R&s)?K73QfF%E=jXgrQO6QfXOCy4F#xY
zqUTpO4&c;G9d%bfPF5fFn2R+2^9POVcxD<r{y=$DwY;r68%!URQ(7ZoG*BNphrt48
zJLFg=*@IM7g+l?tAT+i=k1B(sR1J5)?BRSibBWj#Q)9oGZSK1eyh@$*t^2xzSIsBQ
zo@p--iHRp~df=r(NOx{wpA@QbBD#WCJy2K~ymaWf-rim6#@jVyWDgHzGYGCSIAC4~
zj*{Q*;1DV~H&8~&oAV#YWi~i!b{!rSG6}!J^2_d+dAh9xN6750aMbtiFE9J<Zl<<Q
zT#8)LRETuE3_%Wp;D~GT`rv&vdF8sBN6eOAnmv<WB~NB|^Qh^0nWib2?V!1{(%Tgr
zK?6os28RwkftGW+x`wJvl`<siA_N!Bp6UB!h;lV@3{_l@f2LH(zKCnN)7gw!j%KVx
zEjJ|D?5?0nDWvC7!DhPxy77qoJ0pdH9pI%yXGxc`_Z{$iU!oPcy<WFh=+x3uDP^FP
zfl>zk=Na(muLV12?p*q%l{M|ghK==`jF~&D+jrdpR&D!yF`p<w3QI^~YxCN*!L3`z
z>f8%y4M~GrwrtuWWh-#LZg1Qt(pD<HA&zfZ9f`JZI+lnps9GuN@fpgB7~xdPN)peA
zCXy?R1#{b?3k{sA7CLZa(R~hTCCPPz=j0SJC!nDEN^<m=fAz-+jG88cUl&}{yk0eI
zOz?1uk6$`)t05M*Rv3}(HAXld0bF9VXHsb+5{pKy<SZ3{xO=Y^x02zswH4ai!fhEV
zwZhPiO(zWbnGI5--D*!H4^XQO(lk~lli>rC1bx1gWq&f3w$>#w`0OV|D~uq0t~2(x
zB5wwLvx~*|8c`zV&rQwuPZD#Z^zviIqIFFl&bmF|nOU}yzNOQ*lfXNL!YT6lF774(
zR{;9~e>Yqx3<JLU`$C}t_n{wME)*65rf|>O4EPve2jD*3!~X(s2bMpl0BOzF2e=z6
zzC7R;@OV&8Z+=kd+cg!>!7akmS>c&E`R4Lo=n$XwbhILtd8DxJ<_ZFHaMa@HeY;S|
z0g9;=^;2g&;JbBy`B8EI9V_o%I%hVq$@W1Uvk?0`#8^GG;%nZuQ*QDOq6Q^D1opc{
z<TD_98pm0*S$=C}d+-b#v?lv1@{pr;uq@bBR`#sd_Wx^$X5x-)0NOe`rC5;vt2jvh
zhEoolWHXS1*T<244zgOvzKdKx<+SrYZr9rcIjy7G#>n^D@-*ZfkkdNBj?4QqTmDy&
ze+Jj=dt>C^R)G_*R7x2rWuTOSQU*#HC}p6Ofl>zk-)DgLsqwxvEb(Nay>e8T&PhrL
zd@WH$uDaE80WnU`RP5Uou=m4hw!ODYZaRey?~9{#D-|rgWqJQvAwftEA0DXi-m)h(
zf)`%0<h7ihRj9nBfW0S-*2+|-D`20snaJyXaQbmWl<S7_a|eA!pz>J-M4x6bzg#06
zcTVfkhh!>TKVIw0_JAc-@!>$yr*uQ`(Mr-B2OcXVeN4;oNF-^FTRQ(gIPv~*ch$>v
zKd9Gmn}$&h_i6Yw4IkI=I~qQ#;for+rXkm}JD{-jYuDnj9y^(<jG9{P1+6ibR4=Km
zuB*nT(D6n0)Ks&RaUfmp3F7w!w}`Iaczn_uk5B!@gI?y&5A7?a9n`!(U6MFaOdj+q
z?tjd$2XY$|4%kwJp8?-+=rKAO4KvXDwg;G_V6lDC--5UT`mHYg?}9&GJdeQtc>9fi
zO8Q=(Jr97C<$0XP^Ya_YpDHeBp?nqxkDhPG(|=Rz2lNDT9yFCv92C!Edf>|EZty0-
z54FpkhtC%A|3&Z(*tyFJ@GHGD1ouOhua|zNifwwJzo7*iNuPEM8;U_2RjBNe`U-cR
z>AeWWpVAWz``HhEg}7aG1f1GJ94g|cLi-ue<7b{Wqjx6c$9;XD*7{p@y@uBRrskLF
ziCaGFL-s?>{}(;M%XbLizo7Z|Xuf>+0RF3zubn!;?~3q0D8i@rL6m=Y{PY>#aQJ1^
z5B_W>9B>2pw|Q~*()nT0CgLOi|El+8b3c@C8&m!k(b%|kN$sK~%YcT~?bx`wal=|{
zZ^uLw3Yww%Evp-YA#>dp5nA7{Y0c_}(5CwOt-+>H)9N)1LF}w|?}HEU4tDHYmz(Qt
z25oQus6G3%0e`Hr6$z)qdguID3AWY8ZrH4EtJI2ArlrM7rNoWz7;k*Z3O`hj+RM+|
z_-$V@YS+IBLXlJ=)QWeT^iP9`wI{O|4ci0<kyt#GNm&v7FM+X6)V~Nq7XELb|06I~
z9@>*i=^*4E3fwO?Mbb*Ogxj3&I^8c_)qfPmxW!)^#?bVKhB1s#2>+Z21vjqy$p1-<
z{Wrr{2e=k#L=weq+@Ctoo(}H;OeYm=WlX&v7XEl5ZTa`cGyaYwUd|@d2Q+K-nhlH6
z;k^!pJ((E(T@;H5sook+wF-aaKpa&lm`*CmJ}a5Ri}5i;2(qNr7A6CYb+o00U&iWB
zTb;mVC-Nr~vX}X-)=)GVZnr|M5!f=NEJ9&=>8vamhvRU2tOd0u;0q_etUx)#r4;^_
zM0>jxPxC<hgpj5yOd<Mf1nwhc${|BuXL2>9>nvXYK}$MJ2eQ;@dxincXTa+|rk4Y8
z+IMPuh6XEY$oB<iV@ux{`zN(MLwdjMEN=U+K!{hj367uF&kWgKFDvc!Hv0?%!R9Bs
z{k*<rcv>fn+DIt3pVzrD@F--ipV#XQd3_EQS={-51YBARvpuhW8M@<VJ%;o@G+J9T
z&FgE1yJfA{bIkE;w&(#uIc0l$eXrN~Ij!i<AGiMrZNEnA@$&*he!d`ky5(~HZ2wL0
zC@1XSzEA1<6+cHoMHct<e;Qn;y`k+G=9HOz=G!2{zjfIUXaU1tOtMWS!gkF4zRRAU
z2N}}$BWH2P|04(~e$Kxzw6=x=By<$cAG3Y}#A(m#Zieo8#$EqQF8l4Xm1l;}{5SRg
zGt{Tz(8Jb#{^kFTy5nV<;csE<bJ+9qao-{^QCLq=*J<X{_@chU_4D%x{||Ur$FIAy
z?T@eH2y97(?fHLbxy4|jEKtBcbG;0&K;pFL=bOB?@31Wu!gdV3Fm~GWx_x`K(&sqY
z4cpl?jETqb^LnvW+q?CtE@zpd>(}bGi0x^9cNVsHLc8xe6i;f0b-KgReTc3B$HU{F
jo;x{qhq<^|tEeWIbFP<lY56*y{lYIPjX!Z2xJdC|t|{d6

literal 0
HcmV?d00001

diff --git a/tests.c b/tests.c
new file mode 100644
index 0000000..f45de54
--- /dev/null
+++ b/tests.c
@@ -0,0 +1,77 @@
+//
+// Created by frederik on 1/11/25.
+//
+
+#include "MMS24-25.h"
+#include <stdio.h>
+#include <stdlib.h>
+#include <stdarg.h>
+#include <math.h>
+
+void info(char* fmt,...) {
+    printf("[\x1B[37;1m INFO \x1B[0m]\t\x1B[34m'");
+
+    va_list args;
+    va_start(args, fmt);
+    vprintf(fmt, args);
+    va_end(args);
+
+    printf("'\x1B[0m\n");
+}
+
+void success(char* fmt,...) {
+    printf("[\x1B[32;1m SUCCESS \x1B[0m]\t\x1B[34m'");
+
+    va_list args;
+    va_start(args, fmt);
+    vprintf(fmt, args);
+    va_end(args);
+
+    printf("'\x1B[0m\n");
+}
+
+void error(const char* fmt, ...) {
+    va_list args;
+
+    printf("<\x1B[31;1m ERROR \x1B[0m: \x1B[34m'");
+    va_start(args, fmt);
+    vprintf(fmt, args);
+    va_end(args);
+    printf("'\x1B[0m >\n");
+    exit(-1);
+}
+
+void testInterpolation() {
+    info("testInterpolation();");
+
+    double dRangeMin = -10;
+    double dRangeMax = 10;
+    double dRangeStep = 0.5;
+
+    for (double dX1=dRangeMin;dX1<dRangeMax;dX1+=dRangeStep) {
+        for (double dX2=dRangeMin;dX2<dRangeMax;dX2+=dRangeStep) {
+            for (double dY1=dRangeMin;dY1<dRangeMax;dY1+=dRangeStep) {
+                for (double dY2=dRangeMin;dY2<dRangeMax;dY2+=dRangeStep) {
+                    for (double dX=dRangeMin;dX<dRangeMax;dX+=dRangeStep) {
+                        if (dX1 == dX2 || dX1 == dX || dX2 == dX)
+                            continue;
+
+                        double dR1 = interpolateLine(dX1, dY1, dX2, dY2, dX);
+                        double dR2 = interpolateLine(dX1, dY1, dX, dR1, dX2);
+
+                        if (round(dR2*1e6) != round(dY2*1e6))
+                            error("Failed assertion for: (%lf, %lf), (%lf, %lf) -> (%lf, %lf); %lf", dX1, dY1, dX2, dY2, dX, dR1, dR2);
+                    }
+                }
+            }
+        }
+    }
+
+    success("testInterpolation();");
+}
+
+int main (int iArgC, char** ppcArgV){
+    testInterpolation();
+
+    return 0;
+}