plplot 0.0.0

data/examples/x03.rb

@@ -0,0 +1,81 @@
+#!/usr/bin/env ruby
+$-w = true if $0 == __FILE__
+
+require 'rubygems'
+require 'plplot'
+include PLplot
+include NMath
+
+# Polar plot demo.
+#
+# Generates polar plot, with 1-1 scaling.
+
+n = 361
+
+dtr = PI / 180.0
+x0 = cos(NArray.float(n).indgen!.mul!(dtr))
+y0 = sin(NArray.float(n).indgen!.mul!(dtr))
+r  = sin(NArray.float(n).indgen!.mul!(5*dtr))
+
+# Parse and process command line arguments
+
+PLOptionParser.parse!
+
+# Set orientation to landscape - note not all device drivers 
+# support this, in particular most interactive drivers do not
+plsori(1)
+
+# Initialize plplot
+
+plinit
+
+# Set up viewport and window, but do not draw box
+
+plenv(-1.3, 1.3, -1.3, 1.3, 1, -2)
+
+# Draw circles for polar grid
+#
+1.upto(10) {|i| plarc(0.0, 0.0, 0.1 * i, 0.1 * i, 0.0, 360.0, 0)}
+
+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.to_s
+
+  #  Write labels for angle
+
+  if (theta < 9.99)
+    offset = 0.45
+  elsif (theta < 99.9)
+    offset = 0.30
+  else
+    offset = 0.15
+  end
+
+  #  Slightly off zero to avoid floating point logic flips at 90 and 270 deg.
+  if (dx >= -0.00001)
+    plptex(dx, dy, dx, dy, -offset, text)
+  else
+    plptex(dx, dy, -dx, -dy, offset+1, text)
+  end
+end
+
+#  Draw the graph
+
+x = r * x0
+y = r * y0
+plcol0(3)
+plline(x, y)
+
+plcol0(4)
+plmtex("t", 2.0, 0.5, 0.5, "#frPLplot Example 3 - r(#gh)=sin 5#gh")
+
+#  Close the plot at end
+
+plend

data/examples/x04.rb

@@ -0,0 +1,98 @@
+#!/usr/bin/env ruby
+$-w = true if $0 == __FILE__
+
+require 'rubygems'
+require 'plplot'
+include PLplot
+include NMath
+
+# Log plot demo.
+#
+# Illustration of logarithmic axes, and redefinition of window.
+
+# Makes transliteration from C example easier
+def pow(a,b)
+  a ** b
+end
+
+#--------------------------------------------------------------------------*\
+# plot1
+#
+# Log-linear plot.
+#--------------------------------------------------------------------------*/
+
+def plot1(type)
+
+  n = 101
+  freql = NArray.float(n)
+  ampl  = NArray.float(n)
+  phase = NArray.float(n)
+
+  pladv(0)
+
+  #  Set up data for log plot
+
+  f0 = 1.0
+  n.times do |i|
+    freql[i] = -2.0 + i / 20.0
+    freq = pow(10.0, freql[i])
+    ampl[i] = 20.0 * log10(1.0 / sqrt(1.0 + pow((freq / f0), 2.0)))
+    phase[i] = -(180.0 / PI) * atan(freq / f0)
+  end
+
+  plvpor(0.15, 0.85, 0.1, 0.9)
+  plwind(-2.0, 3.0, -80.0, 0.0)
+
+  #  Try different axis and labelling styles.
+
+  plcol0(1)
+  case (type)
+  when 0
+    plbox("bclnst", 0.0, 0, "bnstv", 0.0, 0)
+  when 1
+    plbox("bcfghlnst", 0.0, 0, "bcghnstv", 0.0, 0)
+  end
+
+  #  Plot ampl vs freq
+
+  plcol0(2)
+  plline(freql, ampl)
+  plcol0(1)
+  plptex(1.6, -30.0, 1.0, -20.0, 0.5, "-20 dB/decade")
+
+  #  Put labels on
+
+  plcol0(1)
+  plmtex("b", 3.2, 0.5, 0.5, "Frequency")
+  plmtex("t", 2.0, 0.5, 0.5, "Single Pole Low-Pass Filter")
+  plcol0(2)
+  plmtex("l", 5.0, 0.5, 0.5, "Amplitude (dB)")
+
+  #  For the gridless case, put phase vs freq on same plot
+
+  if (type == 0)
+    plcol0(1)
+    plwind(-2.0, 3.0, -100.0, 0.0)
+    plbox("", 0.0, 0, "cmstv", 30.0, 3)
+    plcol0(3)
+    plline(freql, phase)
+    plcol0(3)
+    plmtex("r", 5.0, 0.5, 0.5, "Phase shift (degrees)")
+  end
+end
+
+# Parse and process command line arguments
+
+PLOptionParser.parse!
+
+#  Initialize plplot
+
+plinit
+plfont(2)
+
+#  Make log plots using two different styles.
+
+plot1(0)
+plot1(1)
+
+plend

data/examples/x05.rb

@@ -0,0 +1,33 @@
+#!/usr/bin/env ruby
+$-w = true if $0 == __FILE__
+
+require 'plplot'
+include PLplot
+include NMath
+
+# Histogram demo.
+#
+# Draws a histogram from sample data.
+
+NPTS = 2047
+
+# Parse and process command line arguments
+
+PLOptionParser.parse!
+
+# Initialize plplot
+
+plinit
+
+# Fill up data points
+
+delta = 2.0 * PI / NPTS
+data = sin(NArray.float(NPTS).indgen!.mul!(delta))
+
+plcol0(1)
+plhist(data, 44, -1.1 .. 1.1)
+plcol0(2)
+pllab("#frValue", "#frFrequency",
+      "#frPLplot Example 5 - Probability function of Oscillator")
+
+plend

data/examples/x06.rb

@@ -0,0 +1,62 @@
+#!/usr/bin/env ruby
+$-w = true if $0 == __FILE__
+
+require 'rubygems'
+require 'plplot'
+include PLplot
+
+# --------------------------------------------------------------------------
+# Font demo.
+#
+# Displays the entire "plpoin" symbol (font) set.
+# --------------------------------------------------------------------------
+
+# Parse and process command line arguments
+
+PLOptionParser.parse!
+
+#  Initialize plplot
+
+plinit
+
+pladv(0)
+
+#  Set up viewport and window
+
+plcol0(2)
+plvpor(0.1, 1.0, 0.1, 0.9)
+plwind(0.0, 1.0, 0.0, 1.3)
+
+#  Draw the grid using plbox
+
+plbox("bcg", 0.1, 0, "bcg", 0.1, 0)
+
+#  Write the digits below the frame
+
+plcol0(15)
+10.times do |i|
+  plmtex("b", 1.5, (0.1 * i + 0.05), 0.5, i)
+end
+
+k = 0
+13.times do |i|
+
+  #  Write the digits to the left of the frame
+
+  plmtex("lv", 1.0, (1.0 - (2 * i + 1) / 26.0), 1.0, 10*i)
+
+  10.times do |j|
+    x = 0.1 * j + 0.05
+    y = 1.25 - 0.1 * i
+
+    #  Display the symbols (plpoin accepts scalars as well as NArrays
+
+    if (k < 128)
+      plpoin(x, y, k)
+    end
+    k = k + 1
+  end
+end
+
+plmtex("t", 1.5, 0.5, 0.5, "PLplot Example 6 - plpoin symbols")
+plend

data/examples/x07.rb

@@ -0,0 +1,69 @@
+#!/usr/bin/env ruby
+$-w = true if $0 == __FILE__
+
+require 'rubygems'
+require 'plplot'
+include PLplot
+include NMath
+
+# --------------------------------------------------------------------------*\
+#  Font demo.
+# 
+#  Displays the entire "plsym" symbol (font) set.
+# --------------------------------------------------------------------------*/
+
+base = [
+     0,  200,  500,  600,  700,  800,  900,
+  2000, 2100, 2200, 2300, 2400, 2500, 2600, 2700, 2800, 2900
+]
+
+# Parse and process command line arguments
+
+PLOptionParser.parse!
+
+#  Initialize plplot
+
+plinit
+
+plfontld(1)
+
+17.times do |l|
+  pladv(0)
+
+  #  Set up viewport and window
+
+  plcol0(2)
+  plvpor(0.15, 0.95, 0.1, 0.9)
+  plwind(0.0, 1.0, 0.0, 1.0)
+
+  #  Draw the grid using plbox
+
+  plbox("bcg", 0.1, 0, "bcg", 0.1, 0)
+
+  #  Write the digits below the frame
+
+  plcol0(15)
+  10.times do |i|
+    plmtex("b", 1.5, (0.1 * i + 0.05), 0.5, i)
+  end
+
+  k = 0
+  10.times do |i|
+
+    #  Write the digits to the left of the frame
+
+    plmtex("lv", 1.0, (0.95 - 0.1 * i), 1.0, base[l] + 10 * i)
+    10.times do |j|
+      x = 0.1 * j + 0.05
+      y = 0.95 - 0.1 * i
+
+      #  Display the symbols
+
+      plsym(x, y, base[l] + k)
+      k = k + 1
+    end
+  end
+
+  plmtex("t", 1.5, 0.5, 0.5, "PLplot Example 7 - PLSYM symbols")
+end
+plend

data/examples/x08.rb

@@ -0,0 +1,179 @@
+#!/usr/bin/env ruby
+$-w = true if $0 == __FILE__
+
+require 'rubygems'
+require 'plplot'
+include PLplot
+include NMath
+
+# 3-d plot demo.
+#
+# Copyright (C) 2004  Alan W. Irwin
+# Copyright (C) 2004  Rafael Laboissiere
+#
+# 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
+
+XPTS = 35               #  Data points in x
+YPTS = 46               #  Data points in y
+
+alt = [60.0, 20.0]
+az  = [30.0, 60.0]
+
+title = [
+    "#frPLplot Example 8 - Alt=60, Az=30",
+    "#frPLplot Example 8 - Alt=20, Az=60",
+]
+
+# Makes transliteration from C example easier
+def pow(a,b)
+  a ** b
+end
+
+#--------------------------------------------------------------------------
+# cmap1_init1
+#
+# Initializes color map 1 in HLS space.
+# Basic grayscale variation from half-dark (which makes more interesting
+# looking plot compared to dark) to light.
+# An interesting variation on this:
+#      s[1] = 1.0
+#--------------------------------------------------------------------------
+
+def cmap1_init(gray)
+
+  i = [0.0, 1.0];       #  [left, right] boundaries
+
+  if (gray==1)
+    h = [0.0, 0.0]     #  hue -- low: red (arbitrary if s=0)
+    l = [0.5, 1.0]     #  lightness -- low: half-dark
+    s = [0.0, 0.0]     #  minimum saturation
+  else
+    h = [240, 0]       #  blue -> green -> yellow -> red
+    l = [0.6, 0.6]
+    s = [0.8, 0.8]
+  end
+
+  plscmap1n(256)
+  plscmap1l(PL::CMAP_HLS, i, h, l, s)
+end
+
+# Use PLOptionParser to parse options.  Use global variables to track user
+# options since that matches the C example most closely.
+
+# Defaults for user options
+$sombrero = false
+
+PLOP = PLOptionParser.new do |op|
+  op.separator('')
+  op.separator('x08 options:')
+  op.on('--sombrero', 'Use the "sombrero" function') do
+    $sombrero = true
+  end
+end
+
+LEVELS = 10
+
+#--------------------------------------------------------------------------
+# main
+#
+# Does a series of 3-d plots for a given data set, with different
+# viewing options in each plot.
+#--------------------------------------------------------------------------
+
+PLOP.parse!
+rosen = !$sombrero
+
+plinit
+
+#  Allocate data structures
+
+x = NArray.float(XPTS).indgen!.div!(XPTS/2.0).sbt!(1)
+y = NArray.float(YPTS).indgen!.div!(YPTS/2.0).sbt!(1)
+z = NArray.float(YPTS, XPTS)
+
+if(rosen)
+  x.mul!(1.5)
+  y.add!(0.5)
+end
+
+XPTS.times do |i|
+  xx = x[i]
+  YPTS.times do |j|
+    yy = y[j]
+    if (rosen)
+      z[j,i] = pow(1 - xx, 2) + 100 * pow(yy - pow(xx, 2), 2)
+      #  The log argument may be zero for just the right grid. 
+      if (z[j,i] > 0)
+        z[j,i] = log(z[j,i])
+      else
+        z[j,i] = -5; #  -MAXFLOAT would mess-up up the scale
+      end
+    else
+      r = sqrt(xx * xx + yy * yy)
+      z[j,i] = exp(-r * r) * cos(2.0 * PI * r)
+    end
+  end
+end
+
+# It would be convenient if narray had a minmax function
+zmin = z.min
+zmax = z.max
+
+nlevel = LEVELS
+step = (zmax-zmin)/(nlevel+1)
+clevel = NArray.float(nlevel).indgen!(1).mul!(step).add!(zmin)
+
+pllightsource(1,1,1)
+
+2.times do |k|
+  4.times do |ifshade|
+    pladv(0)
+    plvpor(0.0, 1.0, 0.0, 0.9)
+    plwind(-1.0, 1.0, -0.9, 1.1)
+    plcol0(3)
+    plmtex("t", 1.0, 0.5, 0.5, title[k])
+    plcol0(1)
+    if (rosen)
+      plw3d(1.0, 1.0, 1.0, -1.5, 1.5, -0.5, 1.5, zmin, zmax, alt[k], az[k])
+    else
+      plw3d(1.0, 1.0, 1.0, -1.0, 1.0, -1.0, 1.0, zmin, zmax, alt[k], az[k])
+    end
+
+    plbox3("bnstu", "x axis", 0.0, 0,
+           "bnstu", "y axis", 0.0, 0,
+           "bcdmnstuv", "z axis", 0.0, 0)
+    plcol0(2)
+
+    if (ifshade == 0) #  diffuse light surface plot
+      cmap1_init(1)
+      plsurf3d(x, y, z)
+    elsif (ifshade == 1) #  magnitude colored plot
+      cmap1_init(0)
+      plsurf3d(x, y, z, PL::MAG_COLOR)
+    elsif (ifshade == 2) #   magnitude colored plot with faceted squares
+      cmap1_init(0)
+      plsurf3d(x, y, z, PL::MAG_COLOR | PL::FACETED)
+    else                    #  magnitude colored plot with contours
+      cmap1_init(0)
+      plsurf3d(x, y, z, PL::MAG_COLOR | PL::SURF_CONT | PL::BASE_CONT, clevel)
+    end
+  end
+end
+
+#  Clean up
+
+plend