From 2806fbe59aea0dbd639043c818e2f80e5db5fe78 Mon Sep 17 00:00:00 2001
From: adrianaghiozzi <67669644+adrianaghiozzi@users.noreply.github.com>
Date: Thu, 24 Oct 2024 14:17:25 -0700
Subject: [PATCH 1/2] Read neutron rate output file in read_outputs

---
 src/outputs.jl | 31 ++++++++++++++++++++++---------
 1 file changed, 22 insertions(+), 9 deletions(-)

diff --git a/src/outputs.jl b/src/outputs.jl
index d2b10c9..1175635 100644
--- a/src/outputs.jl
+++ b/src/outputs.jl
@@ -6,14 +6,21 @@ Base.@kwdef mutable struct RABBIToutput
     torqdepo_data::Union{Array{<:Real},Missing} = missing
     rho_data::Union{Vector{<:Real},Missing} = missing
     time_data::Union{Vector{<:Real},Missing} = missing 
+    nrate_data::Union{Array{<:Real}, Missing} = missing
 
 end
 
 function read_outputs(path::String; filename::String ="run")
     struct_len = 4 
     result_path = abspath(joinpath(path, "$filename/beam1/rtfi_result_oav.bin"))
+    nrate_path = abspath(joinpath(path, "$filename/rtfi_nrate_oav.bin"))
 
     output = RABBIToutput()
+
+    function read_and_reshape(f, element_count; dims)
+        data = Float64[struct_unpack(read(f, struct_len)) for _ in 1:element_count]
+        return reshape(data, dims)
+    end
     
     open(string(result_path), "r") do f
         seekend(f)
@@ -28,11 +35,6 @@ function read_outputs(path::String; filename::String ="run")
             output.time_data = Float64[struct_unpack(read(f, struct_len)) for _ in 1:ntime]
             output.rho_data = Float64[struct_unpack(read(f, struct_len)) for _ in 1:nrho]
     
-            function read_and_reshape(f, element_count; dims)
-                data = Float64[struct_unpack(read(f, struct_len)) for _ in 1:element_count]
-                return reshape(data, dims)
-            end
-    
             bdens_data = read_and_reshape(f, nrho * ntime, dims=(nrho, ntime))
             press_data = read_and_reshape(f, nrho * ntime, dims=(nrho, ntime))
             output.powe_data = read_and_reshape(f, nrho * ntime, dims=(nrho, ntime))            
@@ -61,11 +63,22 @@ function read_outputs(path::String; filename::String ="run")
                 dArea = read_and_reshape(f, nrho * ntime, dims=(nrho, ntime))
                 output.torqdepo_data = read_and_reshape(f, nrho * nv * ntime, dims=(nrho, ntime, nv))
                 torqjxb_data = read_and_reshape(f, nrho * nv * ntime, dims=(nrho, ntime, nv))
-            
-            return output
-        
+                    
             break
             end
         end
-    end   
+    end 
+    
+    open(string(nrate_path), "r") do f
+        seekstart(f)
+
+        nrho = length(output.rho_data)
+        ntime = length(output.time_data)
+        while !eof(f)
+            output.nrate_data = read_and_reshape(f, nrho * ntime, dims = (nrho, ntime))
+        end
+    end
+
+    return output
+
 end
\ No newline at end of file

From c7253a21fa19cce52fbf256bd52128c74a3adcd3 Mon Sep 17 00:00:00 2001
From: adrianaghiozzi <67669644+adrianaghiozzi@users.noreply.github.com>
Date: Mon, 28 Oct 2024 13:29:48 -0700
Subject: [PATCH 2/2] Add runtests for RABBIT

---
 .github/workflows/runtests.yml     |  49 +++++++++++++++++++++++++++++
 test/Project.toml                  |   2 ++
 test/run/beam1/rtfi_result_oav.bin | Bin 0 -> 5134 bytes
 test/run/rtfi_nrate_oav.bin        | Bin 0 -> 800 bytes
 test/runtests.jl                   |  17 ++++++++++
 5 files changed, 68 insertions(+)
 create mode 100644 .github/workflows/runtests.yml
 create mode 100644 test/Project.toml
 create mode 100644 test/run/beam1/rtfi_result_oav.bin
 create mode 100644 test/run/rtfi_nrate_oav.bin
 create mode 100644 test/runtests.jl

diff --git a/.github/workflows/runtests.yml b/.github/workflows/runtests.yml
new file mode 100644
index 0000000..192e0d6
--- /dev/null
+++ b/.github/workflows/runtests.yml
@@ -0,0 +1,49 @@
+name: CI
+
+on:
+  push:
+    branches:
+      - master
+  pull_request:
+
+jobs:
+  test:
+    name: Julia ${{ matrix.version }} - ${{ matrix.os }} - ${{ matrix.arch }}
+    runs-on: ${{ matrix.os }}
+    strategy:
+      fail-fast: false
+      matrix:
+        version:
+          - "1.x"
+        os:
+          - ubuntu-latest
+        arch:
+          - x64
+
+    steps:
+      - uses: actions/checkout@v4
+
+      - uses: julia-actions/setup-julia@latest
+
+      - uses: julia-actions/cache@v1
+
+      - name: Add FuseRegistry
+        run: |
+          rm -rf ~/.julia/registries/FuseRegistry
+          julia -e 'using Pkg; Pkg.Registry.add(RegistrySpec(url="https://github.com/ProjectTorreyPines/FuseRegistry.jl.git")); Pkg.Registry.add("General"); Pkg.Registry.update()'
+
+      - name: Replace git@github.com with https in Package.toml files
+        run: |
+          find ~/.julia/registries/FuseRegistry -type f -name 'Package.toml' -exec sed -i 's|git@github.com:|https://project-torrey-pines:${{secrets.PTP_READ_TOKEN}}@github.com/|g' {} +
+
+      - name: Install dependencies
+        run: |
+          julia -e '
+            using Pkg
+            Pkg.develop(PackageSpec(path=pwd()))
+            Pkg.instantiate()
+          '
+
+      - uses: julia-actions/julia-buildpkg@v1
+
+      - uses: julia-actions/julia-runtest@v1
diff --git a/test/Project.toml b/test/Project.toml
new file mode 100644
index 0000000..0c36332
--- /dev/null
+++ b/test/Project.toml
@@ -0,0 +1,2 @@
+[deps]
+Test = "8dfed614-e22c-5e08-85e1-65c5234f0b40"
diff --git a/test/run/beam1/rtfi_result_oav.bin b/test/run/beam1/rtfi_result_oav.bin
new file mode 100644
index 0000000000000000000000000000000000000000..21c839b8fa9a0767ddbd514ced8a64c6e6f4d715
GIT binary patch
literal 5134
zcmeHJ`8$<c8|E02DMN-Rp(shl4EtScy|(=hn+BCK?)`2hk~E^C$y_2snT3)J6%jIo
z5K)FgN}M_xlujy9j#K1w{)g|ouJ4C)U)OUz>$!ipuV>vq+-nI530VsXU1HBQYoX=-
z4}|8>+0g;g)6(GF5t*8rA}=p5)Y%a;*Yl8&Py<R&dp^fc%+2Dssi`6^B&3OVcC5zJ
z)2whWFKUi+nB&~%INv!g5RYj4!dEq!;9r+-<bRRk@MUX_`S>Fheg$uq*QD3S+jhU4
z7twx&XOPsvvkkn(OYu+VRdj{$F4>d3%JnmD?#h{Nb(Ji)#K;m?7c}3srhc41^zA!e
zwzH9c)SJ!ceK+LCA5rF`{j)p`g?=9Mv683edyRJ_y`87Dyp%_gX*^@65FSz_c$=ku
zxE;Nk<_4$C-2M<RbZv?G<uWDJOm=p*l4n=kAk%hjB12oPNPQ1=(xd4oarXTKg4~x&
zG*3nmvWp%Pbz#?u_3<%8uevjFq)eIccCX+FlI=O}v2ShCbnR`dC9>HDYAxid$TreE
z@H%<GW&_DOXHLGFR3+E?{vZZ&+KJiSc|^P6Ibx6NL*h|OA>pDEO;8*s!bx12;QE(w
zj_IxAY+e4srZR=Ji4BftFYH{$^?kOLyF!cL4qX!A?g|m%Htd|HHB~>;QlT&DE#95<
z_nocuWffQY#egOK&Tc6k=O;?PtnQ>fUOq?dIR>Z&%g4!-lR>0^M=z22b3ONJm?L+r
zo#1}?_8Zss_ushZmrm0ov17Dl@f-TmKsRl@w2OAl+D2b0v7ocNHRxp<#c12eM^yK^
zSnAzzj4B-bMBev0NT&BcAf#G8a&OQ3le=weDR)S085b2PatG>v&~B_=I=HHYZthB=
z%hV6kxmlI;M#&u7Aw7n^#2lbi*I3bJ-#$}5SxJ<RmMN84R6stco+Jd~Bi!@T!`wrn
z#oXD^rQES6a$NIMleA(@H=Q<KN{^jLp|x-CrDv+E>BOhmbZgUDTH1Uc9iwhRFPr$2
zYMDr&wrw$>;<j8SO<q1HjtAbMW5yM^3#=D&<F0qm?JFy2cTOVh`EoZMU~EPYno7_;
zEsv=cVPB~bzgDUsIhl%GW>2jzm!n)1pOI6MDP$<Vkz60tM_ku3A({`>(7CdT+`nF`
zayu8c(V+`V=`Onj+U~m-J@~gNeX4Li{o9*oRI2Y7mC$mJGT)U*?UCF}wbn1B{3ZrS
z8R--<mTgNKinkL64qAktIh!82a+L;ix?4w1(X(DYv`#5O4^J$ioeCx?@r+tZSSNwn
zyrYu3awv@w&E8F&_S2#kzJEmu&gGNUPkhM9;~JzUKZh_279(PGt?9+-xpd-|+qC>l
zAie*?PI@HDnm*H^P8%yuQig%`lvHIBrP+6z8cIy2-fDYM0Z+9lNr@41xA!%2sCPH1
z9;8M-qSJ}E@G*{cUIGL?i_lJuSX6bi1sPRMqvCywaJOMoyfS+O{_dm~KIM55|D1Lf
zukuX74@PF=-J=CK8C8l)+E(L0t^psp-GZxCKEQh)w&MS7y*u6!0!`5vL3(Sd0JMb(
zV)Bm)PKs|A<p1Ib_7vy|*4N7mT>Yn*ai<PuRx_7*bU%=>a?oQ!njiAN3_9}5BJc9l
zm&)-vcWig#U;O0SmZj)^!C%3>b6c~ZAh?k+{Zz>;=()<YJkDd1MJ_YrOctXooX32g
zEn>8PD`C_^Z!!9KA+!A61!fewGUf->n3R-uzQsbCe=0AQmn$yqUVdD{y=3JTLEdL0
zCc4>#$+J;njIKy9^PfvFQj?O*V|p=T|CD0R{W4{?Z!}}<o!prtIa`=x`{fvFc?W;C
zJDUG;Op!ku9Kb7W3l#*|L<wA)k_0kZA%dGX{RC5Y_=0&o9KqYK`hxyUMM1FR4@Mhy
zFnUieGd3zmn71OUm}1>7e#oFJU!L8_E2)y<J=1b<tIFtcO_5Y|U-MPoJ>0NaFc#jz
zcy6y@48rpnndkzBKb*_xuFYW_6!RI!4@Hc6e=#E`Rm!Mb%V+xcB`|)Tj!fEwGPCbz
z3;*Ir#8;WiZHr_n_ZPAf?iEovf~?&pOkIj4qmZ$PxiTTe*l(3)0vhF+j1&##%nOQ1
zy=TI_|760X6funE)lJN%JQ-&8NGqTG^E5w4b^%`_+mja}uYjL0RKwMK)bM^uRb2A5
z0-m}~1`k>-f$y*r#<%mwk;aC1NG9<WT5S6Uxn6sR4AnlPh2Lk8?bv+WX_Yem!cPx(
zT0`M-k6iIxd4B>1t4Hfqo})uuLr6|@3>Ba4M=P$~Mon*W(ST7X+Cb@{72>k!fYcOp
z2+xPvyK3O$_o65wb1T}TlY+{3w;*}T*Jxw-U+Bp;rMf8)TKC$LgLa)eg|dHm&qd&V
zl;#nM%tAa+!w(FpzAuBtk3%qzwc&87wh43?Syu-d+iI5@^}xOn9{R&B0}VdEiV{U*
z(WtfpUIEnbS3_#JkBTbZHmQL73d`Uz((`a*-C0z2@+*pK8Ac&OuaJk)D|GPGF!II6
zP~3W9{H3Wh{?1Va|CqiKm#^aB=2x8Y*XP~?%`bP*wTg%6W<?h|CH)r3`gJ28iE<>%
z<f8pPVJMO|Kz`Elh)nn!_TS2fnf*0z-A@TLcGVdf6KSYKsTD<;yg>`+?prdcRL5sU
z*S!gXNZWlkddqOp>3c57<U<H*uk}QQWdKbSRzXX-9&D3)IQ&vz4A*VmSeLN-LG8-B
zgYa|IcI2FJ5ji&IAU~T(RDCMN)~PDVR_oFM+wJT_w!`n|*@gSgvs=Hq)VjMb*Lv;g
z6D5h+y(MC#o}G7tp4~rxlz$*OR{{HWdU|R8gNT(+v6MYKz-o_faHv)u7EG|<%9|dr
z?^^^Un$lpEaS=2$Yk(eAJ#ca9Tlly1Bn)Om&}SDZwEC+8N?)XgI@>f+YmYAKmD5N6
zZN1CtTaain52fBefJ{~GkzBtGdRJwL{QXr?x3dVEJ~;%Zch*Ct;SAUq8VHxr5Iz<a
zg4!x+z)4pdyvmQq0{jP9tTBDoGmYC;san}6<3%RIW*kuO>{e*M(+*aU=5Tp{IW%%G
zgawKEFn`zt?r&tl{w)CBJdR*5X$G~;m%wG>BY>lk3zXJ5116vsi@l(NT@yZxd>`yb
zvU?Yx_|8m_W}XcsG$MekuODddIs}w227uY5D3Bgd41#?#f$)}0AeVj*$a_@+u_r-5
zX|*v3y*7?Tcm`uKEBaVlr#7H%&>77W@kIwl?a;_LLZ7t^&}mW?waba2#%*t4QqNs@
z;&BFS{}}+2tuVMCP6(QoXMh=?17d{FV<(<GW%(#-v7AFLTaE6`La*4F$a%m4)h>2{
z=lj<~@vjyz!^{$52TkFPTZT|!ff=;A!Gf<BVlaFZK;2{$*m_3|N`H9+UhTRBx*WEE
zWo-|!KR(RI^sEBVS+V`7?xidWoXG<5#<`&V@fjc~6A08hj{>S87zk%Yg9RsxLAJ|9
z@Lb~}IKH9@Y`k6p+-CxS&k_UB`uZc*w)rS_f@x*3V~26xntLqQQ7(uqm%@hlvS7f?
z7_e^kVTQ>Tpx40+h(3}7@6M=##M9OQbiTlpRFg2J-1C_2yxkaYxh2a$`ZDXUoF>+D
z|4XbRcUjg?<Rn19VKsJST?kmy1VC|YFt85N0y=r3;J(>L5L9mjZg?I7JtFboyG1x?
z_DurybSfy)hywA?c7dd)Yd~|(?^r8li*@@-S*utSVWlcH;MRyV#1@3rWY&JfYKC`W
zsT=iy+`bX)npQh@N9TQ2vEB|4sWF26G(Q70J>Oy9My)`mk1|-yOk!!tF<6fOFV@KO
zxA>;`TGr$p8iYwoU}K-8L3yzu_~_M#-95(!)67QTX}k#hTD=&ID6>KB+d-^kCILIf
zkHb1HcwzO)R;(3*9M-NUO{@Vji{;<KvR-Z?0R;Kg*x31S5OtUagGVDkpv`K4FO>t|
z%=UxA?<^2_>kz2r$ARJKFfdb<2-I9tK>fkfAgXL9&<xTCz72htuA>!pUsuYSYfyyQ
zH&%mq5gF+1hSi*m{u>(~55q(hOo1al59kjKVx=djs#Sg4fu;X2mdZX2qDx<6`E?ec
zZ$J@b_kY36$|JCIgYQ}8&Z}UIX%t*Dk_zR@UEo?}C)h7#3qO3ffU`NvVbWb`xU^sj
zn15uzGJ8w##!~@kg8|I5rVyj*;xOwHA8cNs6?SdY0E_KHv-m-JP<ZqV)G17c^`9J}
z=-^hU{yPnu4qHO~R$VBTv;g{)3B!XiE&yCu3)Wc6gClM)u`eklSiDanW*xN`GkIu%
z9WCu<!8X7$ZJ3Wv6vd;_*~>^`$N&i|Xdvxr5meOL506+D!;9-q!_WIUFac38ZbA>v
z-k1iF3kyJXSp;~N6AbKQc)(CiA2h$I#N;$&u_w4VVrR#oG=&TlsId|)GGBzs9||GC
z*FGq=uNc}*pMjNt4Mo+t@Nlv&R7jlyKJtZNKsO2u1crhQ0d63SSOJ!}+{A1$ez1O<
zOKYz(9g^^}Aw5J*$SCq9ajd?ch%8Da8U(&X{fIGfzfF`lkk!u7_^M27-#g8D1v@w^
zcVuwFJ)Agu_6c!vS0vh0NE+CrT)n~G%vNFdtkfmly=};@d?PX_WPmUdyF+-VClSGQ
z-o(k9Mnr0~81Z4Gi{r1VLezfwiz6BNfHTvO!P#`znWJ-9h=X*`+k6*UY4c7$mpxoL
ZkNy9?-#@S3zxdaIe;xSOf&YaA{{g{~yJ-Lb

literal 0
HcmV?d00001

diff --git a/test/run/rtfi_nrate_oav.bin b/test/run/rtfi_nrate_oav.bin
new file mode 100644
index 0000000000000000000000000000000000000000..2eb68f4e46f15f516548336274de82814fc75054
GIT binary patch
literal 800
zcmV+*1K<4N--cETl7v?MrEpdS^d(jc@(NZJklIyaLZns4hjdkm7CBYvbM;gY61G(c
zNt;zvplMan$Q@O<*1c4xgh*5ej=)oGR4-GvVwX|{=IBwRE{0ZjdV^LJ2yIr9<{?%e
zCI?nW57ku)_n}o{q;yrK&N)>$>f=<znZi|bs+d)JKV?;SnHyEaN4r$kjYL#IX1Y_3
zUnf(b)Q(avC8<&M6@yf{qI6W|E+teWu2@sjeKAw2AOTaW1F=%}Q*Tlp#VArifBsQ%
zBZX2|076nj4hB*uSff#MnjcX_=wDFPP>WBJlABIy@or5i0tif8)Pz*yqjyw(B_>oM
zLR3?{z${Zx1N>41kg!r)k9AVuU@cPm%>`1!GlfzOlRi?~i33t@;-OJpyB$%pVOmgy
z0EJKRNRduj)LBj5`{GN0?|@WQi&s?om?cyd{^?W7fRR&SZ6s4T5UW!CPAgIYPq$Io
zX{S&R;AvE#uWD4twKG)CP5D#snVeJhb}~~Qrn*vBpgK~Qr^HbJXe&`=O@36mSyEId
z9w1cyBH&Yw|A|wXfErWMsGm~d^dC~FzO7M@t3pt#evVWEbY)Z=11?mC|L#)&q?A+E
z&?i$)O|DW}I4)A|A-YktXTnena;-v$d#Xb2yP86FX;nf$t};T2hzvr}JJUh!I-Nn*
z1a3jX(j7r_BkDnTvByCUIF>;;=t@D+AoM^??tnn}p7%dfcws-Ot-?M-z!g3pj;ule
z6sba$XqZA%Gf+Zc%q>E({s=;VNYFu2O`SoQ@oquzixNS6K>9%*Q^Y}CXOltb7)U{o
z*6=`)!+b!HwC_LlJXb$Kt+zf-gxNhi*yus?-Q+>3xz9nyKaW9X-)})gsXRe=$qhk%
z6VpH+#GgP?-%~(TDIGy+3I;*&smVaAGKfIy-6KHL6|FxSHYPvFSC>AmtO7nC2}wOa
zrsqLQ`{F@un#@5$z>7hBQE5Rm={G@n&I>^=EYm;>W1v9XyF@@S`X@n>#{)rc;>AGM
ePliB}s3Abusi{Be<sd(kd5}JAl=eL%0TMlreqxaT

literal 0
HcmV?d00001

diff --git a/test/runtests.jl b/test/runtests.jl
new file mode 100644
index 0000000..3566742
--- /dev/null
+++ b/test/runtests.jl
@@ -0,0 +1,17 @@
+using RABBIT
+using Test 
+
+@testset "RABBIT" begin 
+
+    works = try
+        output = RABBIT.read_outputs(@__DIR__)
+        @test isapprox(output.powe_data[:,1][1], 208692, rtol=0.10)
+        @test isapprox(output.powi_data[:,1][1], 482856, rtol=0.10)
+        @test isapprox(output.torqdepo_data[:, :, 1][1], 0.20518352, rtol=0.10)
+        @test isapprox(output.nrate_data[:, :, 1][1], 29.10558, rtol=0.10)
+        true
+    catch
+        false
+    end
+    @test works==true
+end
\ No newline at end of file