Skip to content

Commit

Permalink
Add test for log(ft) to test driver
Browse files Browse the repository at this point in the history
  • Loading branch information
@jaharris87
jaharris87 committed Jul 2, 2024
1 parent bd5d806 commit 76aab72
Show file tree
Hide file tree
Showing 4 changed files with 126 additions and 0 deletions.
59 changes: 59 additions & 0 deletions test/Data_SN231_logft/ab_co
View file
Original file line number Diff line number Diff line change
@@ -0,0 +1,59 @@
50% C12, 49% O16, 1% Ne22 (Ye=0.4995)
n 0.0000000E+00 p 0.0000000E+00 d 0.0000000E+00 t 0.0000000E+00
he3 0.0000000E+00 he4 0.0000000E+00 li6 0.0000000E+00 li7 0.0000000E+00
be7 0.0000000E+00 be9 0.0000000E+00 b8 0.0000000E+00 b10 0.0000000E+00
b11 0.0000000E+00 c11 0.0000000E+00 c12 4.1250000E-02 c13 0.0000000E+00
c14 0.0000000E+00 n12 0.0000000E+00 n13 0.0000000E+00 n14 0.0000000E+00
n15 0.0000000E+00 o14 0.0000000E+00 o15 0.0000000E+00 o16 3.0937500E-02
o17 0.0000000E+00 o18 0.0000000E+00 o19 0.0000000E+00 f17 0.0000000E+00
f18 0.0000000E+00 f19 0.0000000E+00 f20 0.0000000E+00 f21 0.0000000E+00
ne18 0.0000000E+00 ne19 0.0000000E+00 ne20 0.0000000E+00 ne21 0.0000000E+00
ne22 4.5454545E-04 ne23 0.0000000E+00 ne24 0.0000000E+00 na20 0.0000000E+00
na21 0.0000000E+00 na22 0.0000000E+00 na23 0.0000000E+00 na24 0.0000000E+00
na25 0.0000000E+00 na26 0.0000000E+00 mg22 0.0000000E+00 mg23 0.0000000E+00
mg24 0.0000000E+00 mg25 0.0000000E+00 mg26 0.0000000E+00 mg27 0.0000000E+00
mg28 0.0000000E+00 mg29 0.0000000E+00 al24 0.0000000E+00 al25 0.0000000E+00
al26 0.0000000E+00 al27 0.0000000E+00 al28 0.0000000E+00 al29 0.0000000E+00
al30 0.0000000E+00 si26 0.0000000E+00 si27 0.0000000E+00 si28 0.0000000E+00
si29 0.0000000E+00 si30 0.0000000E+00 si31 0.0000000E+00 si32 0.0000000E+00
si33 0.0000000E+00 p28 0.0000000E+00 p29 0.0000000E+00 p30 0.0000000E+00
p31 0.0000000E+00 p32 0.0000000E+00 p33 0.0000000E+00 p34 0.0000000E+00
p35 0.0000000E+00 s30 0.0000000E+00 s31 0.0000000E+00 s32 0.0000000E+00
s33 0.0000000E+00 s34 0.0000000E+00 s35 0.0000000E+00 s36 0.0000000E+00
s37 0.0000000E+00 cl32 0.0000000E+00 cl33 0.0000000E+00 cl34 0.0000000E+00
cl35 0.0000000E+00 cl36 0.0000000E+00 cl37 0.0000000E+00 cl38 0.0000000E+00
cl39 0.0000000E+00 ar34 0.0000000E+00 ar35 0.0000000E+00 ar36 0.0000000E+00
ar37 0.0000000E+00 ar38 0.0000000E+00 ar39 0.0000000E+00 ar40 0.0000000E+00
ar41 0.0000000E+00 ar42 0.0000000E+00 k36 0.0000000E+00 k37 0.0000000E+00
k38 0.0000000E+00 k39 0.0000000E+00 k40 0.0000000E+00 k41 0.0000000E+00
k42 0.0000000E+00 k43 0.0000000E+00 k44 0.0000000E+00 ca38 0.0000000E+00
ca39 0.0000000E+00 ca40 0.0000000E+00 ca41 0.0000000E+00 ca42 0.0000000E+00
ca43 0.0000000E+00 ca44 0.0000000E+00 ca45 0.0000000E+00 ca46 0.0000000E+00
ca47 0.0000000E+00 ca48 0.0000000E+00 sc41 0.0000000E+00 sc42 0.0000000E+00
sc43 0.0000000E+00 sc44 0.0000000E+00 sc45 0.0000000E+00 sc46 0.0000000E+00
sc47 0.0000000E+00 sc48 0.0000000E+00 sc49 0.0000000E+00 ti42 0.0000000E+00
ti43 0.0000000E+00 ti44 0.0000000E+00 ti45 0.0000000E+00 ti46 0.0000000E+00
ti47 0.0000000E+00 ti48 0.0000000E+00 ti49 0.0000000E+00 ti50 0.0000000E+00
ti51 0.0000000E+00 v45 0.0000000E+00 v46 0.0000000E+00 v47 0.0000000E+00
v48 0.0000000E+00 v49 0.0000000E+00 v50 0.0000000E+00 v51 0.0000000E+00
v52 0.0000000E+00 v53 0.0000000E+00 cr47 0.0000000E+00 cr48 0.0000000E+00
cr49 0.0000000E+00 cr50 0.0000000E+00 cr51 0.0000000E+00 cr52 0.0000000E+00
cr53 0.0000000E+00 cr54 0.0000000E+00 cr55 0.0000000E+00 mn49 0.0000000E+00
mn50 0.0000000E+00 mn51 0.0000000E+00 mn52 0.0000000E+00 mn53 0.0000000E+00
mn54 0.0000000E+00 mn55 0.0000000E+00 mn56 0.0000000E+00 mn57 0.0000000E+00
fe51 0.0000000E+00 fe52 0.0000000E+00 fe53 0.0000000E+00 fe54 0.0000000E+00
fe55 0.0000000E+00 fe56 0.0000000E+00 fe57 0.0000000E+00 fe58 0.0000000E+00
fe59 0.0000000E+00 fe60 0.0000000E+00 co52 0.0000000E+00 co53 0.0000000E+00
co54 0.0000000E+00 co55 0.0000000E+00 co56 0.0000000E+00 co57 0.0000000E+00
co58 0.0000000E+00 co59 0.0000000E+00 co60 0.0000000E+00 co61 0.0000000E+00
ni55 0.0000000E+00 ni56 0.0000000E+00 ni57 0.0000000E+00 ni58 0.0000000E+00
ni59 0.0000000E+00 ni60 0.0000000E+00 ni61 0.0000000E+00 ni62 0.0000000E+00
ni63 0.0000000E+00 ni64 0.0000000E+00 cu57 0.0000000E+00 cu58 0.0000000E+00
cu59 0.0000000E+00 cu60 0.0000000E+00 cu61 0.0000000E+00 cu62 0.0000000E+00
cu63 0.0000000E+00 cu64 0.0000000E+00 cu65 0.0000000E+00 zn59 0.0000000E+00
zn60 0.0000000E+00 zn61 0.0000000E+00 zn62 0.0000000E+00 zn63 0.0000000E+00
zn64 0.0000000E+00 zn65 0.0000000E+00 zn66 0.0000000E+00 zn67 0.0000000E+00
zn68 0.0000000E+00 ga62 0.0000000E+00 ga63 0.0000000E+00 ga64 0.0000000E+00
ga65 0.0000000E+00 ga66 0.0000000E+00 ga67 0.0000000E+00 ga68 0.0000000E+00
ga69 0.0000000E+00 ge65 0.0000000E+00 ge66 0.0000000E+00 ge67 0.0000000E+00
ge68 0.0000000E+00 ge69 0.0000000E+00 ge70 0.0000000E+00
20 changes: 20 additions & 0 deletions test/Test_Problems/setup_logft_sn231
View file
Original file line number Diff line number Diff line change
@@ -0,0 +1,20 @@

## Output Controls
3 Diagnostic Output Level
2 Per Timestep Output Level
# ASCII output filename root, network will append zone number
Test_Results/ev_logft_sn160_
# Binary output filename root, network will append zone number
Test_Results/ts_logft_sn160_
# Species to output in ASCII output (format 14a5): 14
p he4 c12 o16 ne20 mg24 si28 ca40 ti44 cr48 fe52 fe54 ni56 zn60
## Input Controls
# Nuclear Data Directory
Data_SN231_logft
# Initial Abundance and Thermodynamic Trajectory Files
Data_SN231_logft/ab_co
Test_Problems/th_co_burn_1
Data_SN231_logft/ab_co
Test_Problems/th_co_burn_2
Data_SN231_logft/ab_co
Test_Problems/th_co_burn_3
32 changes: 32 additions & 0 deletions test/test_settings_logft
View file
Original file line number Diff line number Diff line change
@@ -0,0 +1,32 @@
## Problem Description
2000 timesteps
diagnostic level 1
New parameters
## Job Controls
1 Initial Zone
3 # of Zones
2 Include Weak Reactions (yes=1,no=0,only=-1)
1 Include Screening (yes=1)
1 Process Nuclear Data at Run Time (yes=1,no=0)
## Neutrinos
0 Include Neutrino Reactions (yes=1, no=0)
## NSE Initial Conditions
11.0 Temperature in GK to use NSE initial conditions instead of file
## Integration Controls
1 Choice of integration Scheme (1=Backward Euler, 2= Bader-Deufelhard)
9999 Max. number of timesteps before quit
5 Max. iterations per step
1 Rebuild the jacobian every ijac iterations after the first
0 Convergence Condition (Mass Cons.=0, (dY/Y small)=1)
1.00E-01 Max. Abundance Change per timestep
1.00E-07 Smallest Abundance used in timestep calculation
1.00E-06 Mass Conservation Limit
1.00E-04 Convergence Criterion
1.00E-30 Lower Abundance limit, smaller abundances = 0
2.00E+00 Max. Factor to change dt in a timestep
## Self-heating Controls
1 Include self-heating (yes=1,no=0)
1.00E-02 Max. Temperature Change per timestep
1.00E-04 Temperature Convergence Criterion
## Zone Batching Controls
2 Blocking size for zone loop
15 changes: 15 additions & 0 deletions test/test_xnet.sh
View file
Original file line number Diff line number Diff line change
Expand Up @@ -45,6 +45,7 @@
## torch47 (47) 52
## SN160 (160) 53
## SN160 (160) (BDF integrator) 54
## SN231 (231) (log-ft rates) 55
## =====================================================================
## All parallel test problems 30
## 4 different SN alpha (16) 31
Expand Down Expand Up @@ -109,6 +110,13 @@ function do_test_bdf {
test_diff Test_Results/net_diag_$2 Test_Problems/Results/net_diag_$2 $2
}

function do_test_logft {
cat test_settings_logft Test_Problems/setup_$2 >| control
$1
mv -f net_diag01 Test_Results/net_diag_$2
test_diff Test_Results/net_diag_$2 Test_Problems/Results/net_diag_$2 $2
}

function do_test_batch {
cat test_settings_batch Test_Problems/setup_$2 >| control
$1
Expand Down Expand Up @@ -297,6 +305,13 @@ for xnet in ${xnet_list[@]}; do
do_test_bdf $xnet $test_name
fi

# Self-heating test using SN160 (explosive burning of degenerate C/O) (with BDF integrator)
if [ $itest -eq 50 -o $itest -eq 55 ]; then
echo "Test: Self-heating from explosive burning of degenerate C/O with SN231 network and log(ft) rates"
test_th="logft"; test_net="sn231"; test_name=${test_th}_${test_net}
do_test_logft $xnet $test_name
fi

# Zone batching test using alpha (explosive burning of degenerate C/O)
if [ $itest -eq 60 -o $itest -eq 61 ]; then
echo "Test: Zone batching w/ self-heating from explosive burning of degenerate C/O with alpha network"
Expand Down

0 comments on commit 76aab72

Please sign in to comment.