# 18 Aug 2017
# Each of these plot commands failed on some version between 5.0 and 5.3.
# I collect them here to use as unit tests for continued success
# of the code added to repair breakage.
#
# Usage:  call 'sampling.dem' <testno>
#

set urange [1:100]
set vrange [1:100]
set xrange [1:100]
set key opaque box

if (ARGC < 1) ARG1 = 0

set trange [20:50]

if (ARG1 <= 1) {
    print "test 1: explicit trange distinct from xrange"
    plot '+' using 1:(10. + sin($1)) title "trange [20:50]"
    pause -1 "Hit <cr> to continue"
}

if (ARG1 <= 2) {
    print "test 2: range set by 'sample' keyword, linear x axis"
    plot sample [t=20:50] '+' using (t):(10. + sin(t))
    pause -1 "Hit <cr> to continue"
}


set log x

if (ARG1 <= 3) {
    print "test 3: range set by 'sample' keyword, logscale x axis"
    plot sample [t=20:50] '+' using (t):(10. + sin(t))
    pause -1 "Hit <cr> to continue"
}

unset log

if (ARG1 <= 4) {
    print "test 4: splot '++' with autoscaled y (linear xy)"
    splot '++' using 1:2:($1*$2)
    pause -1 "Hit <cr> to continue"
}

set log xy

if (ARG1 <= 5) {
    print "test 5: splot '++' with autoscaled y (logscale xy)"
    splot '++' using 1:2:($1*$2)
    pause -1 "Hit <cr> to continue"
}

unset log

if (ARG1 <= 6) {
    print "test 6: plot '++' with image (linear xy)"
    plot '++' using 1:2:($1*$2) with image
    pause -1 "Hit <cr> to continue"
}


set log xy

if (ARG1 <= 7) {
    print "test 7: plot '++' with image (logscale xy)"
    plot '++' using 1:2:($1*$2) with image
    pause -1 "Hit <cr> to continue"
}

reset 
#
#  piecewise functions along nonlinear x axis
#
unset border
set key center top reverse Left
set xzeroaxis
set yzeroaxis
set xtics axis out scale 1,8
set xtics add  (1.00000 1, 6.28319 1)
set ytics axis
set xrange [ .1 : 12 ] noreverse nowriteback



if (ARG1 <= 8) {
    print "test 8: multiple sampling ranges in one 2D plot command"
    set multiplot layout 2,1
    set title "Piecewise function sampling along linear x" font ",13"
    plot sample [*:1] x, [1:2.*pi] cos(x), [2.*pi:10] (x-8)**2
    set nonlinear x via x inv x
    set title "nonlinear (identity mapped) x" font ",13"
    plot sample [*:1] x, [1:2.*pi] cos(x), [2.*pi:10] (x-8)**2
    unset multiplot
    pause -1 "Hit <cr> to continue"
}

reset

set xrange [1:100]
set yrange [1:100]
set urange [0:100]
set vrange [70:90]
set style data lines
set xyplane at 0

set hidden3d
set title "3D sampling range distinct from plot x/y range"
print "test 9: 3D sampling range distinct from plot x/y range"

splot '++' using 1:2:($1*25.*sin($2/10)), \
      [u=30:70][v=0:50] '++' using 1:2:(u*v), \
      [u=40:80][v=30:60] '++' using (u):(v):(u*sqrt(v)) lt 4, \
      [u=1:100][v=500:1000] '++' using (90):(u):(v) lt 6

pause -1 "Hit <cr> to continue"

reset

unset key
set view map
set xrange [ 20:100 ]
set yrange [ 0:80 ]
set size ratio 1.0
radius(x,y) = 10 * sqrt((x-50)**2 + (y-50)**2)
alpha(x,y) = radius(x,y) > 255 ? 0 : 255-radius(x,y)
set title "3D custom sampling on u and v using pseudofile '++'"
print "test 10: splot '++' with explicit sampling intervals"

splot sample [u=25:75:1][v=25:75:4] '++' \
      using 1:2:0:(u*4-50):(v*4-50):(255*sin(u*v/2500.)):(alpha(u,v)) with rgbalpha, \
      [u=25:75:4][v=25:75:1] '++' \
      using (u+20):(v-20):0:(u*4-50):(v*4-50):(255*sin(u*v/2500.)):(alpha(u,v)) with rgbalpha

pause -1 "Hit <cr> to continue"

set title "2D custom sampling on u and v using pseudofile '++'"
print "test 10: plot '++' with explicit sampling intervals"

plot sample [u=25:75:1][v=25:75:4] '++' \
      using 1:2:(u*4-50):(v*4-50):(255*sin(u*v/2500.)):(alpha(u,v)) with rgbalpha, \
      [u=25:75:4][v=25:75:1] '++' \
      using (u+20):(v-20):(u*4-50):(v*4-50):(255*sin(u*v/2500.)):(alpha(u,v)) with rgbalpha

pause -1 "Hit <cr> to continue"

reset

