TODO.lisp

   1 ;;; -*- mode: lisp -*-
   2
   3 ;;; Time-stamp: <2009-06-29 21:19:00 tony>
   4 ;;; Creation:   <2008-09-08 08:06:30 tony>
   5 ;;; File:       TODO.lisp
   6 ;;; Author:     AJ Rossini <blindglobe@gmail.com>
   7 ;;; Copyright:  (c) 2007-2008, AJ Rossini <blindglobe@gmail.com>.  BSD.
   8 ;;; Purpose: Stuff that needs to be made working sits inside the
   9 ;;;          progns... This file contains the current challenges to
  10 ;;;          solve, including a description of the setup and the work
  11 ;;;          to solve....
  12
  13 ;;; What is this talk of 'release'? Klingons do not make software
  14 ;;; 'releases'.  Our software 'escapes', leaving a bloody trail of
  15 ;;; designers and quality assurance people in its wake.
  16
  17 ;;; SET UP
  18
  19 (in-package :cl-user)
  20 ;;(asdf:oos 'asdf:load-op 'lisp-matrix)
  21 ;;(asdf:oos 'asdf:compile-op 'lispstat :force t)
  22 ;;(asdf:oos 'asdf:load-op 'lispstat)
  23
  24 (in-package :lisp-stat-unittests)
  25
  26 ;; tests = 80, failures = 8, errors = 15
  27 (run-tests :suite 'lisp-stat-ut)
  28 (describe (run-tests :suite 'lisp-stat-ut))
  29
  30 ;; FIXME: Example: currently not relevant, yet
  31 ;;   (describe (lift::run-test :test-case  'lisp-stat-unittests::create-proto
  32 ;;                             :suite 'lisp-stat-unittests::lisp-stat-ut-proto))
  33
  34 (describe (lift::run-tests :suite 'lisp-stat-ut-dataframe))
  35 (lift::run-tests :suite 'lisp-stat-ut-dataframe)
  36
  37 (describe
  38  (lift::run-test
  39   :test-case  'lisp-stat-unittests::create-proto
  40   :suite 'lisp-stat-unittests::lisp-stat-ut-proto))
  41
  42 (describe 'lisp-stat-ut)
  43
  44 (in-package :ls-user)
  45
  46
  47 #+nil
  48 (progn
  49   (asdf:oos 'asdf:load-op 'xarray)
  50
  51   *features*
  52   *package*
  53
  54   ;; need to consider CLS features:  CFFI, CL-BLAPACK, LISP-MATRIX, XARRAY
  55   ;; CL-PLPLOT,  CL-CAIRO2/CL-2D.
  56
  57   ;; SBCL/CCL -> ...
  58   ;; CLISP/... -> ...
  59   )
  60
  61 #+nil
  62 (progn
  63   ;; Plotting -- need to figure out the core-dump, or change libraries.
  64   ;; (asdf:oos 'asdf:load-op 'cl-plplot)
  65
  66   ;; (defparameter *gdev* "xwin")
  67   (defparameter *gdev* "xcairo")
  68   ;; (cl-plplot::plsdev *gdev*) ; -- usually handled within call.
  69   (plot-ex)
  70   (plot-ex)
  71   ;; Boom! -- there is currently a loose pointer floating around that
  72   ;; causes errors the 3rd time that we create a plot (and crashes
  73   ;; SBCL the 4th time).  Order independent.
  74   (plot-ex)
  75
  76   (contour-plot-ex)
  77   (fn-contour-plot-ex)
  78   (shade-plot-ex)
  79   (3D-plot-ex))
  80
  81
  82 (progn
  83   ;; REVIEW: general Lisp use guidance
  84
  85   (fdefinition 'make-matrix)
  86   (documentation 'make-matrix 'function)
  87
  88 #| Examples from CLHS, a bit of guidance.
  89
  90   ;; This function assumes its callers have checked the types of the
  91   ;; arguments, and authorizes the compiler to build in that assumption.
  92   (defun discriminant (a b c)
  93    (declare (number a b c))
  94    "Compute the discriminant for a quadratic equation."
  95    (- (* b b) (* 4 a c))) =>  DISCRIMINANT
  96   (discriminant 1 2/3 -2) =>  76/9
  97
  98   ;; This function assumes its callers have not checked the types of the
  99   ;; arguments, and performs explicit type checks before making any assumptions.
 100  (defun careful-discriminant (a b c)
 101    "Compute the discriminant for a quadratic equation."
 102    (check-type a number)
 103    (check-type b number)
 104    (check-type c number)
 105    (locally (declare (number a b c))
 106      (- (* b b) (* 4 a c)))) =>  CAREFUL-DISCRIMINANT
 107  (careful-discriminant 1 2/3 -2) =>  76/9
 108 |#
 109   )
 110
 111
 112 #+nil
 113 (progn ;; experiments with GSL and the Lisp interface.
 114   (asdf:oos 'asdf:load-op 'gsll)
 115   (asdf:oos 'asdf:load-op 'gsll-tests) ; requires lisp-unit
 116
 117   ;; the following should be equivalent
 118   (defparameter *t1*  (LIST 6.18d0 6.647777777777779d0 6.18d0))
 119   (defparameter *t2*  (MULTIPLE-VALUE-LIST
 120                (LET ((VEC
 121                       (gsll:make-marray 'DOUBLE-FLOAT
 122                                         :INITIAL-CONTENTS '(-3.21d0 1.0d0 12.8d0)))
 123                      (WEIGHTS
 124                       (gsll:MAKE-MARRAY 'DOUBLE-FLOAT
 125                                         :INITIAL-CONTENTS '(3.0d0 1.0d0 2.0d0))))
 126                  (LET ((MEAN (gsll:MEAN VEC)))
 127                    (LIST (gsll:ABSOLUTE-DEVIATION VEC)
 128                          (gsll:WEIGHTED-ABSOLUTE-DEVIATION VEC WEIGHTS)
 129                          (gsll:ABSOLUTE-DEVIATION VEC MEAN))))))
 130   (eql *t1* *t2*)
 131   (equal *t1* *t2*)
 132
 133   ;; from (gsll:examples 'gsll::numerical-integration) ...
 134   (gsll:integration-qng gsll::one-sine 0.0d0 PI)
 135
 136   (gsll:defun-single axpb (x) (+ (* 2 x) 3)) ;; a<-2, b<-3
 137   (gsll:integration-qng axpb 1d0 2d0)
 138
 139   (let ((a 2)
 140         (b 3))
 141     (defun-single axpb2 (x) (+ (* a x) b)))
 142   (gsll:integration-qng axpb2 1d0 2d0)
 143
 144   ;;   BAD
 145   ;;   (gsll:integration-qng
 146   ;;    (let ((a 2)
 147   ;;     (b 3))
 148   ;;      (defun-single axpb2 (x) (+ (* a x) b)))
 149   ;;    1d0 2d0)
 150
 151   ;; right, but weird expansion...
 152   (gsll:integration-qng
 153    (let ((a 2)
 154          (b 3))
 155      (defun axpb2 (x) (+ (* a x) b))
 156      (gsll:def-single-function axpb2)
 157      axpb2)
 158    1d0 2d0)
 159
 160   ;; Linear least squares
 161
 162   (gsll:gsl-lookup "gsl_linalg_LU_decomp") ; => gsll:lu-decomposition
 163   (gsll:gsl-lookup "gsl_linalg_LU_solve") ; => gsll:lu-solve
 164   )
 165
 166
 167 #+nil
 168 (progn
 169   (asdf:oos 'asdf:load-op 'versioned-objects)
 170   (asdf:oos 'asdf:load-op 'validations)
 171
 172   )
 173
 174
 175 ;; SETUP FOR PLOT EXAMPLE:
 176
 177 (asdf:oos 'asdf:load-op 'lispstat)
 178 (asdf:oos 'asdf:load-op 'cl-cairo2-x11)
 179 (asdf:oos 'asdf:load-op 'cl-2d)
 180
 181 (defpackage :cl-2d-user-x11
 182   (:use :cl :cl-cairo2  :cl-2d :cl-numlib :cl-colors :bind))
 183
 184 (in-package :cl-2d-user-x11)
 185
 186 ;; PLOT EXAMPLE
 187 #+nil
 188 (progn
 189
 190
 191   ;; this is how you create an X11 frame.  If you supply a
 192   ;; background-color to as-frame, each plot will clear the frame with
 193   ;; this color.
 194
 195   (defparameter *frame1* (as-frame (create-xlib-image-context 300 300)
 196                                   :background-color +white+))
 197
 198   ;; or netbook size, picture is similar but on a lower-res display window.
 199   (defparameter *frame2* (as-frame (create-xlib-image-context 200 200)
 200                                   :background-color +white+))
 201
 202   (plot-function *frame1*
 203                  #'exp (interval-of 0 2)
 204                  :x-title "x"
 205                  :y-title "exp(x)")
 206
 207   ;; split the frame, and you can draw on the subframes independently.
 208   ;; I do this a lot.
 209
 210   (bind ((#2A((f1 f2) (f3 f4))
 211              (split-frame *frame2* (percent 50) (percent 50))))
 212         (defparameter *f1* f1)
 213         (defparameter *f2* f2)
 214         (defparameter *f3* f3)
 215         (defparameter *f4* f4))
 216
 217   (plot-function *f1* #'sin (interval-of 0 2) :x-title "x" :y-title "sin(x)")
 218   (plot-function *f2* #'cos (interval-of 0 2) :x-title "x" :y-title "cos(x)")
 219   (plot-function *f3* #'tan (interval-of 0 2) :x-title "x" :y-title "tan(x)")
 220   (plot-function *f4* #'/ (interval-of 0 2) :x-title "x" :y-title "1/x")
 221
 222   (clear *frame1*)
 223
 224
 225   (let* ((n 500)
 226          (xs (num-sequence :from 0 :to 10 :length n))
 227          (ys (map 'vector #'(lambda (x) (+ x 8 (random 4.0))) xs))
 228          (weights (replicate #'(lambda () (1+ (random 10))) n 'fixnum))
 229          (da (plot-simple *frame1* (interval-of 0 10) (interval-of 10 20)
 230                           :x-title "x" :y-title "y")))
 231     (draw-symbols da xs ys :weights weights))
 232
 233   (xlib-image-context-to-png (context *frame1*) "/home/tony/test1.png")
 234   (xlib-image-context-to-png (context *frame2*) "/home/tony/test2.png")
 235   )
 236
 237
 238 ;;; EXAMPLE FOR DSC2009
 239   (defparameter *frame2* (as-frame (create-xlib-image-context 200 200)
 240                                   :background-color +white+))
 241   (bind ((#2A((f1 f2) (f3 f4))
 242              (split-frame *frame2* (percent 50) (percent 50))))
 243         (defparameter *f1* f1)
 244         (defparameter *f2* f2)
 245         (defparameter *f3* f3)
 246         (defparameter *f4* f4))
 247   (plot-function *f1* #'sin (interval-of 0 2) :x-title "x" :y-title "sin(x)")
 248   (plot-function *f2* #'cos (interval-of 0 2) :x-title "x" :y-title "cos(x)")
 249   (plot-function *f3* #'tan (interval-of 0 2) :x-title "x" :y-title "tan(x)")
 250   (let* ((n 500)
 251          (xs (num-sequence :from 0 :to 10 :length n))
 252          (ys (map 'vector #'(lambda (x) (+ x 8 (random 4.0))) xs))
 253          (weights (replicate #'(lambda () (1+ (random 10))) n 'fixnum))
 254          (da (plot-simple *f4* (interval-of 0 10) (interval-of 10 20)
 255                           :x-title "x" :y-title "y")))
 256     (draw-symbols da xs ys :weights weights))
 257   (xlib-image-context-to-png (context *f1*) "/home/tony/test1.png")
 258   (xlib-image-context-to-png (context *frame2*) "/home/tony/test2.png")
 259
 260
 261
 262 ;; back to normal application
 263 (in-package :ls-user)
 264
 265
 266 (defun testme (&key (a 3) (b (+ a 3)))
 267   b)
 268
 269 (testme)
 270 (testme :a 2)
 271 (testme :b 4)
 272 (testme :a 2 :b (* a 5))