1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
|
#!/bin/bash
if [ $# != 2 ]; then
echo "Usage: ./make_galrep.sh <range> <src_dir>"
exit 1
else
r=$1
d=$2
fi
echo 'running Magma galrep script on '${r} in ${d}
pushd ${HOME}/galrep
# Here we need the complete label then the a-invariant list, ":"-separated
# - in curvedata files this is fields 1,11
# - in allcurves files (as output by process_raw_curves): 1+2+4,6
awk '{print $1$2$4":"$6;}' ${d}/allcurves/allcurves.${r} > temp.${r}
#awk '{print $1":"$11;}' ${d}/curvedata/curvedata.${r} > temp.${r}
echo "ComputeQGaloisImages(\"temp.$r\", \"galrep.$r\"); quit;" > magma_temp.m
cat magma_temp.m
magma $HOME/galrep/nfgalrep.m magma_temp.m > /dev/null
mv galrep.${r} ${d}/galrep/galrep.${r}.copy
/bin/rm temp.${r} magma_temp.m
popd
|
