superinv.txt

射线追踪程序

#!/bin/sh
test $# -lt 1 && exec echo $0 name

RCP="scp -q -C -oProtocol=1"
RSH="ssh -q -1"

rm -rf global
mkdir global

for m in `grep -v "^#" MACHINES`
do  $RCP /usr/local/share/proj/proj_def.dat inv $1.dat $1.lay $1.cnf $m:tomo3d-data/.
    $RSH $m "(cd tomo3d-data; ./inv $1 clean)"
    touch $m.log
    xterm -fn 5x7 -T $m -e tail -f $m.log&
    rm -f $m.lock
done

echo 10 > oa.tmp
echo 1 > layer.tmp
ENDLAYER=`cat $1.lay | wc -l | awk '{ print $1-2 }'`

while [ -f oa.tmp ]
do  for m in `grep -v "^#" MACHINES`
    do  test -f $m.lock || (\
            touch $m.lock
            oa=`cat oa.tmp`
            layer=`cat layer.tmp`
            echo Running $layer/$ENDLAYER dir=$oa,layer=$layer on $m
            $RCP layer.tmp $m:tomo3d-data/LAYER.file
            $RCP oa.tmp    $m:tomo3d-data/OA.file
            $RCP inv $1.dat $1.lay $1.cnf $m:tomo3d-data/.

            echo `expr $layer + 1` > layer.tmp
            if [ $layer -ge $ENDLAYER ]; then
                echo `expr $oa + 1` > oa.tmp
                echo 1 > layer.tmp
            fi
            ( echo "$RSH $m \"(cd tomo3d-data; ./inv $1 A RT B $oa $layer)\" > $m.log 2>&1"
              echo "$RCP $m:tomo3d-data/all.geo global/$1.geo-$layer-$oa-$m"
              echo echo $m finished
              echo rm -f $m.lock
            ) > $m.script
            sh $m.script &
        )
        sleep 1
    done
    sleep 5
    test `cat oa.tmp` -gt 17 && rm -f oa.tmp
done

while `ls *lock`; do sleep 1; done