plplot 0.0.0

data/examples/x15.rb

@@ -0,0 +1,266 @@
+#!/usr/bin/env ruby
+$-w = true if $0 == __FILE__
+
+require 'rubygems'
+require 'plplot'
+include PLplot
+include Math
+
+# Shade plot demo.
+#
+# Maurice LeBrun
+# IFS, University of Texas at Austin
+# 31 Aug 1993
+
+XPTS = 35 # Data points in x
+YPTS = 46 # Data points in y
+
+
+#--------------------------------------------------------------------------
+# cmap1_init1
+#
+# Initializes color map 1 in HLS space.
+#--------------------------------------------------------------------------
+
+def cmap1_init1
+  i = [
+    0,           # left boundary
+    0.45,        # just before center
+    0.55,        # just after center
+    1            # right boundary
+  ]
+
+  h = [
+    260,         # hue -- low: blue-violet
+    260,         # only change as we go over vertex
+    20,          # hue -- high: red
+    20           # keep fixed
+  ]
+
+  l = [
+    0.5,         # lightness -- low
+    0.0,         # lightness -- center
+    0.0,         # lightness -- center
+    0.5          # lightness -- high
+  ]
+
+  s = [
+    1,           # maximum saturation
+    1,           # maximum saturation
+    1,           # maximum saturation
+    1            # maximum saturation
+  ]
+
+  plscmap1l( 0, i, h, l, s )
+end
+
+#--------------------------------------------------------------------------\
+# cmap1_init2
+#
+# Initializes color map 1 in HLS space.
+#--------------------------------------------------------------------------/
+
+def cmap1_init2
+  i = [
+    0,           # left boundary
+    0.45,        # just before center
+    0.55,        # just after center
+    1            # right boundary
+  ]
+
+  h = [
+    260,         # hue -- low: blue-violet
+    260,         # only change as we go over vertex
+    20,          # hue -- high: red
+    20           # keep fixed
+  ]
+
+  l = [
+    0.6,         # lightness -- low
+    0.0,         # lightness -- center
+    0.0,         # lightness -- center
+    0.6          # lightness -- high
+  ]
+
+  s = [
+    1,           # saturation -- low
+    0.5,         # saturation -- center
+    0.5,         # saturation -- center
+    1            # saturation -- high
+  ]
+
+  plscmap1l( 0, i, h, l, s )
+end
+
+#--------------------------------------------------------------------------\
+# plot1
+#
+# Illustrates a single shaded region.
+#--------------------------------------------------------------------------/
+
+def plot1(z)
+
+  zmin, zmax = plMinMax2dGrid(z)
+
+  pladv( 0 )
+  plvpor( 0.1, 0.9, 0.1, 0.9 )
+  plwind( -1.0, 1.0, -1.0, 1.0 )
+
+  # Plot using identity transform
+
+  shade_min = zmin + ( zmax - zmin ) * 0.4
+  shade_max = zmin + ( zmax - zmin ) * 0.6
+  sh_cmap   = 0
+  sh_color  = 7
+  sh_width  = 2
+  min_color = 9
+  max_color = 2
+  min_width = 2
+  max_width = 2
+
+  plpsty( 8 )
+  plshade(z, -1..1, -1..1,
+          shade_min..shade_max,
+          sh_color, sh_width,
+          min_color..max_color,
+          min_width..max_width,
+          1)
+
+  plcol0( 1 )
+  plbox( "bcnst", 0.0, 0, "bcnstv", 0.0, 0 )
+  plcol0( 2 )
+  pllab( "distance", "altitude", "Bogon flux" )
+end
+
+#--------------------------------------------------------------------------
+# plot2
+#
+# Illustrates multiple adjacent shaded regions, using different fill
+# patterns for each region.
+#--------------------------------------------------------------------------
+
+def plot2(z)
+
+  inc = [
+    [ 450       ], [ -450       ], [ 0      ], [ 900      ],
+    [ 300       ], [  450, -450 ], [ 0, 900 ], [   0, 450 ],
+    [ 450, -450 ], [    0,  900 ]
+  ]
+
+  del = [
+    [ 2000       ], [ 2000       ], [ 2000       ],
+    [ 2000       ], [ 2000       ], [ 2000, 2000 ],
+    [ 2000, 2000 ], [ 2000, 2000 ], [ 4000, 4000 ],
+    [ 4000, 2000 ]
+  ]
+
+  sh_width = 2
+
+  pladv( 0 )
+  plvpor( 0.1, 0.9, 0.1, 0.9 )
+  plwind( -1.0, 1.0, -1.0, 1.0 )
+
+  # Plot using identity transform
+
+  zmin, zmax = plMinMax2dGrid(z)
+
+  10.times do |i|
+    shade_min = zmin + ( zmax - zmin ) * i / 10.0
+    shade_max = zmin + ( zmax - zmin ) * ( i + 1 ) / 10.0
+    sh_color  = i + 6
+    plpat( inc[i], del[i] )
+
+    plshade(z, -1..1, -1..1,
+            shade_min..shade_max,
+            sh_color, sh_width,
+            nil, nil,
+            1)
+  end
+
+  plcol0( 1 )
+  plbox( "bcnst", 0.0, 0, "bcnstv", 0.0, 0 )
+  plcol0( 2 )
+  pllab( "distance", "altitude", "Bogon flux" )
+end
+
+#--------------------------------------------------------------------------
+# plot3
+#
+# Illustrates shaded regions in 3d, using a different fill pattern for
+# each region.
+#--------------------------------------------------------------------------
+
+def plot3
+  xx = [
+    [ -1.0, 1.0, 1.0, -1.0, -1.0 ],
+    [ -1.0, 1.0, 1.0, -1.0, -1.0 ]
+  ]
+  yy = [
+    [  1.0,  1.0, 0.0, 0.0,  1.0 ],
+    [ -1.0, -1.0, 0.0, 0.0, -1.0 ]
+  ]
+  zz = [
+    [ 0.0, 0.0, 1.0, 1.0, 0.0 ],
+    [ 0.0, 0.0, 1.0, 1.0, 0.0 ]
+  ]
+
+  pladv( 0 )
+  plvpor( 0.1, 0.9, 0.1, 0.9 )
+  plwind( -1.0, 1.0, -1.0, 1.0 )
+  plw3d( 1.0, 1.0, 1.0, -1.0, 1.0, -1.0, 1.0, 0.0, 1.5, 30, -40 )
+
+  # Plot using identity transform
+
+  plcol0( 1 )
+  plbox3( "bntu", "X", 0.0, 0, "bntu", "Y", 0.0, 0, "bcdfntu", "Z", 0.5, 0 )
+  plcol0( 2 )
+  pllab( "", "", "3-d polygon filling" )
+
+  plcol0( 3 )
+  plpsty( 1 )
+  plline3( xx[0], yy[0], zz[0] )
+  plfill3( xx[0][0,4], yy[0][0,4], zz[0][0,4] )
+  plpsty( 2 )
+  plline3( xx[1], yy[1], zz[1] )
+  plfill3( xx[1][0,4], yy[1][0,4], zz[1][0,4] )
+end
+
+#--------------------------------------------------------------------------
+# Does a variety of shade plots.
+#--------------------------------------------------------------------------
+
+
+# Parse and process command line arguments
+
+PLOptionParser.parse!
+
+# Set up color map 0
+
+#plscmap0n(3)
+
+# Set up color map 1
+
+cmap1_init2
+
+# Initialize plplot
+
+plinit
+
+# Set up data array
+
+z = NArray.float(YPTS, XPTS)
+
+XPTS.times do |i|
+  xx = ( i - ( XPTS.to_f / 2 ) ) / ( XPTS.to_f / 2 )
+  YPTS.times do |j|
+    yy = ( j - ( YPTS.to_f / 2 ) ) / ( YPTS.to_f / 2 ) - 1.0
+
+    z[j,i] = xx * xx - yy * yy + ( xx - yy ) / ( xx * xx + yy * yy + 0.1 )
+  end
+end
+
+plot1(z)
+plot2(z)
+plot3
+
+plend

data/examples/x16.rb

@@ -0,0 +1,349 @@
+#!/usr/bin/env ruby
+$-w = true if $0 == __FILE__
+
+require 'rubygems'
+require 'plplot'
+include PLplot
+include Math
+
+# plshade demo, using color fill.
+#
+# Maurice LeBrun
+# IFS, University of Texas at Austin
+# 20 Mar 1994
+
+# Fundamental settings.
+#
+# To get smoother color variation, increase ns, nx, and ny.  To get faster
+# response (especially on a serial link), decrease them.  A decent but quick
+# test results from ns around 5 and nx, ny around 25.
+
+# Defaults for user options.  Use global variables to track user options since
+# that matches the C example most closely.
+
+$ns      = 20      # Default number of shade levels
+$nx      = 35      # Default number of data points in x
+$ny      = 46      # Default number of data points in y
+$exclude = false   # By default do not plot a page illustrating
+                   # exclusion.  API is probably going to change
+                   # anyway, but can only be reproduced by
+                   # the C and Ruby (and ocaml?) front ends
+
+# polar plot data
+PERIMETERPTS = 100
+
+# Transformation function (modified below)
+
+TR = [
+  1.0, 0.0, -1.0,
+  0.0, 1.0, -1.0
+]
+
+def mypltr(x, y, tr)
+    tx = tr[0] * x.to_f + tr[1] * y.to_f + tr[2]
+    ty = tr[3] * x.to_f + tr[4] * y.to_f + tr[5]
+    [tx, ty]
+end
+
+PLOP = PLOptionParser.new do |op|
+  op.separator('')
+  op.separator('x16 options:')
+  op.on('--exclude', "Plot the \"exclusion\" page [#{$exclude}]") do
+    $exclude = true
+  end
+  op.on('--ns LEVELS', Integer, "Sets number of shade levels [#{$ns}]") do |o|
+    $ns = o
+  end
+  op.on('--nx XPTS', Integer, "Sets number of data points in x [#{$nx}]") do |o|
+    $nx = o
+  end
+  op.on('--ny YPTS', Integer, "Sets number of data points in y [#{$ny}]") do |o|
+    $nx = o
+  end
+  op.separator('')
+  op.separator('To get smoother color variation, increase ns, nx, and ny.  To get faster')
+  op.separator('response (especially on a serial link), decrease them.  A decent but quick')
+  op.separator('test results from ns around 5 and nx, ny around 25.')
+end
+
+zdefined = lambda do |x, y|
+    z = sqrt(x * x + y * y)
+    z < 0.4 || z > 0.6
+end
+
+#--------------------------------------------------------------------------
+# Does several shade plots using different coordinate mappings.
+#--------------------------------------------------------------------------
+
+fill_width = 2
+cont_color = 0
+cont_width = 0
+
+# Parse and process command line arguments
+
+PLOP.parse!
+
+# Load colour palettes
+
+plspal0("cmap0_black_on_white.pal")
+plspal1("cmap1_gray.pal", 1)
+
+# Reduce colors in cmap 0 so that cmap 1 is useful on a 16-color display
+
+plscmap0n(3)
+
+# Initialize plplot
+
+plinit
+
+# Allocate data structures
+
+TR[0] = 2.0 / ($nx - 1)
+TR[4] = 2.0 / ($ny - 1)
+
+clevel = NArray.float($ns)
+shedge = NArray.float($ns+1)
+xg1    = NArray.float($nx)
+yg1    = NArray.float($ny)
+
+xg2 = NArray.float($ny, $nx)
+yg2 = NArray.float($ny, $nx)
+
+z = NArray.float($ny, $nx)
+w = NArray.float($ny, $nx)
+
+px = NArray.float(PERIMETERPTS)
+py = NArray.float(PERIMETERPTS)
+
+# Set up data array
+
+$nx.times do |i|
+  x = (i - ($nx / 2.0)) / ($nx / 2.0)
+
+  $ny.times do |j|
+    y =  (j - ($ny / 2.0)) / ($ny / 2.0) - 1
+
+    z[j,i] = -sin(7.0 * x) * cos(7.0 * y) + x * x - y * y
+    w[j,i] = -cos(7.0 * x) * sin(7.0 * y) + 2 * x * y
+  end
+end
+
+zmin, zmax = plMinMax2dGrid(z)
+
+$ns.times do |i|
+  clevel[i] = zmin + (zmax - zmin) * (i + 0.5) / $ns
+end
+
+($ns+1).times do |i|
+  shedge[i] = zmin + (zmax - zmin) * i / $ns
+end
+
+# Set up coordinate grids
+
+$nx.times do |i|
+  $ny.times do |j|
+    x, y = mypltr(i, j, TR)
+
+    argx    = x * PI / 2
+    argy    = y * PI / 2
+    distort = 0.4
+
+    xg1[i] = x + distort * cos(argx)
+    yg1[j] = y - distort * cos(argy)
+
+    xg2[j,i] = x + distort * cos(argx) * cos(argy)
+    yg2[j,i] = y - distort * cos(argx) * cos(argy)
+  end
+end
+
+# Plot using identity transform
+
+pladv(0)
+plvpor(0.1, 0.9, 0.1, 0.9)
+plwind(-1.0, 1.0, -1.0, 1.0)
+
+plpsty(0)
+
+plshade(z, -1..1, -1..1,
+        shedge, fill_width,
+        cont_color, cont_width,
+        1)
+
+plcol0(1)
+plbox("bcnst", 0.0, 0, "bcnstv", 0.0, 0)
+plcol0(2)
+#plcont(w, nil, nil, clevel) {|x,y| mypltr(x,y,TR)}
+
+pllab("distance", "altitude", "Bogon density")
+
+# Plot using 1d coordinate transform
+
+# Load colour palettes
+plspal0("cmap0_black_on_white.pal")
+plspal1("cmap1_blue_yellow.pal", 1)
+# Reduce colors in cmap 0 so that cmap 1 is useful on a 16-color display
+plscmap0n(3)
+
+pladv(0)
+plvpor(0.1, 0.9, 0.1, 0.9)
+plwind(-1.0, 1.0, -1.0, 1.0)
+
+plpsty(0)
+
+plshade(z, -1..1, -1..1,
+        shedge, fill_width,
+        cont_color, cont_width,
+        1, xg1, yg1)
+
+plcol0(1)
+plbox("bcnst", 0.0, 0, "bcnstv", 0.0, 0)
+plcol0(2)
+
+#plcont(w, nil, nil, clevel, xg1, xg2)
+
+pllab("distance", "altitude", "Bogon density")
+
+# Plot using 2d coordinate transform
+
+# Load colour palettes
+plspal0("cmap0_black_on_white.pal")
+plspal1("cmap1_blue_red.pal", 1)
+# Reduce colors in cmap 0 so that cmap 1 is useful on a 16-color display 
+plscmap0n(3)
+
+pladv(0)
+plvpor(0.1, 0.9, 0.1, 0.9)
+plwind(-1.0, 1.0, -1.0, 1.0)
+
+plpsty(0)
+
+plshade(z, -1..1, -1..1,
+        shedge, fill_width,
+        cont_color, cont_width,
+        0, xg2, yg2)
+
+plcol0(1)
+plbox("bcnst", 0.0, 0, "bcnstv", 0.0, 0)
+plcol0(2)
+plcont(w, nil, nil, clevel, xg2, yg2)
+
+pllab("distance", "altitude", "Bogon density, with streamlines")
+
+# Plot using 2d coordinate transform 
+
+# Load colour palettes
+plspal0(nil)
+plspal1(nil, 1)
+# Reduce colors in cmap 0 so that cmap 1 is useful on a 16-color display 
+plscmap0n(3)
+
+pladv(0)
+plvpor(0.1, 0.9, 0.1, 0.9)
+plwind(-1.0, 1.0, -1.0, 1.0)
+
+plpsty(0)
+
+plshade(z, -1..1, -1..1,
+        shedge, fill_width,
+        2, 3,
+        0, xg2, yg2)
+
+plcol0(1)
+plbox("bcnst", 0.0, 0, "bcnstv", 0.0, 0)
+plcol0(2)
+#plcont(w, nil, nil, clevel, xg2, yg2)
+
+pllab("distance", "altitude", "Bogon density")
+
+# Note this exclusion API will probably change. 
+
+# Plot using 2d coordinate transform and exclusion
+
+if $exclude
+  # Load colour palettes
+  plspal0("cmap0_black_on_white.pal")
+  plspal1("cmap1_gray.pal", 1)
+  # Reduce colors in cmap 0 so that cmap 1 is useful on a 16-color display 
+  plscmap0n(3)
+
+  pladv(0)
+  plvpor(0.1, 0.9, 0.1, 0.9)
+  plwind(-1.0, 1.0, -1.0, 1.0)
+
+  plpsty(0)
+
+  pldefined(zdefined)
+
+  plshade(z, -1..1, -1..1,
+           shedge, fill_width,
+           cont_color, cont_width,
+           0, xg2, yg2)
+
+  pldefined(nil)
+
+  plcol0(1)
+  plbox("bcnst", 0.0, 0, "bcnstv", 0.0, 0)
+
+  pllab("distance", "altitude", "Bogon density with exclusion")
+end
+
+# Example with polar coordinates. 
+
+# Load colour palettes
+plspal0("cmap0_black_on_white.pal")
+plspal1("cmap1_gray.pal", 1)
+# Reduce colors in cmap 0 so that cmap 1 is useful on a 16-color display 
+plscmap0n(3)
+
+pladv(0)
+plvpor(0.1, 0.9, 0.1, 0.9)
+plwind(-1, 1, -1, 1)
+
+
+plpsty(0)
+
+# Build new coordinate matrices. 
+
+$nx.times do |i|
+  r = i.to_f / ($nx - 1)
+  $ny.times do |j|
+    t = 2 * PI / ($ny - 1) * j
+    xg2[j,i] = r * cos(t)
+    yg2[j,i] = r * sin(t)
+    z[j,i]   = exp(-r * r) * cos(5 * PI * r) * cos(5 * t)
+  end
+end
+
+# Need a new shedge to go along with the new data set. 
+
+zmin, zmax = plMinMax2dGrid(z)
+
+($ns+1).times do |i|
+  shedge[i] = zmin + (zmax - zmin) * i / $ns
+end
+
+#  Now we can shade the interior region. 
+
+plshade(z, -1..1, -1..1,
+        shedge, fill_width,
+        cont_color, cont_width,
+        0, xg2, yg2)
+
+# Now we can draw the perimeter.  (If do before, shade stuff may overlap.) 
+
+PERIMETERPTS.times do |i|
+  t     = 2 * PI / (PERIMETERPTS - 1) * i
+  px[i] = cos(t)
+  py[i] = sin(t)
+end
+plcol0(1)
+plline(px, py)
+
+# And label the plot.
+
+plcol0(2)
+pllab("", "", "Tokamak Bogon Instability")
+
+# Clean up 
+
+plend