Subversion Repositories wimsdev

Rev

Rev 7152 | Go to most recent revision | Details | Compare with Previous | Last modification | View Log | RSS feed

Rev Author Line No. Line
14 reyssat 1
#! /bin/sh
2
#
3
# Mandatary parms: $1=plot number. $2=session number.
4
# environment parms: ins_source.
5
# reserved words: insplot_set (setup), insplot_data, ins_border (border size).
6
 
7
if [ -z "$2" ]; then exit; fi
8
 
9
CONVERT_DIS=2
10
CONVERT_OPT="-quality 90 -dispose $CONVERT_DIS"
11
 
12
GNUPLOT_HEADER='Pi=pi
13
PI=pi
14
e=exp(1)
15
E=exp(1)
16
sign(x) = sgn(x)
17
ln(x) = log(x)
18
lg(x) = log10(x)
19
ch(x) = cosh(x)
20
sh(x) = sinh(x)
21
th(x) = tanh(x)
22
arccos(x) = acos(x)
23
arcsin(x) = asin(x)
24
tg(x)	  = tan(x)
25
arctan(x) = atan(x)
26
arctg(x)  = atan(x)
27
Argch(x)  = acosh(x)
28
Argsh(x)  = asinh(x)
29
Argth(x)  = atanh(x)
30
argch(x)  = acosh(x)
31
argsh(x)  = asinh(x)
32
argth(x)  = atanh(x)
33
cotan(x)  = (1/tan(x))
34
cot(x)	  = (1/tan(x))
35
ctg(x)    = (1/tan(x))
36
sec(x)	  = (1/cos(x))
37
csc(x)	  = (1/sin(x))
38
max(x,y)  = (x>y? x : y)
39
min(x,y)  = (x<y? x : y)
40
round(x)  = floor(x+0.5)
41
rint(x)	  = floor(x+0.5)
42
set samples 400
43
'
44
 
45
plot_cmd_2D=plot
46
plot_cmd_3D=splot
47
insplot_cmd=
48
if [ "$insplot_method" = "2D" ]; then insplot_cmd=$plot_cmd_2D; fi
49
if [ "$insplot_method" = "3D" ]; then insplot_cmd=$plot_cmd_3D; fi
50
if [ -z $insplot_cmd ]; then exit; fi
51
 
52
rm -f $2/insert$w_wims_inssub-$1* 2>/dev/null
53
 
54
# data file names
55
plot_parm=`echo "$ins_source" | sed "s!insplot_data!$2/insplot_data!g"`
56
insplot_parm=`echo "$insplot_source" | sed "s!insplot_data!$2/insplot_data!g"`
57
 
58
if [ ! -z "$w_insplot_data" ]; then
59
  cat >$2/insplot_data <<@
60
$w_insplot_data
61
@
12023 bpr 62
fi
14 reyssat 63
if [ -z "$w_ins_anim_loop" ]; then w_ins_anim_loop=0; fi
64
insplot_colors=`echo $w_insplot_colors | tr ';\12' '  '`
65
 
66
if [ ! -z $w_insplot_transparent ]; then
67
 gnuplot_transp=transparent
68
# A bad fix for convert transparency incompatibility trouble
69
# convert_transp="-transparent $w_insplot_transparent"
1058 bpr 70
 convert_transp="-transparent rgb(255,255,255)"
14 reyssat 71
else
72
 convert_transp=
73
fi
74
 
75
if [ "$w_anim_format" != "mng" ]; then
76
 w_anim_format=gif
77
fi
78
 
79
if [ -z "$w_gnuplot_format" ] || [ "$w_gnuplot_format" = "default" ]; then
80
 intermediate_format=$w_ins_format
81
else
82
 intermediate_format=$w_gnuplot_format
83
fi
84
 
12023 bpr 85
###fix compatibility with gnuplot 4
4303 bpr 86
### and adress new chroot...
87
chroot_version=`cat $w_wims_home/chroot/version`
88
if [ "$chroot_version" = "4" ] ; then
89
    gnuplot_version=4
90
else
91
    gnuplot_version=`gnuplot --version | awk '{ print $2 }' | awk -F '.'  '{ print $1 }'`
92
fi
892 bpr 93
if [ "$gnuplot_version" -lt "4" ] || [ -z "$gnuplot_version" ]; then
94
  insplot_termset="set terminal $intermediate_format $gnuplot_transp $w_insplot_font size 500,400 $insplot_colors"
95
else
96
  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"}'`
97
  insplot_termset="set terminal $intermediate_format $gnuplot_transp $w_insplot_font $testsize $insplot_colors"
98
# 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;'`
7152 bpr 99
  insplot_set=`echo "$insplot_set" | sed "s/set[ 	]*size *[0-9]\.[0-9]*,[0-9]\.[0-9]*//g"`
12023 bpr 100
  insplot_set=`echo "$insplot_set" | sed "s/set +no/unset /g"`
892 bpr 101
fi
102
 
1059 bpr 103
### todo : dans les nouvelles versions de gnuplot, transformer "set noxxx" en "unset xxx"
14 reyssat 104
if [ ! -z "$multiplot" ]; then
105
 plotcmd="set multiplot"
106
 i=1
107
 until [ $i -gt $multiplot ]; do
108
  tt=`echo "$insplot_parm" | awk 'FNR=='$i' {print;exit}'`
109
  t2=`echo "$w_insplot_split" \
110
	| awk -F ';' 'FNR=='$i' {for(j=1;j<=NF;j++) print "set "$j;exit}'`
12023 bpr 111
  t2=`echo "$t2" | sed "s/set +no/unset /g"`
14 reyssat 112
  c1=`echo $tt | awk '{print $1}'`
113
  c2=`echo $tt | awk '{print substr($0,3)}'`
12023 bpr 114
  if [ $c1 = "2D" ]; then
14 reyssat 115
    cc="$plot_cmd_2D $c2"
116
  else
117
    if [ $c1 = "3D" ]; then
118
      cc="$plot_cmd_3D $c2"
119
    else
120
      cc="$insplot_cmd $tt"
121
    fi
122
  fi
123
  plotcmd="$plotcmd
124
	$t2
125
	$cc"
126
  i=$(( $i+1 ))
127
 done
128
 plotcmd="$plotcmd
12023 bpr 129
	unset multiplot"
14 reyssat 130
else
131
 plotcmd="$insplot_cmd $plot_parm"
132
fi
133
 
134
plot_content=
135
anim_list=
136
i=0
137
while [ $i -lt $ins_anim_frames ]; do
138
plot_content="$plot_content
139
 step=$i
140
 animstep=$i
141
 anim_step=$i
142
 s=$i/$ins_anim_frames.0
143
 $insplot_set
144
 $insplot_termset
145
 set output \"$2/insert$w_wims_inssub-$1-$i.$intermediate_format\"
146
 $plotcmd
147
"
148
anim_list="$anim_list $2/insert$w_wims_inssub-$1-$i.$intermediate_format"
149
i=$(( $i+1 ))
150
done
151
 
152
bin/ch..root gnuplot <<@
153
$GNUPLOT_HEADER
154
$plot_content
155
@
156
 
157
echo "$GNUPLOT_HEADER
158
$plot_content" >$2/gnuplot.cmd
159
 
160
alist=
161
for i in $anim_list
162
do
163
 if [ -s $i ]; then alist="$alist $i"; fi
164
done
165
 
166
if [ $ins_anim_frames -lt 2 ]; then
167
 if [ $w_ins_format = $intermediate_format ]; then
168
  bin/ch..root mv $2/insert$w_wims_inssub-$1-0.$w_ins_format $2/insert$w_wims_inssub-$1.$w_ins_format
169
 else
170
  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
171
 fi
172
else
173
 w_ins_format=$w_anim_format
174
 if [ "$intermediate_format" = "gif" ] && [ "$w_ins_format" = "gif" ]; then
175
  whirlgif -disp prev -time $ins_anim_delay -loop $w_ins_anim_loop -o $2/insert$w_wims_inssub-$1.gif $alist
176
 else
177
  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
178
 fi
179
fi
180
 
181
# Remove zero-sized jpeg file, for otherwise netscape doesn't understand,
182
# and waits forever trying to interprete it.
183
if [ ! -s $2/insert$w_wims_inssub-$1.$w_ins_format ]; then rm -f $2/insert$w_wims_inssub-$1.$w_ins_format; fi
184
 
185
echo "$plot_content" >$2/insplot_cmd
186