#! /bin/bash 
# Last edited on 2020-01-17 14:23:02 by jstolfi

# Plots the phase diagram of the mean-field evolution map for mu = 0, WE = WI.

export PATH=".:..:${PATH}"

dead="$1"; shift;   # Dead step option.
dense="$1"; shift;  # Fill plot with F-plots in a regular array.
region="$1"; shift; # Region to highlight, or -1 for all.

name=`printf "dead_%d_phases" "${dead}"`
echo "dead = ${dead} name = ${name}"

if [[ ${dead} -eq 0 ]]; then
  # No dead step:
  kappa_min="-2.7"
  kappa_max="+2.7"
  lambda_min="-2.2"
  lambda_max="+3.2"
else
  # Dead step:
  if [[ ${region} -lt 0 ]]; then
    # All regions:
    kappa_min="-4.7"
    kappa_max="+8.7"
    lambda_min="-4.7"
    lambda_max="+8.7"
  elif [[ ${region} -eq 4 ]]; then
    kappa_min="-2.5"
    kappa_max="+0.5"
    lambda_min="-0.5"
    lambda_max="+2.5"
  elif [[ ${region} -eq 5 ]]; then
    kappa_min="-0.5"
    kappa_max="+1.5"
    lambda_min="-0.5"
    lambda_max="+1.5"
  elif [[ ${region} -eq 6 ]]; then
    kappa_min="-0.5"
    kappa_max="+2.0"
    lambda_min="-0.5"
    lambda_max="+2.0"
  elif [[ ${region} -eq 7 ]]; then
    kappa_min="-0.5"
    kappa_max="+2.5"
    lambda_min="-0.5"
    lambda_max="+2.5"
  elif [[ ${region} -eq 13 ]]; then
    kappa_min="+0.5"
    kappa_max="+4.5"
    lambda_min="-3.5"
    lambda_max="+0.5"
  else
    # Generic region:
    kappa_min="-4.7"
    kappa_max="+8.7"
    lambda_min="-4.7"
    lambda_max="+8.7"
  fi
fi

compute_mean_phases.gawk \
    -f neuromat_mean_field_lib.gawk \
    -v kappa_min="${kappa_min}" \
    -v kappa_max="${kappa_max}" \
    -v lambda_min="${lambda_min}" \
    -v lambda_max="${lambda_max}" \
    -v region="${region}" \
    -v dense="${dense}" \
    -v dead="${dead}" \
    > out/${name}.txt
  plot_mean_phases.sh ${kappa_min} ${kappa_max} ${lambda_min} ${lambda_max} ${name}