plplot 0.0.0

data/examples/x12.rb

@@ -0,0 +1,65 @@
+#!/usr/bin/env ruby
+$-w = true if $0 == __FILE__
+
+require 'rubygems'
+require 'plplot'
+include PLplot
+
+# Bar chart demo.
+
+def plfbox( x0, y0 )
+    x = [x0, x0, 1+x0, 1+x0]
+    y = [ 0, y0,   y0,    0]
+    plfill( x, y )
+    plcol0( 1 )
+    pllsty( 1 )
+    plline( x, y )
+end
+
+#--------------------------------------------------------------------------
+# Does a simple bar chart, using color fill.  If color fill is
+# unavailable, pattern fill is used instead (automatic).
+#--------------------------------------------------------------------------
+
+#    int          i
+#    char         string[20]
+#    PLFLT        y0[10]
+
+pos   = [ 0.0, 0.25, 0.5, 0.75, 1.0 ]
+red   = [ 0.0, 0.25, 0.5, 1.0, 1.0 ]
+green = [ 1.0, 0.5, 0.5, 0.5, 1.0 ]
+blue  = [ 1.0, 1.0, 0.5, 0.25, 0.0 ]
+
+# Parse and process command line arguments
+
+PLOptionParser.parse!
+
+# Initialize plplot
+
+plinit
+
+pladv( 0 )
+plvsta
+plwind( 1980.0, 1990.0, 0.0, 35.0 )
+plbox( "bc", 1.0, 0, "bcnv", 10.0, 0 )
+plcol0( 2 )
+pllab( "Year", "Widget Sales (millions)", "#frPLplot Example 12" )
+
+y0 = [5, 15, 12, 24, 28, 30, 20, 8, 12, 3]
+
+plscmap1l( PL::CMAP_RGB, pos, red, green, blue)
+
+10.times do |i|
+  #plcol0(i + 1)
+  plcol1( i / 9.0 )
+  plpsty( 0 )
+  plfbox( ( 1980.0 + i ), y0[i] )
+  string = "%.0f" % y0[i]
+  plptex( ( 1980.0 + i + 0.5 ), ( y0[i] + 1.0 ), 1.0, 0.0, 0.5, string )
+  string = "#{1980 + i}"
+  plmtex( "b", 1.0, ( ( i + 1 ) * 0.1 - 0.05 ), 0.5, string )
+end
+
+# Don't forget to call plend to finish off!
+
+plend

data/examples/x13.rb

@@ -0,0 +1,86 @@
+#!/usr/bin/env ruby
+$-w = true if $0 == __FILE__
+
+require 'rubygems'
+require 'plplot'
+include PLplot
+include Math
+
+# Pie chart demo.
+
+text = [
+    "Maurice",
+    "Geoffrey",
+    "Alan",
+    "Rafael",
+    "Vince"
+]
+
+#--------------------------------------------------------------------------
+# Does a simple pie chart.
+#--------------------------------------------------------------------------
+
+#    int          i, j, dthet, theta0, theta1, theta
+#    PLFLT        just, dx, dy
+#    static PLFLT x[500], y[500], per[5]
+
+per = [10.0, 32.0, 12.0, 30.0, 16.0]
+
+# Parse and process command line arguments
+
+PLOptionParser.parse!
+
+# Initialize plplot
+
+plinit
+
+pladv( 0 )
+# Ensure window has aspect ratio of one so circle is
+# plotted as a circle.
+plvasp( 1.0 )
+plwind( 0.0, 10.0, 0.0, 10.0 )
+# plenv(0.0, 10.0, 0.0, 10.0, 1, -2)
+plcol0( 2 )
+# n.b. all theta quantities scaled by 2*PI/500 to be integers to avoid
+# floating point logic problems.
+theta0 = 0
+dthet  = 1
+5.times do |i|
+  x = [5.0]
+  y = [5.0]
+  j = 1
+  # n.b. the theta quantities multiplied by 2*PI/500 afterward so
+  # in fact per is interpreted as a percentage.
+  theta1 = ( theta0 + 5 * per[i] )
+  if ( i == 4 )
+    theta1 = 500
+  end
+  theta0.step(theta1, dthet) do |theta|
+    x[j] = 5 + 3 * cos( ( 2.0 * PI / 500.0 ) * theta )
+    y[j] = 5 + 3 * sin( ( 2.0 * PI / 500.0 ) * theta )
+    j += 1
+  end
+  plcol0( i + 1 )
+  plpsty( ( i + 3 ) % 8 + 1 )
+  plfill( x, y )
+  plcol0( 1 )
+  plline( x, y )
+  just = ( 2.0 * PI / 500.0 ) * ( theta0 + theta1 ) / 2.0
+  dx   = 0.25 * cos( just )
+  dy   = 0.25 * sin( just )
+  if ( ( theta0 + theta1 ) < 250 || ( theta0 + theta1 ) > 750 )
+    just = 0.0
+  else
+    just = 1.0
+  end
+
+  plptex( ( x[j / 2] + dx ), ( y[j / 2] + dy ), 1.0, 0.0, just, text[i] )
+  theta0 = theta1 #- dthet
+end
+plfont( 2 )
+plschr( 0.0, 1.3 )
+plptex( 5.0, 9.0, 1.0, 0.0, 0.5, "Percentage of Sales" )
+
+# Don't forget to call PLEND to finish off!
+
+plend

data/examples/x14.rb

@@ -0,0 +1,391 @@
+#!/usr/bin/env ruby
+$-w = true if $0 == __FILE__
+
+require 'rubygems'
+require 'plplot'
+include PLplot
+include Math
+
+# Demo of multiple stream/window capability (requires Tk or Tcl-DP).
+#
+# Maurice LeBrun
+# IFS, University of Texas at Austin
+#
+# Copyright (C) 2004  Alan W. Irwin
+#
+# This file is part of PLplot.
+#
+# PLplot is free software; you can redistribute it and/or modify
+# it under the terms of the GNU General Library Public License as published
+# by the Free Software Foundation; either version 2 of the License, or
+# (at your option) any later version.
+#
+# PLplot is distributed in the hope that it will be useful,
+# but WITHOUT ANY WARRANTY; without even the implied warranty of
+# MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
+# GNU Library General Public License for more details.
+#
+# You should have received a copy of the GNU Library General Public License
+# along with PLplot; if not, write to the Free Software
+# Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA
+
+#static PLFLT x[101], y[101]
+#static PLFLT xscale, yscale, xoff, yoff, xs[6], ys[6]
+#static PLINT space0 = 0, mark0 = 0, space1 = 1500, mark1 = 1500
+
+# ===============================================================
+
+# Makes transliteration from C example easier
+def pow(a,b)
+  a ** b
+end
+
+def plot1(xscale, yscale, xoff, yoff)
+  x = NArray.float(60)
+  y = NArray.float(60)
+  60.times do |i|
+    x[i] = xoff + xscale * (i + 1) / 60.0
+    y[i] = yoff + yscale * pow(x[i], 2.0)
+  end
+
+  xmin = x[0]
+  xmax = x[59]
+  ymin = y[0]
+  ymax = y[59]
+
+  xs = []
+  ys = []
+  6.times do |i|
+    xs[i] = x[i * 10 + 3]
+    ys[i] = y[i * 10 + 3]
+  end
+
+  # Set up the viewport and window using PLENV. The range in X is
+  # 0.0 to 6.0, and the range in Y is 0.0 to 30.0. The axes are
+  # scaled separately (just = 0), and we just draw a labelled
+  # box (axis = 0).
+
+  plcol0(1)
+  plenv(xmin, xmax, ymin, ymax, 0, 0)
+  plcol0(6)
+  pllab("(x)", "(y)", "#frPLplot Example 1 - y=x#u2")
+
+  # Plot the data points
+
+  plcol0(9)
+  plpoin(xs, ys, 9)
+
+  # Draw the line through the data
+
+  plcol0(4)
+  plline(x, y)
+  plflush
+end
+
+
+# ===============================================================
+
+def plot2
+  # Set up the viewport and window using PLENV. The range in X is -2.0 to
+  # 10.0, and the range in Y is -0.4 to 2.0. The axes are scaled separately
+  # (just = 0), and we draw a box with axes (axis = 1).
+
+  plcol0(1)
+  plenv(-2.0, 10.0, -0.4, 1.2, 0, 1)
+  plcol0(2)
+  pllab("(x)", "sin(x)/x", "#frPLplot Example 1 - Sinc Function")
+
+  # Fill up the arrays
+
+  x = NArray.float(100)
+  y = NArray.float(100)
+  100.times do |i|
+    x[i] = (i - 19.0) / 6.0
+    y[i] = 1.0
+    if (x[i] != 0.0)
+      y[i] = sin(x[i]) / x[i]
+    end
+  end
+
+  # Draw the line
+
+  plcol0(3)
+  plline(x, y)
+  plflush
+end
+
+# ===============================================================
+
+def plot3
+  # For the final graph we wish to override the default tick intervals, and
+  # so do not use PLENV
+
+  pladv(0)
+
+  # Use standard viewport, and define X range from 0 to 360 degrees, Y range
+  # from -1.2 to 1.2.
+
+  plvsta
+  plwind(0.0, 360.0, -1.2, 1.2)
+
+  # Draw a box with ticks spaced 60 degrees apart in X, and 0.2 in Y.
+
+  plcol0(1)
+  plbox("bcnst", 60.0, 2, "bcnstv", 0.2, 2)
+
+  # Superimpose a dashed line grid, with 1.5 mm marks and spaces. plstyl
+  # expects a pointer!!
+
+  plstyl(1500, 1500)
+  plcol0(2)
+  plbox("g", 30.0, 0, "g", 0.2, 0)
+  plstyl
+
+  plcol0(3)
+  pllab("Angle (degrees)", "sine", "#frPLplot Example 1 - Sine function")
+
+  x = NArray.float(101)
+  y = NArray.float(101)
+  101.times do |i|
+    x[i] = 3.6 * i
+    y[i] = sin(x[i] * PI / 180.0)
+  end
+
+  plcol0(4)
+  plline(x, y)
+  plflush
+end
+
+# ===============================================================
+
+def plot4
+  x0 = NArray.float(361)
+  y0 = NArray.float(361)
+  x  = NArray.float(361)
+  y  = NArray.float(361)
+  dtr = PI / 180.0
+  361.times do |i|
+    x0[i] = cos(dtr * i)
+    y0[i] = sin(dtr * i)
+  end
+
+  # Set up viewport and window, but do not draw box
+
+  plenv(-1.3, 1.3, -1.3, 1.3, 1, -2)
+  (1..10).each do |i|
+    361.times do |j|
+      x[j] = 0.1 * i * x0[j]
+      y[j] = 0.1 * i * y0[j]
+    end
+
+    # Draw circles for polar grid
+
+    plline(x, y)
+  end
+
+  plcol0(2)
+  12.times do |i|
+    theta = 30.0 * i
+    dx    = cos(dtr * theta)
+    dy    = sin(dtr * theta)
+
+    # Draw radial spokes for polar grid
+
+    pljoin(0.0, 0.0, dx, dy)
+    text = theta.round
+
+    # Write labels for angle
+
+    # Slightly off zero to avoid floating point logic flips at 90 and 270 deg.
+    if (dx >= -0.00001)
+      plptex(dx, dy, dx, dy, -0.15, text)
+    else
+      plptex(dx, dy, -dx, -dy, 1.15, text)
+    end
+  end
+
+  # Draw the graph
+
+  361.times do |i|
+    r    = sin(dtr * (5 * i))
+    x[i] = x0[i] * r
+    y[i] = y0[i] * r
+  end
+  plcol0(3)
+  plline(x, y)
+
+  plcol0(4)
+  plmtex("t", 2.0, 0.5, 0.5,
+         "#frPLplot Example 3 - r(#gh)=sin 5#gh")
+  plflush
+end
+
+# ===============================================================
+
+# Demonstration of contour plotting
+
+XPTS = 35
+YPTS = 46
+XSPA = 2.0 / (XPTS - 1.0)
+YSPA = 2.0 / (YPTS - 1.0)
+
+TR = [ XSPA, 0.0, -1.0, 0.0, YSPA, -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
+
+
+CLEVEL = [ -1.0, -0.8, -0.6, -0.4, -0.2, 0.0, 0.2, 0.4, 0.6, 0.8, 1.0 ]
+
+def plot5
+  mark  = 1500
+  space = 1500
+
+  # Set up function arrays
+
+  z = NArray.float(YPTS, XPTS);
+  w = NArray.float(YPTS, XPTS);
+
+  XPTS.times do |i|
+    xx = (i - (XPTS / 2.0)) / (XPTS / 2.0);
+    YPTS.times do |j|
+      yy = (j - (YPTS / 2.0)) / (YPTS / 2.0) - 1.0;
+      z[j,i] = xx * xx - yy * yy;
+      w[j,i] = 2 * xx * yy;
+    end
+  end
+
+  plenv(-1.0, 1.0, -1.0, 1.0, 0, 0)
+  plcol0(2)
+  plcont(z, nil, nil, CLEVEL) {|x,y| mypltr(x, y, TR)}
+  plstyl(mark, space)
+  plcol0(3)
+  plcont(w, nil, nil, CLEVEL) {|x,y| mypltr(x, y, TR)}
+  plcol0(1)
+  pllab("X Coordinate", "Y Coordinate", "Streamlines of flow")
+  plflush
+end
+
+#-------------------------------------------------------------------------
+# Plots several simple functions from other example programs.
+#
+# This version sends the output of the first 4 plots (one page) to two
+# independent streams.
+#--------------------------------------------------------------------------
+
+# Select either TK or DP driver and use a small window
+# Using DP results in a crash at the end due to some odd cleanup problems
+# The geometry strings MUST be in writable memory
+
+geometry_master = "500x410+100+200"
+geometry_slave  = "500x410+650+200"
+
+# Parse and process command line arguments
+
+PLOptionParser.parse!
+
+# If valid geometry specified on command line, use it for both streams.
+xp0, yp0, xleng0, yleng0, xoff0, yoff0 = plgpage
+valid_geometry = (xleng0 > 0 && yleng0 > 0)
+
+# Set up first stream
+
+if (valid_geometry)
+  plspage(xp0, yp0, xleng0, yleng0, xoff0, yoff0)
+else
+  plsetopt("geometry", geometry_master)
+end
+
+plssub(2, 2)
+plinit
+
+driver = plgdev
+fam, num, bmax = plgfam
+
+printf("Demo of multiple output streams via the %s driver.\n", driver)
+printf("Running with the second stream as slave to the first.\n")
+printf("\n")
+
+# Start next stream
+
+plsstrm(1)
+
+if (valid_geometry)
+  plspage(xp0, yp0, xleng0, yleng0, xoff0, yoff0)
+else
+  plsetopt("geometry", geometry_slave)
+end
+
+# Turn off pause to make this a slave (must follow master)
+plspause(0)
+plsdev(driver)
+plsfam(fam, num, bmax)
+# Currently number of digits in format number can only be
+# set via the command line option
+plsetopt("fflen", "2")
+plinit
+
+# Set up the data & plot
+# Original case
+
+plsstrm(0)
+
+xscale = 6.0
+yscale = 1.0
+xoff   = 0.0
+yoff   = 0.0
+plot1(xscale, yscale, xoff, yoff)
+
+# Set up the data & plot
+
+xscale = 1.0
+yscale = 1.0e+6
+plot1(xscale, yscale, xoff, yoff)
+
+# Set up the data & plot
+
+xscale = 1.0
+yscale = 1.0e-6
+digmax = 2
+plsyax(digmax, 0)
+plot1(xscale, yscale, xoff, yoff)
+
+# Set up the data & plot
+
+xscale = 1.0
+yscale = 0.0014
+yoff   = 0.0185
+digmax = 5
+plsyax(digmax, 0)
+plot1(xscale, yscale, xoff, yoff)
+
+# To slave
+# The pleop ensures the eop indicator gets lit.
+
+plsstrm(1)
+plot4
+pleop
+
+# Back to master
+
+plsstrm(0)
+plot2
+plot3
+
+# To slave
+
+plsstrm(1)
+plot5
+pleop
+
+# Back to master to wait for user to advance
+
+plsstrm(0)
+pleop
+
+# Call plend to finish off.
+
+plend