96 |
|
|
97 |
testoutput_var() |
testoutput_var() |
98 |
{ |
{ |
99 |
# testoutput_var dir s1 label subdir extension |
# testoutput_var dir s1 label subdir reference_output |
100 |
# |
# |
101 |
# compares 1 variable output selected from file $dir/$subdir/$OUTPUTFILE |
# compares 1 variable output selected from file $dir/$subdir/$OUTPUTFILE |
102 |
# with same output from reference file $dir/results/output.$extension |
# with same output from reference file $dir/results/$reference_output |
103 |
# using search strings s1 and text label |
# using search strings s1 and text label |
104 |
|
|
105 |
if [ $debug -gt 0 ]; then |
if [ $debug -gt 0 ]; then |
108 |
if [ -r $1/$4/$OUTPUTFILE ]; then |
if [ -r $1/$4/$OUTPUTFILE ]; then |
109 |
grep "$2" $1/$4/$OUTPUTFILE | sed 's/.*=//' | cat -n > tmp1.txt |
grep "$2" $1/$4/$OUTPUTFILE | sed 's/.*=//' | cat -n > tmp1.txt |
110 |
lncntA=`wc -l tmp1.txt | awk '{print $1}' ` |
lncntA=`wc -l tmp1.txt | awk '{print $1}' ` |
111 |
if [ $lncntA -lt 3 ]; then |
if [ $lncntA -lt 2 ]; then |
112 |
if [ $verbose -gt 0 ]; then |
if [ $verbose -gt 0 ]; then |
113 |
echo Not enough lines of output when searching for "$2" 1>&2 |
echo Not enough lines of output when searching for "$2" 1>&2 |
114 |
fi |
fi |
119 |
return 99 |
return 99 |
120 |
fi |
fi |
121 |
if [ $debug -gt 0 ]; then |
if [ $debug -gt 0 ]; then |
122 |
echo testoutput_var: grep "$2" $1/results/output.$5 1>&2 |
echo testoutput_var: grep "$2" $1/results/$5 1>&2 |
123 |
fi |
fi |
124 |
grep "$2" $1/results/output.$5 | sed 's/.*=//' | cat -n > tmp2.txt |
grep "$2" $1/results/$5 | sed 's/.*=//' | cat -n > tmp2.txt |
125 |
lncntB=`wc -l tmp2.txt | awk '{print $1}' ` |
lncntB=`wc -l tmp2.txt | awk '{print $1}' ` |
126 |
if [ $lncntB -lt 3 ]; then |
if [ $lncntB -lt 2 ]; then |
127 |
if [ $verbose -gt 0 ]; then |
if [ $verbose -gt 0 ]; then |
128 |
echo Not enough lines of output when searching for "$2" 1>&2 |
echo Not enough lines of output when searching for "$2" 1>&2 |
129 |
fi |
fi |
161 |
digits_of_similarity=`./tr_cmpnum < tmp4.txt` |
digits_of_similarity=`./tr_cmpnum < tmp4.txt` |
162 |
if [ $digits_of_similarity -eq 99 ]; then |
if [ $digits_of_similarity -eq 99 ]; then |
163 |
if [ $verbose -gt 0 ]; then |
if [ $verbose -gt 0 ]; then |
164 |
echo testoutput_var: No comparison was available for \"$2\" 1>&2 |
echo testoutput_var: No comparison was available for \"$3\" 1>&2 |
165 |
fi |
fi |
166 |
digits_of_similarity=99 |
digits_of_similarity=99 |
167 |
else |
else |
168 |
if [ $verbose -gt 0 ]; then |
if [ $verbose -gt 0 ]; then |
169 |
echo There were $digits_of_similarity decimal places of similarity for \"$2\" 1>&2 |
echo There were $digits_of_similarity decimal places of similarity for \"$3\" 1>&2 |
170 |
fi |
fi |
171 |
fi |
fi |
172 |
rm -f tmp1.txt tmp2.txt tmp3.txt tmp4.txt |
rm -f tmp1.txt tmp2.txt tmp3.txt tmp4.txt |
189 |
done |
done |
190 |
} |
} |
191 |
|
|
|
testoutput_ad() |
|
|
{ |
|
|
grep $3 $1/results_ad/output_adm.txt | awk '{print NR " " $5}' > t05.txt |
|
|
grep $3 $1/$2/$OUTPUTFILE | awk '{print NR " " $5}' > t15.txt |
|
|
grep $3 $1/results_ad/output_adm.txt | awk '{print NR " " $6}' > t06.txt |
|
|
grep $3 $1/$2/$OUTPUTFILE | awk '{print NR " " $6}' > t16.txt |
|
|
join t05.txt t15.txt > t5.txt |
|
|
join t06.txt t16.txt > t6.txt |
|
|
echo "-1" >> t5.txt |
|
|
echo "-1" >> t6.txt |
|
|
digits_5=`./tr_cmpnum < t5.txt` |
|
|
digits_6=`./tr_cmpnum < t6.txt` |
|
|
dashnum $digits_5 $digits_6 |
|
|
rm -f t[01][56].txt t[56].txt |
|
|
} |
|
|
|
|
192 |
check_for_add_mon_output() |
check_for_add_mon_output() |
193 |
{ |
{ |
194 |
# Check for additional types of monitor output |
# Check for additional types of monitor output |
216 |
|
|
217 |
testoutput_run() |
testoutput_run() |
218 |
{ |
{ |
219 |
# testoutput_run directory subdir extension |
# testoutput_run directory subdir reference_output |
220 |
# |
# |
221 |
# test output from 1 run in "directory" |
# test output from 1 run in "directory" |
222 |
if test "x$ADM" = x ; then |
# --> same processing for adjoint & forward test |
223 |
# default list of output variables to be checked: |
# default list of output variables to be checked: |
224 |
# 1rst : main variable used to decide if it pass or FAIL |
# 1rst : main variable used to decide if it pass or FAIL |
225 |
# others : number of matching digits to be printed in summary.txt |
# others : number of matching digits to be printed in summary.txt |
226 |
listChk=$DEF_CHECK_LIST |
listChk=$DEF_CHECK_LIST |
|
nbDef=`echo $listChk | awk '{print NF-1+3*gsub("+","E")}'` |
|
227 |
# load experiment-specific list from file "tr_checklist" (if it exist) |
# load experiment-specific list from file "tr_checklist" (if it exist) |
228 |
if test -r $1/$2/tr_checklist ; then listChk=`cat $1/$2/tr_checklist` ; fi |
if test -r $1/$2/tr_checklist ; then listChk=`cat $1/$2/tr_checklist` ; fi |
229 |
sVar=`echo $listChk | awk '{print $1}'` |
sVar=`echo $listChk | awk '{print $1}'` |
230 |
# remove 1rst var and expand the list: + => min max mean s.d |
# remove 1rst var and expand the list: + => min max mean s.d |
231 |
listVar=`echo $listChk | awk '{ for(i=2;i<=NF;i++){if (sub("+","",$i)) \ |
listVar=`echo $listChk | sed 's/ [a-zA-Z0-9]*+/&mn &mx &av &sd/g' \ |
232 |
{printf " %s %s %s %s",$i"mn",$i"mx",$i"av",$i"sd"} else {printf " %s",$i} }}'` |
| sed 's/+//g' | sed "s/^$sVar//"` |
233 |
if [ $debug -gt 0 ]; then echo "testoutput_run: listVar(I)='$listVar'" 1>&2 ; fi |
if [ $debug -gt 0 ]; then echo "testoutput_run: listVar(I)='$listVar'" 1>&2 ; fi |
234 |
for ii in 1 2 3 4 5 6 7 8 9 ; do |
for ii in 1 2 3 4 5 6 7 8 9 ; do |
235 |
tst=`eval 'echo "$HAVE_PTR0'$ii'"'` |
tst=`eval 'echo "$HAVE_PTR0'$ii'"'` |
252 |
case $xx in |
case $xx in |
253 |
'PS') if [ $debug -gt 0 ] |
'PS') if [ $debug -gt 0 ] |
254 |
then echo testoutput_run: testoutput_var $1 cg2d_init_res 1>&2 ; fi |
then echo testoutput_run: testoutput_var $1 cg2d_init_res 1>&2 ; fi |
255 |
testoutput_var $1 "cg2d_init_res" "cg2d init. residual" $2 $3; yy=$? |
testoutput_var $1 "cg2d_init_res" "Press. Solver (cg2d)" $2 $3; yy=$? |
256 |
if [ $debug -gt 0 ] ; then echo testoutput_run: cg2dres=$yy 1>&2 ; fi ;; |
if [ $debug -gt 0 ] ; then echo testoutput_run: cg2dres=$yy 1>&2 ; fi ;; |
257 |
'Tmn') testoutput_var $1 "dynstat_theta_min" "theta minimum" $2 $3; yy=$? ;; |
'Cost') testoutput_var $1 "ADM precision_derivative_cost" "ADM Cost" $2 $3; yy=$? ;; |
258 |
'Tmx') testoutput_var $1 "dynstat_theta_max" "theta maximum" $2 $3; yy=$? ;; |
'Grad') testoutput_var $1 "ADM precision_derivative_grad" "ADM Grad" $2 $3; yy=$? ;; |
259 |
'Tav') testoutput_var $1 "dynstat_theta_mean" "theta mean" $2 $3; yy=$? ;; |
'Tmn') testoutput_var $1 "dynstat_theta_min" "Theta minimum" $2 $3; yy=$? ;; |
260 |
'Tsd') testoutput_var $1 "dynstat_theta_sd" "theta s.d." $2 $3; yy=$? ;; |
'Tmx') testoutput_var $1 "dynstat_theta_max" "Theta maximum" $2 $3; yy=$? ;; |
261 |
'Smn') testoutput_var $1 "dynstat_salt_min" "salt minimum" $2 $3; yy=$? ;; |
'Tav') testoutput_var $1 "dynstat_theta_mean" "Theta mean" $2 $3; yy=$? ;; |
262 |
'Smx') testoutput_var $1 "dynstat_salt_max" "salt maximum" $2 $3; yy=$? ;; |
'Tsd') testoutput_var $1 "dynstat_theta_sd" "Theta Std.Dev" $2 $3; yy=$? ;; |
263 |
'Sav') testoutput_var $1 "dynstat_salt_mean" "salt mean" $2 $3; yy=$? ;; |
'Smn') testoutput_var $1 "dynstat_salt_min" "Salt minimum" $2 $3; yy=$? ;; |
264 |
'Ssd') testoutput_var $1 "dynstat_salt_sd" "salt s.d." $2 $3; yy=$? ;; |
'Smx') testoutput_var $1 "dynstat_salt_max" "Salt maximum" $2 $3; yy=$? ;; |
265 |
|
'Sav') testoutput_var $1 "dynstat_salt_mean" "Salt mean" $2 $3; yy=$? ;; |
266 |
|
'Ssd') testoutput_var $1 "dynstat_salt_sd" "Salt Std.Dev" $2 $3; yy=$? ;; |
267 |
'Umn') testoutput_var $1 "dynstat_uvel_min" "U minimum" $2 $3; yy=$? ;; |
'Umn') testoutput_var $1 "dynstat_uvel_min" "U minimum" $2 $3; yy=$? ;; |
268 |
'Umx') testoutput_var $1 "dynstat_uvel_max" "U maximum" $2 $3; yy=$? ;; |
'Umx') testoutput_var $1 "dynstat_uvel_max" "U maximum" $2 $3; yy=$? ;; |
269 |
'Uav') testoutput_var $1 "dynstat_uvel_mean" "U mean" $2 $3; yy=$? ;; |
'Uav') testoutput_var $1 "dynstat_uvel_mean" "U mean" $2 $3; yy=$? ;; |
270 |
'Usd') testoutput_var $1 "dynstat_uvel_sd" "U s.d." $2 $3; yy=$? ;; |
'Usd') testoutput_var $1 "dynstat_uvel_sd" "U Std.Dev" $2 $3; yy=$? ;; |
271 |
'Vmn') testoutput_var $1 "dynstat_vvel_min" "V minimum" $2 $3; yy=$? ;; |
'Vmn') testoutput_var $1 "dynstat_vvel_min" "V minimum" $2 $3; yy=$? ;; |
272 |
'Vmx') testoutput_var $1 "dynstat_vvel_max" "V maximum" $2 $3; yy=$? ;; |
'Vmx') testoutput_var $1 "dynstat_vvel_max" "V maximum" $2 $3; yy=$? ;; |
273 |
'Vav') testoutput_var $1 "dynstat_vvel_mean" "V mean" $2 $3; yy=$? ;; |
'Vav') testoutput_var $1 "dynstat_vvel_mean" "V mean" $2 $3; yy=$? ;; |
274 |
'Vsd') testoutput_var $1 "dynstat_vvel_sd" "V s.d." $2 $3; yy=$? ;; |
'Vsd') testoutput_var $1 "dynstat_vvel_sd" "V Std.Dev" $2 $3; yy=$? ;; |
275 |
'pt1mn'|'pt2mn'|'pt3mn'|'pt4mn'|'pt5mn') ii=`echo $xx | sed 's/pt//' | sed 's/..$//'` |
'pt1mn'|'pt2mn'|'pt3mn'|'pt4mn'|'pt5mn') ii=`echo $xx | sed 's/pt//' | sed 's/..$//'` |
276 |
testoutput_var $1 "trcstat_ptracer0"$ii"_min" "p0"$ii"_min" $2 $3; yy=$? ;; |
testoutput_var $1 "trcstat_ptracer0"$ii"_min" "p0"$ii"_min" $2 $3; yy=$? ;; |
277 |
'pt1mx'|'pt2mx'|'pt3mx'|'pt4mx'|'pt5mx') ii=`echo $xx | sed 's/pt//' | sed 's/..$//'` |
'pt1mx'|'pt2mx'|'pt3mx'|'pt4mx'|'pt5mx') ii=`echo $xx | sed 's/pt//' | sed 's/..$//'` |
279 |
'pt1av'|'pt2av'|'pt3av'|'pt4av'|'pt5av') ii=`echo $xx | sed 's/pt//' | sed 's/..$//'` |
'pt1av'|'pt2av'|'pt3av'|'pt4av'|'pt5av') ii=`echo $xx | sed 's/pt//' | sed 's/..$//'` |
280 |
testoutput_var $1 "trcstat_ptracer0"$ii"_mean" "p0"$ii"_mean" $2 $3; yy=$? ;; |
testoutput_var $1 "trcstat_ptracer0"$ii"_mean" "p0"$ii"_mean" $2 $3; yy=$? ;; |
281 |
'pt1sd'|'pt2sd'|'pt3sd'|'pt4sd'|'pt5sd') ii=`echo $xx | sed 's/pt//' | sed 's/..$//'` |
'pt1sd'|'pt2sd'|'pt3sd'|'pt4sd'|'pt5sd') ii=`echo $xx | sed 's/pt//' | sed 's/..$//'` |
282 |
testoutput_var $1 "trcstat_ptracer0"$ii"_sd" "p0"$ii"_sd" $2 $3; yy=$? ;; |
testoutput_var $1 "trcstat_ptracer0"$ii"_sd" "p0"$ii"_StDv" $2 $3; yy=$? ;; |
283 |
|
'Qntmn') testoutput_var $1 "extforcing_qnet_min" "Qnet minimum" $2 $3; yy=$? ;; |
284 |
|
'Qntmx') testoutput_var $1 "extforcing_qnet_max" "Qnet maximum" $2 $3; yy=$? ;; |
285 |
|
'Qntav') testoutput_var $1 "extforcing_qnet_mean" "Qnet mean" $2 $3; yy=$? ;; |
286 |
|
'Qntsd') testoutput_var $1 "extforcing_qnet_sd" "Qnet Std.Dev" $2 $3; yy=$? ;; |
287 |
|
'aSImn') testoutput_var $1 "seaice_area_min" "SIce Area min" $2 $3; yy=$? ;; |
288 |
|
'aSImx') testoutput_var $1 "seaice_area_max" "SIce Area max" $2 $3; yy=$? ;; |
289 |
|
'aSIav') testoutput_var $1 "seaice_area_mean" "SIce Area mean" $2 $3; yy=$? ;; |
290 |
|
'aSIsd') testoutput_var $1 "seaice_area_sd" "SIce Area StDv" $2 $3; yy=$? ;; |
291 |
|
'hSImn') testoutput_var $1 "seaice_heff_min" "SIce Heff min" $2 $3; yy=$? ;; |
292 |
|
'hSImx') testoutput_var $1 "seaice_heff_max" "SIce Hell max" $2 $3; yy=$? ;; |
293 |
|
'hSIav') testoutput_var $1 "seaice_heff_mean" "SIce Hell mean" $2 $3; yy=$? ;; |
294 |
|
'hSIsd') testoutput_var $1 "seaice_heff_sd" "SIce Hell StDv" $2 $3; yy=$? ;; |
295 |
|
'AthSiG') testoutput_var $1 "thSI_Ice_Area_G" "thSIc Area Global" $2 $3; yy=$? ;; |
296 |
|
'AthSiS') testoutput_var $1 "thSI_Ice_Area_S" "thSIc Area South" $2 $3; yy=$? ;; |
297 |
|
'AthSiN') testoutput_var $1 "thSI_Ice_Area_N" "thSIc Area North" $2 $3; yy=$? ;; |
298 |
|
'HthSiG') testoutput_var $1 "thSI_IceH_ave_G" "thSIc H Global" $2 $3; yy=$? ;; |
299 |
|
'HthSiS') testoutput_var $1 "thSI_IceH_ave_S" "thSIc H South" $2 $3; yy=$? ;; |
300 |
|
'HthSiN') testoutput_var $1 "thSI_IceH_ave_N" "thSIc H North" $2 $3; yy=$? ;; |
301 |
*) yy=99; echo "WARNING: asking for var=$xx : not recognized !" 1>&2 ;; |
*) yy=99; echo "WARNING: asking for var=$xx : not recognized !" 1>&2 ;; |
302 |
esac |
esac |
303 |
if test $xx = $sVar |
if test $xx = $sVar |
307 |
done |
done |
308 |
|
|
309 |
nbVar=`echo $listVar | awk '{print NF}'` |
nbVar=`echo $listVar | awk '{print NF}'` |
310 |
if [ $nbVar -lt $nbDef ] ; then |
if [ $nbVar -lt $LEN_CHECK_LIST ] ; then |
311 |
#-- fill line (up to standard length) with dot: |
#-- fill line (up to standard length) with dot: |
312 |
adNul=`expr $nbDef - $nbVar | awk '{for(i=1;i<=$1;i++){print "."}}'` |
adNul=`expr $LEN_CHECK_LIST - $nbVar | awk '{for(i=1;i<=$1;i++){print "."}}'` |
313 |
echo $allargs $adNul |
echo $allargs $adNul |
314 |
else |
else |
315 |
echo $allargs |
echo $allargs |
316 |
fi |
fi |
317 |
|
# <-- same processing for adjoint & forward test |
|
else |
|
|
testoutput_ad $1 $2 "precision_grdchk_result" |
|
|
fi |
|
318 |
} |
} |
319 |
|
|
320 |
genmakemodel() |
genmakemodel() |
586 |
RETVAL=8 |
RETVAL=8 |
587 |
ENDVAL=-1 |
ENDVAL=-1 |
588 |
else |
else |
589 |
if test $OUTPUTFILE -ot $EXECUTABLE ; then |
if test ! -f $OUTPUTFILE -o $OUTPUTFILE -ot $EXECUTABLE ; then |
590 |
|
# output do not exist or is older than executable: |
591 |
rm -f run.log ; touch run.log |
rm -f run.log ; touch run.log |
592 |
if test -f run.log_00 ; then cat run.log_00 >> run.log ; fi |
if test -f run.log_00 ; then cat run.log_00 >> run.log ; fi |
593 |
( eval $COMMAND ) >> run.log 2>&1 |
( eval $COMMAND ) >> run.log 2>&1 |
638 |
double a,b,abave,relerr; |
double a,b,abave,relerr; |
639 |
best = -22; |
best = -22; |
640 |
lncnt = 0; |
lncnt = 0; |
641 |
while( 1 & (lncnt+=1) < 999 ) { |
while( 1 & ( (lncnt+=1) < 999 ) ) { |
642 |
scanf("%d", &linnum); |
scanf("%d", &linnum); |
643 |
if (linnum == -1) break; |
if (linnum == -1) break; |
644 |
scanf("%lf", &a); scanf("%lf", &b); |
scanf("%lf", &a); scanf("%lf", &b); |
685 |
listPrt=$@ |
listPrt=$@ |
686 |
listRes=`echo $listPrt | sed 's/>//' | sed 's/<//'` |
listRes=`echo $listPrt | sed 's/>//' | sed 's/<//'` |
687 |
xx=`echo $listPrt | sed 's/.*>//' | sed 's/<.*//' | awk '{print $1}'` |
xx=`echo $listPrt | sed 's/.*>//' | sed 's/<.*//' | awk '{print $1}'` |
688 |
|
printf '%3s' $listPrt |
689 |
|
# line below does not work on hp-ux_ia64 : do those substitutions later on |
690 |
# printf '%3s' $listPrt | sed 's/ 99/ --/g' | sed 's/ > />/' | sed 's/ < /</' |
# printf '%3s' $listPrt | sed 's/ 99/ --/g' | sed 's/ > />/' | sed 's/ < /</' |
|
printf '%3s' $listPrt | sed 's/ 99/ --/g' | sed 's/ > / /' | sed 's/ < / /' |
|
691 |
|
|
692 |
if [ $xx = '..' ]; then |
if [ $xx = '..' ]; then |
693 |
printf ' N/O ' |
printf ' N/O ' |
711 |
{ |
{ |
712 |
if [ $# -eq 1 ]; then |
if [ $# -eq 1 ]; then |
713 |
for arg in * ; do |
for arg in * ; do |
714 |
test -d $arg/$1 && echo $arg |
test -f $arg/$1 && echo $arg |
715 |
done |
done |
716 |
else |
else |
717 |
echo $* |
echo $* |
896 |
NODEPEND=t |
NODEPEND=t |
897 |
fi |
fi |
898 |
|
|
899 |
|
#- setting for forward or ADM testing |
900 |
|
if test "x$ADM" = xt ; then |
901 |
|
code_dir=code_ad |
902 |
|
input_dirs='input_ad input' |
903 |
|
ref_outp="output_adm.txt" |
904 |
|
EXECUTABLE="mitgcmuv_ad" |
905 |
|
else |
906 |
|
code_dir=code |
907 |
|
input_dirs='input' |
908 |
|
ref_outp="output.txt" |
909 |
|
EXECUTABLE="mitgcmuv" |
910 |
|
fi |
911 |
|
|
912 |
if test "x$TESTDIRS" = x ; then |
if test "x$TESTDIRS" = x ; then |
913 |
if test "x$ADM" = xt ; then |
LIST=`scandirs results/$ref_outp` |
|
LIST=`scandirs results_ad` |
|
|
else |
|
|
LIST=`scandirs results` |
|
|
fi |
|
914 |
else |
else |
915 |
#- expand group of experiments: |
#- expand group of experiments: |
916 |
LIST=" " |
LIST=" " |
954 |
OPTFILE=$MITGCM_OF |
OPTFILE=$MITGCM_OF |
955 |
fi |
fi |
956 |
|
|
957 |
if test "x$ADM" = xt ; then |
OUTPUTFILE=$ref_outp |
|
EXECUTABLE="mitgcmuv_ad" |
|
|
OUTPUTFILE="output_adm.txt" |
|
|
else |
|
|
EXECUTABLE="mitgcmuv" |
|
|
OUTPUTFILE="output.txt" |
|
|
fi |
|
|
|
|
958 |
if test "x$COMMAND" = x ; then |
if test "x$COMMAND" = x ; then |
959 |
COMMAND="./$EXECUTABLE > $OUTPUTFILE" |
COMMAND="./$EXECUTABLE > $OUTPUTFILE" |
960 |
fi |
fi |
962 |
OUTPUTFILE="STDOUT.0000" |
OUTPUTFILE="STDOUT.0000" |
963 |
fi |
fi |
964 |
|
|
|
#echo "OK" |
|
965 |
echo "OK (COMMAND= $COMMAND )" |
echo "OK (COMMAND= $COMMAND )" |
966 |
|
|
967 |
# set the Default List of output variables to be checked: |
# set the Default List of output variables to be checked: |
968 |
# (use default or load experiment-specific list from file "tr_checklist") |
# (use default or load experiment-specific list from file "tr_checklist") |
969 |
# content : 1rst = main variable used to decide if it pass or FAIL |
# content : 1rst = main variable used to decide if it pass or FAIL |
970 |
# others = number of matching digits to be printed in summary.txt |
# others = number of matching digits to be printed in summary.txt |
971 |
if test "x$ADM" != xt ; then |
if test "x$ADM" = x ; then |
972 |
DEF_CHECK_LIST='PS PS T+ S+ U+ V+ pt1+ pt2+ pt3+ pt4+ pt5+' |
DEF_CHECK_LIST='PS PS T+ S+ U+ V+ pt1+ pt2+ pt3+ pt4+ pt5+' |
973 |
EMPTY_RESULTS='.. .. .. .. .. .. .. .. .. .. .. .. .. .. .. .. ..' |
EMPTY_RESULTS='.. .. .. .. .. .. .. .. .. .. .. .. .. .. .. .. ..' |
974 |
nbDef=`echo $DEF_CHECK_LIST | awk '{print NF-1+3*gsub("+","E")}'` |
LEN_CHECK_LIST=`echo $DEF_CHECK_LIST | sed 's/ [a-zA-Z0-9]*+/&mn &mx &av &sd/g' | awk '{print NF-1}'` |
975 |
ii=`echo $EMPTY_RESULTS | awk '{print NF}'` |
ii=`echo $EMPTY_RESULTS | awk '{print NF}'` |
976 |
EMPTY_RESULTS=$EMPTY_RESULTS`expr $nbDef - $ii | awk 'BEGIN{FS=":"}{for(i=1;i<=$1;i++){printf " ."}}'` |
EMPTY_RESULTS=$EMPTY_RESULTS`expr $LEN_CHECK_LIST - $ii | awk 'BEGIN{FS=":"}{for(i=1;i<=$1;i++){printf " ."}}'` |
977 |
|
else |
978 |
|
DEF_CHECK_LIST='Grad Cost Grad' |
979 |
|
EMPTY_RESULTS='.. ..' |
980 |
|
LEN_CHECK_LIST=`echo $DEF_CHECK_LIST | sed 's/ [a-zA-Z0-9]*+/&mn &mx &av &sd/g' | awk '{print NF-1}'` |
981 |
fi |
fi |
982 |
|
|
983 |
# create the FORTRAN comparison code |
# create the FORTRAN comparison code |
1071 |
else |
else |
1072 |
echo "ADJOINT=true" >> $SUMMARY |
echo "ADJOINT=true" >> $SUMMARY |
1073 |
echo >> $SUMMARY |
echo >> $SUMMARY |
1074 |
|
if [ $MATCH_CRIT -lt 10 ] ; |
1075 |
|
then line_0="default "$MATCH_CRIT ; |
1076 |
|
else line_0="default "$MATCH_CRIT ; fi |
1077 |
|
echo "$line_0" | tee -a $SUMMARY |
1078 |
cat << EOF | tee -a $SUMMARY |
cat << EOF | tee -a $SUMMARY |
1079 |
G D M C G |
G D M C G |
1080 |
E p a R o r |
E p a R o r |
1108 |
# results in the correct location--or skip this directory! |
# results in the correct location--or skip this directory! |
1109 |
fout= |
fout= |
1110 |
if test "x$ADM" = x ; then |
if test "x$ADM" = x ; then |
1111 |
fout=$dir"/results/output.txt" |
fout=$dir"/results/"$ref_outp |
1112 |
else |
else |
1113 |
fout=$dir"/results_ad/output_adm.txt" |
fout=$dir"/results_ad/"$ref_outp |
1114 |
fi |
fi |
1115 |
if test ! -r $fout ; then |
if test ! -r $fout ; then |
1116 |
echo "can't read \"$fout\" -- skipping $dir" |
echo "can't read \"$fout\" -- skipping $dir" |
1125 |
builddir="build" |
builddir="build" |
1126 |
if test ! -d $dir/$builddir ; then mkdir $dir/$builddir ; fi |
if test ! -d $dir/$builddir ; then mkdir $dir/$builddir ; fi |
1127 |
rundir="run" |
rundir="run" |
|
#rundir=$builddir |
|
1128 |
if test ! -d $dir/$rundir ; then |
if test ! -d $dir/$rundir ; then |
1129 |
rundir=$builddir |
rundir=$builddir |
1130 |
fi |
fi |
1131 |
|
CODE_DIR=$dir/$code_dir |
|
if test "x$ADM" = x ; then |
|
|
code_dir=code |
|
|
CODE_DIR=$dir/code |
|
|
input_dirs='input' |
|
|
else |
|
|
code_dir=code_ad |
|
|
CODE_DIR=$dir/code_ad |
|
|
input_dirs='input_ad input' |
|
|
fi |
|
1132 |
BUILD_DIR=$dir/$builddir |
BUILD_DIR=$dir/$builddir |
1133 |
|
|
1134 |
if test ! -r $CODE_DIR"/SIZE.h_mpi" -a "x$MPI" = "xt" ; then |
if test ! -r $CODE_DIR"/SIZE.h_mpi" -a "x$MPI" = "xt" ; then |
1190 |
&& makemodel $dir/$builddir && make=Y \ |
&& makemodel $dir/$builddir && make=Y \ |
1191 |
&& linkdata $dir/$rundir $input_dirs \ |
&& linkdata $dir/$rundir $input_dirs \ |
1192 |
&& runmodel $dir/$rundir && run=Y \ |
&& runmodel $dir/$rundir && run=Y \ |
1193 |
&& results=`testoutput_run $dir $rundir "txt"` |
&& results=`testoutput_run $dir $rundir $ref_outp` |
1194 |
fi |
fi |
1195 |
|
|
1196 |
echo |
echo |
1197 |
if test "x$ADM" = x ; then |
# --> same processing for adjoint & forward test |
1198 |
fres=`formatresults $dir ${genmake:-N} ${makedepend:-N} ${make:-N} ${run:-N} $results` |
fres=`formatresults $dir ${genmake:-N} ${makedepend:-N} ${make:-N} ${run:-N} $results` |
1199 |
echo |
echo |
1200 |
echo "$fres" >> $SUMMARY |
echo "$fres" | sed 's/ 99/ --/g' | sed 's/ > />/' | sed 's/ < /</' >> $SUMMARY |
1201 |
touch $CDIR"/summary.txt" |
touch $CDIR"/summary.txt" |
1202 |
echo "fresults='$fres'" >> $CDIR"/summary.txt" |
echo "fresults='$fres'" | sed 's/ 99/ --/g' >> $CDIR"/summary.txt" |
1203 |
echo "MACH='$MACH'" >> $CDIR"/summary.txt" |
echo "MACH='$MACH'" >> $CDIR"/summary.txt" |
1204 |
echo "UNAMEA='$UNAMEA'" >> $CDIR"/summary.txt" |
echo "UNAMEA='$UNAMEA'" >> $CDIR"/summary.txt" |
1205 |
echo "DATE='$DATE'" >> $CDIR"/summary.txt" |
echo "DATE='$DATE'" >> $CDIR"/summary.txt" |
1216 |
run_clean $dir/tr_run.$ex |
run_clean $dir/tr_run.$ex |
1217 |
linkdata $dir/tr_run.$ex input.$ex input |
linkdata $dir/tr_run.$ex input.$ex input |
1218 |
runmodel $dir/tr_run.$ex && run=Y \ |
runmodel $dir/tr_run.$ex && run=Y \ |
1219 |
&& results=`testoutput_run $dir tr_run.$ex ${ex}".txt"` |
&& results=`testoutput_run $dir tr_run.$ex "output.${ex}.txt"` |
|
fres=`printf '%s %s %s %s' ${genmake:-N} ${makedepend:-N} ${make:-N} ${run:-N}` |
|
1220 |
fres=`formatresults $dir ${genmake:-N} ${makedepend:-N} ${make:-N} ${run:-N} $results` |
fres=`formatresults $dir ${genmake:-N} ${makedepend:-N} ${make:-N} ${run:-N} $results` |
1221 |
fres="$fres.$ex" |
fres="$fres.$ex" |
1222 |
echo |
echo |
1223 |
echo "$fres" >> $SUMMARY |
echo "$fres" | sed 's/ 99/ --/g' | sed 's/ > />/' | sed 's/ < /</' >> $SUMMARY |
1224 |
touch $CDIR"/summary.txt" |
touch $CDIR"/summary.txt" |
1225 |
echo "fresults='$fres'" >> $CDIR"/summary.txt" |
echo "fresults='$fres'" | sed 's/ 99/ --/g' >> $CDIR"/summary.txt" |
1226 |
echo "MACH='$MACH'" >> $CDIR"/summary.txt" |
echo "MACH='$MACH'" >> $CDIR"/summary.txt" |
1227 |
echo "UNAMEA='$UNAMEA'" >> $CDIR"/summary.txt" |
echo "UNAMEA='$UNAMEA'" >> $CDIR"/summary.txt" |
1228 |
echo "DATE='$DATE'" >> $CDIR"/summary.txt" |
echo "DATE='$DATE'" >> $CDIR"/summary.txt" |
1229 |
echo "tdir='$dir.$ex'" >> $CDIR"/summary.txt" |
echo "tdir='$dir.$ex'" >> $CDIR"/summary.txt" |
1230 |
if test "x$POSTCLEAN" = xt ; then |
if test "x$POSTCLEAN" = xt ; then |
1231 |
run_clean $dir/tr_run.$ex |
run_clean $dir/tr_run.$ex |
1232 |
fi |
fi |
1233 |
done |
done |
1234 |
else |
# <-- same processing for adjoint & forward test |
1235 |
fres=`printf '%s %s %s %s' ${genmake:-N} ${makedepend:-N} ${make:-N} ${run:-N}` |
if test "x$ADM" = xt ; then |
1236 |
fres=$fres"$results $dir" |
grep -A3 'Seconds in section "ALL' $dir/$rundir/$OUTPUTFILE \ |
|
echo |
|
|
echo "$fres" >> $SUMMARY |
|
|
touch $CDIR"/summary.txt" |
|
|
echo "fresults='$fres'" >> $CDIR"/summary.txt" |
|
|
echo "MACH='$MACH'" >> $CDIR"/summary.txt" |
|
|
echo "UNAMEA='$UNAMEA'" >> $CDIR"/summary.txt" |
|
|
echo "DATE='$DATE'" >> $CDIR"/summary.txt" |
|
|
echo "tdir='$dir'" >> $CDIR"/summary.txt" |
|
|
grep -A3 'Seconds in section "ALL' $dir/$rundir/$OUTPUTFILE \ |
|
1237 |
>> $CDIR"/summary.txt" |
>> $CDIR"/summary.txt" |
1238 |
fi |
fi |
1239 |
|
|
1282 |
rm -f tr_cmpnum.c tr_cmpnum |
rm -f tr_cmpnum.c tr_cmpnum |
1283 |
|
|
1284 |
if test "x$CLEANUP" != xt ; then |
if test "x$CLEANUP" != xt ; then |
1285 |
cat $SUMMARY | sed 's/ \. \. \. \. \. \. \. \. \. \. \. \.//' |
cat $SUMMARY | sed 's/ \. \. \. \. \. \. \. \. \. \. \. \. //' |
1286 |
if test -e tr_out.txt ; then |
if test -e tr_out.txt ; then |
1287 |
mv tr_out.txt tr_out.txt.old |
mv tr_out.txt tr_out.txt.old |
1288 |
fi |
fi |
1289 |
cat $SUMMARY | sed '11,$ s/ \.//g' > tr_out.txt |
cat $SUMMARY | sed '/^[YN] [YN] [YN] [YN]/ s/ \. //g' > tr_out.txt |
1290 |
fi |
fi |
1291 |
|
|
1292 |
if test "x$DELDIR" = xt ; then |
if test "x$DELDIR" = xt ; then |