script_annulus_benchmark

#! /bin/bash
#   nelr        = int(sys.argv[1])
#   visu        = int(sys.argv[2])
#   nqperdim    = int(sys.argv[3])
#   mapping     = int(sys.argv[4])
#   xi          = int(sys.argv[5])
#   etablobstar = float(sys.argv[6])
#   rhoblobstar = float(sys.argv[7])
#   zblob       = float(sys.argv[8])
#   Rblob       = float(sys.argv[9])
#   etalithstar = float(sys.argv[10])

echo '#did you choose the right planet in flaps.py?'

nqperdim=3
xi=6

rm errors_map*
rm timings_*

for mapping in 2 3 4; do
    for l in 16 24 32 40 48 56 64 72 80 96 128; do
        echo nelr= $l - nqperdim= $nqperdim - mapping= $mapping - xi= $xi
        python3 flaps.py $l 1 $nqperdim $mapping $xi 1 1 1 1 1 > opla_${l}_${xi}
        grep errv opla_${l}_${xi} >> errors_mappingQ${mapping}
        grep sanity opla_${l}_${xi} >> timings_sanity_mappingQ${mapping}
        grep build opla_${l}_${xi} | grep FE >> timings_build_mappingQ${mapping}
        grep solving opla_${l}_${xi} >> timings_solving_mappingQ${mapping}
    done
done