Subversion Repositories wimsdev

#! /bin/sh
#
# Mandatary parms: $1=plot number. $2=session number.
# environment parms: ins_source.
# reserved words: insplot_set (setup), insplot_data, ins_border (border size).

if [ -z "$2" ]; then exit; fi

CONVERT_DIS=2
CONVERT_OPT="-quality 90 -dispose $CONVERT_DIS"

GNUPLOT_HEADER='Pi=pi
PI=pi
e=exp(1)
E=exp(1)
sign(x) = sgn(x)
ln(x) = log(x)
lg(x) = log10(x)
ch(x) = cosh(x)
sh(x) = sinh(x)
th(x) = tanh(x)
arccos(x) = acos(x)
arcsin(x) = asin(x)
tg(x)     = tan(x)
arctan(x) = atan(x)
arctg(x)  = atan(x)
Argch(x)  = acosh(x)
Argsh(x)  = asinh(x)
Argth(x)  = atanh(x)
argch(x)  = acosh(x)
argsh(x)  = asinh(x)
argth(x)  = atanh(x)
cotan(x)  = (1/tan(x))
cot(x)    = (1/tan(x))
ctg(x)    = (1/tan(x))
sec(x)    = (1/cos(x))
csc(x)    = (1/sin(x))
max(x,y)  = (x>y? x : y)
min(x,y)  = (x<y? x : y)
round(x)  = floor(x+0.5)
rint(x)   = floor(x+0.5)
set samples 400
set colors classic
'

plot_cmd_2D=plot
plot_cmd_3D=splot
insplot_cmd=
if [ "$insplot_method" = "2D" ]; then insplot_cmd=$plot_cmd_2D; fi
if [ "$insplot_method" = "3D" ]; then insplot_cmd=$plot_cmd_3D; fi
if [ -z $insplot_cmd ]; then exit; fi

rm -f $2/insert$w_wims_inssub-$1* 2>/dev/null

# data file names
plot_parm=`echo "$ins_source" | sed "s!insplot_data!$2/insplot_data!g"`
insplot_parm=`echo "$insplot_source" | sed "s!insplot_data!$2/insplot_data!g"`

if [ ! -z "$w_insplot_data" ]; then
  cat >$2/insplot_data <<@
$w_insplot_data
@
fi
if [ -z "$w_ins_anim_loop" ]; then w_ins_anim_loop=0; fi
insplot_colors=`echo $w_insplot_colors | tr ';\12' '  '`

if [ ! -z $w_insplot_transparent ]; then
 gnuplot_transp=transparent
# A bad fix for convert transparency incompatibility trouble
# convert_transp="-transparent $w_insplot_transparent"
 convert_transp="-transparent rgb(255,255,255)"
else
 convert_transp=
fi

if [ "$w_anim_format" != "mng" ]; then
 w_anim_format=gif
fi

if [ -z "$w_gnuplot_format" ] || [ "$w_gnuplot_format" = "default" ]; then
 intermediate_format=$w_ins_format
else
 intermediate_format=$w_gnuplot_format
fi

###fix compatibility with gnuplot 4
### and adress new chroot...
chroot_version=`cat $w_wims_home/chroot/version`
if [ "$chroot_version" = "4" ] ; then
    gnuplot_version=4
else
    gnuplot_version=`gnuplot --version | awk '{ print $2 }' | awk -F '.'  '{ print $1 }'`
fi
if [ "$gnuplot_version" -lt "4" ] || [ -z "$gnuplot_version" ]; then
  insplot_termset="set terminal $intermediate_format $gnuplot_transp $w_insplot_font size 500,400 $insplot_colors"
else
  testsize=`echo $insplot_set | perl -e 'if (<> =~ /size (\d+\.\d+),(\d+\.\d+)/) {$sx=$1*500 ; $sy=$2*400 ; print "size $sx,$sy" ; } else {print "size 500,400"}'`
  insplot_termset="set terminal $intermediate_format $gnuplot_transp $w_insplot_font $testsize $insplot_colors"
# insplot_set=`echo $insplot_set | perl -e '$ins_set=<> ; $ins_set =~ s/(set +size +\d+\.\d+,\d+\.\d+)// ; $ins_set =~ s/set/\nset/g ; print $ins_set;'`
  insplot_set=`echo "$insplot_set" | sed "s/set[        ]*size *[0-9]\.[0-9]*,[0-9]\.[0-9]*//g"`
  insplot_set=`echo "$insplot_set" | sed "s/set +no/unset /g"`
fi

### todo : dans les nouvelles versions de gnuplot, transformer "set noxxx" en "unset xxx"
if [ ! -z "$multiplot" ]; then
 plotcmd="set multiplot"
 i=1
 until [ $i -gt $multiplot ]; do
  tt=`echo "$insplot_parm" | awk 'FNR=='$i' {print;exit}'`
  t2=`echo "$w_insplot_split" \
        | awk -F ';' 'FNR=='$i' {for(j=1;j<=NF;j++) print "set "$j;exit}'`
  t2=`echo "$t2" | sed "s/set +no/unset /g"`
  c1=`echo $tt | awk '{print $1}'`
  c2=`echo $tt | awk '{print substr($0,3)}'`
  if [ $c1 = "2D" ]; then
    cc="$plot_cmd_2D $c2"
  else
    if [ $c1 = "3D" ]; then
      cc="$plot_cmd_3D $c2"
    else
      cc="$insplot_cmd $tt"
    fi
  fi
  plotcmd="$plotcmd
        $t2
        $cc"
  i=$(( $i+1 ))
 done
 plotcmd="$plotcmd
        unset multiplot"
else
 plotcmd="$insplot_cmd $plot_parm"
fi

plot_content=
anim_list=
i=0
while [ $i -lt $ins_anim_frames ]; do
plot_content="$plot_content
 step=$i
 animstep=$i
 anim_step=$i
 s=$i/$ins_anim_frames.0
 $insplot_set
 $insplot_termset
 set output \"$2/insert$w_wims_inssub-$1-$i.$intermediate_format\"
 $plotcmd
"
anim_list="$anim_list $2/insert$w_wims_inssub-$1-$i.$intermediate_format"
i=$(( $i+1 ))
done

bin/ch..root gnuplot <<@
$GNUPLOT_HEADER
$plot_content
@

echo "$GNUPLOT_HEADER
$plot_content" >$2/gnuplot.cmd

alist=
for i in $anim_list
do
 if [ -s $i ]; then alist="$alist $i"; fi
done

if [ $ins_anim_frames -lt 2 ]; then
 if [ $w_ins_format = $intermediate_format ]; then
  bin/ch..root mv $2/insert$w_wims_inssub-$1-0.$w_ins_format $2/insert$w_wims_inssub-$1.$w_ins_format
 else
  bin/ch..root convert $CONVERT_OPT $convert_transp -density $w_ins_density $2/insert$w_wims_inssub-$1-0.$intermediate_format $2/insert$w_wims_inssub-$1.$w_ins_format 2>$2/convert.err || bin/convert..switch $2/convert.err
 fi
else
 w_ins_format=$w_anim_format
 if [ "$intermediate_format" = "gif" ] && [ "$w_ins_format" = "gif" ]; then
  whirlgif -disp prev -time $ins_anim_delay -loop $w_ins_anim_loop -o $2/insert$w_wims_inssub-$1.gif $alist
 else
  bin/ch..root convert $CONVERT_OPT $convert_transp -delay $ins_anim_delay -loop $w_ins_anim_loop $alist $2/insert$w_wims_inssub-$1.$w_ins_format 2>$2/convert.err || bin/convert..switch $2/convert.err
 fi
fi

# Remove zero-sized jpeg file, for otherwise netscape doesn't understand,
# and waits forever trying to interprete it.
if [ ! -s $2/insert$w_wims_inssub-$1.$w_ins_format ]; then rm -f $2/insert$w_wims_inssub-$1.$w_ins_format; fi

echo "$plot_content" >$2/insplot_cmd