#!/bin/bash
# Reconstruction script for pol2 example.

# Initial stage.
gmm-rec.py pol2_input.mrc -v -o coarse.mrc -s 2000 50 50 -g -r 5000 -d 32 -a 4 -A 2

gmm-plot-stack.py -i coarse_input.mrc coarse_initial_projections.mrc coarse_final_projections.mrc -r 4 -c 6
gmm-plot-stats.py -i coarse_stats.csv -o coarse

# Refinement stage.
gmm-rec.py pol2_input.mrc -v -o refine.mrc -s 1000 -K 500 -t coarse_transformations.csv -N 10000 -a 4 -A 2

gmm-plot-stack.py -i refine_input.mrc refine_initial_projections.mrc refine_final_projections.mrc -r 4 -c 6
gmm-plot-stats.py -i refine_stats.csv -o refine
