densities

#!/bin/tcsh -f
unlimit stack
setenv SCMAXHEAP 1700
rm -f common-adifor.s common-tapenade.s g_common-adifor.s g_particle-adifor1.s
rm -f g_saddle-adifor1.s h_common-adifor.s h_g_common-adifor.s
rm -f h_g_particle-adifor1.s h_g_saddle-adifor1.s h_particle-adifor2.s
rm -f h_saddle-adifor1.s h_saddle-adifor2.s particle-adifor1000.s
rm -f particle-adifor1.s particle-adifor2.s particle-fadbad++.s
rm -f particle-tapenade1000.s particle-tapenade_gv.s particle-tapenade_hv.s
rm -f saddle-adifor1000.s saddle-adifor1.s saddle-adifor2.s
rm -f saddle-fadbad++.s saddle-tapenade1000.s saddle-tapenade_gv.s
rm -f saddle-tapenade_hv.s foo.s
rm -f saddle-stalingrad1000.{c,s} particle-stalingrad1000.{c,s}

g++ -S\
    -I ~/pkg/FADBAD++\
    -O3\
    -march=k8\
    -mfpmath=sse\
    -fomit-frame-pointer\
    -ffast-math\
    saddle-fadbad++.cpp
echo saddle FADBAD++
./density saddle-fadbad++

Adifor2.1 AD_SCRIPT=saddle-adifor1.adf
Adifor2.1 AD_SCRIPT=saddle-adifor2.adf
gfortran -S\
    -O3\
    -march=k8\
    -mfpmath=sse\
    -fomit-frame-pointer\
    -ffast-math\
    saddle-adifor1000.f saddle-adifor1.f saddle-adifor2.f common-adifor.f\
    g_saddle-adifor1.f h_common-adifor.f h_g_saddle-adifor1.f\
    h_saddle-adifor1.f h_saddle-adifor2.f
rm -f *.{A,aux} *~ .*.f
rm -rf AD_cache
rm -f g_saddle-adifor1.f h_common-adifor.f h_g_saddle-adifor1.f\
      h_saddle-adifor1.f h_saddle-adifor2.f
rm -f foo.s
cat saddle-adifor1000.s saddle-adifor1.s saddle-adifor2.s common-adifor.s\
    g_saddle-adifor1.s h_common-adifor.s h_g_saddle-adifor1.s\
    h_saddle-adifor1.s h_saddle-adifor2.s >foo.s
echo saddle ADIFOR
./density foo

rm -f DIFFSIZES.inc
tapenade -root inner\
         -d -multi -o saddle-tapenade\
         -diffvarname "_g" -difffuncname "_g"\
         common-tapenade.f saddle-tapenade.f
cp saddle-tapenade-DIFFSIZES1.inc DIFFSIZES.inc
tapenade -root outer\
         -d -multi -o saddle-tapenade\
         -diffvarname "_h" -difffuncname "_h"\
         common-tapenade.f saddle-tapenade.f saddle-tapenade_gv.f
cp saddle-tapenade-DIFFSIZES2.inc DIFFSIZES.inc
rm -f foo.f
sed -f saddle-tapenade.sed <saddle-tapenade_hv.f >foo.f
mv foo.f saddle-tapenade_hv.f
gfortran -S\
    -O3\
    -march=k8\
    -mfpmath=sse\
    -fomit-frame-pointer\
    -ffast-math\
    common-tapenade.f\
    saddle-tapenade1000.f saddle-tapenade_gv.f saddle-tapenade_hv.f
rm -f DIFFSIZES.inc saddle-tapenade_{g,h}v.{f,msg}
rm -f foo.s
cat common-tapenade.s\
    saddle-tapenade1000.s saddle-tapenade_gv.s saddle-tapenade_hv.s >foo.s
echo saddle Tapenade
./density foo

stalingrad -no-assert -compile -c -sra -inline -alias saddle-stalingrad1000
gcc -S\
    -O3\
    -march=k8\
    -mfpmath=sse\
    -fomit-frame-pointer\
    -ffast-math\
    saddle-stalingrad1000.c
echo saddle Stalingrad
./density saddle-stalingrad1000

g++ -S\
    -I ~/pkg/FADBAD++\
    -O3\
    -march=k8\
    -mfpmath=sse\
    -fomit-frame-pointer\
    -ffast-math\
    particle-fadbad++.cpp
echo particle FADBAD++
./density particle-fadbad++

Adifor2.1 AD_SCRIPT=particle-adifor1.adf
Adifor2.1 AD_SCRIPT=particle-adifor2.adf
gfortran -S\
    -O3\
    -march=k8\
    -mfpmath=sse\
    -fomit-frame-pointer\
    -ffast-math\
    particle-adifor1000.f particle-adifor1.f particle-adifor2.f\
    common-adifor.f\
    g_common-adifor.f h_common-adifor.f h_g_common-adifor.f\
    g_particle-adifor1.f h_g_particle-adifor1.f h_particle-adifor2.f
rm -f *.{A,aux} *~ .*.f
rm -rf AD_cache
rm -f g_common-adifor.f h_common-adifor.f h_g_common-adifor.f\
      g_particle-adifor1.f h_g_particle-adifor1.f h_particle-adifor2.f
rm -f foo.s
cat particle-adifor1000.s particle-adifor1.s particle-adifor2.s\
    common-adifor.s\
    g_common-adifor.s h_common-adifor.s h_g_common-adifor.s\
    g_particle-adifor1.s h_g_particle-adifor1.s h_particle-adifor2.s >foo.s
echo particle ADIFOR
./density foo

rm -f DIFFSIZES.inc
tapenade -root p\
         -d -multi -o particle-tapenade\
         -diffvarname "_g" -difffuncname "_g"\
         common-tapenade.f particle-tapenade.f
cp particle-tapenade-DIFFSIZES1.inc DIFFSIZES.inc
tapenade -root naive_euler\
         -d -multi -o particle-tapenade\
         -diffvarname "_h" -difffuncname "_h"\
         common-tapenade.f particle-tapenade.f particle-tapenade_gv.f
cp particle-tapenade-DIFFSIZES2.inc DIFFSIZES.inc
rm -f foo.f
sed -f particle-tapenade.sed <particle-tapenade_hv.f >foo.f
mv foo.f particle-tapenade_hv.f
gfortran -S\
    -O3\
    -march=k8\
    -mfpmath=sse\
    -fomit-frame-pointer\
    -ffast-math\
    common-tapenade.f\
    particle-tapenade1000.f particle-tapenade_gv.f particle-tapenade_hv.f
rm -f DIFFSIZES.inc particle-tapenade_{g,h}v.{f,msg}
rm -f foo.s
cat common-tapenade.s\
    particle-tapenade1000.s particle-tapenade_gv.s particle-tapenade_hv.s\
    >foo.s
echo particle Tapenade
./density foo

stalingrad -no-assert -compile -c -sra -inline -alias particle-stalingrad1000
gcc -S\
    -O3\
    -march=k8\
    -mfpmath=sse\
    -fomit-frame-pointer\
    -ffast-math\
    particle-stalingrad1000.c
echo particle Stalingrad
./density particle-stalingrad1000

rm -f common-adifor.s common-tapenade.s g_common-adifor.s g_particle-adifor1.s
rm -f g_saddle-adifor1.s h_common-adifor.s h_g_common-adifor.s
rm -f h_g_particle-adifor1.s h_g_saddle-adifor1.s h_particle-adifor2.s
rm -f h_saddle-adifor1.s h_saddle-adifor2.s particle-adifor1000.s
rm -f particle-adifor1.s particle-adifor2.s particle-fadbad++.s
rm -f particle-tapenade1000.s particle-tapenade_gv.s particle-tapenade_hv.s
rm -f saddle-adifor1000.s saddle-adifor1.s saddle-adifor2.s
rm -f saddle-fadbad++.s saddle-tapenade1000.s saddle-tapenade_gv.s
rm -f saddle-tapenade_hv.s foo.s
rm -f saddle-stalingrad1000.{c,s} particle-stalingrad1000.{c,s}