#!/bin/bash
# SIM-003 v6c WAVE RE-INJECTION — Hetzner deployment
# Date: 2026-04-01
# Pockets + voids + emergent global. Cascade only at global eq.
set -e
cd ~/prometheus_sim
mkdir -p results_sim003_v6c logs

TS=$(date -u +%Y%m%d_%H%M%S)
LOG="logs/v6c_wave_reinjection_${TS}.log"

echo "============================================" | tee "$LOG"
echo "SIM-003 v6c: WAVE RE-INJECTION"               | tee -a "$LOG"
echo "$(date -u)"                                    | tee -a "$LOG"
echo "Pockets + voids, emergent global_r"            | tee -a "$LOG"
echo "Local overflow → wave continues (no cascade)" | tee -a "$LOG"
echo "Global overflow (mean≥0.5) → cascade"         | tee -a "$LOG"
echo "============================================" | tee -a "$LOG"

# Phase 1: Fe at 5000 cycles — enough for 2D to reach global eq?
echo "" | tee -a "$LOG"
echo "=== PHASE 1: Fe 5000 cycles ===" | tee -a "$LOG"
python3 SIM-003_v6c_wave_reinjection.py \
    --element Fe \
    --cycles 5000 \
    --max-dim 3 \
    --delta-r 0.0001 \
    --grid-1d 2048 \
    --grid-2d 256 \
    --grid-3d 128 \
    2>&1 | tee -a "$LOG"

# Phase 2: Fe extended at 20000 cycles if 2D hasn't equilibrated
echo "" | tee -a "$LOG"
echo "=== PHASE 2: Fe 20000 cycles ===" | tee -a "$LOG"
python3 SIM-003_v6c_wave_reinjection.py \
    --element Fe \
    --cycles 20000 \
    --max-dim 3 \
    --delta-r 0.0001 \
    --grid-1d 2048 \
    --grid-2d 256 \
    --grid-3d 128 \
    2>&1 | tee -a "$LOG"

# Phase 3: Key elements at 5000 cycles
echo "" | tee -a "$LOG"
echo "=== PHASE 3: KEY ELEMENTS 5000 cycles ===" | tee -a "$LOG"
for elem in H He Li C Na Cu W Au; do
    echo "--- $elem ---" | tee -a "$LOG"
    python3 SIM-003_v6c_wave_reinjection.py \
        --element "$elem" \
        --cycles 5000 \
        --max-dim 3 \
        --delta-r 0.0001 \
        --grid-1d 2048 \
        --grid-2d 256 \
        --grid-3d 128 \
        2>&1 | tee -a "$LOG"
done

echo "" | tee -a "$LOG"
echo "ALL COMPLETE: $(date -u)" | tee -a "$LOG"