From e8d6bb189ce835c27bf041ad47655cbf3f03cb5e Mon Sep 17 00:00:00 2001 From: connorhhmoore Date: Wed, 2 Sep 2026 16:26:20 +0200 Subject: Added driver script for dynamic enrichment --- triso_enrichment_driver.sh | 26 ++++++++++++++++++++++++++ 1 file changed, 26 insertions(+) create mode 100755 triso_enrichment_driver.sh (limited to 'triso_enrichment_driver.sh') 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 -- cgit v1.2.3