Doc/talks/Rossini-DSC-July2009.lisp

   1 ;;; -*- mode: lisp -*-
   2
   3 ;;; Time-stamp: <2009-07-14 15:09:33 tony>
   4 ;;; Creation:   <2009-03-12 17:14:56 tony>
   5 ;;; File:       plotting-data.lisp
   6 ;;; Author:     AJ Rossini <blindglobe@gmail.com>
   7 ;;; Copyright:  (c)2009--, AJ Rossini.  BSD, LLGPL, or GPLv2, depending
   8 ;;;             on how it arrives.
   9 ;;; Purpose:    Example of generating plots.
  10
  11 ;;; What is this talk of 'release'? Klingons do not make software
  12 ;;; 'releases'.  Our software 'escapes', leaving a bloody trail of
  13 ;;; designers and quality assurance people in its wake.
  14
  15 ;;; This organization and structure is new to the 21st Century
  16 ;;; version..   Think, "21st Century Schizoid Man".
  17
  18 ;; SETUP FOR PLOT EXAMPLE:
  19 (in-package :cl-user)
  20
  21 (defpackage :cl-2d-user-x11
  22   (:use :cl :cl-cairo2  :cl-2d :cl-numlib :cl-colors :bind :cls)
  23   (:shadowing-import-from :lisp-stat call-method call-next-method
  24
  25       expt + - * / ** mod rem abs 1+ 1- log exp sqrt sin cos tan
  26       asin acos atan sinh cosh tanh asinh acosh atanh float random
  27       truncate floor ceiling round minusp zerop plusp evenp oddp
  28       < <= = /= >= > > ;; complex
  29       conjugate realpart imagpart phase
  30       min max logand logior logxor lognot ffloor fceiling
  31       ftruncate fround signum cis
  32
  33       <= float imagpart))
  34
  35 (in-package :cl-2d-user-x11)
  36
  37 ;; PLOT EXAMPLE
  38 #+nil
  39 (progn
  40
  41
  42   ;; this is how you create an X11 frame.  If you supply a
  43   ;; background-color to as-frame, each plot will clear the frame with
  44   ;; this color.
  45
  46   (defparameter *frame1* (as-frame (create-xlib-image-context 300 300)
  47                                   :background-color +white+))
  48
  49   ;; or netbook size, picture is similar but on a lower-res display window.
  50   (defparameter *frame2* (as-frame (create-xlib-image-context 200 200)
  51                                   :background-color +white+))
  52
  53   (plot-function *frame1*
  54                  #'exp (interval-of 0 2)
  55                  :x-title "x"
  56                  :y-title "exp(x)")
  57
  58   ;; split the frame, and you can draw on the subframes independently.
  59   ;; I do this a lot.
  60
  61   (bind ((#2A((f1 f2) (f3 f4))
  62              (split-frame *frame2* (percent 50) (percent 50))))
  63         (defparameter *f1* f1)
  64         (defparameter *f2* f2)
  65         (defparameter *f3* f3)
  66         (defparameter *f4* f4))
  67
  68
  69   (bind ((#2A((f1 f2) (f3 f4))
  70              (split-frame *frame2* (percent 75) (percent 25))))
  71         (defparameter *f1* f1)
  72         (defparameter *f2* f2)
  73         (defparameter *f3* f3)
  74         (defparameter *f4* f4))
  75
  76   (plot-function *f1* #'sin (interval-of 0 2) :x-title "x" :y-title "sin(x)")
  77   (plot-function *f2* #'cos (interval-of 0 2) :x-title "x" :y-title "cos(x)")
  78   (plot-function *f3* #'tan (interval-of 0 2) :x-title "x" :y-title "tan(x)")
  79   (plot-function *f4* #'/ (interval-of 0 2) :x-title "x" :y-title "1/x")
  80
  81   (clear *frame1*)
  82
  83
  84   (let* ((n 500)
  85          (xs (num-sequence :from 0 :to 10 :length n))
  86          (ys (map 'vector #'(lambda (x) (+ x 8 (random 4.0))) xs))
  87          (weights (replicate #'(lambda () (1+ (random 10))) n 'fixnum))
  88          (da (plot-simple *frame1* (interval-of 0 10) (interval-of 10 20)
  89                           :x-title "x" :y-title "y")))
  90     (draw-symbols da xs ys :weights weights))
  91
  92   (xlib-image-context-to-png (context *frame1*) "/home/tony/test1.png")
  93   (xlib-image-context-to-png (context *frame2*) "/home/tony/test2.png")
  94   )
  95
  96
  97 ;;; EXAMPLE FOR DSC2009
  98   (defparameter *frame2* (as-frame (create-xlib-image-context 400 400)
  99                                   :background-color +white+))
 100
 101   (bind ((#2A((f1 f2) (f3 f4))
 102              (split-frame *frame2* (percent 50) (percent 50))))
 103         (defparameter *f1* f1)
 104         (defparameter *f2* f2)
 105         (defparameter *f3* f3)
 106         (defparameter *f4* f4))
 107   (plot-function *f1* #'sin (interval-of 0 2) :x-title "x" :y-title "sin(x)")
 108   (plot-function *f2* #'cos (interval-of 0 2) :x-title "x" :y-title "cos(x)")
 109   (plot-function *f3* #'tan (interval-of 0 2) :x-title "x" :y-title "tan(x)")
 110
 111  (num-sequence :from 0 :to 10 :length 30)
 112
 113   (let* ((n 500)
 114          (xs (num-sequence :from 0 :to 10 :length n))
 115          (ys (map 'vector #'(lambda (x) (+ x 8 (random 4.0))) xs))
 116          (weights (replicate #'(lambda () (1+ (random 10))) n 'fixnum))
 117          (da (plot-simple *f4* (interval-of 0 10) (interval-of 10 20)
 118                           :x-title "x" :y-title "y")))
 119     (draw-symbols da xs ys :weights weights))
 120   (xlib-image-context-to-png (context *f1*) "/home/tony/test1.png")
 121   (xlib-image-context-to-png (context *frame2*) "/home/tony/test2.png")
 122   (destroy (context  *frame2*))
 123
 124