#!/bin/bash

CORES=16;
export EXEC=examples/bin/fmm_cheb

# Set run parameters
declare -a    nodes=();
declare -a    cores=();
declare -a mpi_proc=();
declare -a  threads=();
declare -a testcase=();
declare -a    n_pts=();
declare -a        m=();
declare -a        q=();
declare -a      tol=();
declare -a    depth=();
declare -a     unif=();
declare -a     adap=();
declare -a max_time=();


###################################################################################################
#                 NON-UNIFORM OCTREE, LAPLACE KERNEL, STRONG SCALABILITY                          #
###################################################################################################

# m=10, q=14, octants=
nodes+=(            1         8        64       512      4096     32768 :)
cores+=(     ${CORES}  ${CORES}  ${CORES}  ${CORES}  ${CORES}  ${CORES} :)
mpi_proc+=(         4        32       256      2048     16384     16384 :)
threads+=(   ${CORES}  ${CORES}  ${CORES}  ${CORES}  ${CORES}  ${CORES} :)
testcase+=(         1         1         1         1         1         1 :)
n_pts+=(    $((8**7)) $((8**7)) $((8**7)) $((8**7)) $((8**7)) $((8**7)) :)
m_pts+=(            1         1         1         1         1         1 :)
m+=(               10        10        10        10        10        10 :)
q+=(               14        14        14        14        14        14 :)
tol+=(           1e-0      1e-0      1e-0      1e-0      1e-0      1e-0 :)
depth+=(           15        15        15        15        15        15 :)
unif+=(             0         0         0         0         0         0 :)
adap+=(             0         0         0         0         0         0 :)
max_time+=(       800       800       800       800       800       800 :)


###################################################################################################
#                 NON-UNIFORM OCTREE, STOKES KERNEL, STRONG SCALABILITY                           #
###################################################################################################

# m=10, q=14, octants=
nodes+=(            1         8        64       512      4096     32768 :)
cores+=(     ${CORES}  ${CORES}  ${CORES}  ${CORES}  ${CORES}  ${CORES} :)
mpi_proc+=(         1         8        64       512      4096     32768 :)
threads+=(   ${CORES}  ${CORES}  ${CORES}  ${CORES}  ${CORES}  ${CORES} :)
testcase+=(         3         3         3         3         3         3 :)
n_pts+=(    $((8**7)) $((8**7)) $((8**7)) $((8**7)) $((8**7)) $((8**7)) :)
m_pts+=(            1         1         1         1         1         1 :)
m+=(               10        10        10        10        10        10 :)
q+=(               14        14        14        14        14        14 :)
tol+=(           1e-0      1e-0      1e-0      1e-0      1e-0      1e-0 :)
depth+=(           15        15        15        15        15        15 :)
unif+=(             0         0         0         0         0         0 :)
adap+=(             0         0         0         0         0         0 :)
max_time+=(      2400      2400      2400      2400      2400      2400 :)


###################################################################################################
#                  UNIFORM OCTREE, HELMHOLTZ KERNEL, STRONG SCALABILITY                           #
###################################################################################################

# m=10, q=14, octants=262144
nodes+=(            4         8        16        32        64        128        256        512       1024 :)
cores+=(     ${CORES}  ${CORES}  ${CORES}  ${CORES}  ${CORES}   ${CORES}   ${CORES}   ${CORES}   ${CORES} :)
mpi_proc+=(         4         8        16        32        64        128        256        512       1024 :)
threads+=(   ${CORES}  ${CORES}  ${CORES}  ${CORES}  ${CORES}   ${CORES}   ${CORES}   ${CORES}   ${CORES} :)
testcase+=(         5         5         5         5         5          5          5          5          5 :)
n_pts+=(    $((8**6)) $((8**6)) $((8**6)) $((8**6)) $((8**6))  $((8**6))  $((8**6))  $((8**6))  $((8**6)) :)
m_pts+=(            1         1         1         1         1          1          1          1          1 :)
m+=(               10        10        10        10        10         10         10         10         10 :)
q+=(               14        14        14        14        14         14         14         14         14 :)
tol+=(           1e-0      1e-0      1e-0      1e-0      1e-0       1e-0       1e-0       1e-0       1e-0 :)
depth+=(           15        15        15        15        15         15         15         15         15 :)
unif+=(             1         1         1         1         1          1          1          1          1 :)
adap+=(             0         0         0         0         0          0          0          0          0 :)
max_time+=(      2400      2400      2400      2400      2400       2400       2400       2400       2400 :)


###################################################################################################

# Export arrays
export    nodes_="$(declare -p    nodes)";
export    cores_="$(declare -p    cores)";
export mpi_proc_="$(declare -p mpi_proc)";
export  threads_="$(declare -p  threads)";
export testcase_="$(declare -p testcase)";
export    n_pts_="$(declare -p    n_pts)";
export    m_pts_="$(declare -p    m_pts)";
export        m_="$(declare -p        m)";
export        q_="$(declare -p        q)";
export      tol_="$(declare -p      tol)";
export    depth_="$(declare -p    depth)";
export     unif_="$(declare -p     unif)";
export     adap_="$(declare -p     adap)";
export max_time_="$(declare -p max_time)";

export RESULT_FNAME=$(basename ${0%.*}).out;
export WORK_DIR=$(dirname ${PWD}/$0)/..
cd ${WORK_DIR}

TERM_WIDTH=$(stty size | cut -d ' ' -f 2)
./scripts/.submit_jobs.sh | cut -b -${TERM_WIDTH}