summaryrefslogtreecommitdiff
path: root/triso_enrichment_driver.sh
diff options
context:
space:
mode:
Diffstat (limited to 'triso_enrichment_driver.sh')
-rwxr-xr-xtriso_enrichment_driver.sh26
1 files changed, 26 insertions, 0 deletions
diff --git a/triso_enrichment_driver.sh b/triso_enrichment_driver.sh
new file mode 100755
index 0000000..fc72c4f
--- /dev/null
+++ b/triso_enrichment_driver.sh
@@ -0,0 +1,26 @@
+#!/bin/bash
+
+prog=(./pin.py -e -f ucodyn)
+mkdir -p results/
+touch results/vwh.out
+touch results/exp.out
+
+for enr in $(seq 0.5 0.5 20.0); do
+
+ echo "Running for enrichment: ${enr}..."
+
+ # 1. Basic calculation without homogenization
+ "${prog[@]}" -t none -en $enr > results/explicit.out
+ exp_keff=$(cat results/explicit.out| grep "Combined k-effective" | cut -d "=" -f 2 | cut -d "+" -f 1)
+ exp_pm=$(cat results/explicit.out | grep "Combined k-effective" | cut -d "=" -f 2 | cut -d "-" -f 2)
+ echo $enr $exp_keff $exp_pm >> results/exp.out
+
+ # 2. VWH calculation
+ "${prog[@]}" -t vwh -en $enr > results/vwh.out
+ vwh_keff=$(cat results/vwh.out| grep "Combined k-effective" | cut -d "=" -f 2 | cut -d "+" -f 1)
+ vwh_pm=$(cat results/vwh.out | grep "Combined k-effective" | cut -d "=" -f 2 | cut -d "-" -f 2)
+ echo $enr $vwh_keff $vwh_pm >> results/vwh.out
+
+ echo "Found ${exp_keff} and ${vwh_keff}!"
+
+done