Thread: [Plplot-cvs] SF.net SVN: plplot:[8535] trunk

SourceForge Headquarters 225 Broadway Suite 1600 San Diego, CA 92101 +1 (858) 454-5900

Revision: 8535
          http://plplot.svn.sourceforge.net/plplot/?rev=8535&view=rev
Author:   airwin
Date:     2008-07-19 17:32:27 +0000 (Sat, 19 Jul 2008)

Log Message:
-----------
AWI for Hezekiah M. Carty.  Add ocaml standard examples 1, 2, and 3.

Modified Paths:
--------------
    trunk/examples/ocaml/CMakeLists.txt
    trunk/examples/ocaml/Makefile.examples.in
    trunk/plplot_test/test_ocaml.sh.in

Added Paths:
-----------
    trunk/examples/ocaml/x01.ml
    trunk/examples/ocaml/x02.ml
    trunk/examples/ocaml/x03.ml

Modified: trunk/examples/ocaml/CMakeLists.txt
===================================================================

--- trunk/examples/ocaml/CMakeLists.txt	2008-07-19 17:21:37 UTC (rev 8534)
+++ trunk/examples/ocaml/CMakeLists.txt	2008-07-19 17:32:27 UTC (rev 8535)
@@ -21,6 +21,9 @@
 # Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301  USA
 
 set(ocaml_STRING_INDICES
+  "01"
+  "02"
+  "03"
   "11"
   "19"
   )

Modified: trunk/examples/ocaml/Makefile.examples.in
===================================================================
--- trunk/examples/ocaml/Makefile.examples.in	2008-07-19 17:21:37 UTC (rev 8534)
+++ trunk/examples/ocaml/Makefile.examples.in	2008-07-19 17:32:27 UTC (rev 8535)
@@ -28,6 +28,9 @@
 @pkg_config_true@PKG_CONFIG_ENV = @PKG_CONFIG_ENV@
 
 EXECUTABLES_list = \
+		 x01$(EXEEXT) \
+		 x02$(EXEEXT) \
+		 x03$(EXEEXT) \
 		 x11$(EXEEXT) \
 		 x19$(EXEEXT) \
 

Added: trunk/examples/ocaml/x01.ml
===================================================================
--- trunk/examples/ocaml/x01.ml	                        (rev 0)
+++ trunk/examples/ocaml/x01.ml	2008-07-19 17:32:27 UTC (rev 8535)
@@ -0,0 +1,337 @@
+(* $Id$
+
+   Simple line plot and multiple windows demo.  A fairly straightforward
+   translation from the original C to OCaml.
+
+   Copyright (C) 2008  Hezekiah M. Carty
+
+  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
+
+*)
+
+open Plplot
+open Printf
+
+type ex1_params = {
+  xscale : float;
+  yscale : float;
+  xoff : float;
+  yoff : float;
+}
+
+let pi = atan 1.0 *. 4.0
+
+(*
+/* Variables and data arrays used by plot generators */
+
+static PLFLT x[101], y[101];
+static PLFLT xscale, yscale, xoff, yoff, xs[6], ys[6];
+static PLGraphicsIn gin;
+
+static int locate_mode;
+static int test_xor;
+static int fontset = 1;
+static char *f_name = NULL;
+
+/* Options data structure definition. */
+
+static PLOptionTable options[] = {
+{
+    "locate",                   /* Turns on test of API locate function */
+    NULL,
+    NULL,
+    &locate_mode,
+    PL_OPT_BOOL,
+    "-locate",
+    "Turns on test of API locate function" },
+{
+    "xor",                      /* Turns on test of xor function */
+    NULL,
+    NULL,
+    &test_xor,
+    PL_OPT_BOOL,
+    "-xor",
+    "Turns on test of XOR" },
+{
+    "font",                     /* For switching between font set 1 & 2 */
+    NULL,
+    NULL,
+    &fontset,
+    PL_OPT_INT,
+    "-font number",
+    "Selects stroke font set (0 or 1, def:1)" },
+{
+    "save",                     /* For saving in postscript */
+    NULL,
+    NULL,
+    &f_name,
+    PL_OPT_STRING,
+    "-save filename",
+      "Save plot in color postscript `file'" },
+{
+    NULL,                       /* option */
+    NULL,                       /* handler */
+    NULL,                       /* client data */
+    NULL,                       /* address of variable to set */
+    0,                          /* mode flag */
+    NULL,                       /* short syntax */
+    NULL }                      /* long syntax */
+};
+
+const char *notes[] = {"Make sure you get it right!", NULL};
+
+*)
+
+let plot1 (do_test, test_xor) params =
+  let x =
+    Array.init 60 (
+      fun i ->
+        params.xoff +. params.xscale *. (float_of_int i +. 1.0) /. 60.0
+    )
+  in
+  let y =
+    Array.init 60 (
+      fun i ->
+        params.yoff +. params.yscale *. x.(i) ** 2.0
+    )
+  in
+
+  let xmin = x.(0) in
+  let xmax = x.(59) in
+  let ymin = y.(0) in
+  let ymax = y.(59) in
+
+  let xs = Array.init 6 (fun i -> x.(i * 10 + 3)) in
+  let ys = Array.init 6 (fun i -> y.(i * 10 + 3)) in
+
+  (* 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 2;
+  pllab "(x)" "(y)" "#frPLplot Example 1 - y=x#u2";
+
+  (* Plot the data points *)
+  plcol0 4;
+  plpoin xs ys 9;
+
+  (* Draw the line through the data *)
+  plcol0 3;
+  plline x y;
+
+  (* COMMENT THIS SECTION OUT FOR NOW, AS IT DOES NOT FUNCTION EXACTLY
+     AS THE C EXAMPLE.
+  (* xor mode enable erasing a line/point/text by replotting it again *)
+  (* it does not work in double buffering mode, however *)
+  let () =
+    if do_test && test_xor then
+      let st = plxormod 1 in (* enter xor mode *)
+      if st > 0 then begin
+        for i = 1 to 10 do
+          let i' = i * 6 - 6 in
+          plpoin (Array.sub x i' 1) (Array.sub y i' 1) 9;   (* draw a point *)
+          Unix.sleep 1;                                     (* wait a little *)
+          plflush ();                               (* force an update of the driver *)
+          plpoin (Array.sub x i' 1) (Array.sub y i' 1) 9;   (* erase point *)
+        done;
+        ignore (plxormod 0)                                 (* leave xor mode *)
+      end
+      else
+        print_endline "xor mode not supported."
+  in
+  END BIG BLOCK COMMENT *)
+  (* All done. *)
+  ()
+
+let 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 *)
+  let x = Array.init 100 (fun i -> (float_of_int i -. 19.0) /. 6.0) in
+  let y =
+    Array.init 100 (
+      fun i ->
+        if x.(i) <> 0.0 then
+          sin x.(i) /. x.(i)
+        else
+          1.0
+    )
+  in
+
+  (* Draw the line *)
+  plcol0 3;
+  plwid 2;
+  plline x y;
+  plwid 1;
+
+  (* All done. *)
+  ()
+
+let plot3 () =
+  let space1 = 1500 in
+  let mark1 = 1500 in
+
+  (* 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 [|mark1|] [|space1|];
+  plcol0 2;
+  plbox "g" 30.0 0 "g" 0.2 0;
+  plstyl [||] [||];
+
+  plcol0 3;
+  pllab "Angle (degrees)" "sine" "#frPLplot Example 1 - Sine function";
+
+  let x = Array.init 101 (fun i -> 3.6 *. float_of_int i) in
+  let y = Array.init 101 (fun i -> sin (x.(i) *. pi /. 180.0)) in
+
+  plcol0 4;
+  plline x y;
+
+  (* All done. *)
+  ()
+
+(*--------------------------------------------------------------------------*\
+ * main
+ *
+ * Generates several simple line plots.  Demonstrates:
+ *   - subwindow capability
+ *   - setting up the window, drawing plot, and labelling
+ *   - changing the color
+ *   - automatic axis rescaling to exponential notation
+ *   - placing the axes in the middle of the box
+ *   - gridded coordinate axes
+\*--------------------------------------------------------------------------*)
+let main fontset =
+  (* plplot initialization *)
+  (* Divide page into 2x2 plots unless user overrides *)
+  plssub 2 2;
+
+  (* NOT IMPLEMENTED YET
+  (* Parse and process command line arguments *)
+  plMergeOpts(options, "x01c options", notes);
+  plparseopts(&argc, argv, PL_PARSE_FULL);
+  *)
+
+  (* Get version number, just for kicks *)
+  let plplot_version = plgver () in
+  fprintf stdout "PLplot library version: %s\n" plplot_version;
+
+  (* Initialize plplot *)
+  plinit ();
+
+  (* Select font set as per input flag *)
+  plfontld (if fontset then 1 else 0);
+
+  (* Set up the data *)
+  (* Original case *)
+  let plot1_params =
+    {
+      xscale = 6.0;
+      yscale = 1.0;
+      xoff = 0.0;
+      yoff = 0.0;
+    }
+  in
+  (* Do a plot *)
+  plot1 (false, false) plot1_params;
+
+  (* Set up the data *)
+  let plot1_params =
+    {
+      plot1_params with
+        xscale = 1.0;
+        yscale = 0.0014;
+        yoff = 0.0185;
+    }
+  in
+  (* Do a plot *)
+  let digmax = 5 in
+  plsyax digmax 0;
+  plot1 (true, true) plot1_params;
+
+  plot2 ();
+
+  plot3 ();
+
+  (* NOT IMPLEMENTED YET
+  (*
+   * Show how to save a plot:
+   * Open a new device, make it current, copy parameters,
+   * and replay the plot buffer
+   *)
+    if (f_name) { /* command line option '-save filename' */
+
+      printf("The current plot was saved in color Postscript under the name `%s'.\n", f_name);
+      plgstrm(&cur_strm);    /* get current stream */
+      plmkstrm(&new_strm);   /* create a new one */
+
+      plsfnam(f_name);       /* file name */
+      plsdev("psc");         /* device type */
+
+      plcpstrm(cur_strm, 0); /* copy old stream parameters to new stream */
+      plreplot();            /* do the save by replaying the plot buffer */
+      plend1();              /* finish the device */
+
+      plsstrm(cur_strm);     /* return to previous stream */
+    }
+
+  /* Let's get some user input */
+
+    if (locate_mode) {
+        for (;;) {
+          if (! plGetCursor(&gin)) break;
+          if (gin.keysym == PLK_Escape) break;
+
+            pltext();
+            if (gin.keysym < 0xFF && isprint(gin.keysym))
+                printf("subwin = %d, wx = %f,  wy = %f, dx = %f,  dy = %f,  c = '%c'\n",
+                       gin.subwindow, gin.wX, gin.wY, gin.dX, gin.dY, gin.keysym);
+            else
+                printf("subwin = %d, wx = %f,  wy = %f, dx = %f,  dy = %f,  c = 0x%02x\n",
+                       gin.subwindow, gin.wX, gin.wY, gin.dX, gin.dY, gin.keysym);
+
+            plgra();
+        }
+    }
+  *)
+
+  (* Don't forget to call plend() to finish off! *)
+  plend ();
+  ()
+
+let () = main true


Property changes on: trunk/examples/ocaml/x01.ml
___________________________________________________________________
Added: svn:eol-style
   + native

Added: trunk/examples/ocaml/x02.ml
===================================================================
--- trunk/examples/ocaml/x02.ml	                        (rev 0)
+++ trunk/examples/ocaml/x02.ml	2008-07-19 17:32:27 UTC (rev 8535)
@@ -0,0 +1,144 @@
+(* $Id$
+
+        Multiple window and color map 0 demo.
+   Copyright (C) 2008  Hezekiah M. Carty
+
+  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
+*)
+
+open Plplot
+
+(*--------------------------------------------------------------------------*\
+ * draw_windows
+ *
+ * Draws a set of numbered boxes with colors according to cmap0 entry.
+\*--------------------------------------------------------------------------*)
+
+let draw_windows nw cmap0_offset =
+  plschr 0.0 3.5;
+  plfont 4;
+
+  for i = 0 to nw - 1 do
+    plcol0 (i + cmap0_offset);
+    let text = string_of_int i in
+    pladv 0;
+    let vmin = 0.1 in
+    let vmax = 0.9 in
+    for j = 0 to 2 do
+      plwid (j + 1);
+      let j' = float_of_int j *. 0.1 in
+      plvpor (vmin +. j') (vmax -. j') (vmin +. j') (vmax -. j');
+      plwind 0.0 1.0 0.0 1.0;
+      plbox "bc" 0.0 0 "bc" 0.0 0;
+    done;
+    plwid 1;
+    plptex 0.5 0.5 1.0 0.0 0.5 text;
+  done;
+  ()
+
+(*--------------------------------------------------------------------------*\
+ * demo1
+ *
+ * Demonstrates multiple windows and default color map 0 palette.
+\*--------------------------------------------------------------------------*)
+
+let demo1 () =
+  plbop ();
+  (* Divide screen into 16 regions *)
+  plssub 4 4;
+  draw_windows 16 0;
+  pleop ();
+  ()
+
+(*--------------------------------------------------------------------------*\
+ * demo2
+ *
+ * Demonstrates multiple windows, user-modified color map 0 palette, and
+ * HLS -> RGB translation.
+\*--------------------------------------------------------------------------*)
+
+let demo2 () =
+  (* Set up cmap0 *)
+  (* Use 100 custom colors in addition to base 16 *)
+  let r = Array.make 116 0 in
+  let g = Array.make 116 0 in
+  let b = Array.make 116 0 in
+
+  (* Min & max lightness values *)
+  let lmin = 0.15 and lmax = 0.85 in
+
+  plbop ();
+
+  (* Divide screen into 100 regions *)
+  plssub 10 10;
+
+  for i = 0 to 99 do
+    (* Bounds on HLS, from plhlsrgb() commentary --
+     *  hue             [0., 360.]      degrees
+     *  lightness       [0., 1.]        magnitude
+     *  saturation      [0., 1.]        magnitude
+     *)
+    (* Vary hue uniformly from left to right *)
+    let h = (360. /. 10. ) *. float_of_int ( i mod 10 ) in
+    (* Vary lightness uniformly from top to bottom, between min & max *)
+    let l = lmin +. (lmax -. lmin) *. float_of_int (i / 10) /. 9. in
+    (* Use max saturation *)
+    let s = 1.0 in
+
+    let (r1, g1, b1) = plhlsrgb h l s in
+
+    r.(i+16) <- int_of_float (r1 *. 255.);
+    g.(i+16) <- int_of_float (g1 *. 255.);
+    b.(i+16) <- int_of_float (b1 *. 255.);
+  done;
+
+  (* Load default cmap0 colors into our custom set *)
+  for i = 0 to 15 do
+    let (r1, g1, b1) = plgcol0 i in
+    r.(i) <- r1;
+    g.(i) <- g1;
+    b.(i) <- b1;
+  done;
+
+  (* Now set cmap0 all at once (faster, since fewer driver calls) *)
+  plscmap0 r g b;
+  draw_windows 100 16;
+  pleop ();
+  ()
+
+(*--------------------------------------------------------------------------*\
+ * main
+ *
+ * Demonstrates multiple windows and color map 0 palette, both default and
+ * user-modified.
+\*--------------------------------------------------------------------------*)
+
+let main () =
+  (* NOT IMPLEMENTED
+  (* Parse and process command line arguments *)
+    (void) plparseopts(&argc, argv, PL_PARSE_FULL);
+  *)
+
+  (* Initialize plplot *)
+  plinit ();
+  (* Run demos *)
+  demo1 ();
+  demo2 ();
+  plend ();
+  ()
+
+let () = main ()


Property changes on: trunk/examples/ocaml/x02.ml
___________________________________________________________________
Added: svn:eol-style
   + native

Added: trunk/examples/ocaml/x03.ml
===================================================================
--- trunk/examples/ocaml/x03.ml	                        (rev 0)
+++ trunk/examples/ocaml/x03.ml	2008-07-19 17:32:27 UTC (rev 8535)
@@ -0,0 +1,92 @@
+(* $Id$
+
+  Polar plot demo.
+   Copyright (C) 2008  Hezekiah M. Carty
+
+  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
+*)
+
+open Plplot
+
+let pi = atan 1.0 *. 4.0
+
+(*--------------------------------------------------------------------------*\
+ * main
+ *
+ * Generates polar plot, with 1-1 scaling.
+\*--------------------------------------------------------------------------*)
+
+let main () =
+  let dtr = pi /. 180.0 in
+  let x0 = Array.init 361 (fun i -> cos (dtr *. float_of_int i)) in
+  let y0 = Array.init 361 (fun i -> sin (dtr *. float_of_int i)) in
+
+  (* NOT SUPPORTED YET
+  (* Parse and process command line arguments *)
+  (void) plparseopts(&argc, argv, PL_PARSE_FULL);
+  *)
+
+  (* Initialize plplot *)
+  plinit ();
+
+  (* Set up viewport and window, but do not draw box *)
+  plenv (-1.3) 1.3 (-1.3) 1.3 1 (-2);
+  for i = 1 to 10 do
+    let x = Array.init 361 (fun j -> 0.1 *. float_of_int i *. x0.(j)) in
+    let y = Array.init 361 (fun j -> 0.1 *. float_of_int i *. y0.(j)) in
+    (* Draw circles for polar grid *)
+    plline x y;
+  done;
+
+  plcol0 2;
+  for i = 0 to 11 do
+    let theta = 30.0 *. float_of_int i in
+    let dx = cos (dtr *. theta) in
+    let dy = sin (dtr *. theta) in
+
+    (* Draw radial spokes for polar grid *)
+    pljoin 0.0 0.0 dx dy;
+    let text = string_of_int (int_of_float (floor (theta +. 0.5))) in
+
+    (* Write labels for angle *)
+
+    (* Slightly off zero to avoid floating point logic flips at 90 and
+       270 deg.*)
+    if dx >= -0.00001 then
+      plptex dx dy dx dy (-0.15) text
+    else
+      plptex dx dy (-.dx) (-.dy) 1.15 text
+  done;
+
+  (* Draw the graph *)
+  let x =
+    Array.init 361 (fun i -> x0.(i) *. sin (dtr *. float_of_int (5 * i)))
+  in
+  let y =
+    Array.init 361 (fun i -> y0.(i) *. sin (dtr *. float_of_int (5 * i)))
+  in
+  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 ();
+  ()
+
+let () = main ()


Property changes on: trunk/examples/ocaml/x03.ml
___________________________________________________________________
Added: svn:eol-style
   + native

Modified: trunk/plplot_test/test_ocaml.sh.in
===================================================================
--- trunk/plplot_test/test_ocaml.sh.in	2008-07-19 17:21:37 UTC (rev 8534)
+++ trunk/plplot_test/test_ocaml.sh.in	2008-07-19 17:32:27 UTC (rev 8535)
@@ -24,7 +24,7 @@
 
 # Do the standard non-interactive examples.
 # Currently only 11 and 19 are implemented
-for index in 11 19; do
+for index in 01 02 03 11 19; do
   $ocamldir/x${index}ocaml -dev $device -o ${OUTPUT_DIR}/x${index}ocaml.$dsuffix \
     $options 2> test.error
   # Look for any status codes (segfaults, plexit) from the examples themselves.


This was sent by the SourceForge.net collaborative development platform, the world's largest Open Source development site.




Thread: [Plplot-cvs] SF.net SVN: plplot:[8535] trunk

Cross-platform, scientific graphics plotting library

plplot-cvs