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 |
if [ $debug -gt 0 ]; then |
# default list of output variables to be checked: |
224 |
echo testoutput_run: testoutput_var $1 cg2d_init_res 1>&2 |
# 1rst : main variable used to decide if it pass or FAIL |
225 |
fi |
# others : number of matching digits to be printed in summary.txt |
226 |
testoutput_var $1 "cg2d_init_res" "cg2d init. residual" $2 $3; cg2dres=$? |
listChk=$DEF_CHECK_LIST |
227 |
if [ $debug -gt 0 ]; then |
# load experiment-specific list from file "tr_checklist" (if it exist) |
228 |
echo testoutput_run: cg2dres=$cg2dres 1>&2 |
if test -r $1/$2/tr_checklist ; then listChk=`cat $1/$2/tr_checklist` ; fi |
229 |
fi |
sVar=`echo $listChk | awk '{print $1}'` |
230 |
testoutput_var $1 "dynstat_theta_min" "theta minimum" $2 $3; tmin=$? |
# remove 1rst var and expand the list: + => min max mean s.d |
231 |
testoutput_var $1 "dynstat_theta_max" "theta maximum" $2 $3; tmax=$? |
listVar=`echo $listChk | sed 's/ [a-zA-Z0-9]*+/&mn &mx &av &sd/g' \ |
232 |
testoutput_var $1 "dynstat_theta_mean" "theta mean" $2 $3; tmean=$? |
| sed 's/+//g' | sed "s/^$sVar//"` |
233 |
testoutput_var $1 "dynstat_theta_sd" "theta s.d." $2 $3; tsd=$? |
if [ $debug -gt 0 ]; then echo "testoutput_run: listVar(I)='$listVar'" 1>&2 ; fi |
234 |
testoutput_var $1 "dynstat_salt_min" "salt minimum" $2 $3; smin=$? |
for ii in 1 2 3 4 5 6 7 8 9 ; do |
235 |
testoutput_var $1 "dynstat_salt_max" "salt maximum" $2 $3; smax=$? |
tst=`eval 'echo "$HAVE_PTR0'$ii'"'` |
236 |
testoutput_var $1 "dynstat_salt_mean" "salt mean" $2 $3; smean=$? |
#echo "-- ptr test=" $tst "number of var=" `echo $listVar | awk '{print NF}'` 1>&2 |
237 |
testoutput_var $1 "dynstat_salt_sd" "salt s.d." $2 $3; ssd=$? |
if test "x$tst" != xt ; then listVar=`echo "$listVar" | sed "s/ pt$ii..//g"` ; fi |
238 |
testoutput_var $1 "dynstat_uvel_min" "U minimum" $2 $3; umin=$? |
done |
239 |
testoutput_var $1 "dynstat_uvel_max" "U maximum" $2 $3; umax=$? |
tst=`echo $sVar $listVar | awk '{ for(i=2;i<=NF;i++){t+=($i==$1)}; print t }'` |
240 |
testoutput_var $1 "dynstat_uvel_mean" "U mean" $2 $3; umean=$? |
if test $tst != 1 ; then |
241 |
testoutput_var $1 "dynstat_uvel_sd" "U s.d." $2 $3; usd=$? |
if test $tst = 0 ; then echo "==> WARNING: selected var >$sVar< not found" 1>&2 |
242 |
testoutput_var $1 "dynstat_vvel_min" "V minimum" $2 $3; vmin=$? |
else echo "==> WARNING: found selected var >$sVar< $tst times" 1>&2 ; fi |
243 |
testoutput_var $1 "dynstat_vvel_max" "V maximum" $2 $3; vmax=$? |
echo "==> WARNING: in checked list:" $listVar 1>&2 |
244 |
testoutput_var $1 "dynstat_vvel_mean" "V mean" $2 $3; vmean=$? |
#- put it back once: |
245 |
testoutput_var $1 "dynstat_vvel_sd" "V s.d." $2 $3; vsd=$? |
listVar=" $sVar "`echo "$listVar " | sed "s/ $sVar / /g"` |
246 |
|
fi |
247 |
# This is for PTRACERS |
if [ $debug -gt 0 ]; then echo "testoutput_run: listVar(M)='$listVar'" 1>&2 ; fi |
248 |
for ii in $PTRACERS_NUM ; do |
echo "listVar='$listVar'" > $CDIR"/summary.txt" |
249 |
eval `echo "p0"$ii"_min=99"` |
allargs="" |
250 |
eval `echo "p0"$ii"_max=99"` |
for xx in $listVar |
251 |
eval `echo "p0"$ii"_mean=99"` |
do |
252 |
eval `echo "p0"$ii"_sd=99"` |
case $xx in |
253 |
tst=`eval 'echo "$HAVE_PTR0'$ii'"'` |
'PS') if [ $debug -gt 0 ] |
254 |
#echo 'tst = '$tst |
then echo testoutput_run: testoutput_var $1 cg2d_init_res 1>&2 ; fi |
255 |
if test "x$tst" = xt ; then |
testoutput_var $1 "cg2d_init_res" "Press. Solver (cg2d)" $2 $3; yy=$? |
256 |
a="trcstat_ptracer0" |
if [ $debug -gt 0 ] ; then echo testoutput_run: cg2dres=$yy 1>&2 ; fi ;; |
257 |
testoutput_var $1 "$a"$ii"_min" "p0"$ii"_min" $2 $3 |
'Cost') testoutput_var $1 "ADM precision_derivative_cost" "ADM Cost" $2 $3; yy=$? ;; |
258 |
RETVAL=$? ; eval `echo "p0"$ii"_min="$RETVAL` |
'Grad') testoutput_var $1 "ADM precision_derivative_grad" "ADM Grad" $2 $3; yy=$? ;; |
259 |
testoutput_var $1 "$a"$ii"_max" "p0"$ii"_max" $2 $3 |
'Tmn') testoutput_var $1 "dynstat_theta_min" "Theta minimum" $2 $3; yy=$? ;; |
260 |
RETVAL=$? ; eval `echo "p0"$ii"_max="$RETVAL` |
'Tmx') testoutput_var $1 "dynstat_theta_max" "Theta maximum" $2 $3; yy=$? ;; |
261 |
testoutput_var $1 "$a"$ii"_mean" "p0"$ii"_mean" $2 $3 |
'Tav') testoutput_var $1 "dynstat_theta_mean" "Theta mean" $2 $3; yy=$? ;; |
262 |
RETVAL=$? ; eval `echo "p0"$ii"_mean="$RETVAL` |
'Tsd') testoutput_var $1 "dynstat_theta_sd" "Theta Std.Dev" $2 $3; yy=$? ;; |
263 |
testoutput_var $1 "$a"$ii"_sd" "p0"$ii"_sd" $2 $3 |
'Smn') testoutput_var $1 "dynstat_salt_min" "Salt minimum" $2 $3; yy=$? ;; |
264 |
RETVAL=$? ; eval `echo "p0"$ii"_sd="$RETVAL` |
'Smx') testoutput_var $1 "dynstat_salt_max" "Salt maximum" $2 $3; yy=$? ;; |
265 |
fi |
'Sav') testoutput_var $1 "dynstat_salt_mean" "Salt mean" $2 $3; yy=$? ;; |
266 |
done |
'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=$? ;; |
268 |
allargs="$cg2dres $tmin $tmax $tmean $tsd $smin $smax $smean $ssd" |
'Umx') testoutput_var $1 "dynstat_uvel_max" "U maximum" $2 $3; yy=$? ;; |
269 |
allargs="$allargs $umin $umax $umean $usd $vmin $vmax $vmean $vsd" |
'Uav') testoutput_var $1 "dynstat_uvel_mean" "U mean" $2 $3; yy=$? ;; |
270 |
allargs="$allargs $p01_min $p01_max $p01_mean $p01_sd" |
'Usd') testoutput_var $1 "dynstat_uvel_sd" "U Std.Dev" $2 $3; yy=$? ;; |
271 |
allargs="$allargs $p02_min $p02_max $p02_mean $p02_sd" |
'Vmn') testoutput_var $1 "dynstat_vvel_min" "V minimum" $2 $3; yy=$? ;; |
272 |
allargs="$allargs $p03_min $p03_max $p03_mean $p03_sd" |
'Vmx') testoutput_var $1 "dynstat_vvel_max" "V maximum" $2 $3; yy=$? ;; |
273 |
allargs="$allargs $p04_min $p04_max $p04_mean $p04_sd" |
'Vav') testoutput_var $1 "dynstat_vvel_mean" "V mean" $2 $3; yy=$? ;; |
274 |
allargs="$allargs $p05_min $p05_max $p05_mean $p05_sd" |
'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/..$//'` |
276 |
eval "dashnum $allargs" |
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/..$//'` |
278 |
else |
testoutput_var $1 "trcstat_ptracer0"$ii"_max" "p0"$ii"_max" $2 $3; yy=$? ;; |
279 |
testoutput_ad $1 $2 "precision_grdchk_result" |
'pt1av'|'pt2av'|'pt3av'|'pt4av'|'pt5av') ii=`echo $xx | sed 's/pt//' | sed 's/..$//'` |
280 |
fi |
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/..$//'` |
282 |
|
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 ;; |
302 |
|
esac |
303 |
|
if test $xx = $sVar |
304 |
|
then allargs="$allargs > $yy <" |
305 |
|
else allargs="$allargs $yy" |
306 |
|
fi |
307 |
|
done |
308 |
|
|
309 |
|
nbVar=`echo $listVar | awk '{print NF}'` |
310 |
|
if [ $nbVar -lt $LEN_CHECK_LIST ] ; then |
311 |
|
#-- fill line (up to standard length) with dot: |
312 |
|
adNul=`expr $LEN_CHECK_LIST - $nbVar | awk '{for(i=1;i<=$1;i++){print "."}}'` |
313 |
|
echo $allargs $adNul |
314 |
|
else |
315 |
|
echo $allargs |
316 |
|
fi |
317 |
|
# <-- same processing for adjoint & forward test |
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); |
682 |
nm=$1 |
nm=$1 |
683 |
printf '%s %s %s %s' $2 $3 $4 $5 |
printf '%s %s %s %s' $2 $3 $4 $5 |
684 |
shift; shift; shift; shift; shift; |
shift; shift; shift; shift; shift; |
685 |
printf '%3s' $@ |
listPrt=$@ |
686 |
|
listRes=`echo $listPrt | sed 's/>//' | sed 's/<//'` |
687 |
|
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/ < /</' |
691 |
|
|
692 |
if [ $1 = '--' ]; then |
if [ $xx = '..' ]; then |
693 |
|
printf ' N/O ' |
694 |
|
elif [ $xx = '--' ]; then |
695 |
|
printf ' N/O ' |
696 |
|
elif [ $xx = 99 ]; then |
697 |
printf ' N/O ' |
printf ' N/O ' |
698 |
else |
else |
699 |
if [ $1 -ge $MATCH_CRIT ]; then |
if [ $xx -ge $MATCH_CRIT ]; then |
700 |
printf ' pass' |
printf ' pass' |
701 |
else |
else |
702 |
printf ' FAIL' |
printf ' FAIL' |
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: |
968 |
|
# (use default or load experiment-specific list from file "tr_checklist") |
969 |
|
# content : 1rst = main variable used to decide if it pass or FAIL |
970 |
|
# others = number of matching digits to be printed in summary.txt |
971 |
|
if test "x$ADM" = x ; then |
972 |
|
DEF_CHECK_LIST='PS PS T+ S+ U+ V+ pt1+ pt2+ pt3+ pt4+ pt5+' |
973 |
|
EMPTY_RESULTS='.. .. .. .. .. .. .. .. .. .. .. .. .. .. .. .. ..' |
974 |
|
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}'` |
976 |
|
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 |
982 |
|
|
983 |
# create the FORTRAN comparison code |
# create the FORTRAN comparison code |
984 |
createcodelet |
createcodelet |
985 |
|
|
1044 |
echo >> $SUMMARY |
echo >> $SUMMARY |
1045 |
if test "x$ADM" = x ; then |
if test "x$ADM" = x ; then |
1046 |
if [ $MATCH_CRIT -lt 10 ] ; |
if [ $MATCH_CRIT -lt 10 ] ; |
1047 |
then line_0=" "$MATCH_CRIT ; |
then line_0="default "$MATCH_CRIT ; |
1048 |
else line_0=" "$MATCH_CRIT ; fi |
else line_0="default "$MATCH_CRIT ; fi |
1049 |
line_0=" $line_0 ----T----- ----S----- ----U----- ----V-----" |
line_0="$line_0 ----T----- ----S----- ----U----- ----V-----" |
1050 |
# line_0=" ----T----- ----S----- ----U----- ----V-----" |
# line_0=" ----T----- ----S----- ----U----- ----V-----" |
1051 |
line_1="G D M c m s m s m s m s" |
line_1="G D M c m s m s m s m s" |
1052 |
line_2="E p a R g m m e . m m e . m m e . m m e ." |
line_2="E p a R g m m e . m m e . m m e . m m e ." |
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 |
1106 |
|
|
1107 |
# Verify that the testdir exists and contains previous |
# Verify that the testdir exists and contains previous |
1108 |
# results in the correct location--or skip this directory! |
# results in the correct location--or skip this directory! |
1109 |
fout= |
fout=$dir"/results/"$ref_outp |
|
if test "x$ADM" = x ; then |
|
|
fout=$dir"/results/output.txt" |
|
|
else |
|
|
fout=$dir"/results_ad/output_adm.txt" |
|
|
fi |
|
1110 |
if test ! -r $fout ; then |
if test ! -r $fout ; then |
1111 |
echo "can't read \"$fout\" -- skipping $dir" |
echo "can't read \"$fout\" -- skipping $dir" |
1112 |
continue |
continue |
1120 |
builddir="build" |
builddir="build" |
1121 |
if test ! -d $dir/$builddir ; then mkdir $dir/$builddir ; fi |
if test ! -d $dir/$builddir ; then mkdir $dir/$builddir ; fi |
1122 |
rundir="run" |
rundir="run" |
|
#rundir=$builddir |
|
1123 |
if test ! -d $dir/$rundir ; then |
if test ! -d $dir/$rundir ; then |
1124 |
rundir=$builddir |
rundir=$builddir |
1125 |
fi |
fi |
1126 |
|
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 |
|
1127 |
BUILD_DIR=$dir/$builddir |
BUILD_DIR=$dir/$builddir |
1128 |
|
|
1129 |
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 |
1165 |
fi |
fi |
1166 |
echo |
echo |
1167 |
unset genmake makedepend make run |
unset genmake makedepend make run |
1168 |
results='-- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- --' |
results=$EMPTY_RESULTS |
1169 |
|
|
1170 |
# Create an output dir for each OPTFILE/tdir combination |
# Create an output dir for each OPTFILE/tdir combination |
1171 |
rel_CDIR=$DRESULTS"/"$dir |
rel_CDIR=$DRESULTS"/"$dir |
1185 |
&& makemodel $dir/$builddir && make=Y \ |
&& makemodel $dir/$builddir && make=Y \ |
1186 |
&& linkdata $dir/$rundir $input_dirs \ |
&& linkdata $dir/$rundir $input_dirs \ |
1187 |
&& runmodel $dir/$rundir && run=Y \ |
&& runmodel $dir/$rundir && run=Y \ |
1188 |
&& results=`testoutput_run $dir $rundir "txt"` |
&& results=`testoutput_run $dir $rundir $ref_outp` |
1189 |
fi |
fi |
1190 |
|
|
1191 |
echo |
echo |
1192 |
if test "x$ADM" = x ; then |
# --> same processing for adjoint & forward test |
1193 |
fres=`formatresults $dir ${genmake:-N} ${makedepend:-N} ${make:-N} ${run:-N} $results` |
fres=`formatresults $dir ${genmake:-N} ${makedepend:-N} ${make:-N} ${run:-N} $results` |
1194 |
echo |
echo |
1195 |
echo "$fres" >> $SUMMARY |
echo "$fres" | sed 's/ 99/ --/g' | sed 's/ > />/' | sed 's/ < /</' >> $SUMMARY |
1196 |
echo "fresults='$fres'" > $CDIR"/summary.txt" |
touch $CDIR"/summary.txt" |
1197 |
|
echo "fresults='$fres'" | sed 's/ 99/ --/g' >> $CDIR"/summary.txt" |
1198 |
echo "MACH='$MACH'" >> $CDIR"/summary.txt" |
echo "MACH='$MACH'" >> $CDIR"/summary.txt" |
1199 |
echo "UNAMEA='$UNAMEA'" >> $CDIR"/summary.txt" |
echo "UNAMEA='$UNAMEA'" >> $CDIR"/summary.txt" |
1200 |
echo "DATE='$DATE'" >> $CDIR"/summary.txt" |
echo "DATE='$DATE'" >> $CDIR"/summary.txt" |
1202 |
|
|
1203 |
for ex in $extra_runs ; do |
for ex in $extra_runs ; do |
1204 |
unset run |
unset run |
1205 |
results='-- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- -- --' |
results=$EMPTY_RESULTS |
1206 |
# Create an output dir for each OPTFILE/tdir.ex combination |
# Create an output dir for each OPTFILE/tdir.ex combination |
1207 |
rel_CDIR=$DRESULTS"/"$dir"."$ex |
rel_CDIR=$DRESULTS"/"$dir"."$ex |
1208 |
mkdir $rel_CDIR |
mkdir $rel_CDIR |
1211 |
run_clean $dir/tr_run.$ex |
run_clean $dir/tr_run.$ex |
1212 |
linkdata $dir/tr_run.$ex input.$ex input |
linkdata $dir/tr_run.$ex input.$ex input |
1213 |
runmodel $dir/tr_run.$ex && run=Y \ |
runmodel $dir/tr_run.$ex && run=Y \ |
1214 |
&& 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}` |
|
1215 |
fres=`formatresults $dir ${genmake:-N} ${makedepend:-N} ${make:-N} ${run:-N} $results` |
fres=`formatresults $dir ${genmake:-N} ${makedepend:-N} ${make:-N} ${run:-N} $results` |
1216 |
fres="$fres.$ex" |
fres="$fres.$ex" |
1217 |
echo |
echo |
1218 |
echo "$fres" >> $SUMMARY |
echo "$fres" | sed 's/ 99/ --/g' | sed 's/ > />/' | sed 's/ < /</' >> $SUMMARY |
1219 |
echo "fresults='$fres'" > $CDIR"/summary.txt" |
touch $CDIR"/summary.txt" |
1220 |
|
echo "fresults='$fres'" | sed 's/ 99/ --/g' >> $CDIR"/summary.txt" |
1221 |
echo "MACH='$MACH'" >> $CDIR"/summary.txt" |
echo "MACH='$MACH'" >> $CDIR"/summary.txt" |
1222 |
echo "UNAMEA='$UNAMEA'" >> $CDIR"/summary.txt" |
echo "UNAMEA='$UNAMEA'" >> $CDIR"/summary.txt" |
1223 |
echo "DATE='$DATE'" >> $CDIR"/summary.txt" |
echo "DATE='$DATE'" >> $CDIR"/summary.txt" |
1224 |
echo "tdir='$dir.$ex'" >> $CDIR"/summary.txt" |
echo "tdir='$dir.$ex'" >> $CDIR"/summary.txt" |
1225 |
if test "x$POSTCLEAN" = xt ; then |
if test "x$POSTCLEAN" = xt ; then |
1226 |
run_clean $dir/tr_run.$ex |
run_clean $dir/tr_run.$ex |
1227 |
fi |
fi |
1228 |
done |
done |
1229 |
else |
# <-- same processing for adjoint & forward test |
1230 |
fres=`printf '%s %s %s %s' ${genmake:-N} ${makedepend:-N} ${make:-N} ${run:-N}` |
if test "x$ADM" = xt ; then |
1231 |
fres=$fres"$results $dir" |
head -1 $dir/$builddir/taf_ad.log >> $CDIR"/summary.txt" |
1232 |
echo |
grep -A3 'Seconds in section "ALL' $dir/$rundir/$OUTPUTFILE \ |
|
echo "$fres" >> $SUMMARY |
|
|
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 \ |
|
1233 |
>> $CDIR"/summary.txt" |
>> $CDIR"/summary.txt" |
1234 |
fi |
fi |
1235 |
|
|
1278 |
rm -f tr_cmpnum.c tr_cmpnum |
rm -f tr_cmpnum.c tr_cmpnum |
1279 |
|
|
1280 |
if test "x$CLEANUP" != xt ; then |
if test "x$CLEANUP" != xt ; then |
1281 |
cat $SUMMARY | sed 's/ -- -- -- --//g' |
cat $SUMMARY | sed 's/ \. \. \. \. \. \. \. \. \. \. \. \. //' |
1282 |
if test -e tr_out.txt ; then |
if test -e tr_out.txt ; then |
1283 |
mv tr_out.txt tr_out.txt.old |
mv tr_out.txt tr_out.txt.old |
1284 |
fi |
fi |
1285 |
cat $SUMMARY | sed 's/ -- -- -- --//g' > tr_out.txt |
cat $SUMMARY | sed '/^[YN] [YN] [YN] [YN]/ s/ \. //g' > tr_out.txt |
1286 |
fi |
fi |
1287 |
|
|
1288 |
if test "x$DELDIR" = xt ; then |
if test "x$DELDIR" = xt ; then |