From 3f36a2a4c3e9c6b0754ff08acfbbc7d5e1137c18 Mon Sep 17 00:00:00 2001
From: Jimbles <j.avery.10@ucl.ac.uk>
Date: Tue, 29 Jun 2021 11:06:52 +0100
Subject: [PATCH 1/3] Removed windows specific file separators

---
 example/KHU_example.m       |   4 ++--
 example/ScouseTom_example.m |   6 +++---
 example/SwissTom_example.m  |   2 +-
 src/matlab/KHU_Load.m       |   2 +-
 test/S3b_MF1_log.mat        | Bin 111210 -> 0 bytes
 5 files changed, 7 insertions(+), 7 deletions(-)
 delete mode 100644 test/S3b_MF1_log.mat

diff --git a/example/KHU_example.m b/example/KHU_example.m
index 46aea67..f8a02ab 100644
--- a/example/KHU_example.m
+++ b/example/KHU_example.m
@@ -7,11 +7,11 @@
 
 %load boundary voltages of baseline recording, cleaning the data and
 %plotting
-[BV_baseline, KHUs_baseline]=KHU_Load('..\resources\data\KHU\B1',1,1);
+[BV_baseline, KHUs_baseline]=KHU_Load('../resources/data/KHU/B1',1,1);
 
 %% load Perturbation
 
 %load boundary voltages of baseline recording, cleaning the data but not
 %plotting
-[BV_perturbation, KHUs_perturbation]=KHU_Load('..\resources\data\KHU\P1',1,0);
+[BV_perturbation, KHUs_perturbation]=KHU_Load('../resources/data/KHU/P1',1,0);
 
diff --git a/example/ScouseTom_example.m b/example/ScouseTom_example.m
index 5aa028f..38bdf4f 100644
--- a/example/ScouseTom_example.m
+++ b/example/ScouseTom_example.m
@@ -7,14 +7,14 @@
 
 %load boundary voltages of resistor phantom recording. The log files
 %generated by the ScouseTom ssytem are located in same directory
-[BVstruc_TD]=ScouseTom_Load('..\resources\data\ScouseTom\forreal_2_clean.bdf');
+[BVstruc_TD]=ScouseTom_Load('../resources/data/ScouseTom/forreal_2_clean.bdf');
 
 %% load Multifrequency data
 
 %load boundary voltages of resistor phantom recording using multiple
 %frequecies. The log files generated by the ScouseTom system is located in
 %the same directory
-[BVstruc_MF]=ScouseTom_Load('..\resources\data\ScouseTom\MultiFreq.bdf');
+[BVstruc_MF]=ScouseTom_Load('../resources/data/ScouseTom/MultiFreq.bdf');
 
 %% load contact impedance check
 
@@ -22,4 +22,4 @@
 % electrodes. This code is a bit different as we want to estimate the
 % contact impedance and plot it quickly. The protocol is also independent
 % of that set by the ExpSetup
-[Zout]=ScouseTom_Load('..\resources\data\ScouseTom\S1a_Z4.bdf');
+[Zout]=ScouseTom_Load('../resources/data/ScouseTom/S1a_Z4.bdf');
diff --git a/example/SwissTom_example.m b/example/SwissTom_example.m
index 4ad9676..53e8d02 100644
--- a/example/SwissTom_example.m
+++ b/example/SwissTom_example.m
@@ -7,4 +7,4 @@
 
 %load boundary voltages of resistor phantom recording. Removing injected
 %channels, saturated channels, and plotting data
-[STout]=SwissTom_Load('..\resources\data\Swisstom\r-phantom.eit',1,1,1);
+[STout]=SwissTom_Load('../resources/data/Swisstom/r-phantom.eit',1,1,1);
diff --git a/src/matlab/KHU_Load.m b/src/matlab/KHU_Load.m
index 65a4fe4..6186da8 100644
--- a/src/matlab/KHU_Load.m
+++ b/src/matlab/KHU_Load.m
@@ -359,7 +359,7 @@
 KHU.info=info;
 
 %get the name for the new files
-k=strfind(dirname,filesep);
+k=strfind(dirname,'/');
 newnamestr=dirname(k(end)+1:end);
 
 
diff --git a/test/S3b_MF1_log.mat b/test/S3b_MF1_log.mat
deleted file mode 100644
index 220b8884ce6345c1da99b6d029b9428c09ddb399..0000000000000000000000000000000000000000
GIT binary patch
literal 0
HcmV?d00001

literal 111210
zcmeI53y>XGdFRi}NcQa52H7~E6*gS1!5IX$XGRZ?6~K2!507L^8cVV*t8k=wT)jNf
zjAusnNP$98fXGg=5hiSm2|+et5n2{hW~-ug!mjWZtXEVPuS2$kw~`XBB&9eTR;whn
ziYps$pZ@;8nZ7;KjbyC84e3_tetjN&PWL&#@B7X<{g~BV8}91rF5P<bg3{`)4L5Dw
zF;snjsduPveCzP&t|g`3rT48|vv}eAOG`(qedE=^((vvjr4^$)N~`wlE_Ey{wJ%!I
zzHmvf)6v$xxYV<3#iG*K!1n5{zEb<mZEdBSOEYbbt#RdQb_-db7I-x1v#VO>wS=`d
zH@jJGj-5!OiPyC4e5N6cAZ)45ZPAXUI(JavBRamp`G3yMRs7W|Kb~*Hy5-AVkdyL}
zI?SD(ZF`hMd?EjNv_Q}Du3mqo8+ExolFy4X1>N$Ak@eN_JtKq%1FPltZy2p^9lOD~
z6{FRM){YKVN0p&>d*4_!5&_53f`G_bHT=#Izn%A40S9$}@!Ocl?}=v1Pq`=FX5k|`
zk3#%xw*dTZ0WireZ?)@gaF@-<*{N1b58;0uO|*W^+GP!Hwkzq_^Y_~hr%~Eu;mU_i
z3R72Vs?FD<@OIf@otE3a+43d*b}C{2U{Vf?$F841!PB#6p4#A=!~Q|v8=SR3L!-4M
z?RmA%EFNwOJfM%ex^#O`lHV$B;T3M-YiifoD!)+$HoHS%t|mO{{uPxlB2nI}^dBet
zlrKTLK}$ce&E%F{%KuO(uhH52YpGOr-FKyH3--L{zA8v*cEKPr7*B2<80s4v^YL#9
z{$CcR`;c>k!+ZLN0*gF#+~whMJ&Ss~HuPBje-Y%P{cedHbk!dg7W^3Bz6{%1FDGSY
z`Mjb#-#&9Xv8S{%avd=<%d?6u-8*wSvFO3^Hxo0nY*qB(H_e<*>>TYh+sKXR*DdSL
zx4X7tC;x8Ag_8G6?l5^`RN=iQD<4z*0h6acVbVQnvh?dFE5B*-#N#ICe#YdXCzbw%
zQ09$q`sM88c-jcT?lAm$r?V{{zIeNHK2*IXcUCG$zXS#=+|J%^a<0O2B*}-2=TAD0
zeDr%+<uyxQC&}^gc<SeTO3Hs@?6~;~-z-VJ?^n1@vQ@Hjf%10c`;g=>NPa-_HpwND
zAC$aRlH({RUT>%JEt2ezTrAlx$?=qj{>0-K#PU;~m7iDsn?+-i-^E|V<;C+q$N8;k
z(xX^^6zvjg7unsswM%@QJy$?Xw1wRo^y{2(Y~in_F^sb7N*j+{7371f$??d_-CKvP
z$e-vJ`(C$kWy^o4gcoq_S4_5k*5u*wTIq%?T&YTG{8)Lw!u#)+Wc;~d=6sb+%6E_C
zx=`kgBl@L&cE^QV(JrQ+Q}(!M7vBQST8GMeyI5~I*R^XT-1X~key4YIczk$Zc!>1n
zL)C%JyZ7uOymWZ?c;CSI((ToO2UqUiSsfVPuwz#+Z(LpN8{0El-BsOf5$nfy?AjCD
zt{5KOFuG&gw(6*5bJOD2ZrwUo9goEK?im`Z-nV0Le7mo6D*wI1L*sqhf~vx8U3Db5
z@iFGQ`vy1rb5yijI=pLSXh+`w`EUn{-%SF3V+#kq!P4fx&F+lKJviHTsQ2LY7Jk9h
z$n`4UeikYej*Gd9IDADsr8ay*!P30i_;~1$B6@x<o>Cj0P_XUV+W2^Aiz4Rw1%{OE
zh|-$>WlE+Xsh;-!_&YxF*LNSERq<z7*n0zY%C6abFR$&=Dn)D19j=B3i?FdJV@yAv
zm>PFB`;K;%@OS#UpJ|@Y&xS7haJ}2%O4k{lmS69?GVnrL@2l<_8Q<#@{8`)i+EsGy
zQ_|K^>o%@B_W6oG&w96`_u$dZ?jOxczJFA*+*@=|vsPhx6rOvt?VM7$wA6Oy(C~^c
z#7}1n)YCs%i`hC-7J8FD->%gg?r)t@w;nh1l-S8HD5ZB{VwN!Bypqq}mY(mpl25g#
z=i8#>&iAJ0t5<40`TE0^S2n)yPxe;)wdaj5+NB~p8ee~eM$BwH9QFIr^`<nxzrc04
zwZTaJW;gC8g8#3#U1kLba$W=5`@%;auGF}rA!zA{|K72CeDh#+V0e)E%$(pt{PE6}
z!e7c^T=RO@%V^Jf`2W?nKH7D9^QSBRg7V5)ai6ES&l4(~>j*h@Zn5M~-+i8o@#9H$
z^7Yyop3~aC%v$rf4A@k87xuiu*!bw40W-l?kbdiF$Lti0hrHUcgKF=%zr>HQ{1gvg
z{QbH^VRKAAClfrCV(pV8dzJor<rW>EF7w(i_3*paANE{8(Jp?4Bzk|@S|@MqV)Lq5
zldBX*-f8EV^xLY0(hW&DEPk$%o6G6>PAU2LqV#-uPp!Y77_Rs>${Szw8~1_ks_x$A
z@5B(+ZJf#Qnm&K`$GHu|<9$OU2H~~AU7sDh`-V7<a{Ya#HNowq;NH>ZrF#d4s$-Ti
z72d2nLR0a}s#|yLu5RAtvrgR)x^Krwbu@UUe0<05ZO$#Lj)nJt{3)y-@6%naAU=4`
zeEs<7<oRx1AKak|%Hh5d;bqnSJ==neVPQe~;BLFHt&;e|9y_hq{R?-!34YkD(A4}e
z9!?sMO`o@#<~O;L>v3Ify}QS4cB|bAx4|t7_SOe-x!dIYoYrrzhr?nXTYJBE?BzfH
z$C14^cD>&Ayx8Y)FJ>GMcN-J9yEdd?io1CD!ntb;OKJ+9!|n*4LB1(qZ^VrT(udmJ
z;(*8Y;NKFrDEw`63&Z?BYy1`K1;u)SsZ>YB<YLL6?ezleOzWhbv61D^!u{7Yd;62$
ze{BiZgqtUy6EXkI^W+VUMv<D;#lx(toYjD#*Rlt17}<T@srOf!!{y^9*Es6~f%k8B
zA97xHgtFb;s<4k;6j<yw#oso$Z!yu|j=+X@hWQqSzm(<6SRBSIjHFJ_W9JafNeNcV
zZQ;=i!e2W-&f3>tF;;_>Vdvt#j8Es+j&BdH%EwZJ9q0X$<+XU)<Ht+R#Ll!sXVanf
z=$Nd{%I|P?EgeCN`}!8BedGN=y}o7bY`F@%wvNg3r@!LoVdpbdu7%rGKk@C8=V`}V
z+qCPfyvMdyfdd~y$?PqhjGN+R@nqbslY5<$ds8g=nkc6|93fESf{;F|eX0zr-|{R}
z+a^h0i$DIy$@ZX)sWpAipNE|hW%~ZHIN#&8w>XHS1${rV?-^r%L8Rwl!1pOMPSvzP
zJ>+SI{WkSJXVIx}dr;~W`%}5CS3AY{t*XVJZT`mn@#imC?+B)#;nq7}I|*V!@n_zb
zkU`OhNWO`knI-eRgYTI+{ih`sznnDhOj`HF-(2|kpLu85Ee%Hb9}Z@jD}%Y_O@W*U
zW&-P7HMn!SC)gVc((=dO)P0EQ??0yb{|0CNd&}gwr`CUi8SNPl<G&BuVUy2;PWRti
zLgbo)xqG}tH~-j~{^r1nEDrj@&$?_z=40A}eLox99NgL7<@&-q+j*Oh&1@HDyTf5?
zIJe~%?Ieg3itk*Ikb!i`k$0!(W4)fc6-l{y&g)oLQjYq@Pk#3!_B?W$9;Ceoe5ErD
zXmD=|o;jS$bBJ4aREGxr`kjB=*AG$O68!hNz`Wp^VBc@ZHwF8p6=77E@5)Hn4~3)Q
zrXcLkiTr+m>+2xwht&S{W5Rx&r6m&fY2N;k-W=?ET_^0zerqJ`kM(ts?jOaD*6+_P
zTAwc~TAshQXxP7c?8~4$fA0RA3Hx)4hJ72|v2^WSp)K*5gAKdl%k~@VuB_i``Vh8T
zG@T_|--dqk=k7O-XMH!EGz9VfI{g-3@Ac>BH=WM3^<UrSe)H|(H`aj(`!@C4)H<+*
zeP8$u|1_0eyH3n``F8c2&lg=crhMOzeoI<M=COLNmvv>@g?i%a%+$;Eu<lHMkWRa@
z?#yfFq{r8#X;028zD`ZM(EhAjqesNEj_vi%Z}d0sf3V)ib@M(5uLaN#oImYIJnxCT
z&5958#qTD)FT(vMu0Q_%2(JZD|5Sg~zE8sI4b(d=pYA5SUFXGhPUWk8-{d;wqg|)+
z)xLisb2-T~@1wj!$I~wKC%5U~C;pyFNypLtytnc$#lvs>Jr*Uq^8i9bA&V9$TA*lw
zq6LZ;C|aOsfuaS97ARWaXQl;M&t)Ad>3Lq(1>*AJpXari=;ZsqIelM~^{n`J5vSI(
zF7EmA;JC2qE_OYOc`ee<w^`5eb8tV;Bkbp-e!fW9&v~O^Ki7_~Kl!<?FPHp&F71yY
z?B~ut4`Dw?_H_{Ub74QfCG6+C(e)}nw~f~C=d^y_N`60=jh=rg2>bad>3+@%pHt5P
zv<9(3yuQoGJRC0L;S0at&=!_!&l~ww?N!0uh22`+=6Rz<!JP<ypKnQckIf&DzxD0#
zGJ<w%3jR-hmznzl$NkNcVCymA!uOJchXD8<Gg8T}`^(m{hZIr%l7)|HfuZI5mT*D~
zLDwv_^i#^e?vHHejKb%(TKajyPuULkr{yY3XJKjHn1xG<?mlYyO`Dv1`C3b8RsB1E
z$-+I#KkpgK&+}>>*IPob@^}A=g-4Ws{x>cEeuZ1!V>=Vdf9@YF{RyQvep~sqK(^|$
zmcV*j%X=+cvTUyVBer9-Id`DT(kn`T^?>d4D!fD&vPI!%9<`kjh0m?A<M%6k@>8~R
zNa3URSo&dw4?JNzhqaJ<?(Z#pOyPN(EkFIs_qM`?<sj=BUFYM!X8BZXaJR>{e#rsJ
zLCLD*R>^IW+a)P&j>6dA`1fnC<h!harSab*jeb|u+aAifMIVI5esKTfQhvulryh#w
zcA+zjf4$#Th&~sT^c_}j%enr6zTLyScMS9mmH2c6EgL<{sqd?v)q$h`qc%MrIxRiu
zxQ&|PIV}8wIGOj=^n6E@{PKU9xg27jP`$6bFTFgzkMsDG>G_T+x$(bJ!%rsR?2OJc
zK8|ADR>>|N_TwRcO5Wc#2mfaWcR+mDKX2v3?g9Vf{iVf6<9vLxvnNHH+?mFm{svjF
zz5G8-zH+-=Z24geKlx{hUuy9$O;~u{4^2M3!r~wMsD)epT=A<czI(rgpZiObn^r6S
z;}+icSBk&e;$Quwg**P2$!j(!{xJ)`@*~B|NI&x%7JmBwm^``3@?Z053vc>6lLsEK
z_>-SeyqwB`KE;35;vbh*EE%-;1HWVOtC}plWV_@SEWYtF3!mGm_%B-gv&|Mhx=ZO_
zviQ*}EZi?Aa`cGe->mp?9rsm>pMQ;o=RKnIuPgqo7JhNB;=f_>FV9u{M@_zT+~S{n
zyM-UyZ}Ra!Q2aX;|8d3tCyQTKws7~anq2jq;%~I@tB)%FJBpug;b(r^B)Zaalj1*T
z^4afN{^xG7@Xp^={P!$=Ux$UuPnjHjQSpl{-2WBD%l2SfupRS%VEH>fX!&0Jnn@oU
zwX?S^zp(FB`=JXZ3{;1Ptn@Y$7x|0c7ye~t|3do+;h^2+2c>)5AB)?=r1(TT*<^Nd
z%Yb!Y)J_uTQ?78Ou<!L5r(oaP-3B*U-5PvluF9}8h_!r8jkAtu?BY+!pQbnAq4^DD
z4SyH<z4V#*^W$)T9*n<{jHi5EY1^5HYEzTiZ}P-zvlpSF1&S6ZTA*lwq6LZ;C|aOs
zfuaS97WgS@fr(%Io9}-sxVh5Ycw)uRz4WompIY?nhyH2lEx+){r~c{i<qw}|c+<%K
zrt40-%R4WA_j}j3{@|nkK3Gez6{>p9bXSFckxDbj*ZGCoOSO_p`yABLLVN{~85t`N
z@$kjpNj((uIrU|va>amszvRu5^Cg{RNph~_yCvHsua%r5*)GYO*f%QND%mV~T_~d~
zGS*+VvU17(f|#xJ-!kafYUjB0OFg;0ui!RZ$@V=o?JwUsID<dqD7bf7u^xrzvcm8S
zckmnM9FE}^eTegGTt~UJ{>p~it#*;YZ~nNg7Tl(A3x9A5U)V*shF|!`4#Ov0V_)GP
z&fyx4;oUfna%=sS4Y%8K%<Wi?xux0qw9_xCg=>c&`*C;7cj&mCk`GGic1G|9rfs_=
zhb2cOSuh?|m~qKmg%S4V^kr^Nqs4Gnk@RVDmc?+_aay_}S|uCzB;~O9xvF>HpQh)t
z?@>5+vYRWJIlpO#bKDJ`IoBH>%Q&Qd<3sEcb_e@~UBOP3D|YqR9qbTx2Rnwn!_Hv`
zosP%P*4nkGy)r$`#$KT(dGxdUb8N5b&uyG5xP#L;XK)PHaqi(4zVWy4jb6aD@gL>Z
z`YRi5vFCZ@_Fv_gTl`j={0C3Bm~ET>Taw@6>6hs_ET%^Zqx+I_SbV4IedXo!e66bY
zz{+&L)vIE7Dkr_vH!ijQ`j7si*XS>Ln`94W*5j#pCHjt?a20Dc{3CwL>_k-mjr(l$
zpYdBB_5XKrtpD}r7M|c0PUHRx-b&&R&fy#W;Tc}b%8!4<pTV{HuPC?HU)gXwV#*Z!
zHc#&P9r~r7>v2bO%q@Pax6zu7rz~q4F%qcAfAAD$B&lY1yFxGha!vHKz(&=(>CCjG
znus3NyQDYWZ^1$L$7>G1<N-LV-+T*wM?cFIyJ++j`+(k~=Zrh!`i@<|FQMny6YL0f
z!|YB}|Bd@>^nY}+hx4cZzn^3MuRphN6z3Ao@IUYeuW$#~#^2O<3x3f9I7T1f-Sj5P
zt@T$n+>R%3OLU)h`X%eVZa1d$;JjH*dsy+0NKQ!Zm1O$TrEs_8QpshK%OzJx-XYl|
zxl(eK<eidtNv@V$BZ;?&uQS4ZJ;w##&aqu!TrkpT4awUSYSw9D%PfoHZHs9c6md}U
zadr|iPN?3Kf0K|w(c4w;xqIp9ne!h~^sA4}oPMcd;jDh+0_+_24||He#ExQb89%s+
zwF&mNtT1*JdyM^N{83Uq>}+coj@nn_J{$XrUCQIU;E!{x|Mlk<PT|csat>bM3*O)t
z-r*aL;ToRd6}^COcs9PH+**HS!!7<ZkKFz+$K0CVy4)RVbm7KZeHNPjyI=l;xBjN*
zu$UuCIEo`l$oPcnUGldH85Dg`_4a=_-ESRJe%?k*DwQ{^d0R0l--X84U;oi_{0w@H
zKBK4D2lN=bfPP~Kupgz0wE=bjyMX<`?v%rPQT;dWv(bOX%X!rQALUs8>(4FR!4X`+
zn{gE|81fbLuOwW;Gu*>9oWec&0q^i_yhpjU{>p~i?%={~Yi1{u^h-VGTmM*aoAP`a
zJi;Hml`Gbf@CBFf4ClsO@Xa3OhfDMUUg0^;z3EAmTkEfExLsx!8T{st+m<f>rg>+0
zhWY%W$7R8HAglqOygYq8b8435du%8@UvT3l+?IUNWYX_h)yvlslJcEay?jj}DPObd
zeew5_a$K5t>|_1b3D9Ts7`;XBZ>-pr$MqOJNAIyCW#z|SU?<{wkDW2Q6V-p?J{$dC
zVHX$t=1>3MmSg>|KeuoRr|=7B@EYe7zTpzT1g~&e631mp{3D#3-bA^z{>p~io`6`}
z^2hCa3T{)zTX0vdSpUHv9K$s{!Yf?DAAG|*oTCTuYMf4u&y4FRx7J_TaJ$MbI{2;S
zmirEKnk>VpzG^g?J-##BBKq~!qj)mmR{8bE{w*;x%ektz<?8hD34ihM)jjF?+LV0s
zt25VksbaCq^_*`lDD*!?|Iu&s6#Zo!g*`xT@k8i6_5l6H4wNg#6m|o9fxW>#n4O90
zzx7u(`hV9X__g|<2X4CxZd160D|mumID}*PhEq6%Ke&$bZd^JC?{Th;|0uWCU)gZG
zCdb^a$uYP1^-lTM=iZj?H(TY`pZQ37zFsBwe`B&^X56W$UcPQPW2PXoP4&L0uOBAm
z)5jy-v%jB|<I==qAM3YHfL`LC&{ymLdX1lo>pl8ju2>6UH?SAj59|x}20LK(B&z?$
zeKz{vJK3$Zev9Zn?et3yNDfL?CAUh(zfAl_zh^n7dVX4-*IKLdrJe6qc%9^W$qkbC
zNZu=XpX5f#`z1F?wo7(ME|6R(xkz%c<gJpn`s#we%WRv`=3JjTx?gfYa!|4=xm9wT
z<aWs&k~<|IlpK=WCAnL2SaL-2A<0q6G0AbsJ(3SgVom4jIuO{Edd?R<oMZn*`_#|w
zmMhj*N%j@{iCxBSw&^(RC-xOP%y<lYi#=!DR?_iir=xb)`YRi|i~h&g3+R`6a=WwO
zHf6p5cW?!F<r=P>gTpw-PWj*yp5YkI8E3#NdNO^yVEvU1xA>_zw~U9kXk7ZFzM}0C
zwqA`(AJbR16SBDp8Tirvn3Tifx2xWcci|Cc&c9LhmOq)EZ@-cs*Vn&iF71tv&4d0=
z6#Ac{|L8q>iaw*K=rwweeq$%F3-}}KNx5P)VF&Pg%qOr1W(T5v%la!D{l9k-`P%tH
z9=QEP!EFk+_z$>(Gq{AuILGh~$8Z|w65in%y+BWl`zW{8U)gZW{3-6Y^1$t<bIh&z
z_1W&A{Oc<}pYAs&<X<2E*Xj8VD|yv-CJoArJBL*7I{h_`88Zcu`&BQ0Z6hgPz2f1l
ze)BE#8@<F1pvU+n>;iTGy{`2?!55$^)|U7?>;v`$y~iGy9f|6{ai5L;vtE)%{eLRQ
z`d@!;;R)X03SQ&<!8M%Yzu*~;(Fgd%Kf*hDVfql|*7_?OZaYlLg5Ug&x4u?zn=;;l
zH#ovR!W|r9Z{Zsb;S!$V7~bI<UgLUUdK2Z=`YRi5S$~dm%R0gHjRy97{Uy>#%ws#p
zXIaFH)%5YiX&uL3CyNy`8GJ$Y^4H5I6Jk3jR4;$sES5YOJfppNzdIR~$qsh8p6i5v
zSm=L>{-dwxCwhy%qu<y8^toIyfUpDjKkNYZ1OHc2e(VeO#_T~<|E<5W(f`GE5y5Z%
z{MN~W+Z1l`Bk&7faZceH9^o5a<2;tdDcr+1T%#Apb(CA{uWY#GzG$3V^IMm@W15fh
z*BU1=kL?`Re3ZZD7)zcEJ}<ZS<bRutnsMj2>fO|y?gx*k-bN<aGnVGu3o7PW{gugt
zZ0>~W<*!gC<-5@MJov3sh5o1LKl+P)<CoA=^cH<bZ?OaTAM6GCjy=GxU}vx|__K1w
zvlP{T>#uC|-)iw^kUzim7X`N|+`=Jz!DqSRd(%0%ggbbHKlp@KxJMu09z8%m;N19+
za%=td#&g>eTn%18f2sWzp4Z7Ew|`r3o5C%;l`BEF1{-!1{@@zk;2SRC89w12PT>>&
z;T`^s`zW{8U)gZWb7gttwxK(}?<q5ms^9&DIG1o7=MJ9X5<3fj@Cfhair0dv_qyO6
z&f(a2k8*4Ml?}H%|CmQ^-&}B;VsGIPF5wD(;V{lI947ILeXiwmievPpT=8^AxwZbv
zhFhM~%_Fzh72Kw9i=Bi^_=RWuJp5u`;nM8o6pwIR%cXO04#&0n5#`qUD;sWko;;7-
z&M&x4;TA3#ug3X>H~7R3!!;bkKRm-Z{Kk1TI~%+KHO<~ye`Uih?>*#^+s=a96mH=W
z`v+%m3}5gJm)K?OD_p`cT*EK?!#Dh!-bC%K^;b6B@_tSpxm{6ko5C#|!Y|yxA-uvV
z9Kt#F8ZO}x?%^6<(F3?PJ&AH_{gn;3yhoNtZr2ywrf>^Ca0q8`Xk1OL<H9NY!72R1
zIsC(GoO|@d_>OXG{gn;3yziJtZu<&uQ@D-u2glgUIEUC_c!pQ_gh%W$oG0mp@f_vW
z`YRi5c`rAQ+zuDqrf`d0g-`f{Pq@RbVmILv&au<z0rnQY;TV3=3*$b@t@T$n-17c-
zoLinxJT%LN+kpqT(`MoL?{P&O{jGRPZTLu|g_r!_+W2_r3Dw)bCB5DgI`HIY;w99E
zkEz~G=WF9Hb%^^D^?W|@BZdB_=s$Xno}$O-H+qU)K(DbE_$~Ardw{(_@5>cq4SRzf
zG5ZqLf9tPo^naos4QOU3l)Ue>UC&Xy#0|EZTvLHVQqxt>O~uW$x1aZzZcNwf5!L(Z
zqv`pcQ1Z+Fd8XRxA6I<$#QK2szc>23Cy$6<INq!9KFKc0Zpo#R%Osaeu8_P#vPW{I
z<SNNKCGV14ExAT=tz@s{-ID7h*Gq1Yyhrk0$@?TXO5QKINwQtCLvn%SLdivvizRQB
z?3BDsvQ4sY!1||OazJuWvMRY%a+~CK$sLkAB_EU=lH4V^TXI-(MDiiYQOPmMamhWB
z4@*8GIU%`Ma-U?EWVhr}$z_tuC09t^A=x9jQgW5#osxG+u9jRQxmL1Q@@~m>lItZm
zNZuoPuO!xOZp>8WVchUoF>Xj1H{iz^2QYq!j~5sRlq=T7_<hC^%r6*kF#cd1!T7}K
zc*Ysi=MC0h*^C>oukm?99=Ls?;5LQZIFE1zhj0z2NnFDr{Kxr)Z^p;9JV(d7)?e9h
z%Qzs9-2R(_+Z1l$4X)r4`v`||p5Yj~nq;5h6ppdi=!5Ya<<|Nu8*Uk|=8@ZP6x^n8
z3y*Mw{e?sLg*!Md$=|~z_7g7SdI0CxarDKwk8*4Ml?}JdtMbV0w+e1kxP>RU!+ydm
zJjZ#1Z}^34c!hg7hi~jT`eK|%xwZbvhFj+Uac*s0Y?eFU=)!GX3o|$i$A8bvvY77S
z^n9!jpZik0gjl%A!lPHDmlr!KcfmZ&6TVmIe~SL2ujnayi(aFrar=PZLGQ5>=r{JF
zTrnuIC)f?_j@gB%{#$=#qyMb;#q~cA-2Q37Z3?$>zTg+m;1@pO6`tV`UeN>i$6vuU
zyy73H&l9Y_vf-9>={$1#N{+e3ub1Rs=U<cVU+3ul)VxR1^K~k@`)lbgxK;H&_m=eX
z(o3m-D)y1vg$=GXJgEkK-{@+cB3Bs?H@UgtJH1zzZr7XV6}RvTx9~N!&+H)9AKRz`
zo86(X#3px8C8O8qCH@ILN1w3=*az$edK|X{<%(5+{fOHE?1$NjsNXW~v(bO<L*-Hb
z&*xbG>(4D*;#c4hj^P!~@K5jym++3if@^q1Kj0g#O+TXCT7PB3E%!_F$n9lIi}#<C
z?JZowYZ8y}2ItsgIE8C?g-bX`FW?uRO+TXCT7PB3E%)`~+%lf%k>A=k507VFB>vl}
z@x;!LrRU>uxR-d~E?$try{h-gx2KnvaCmN_vCFKFGz8-Vzcn=2bv5ctkxgz+I6jyY
z#5ad&mGF1^xWUTt`DfaqJiVpR{|gvzp`Yk8dW*i}$FL9RId%ZQRjycnpx@XR>;rbB
zq<m&~ruSDi`p<Jb&}QfEze{7|)OvFJu7cYXZt+)e3XgCIr}!H<g+I6~i$nN^fAj%P
z(F3?Po}=7ae`Uih&zt3u+gozXEq;AO{&ihB-M?;;U+@0a^nA>h=YI<eGP9gey|3Jm
zUf#?{=O-5X$aQDGPJn)*@Awz=9(|4LJ^EX&SZiV*um{)`><RWIZXe7pMD^de&qn`w
zzBG^e|H0hr|K-5~l=c4nADB$~Ju7|f_+Yxeo|c_>@oVY%cusWHXi|+9->ma{?xpm6
z=XHK3Ka^f?NyY3OOE2$2Q|qt)=q>tb`sf_`jy{(wc1_rYxZPkpfPOpWXB>bX!7gAw
z%sxc*-?-04|9PG~kNUqV$NFD?Zs83+;R*iW6rSUp!!P~`F5w)1h#t($bCg@_uWY#G
zy@xz<yRqOlWxfTMajxJPezC`J3GZ+Tr`T8chjX}wbL_9_LzG+VuWY#G{hT~<yS3mp
zg<E)pGq{8=_=P_>$Nt9o#h$_`e8MyO0q3R{QEsijvf-BZ$l~1EJb1P{rg3Zesr2#V
zDJ8G_?=E$FUgObCi_*(Gu6mz+Iz8WMC6B(CT5a7D8Bf%6op8L+{}la4U(s{)zFaYA
z&{y;vJAhut{Tg-zJAvNg->@&(hv|OH`YRj#=e^Om{^x<)UoN;!;THbloWgUQLwJMl
zIG6BkJn}ussqeJIJ$hn(GwQdjzp~+$_nYI~n%}z29hRSIxhdU$9Fm`T=5y)!jw|`e
zlV%l?evhc$r*BR#?+Mkr>htOOPAK`=?<UoFY2x#s|G!!2e~SL2zvwf1kA9-h<%$tv
z`sy5Z13kx1;K$<r2|tD(#E+Rhi0Z%fS2p_3_YC6tp9gL~UvQhkEu6v|JmH_>yuu%x
z!!<m@FWjRSaE{-CcjG(Ct@T$n-16OxIJfw%Qln|>lYeG1>37a7OL=8!y8q}@0na{?
zRFcKFs@{&%>G|fW-UG|i>#eAmmnPE7d!tkHp#Mh;{ZG+<^c+1#Z_5>90sT#~6W9Ut
z8oy?GJ>9=yKWgnoRR68NveAFOvlG|<JaBuw;5LO@_=I2hgfBQQiDUc{Jmcr!6#ns3
za15t#j-NE%qug45Wy3AsgNk#D->S%ewEVeganf&_{KwO)(*4Iq6;S?IQb`uyqk5nI
zi}ZY*s&~m<>Gf_=G2Q#q%X_0!^PvCVDfB-@|ItJA5`D%V#QjmZVvJxP;(Cu=K;Pr`
z1G|8oF?$l#f9tPo^q=p)#q~cA+<w2{HicU_!oR>BJmat67{1{de&N|T4A>UmaemPY
z{F`wf<<|Nu8*cedVVqn1))x7XeSc-zn)KT%|FP-rbpOG0=cP|1m1ObTRqr)_ot|%_
z>OHqEz25yQ=G9N8m-j}e=0X2|Sm=L>{-fvkAM_l5gPx-2*a`F+dw?CmUSJ=}6=N!H
z7qByC52E^S{gsXW^Ig)o{^x<)|6OpK!YzIz&KDfRHJrgSehHr88(!fbzTq0q;oI~e
z%B}TRHr(?4+&H)Rt^M*Jul&fgFX?wg{^OdB>Hgyh74Z14CzWLJ2UYKZv+4OJRPWLI
z)9XE~VxIYp^zz>5)I8|_|10!AMgPkc>+`stqQ6P{j{c(W_$BNB_5%BXzr)^O2h7ey
z_22p{8~x{d*>U~P1GiT!E8dq+eqSEG;1llfTlgFJgg<zMSNMnbIG=EiemLoi`MoH&
z)?e9h%Xi)L$nD&M+Z1~Xmv9P?a0st(Ju`Pfzl%dShF^G%bBw;2o<zB|{>p}1{{BFm
zTN`g(?hb1_bD%GMJoCJ>t^c>RlOT>aTEggS2^kc9MD@-aNU!$=<v;n^^zu%q-p9kQ
zC?w}Srh1nQrq`QXLTySO{ML<y{-@}Ft^aTieZ?N2&*&?9jDBMuup`(3><IP(e~29^
zSBw-JcSZHz`YRj#=kJZg^*;~XE-biB;TC@2EzT>Pu^s@2@E+$EUf~;_(Fgc9UY(on
zx2(Ui;g-LPlSgj5bIdL4D5o0D`n)_V-OoKQfB52#^n7QnIJe{r+DQ<nRd3^E>E*qk
zde7}l&v#z+ym&aF!b=p5eXQSi9?;8j#V!$j#~z@s=rMi@{l#8j2b{{m4&e8&2ly}S
z$n^P^ai5L;^LMcFsQ+totpD}r77pPKuHX%>;TQhk4!-eg@Ez9!IEQC=H_oHuE$gpr
zxaIG)<&oP53T{)zTkwdTgjcwPKe&cVIE>q8_=QV2g?ID>{V=_Va%=sS4Y&MV!8~&N
zV8Lw)w{QtZ*iZO^U-*VgIF9oP_t<H8h3`1`)9tPGS2o=8_b>Cv?L>~bwRys9cXpNw
zxBijzarj&%KlWwIpY+?T_XO9yDZRY&s`u50(({#+eDp|CjrEPM|9C6zf65hW5%d~;
z#U7yd*ah?*`+<GHU*Qk2E7$|<$#lPE+-IZz{Jq;e>i;KmtpD}r7OvnC{^ESWA)LZ{
zoNsu>PoW3!Y@9lW9vIKj@s{;hHr(=ene)i)Zx!68jJMzszTgqg;2C=f*Kh~-@Cw)1
zY50d@>@2)bx3|_`*>KC>507)pcw#Q!9hgY>>lG#Uf8D%Q(r=r_6Z7AiUS3J{zPL9%
zUyqVsJ(g5sedDpq_57~R(}n)0=s)_4-lE^=J@x_pMz664*o$(-NJFpDd+Y>ykKe)`
zn4O9GE$gpr^uKoi{b*(<lzHIx>jk$d+`?bEqK=b<H#o+R!7&`dH$1{M{K6;P!#P}=
zzC^jT{>p~ibvfqtn>pqdzh04leQbZaf8DO+(Pu?jg4ig(zV99B<@KoE@`3bxEMGkS
z2MHC{IU4&|zx#>hiZuaxi=N`g@N4LO+#aCE*b(do_5gjy9$+`H3uZr}`fuE4qyHOn
ztpDH5vHsVeTX>55A9#i@IEHKZ#Gk<@T*5K>ki@@nALZ8iD;sX_%`vw>D7Z}-Z^0`Z
zVn5+IZa3pR!ZBRJd)z+5GyJ0u)8`4+U)gZGF@an2LC*DQr(crng9=wAw@Pl4+%CC8
za;M~jl0%ZaBzH>=OO8lBBsnTMCOIy-NAh7wJ~PF8(0diXPqIt0TXL!7GRftVD<tob
z?2%k4xk~a*$-5+1ORkYzE7>b~x8yp>^^zMT?~%M$@;=FplJ`q)l5CgkkX#_SP;!yv
zV#!-2J0)+EG}?mxR!8?s4oD74RwcJeZj;<DxkGZN<b#q!lDi~#OAbqpNIoPvDmf-O
zF1bhYVaZ1%CnWbu?vw12?3P?AxlD4o<O<0<Bzq)RO7i#G?^Kw-gS}c|{x0rXh537<
zcPq@_aa^x3f7k6Eh536~_bJTZ&$(Y={+>p=!u*|r1q$>1_eBcxUF};H<~zH$Da`jv
z`=ro(&u~Ctz7tnfnD0hyQ<(4Z>`<8Rb3CXp-#yr+Fz-PRE6n?)4=K$1W@8HT9?u?y
zc|YP2g?V0jufjZs+NChhr7cyM=Y5td%ySfXD9ruol?rpe>Q05Z53pKc)>qdm%zDe+
z3Nw#guQ21<dlY88aG%21;rkWdB-t+6A-O<uq2waT#gey5c1qqR*(TX1!`&~*Y>o9`
znS(+UEl{*T(E|T$TY&j=J@4=QWif9{nYS@NWB%qU25U)@^#tZ!%+uy5oprOa!Z%8`
zN-_^*p2mERbqJdWM%Npxzp|OPF+NDTzuWAzJ;2V)^1KEBEvIHqC-zK}?OpQ$Mhv9X
z%2Qege&*k$=bNL$wrE^ZTT1;ydC>pAEA&4_|IyF5eq#r)52mNQC+8gYqFgb`uoLJt
z_5*##E?_sX4`vsl`fvS}js7zqi0gkIxNTfse4iuv`yBWoIE6F#f=4)nV>pCgIK+>k
z4{(fL;J@H~y5F+?%7$C!t$F13%>}n9_7+a#yul%y!72R0D?GwCT*E#3fPNUKQ{PiJ
zexvr*`YRi5SzpN`x7X*GTi$biLF31Li_^!2Cp3QS|4MqkGfH0fC+5YIeorxO?MyH4
zdDZ*s(e!)_llwoHRO6+I->`0ZaBE_g>;AHx-)q_qsk)>0So~py4?JNzhZ{|v`+Jke
z6rQ))^7EaQb$?{(XH2+rTP^)O;iqh;*`i&`RoYSe&Kt8Z?}2n5wft?$_VTs1)2jSC
zf62l<N}u<P<=?1q$Mv?;tNh)+Vqu;ing318&tonv@3Hg=<v;fimi~mwYy7tItDjeW
z))MAev}<{<g-Z%|f5djox;S^B%hD@KfAxUv^eVhW7qUg+XCAek5rxmKvE%nEeDYJa
zb4cNzwEv8c>m6K_zhw2Pcvt#G+dHP{lS^#pgu*S~xAYTgziSp+!fB<y^hMi&|AH55
zkCGZVYhEr^znK>tRDU+O*6_?58ZF{#U5AogJly2wI-hoR>2}u~=B>DeSGa|*iI>?y
ztWQrK_c9JBSHi;`<A(Tnf$;(UpYaCcTE-!a8yI&mzF-{8ID&D8jW43(2lJbz;SKJp
z;0l`E>aGoUb#-qZ80s4v^XK`N;QwX8f7)emcu)UO)wwGj$6X#C*R!a%YeSFaZ;9l;
z<%ZGf*0CFI;eJdW#}Dtz@%W+s+`<uD!WaA*S5xZ>@K~w@7z7(!!ZBQzD?ZgZ_=j(J
zg>%!BD7V&M@ONps<-TDaxxKC6HpSk;Cw##l{2Euz!8N?X9UQ}Dx#DSc4t;>*IREf#
R+()^!{>p}1?r+Dr{r@(`Uv2;Z


From c51de3149403692c5793c10323978b2bd3744052 Mon Sep 17 00:00:00 2001
From: Jimbles <j.avery.10@ucl.ac.uk>
Date: Tue, 29 Jun 2021 11:09:30 +0100
Subject: [PATCH 2/3] Added tests and data non windows specific

---
 tests/Test_ScouseTom.m          | 163 +++++++++++++++++++++++++
 tests/data/readme.md            |   2 +
 tests/filters/check_acc.m       | 121 ++++++++++++++++++
 tests/filters/test_accuracy.m   |   0
 tests/filters/test_decimate.m   |  55 +++++++++
 tests/filters/test_detrend.m    | 106 ++++++++++++++++
 tests/filters/test_filters.m    | 209 ++++++++++++++++++++++++++++++++
 tests/filters/test_singlecase.m |  90 ++++++++++++++
 tests/readme.md                 |   6 +
 9 files changed, 752 insertions(+)
 create mode 100644 tests/Test_ScouseTom.m
 create mode 100644 tests/data/readme.md
 create mode 100644 tests/filters/check_acc.m
 create mode 100644 tests/filters/test_accuracy.m
 create mode 100644 tests/filters/test_decimate.m
 create mode 100644 tests/filters/test_detrend.m
 create mode 100644 tests/filters/test_filters.m
 create mode 100644 tests/filters/test_singlecase.m
 create mode 100644 tests/readme.md

diff --git a/tests/Test_ScouseTom.m b/tests/Test_ScouseTom.m
new file mode 100644
index 0000000..ef17743
--- /dev/null
+++ b/tests/Test_ScouseTom.m
@@ -0,0 +1,163 @@
+% test code is functioning on real examples. Requires files from https://doi.org/10.5281/zenodo.4783547
+
+datadir='I:/Load_data_testing/';
+%datadir='data/';
+%% load bdfs
+
+somethingwentwrong=0;
+
+
+try
+    BVs=ScouseTom_Load([datadir 'Baseline1.bdf']);
+catch err
+    fprintf(2,'error on initial test\n');
+    somethingwentwrong=1;
+    fprintf(2, '%s\n', getReport(err, 'extended'));
+end
+
+disp('=================================================');
+disp('=================================================');
+disp('=================================================');
+
+try
+    BVs=ScouseTom_Load([datadir 'Baseline1_missingstart.bdf']);
+catch err
+    fprintf(2,'error on missing start test\n');
+    somethingwentwrong=1;
+    fprintf(2, '%s\n', getReport(err, 'extended'));
+end
+
+disp('=================================================');
+disp('=================================================');
+disp('=================================================');
+
+try
+    BVs=ScouseTom_Load([datadir 'Baseline1_missingend.bdf']);
+catch err
+    fprintf(2,'error on missing end test\n');
+    somethingwentwrong=1;
+    fprintf(2, '%s\n', getReport(err, 'extended'));
+end
+
+disp('=================================================');
+disp('=================================================');
+disp('=================================================');
+
+try
+    BVs=ScouseTom_Load([datadir 'Baseline1_missingboth.bdf']);
+catch err
+    fprintf(2,'error on missing both test\n');
+    somethingwentwrong=1;
+    fprintf(2, '%s\n', getReport(err, 'extended'));
+end
+
+disp('=================================================');
+disp('=================================================');
+disp('=================================================');
+%% multifreq
+
+try
+    BVs=ScouseTom_Load([datadir 'MultiFreq.bdf']);
+catch err
+    fprintf(2,'error on Multifreq test\n');
+    somethingwentwrong=1;
+    fprintf(2, '%s\n', getReport(err, 'extended'));
+end
+
+disp('=================================================');
+disp('=================================================');
+disp('=================================================');
+
+try
+    BVs=ScouseTom_Load([datadir 'MultiFreq_missingstart.bdf']);
+catch err
+    fprintf(2,'error on Multifreq missing start test\n');
+    somethingwentwrong=1;
+    fprintf(2, '%s\n', getReport(err, 'extended'));
+end
+
+disp('=================================================');
+disp('=================================================');
+disp('=================================================');
+
+try
+    BVs=ScouseTom_Load([datadir 'MultiFreq_missingend.bdf']);
+catch err
+    fprintf(2,'error on Multifreq missing end test\n');
+    somethingwentwrong=1;
+    fprintf(2, '%s\n', getReport(err, 'extended'));
+end
+
+disp('=================================================');
+disp('=================================================');
+disp('=================================================');
+
+try
+    BVs=ScouseTom_Load([datadir 'MultiFreq_missingboth.bdf']);
+catch err
+    fprintf(2,'error on Multifreq missing both test\n');
+    somethingwentwrong=1;
+end
+
+disp('=================================================');
+disp('=================================================');
+disp('=================================================');
+
+%% Z check
+try
+    Z=ScouseTom_Load([datadir 'firstinnir_zcheck_1.bdf']);
+catch err
+    fprintf(2,'error on Z load \n');
+    somethingwentwrong=1;
+    fprintf(2, '%s\n', getReport(err, 'extended'));
+end
+
+
+%% AC check
+
+try
+    BVs=ScouseTom_Load([datadir 'rptest2.eeg']);
+catch err
+    fprintf(2,'error on Actichamp test 1\n');
+    somethingwentwrong=1;
+end
+
+disp('=================================================');
+disp('=================================================');
+disp('=================================================');
+
+try
+    BVs=ScouseTom_Load([datadir 'rptest3.eeg']);
+catch err
+    fprintf(2,'error on Actichamp test 2\n');
+    somethingwentwrong=1;
+end
+
+disp('=================================================');
+disp('=================================================');
+disp('=================================================');
+
+
+%% Big file
+
+try
+    BVs=ScouseTom_Load([datadir 'RealThing_NIRFACE_2.bdf']);
+catch err
+    fprintf(2,'error on Multifreq big file test\n');
+    somethingwentwrong=1;
+end
+
+disp('=================================================');
+disp('=================================================');
+disp('=================================================');
+
+
+
+%% did it work?
+
+if somethingwentwrong
+    fprintf(2,'Oh No! SOMETHING WENT WRONG!\n');
+else
+    fprintf('YAY! IT WORKED :)\n');
+end
+
diff --git a/tests/data/readme.md b/tests/data/readme.md
new file mode 100644
index 0000000..e446dff
--- /dev/null
+++ b/tests/data/readme.md
@@ -0,0 +1,2 @@
+# Test dataset
+Put files from [![DOI](https://zenodo.org/badge/DOI/10.5281/zenodo.4783527.svg)](https://doi.org/10.5281/zenodo.4783527) in this directory to run tests
diff --git a/tests/filters/check_acc.m b/tests/filters/check_acc.m
new file mode 100644
index 0000000..fb08a3b
--- /dev/null
+++ b/tests/filters/check_acc.m
@@ -0,0 +1,121 @@
+function [ Amp_error, Phase_error,Vsig,Vsigdemod,Filt,trim_demod] = check_acc( Fc,InjTime,Amp_Inj,Amp_Meas,InjPhase,MeasPhaseDiff,DCoffset,DCoffsetInj, Padsec, Fs,decimate_factor)
+    %UNTITLED Summary of this function goes here
+    %   Detailed explanation goes here
+    
+    %%
+    inj = [ 2 4];
+    
+    if exist('DCoffset','var') == 0 || isempty(DCoffset)
+        DCoffset=0;
+    end
+    
+    if exist('DCoffsetInj','var') == 0 || isempty(DCoffsetInj)
+        DCoffsetInj=0;
+    end
+    
+    
+    BW=100;
+    chn=5;
+    
+    if exist('Padsec','var') == 0 || isempty(Padsec)
+        Padsec=4;
+    end
+    
+    if exist('Fs','var') == 0 || isempty(Fs)
+        Fs=16384;
+    end
+    
+    if exist('decimate_factor','var') == 0 || isempty(decimate_factor)
+        decimate_factor=0;
+    end
+    
+    
+    %% Create ideal values, and voltages
+    
+    AmpActual= repmat(Amp_Meas,chn,1);
+    AmpActual(inj)=Amp_Inj;
+    
+    a=MeasPhaseDiff;
+    
+    %force normal valus of deg
+    a = mod(a,360); % [0 2pi)
+    
+    % shift
+    jj = a > 180;
+    a(jj) = a(jj) - 360;
+    jj = a < 0 - 180;
+    a(jj) = a(jj) + 360;
+    
+    MeasPhaseDiff_corr = a;
+    
+    MeasPhase = InjPhase + MeasPhaseDiff_corr; %deg
+    PhaseActual= repmat(MeasPhaseDiff_corr,chn,1);
+    PhaseActual(inj)=0;
+    
+    Totaltime=ceil(InjTime + 2*Padsec);
+    t = 0:1/Fs:Totaltime-1/Fs;
+    %make sin wave
+    
+    v_m= Amp_Meas*sin(2*pi*Fc*t+(pi*MeasPhase/180))+DCoffset;
+    
+    % change amplitude
+    v_i=Amp_Inj*sin(2*pi*Fc*t+(pi*InjPhase/180))+DCoffsetInj;
+    
+    
+    %%
+    V=repmat(v_m,chn,1);
+    
+    V(inj(1),:) = v_i;
+    V(inj(2),:) = v_i;
+    
+    V=V';
+    
+    %% decimate now if we want to
+    
+    if decimate_factor
+        for iChn = 1 :size(V,2)
+            Vtmp(:,iChn) = decimate(V(:,iChn),decimate_factor);
+        end
+        Fs=Fs/decimate_factor;
+        V=Vtmp;
+    end
+    
+    %%
+    
+    
+    %pad with a second of data either side, so the hilbert is more realistic
+    datastart = round(Padsec*Fs);
+    dataend = round((Padsec+InjTime)*Fs);
+    
+    V(1:datastart,:)=V(1:datastart,:)*0.1;
+    V(dataend:end,:)=V(dataend:end,:)*0.1;
+    
+    InjectionWindows =[datastart dataend];
+    
+    
+    %%
+    %find the corresponding filter settings
+    [trim_demod,Filt,Fc_found]=ScouseTom_data_GetFilterTrim(V(datastart:dataend,inj(1)),Fs,[],[],[]);
+    
+    Vsig = V(datastart:dataend,inj(1));
+    
+    [ Vdata_demod,Pdata_demod ] = ScouseTom_data_DemodHilbert( V,Filt);
+    Vsigdemod = Vdata_demod(datastart:dataend,inj(1));
+    
+    %%
+    [Vmag,PhaseRaw,VmagSTD,PhaseRawSTD]=ScouseTom_data_getBV(Vdata_demod,Pdata_demod,trim_demod,InjectionWindows);
+    % [Vmag,PhaseRaw,VmagSTD,PhaseRawSTD]=ScouseTom_data_getBVChunk(Datafilt,trim_demod,InjectionWindows);
+    
+    
+    [Phase]=ScouseTom_data_PhaseEst(PhaseRaw,inj,1);
+    
+    Phase_deg=(180/pi) * (Phase);
+    
+    Phase_error = PhaseActual - Phase_deg';
+    
+    Amp_error = AmpActual - Vmag';
+    
+    fprintf('Amp error : %.6f, Phase error : %.6f\n',mean(Amp_error),mean(Phase_error));
+    end
+    
+    
\ No newline at end of file
diff --git a/tests/filters/test_accuracy.m b/tests/filters/test_accuracy.m
new file mode 100644
index 0000000..e69de29
diff --git a/tests/filters/test_decimate.m b/tests/filters/test_decimate.m
new file mode 100644
index 0000000..b3e8a05
--- /dev/null
+++ b/tests/filters/test_decimate.m
@@ -0,0 +1,55 @@
+Fcur = 1000;
+FreqNum = size(Fcur,2);
+
+% Cycles = 2*6000;
+% T=(1./Fcur); %Period in s
+% InjTime=(T.*Cycles);
+
+InjTime=2;
+
+
+Amp_Inj = 500;
+Amp_Meas = 150;
+InjPhase=0;
+MeasPhaseDiff=0;
+DCoffset = 0;
+DCoffsetinj = 0;
+
+Fs=100000;
+
+
+[Amp_error1, Phase_error1,V1,Vd1,Filt1,tr1] = check_acc( Fcur,InjTime,Amp_Inj,Amp_Meas,InjPhase,MeasPhaseDiff,DCoffset,DCoffsetinj,[],Fs);
+%%
+dec_factor=2;
+
+Fs2=Fs/dec_factor;
+
+[Amp_error2, Phase_error2,V2,Vd2,Filt2,tr2] = check_acc( Fcur,InjTime,Amp_Inj,Amp_Meas,InjPhase,MeasPhaseDiff,DCoffset,DCoffsetinj,[],Fs2);
+%%
+[Amp_error3, Phase_error3,V3,Vd3,Filt3,tr3] = check_acc( Fcur,InjTime,Amp_Inj,Amp_Meas,InjPhase,MeasPhaseDiff,DCoffset,DCoffsetinj,[],Fs,dec_factor);
+
+
+
+%%
+
+t1 = (0:1:length(Vd1)-1)/Fs;
+t2 = (0:1:length(Vd2)-1)/Fs2;
+t3 = (0:1:length(Vd3)-1)/Fs2;
+
+figure;
+hold on
+plot(t1,Vd1);
+plot(t2,Vd2);
+plot(t3,Vd3);
+hold off
+legend('Original','Ideal downsample','Decimated');
+
+figure;
+hold on
+plot(t1,V1);
+plot(t2,V2);
+plot(t3,V3);
+hold off
+legend('Original','Ideal downsample','Decimated');
+xlim([1, 1.01])
+
diff --git a/tests/filters/test_detrend.m b/tests/filters/test_detrend.m
new file mode 100644
index 0000000..9e561fb
--- /dev/null
+++ b/tests/filters/test_detrend.m
@@ -0,0 +1,106 @@
+
+Fc=20;
+
+Fcur = 15;
+FreqNum = size(Fcur,2);
+
+
+NumInj= 5;
+
+Cycles = 128;
+T=(1./Fcur); %Period in s
+InjTime=(T.*Cycles)*NumInj;
+InjSamples=(T.*Cycles)*Fs;
+
+% InjTime=10;
+
+Amp_Inj = 500;
+Amp_Meas = 150;
+InjPhase=0;
+MeasPhaseDiff=-30;
+
+inj = [ 2 4];
+
+DCoffset=10+10*(1:NumInj);
+DCoffsetInj=100 + 50*(1:NumInj);
+
+
+Fs=16384;
+BW=100;
+chn=5;
+
+Padsec=4;
+
+%% Create ideal values, and voltages
+
+AmpActual= repmat(Amp_Meas,chn,1);
+AmpActual(inj)=Amp_Inj;
+
+a=MeasPhaseDiff;
+
+%force normal valus of deg
+a = mod(a,360); % [0 2pi)
+
+% shift
+jj = a > 180;
+a(jj) = a(jj) - 360;
+jj = a < 0 - 180;
+a(jj) = a(jj) + 360;
+
+MeasPhaseDiff_corr = a;
+
+MeasPhase = InjPhase + MeasPhaseDiff_corr; %deg
+PhaseActual= repmat(MeasPhaseDiff_corr,chn,1);
+PhaseActual(inj)=0;
+
+Totaltime=ceil(InjTime + 2*Padsec);
+t = 0:1/Fs:Totaltime-1/Fs;
+%make sin wave
+
+v_m= Amp_Meas*sin(2*pi*Fc*t+(pi*MeasPhase/180));
+
+% change amplitude
+v_i=Amp_Inj*sin(2*pi*Fc*t+(pi*InjPhase/180));
+
+%%
+V=repmat(v_m,chn,1);
+
+V(inj(1),:) = v_i;
+V(inj(2),:) = v_i;
+
+V=V';
+
+%pad with a second of data either side, so the hilbert is more realistic
+datastart = round(Padsec*Fs);
+dataend = round((Padsec+InjTime)*Fs);
+
+V(1:datastart,:)=V(1:datastart,:)*0.1;
+V(dataend:end,:)=V(dataend:end,:)*0.1;
+
+InjectionWindows(:,1) = (datastart + ceil(InjSamples * (0:NumInj-1)))';
+InjectionWindows(:,2) = [InjectionWindows(2:end,1) ; dataend];
+
+
+DCoffsetVec = ones(chn,1)*DCoffset;
+DCoffsetVec(inj(1),:) = DCoffsetInj;
+DCoffsetVec(inj(2),:) = DCoffsetInj;
+
+for iInj = 1:NumInj
+
+V(InjectionWindows(iInj,1):InjectionWindows(iInj,2),:)=bsxfun(@plus,V(InjectionWindows(iInj,1):InjectionWindows(iInj,2),:),DCoffsetVec(:,iInj)');
+
+end
+
+
+%% 
+
+Voff =mean(V(datastart:dataend,:));
+
+plot(bsxfun(@minus,V,Voff))
+
+
+
+
+
+
+
diff --git a/tests/filters/test_filters.m b/tests/filters/test_filters.m
new file mode 100644
index 0000000..5a6b8c7
--- /dev/null
+++ b/tests/filters/test_filters.m
@@ -0,0 +1,209 @@
+%%
+
+% Fs=16384;
+Fs=100000;
+
+Fc_all = 5:100:550;
+FreqNum = length(Fc_all);
+
+BWtarget =100;
+
+Fstopomin=0.5;
+Fcomin = 1.5;
+
+
+for iFreq = 1: FreqNum
+    
+    
+    Fc=Fc_all(iFreq);
+    
+    disp(Fc);
+    
+    if Fc == 25;
+        disp('stop');
+    end
+    
+    %% old filter
+    
+    BWold = BWtarget;
+    Fco(2) = Fc+(BWold/2);
+    Fco(1) = max([(Fc -(BWold/2)), Fstopomin]);
+    [B,A] = butter(3,Fco./(Fs/2));
+    
+    stold(iFreq)= isstable(B,A);
+    imp_lenold(iFreq)=impzlength(B,A,0.001);
+    imp_len_msold(iFreq) = (imp_lenold(iFreq) / Fs)*(1000);
+    len_cyc = (1/Fc)*1000;
+    imp_len_cycold(iFreq)=imp_len_msold(iFreq)/len_cyc;
+    
+    [Hold,W] =freqz(B,A,Fco(1):0.2:Fco(2),Fs);
+    fc_idx = find (W > Fc,1);
+    Gainold(iFreq) =abs(Hold(fc_idx));
+    
+    
+    %% New Bandpass
+    
+    BWbp = BWtarget;
+    
+    StopBandDiffbp = 350;
+    
+    Fpass1bp = Fc-BWbp/2;
+    
+    Astop1 = 60;
+    Astop2= 60;
+    Apass  = 0.5;
+    
+    
+    if ~(Fpass1bp > 0)
+        Fpass1bp = max([Fc/2 Fcomin]);
+        Astop1 = 15;
+    end
+    
+    Fstop1bp = Fpass1bp -StopBandDiffbp;
+    Fstop1bp = max([Fstop1bp Fstopomin]);
+    
+    Fpass2bp = Fc+BWbp/2;
+    Fstop2bp = Fpass2bp + StopBandDiffbp;
+    
+    
+    
+    hbp = fdesign.bandpass('fst1,fp1,fp2,fst2,ast1,ap,ast2', Fstop1bp, Fpass1bp, ...
+        Fpass2bp, Fstop2bp, Astop1, Apass, Astop2, Fs);
+    
+    dbp = design(hbp, 'butter', 'MatchExactly', 'passband', 'SOSScaleNorm', 'Linf');
+    
+    %     dbp = designfilt('bandpassiir', ...       % Response type
+    %         'FilterOrder', 20, ...
+    %        'HalfPowerFrequency1',Fpass1bp, ...    % Frequency constraints
+    %        'HalfPowerFrequency2',Fpass2bp, ...
+    %        'SampleRate',Fs)   ;            % Sample rate
+    %
+    
+    
+    
+    stbp(iFreq)= isstable(dbp);
+    imp_lenbp(iFreq)=impzlength(dbp,0.001);
+    imp_len_msbp(iFreq) = (imp_lenbp(iFreq) / Fs)*(1000);
+    len_cyc = (1/Fc)*1000;
+    imp_len_cycbp(iFreq)=imp_len_msbp(iFreq)/len_cyc;
+    
+    [Hbp,W] =freqz(dbp,Fstop1bp:0.2:Fstop2bp,Fs);
+    fc_idx = find (W > Fc,1);
+    Gainbp(iFreq) =abs(Hbp(fc_idx));
+    
+    
+    %% New Cascaded filter
+    
+    BWc = BWtarget;
+    
+    StopBandDiffc = 350;
+    
+    Fpass1c = Fc-BWc/2;
+    
+    Astop1 = 60;
+    Astop2= 60;
+    Apass  = 0.5;
+    
+    if ~(Fpass1c > 0)
+        Fpass1c = max([Fc/2 Fcomin]);
+        Astop1 =30;
+    end
+    
+    Fstop1c = Fpass1c -StopBandDiffc;
+    Fstop1c = max([Fstop1c Fstopomin]);
+    
+    Fpass2c = Fc+BWc/2;
+    Fstop2c = Fpass2c + StopBandDiffc;
+    
+    
+    
+    %     hchp = fdesign.highpass('fst,fp,ast,ap', Fstop1c, Fpass1c, Astop1, Apass, Fs);
+    %
+    %     dchp = design(hchp, 'butter', 'MatchExactly', 'passband', 'SOSScaleNorm', 'Linf');
+    
+    %     hclp = fdesign.lowpass('fp,fst,ap,ast', Fpass2c, Fstop2c, Apass, Astop2, Fs);
+    %
+    %     dclp = design(hclp, 'butter', 'MatchExactly', 'passband', 'SOSScaleNorm', 'Linf');
+    
+    dclp = designfilt('lowpassiir', ...       % Response type
+        'FilterOrder', 10, ...
+        'HalfPowerFrequency',Fpass2c, ...    % Frequency constraints
+        'SampleRate',Fs)   ;            % Sample rate
+    
+    
+    dchp = designfilt('highpassiir', ...       % Response type
+        'FilterOrder', 10, ...
+        'HalfPowerFrequency',Fpass1c, ...    % Frequency constraints
+        'SampleRate',Fs)   ;            % Sample rate
+    
+    
+    stc(iFreq)= all([isstable(dchp) isstable(dclp)]);
+    imp_lenclp(iFreq)=impzlength(dclp,0.001) ;
+    imp_lenchp(iFreq) = impzlength(dchp,0.001);
+    imp_len_msclp(iFreq) = (imp_lenclp(iFreq) / Fs)*(1000);
+    imp_len_mschp(iFreq) = (imp_lenchp(iFreq) / Fs)*(1000);
+    len_cyc = (1/Fc)*1000;
+    imp_len_cycclp(iFreq)=imp_len_msclp(iFreq)/len_cyc;
+    imp_len_cycchp(iFreq)=imp_len_mschp(iFreq)/len_cyc;
+    [Hclp,W] =freqz(dclp,Fstop1c:0.2:Fstop2c,Fs);
+    [Hchp,W] =freqz(dchp,Fstop1c:0.2:Fstop2c,Fs);
+    fc_idx = find (W > Fc,1);
+    Gainclp(iFreq) = abs(Hclp(fc_idx)) ;
+    
+    Gainchp(iFreq) =  abs(Hchp(fc_idx));
+    
+    %% Freq Info
+    
+end
+
+
+%%
+
+figure;
+hold on
+plot(Fc_all,10*log10(Gainold));
+plot(Fc_all,10*log10(Gainbp));
+plot(Fc_all,10*log10(Gainclp));
+plot(Fc_all,10*log10(Gainchp));
+hold off
+legend('Old','Bandp','Casclp','Caschp')
+% legend('Bandp','Casclp','Caschp')
+title('gain')
+xlabel('Freq')
+ylabel('Gain db');
+
+
+figure;
+hold on
+plot(Fc_all,imp_len_msold);
+plot(Fc_all,imp_len_msbp);
+plot(Fc_all,imp_len_msclp);
+plot(Fc_all,imp_len_mschp);
+
+hold off
+legend('Old','Bandp','Casclp','Caschp')
+% legend('Bandp','Casclp','Caschp')
+title('timedelay ms')
+ylabel('Settling time ms')
+xlabel('Freq')
+
+
+
+figure;
+hold on
+plot(Fc_all,imp_len_cycold);
+plot(Fc_all,imp_len_cycbp);
+plot(Fc_all,imp_len_cycclp);
+plot(Fc_all,imp_len_cycchp);
+
+hold off
+legend('Old','Bandp','Casclp','Caschp')
+% legend('Bandp','Casclp','Caschp')
+title('timedelay cycles')
+ylabel('Settling time cycles')
+xlabel('Freq')
+
+
+
+
+
diff --git a/tests/filters/test_singlecase.m b/tests/filters/test_singlecase.m
new file mode 100644
index 0000000..1c8f25a
--- /dev/null
+++ b/tests/filters/test_singlecase.m
@@ -0,0 +1,90 @@
+Fcur = 5;
+FreqNum = size(Fcur,2);
+
+Cycles = 64;
+T=(1./Fcur); %Period in s
+InjTime=(T.*Cycles);
+
+% InjTime=10;
+
+
+Amp_Inj = 500;
+Amp_Meas = 150;
+InjPhase=0;
+MeasPhaseDiff=-30;
+DCoffset = 0;
+DCoffsetinj = 0;
+
+% Fs=100000;
+Fs=16384;
+
+[Amp_error1, Phase_error1,V1,Vd1,Filt1,tr1] = check_acc( Fcur,InjTime,Amp_Inj,Amp_Meas,InjPhase,MeasPhaseDiff,DCoffset,DCoffsetinj,[],Fs);
+
+trim_demod=200;
+
+
+t1 = (0:1:length(V1)-1)/Fs;
+
+
+%%
+
+Fcur = 5;
+FreqNum = size(Fcur,2);
+
+Cycles = 32;
+T=(1./Fcur); %Period in s
+InjTime=(T.*Cycles);
+
+% InjTime=10;
+
+
+Amp_Inj = 500;
+Amp_Meas = 150;
+InjPhase=0;
+MeasPhaseDiff=-30;
+DCoffset = 100;
+DCoffsetinj = 500;
+
+% Fs2=100000;
+Fs2=16384;
+
+
+[Amp_error2, Phase_error2,V2,Vd2,Filt2,tr2] = check_acc( Fcur,InjTime,Amp_Inj,Amp_Meas,InjPhase,MeasPhaseDiff,DCoffset,DCoffsetinj,[],Fs2);
+
+
+t2 = (0:1:length(V2)-1)/Fs2;
+
+%
+
+%%
+
+fvtool(Filt1,Filt2)
+%%
+[H1,F1]=freqz(Filt1,Fcur-500:Fcur+500,Fs);
+[H2,F2]=freqz(Filt2,Fcur-500:Fcur+500,Fs2);
+
+figure;
+hold on
+plot(F1,10*log10(abs(H1)));
+plot(F2,10*log10(abs(H2)));
+hold off
+
+%%
+
+figure
+hold on
+plot(t1,V1)
+plot(t2,V2)
+% plot(V8)
+hold off
+
+
+figure
+hold on
+% plot(Vd1(tr1:end-tr1))
+plot(t1,Vd1)
+plot(t2,Vd2)
+% plot(Vd2(tr2:end-tr2))
+% plot(Vd8)
+hold off
+% ylim(Amp_Inj + [-1 +1])
diff --git a/tests/readme.md b/tests/readme.md
new file mode 100644
index 0000000..8f5515a
--- /dev/null
+++ b/tests/readme.md
@@ -0,0 +1,6 @@
+# Functional Tests
+
+`Test_ScouseTom.m` tests require the files from [![DOI](https://zenodo.org/badge/DOI/10.5281/zenodo.4783547.svg)](https://doi.org/10.5281/zenodo.4783547) to be placed in the data directory as the EEG system files get very big very quickly. These tests only check that the code is functional and are by no means unit tests.
+
+### Filter tests
+Checks to characterise frequency dependence of filtering method. Historical curiosity more than anything else.

From 41eb965dd17013d3d5c0c4cdf6746ad09252f184 Mon Sep 17 00:00:00 2001
From: Jimbles <j.avery.10@ucl.ac.uk>
Date: Tue, 29 Jun 2021 11:09:30 +0100
Subject: [PATCH 3/3] Added tests and data non windows specific

---
 tests/Test_ScouseTom.m          | 163 +++++++
 tests/data/readme.md            |   2 +
 tests/filters/check_acc.m       | 121 +++++
 tests/filters/test_accuracy.m   | 811 ++++++++++++++++++++++++++++++++
 tests/filters/test_decimate.m   |  55 +++
 tests/filters/test_detrend.m    | 106 +++++
 tests/filters/test_filters.m    | 209 ++++++++
 tests/filters/test_singlecase.m |  90 ++++
 tests/readme.md                 |   6 +
 9 files changed, 1563 insertions(+)
 create mode 100644 tests/Test_ScouseTom.m
 create mode 100644 tests/data/readme.md
 create mode 100644 tests/filters/check_acc.m
 create mode 100644 tests/filters/test_accuracy.m
 create mode 100644 tests/filters/test_decimate.m
 create mode 100644 tests/filters/test_detrend.m
 create mode 100644 tests/filters/test_filters.m
 create mode 100644 tests/filters/test_singlecase.m
 create mode 100644 tests/readme.md

diff --git a/tests/Test_ScouseTom.m b/tests/Test_ScouseTom.m
new file mode 100644
index 0000000..ef17743
--- /dev/null
+++ b/tests/Test_ScouseTom.m
@@ -0,0 +1,163 @@
+% test code is functioning on real examples. Requires files from https://doi.org/10.5281/zenodo.4783547
+
+datadir='I:/Load_data_testing/';
+%datadir='data/';
+%% load bdfs
+
+somethingwentwrong=0;
+
+
+try
+    BVs=ScouseTom_Load([datadir 'Baseline1.bdf']);
+catch err
+    fprintf(2,'error on initial test\n');
+    somethingwentwrong=1;
+    fprintf(2, '%s\n', getReport(err, 'extended'));
+end
+
+disp('=================================================');
+disp('=================================================');
+disp('=================================================');
+
+try
+    BVs=ScouseTom_Load([datadir 'Baseline1_missingstart.bdf']);
+catch err
+    fprintf(2,'error on missing start test\n');
+    somethingwentwrong=1;
+    fprintf(2, '%s\n', getReport(err, 'extended'));
+end
+
+disp('=================================================');
+disp('=================================================');
+disp('=================================================');
+
+try
+    BVs=ScouseTom_Load([datadir 'Baseline1_missingend.bdf']);
+catch err
+    fprintf(2,'error on missing end test\n');
+    somethingwentwrong=1;
+    fprintf(2, '%s\n', getReport(err, 'extended'));
+end
+
+disp('=================================================');
+disp('=================================================');
+disp('=================================================');
+
+try
+    BVs=ScouseTom_Load([datadir 'Baseline1_missingboth.bdf']);
+catch err
+    fprintf(2,'error on missing both test\n');
+    somethingwentwrong=1;
+    fprintf(2, '%s\n', getReport(err, 'extended'));
+end
+
+disp('=================================================');
+disp('=================================================');
+disp('=================================================');
+%% multifreq
+
+try
+    BVs=ScouseTom_Load([datadir 'MultiFreq.bdf']);
+catch err
+    fprintf(2,'error on Multifreq test\n');
+    somethingwentwrong=1;
+    fprintf(2, '%s\n', getReport(err, 'extended'));
+end
+
+disp('=================================================');
+disp('=================================================');
+disp('=================================================');
+
+try
+    BVs=ScouseTom_Load([datadir 'MultiFreq_missingstart.bdf']);
+catch err
+    fprintf(2,'error on Multifreq missing start test\n');
+    somethingwentwrong=1;
+    fprintf(2, '%s\n', getReport(err, 'extended'));
+end
+
+disp('=================================================');
+disp('=================================================');
+disp('=================================================');
+
+try
+    BVs=ScouseTom_Load([datadir 'MultiFreq_missingend.bdf']);
+catch err
+    fprintf(2,'error on Multifreq missing end test\n');
+    somethingwentwrong=1;
+    fprintf(2, '%s\n', getReport(err, 'extended'));
+end
+
+disp('=================================================');
+disp('=================================================');
+disp('=================================================');
+
+try
+    BVs=ScouseTom_Load([datadir 'MultiFreq_missingboth.bdf']);
+catch err
+    fprintf(2,'error on Multifreq missing both test\n');
+    somethingwentwrong=1;
+end
+
+disp('=================================================');
+disp('=================================================');
+disp('=================================================');
+
+%% Z check
+try
+    Z=ScouseTom_Load([datadir 'firstinnir_zcheck_1.bdf']);
+catch err
+    fprintf(2,'error on Z load \n');
+    somethingwentwrong=1;
+    fprintf(2, '%s\n', getReport(err, 'extended'));
+end
+
+
+%% AC check
+
+try
+    BVs=ScouseTom_Load([datadir 'rptest2.eeg']);
+catch err
+    fprintf(2,'error on Actichamp test 1\n');
+    somethingwentwrong=1;
+end
+
+disp('=================================================');
+disp('=================================================');
+disp('=================================================');
+
+try
+    BVs=ScouseTom_Load([datadir 'rptest3.eeg']);
+catch err
+    fprintf(2,'error on Actichamp test 2\n');
+    somethingwentwrong=1;
+end
+
+disp('=================================================');
+disp('=================================================');
+disp('=================================================');
+
+
+%% Big file
+
+try
+    BVs=ScouseTom_Load([datadir 'RealThing_NIRFACE_2.bdf']);
+catch err
+    fprintf(2,'error on Multifreq big file test\n');
+    somethingwentwrong=1;
+end
+
+disp('=================================================');
+disp('=================================================');
+disp('=================================================');
+
+
+
+%% did it work?
+
+if somethingwentwrong
+    fprintf(2,'Oh No! SOMETHING WENT WRONG!\n');
+else
+    fprintf('YAY! IT WORKED :)\n');
+end
+
diff --git a/tests/data/readme.md b/tests/data/readme.md
new file mode 100644
index 0000000..e446dff
--- /dev/null
+++ b/tests/data/readme.md
@@ -0,0 +1,2 @@
+# Test dataset
+Put files from [![DOI](https://zenodo.org/badge/DOI/10.5281/zenodo.4783527.svg)](https://doi.org/10.5281/zenodo.4783527) in this directory to run tests
diff --git a/tests/filters/check_acc.m b/tests/filters/check_acc.m
new file mode 100644
index 0000000..fb08a3b
--- /dev/null
+++ b/tests/filters/check_acc.m
@@ -0,0 +1,121 @@
+function [ Amp_error, Phase_error,Vsig,Vsigdemod,Filt,trim_demod] = check_acc( Fc,InjTime,Amp_Inj,Amp_Meas,InjPhase,MeasPhaseDiff,DCoffset,DCoffsetInj, Padsec, Fs,decimate_factor)
+    %UNTITLED Summary of this function goes here
+    %   Detailed explanation goes here
+    
+    %%
+    inj = [ 2 4];
+    
+    if exist('DCoffset','var') == 0 || isempty(DCoffset)
+        DCoffset=0;
+    end
+    
+    if exist('DCoffsetInj','var') == 0 || isempty(DCoffsetInj)
+        DCoffsetInj=0;
+    end
+    
+    
+    BW=100;
+    chn=5;
+    
+    if exist('Padsec','var') == 0 || isempty(Padsec)
+        Padsec=4;
+    end
+    
+    if exist('Fs','var') == 0 || isempty(Fs)
+        Fs=16384;
+    end
+    
+    if exist('decimate_factor','var') == 0 || isempty(decimate_factor)
+        decimate_factor=0;
+    end
+    
+    
+    %% Create ideal values, and voltages
+    
+    AmpActual= repmat(Amp_Meas,chn,1);
+    AmpActual(inj)=Amp_Inj;
+    
+    a=MeasPhaseDiff;
+    
+    %force normal valus of deg
+    a = mod(a,360); % [0 2pi)
+    
+    % shift
+    jj = a > 180;
+    a(jj) = a(jj) - 360;
+    jj = a < 0 - 180;
+    a(jj) = a(jj) + 360;
+    
+    MeasPhaseDiff_corr = a;
+    
+    MeasPhase = InjPhase + MeasPhaseDiff_corr; %deg
+    PhaseActual= repmat(MeasPhaseDiff_corr,chn,1);
+    PhaseActual(inj)=0;
+    
+    Totaltime=ceil(InjTime + 2*Padsec);
+    t = 0:1/Fs:Totaltime-1/Fs;
+    %make sin wave
+    
+    v_m= Amp_Meas*sin(2*pi*Fc*t+(pi*MeasPhase/180))+DCoffset;
+    
+    % change amplitude
+    v_i=Amp_Inj*sin(2*pi*Fc*t+(pi*InjPhase/180))+DCoffsetInj;
+    
+    
+    %%
+    V=repmat(v_m,chn,1);
+    
+    V(inj(1),:) = v_i;
+    V(inj(2),:) = v_i;
+    
+    V=V';
+    
+    %% decimate now if we want to
+    
+    if decimate_factor
+        for iChn = 1 :size(V,2)
+            Vtmp(:,iChn) = decimate(V(:,iChn),decimate_factor);
+        end
+        Fs=Fs/decimate_factor;
+        V=Vtmp;
+    end
+    
+    %%
+    
+    
+    %pad with a second of data either side, so the hilbert is more realistic
+    datastart = round(Padsec*Fs);
+    dataend = round((Padsec+InjTime)*Fs);
+    
+    V(1:datastart,:)=V(1:datastart,:)*0.1;
+    V(dataend:end,:)=V(dataend:end,:)*0.1;
+    
+    InjectionWindows =[datastart dataend];
+    
+    
+    %%
+    %find the corresponding filter settings
+    [trim_demod,Filt,Fc_found]=ScouseTom_data_GetFilterTrim(V(datastart:dataend,inj(1)),Fs,[],[],[]);
+    
+    Vsig = V(datastart:dataend,inj(1));
+    
+    [ Vdata_demod,Pdata_demod ] = ScouseTom_data_DemodHilbert( V,Filt);
+    Vsigdemod = Vdata_demod(datastart:dataend,inj(1));
+    
+    %%
+    [Vmag,PhaseRaw,VmagSTD,PhaseRawSTD]=ScouseTom_data_getBV(Vdata_demod,Pdata_demod,trim_demod,InjectionWindows);
+    % [Vmag,PhaseRaw,VmagSTD,PhaseRawSTD]=ScouseTom_data_getBVChunk(Datafilt,trim_demod,InjectionWindows);
+    
+    
+    [Phase]=ScouseTom_data_PhaseEst(PhaseRaw,inj,1);
+    
+    Phase_deg=(180/pi) * (Phase);
+    
+    Phase_error = PhaseActual - Phase_deg';
+    
+    Amp_error = AmpActual - Vmag';
+    
+    fprintf('Amp error : %.6f, Phase error : %.6f\n',mean(Amp_error),mean(Phase_error));
+    end
+    
+    
\ No newline at end of file
diff --git a/tests/filters/test_accuracy.m b/tests/filters/test_accuracy.m
new file mode 100644
index 0000000..6251ee1
--- /dev/null
+++ b/tests/filters/test_accuracy.m
@@ -0,0 +1,811 @@
+% script to test the accuracy of the processing code
+
+Failed = 0;
+Amp_tolerance = 5e-3; %how much error is ok for amp use about 10e-5 smaller than amp
+Phase_tolerance = 5e-3;
+
+plotflag = 1;
+
+%% test 1
+% fixed amplitude, time and freq. Sweep through phase
+
+disp('Test 1');
+
+Fc = 1800;
+
+InjTime=2;
+
+Amp_Inj = 500;
+Amp_Meas = 150;
+
+InjPhase = -95:20:95;
+InjPhaseNum = size(InjPhase,2);
+MeasPhaseDiff = -355:40:355;
+PhaseNum = size(MeasPhaseDiff,2);
+
+Amp_error_tot=zeros(InjPhaseNum,PhaseNum);
+Phase_error_tot=zeros(InjPhaseNum,PhaseNum);
+
+for iInjPhase = 1:InjPhaseNum
+    curInjPhase = InjPhase(iInjPhase);
+    for iPhase = 1:PhaseNum
+        
+        evalc('[Amp_error, Phase_error] = check_acc( Fc,InjTime,Amp_Inj,Amp_Meas,curInjPhase,MeasPhaseDiff(iPhase) );');
+        Amp_error_totT1(iInjPhase,iPhase) = mean(Amp_error);
+        Phase_error_totT1(iInjPhase,iPhase) = mean(Phase_error);
+        
+    end
+end
+%%
+if plotflag
+    figure
+    
+    subplot(2,1,1)
+    plot(MeasPhaseDiff,Amp_error_totT1)
+    ylabel('Error');
+    xlabel('Phase difference beetween inj and meas');
+    title('T1: Amp error for different starting phase');
+    
+    subplot(2,1,2)
+    plot(MeasPhaseDiff,Phase_error_totT1)
+    ylabel('Error');
+    xlabel('Phase difference');
+    title('T1: Phase error for different starting phase');
+    
+end
+
+Test1_AmpError=max(max(abs(Amp_error_totT1)));
+Test1_PhaseError=max(max(abs(Phase_error_totT1)));
+
+
+try
+    assert( Test1_AmpError < Amp_tolerance, 'Test1 Amp error failed');
+catch err
+    Failed = 1;
+    fprintf(2, '%s\n', getReport(err, 'extended'));
+end
+
+try
+    assert( Test1_PhaseError < Phase_tolerance, 'Test1 Phase error failed');
+catch err
+    Failed = 1;
+    fprintf(2, '%s\n', getReport(err, 'extended'));
+end
+
+%% test 2
+
+clear Amp_error_tot Phase_error_tot
+
+%changing amplitude
+disp('Test 2');
+
+InjTime=2;
+Fc = 1800;
+Amp_Inj = 500;
+Amp_Meas = 150;
+InjPhase=0;
+MeasPhaseDiff=-30;
+
+Amp_sc=0.5:0.5:10;
+AmpNum = size(Amp_sc,2);
+
+for iAmp = 1:AmpNum
+    
+    evalc('[Amp_error, Phase_error] = check_acc( Fc,InjTime,Amp_Inj*Amp_sc(iAmp),Amp_Meas*Amp_sc(iAmp),InjPhase,MeasPhaseDiff);');
+    Amp_error_totT2(iAmp) = mean(Amp_error);
+    Phase_error_totT2(iAmp) = mean(Phase_error);
+    
+end
+
+%%
+if plotflag
+    figure
+    
+    subplot(2,1,1)
+    plot(Amp_Meas*Amp_sc,Amp_error_totT2)
+    ylabel('Error');
+    xlabel('Amplitude');
+    title('T2: Amp error for different Amplitudes');
+    
+    subplot(2,1,2)
+    plot(Amp_Meas*Amp_sc,Phase_error_totT2)
+    ylabel('Error');
+    xlabel('Amplitude');
+    title('T2: Phase error for different Amplitudes');
+    
+end
+
+Test2_AmpError=max(max(abs(Amp_error_totT2)));
+Test2_PhaseError=max(max(abs(Phase_error_totT2)));
+
+
+try
+    assert( Test2_AmpError < Amp_tolerance, 'Test2 Amp error failed')
+catch err
+    Failed = 1;
+    fprintf(2, '%s\n', getReport(err, 'extended'));
+end
+
+try
+    assert( Test2_PhaseError < Phase_tolerance, 'Test2 Phase error failed')
+catch err
+    Failed = 1;
+    fprintf(2, '%s\n', getReport(err, 'extended'));
+end
+
+%% test 3
+
+clear Amp_error_tot Phase_error_tot
+
+%changing freq - large fixed inj time
+disp('Test 3  - large fixed injection time 10s');
+
+InjTime=10;
+
+Fc = [5:5:150 200:200:2000];
+FreqNum = size(Fc,2);
+
+Amp_Inj = 500;
+Amp_Meas = 150;
+InjPhase=0;
+MeasPhaseDiff=-30;
+
+for iFreq = 1:FreqNum
+    
+    evalc('[Amp_error, Phase_error] = check_acc( Fc(iFreq),InjTime,Amp_Inj,Amp_Meas,InjPhase,MeasPhaseDiff);');
+    Amp_error_totT3(iFreq) = mean(Amp_error);
+    Phase_error_totT3(iFreq) = mean(Phase_error);
+    
+end
+
+%%
+if plotflag
+    figure
+    
+    subplot(2,1,1)
+    plot(Fc,Amp_error_totT3)
+    ylabel('Error');
+    xlabel('Frequency');
+    title('T3: Amp error 10s injection');
+    
+    subplot(2,1,2)
+    plot(Fc,Phase_error_totT3)
+    ylabel('Error');
+    xlabel('Frequency');
+    title('T3: Phase error 10s injection');
+    
+end
+
+Test3_AmpError=max(max(abs(Amp_error_totT3)));
+Test3_PhaseError=max(max(abs(Phase_error_totT3)));
+
+try
+    assert( Test3_AmpError < Amp_tolerance, 'Test3 Amp error failed')
+catch err
+    Failed = 1;
+    fprintf(2, '%s\n', getReport(err, 'extended'));
+end
+
+try
+    assert( Test3_PhaseError < Phase_tolerance, 'Test3 Phase error failed')
+catch err
+    Failed = 1;
+    fprintf(2, '%s\n', getReport(err, 'extended'));
+end
+
+
+%% test 4
+
+clear Amp_error_tot Phase_error_tot
+
+%changing freq - large fixed inj time
+disp('Test 4 - small fixed injection time 0.5s');
+
+InjTime=.5;
+
+Fc = [5:5:150 200:200:2000];
+FreqNum = size(Fc,2);
+
+Amp_Inj = 500;
+Amp_Meas = 150;
+InjPhase=0;
+MeasPhaseDiff=-30;
+
+for iFreq = 1:FreqNum
+    
+    evalc('[Amp_error, Phase_error] = check_acc( Fc(iFreq),InjTime,Amp_Inj,Amp_Meas,InjPhase,MeasPhaseDiff);');
+    Amp_error_totT4(iFreq) = mean(Amp_error);
+    Phase_error_totT4(iFreq) = mean(Phase_error);
+    
+end
+
+%%
+if plotflag
+    figure
+    
+    subplot(2,1,1)
+    plot(Fc,Amp_error_totT4)
+    ylabel('Error');
+    xlabel('Frequency');
+    title('T4: Amp error 0.5s injection');
+    
+    subplot(2,1,2)
+    plot(Fc,Phase_error_totT4)
+    ylabel('Error');
+    xlabel('Frequency');
+    title('T4: Phase error 0.5s injection');
+    
+end
+
+Test4_AmpError=max(max(abs(Amp_error_totT4)));
+Test4_PhaseError=max(max(abs(Phase_error_totT4)));
+
+try
+    assert( Test4_AmpError < Amp_tolerance, 'Test4 Amp error failed')
+catch err
+    Failed = 1;
+    fprintf(2, '%s\n', getReport(err, 'extended'));
+end
+
+try
+    assert( Test4_PhaseError < Phase_tolerance, 'Test4 Phase error failed')
+catch err
+    Failed = 1;
+    fprintf(2, '%s\n', getReport(err, 'extended'));
+end
+
+%% test 5
+
+clear Amp_error_tot Phase_error_tot
+
+%changing freq - fixed cycles of 32
+disp('Test 5 - fixed cycles of 32');
+
+Fc = [5:5:150 200:200:2000];
+FreqNum = size(Fc,2);
+
+Cycles = 32;
+T=(1./Fc); %Period in s
+InjTime=(T.*Cycles);
+
+
+Amp_Inj = 500;
+Amp_Meas = 150;
+InjPhase=0;
+MeasPhaseDiff=-30;
+
+for iFreq = 1:FreqNum
+    
+    evalc('[Amp_error, Phase_error] = check_acc( Fc(iFreq),InjTime(iFreq),Amp_Inj,Amp_Meas,InjPhase,MeasPhaseDiff);');
+    Amp_error_totT5(iFreq) = mean(Amp_error);
+    Phase_error_totT5(iFreq) = mean(Phase_error);
+    
+end
+
+%%
+if plotflag
+    figure
+    
+    subplot(2,1,1)
+    plot(Fc,Amp_error_totT5)
+    ylabel('Error');
+    xlabel('Frequency');
+    title('T5: Amp error 32 cycles');
+    
+    subplot(2,1,2)
+    plot(Fc,Phase_error_totT5)
+    ylabel('Error');
+    xlabel('Frequency');
+    title('T5: Phase error 32 cycles');
+    
+end
+
+Test5_AmpError=max(max(abs(Amp_error_totT5)));
+Test5_PhaseError=max(max(abs(Phase_error_totT5)));
+
+try
+    assert( Test5_AmpError < Amp_tolerance, 'Test5 Amp error failed')
+catch err
+    Failed = 1;
+    fprintf(2, '%s\n', getReport(err, 'extended'));
+end
+
+try
+    assert( Test5_PhaseError < Phase_tolerance, 'Test5 Phase error failed')
+catch err
+    Failed = 1;
+    fprintf(2, '%s\n', getReport(err, 'extended'));
+end
+
+%% test 6
+
+clear Amp_error_tot Phase_error_tot
+
+%changing freq - fixed cycles of 3
+disp('Test 6 - fixed cycles of 3');
+
+Fc = [5:5:150 200:200:2000];
+FreqNum = size(Fc,2);
+
+Cycles = 3;
+T=(1./Fc); %Period in s
+InjTime=(T.*Cycles);
+
+
+Amp_Inj = 500;
+Amp_Meas = 150;
+InjPhase=0;
+MeasPhaseDiff=-30;
+
+for iFreq = 1:FreqNum
+    
+    evalc('[Amp_error, Phase_error] = check_acc( Fc(iFreq),InjTime(iFreq),Amp_Inj,Amp_Meas,InjPhase,MeasPhaseDiff);');
+    Amp_error_totT6(iFreq) = mean(Amp_error);
+    Phase_error_totT6(iFreq) = mean(Phase_error);
+    
+end
+
+%%
+if plotflag
+    figure
+    
+    subplot(2,1,1)
+    plot(Fc,Amp_error_totT6)
+    ylabel('Error');
+    xlabel('Frequency');
+    title('T6: Amp error 3 cycles');
+    
+    subplot(2,1,2)
+    plot(Fc,Phase_error_totT6)
+    ylabel('Error');
+    xlabel('Frequency');
+    title('T6: Phase error 3 cycles');
+    
+end
+
+Test6_AmpError=max(max(abs(Amp_error_totT6)));
+Test6_PhaseError=max(max(abs(Phase_error_totT6)));
+
+try
+    assert( Test6_AmpError < Amp_tolerance, 'Test6 Amp error failed')
+catch err
+    Failed = 1;
+    fprintf(2, '%s\n', getReport(err, 'extended'));
+end
+
+try
+    assert( Test6_PhaseError < Phase_tolerance, 'Test6 Phase error failed')
+catch err
+    Failed = 1;
+    fprintf(2, '%s\n', getReport(err, 'extended'));
+end
+
+
+%% test 7
+
+clear Amp_error_tot Phase_error_tot
+
+%changing freq - fixed cycles of 64
+disp('Test 7 - fixed cycles of 64');
+
+Fc = [5:5:150 200:200:2000];
+FreqNum = size(Fc,2);
+
+Cycles = 64;
+T=(1./Fc); %Period in s
+InjTime=(T.*Cycles);
+
+
+Amp_Inj = 500;
+Amp_Meas = 150;
+InjPhase=0;
+MeasPhaseDiff=-30;
+
+for iFreq = 1:FreqNum
+    
+    evalc('[Amp_error, Phase_error] = check_acc( Fc(iFreq),InjTime(iFreq),Amp_Inj,Amp_Meas,InjPhase,MeasPhaseDiff);');
+    Amp_error_totT7(iFreq) = mean(Amp_error);
+    Phase_error_totT7(iFreq) = mean(Phase_error);
+    
+end
+
+%%
+if plotflag
+    figure
+    
+    subplot(2,1,1)
+    plot(Fc,Amp_error_totT7)
+    ylabel('Error');
+    xlabel('Frequency');
+    title('T7: Amp error 64 cycles');
+    
+    subplot(2,1,2)
+    plot(Fc,Phase_error_totT7)
+    ylabel('Error');
+    xlabel('Frequency');
+    title('T7: Phase error 64 cycles');
+    
+end
+
+Test7_AmpError=max(max(abs(Amp_error_totT7)));
+Test7_PhaseError=max(max(abs(Phase_error_totT7)));
+
+try
+    assert( Test7_AmpError < Amp_tolerance, 'Test7 Amp error failed')
+catch err
+    Failed = 1;
+    fprintf(2, '%s\n', getReport(err, 'extended'));
+end
+
+try
+    assert( Test7_PhaseError < Phase_tolerance, 'Test7 Phase error failed')
+catch err
+    Failed = 1;
+    fprintf(2, '%s\n', getReport(err, 'extended'));
+end
+
+%% test 8
+
+clear Amp_error_tot Phase_error_tot
+
+%changing freq - fixed cycles of 128
+disp('Test 8 - fixed cycles of 128');
+
+Fc = [5:5:150 200:200:2000];
+FreqNum = size(Fc,2);
+
+Cycles = 128;
+T=(1./Fc); %Period in s
+InjTime=(T.*Cycles);
+
+
+Amp_Inj = 500;
+Amp_Meas = 150;
+InjPhase=0;
+MeasPhaseDiff=-30;
+
+for iFreq = 1:FreqNum
+    
+    evalc('[Amp_error, Phase_error] = check_acc( Fc(iFreq),InjTime(iFreq),Amp_Inj,Amp_Meas,InjPhase,MeasPhaseDiff);');
+    Amp_error_totT8(iFreq) = mean(Amp_error);
+    Phase_error_totT8(iFreq) = mean(Phase_error);
+    
+end
+
+%%
+if plotflag
+    figure
+    
+    subplot(2,1,1)
+    plot(Fc,Amp_error_totT8)
+    ylabel('Error');
+    xlabel('Frequency');
+    title('T8: Amp error 128 cycles');
+    
+    subplot(2,1,2)
+    plot(Fc,Phase_error_totT8)
+    ylabel('Error');
+    xlabel('Frequency');
+    title('T8: Phase error 128 cycles');
+    
+end
+
+Test8_AmpError=max(max(abs(Amp_error_totT8)));
+Test8_PhaseError=max(max(abs(Phase_error_totT8)));
+
+try
+    assert( Test8_AmpError < Amp_tolerance, 'Test8 Amp error failed')
+catch err
+    Failed = 1;
+    fprintf(2, '%s\n', getReport(err, 'extended'));
+end
+
+try
+    assert( Test8_PhaseError < Phase_tolerance, 'Test8 Phase error failed')
+catch err
+    Failed = 1;
+    fprintf(2, '%s\n', getReport(err, 'extended'));
+end
+
+
+%% test 9
+
+clear Amp_error_tot Phase_error_tot
+
+%changing freq - fixed cycles of 128
+disp('Test 9 - fixed cycles of 128 with dc offset');
+
+Fc = [5:5:150 200:200:2000];
+FreqNum = size(Fc,2);
+
+Cycles = 128;
+T=(1./Fc); %Period in s
+InjTime=(T.*Cycles);
+
+
+Amp_Inj = 500;
+Amp_Meas = 150;
+InjPhase=0;
+MeasPhaseDiff=-30;
+DCoffset = 100;
+DCoffsetinj = 400;
+
+for iFreq = 1:FreqNum
+    
+    evalc('[Amp_error, Phase_error] = check_acc( Fc(iFreq),InjTime(iFreq),Amp_Inj,Amp_Meas,InjPhase,MeasPhaseDiff,DCoffset,DCoffsetinj );');
+    Amp_error_totT9(iFreq) = mean(Amp_error);
+    Phase_error_totT9(iFreq) = mean(Phase_error);
+    
+end
+
+%%
+if plotflag
+    figure
+    
+    subplot(2,1,1)
+    plot(Fc,Amp_error_totT9)
+    ylabel('Error');
+    xlabel('Frequency');
+    title('T9: Amp error 128 cycles');
+    
+    subplot(2,1,2)
+    plot(Fc,Phase_error_totT9)
+    ylabel('Error');
+    xlabel('Frequency');
+    title('T9: Phase error 128 cycles');
+    
+end
+
+Test9_AmpError=max(max(abs(Amp_error_totT9)));
+Test9_PhaseError=max(max(abs(Phase_error_totT9)));
+
+try
+    assert( Test9_AmpError < Amp_tolerance, 'Test9 Amp error failed')
+catch err
+    Failed = 1;
+    fprintf(2, '%s\n', getReport(err, 'extended'));
+end
+
+try
+    assert( Test9_PhaseError < Phase_tolerance, 'Test9 Phase error failed')
+catch err
+    Failed = 1;
+    fprintf(2, '%s\n', getReport(err, 'extended'));
+end
+%% test 10
+
+clear Amp_error_tot Phase_error_tot
+
+%changing freq - fixed cycles of 32
+disp('Test 10 - fixed cycles of 32 with dc offset');
+
+Fc = [5:5:150 200:200:2000];
+FreqNum = size(Fc,2);
+
+Cycles = 32;
+T=(1./Fc); %Period in s
+InjTime=(T.*Cycles);
+
+
+Amp_Inj = 500;
+Amp_Meas = 150;
+InjPhase=0;
+MeasPhaseDiff=-30;
+DCoffset = 100;
+DCoffsetinj = 400;
+
+for iFreq = 1:FreqNum
+    
+    evalc('[Amp_error, Phase_error] = check_acc( Fc(iFreq),InjTime(iFreq),Amp_Inj,Amp_Meas,InjPhase,MeasPhaseDiff,DCoffset,DCoffsetinj );');
+    Amp_error_totT10(iFreq) = mean(Amp_error);
+    Phase_error_totT10(iFreq) = mean(Phase_error);
+    
+end
+
+%%
+if plotflag
+    figure
+    
+    subplot(2,1,1)
+    plot(Fc,Amp_error_totT10)
+    ylabel('Error');
+    xlabel('Frequency');
+    title('T10: Amp error 32 cycles');
+    
+    subplot(2,1,2)
+    plot(Fc,Phase_error_totT10)
+    ylabel('Error');
+    xlabel('Frequency');
+    title('T10: Phase error 32 cycles');
+    
+end
+
+Test10_AmpError=max(max(abs(Amp_error_totT10)));
+Test10_PhaseError=max(max(abs(Phase_error_totT10)));
+
+try
+    assert( Test10_AmpError < Amp_tolerance, 'Test10 Amp error failed')
+catch err
+    Failed = 1;
+    fprintf(2, '%s\n', getReport(err, 'extended'));
+end
+
+try
+    assert( Test10_PhaseError < Phase_tolerance, 'Test10 Phase error failed')
+catch err
+    Failed = 1;
+    fprintf(2, '%s\n', getReport(err, 'extended'));
+end
+
+%% test 11
+
+clear Amp_error_tot Phase_error_tot
+
+%changing freq - fixed cycles of 3
+disp('Test 11 - fixed cycles of 3 with dc offset');
+
+Fc = [5:5:150 200:200:2000];
+FreqNum = size(Fc,2);
+
+Cycles = 3;
+T=(1./Fc); %Period in s
+InjTime=(T.*Cycles);
+
+
+Amp_Inj = 500;
+Amp_Meas = 150;
+InjPhase=0;
+MeasPhaseDiff=-30;
+DCoffset = 100;
+DCoffsetinj = 400;
+
+for iFreq = 1:FreqNum
+    
+    evalc('[Amp_error, Phase_error] = check_acc( Fc(iFreq),InjTime(iFreq),Amp_Inj,Amp_Meas,InjPhase,MeasPhaseDiff,DCoffset,DCoffsetinj );');
+    Amp_error_totT11(iFreq) = mean(Amp_error);
+    Phase_error_totT11(iFreq) = mean(Phase_error);
+    
+end
+
+%%
+if plotflag
+    figure
+    
+    subplot(2,1,1)
+    plot(Fc,Amp_error_totT11)
+    ylabel('Error');
+    xlabel('Frequency');
+    title('T11: Amp error 3 cycles');
+    
+    subplot(2,1,2)
+    plot(Fc,Phase_error_totT11)
+    ylabel('Error');
+    xlabel('Frequency');
+    title('T11: Phase error 3 cycles');
+    
+end
+
+Test11_AmpError=max(max(abs(Amp_error_totT11)));
+Test11_PhaseError=max(max(abs(Phase_error_totT11)));
+
+try
+    assert( Test11_AmpError < Amp_tolerance, 'Test11 Amp error failed')
+catch err
+    Failed = 1;
+    fprintf(2, '%s\n', getReport(err, 'extended'));
+end
+
+try
+    assert( Test11_PhaseError < Phase_tolerance, 'Test11 Phase error failed')
+catch err
+    Failed = 1;
+    fprintf(2, '%s\n', getReport(err, 'extended'));
+end
+
+% SNR TEST
+
+
+%%
+
+figure
+hold on
+plot(Fc,Amp_error_totT5);
+plot(Fc,Amp_error_totT7);
+plot(Fc,Amp_error_totT8);
+hold off
+ylabel('Error');
+xlabel('Frequency');
+legend('32','64','128')
+
+
+
+
+%% Check a given ExpSetup
+
+load('S3b_MF1_log.mat','ExpSetup');
+
+
+
+clear Amp_error_tot Phase_error_tot
+
+%changing freq - fixed cycles of 3
+disp('Test 12 - ExpSetup');
+
+Fc = ExpSetup.Freq';
+FreqNum = size(Fc,2);
+
+
+InjTime=ExpSetup.MeasurementTime/1000';
+
+
+Amp_Inj = 9500;
+Amp_Meas = 1200 ;
+InjPhase=0;
+MeasPhaseDiff=-30;
+DCoffset = 100;
+DCoffsetinj = 400;
+
+for iFreq = 1:FreqNum
+    
+    %     evalc('[Amp_error, Phase_error] = check_acc( Fc(iFreq),InjTime(iFreq),Amp_Inj,Amp_Meas,InjPhase,MeasPhaseDiff,DCoffset,DCoffsetinj );');
+    [Amp_error, Phase_error] = check_acc( Fc(iFreq),InjTime(iFreq),Amp_Inj,Amp_Meas,InjPhase,MeasPhaseDiff,DCoffset,DCoffsetinj );
+    
+    
+    Amp_error_totT12(iFreq) = mean(Amp_error);
+    Phase_error_totT12(iFreq) = mean(Phase_error);
+    
+end
+
+%%
+
+if plotflag
+    figure
+    
+    subplot(2,1,1)
+    plot(Fc,100*(Amp_error_totT12/Amp_Inj))
+    ylabel('Error %');
+    xlabel('Frequency');
+    title('T12: Amp error 3 cycles');
+    
+    subplot(2,1,2)
+    
+    ylabel('Error');
+    xlabel('Frequency');
+    title('T12: Phase error 3 cycles');
+    
+end
+%%
+
+figure
+hold on
+plot(Fc/1000,100*(Amp_error_totT12/Amp_Inj))
+plot(Fc/1000,100*(Phase_error_totT12/MeasPhaseDiff))
+
+ylabel('Error %');
+xlabel('Frequency (kHz)');
+
+legend('Amplitude','Phase');
+
+
+
+%%
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
diff --git a/tests/filters/test_decimate.m b/tests/filters/test_decimate.m
new file mode 100644
index 0000000..b3e8a05
--- /dev/null
+++ b/tests/filters/test_decimate.m
@@ -0,0 +1,55 @@
+Fcur = 1000;
+FreqNum = size(Fcur,2);
+
+% Cycles = 2*6000;
+% T=(1./Fcur); %Period in s
+% InjTime=(T.*Cycles);
+
+InjTime=2;
+
+
+Amp_Inj = 500;
+Amp_Meas = 150;
+InjPhase=0;
+MeasPhaseDiff=0;
+DCoffset = 0;
+DCoffsetinj = 0;
+
+Fs=100000;
+
+
+[Amp_error1, Phase_error1,V1,Vd1,Filt1,tr1] = check_acc( Fcur,InjTime,Amp_Inj,Amp_Meas,InjPhase,MeasPhaseDiff,DCoffset,DCoffsetinj,[],Fs);
+%%
+dec_factor=2;
+
+Fs2=Fs/dec_factor;
+
+[Amp_error2, Phase_error2,V2,Vd2,Filt2,tr2] = check_acc( Fcur,InjTime,Amp_Inj,Amp_Meas,InjPhase,MeasPhaseDiff,DCoffset,DCoffsetinj,[],Fs2);
+%%
+[Amp_error3, Phase_error3,V3,Vd3,Filt3,tr3] = check_acc( Fcur,InjTime,Amp_Inj,Amp_Meas,InjPhase,MeasPhaseDiff,DCoffset,DCoffsetinj,[],Fs,dec_factor);
+
+
+
+%%
+
+t1 = (0:1:length(Vd1)-1)/Fs;
+t2 = (0:1:length(Vd2)-1)/Fs2;
+t3 = (0:1:length(Vd3)-1)/Fs2;
+
+figure;
+hold on
+plot(t1,Vd1);
+plot(t2,Vd2);
+plot(t3,Vd3);
+hold off
+legend('Original','Ideal downsample','Decimated');
+
+figure;
+hold on
+plot(t1,V1);
+plot(t2,V2);
+plot(t3,V3);
+hold off
+legend('Original','Ideal downsample','Decimated');
+xlim([1, 1.01])
+
diff --git a/tests/filters/test_detrend.m b/tests/filters/test_detrend.m
new file mode 100644
index 0000000..9e561fb
--- /dev/null
+++ b/tests/filters/test_detrend.m
@@ -0,0 +1,106 @@
+
+Fc=20;
+
+Fcur = 15;
+FreqNum = size(Fcur,2);
+
+
+NumInj= 5;
+
+Cycles = 128;
+T=(1./Fcur); %Period in s
+InjTime=(T.*Cycles)*NumInj;
+InjSamples=(T.*Cycles)*Fs;
+
+% InjTime=10;
+
+Amp_Inj = 500;
+Amp_Meas = 150;
+InjPhase=0;
+MeasPhaseDiff=-30;
+
+inj = [ 2 4];
+
+DCoffset=10+10*(1:NumInj);
+DCoffsetInj=100 + 50*(1:NumInj);
+
+
+Fs=16384;
+BW=100;
+chn=5;
+
+Padsec=4;
+
+%% Create ideal values, and voltages
+
+AmpActual= repmat(Amp_Meas,chn,1);
+AmpActual(inj)=Amp_Inj;
+
+a=MeasPhaseDiff;
+
+%force normal valus of deg
+a = mod(a,360); % [0 2pi)
+
+% shift
+jj = a > 180;
+a(jj) = a(jj) - 360;
+jj = a < 0 - 180;
+a(jj) = a(jj) + 360;
+
+MeasPhaseDiff_corr = a;
+
+MeasPhase = InjPhase + MeasPhaseDiff_corr; %deg
+PhaseActual= repmat(MeasPhaseDiff_corr,chn,1);
+PhaseActual(inj)=0;
+
+Totaltime=ceil(InjTime + 2*Padsec);
+t = 0:1/Fs:Totaltime-1/Fs;
+%make sin wave
+
+v_m= Amp_Meas*sin(2*pi*Fc*t+(pi*MeasPhase/180));
+
+% change amplitude
+v_i=Amp_Inj*sin(2*pi*Fc*t+(pi*InjPhase/180));
+
+%%
+V=repmat(v_m,chn,1);
+
+V(inj(1),:) = v_i;
+V(inj(2),:) = v_i;
+
+V=V';
+
+%pad with a second of data either side, so the hilbert is more realistic
+datastart = round(Padsec*Fs);
+dataend = round((Padsec+InjTime)*Fs);
+
+V(1:datastart,:)=V(1:datastart,:)*0.1;
+V(dataend:end,:)=V(dataend:end,:)*0.1;
+
+InjectionWindows(:,1) = (datastart + ceil(InjSamples * (0:NumInj-1)))';
+InjectionWindows(:,2) = [InjectionWindows(2:end,1) ; dataend];
+
+
+DCoffsetVec = ones(chn,1)*DCoffset;
+DCoffsetVec(inj(1),:) = DCoffsetInj;
+DCoffsetVec(inj(2),:) = DCoffsetInj;
+
+for iInj = 1:NumInj
+
+V(InjectionWindows(iInj,1):InjectionWindows(iInj,2),:)=bsxfun(@plus,V(InjectionWindows(iInj,1):InjectionWindows(iInj,2),:),DCoffsetVec(:,iInj)');
+
+end
+
+
+%% 
+
+Voff =mean(V(datastart:dataend,:));
+
+plot(bsxfun(@minus,V,Voff))
+
+
+
+
+
+
+
diff --git a/tests/filters/test_filters.m b/tests/filters/test_filters.m
new file mode 100644
index 0000000..5a6b8c7
--- /dev/null
+++ b/tests/filters/test_filters.m
@@ -0,0 +1,209 @@
+%%
+
+% Fs=16384;
+Fs=100000;
+
+Fc_all = 5:100:550;
+FreqNum = length(Fc_all);
+
+BWtarget =100;
+
+Fstopomin=0.5;
+Fcomin = 1.5;
+
+
+for iFreq = 1: FreqNum
+    
+    
+    Fc=Fc_all(iFreq);
+    
+    disp(Fc);
+    
+    if Fc == 25;
+        disp('stop');
+    end
+    
+    %% old filter
+    
+    BWold = BWtarget;
+    Fco(2) = Fc+(BWold/2);
+    Fco(1) = max([(Fc -(BWold/2)), Fstopomin]);
+    [B,A] = butter(3,Fco./(Fs/2));
+    
+    stold(iFreq)= isstable(B,A);
+    imp_lenold(iFreq)=impzlength(B,A,0.001);
+    imp_len_msold(iFreq) = (imp_lenold(iFreq) / Fs)*(1000);
+    len_cyc = (1/Fc)*1000;
+    imp_len_cycold(iFreq)=imp_len_msold(iFreq)/len_cyc;
+    
+    [Hold,W] =freqz(B,A,Fco(1):0.2:Fco(2),Fs);
+    fc_idx = find (W > Fc,1);
+    Gainold(iFreq) =abs(Hold(fc_idx));
+    
+    
+    %% New Bandpass
+    
+    BWbp = BWtarget;
+    
+    StopBandDiffbp = 350;
+    
+    Fpass1bp = Fc-BWbp/2;
+    
+    Astop1 = 60;
+    Astop2= 60;
+    Apass  = 0.5;
+    
+    
+    if ~(Fpass1bp > 0)
+        Fpass1bp = max([Fc/2 Fcomin]);
+        Astop1 = 15;
+    end
+    
+    Fstop1bp = Fpass1bp -StopBandDiffbp;
+    Fstop1bp = max([Fstop1bp Fstopomin]);
+    
+    Fpass2bp = Fc+BWbp/2;
+    Fstop2bp = Fpass2bp + StopBandDiffbp;
+    
+    
+    
+    hbp = fdesign.bandpass('fst1,fp1,fp2,fst2,ast1,ap,ast2', Fstop1bp, Fpass1bp, ...
+        Fpass2bp, Fstop2bp, Astop1, Apass, Astop2, Fs);
+    
+    dbp = design(hbp, 'butter', 'MatchExactly', 'passband', 'SOSScaleNorm', 'Linf');
+    
+    %     dbp = designfilt('bandpassiir', ...       % Response type
+    %         'FilterOrder', 20, ...
+    %        'HalfPowerFrequency1',Fpass1bp, ...    % Frequency constraints
+    %        'HalfPowerFrequency2',Fpass2bp, ...
+    %        'SampleRate',Fs)   ;            % Sample rate
+    %
+    
+    
+    
+    stbp(iFreq)= isstable(dbp);
+    imp_lenbp(iFreq)=impzlength(dbp,0.001);
+    imp_len_msbp(iFreq) = (imp_lenbp(iFreq) / Fs)*(1000);
+    len_cyc = (1/Fc)*1000;
+    imp_len_cycbp(iFreq)=imp_len_msbp(iFreq)/len_cyc;
+    
+    [Hbp,W] =freqz(dbp,Fstop1bp:0.2:Fstop2bp,Fs);
+    fc_idx = find (W > Fc,1);
+    Gainbp(iFreq) =abs(Hbp(fc_idx));
+    
+    
+    %% New Cascaded filter
+    
+    BWc = BWtarget;
+    
+    StopBandDiffc = 350;
+    
+    Fpass1c = Fc-BWc/2;
+    
+    Astop1 = 60;
+    Astop2= 60;
+    Apass  = 0.5;
+    
+    if ~(Fpass1c > 0)
+        Fpass1c = max([Fc/2 Fcomin]);
+        Astop1 =30;
+    end
+    
+    Fstop1c = Fpass1c -StopBandDiffc;
+    Fstop1c = max([Fstop1c Fstopomin]);
+    
+    Fpass2c = Fc+BWc/2;
+    Fstop2c = Fpass2c + StopBandDiffc;
+    
+    
+    
+    %     hchp = fdesign.highpass('fst,fp,ast,ap', Fstop1c, Fpass1c, Astop1, Apass, Fs);
+    %
+    %     dchp = design(hchp, 'butter', 'MatchExactly', 'passband', 'SOSScaleNorm', 'Linf');
+    
+    %     hclp = fdesign.lowpass('fp,fst,ap,ast', Fpass2c, Fstop2c, Apass, Astop2, Fs);
+    %
+    %     dclp = design(hclp, 'butter', 'MatchExactly', 'passband', 'SOSScaleNorm', 'Linf');
+    
+    dclp = designfilt('lowpassiir', ...       % Response type
+        'FilterOrder', 10, ...
+        'HalfPowerFrequency',Fpass2c, ...    % Frequency constraints
+        'SampleRate',Fs)   ;            % Sample rate
+    
+    
+    dchp = designfilt('highpassiir', ...       % Response type
+        'FilterOrder', 10, ...
+        'HalfPowerFrequency',Fpass1c, ...    % Frequency constraints
+        'SampleRate',Fs)   ;            % Sample rate
+    
+    
+    stc(iFreq)= all([isstable(dchp) isstable(dclp)]);
+    imp_lenclp(iFreq)=impzlength(dclp,0.001) ;
+    imp_lenchp(iFreq) = impzlength(dchp,0.001);
+    imp_len_msclp(iFreq) = (imp_lenclp(iFreq) / Fs)*(1000);
+    imp_len_mschp(iFreq) = (imp_lenchp(iFreq) / Fs)*(1000);
+    len_cyc = (1/Fc)*1000;
+    imp_len_cycclp(iFreq)=imp_len_msclp(iFreq)/len_cyc;
+    imp_len_cycchp(iFreq)=imp_len_mschp(iFreq)/len_cyc;
+    [Hclp,W] =freqz(dclp,Fstop1c:0.2:Fstop2c,Fs);
+    [Hchp,W] =freqz(dchp,Fstop1c:0.2:Fstop2c,Fs);
+    fc_idx = find (W > Fc,1);
+    Gainclp(iFreq) = abs(Hclp(fc_idx)) ;
+    
+    Gainchp(iFreq) =  abs(Hchp(fc_idx));
+    
+    %% Freq Info
+    
+end
+
+
+%%
+
+figure;
+hold on
+plot(Fc_all,10*log10(Gainold));
+plot(Fc_all,10*log10(Gainbp));
+plot(Fc_all,10*log10(Gainclp));
+plot(Fc_all,10*log10(Gainchp));
+hold off
+legend('Old','Bandp','Casclp','Caschp')
+% legend('Bandp','Casclp','Caschp')
+title('gain')
+xlabel('Freq')
+ylabel('Gain db');
+
+
+figure;
+hold on
+plot(Fc_all,imp_len_msold);
+plot(Fc_all,imp_len_msbp);
+plot(Fc_all,imp_len_msclp);
+plot(Fc_all,imp_len_mschp);
+
+hold off
+legend('Old','Bandp','Casclp','Caschp')
+% legend('Bandp','Casclp','Caschp')
+title('timedelay ms')
+ylabel('Settling time ms')
+xlabel('Freq')
+
+
+
+figure;
+hold on
+plot(Fc_all,imp_len_cycold);
+plot(Fc_all,imp_len_cycbp);
+plot(Fc_all,imp_len_cycclp);
+plot(Fc_all,imp_len_cycchp);
+
+hold off
+legend('Old','Bandp','Casclp','Caschp')
+% legend('Bandp','Casclp','Caschp')
+title('timedelay cycles')
+ylabel('Settling time cycles')
+xlabel('Freq')
+
+
+
+
+
diff --git a/tests/filters/test_singlecase.m b/tests/filters/test_singlecase.m
new file mode 100644
index 0000000..1c8f25a
--- /dev/null
+++ b/tests/filters/test_singlecase.m
@@ -0,0 +1,90 @@
+Fcur = 5;
+FreqNum = size(Fcur,2);
+
+Cycles = 64;
+T=(1./Fcur); %Period in s
+InjTime=(T.*Cycles);
+
+% InjTime=10;
+
+
+Amp_Inj = 500;
+Amp_Meas = 150;
+InjPhase=0;
+MeasPhaseDiff=-30;
+DCoffset = 0;
+DCoffsetinj = 0;
+
+% Fs=100000;
+Fs=16384;
+
+[Amp_error1, Phase_error1,V1,Vd1,Filt1,tr1] = check_acc( Fcur,InjTime,Amp_Inj,Amp_Meas,InjPhase,MeasPhaseDiff,DCoffset,DCoffsetinj,[],Fs);
+
+trim_demod=200;
+
+
+t1 = (0:1:length(V1)-1)/Fs;
+
+
+%%
+
+Fcur = 5;
+FreqNum = size(Fcur,2);
+
+Cycles = 32;
+T=(1./Fcur); %Period in s
+InjTime=(T.*Cycles);
+
+% InjTime=10;
+
+
+Amp_Inj = 500;
+Amp_Meas = 150;
+InjPhase=0;
+MeasPhaseDiff=-30;
+DCoffset = 100;
+DCoffsetinj = 500;
+
+% Fs2=100000;
+Fs2=16384;
+
+
+[Amp_error2, Phase_error2,V2,Vd2,Filt2,tr2] = check_acc( Fcur,InjTime,Amp_Inj,Amp_Meas,InjPhase,MeasPhaseDiff,DCoffset,DCoffsetinj,[],Fs2);
+
+
+t2 = (0:1:length(V2)-1)/Fs2;
+
+%
+
+%%
+
+fvtool(Filt1,Filt2)
+%%
+[H1,F1]=freqz(Filt1,Fcur-500:Fcur+500,Fs);
+[H2,F2]=freqz(Filt2,Fcur-500:Fcur+500,Fs2);
+
+figure;
+hold on
+plot(F1,10*log10(abs(H1)));
+plot(F2,10*log10(abs(H2)));
+hold off
+
+%%
+
+figure
+hold on
+plot(t1,V1)
+plot(t2,V2)
+% plot(V8)
+hold off
+
+
+figure
+hold on
+% plot(Vd1(tr1:end-tr1))
+plot(t1,Vd1)
+plot(t2,Vd2)
+% plot(Vd2(tr2:end-tr2))
+% plot(Vd8)
+hold off
+% ylim(Amp_Inj + [-1 +1])
diff --git a/tests/readme.md b/tests/readme.md
new file mode 100644
index 0000000..8f5515a
--- /dev/null
+++ b/tests/readme.md
@@ -0,0 +1,6 @@
+# Functional Tests
+
+`Test_ScouseTom.m` tests require the files from [![DOI](https://zenodo.org/badge/DOI/10.5281/zenodo.4783547.svg)](https://doi.org/10.5281/zenodo.4783547) to be placed in the data directory as the EEG system files get very big very quickly. These tests only check that the code is functional and are by no means unit tests.
+
+### Filter tests
+Checks to characterise frequency dependence of filtering method. Historical curiosity more than anything else.