lisp/calc/calc-graph.el

   1 ;;; calc-graph.el --- graph output functions for Calc
   2
   3 ;; Copyright (C) 1990, 1991, 1992, 1993, 2001, 2005 Free Software Foundation, Inc.
   4
   5 ;; Author: David Gillespie <daveg@synaptics.com>
   6 ;; Maintainer: Jay Belanger <belanger@truman.edu>
   7
   8 ;; This file is part of GNU Emacs.
   9
  10 ;; GNU Emacs is distributed in the hope that it will be useful,
  11 ;; but WITHOUT ANY WARRANTY.  No author or distributor
  12 ;; accepts responsibility to anyone for the consequences of using it
  13 ;; or for whether it serves any particular purpose or works at all,
  14 ;; unless he says so in writing.  Refer to the GNU Emacs General Public
  15 ;; License for full details.
  16
  17 ;; Everyone is granted permission to copy, modify and redistribute
  18 ;; GNU Emacs, but only under the conditions described in the
  19 ;; GNU Emacs General Public License.   A copy of this license is
  20 ;; supposed to have been given to you along with GNU Emacs so you
  21 ;; can know your rights and responsibilities.  It should be in a
  22 ;; file named COPYING.  Among other things, the copyright notice
  23 ;; and this notice must be preserved on all copies.
  24
  25 ;;; Commentary:
  26
  27 ;;; Code:
  28
  29 ;; This file is autoloaded from calc-ext.el.
  30
  31 (require 'calc-ext)
  32 (require 'calc-macs)
  33
  34 ;;; Graphics
  35
  36 ;; The following three variables are customizable and defined in calc.el.
  37 (defvar calc-gnuplot-name)
  38 (defvar calc-gnuplot-plot-command)
  39 (defvar calc-gnuplot-print-command)
  40
  41 (defvar calc-gnuplot-tempfile "calc")
  42
  43 (defvar calc-gnuplot-default-device)
  44 (defvar calc-gnuplot-default-output)
  45 (defvar calc-gnuplot-print-device)
  46 (defvar calc-gnuplot-print-output)
  47 (defvar calc-gnuplot-keep-outfile nil)
  48 (defvar calc-gnuplot-version nil)
  49
  50 (defvar calc-gnuplot-display (getenv "DISPLAY"))
  51 (defvar calc-gnuplot-geometry)
  52
  53 (defvar calc-graph-default-resolution)
  54 (defvar calc-graph-default-resolution-3d)
  55 (defvar calc-graph-default-precision 5)
  56
  57 (defvar calc-gnuplot-buffer nil)
  58 (defvar calc-gnuplot-input nil)
  59
  60 (defvar calc-gnuplot-last-error-pos 1)
  61 (defvar calc-graph-last-device nil)
  62 (defvar calc-graph-last-output nil)
  63 (defvar calc-graph-file-cache nil)
  64 (defvar calc-graph-var-cache nil)
  65 (defvar calc-graph-data-cache nil)
  66 (defvar calc-graph-data-cache-limit 10)
  67 (defvar calc-graph-no-auto-view nil)
  68 (defvar calc-graph-no-wait nil)
  69 (defvar calc-gnuplot-trail-mark)
  70
  71 (defun calc-graph-fast (many)
  72   (interactive "P")
  73   (let ((calc-graph-no-auto-view t))
  74     (calc-graph-delete t)
  75     (calc-graph-add many)
  76     (calc-graph-plot nil)))
  77
  78 (defun calc-graph-fast-3d (many)
  79   (interactive "P")
  80   (let ((calc-graph-no-auto-view t))
  81     (calc-graph-delete t)
  82     (calc-graph-add-3d many)
  83     (calc-graph-plot nil)))
  84
  85 (defun calc-graph-delete (all)
  86   (interactive "P")
  87   (calc-wrapper
  88    (calc-graph-init)
  89    (save-excursion
  90      (set-buffer calc-gnuplot-input)
  91      (and (calc-graph-find-plot t all)
  92           (progn
  93             (if (looking-at "s?plot")
  94                 (progn
  95                   (setq calc-graph-var-cache nil)
  96                   (delete-region (point) (point-max)))
  97               (delete-region (point) (1- (point-max)))))))
  98    (calc-graph-view-commands)))
  99
 100 (defun calc-graph-find-plot (&optional before all)
 101   (goto-char (point-min))
 102   (and (re-search-forward "^s?plot[ \t]+" nil t)
 103        (let ((beg (point)))
 104          (goto-char (point-max))
 105          (if (or all
 106                  (not (search-backward "," nil t))
 107                  (< (point) beg))
 108              (progn
 109                (goto-char beg)
 110                (if before
 111                    (beginning-of-line)))
 112            (or before
 113                (re-search-forward ",[ \t]+")))
 114          t)))
 115
 116 (defun calc-graph-add (many)
 117   (interactive "P")
 118   (calc-wrapper
 119    (calc-graph-init)
 120    (cond ((null many)
 121           (calc-graph-add-curve (calc-graph-lookup (calc-top-n 2))
 122                                 (calc-graph-lookup (calc-top-n 1))))
 123          ((or (consp many) (eq many 0))
 124           (let ((xdata (calc-graph-lookup (calc-top-n 2)))
 125                 (ylist (calc-top-n 1)))
 126             (or (eq (car-safe ylist) 'vec)
 127                 (error "Y argument must be a vector"))
 128             (while (setq ylist (cdr ylist))
 129               (calc-graph-add-curve xdata (calc-graph-lookup (car ylist))))))
 130          ((> (setq many (prefix-numeric-value many)) 0)
 131           (let ((xdata (calc-graph-lookup (calc-top-n (1+ many)))))
 132             (while (> many 0)
 133               (calc-graph-add-curve xdata
 134                                     (calc-graph-lookup (calc-top-n many)))
 135               (setq many (1- many)))))
 136          (t
 137           (let (pair)
 138             (setq many (- many))
 139             (while (> many 0)
 140               (setq pair (calc-top-n many))
 141               (or (and (eq (car-safe pair) 'vec)
 142                        (= (length pair) 3))
 143                   (error "Argument must be an [x,y] vector"))
 144               (calc-graph-add-curve (calc-graph-lookup (nth 1 pair))
 145                                     (calc-graph-lookup (nth 2 pair)))
 146               (setq many (1- many))))))
 147    (calc-graph-view-commands)))
 148
 149 (defun calc-graph-add-3d (many)
 150   (interactive "P")
 151   (calc-wrapper
 152    (calc-graph-init)
 153    (cond ((null many)
 154           (calc-graph-add-curve (calc-graph-lookup (calc-top-n 3))
 155                                 (calc-graph-lookup (calc-top-n 2))
 156                                 (calc-graph-lookup (calc-top-n 1))))
 157          ((or (consp many) (eq many 0))
 158           (let ((xdata (calc-graph-lookup (calc-top-n 3)))
 159                 (ydata (calc-graph-lookup (calc-top-n 2)))
 160                 (zlist (calc-top-n 1)))
 161             (or (eq (car-safe zlist) 'vec)
 162                 (error "Z argument must be a vector"))
 163             (while (setq zlist (cdr zlist))
 164               (calc-graph-add-curve xdata ydata
 165                                     (calc-graph-lookup (car zlist))))))
 166          ((> (setq many (prefix-numeric-value many)) 0)
 167           (let ((xdata (calc-graph-lookup (calc-top-n (+ many 2))))
 168                 (ydata (calc-graph-lookup (calc-top-n (+ many 1)))))
 169             (while (> many 0)
 170               (calc-graph-add-curve xdata ydata
 171                                     (calc-graph-lookup (calc-top-n many)))
 172               (setq many (1- many)))))
 173          (t
 174           (let (curve)
 175             (setq many (- many))
 176             (while (> many 0)
 177               (setq curve (calc-top-n many))
 178               (or (and (eq (car-safe curve) 'vec)
 179                        (= (length curve) 4))
 180                   (error "Argument must be an [x,y,z] vector"))
 181               (calc-graph-add-curve (calc-graph-lookup (nth 1 curve))
 182                                     (calc-graph-lookup (nth 2 curve))
 183                                     (calc-graph-lookup (nth 3 curve)))
 184               (setq many (1- many))))))
 185    (calc-graph-view-commands)))
 186
 187 (defun calc-graph-add-curve (xdata ydata &optional zdata)
 188   (let ((num (calc-graph-count-curves))
 189         (pstyle (calc-var-value 'var-PointStyles))
 190         (lstyle (calc-var-value 'var-LineStyles)))
 191     (save-excursion
 192       (set-buffer calc-gnuplot-input)
 193       (goto-char (point-min))
 194       (if (re-search-forward (if zdata "^plot[ \t]" "^splot[ \t]")
 195                              nil t)
 196           (error "Can't mix 2d and 3d curves on one graph"))
 197       (if (re-search-forward "^s?plot[ \t]" nil t)
 198           (progn
 199             (end-of-line)
 200             (insert ", "))
 201         (goto-char (point-max))
 202         (or (eq (preceding-char) ?\n)
 203             (insert "\n"))
 204         (insert (if zdata "splot" "plot") " \n")
 205         (forward-char -1))
 206       (insert "{" (symbol-name (nth 1 xdata))
 207               ":" (symbol-name (nth 1 ydata)))
 208       (if zdata
 209           (insert ":" (symbol-name (nth 1 zdata))))
 210       (insert "} "
 211               "title \"" (symbol-name (nth 1 ydata)) "\" "
 212               "with dots")
 213       (setq pstyle (and (eq (car-safe pstyle) 'vec) (nth (1+ num) pstyle)))
 214       (setq lstyle (and (eq (car-safe lstyle) 'vec) (nth (1+ num) lstyle))))
 215     (calc-graph-set-styles
 216      (or (and (Math-num-integerp lstyle) (math-trunc lstyle))
 217          0)
 218      (or (and (Math-num-integerp pstyle) (math-trunc pstyle))
 219          (if (eq (car-safe (calc-var-value (nth 2 ydata))) 'vec)
 220              0 -1)))))
 221
 222 (defun calc-graph-lookup (thing)
 223   (if (and (eq (car-safe thing) 'var)
 224            (calc-var-value (nth 2 thing)))
 225       thing
 226     (let ((found (assoc thing calc-graph-var-cache)))
 227       (or found
 228           (let ((varname (concat "PlotData"
 229                                  (int-to-string
 230                                   (1+ (length calc-graph-var-cache))))))
 231             (setq var (list 'var (intern varname)
 232                             (intern (concat "var-" varname)))
 233                   found (cons thing var)
 234                   calc-graph-var-cache (cons found calc-graph-var-cache))
 235             (set (nth 2 var) thing)))
 236       (cdr found))))
 237
 238 (defun calc-graph-juggle (arg)
 239   (interactive "p")
 240   (calc-graph-init)
 241   (save-excursion
 242     (set-buffer calc-gnuplot-input)
 243     (if (< arg 0)
 244         (let ((num (calc-graph-count-curves)))
 245           (if (> num 0)
 246               (while (< arg 0)
 247                 (setq arg (+ arg num))))))
 248     (while (>= (setq arg (1- arg)) 0)
 249       (calc-graph-do-juggle))))
 250
 251 (defun calc-graph-count-curves ()
 252   (save-excursion
 253     (set-buffer calc-gnuplot-input)
 254     (if (re-search-forward "^s?plot[ \t]" nil t)
 255         (let ((num 1))
 256           (goto-char (point-min))
 257           (while (search-forward "," nil t)
 258             (setq num (1+ num)))
 259           num)
 260       0)))
 261
 262 (defun calc-graph-do-juggle ()
 263   (let (base)
 264     (and (calc-graph-find-plot t t)
 265          (progn
 266            (setq base (point))
 267            (calc-graph-find-plot t nil)
 268            (or (eq base (point))
 269                (let ((str (buffer-substring (+ (point) 2) (1- (point-max)))))
 270                  (delete-region (point) (1- (point-max)))
 271                  (goto-char (+ base 5))
 272                  (insert str ", ")))))))
 273
 274 (defun calc-graph-print (flag)
 275   (interactive "P")
 276   (calc-graph-plot flag t))
 277
 278 (defvar var-DUMMY)
 279 (defvar var-DUMMY2)
 280 (defvar var-PlotRejects)
 281
 282 ;; The following variables are local to calc-graph-plot, but are
 283 ;; used in the functions calc-graph-compute-2d, calc-graph-refine-2d,
 284 ;; calc-graph-recompute-2d, calc-graph-compute-3d and
 285 ;; calc-graph-format-data, which are called by calc-graph-plot.
 286 (defvar calc-graph-yvalue)
 287 (defvar calc-graph-yvec)
 288 (defvar calc-graph-numsteps)
 289 (defvar calc-graph-numsteps3)
 290 (defvar calc-graph-xvalue)
 291 (defvar calc-graph-xvec)
 292 (defvar calc-graph-xname)
 293 (defvar calc-graph-yname)
 294 (defvar calc-graph-xstep)
 295 (defvar calc-graph-ycache)
 296 (defvar calc-graph-ycacheptr)
 297 (defvar calc-graph-refine)
 298 (defvar calc-graph-keep-file)
 299 (defvar calc-graph-xval)
 300 (defvar calc-graph-xlow)
 301 (defvar calc-graph-xhigh)
 302 (defvar calc-graph-yval)
 303 (defvar calc-graph-yp)
 304 (defvar calc-graph-xp)
 305 (defvar calc-graph-zp)
 306 (defvar calc-graph-yvector)
 307 (defvar calc-graph-resolution)
 308 (defvar calc-graph-y3value)
 309 (defvar calc-graph-y3name)
 310 (defvar calc-graph-y3step)
 311 (defvar calc-graph-zval)
 312 (defvar calc-graph-stepcount)
 313 (defvar calc-graph-is-splot)
 314 (defvar calc-graph-surprise-splot)
 315 (defvar calc-graph-blank)
 316 (defvar calc-graph-non-blank)
 317 (defvar calc-graph-curve-num)
 318
 319 (defun calc-graph-plot (flag &optional printing)
 320   (interactive "P")
 321   (calc-slow-wrapper
 322    (let ((calcbuf (current-buffer))
 323          (tempbuf (get-buffer-create "*Gnuplot Temp-2*"))
 324          (tempbuftop 1)
 325          (tempoutfile nil)
 326          (calc-graph-curve-num 0)
 327          (calc-graph-refine (and flag (> (prefix-numeric-value flag) 0)))
 328          (recompute (and flag (< (prefix-numeric-value flag) 0)))
 329          (calc-graph-surprise-splot nil)
 330          (tty-output nil)
 331          cache-env calc-graph-is-splot device output calc-graph-resolution precision samples-pos)
 332      (add-hook 'kill-emacs-hook 'calc-graph-kill-hook)
 333      (save-excursion
 334        (calc-graph-init)
 335        (set-buffer tempbuf)
 336        (erase-buffer)
 337        (set-buffer calc-gnuplot-input)
 338        (goto-char (point-min))
 339        (setq calc-graph-is-splot (re-search-forward "^splot[ \t]" nil t))
 340        (let ((str (buffer-string))
 341              (ver calc-gnuplot-version))
 342          (set-buffer (get-buffer-create "*Gnuplot Temp*"))
 343          (erase-buffer)
 344          (insert "# (Note: This is a temporary copy---do not edit!)\n")
 345          (if (>= ver 2)
 346              (insert "set noarrow\nset nolabel\n"
 347                      "set autoscale xy\nset nologscale xy\n"
 348                      "set xlabel\nset ylabel\nset title\n"
 349                      "set noclip points\nset clip one\nset clip two\n"
 350                      "set format \"%g\"\nset tics\nset xtics\nset ytics\n"
 351                      "set data style linespoints\n"
 352                      "set nogrid\nset nokey\nset nopolar\n"))
 353          (if (>= ver 3)
 354              (insert "set surface\nset nocontour\n"
 355                      "set " (if calc-graph-is-splot "" "no") "parametric\n"
 356                      "set notime\nset border\nset ztics\nset zeroaxis\n"
 357                      "set view 60,30,1,1\nset offsets 0,0,0,0\n"))
 358          (setq samples-pos (point))
 359          (insert "\n\n" str))
 360        (goto-char (point-min))
 361        (if calc-graph-is-splot
 362            (if calc-graph-refine
 363                (error "This option works only for 2d plots")
 364              (setq recompute t)))
 365        (let ((calc-gnuplot-input (current-buffer))
 366              (calc-graph-no-auto-view t))
 367          (if printing
 368              (setq device calc-gnuplot-print-device
 369                    output calc-gnuplot-print-output)
 370            (setq device (calc-graph-find-command "terminal")
 371                  output (calc-graph-find-command "output"))
 372            (or device
 373                (setq device calc-gnuplot-default-device))
 374            (if output
 375                (setq output (car (read-from-string output)))
 376              (setq output calc-gnuplot-default-output)))
 377          (if (or (equal device "") (equal device "default"))
 378              (setq device (if printing
 379                               "postscript"
 380                             (if (or (eq window-system 'x) (getenv "DISPLAY"))
 381                                 "x11"
 382                               (if (>= calc-gnuplot-version 3)
 383                                   "dumb" "postscript")))))
 384          (if (equal device "dumb")
 385              (setq device (format "dumb %d %d"
 386                                   (1- (frame-width)) (1- (frame-height)))))
 387          (if (equal device "big")
 388              (setq device (format "dumb %d %d"
 389                                   (* 4 (- (frame-width) 3))
 390                                   (* 4 (- (frame-height) 3)))))
 391          (if (stringp output)
 392              (if (or (equal output "auto")
 393                      (and (equal output "tty") (setq tty-output t)))
 394                  (setq tempoutfile (calc-temp-file-name -1)
 395                        output tempoutfile))
 396            (setq output (eval output)))
 397          (or (equal device calc-graph-last-device)
 398              (progn
 399                (setq calc-graph-last-device device)
 400                (calc-gnuplot-command "set terminal" device)))
 401          (or (equal output calc-graph-last-output)
 402              (progn
 403                (setq calc-graph-last-output output)
 404                (calc-gnuplot-command "set output"
 405                                      (if (equal output "STDOUT")
 406                                          ""
 407                                        (prin1-to-string output)))))
 408          (setq calc-graph-resolution (calc-graph-find-command "samples"))
 409          (if calc-graph-resolution
 410              (setq calc-graph-resolution (string-to-number calc-graph-resolution))
 411            (setq calc-graph-resolution (if calc-graph-is-splot
 412                                 calc-graph-default-resolution-3d
 413                               calc-graph-default-resolution)))
 414          (setq precision (calc-graph-find-command "precision"))
 415          (if precision
 416              (setq precision (string-to-number precision))
 417            (setq precision calc-graph-default-precision))
 418          (calc-graph-set-command "terminal")
 419          (calc-graph-set-command "output")
 420          (calc-graph-set-command "samples")
 421          (calc-graph-set-command "precision"))
 422        (goto-char samples-pos)
 423        (insert "set samples " (int-to-string (max (if calc-graph-is-splot 20 200)
 424                                                   (+ 5 calc-graph-resolution))) "\n")
 425        (while (re-search-forward "{\\*[^}]+}[^,\n]*" nil t)
 426          (delete-region (match-beginning 0) (match-end 0))
 427          (if (looking-at ",")
 428              (delete-char 1)
 429            (while (memq (preceding-char) '(?\s ?\t))
 430              (forward-char -1))
 431            (if (eq (preceding-char) ?\,)
 432                (delete-backward-char 1))))
 433        (save-excursion
 434          (set-buffer calcbuf)
 435          (setq cache-env (list calc-angle-mode
 436                                calc-complex-mode
 437                                calc-simplify-mode
 438                                calc-infinite-mode
 439                                calc-word-size
 440                                precision calc-graph-is-splot))
 441          (if (and (not recompute)
 442                   (equal (cdr (car calc-graph-data-cache)) cache-env))
 443              (while (> (length calc-graph-data-cache)
 444                        calc-graph-data-cache-limit)
 445                (setcdr calc-graph-data-cache
 446                        (cdr (cdr calc-graph-data-cache))))
 447            (setq calc-graph-data-cache (list (cons nil cache-env)))))
 448        (calc-graph-find-plot t t)
 449        (while (re-search-forward
 450                (if calc-graph-is-splot
 451                    "{\\([^{}:\n]+\\):\\([^{}:\n]+\\):\\([^{}:\n]+\\)}"
 452                  "{\\([^{}:\n]+\\)\\(:\\)\\([^{}:\n]+\\)}")
 453                nil t)
 454          (setq calc-graph-curve-num (1+ calc-graph-curve-num))
 455          (let* ((calc-graph-xname (buffer-substring (match-beginning 1) (match-end 1)))
 456                 (xvar (intern (concat "var-" calc-graph-xname)))
 457                 (calc-graph-xvalue (math-evaluate-expr (calc-var-value xvar)))
 458                 (calc-graph-y3name (and calc-graph-is-splot
 459                              (buffer-substring (match-beginning 2)
 460                                                (match-end 2))))
 461                 (y3var (and calc-graph-is-splot (intern (concat "var-" calc-graph-y3name))))
 462                 (calc-graph-y3value (and calc-graph-is-splot (calc-var-value y3var)))
 463                 (calc-graph-yname (buffer-substring (match-beginning 3) (match-end 3)))
 464                 (yvar (intern (concat "var-" calc-graph-yname)))
 465                 (calc-graph-yvalue (calc-var-value yvar))
 466                 filename)
 467            (delete-region (match-beginning 0) (match-end 0))
 468            (setq filename (calc-temp-file-name calc-graph-curve-num))
 469            (save-excursion
 470              (set-buffer calcbuf)
 471              (let (tempbuftop
 472                    (calc-graph-xp calc-graph-xvalue)
 473                    (calc-graph-yp calc-graph-yvalue)
 474                    (calc-graph-zp nil)
 475                    (calc-graph-xlow nil) (calc-graph-xhigh nil) (y3low nil) (y3high nil)
 476                    calc-graph-xvec calc-graph-xval calc-graph-xstep var-DUMMY
 477                    y3val calc-graph-y3step var-DUMMY2 (calc-graph-zval nil)
 478                    calc-graph-yvec calc-graph-yval calc-graph-ycache calc-graph-ycacheptr calc-graph-yvector
 479                    calc-graph-numsteps calc-graph-numsteps3
 480                    (calc-graph-keep-file (and (not calc-graph-is-splot) (file-exists-p filename)))
 481                    (calc-graph-stepcount 0)
 482                    (calc-symbolic-mode nil)
 483                    (calc-prefer-frac nil)
 484                    (calc-internal-prec (max 3 precision))
 485                    (calc-simplify-mode (and (not (memq calc-simplify-mode
 486                                                        '(none num)))
 487                                             calc-simplify-mode))
 488                    (calc-graph-blank t)
 489                    (calc-graph-non-blank nil)
 490                    (math-working-step 0)
 491                    (math-working-step-2 nil))
 492                (save-excursion
 493                  (if calc-graph-is-splot
 494                      (calc-graph-compute-3d)
 495                    (calc-graph-compute-2d))
 496                  (set-buffer tempbuf)
 497                  (goto-char (point-max))
 498                  (insert "\n" calc-graph-xname)
 499                  (if calc-graph-is-splot
 500                      (insert ":" calc-graph-y3name))
 501                  (insert ":" calc-graph-yname "\n\n")
 502                  (setq tempbuftop (point))
 503                  (let ((calc-group-digits nil)
 504                        (calc-leading-zeros nil)
 505                        (calc-number-radix 10)
 506                        (entry (and (not calc-graph-is-splot)
 507                                    (list calc-graph-xp calc-graph-yp calc-graph-xhigh calc-graph-numsteps))))
 508                    (or (equal entry
 509                               (nth 1 (nth (1+ calc-graph-curve-num)
 510                                           calc-graph-file-cache)))
 511                        (setq calc-graph-keep-file nil))
 512                    (setcar (cdr (nth (1+ calc-graph-curve-num) calc-graph-file-cache))
 513                            entry)
 514                    (or calc-graph-keep-file
 515                        (calc-graph-format-data)))
 516                  (or calc-graph-keep-file
 517                      (progn
 518                        (or calc-graph-non-blank
 519                            (error "No valid data points for %s:%s"
 520                                   calc-graph-xname calc-graph-yname))
 521                        (write-region tempbuftop (point-max) filename
 522                                      nil 'quiet))))))
 523            (insert (prin1-to-string filename))))
 524        (if calc-graph-surprise-splot
 525            (setcdr cache-env nil))
 526        (if (= calc-graph-curve-num 0)
 527            (progn
 528              (calc-gnuplot-command "clear")
 529              (calc-clear-command-flag 'clear-message)
 530              (message "No data to plot!"))
 531          (setq calc-graph-data-cache-limit (max calc-graph-curve-num
 532                                                 calc-graph-data-cache-limit)
 533                filename (calc-temp-file-name 0))
 534          (write-region (point-min) (point-max) filename nil 'quiet)
 535          (calc-gnuplot-command "load" (prin1-to-string filename))
 536          (or (equal output "STDOUT")
 537              calc-gnuplot-keep-outfile
 538              (progn   ; need to close the output file before printing/plotting
 539                (setq calc-graph-last-output "STDOUT")
 540                (calc-gnuplot-command "set output")))
 541          (let ((command (if printing
 542                             calc-gnuplot-print-command
 543                           (or calc-gnuplot-plot-command
 544                               (and (string-match "^dumb" device)
 545                                    'calc-graph-show-dumb)
 546                               (and tty-output
 547                                    'calc-graph-show-tty)))))
 548            (if command
 549                (if (stringp command)
 550                    (calc-gnuplot-command
 551                     "!" (format command
 552                                 (or tempoutfile
 553                                     calc-gnuplot-print-output)))
 554                  (if (symbolp command)
 555                      (funcall command output)
 556                    (eval command))))))))))
 557
 558 (defun calc-graph-compute-2d ()
 559   (if (setq calc-graph-yvec (eq (car-safe calc-graph-yvalue) 'vec))
 560       (if (= (setq calc-graph-numsteps (1- (length calc-graph-yvalue))) 0)
 561           (error "Can't plot an empty vector")
 562         (if (setq calc-graph-xvec (eq (car-safe calc-graph-xvalue) 'vec))
 563             (or (= (1- (length calc-graph-xvalue)) calc-graph-numsteps)
 564                 (error "%s and %s have different lengths" calc-graph-xname calc-graph-yname))
 565           (if (and (eq (car-safe calc-graph-xvalue) 'intv)
 566                    (math-constp calc-graph-xvalue))
 567               (setq calc-graph-xstep (math-div (math-sub (nth 3 calc-graph-xvalue)
 568                                               (nth 2 calc-graph-xvalue))
 569                                     (1- calc-graph-numsteps))
 570                     calc-graph-xvalue (nth 2 calc-graph-xvalue))
 571             (if (math-realp calc-graph-xvalue)
 572                 (setq calc-graph-xstep 1)
 573               (error "%s is not a suitable basis for %s" calc-graph-xname calc-graph-yname)))))
 574     (or (math-realp calc-graph-yvalue)
 575         (let ((arglist nil))
 576           (setq calc-graph-yvalue (math-evaluate-expr calc-graph-yvalue))
 577           (calc-default-formula-arglist calc-graph-yvalue)
 578           (or arglist
 579               (error "%s does not contain any unassigned variables" calc-graph-yname))
 580           (and (cdr arglist)
 581                (error "%s contains more than one variable: %s"
 582                       calc-graph-yname arglist))
 583           (setq calc-graph-yvalue (math-expr-subst calc-graph-yvalue
 584                                         (math-build-var-name (car arglist))
 585                                         '(var DUMMY var-DUMMY)))))
 586     (setq calc-graph-ycache (assoc calc-graph-yvalue calc-graph-data-cache))
 587     (delq calc-graph-ycache calc-graph-data-cache)
 588     (nconc calc-graph-data-cache
 589            (list (or calc-graph-ycache (setq calc-graph-ycache (list calc-graph-yvalue)))))
 590     (if (and (not (setq calc-graph-xvec (eq (car-safe calc-graph-xvalue) 'vec)))
 591              calc-graph-refine (cdr (cdr calc-graph-ycache)))
 592         (calc-graph-refine-2d)
 593       (calc-graph-recompute-2d))))
 594
 595 (defun calc-graph-refine-2d ()
 596   (setq calc-graph-keep-file nil
 597         calc-graph-ycacheptr (cdr calc-graph-ycache))
 598   (if (and (setq calc-graph-xval (calc-graph-find-command "xrange"))
 599            (string-match "\\`\\[\\([0-9.eE+-]*\\):\\([0-9.eE+-]*\\)\\]\\'"
 600                          calc-graph-xval))
 601       (let ((b2 (match-beginning 2))
 602             (e2 (match-end 2)))
 603         (setq calc-graph-xlow (math-read-number (substring calc-graph-xval
 604                                                 (match-beginning 1)
 605                                                 (match-end 1)))
 606               calc-graph-xhigh (math-read-number (substring calc-graph-xval b2 e2))))
 607     (if calc-graph-xlow
 608         (while (and (cdr calc-graph-ycacheptr)
 609                     (Math-lessp (car (nth 1 calc-graph-ycacheptr)) calc-graph-xlow))
 610           (setq calc-graph-ycacheptr (cdr calc-graph-ycacheptr)))))
 611   (setq math-working-step-2 (1- (length calc-graph-ycacheptr)))
 612   (while (and (cdr calc-graph-ycacheptr)
 613               (or (not calc-graph-xhigh)
 614                   (Math-lessp (car (car calc-graph-ycacheptr)) calc-graph-xhigh)))
 615     (setq var-DUMMY (math-div (math-add (car (car calc-graph-ycacheptr))
 616                                         (car (nth 1 calc-graph-ycacheptr)))
 617                               2)
 618           math-working-step (1+ math-working-step)
 619           calc-graph-yval (math-evaluate-expr calc-graph-yvalue))
 620     (setcdr calc-graph-ycacheptr (cons (cons var-DUMMY calc-graph-yval)
 621                             (cdr calc-graph-ycacheptr)))
 622     (setq calc-graph-ycacheptr (cdr (cdr calc-graph-ycacheptr))))
 623   (setq calc-graph-yp calc-graph-ycache
 624         calc-graph-numsteps 1000000))
 625
 626 (defun calc-graph-recompute-2d ()
 627   (setq calc-graph-ycacheptr calc-graph-ycache)
 628   (if calc-graph-xvec
 629       (setq calc-graph-numsteps (1- (length calc-graph-xvalue))
 630             calc-graph-yvector nil)
 631     (if (and (eq (car-safe calc-graph-xvalue) 'intv)
 632              (math-constp calc-graph-xvalue))
 633         (setq calc-graph-numsteps calc-graph-resolution
 634               calc-graph-yp nil
 635               calc-graph-xlow (nth 2 calc-graph-xvalue)
 636               calc-graph-xhigh (nth 3 calc-graph-xvalue)
 637               calc-graph-xstep (math-div (math-sub calc-graph-xhigh calc-graph-xlow)
 638                               (1- calc-graph-numsteps))
 639               calc-graph-xvalue (nth 2 calc-graph-xvalue))
 640       (error "%s is not a suitable basis for %s"
 641              calc-graph-xname calc-graph-yname)))
 642   (setq math-working-step-2 calc-graph-numsteps)
 643   (while (>= (setq calc-graph-numsteps (1- calc-graph-numsteps)) 0)
 644     (setq math-working-step (1+ math-working-step))
 645     (if calc-graph-xvec
 646         (progn
 647           (setq calc-graph-xp (cdr calc-graph-xp)
 648                 calc-graph-xval (car calc-graph-xp))
 649           (and (not (eq calc-graph-ycacheptr calc-graph-ycache))
 650                (consp (car calc-graph-ycacheptr))
 651                (not (Math-lessp (car (car calc-graph-ycacheptr)) calc-graph-xval))
 652                (setq calc-graph-ycacheptr calc-graph-ycache)))
 653       (if (= calc-graph-numsteps 0)
 654           (setq calc-graph-xval calc-graph-xhigh)   ; avoid cumulative roundoff
 655         (setq calc-graph-xval calc-graph-xvalue
 656               calc-graph-xvalue (math-add calc-graph-xvalue calc-graph-xstep))))
 657     (while (and (cdr calc-graph-ycacheptr)
 658                 (Math-lessp (car (nth 1 calc-graph-ycacheptr)) calc-graph-xval))
 659       (setq calc-graph-ycacheptr (cdr calc-graph-ycacheptr)))
 660     (or (and (cdr calc-graph-ycacheptr)
 661              (Math-equal (car (nth 1 calc-graph-ycacheptr)) calc-graph-xval))
 662         (progn
 663           (setq calc-graph-keep-file nil
 664                 var-DUMMY calc-graph-xval)
 665           (setcdr calc-graph-ycacheptr (cons (cons calc-graph-xval (math-evaluate-expr calc-graph-yvalue))
 666                                   (cdr calc-graph-ycacheptr)))))
 667     (setq calc-graph-ycacheptr (cdr calc-graph-ycacheptr))
 668     (if calc-graph-xvec
 669         (setq calc-graph-yvector (cons (cdr (car calc-graph-ycacheptr)) calc-graph-yvector))
 670       (or calc-graph-yp (setq calc-graph-yp calc-graph-ycacheptr))))
 671   (if calc-graph-xvec
 672       (setq calc-graph-xp calc-graph-xvalue
 673             calc-graph-yvec t
 674             calc-graph-yp (cons 'vec (nreverse calc-graph-yvector))
 675             calc-graph-numsteps (1- (length calc-graph-xp)))
 676     (setq calc-graph-numsteps 1000000)))
 677
 678 (defun calc-graph-compute-3d ()
 679   (if (setq calc-graph-yvec (eq (car-safe calc-graph-yvalue) 'vec))
 680       (if (math-matrixp calc-graph-yvalue)
 681           (progn
 682             (setq calc-graph-numsteps (1- (length calc-graph-yvalue))
 683                   calc-graph-numsteps3 (1- (length (nth 1 calc-graph-yvalue))))
 684             (if (eq (car-safe calc-graph-xvalue) 'vec)
 685                 (or (= (1- (length calc-graph-xvalue)) calc-graph-numsteps)
 686                     (error "%s has wrong length" calc-graph-xname))
 687               (if (and (eq (car-safe calc-graph-xvalue) 'intv)
 688                        (math-constp calc-graph-xvalue))
 689                   (setq calc-graph-xvalue (calcFunc-index calc-graph-numsteps
 690                                                (nth 2 calc-graph-xvalue)
 691                                                (math-div
 692                                                 (math-sub (nth 3 calc-graph-xvalue)
 693                                                           (nth 2 calc-graph-xvalue))
 694                                                 (1- calc-graph-numsteps))))
 695                 (if (math-realp calc-graph-xvalue)
 696                     (setq calc-graph-xvalue (calcFunc-index calc-graph-numsteps calc-graph-xvalue 1))
 697                   (error "%s is not a suitable basis for %s" calc-graph-xname calc-graph-yname))))
 698             (if (eq (car-safe calc-graph-y3value) 'vec)
 699                 (or (= (1- (length calc-graph-y3value)) calc-graph-numsteps3)
 700                     (error "%s has wrong length" calc-graph-y3name))
 701               (if (and (eq (car-safe calc-graph-y3value) 'intv)
 702                        (math-constp calc-graph-y3value))
 703                   (setq calc-graph-y3value (calcFunc-index calc-graph-numsteps3
 704                                                 (nth 2 calc-graph-y3value)
 705                                                 (math-div
 706                                                  (math-sub (nth 3 calc-graph-y3value)
 707                                                            (nth 2 calc-graph-y3value))
 708                                                  (1- calc-graph-numsteps3))))
 709                 (if (math-realp calc-graph-y3value)
 710                     (setq calc-graph-y3value (calcFunc-index calc-graph-numsteps3 calc-graph-y3value 1))
 711                   (error "%s is not a suitable basis for %s" calc-graph-y3name calc-graph-yname))))
 712             (setq calc-graph-xp nil
 713                   calc-graph-yp nil
 714                   calc-graph-zp nil
 715                   calc-graph-xvec t)
 716             (while (setq calc-graph-xvalue (cdr calc-graph-xvalue) calc-graph-yvalue (cdr calc-graph-yvalue))
 717               (setq calc-graph-xp (nconc calc-graph-xp (make-list (1+ calc-graph-numsteps3) (car calc-graph-xvalue)))
 718                     calc-graph-yp (nconc calc-graph-yp (cons 0 (copy-sequence (cdr calc-graph-y3value))))
 719                     calc-graph-zp (nconc calc-graph-zp (cons '(skip)
 720                                        (copy-sequence (cdr (car calc-graph-yvalue)))))))
 721             (setq calc-graph-numsteps (1- (* calc-graph-numsteps
 722                                              (1+ calc-graph-numsteps3)))))
 723         (if (= (setq calc-graph-numsteps (1- (length calc-graph-yvalue))) 0)
 724             (error "Can't plot an empty vector"))
 725         (or (and (eq (car-safe calc-graph-xvalue) 'vec)
 726                  (= (1- (length calc-graph-xvalue)) calc-graph-numsteps))
 727             (error "%s is not a suitable basis for %s" calc-graph-xname calc-graph-yname))
 728         (or (and (eq (car-safe calc-graph-y3value) 'vec)
 729                  (= (1- (length calc-graph-y3value)) calc-graph-numsteps))
 730             (error "%s is not a suitable basis for %s" calc-graph-y3name calc-graph-yname))
 731         (setq calc-graph-xp calc-graph-xvalue
 732               calc-graph-yp calc-graph-y3value
 733               calc-graph-zp calc-graph-yvalue
 734               calc-graph-xvec t))
 735     (or (math-realp calc-graph-yvalue)
 736         (let ((arglist nil))
 737           (setq calc-graph-yvalue (math-evaluate-expr calc-graph-yvalue))
 738           (calc-default-formula-arglist calc-graph-yvalue)
 739           (setq arglist (sort arglist 'string-lessp))
 740           (or (cdr arglist)
 741               (error "%s does not contain enough unassigned variables" calc-graph-yname))
 742           (and (cdr (cdr arglist))
 743                (error "%s contains too many variables: %s" calc-graph-yname arglist))
 744           (setq calc-graph-yvalue (math-multi-subst calc-graph-yvalue
 745                                          (mapcar 'math-build-var-name
 746                                                  arglist)
 747                                          '((var DUMMY var-DUMMY)
 748                                            (var DUMMY2 var-DUMMY2))))))
 749     (if (setq calc-graph-xvec (eq (car-safe calc-graph-xvalue) 'vec))
 750         (setq calc-graph-numsteps (1- (length calc-graph-xvalue)))
 751       (if (and (eq (car-safe calc-graph-xvalue) 'intv)
 752                (math-constp calc-graph-xvalue))
 753           (setq calc-graph-numsteps calc-graph-resolution
 754                 calc-graph-xvalue (calcFunc-index calc-graph-numsteps
 755                                        (nth 2 calc-graph-xvalue)
 756                                        (math-div (math-sub (nth 3 calc-graph-xvalue)
 757                                                            (nth 2 calc-graph-xvalue))
 758                                                  (1- calc-graph-numsteps))))
 759         (error "%s is not a suitable basis for %s"
 760                calc-graph-xname calc-graph-yname)))
 761     (if (eq (car-safe calc-graph-y3value) 'vec)
 762         (setq calc-graph-numsteps3 (1- (length calc-graph-y3value)))
 763       (if (and (eq (car-safe calc-graph-y3value) 'intv)
 764                (math-constp calc-graph-y3value))
 765           (setq calc-graph-numsteps3 calc-graph-resolution
 766                 calc-graph-y3value (calcFunc-index calc-graph-numsteps3
 767                                         (nth 2 calc-graph-y3value)
 768                                         (math-div (math-sub (nth 3 calc-graph-y3value)
 769                                                             (nth 2 calc-graph-y3value))
 770                                                   (1- calc-graph-numsteps3))))
 771         (error "%s is not a suitable basis for %s"
 772                calc-graph-y3name calc-graph-yname)))
 773     (setq calc-graph-xp nil
 774           calc-graph-yp nil
 775           calc-graph-zp nil
 776           calc-graph-xvec t)
 777     (setq math-working-step 0)
 778     (while (setq calc-graph-xvalue (cdr calc-graph-xvalue))
 779       (setq calc-graph-xp (nconc calc-graph-xp (make-list (1+ calc-graph-numsteps3) (car calc-graph-xvalue)))
 780             calc-graph-yp (nconc calc-graph-yp (cons 0 (copy-sequence (cdr calc-graph-y3value))))
 781             calc-graph-zp (cons '(skip) calc-graph-zp)
 782             calc-graph-y3step calc-graph-y3value
 783             var-DUMMY (car calc-graph-xvalue)
 784             math-working-step-2 0
 785             math-working-step (1+ math-working-step))
 786       (while (setq calc-graph-y3step (cdr calc-graph-y3step))
 787         (setq math-working-step-2 (1+ math-working-step-2)
 788               var-DUMMY2 (car calc-graph-y3step)
 789               calc-graph-zp (cons (math-evaluate-expr calc-graph-yvalue) calc-graph-zp))))
 790     (setq calc-graph-zp (nreverse calc-graph-zp)
 791           calc-graph-numsteps (1- (* calc-graph-numsteps (1+ calc-graph-numsteps3))))))
 792
 793 (defun calc-graph-format-data ()
 794   (while (<= (setq calc-graph-stepcount (1+ calc-graph-stepcount)) calc-graph-numsteps)
 795     (if calc-graph-xvec
 796         (setq calc-graph-xp (cdr calc-graph-xp)
 797               calc-graph-xval (car calc-graph-xp)
 798               calc-graph-yp (cdr calc-graph-yp)
 799               calc-graph-yval (car calc-graph-yp)
 800               calc-graph-zp (cdr calc-graph-zp)
 801               calc-graph-zval (car calc-graph-zp))
 802       (if calc-graph-yvec
 803           (setq calc-graph-xval calc-graph-xvalue
 804                 calc-graph-xvalue (math-add calc-graph-xvalue calc-graph-xstep)
 805                 calc-graph-yp (cdr calc-graph-yp)
 806                 calc-graph-yval (car calc-graph-yp))
 807         (setq calc-graph-xval (car (car calc-graph-yp))
 808               calc-graph-yval (cdr (car calc-graph-yp))
 809               calc-graph-yp (cdr calc-graph-yp))
 810         (if (or (not calc-graph-yp)
 811                 (and calc-graph-xhigh (equal calc-graph-xval calc-graph-xhigh)))
 812             (setq calc-graph-numsteps 0))))
 813     (if calc-graph-is-splot
 814         (if (and (eq (car-safe calc-graph-zval) 'calcFunc-xyz)
 815                  (= (length calc-graph-zval) 4))
 816             (setq calc-graph-xval (nth 1 calc-graph-zval)
 817                   calc-graph-yval (nth 2 calc-graph-zval)
 818                   calc-graph-zval (nth 3 calc-graph-zval)))
 819       (if (and (eq (car-safe calc-graph-yval) 'calcFunc-xyz)
 820                (= (length calc-graph-yval) 4))
 821           (progn
 822             (or calc-graph-surprise-splot
 823                 (save-excursion
 824                   (set-buffer (get-buffer-create "*Gnuplot Temp*"))
 825                   (save-excursion
 826                     (goto-char (point-max))
 827                     (re-search-backward "^plot[ \t]")
 828                     (insert "set parametric\ns")
 829                     (setq calc-graph-surprise-splot t))))
 830             (setq calc-graph-xval (nth 1 calc-graph-yval)
 831                   calc-graph-zval (nth 3 calc-graph-yval)
 832                   calc-graph-yval (nth 2 calc-graph-yval)))
 833         (if (and (eq (car-safe calc-graph-yval) 'calcFunc-xy)
 834                  (= (length calc-graph-yval) 3))
 835             (setq calc-graph-xval (nth 1 calc-graph-yval)
 836                   calc-graph-yval (nth 2 calc-graph-yval)))))
 837     (if (and (Math-realp calc-graph-xval)
 838              (Math-realp calc-graph-yval)
 839              (or (not calc-graph-zval) (Math-realp calc-graph-zval)))
 840         (progn
 841           (setq calc-graph-blank nil
 842                 calc-graph-non-blank t)
 843           (if (Math-integerp calc-graph-xval)
 844               (insert (math-format-number calc-graph-xval))
 845             (if (eq (car calc-graph-xval) 'frac)
 846                 (setq calc-graph-xval (math-float calc-graph-xval)))
 847             (insert (math-format-number (nth 1 calc-graph-xval))
 848                     "e" (int-to-string (nth 2 calc-graph-xval))))
 849           (insert " ")
 850           (if (Math-integerp calc-graph-yval)
 851               (insert (math-format-number calc-graph-yval))
 852             (if (eq (car calc-graph-yval) 'frac)
 853                 (setq calc-graph-yval (math-float calc-graph-yval)))
 854             (insert (math-format-number (nth 1 calc-graph-yval))
 855                     "e" (int-to-string (nth 2 calc-graph-yval))))
 856           (if calc-graph-zval
 857               (progn
 858                 (insert " ")
 859                 (if (Math-integerp calc-graph-zval)
 860                     (insert (math-format-number calc-graph-zval))
 861                   (if (eq (car calc-graph-zval) 'frac)
 862                       (setq calc-graph-zval (math-float calc-graph-zval)))
 863                   (insert (math-format-number (nth 1 calc-graph-zval))
 864                           "e" (int-to-string (nth 2 calc-graph-zval))))))
 865           (insert "\n"))
 866       (and (not (equal calc-graph-zval '(skip)))
 867            (boundp 'var-PlotRejects)
 868            (eq (car-safe var-PlotRejects) 'vec)
 869            (nconc var-PlotRejects
 870                   (list (list 'vec
 871                               calc-graph-curve-num
 872                               calc-graph-stepcount
 873                               calc-graph-xval calc-graph-yval)))
 874            (calc-refresh-evaltos 'var-PlotRejects))
 875       (or calc-graph-blank
 876           (progn
 877             (insert "\n")
 878             (setq calc-graph-blank t))))))
 879
 880 (defun calc-temp-file-name (num)
 881   (while (<= (length calc-graph-file-cache) (1+ num))
 882     (setq calc-graph-file-cache (nconc calc-graph-file-cache (list nil))))
 883   (car (or (nth (1+ num) calc-graph-file-cache)
 884            (setcar (nthcdr (1+ num) calc-graph-file-cache)
 885                    (list (make-temp-file
 886                           (concat calc-gnuplot-tempfile
 887                                   (if (<= num 0)
 888                                       (char-to-string (- ?A num))
 889                                     (int-to-string num))))
 890                          nil)))))
 891
 892 (defun calc-graph-delete-temps ()
 893   (while calc-graph-file-cache
 894     (and (car calc-graph-file-cache)
 895          (file-exists-p (car (car calc-graph-file-cache)))
 896          (condition-case err
 897              (delete-file (car (car calc-graph-file-cache)))
 898            (error nil)))
 899     (setq calc-graph-file-cache (cdr calc-graph-file-cache))))
 900
 901 (defun calc-graph-kill-hook ()
 902   (calc-graph-delete-temps))
 903
 904 (defun calc-graph-show-tty (output)
 905   "Default calc-gnuplot-plot-command for \"tty\" output mode.
 906 This is useful for tek40xx and other graphics-terminal types."
 907   (call-process-region 1 1 shell-file-name
 908                        nil calc-gnuplot-buffer nil
 909                        "-c" (format "cat %s >/dev/tty; rm %s" output output)))
 910
 911 (defvar calc-dumb-map nil
 912   "The keymap for the \"dumb\" terminal plot.")
 913
 914 (defun calc-graph-show-dumb (&optional output)
 915   "Default calc-gnuplot-plot-command for Pinard's \"dumb\" terminal type.
 916 This \"dumb\" driver will be present in Gnuplot 3.0."
 917   (interactive)
 918   (save-window-excursion
 919     (switch-to-buffer calc-gnuplot-buffer)
 920     (delete-other-windows)
 921     (goto-char calc-gnuplot-trail-mark)
 922     (or (search-forward "\f" nil t)
 923         (sleep-for 1))
 924     (goto-char (point-max))
 925     (re-search-backward "\f\\|^[ \t]+\\^$\\|G N U P L O T")
 926     (if (looking-at "\f")
 927         (progn
 928           (forward-char 1)
 929           (if (eolp) (forward-line 1))
 930           (or (calc-graph-find-command "time")
 931               (calc-graph-find-command "title")
 932               (calc-graph-find-command "ylabel")
 933               (let ((pt (point)))
 934                 (insert-before-markers (format "(%s)" (current-time-string)))
 935                 (goto-char pt)))
 936           (set-window-start (selected-window) (point))
 937           (goto-char (point-max)))
 938       (end-of-line)
 939       (backward-char 1)
 940       (recenter '(4)))
 941     (or calc-dumb-map
 942         (progn
 943           (setq calc-dumb-map (make-sparse-keymap))
 944           (define-key calc-dumb-map "\n" 'scroll-up)
 945           (define-key calc-dumb-map " " 'scroll-up)
 946           (define-key calc-dumb-map "\177" 'scroll-down)
 947           (define-key calc-dumb-map "<" 'scroll-left)
 948           (define-key calc-dumb-map ">" 'scroll-right)
 949           (define-key calc-dumb-map "{" 'scroll-down)
 950           (define-key calc-dumb-map "}" 'scroll-up)
 951           (define-key calc-dumb-map "q" 'exit-recursive-edit)
 952           (define-key calc-dumb-map "\C-c\C-c" 'exit-recursive-edit)))
 953     (use-local-map calc-dumb-map)
 954     (setq truncate-lines t)
 955     (message "Type `q' or `C-c C-c' to return to Calc")
 956     (recursive-edit)
 957     (bury-buffer "*Gnuplot Trail*")))
 958
 959 (defun calc-graph-clear ()
 960   (interactive)
 961   (if calc-graph-last-device
 962       (if (or (equal calc-graph-last-device "x11")
 963               (equal calc-graph-last-device "X11"))
 964           (calc-gnuplot-command "set output"
 965                                 (if (equal calc-graph-last-output "STDOUT")
 966                                     ""
 967                                   (prin1-to-string calc-graph-last-output)))
 968         (calc-gnuplot-command "clear"))))
 969
 970 (defun calc-graph-title-x (title)
 971   (interactive "sX axis title: ")
 972   (calc-graph-set-command "xlabel" (if (not (equal title ""))
 973                                        (prin1-to-string title))))
 974
 975 (defun calc-graph-title-y (title)
 976   (interactive "sY axis title: ")
 977   (calc-graph-set-command "ylabel" (if (not (equal title ""))
 978                                        (prin1-to-string title))))
 979
 980 (defun calc-graph-title-z (title)
 981   (interactive "sZ axis title: ")
 982   (calc-graph-set-command "zlabel" (if (not (equal title ""))
 983                                        (prin1-to-string title))))
 984
 985 (defun calc-graph-range-x (range)
 986   (interactive "sX axis range: ")
 987   (calc-graph-set-range "xrange" range))
 988
 989 (defun calc-graph-range-y (range)
 990   (interactive "sY axis range: ")
 991   (calc-graph-set-range "yrange" range))
 992
 993 (defun calc-graph-range-z (range)
 994   (interactive "sZ axis range: ")
 995   (calc-graph-set-range "zrange" range))
 996
 997 (defun calc-graph-set-range (cmd range)
 998   (if (equal range "$")
 999       (calc-wrapper
1000        (let ((val (calc-top-n 1)))
1001          (if (and (eq (car-safe val) 'intv) (math-constp val))
1002              (setq range (concat
1003                           (math-format-number (math-float (nth 2 val))) ":"
1004                           (math-format-number (math-float (nth 3 val)))))
1005            (if (and (eq (car-safe val) 'vec)
1006                     (= (length val) 3))
1007                (setq range (concat
1008                             (math-format-number (math-float (nth 1 val))) ":"
1009                             (math-format-number (math-float (nth 2 val)))))
1010              (error "Range specification must be an interval or 2-vector")))
1011          (calc-pop-stack 1))))
1012   (if (string-match "\\[.+\\]" range)
1013       (setq range (substring range 1 -1)))
1014   (if (and (not (string-match ":" range))
1015            (or (string-match "," range)
1016                (string-match " " range)))
1017       (aset range (match-beginning 0) ?\:))
1018   (calc-graph-set-command cmd (if (not (equal range ""))
1019                                   (concat "[" range "]"))))
1020
1021 (defun calc-graph-log-x (flag)
1022   (interactive "P")
1023   (calc-graph-set-log flag 0 0))
1024
1025 (defun calc-graph-log-y (flag)
1026   (interactive "P")
1027   (calc-graph-set-log 0 flag 0))
1028
1029 (defun calc-graph-log-z (flag)
1030   (interactive "P")
1031   (calc-graph-set-log 0 0 flag))
1032
1033 (defun calc-graph-set-log (xflag yflag zflag)
1034   (let* ((old (or (calc-graph-find-command "logscale") ""))
1035          (xold (string-match "x" old))
1036          (yold (string-match "y" old))
1037          (zold (string-match "z" old))
1038          str)
1039     (setq str (concat (if (if xflag
1040                               (if (eq xflag 0) xold
1041                                 (> (prefix-numeric-value xflag) 0))
1042                             (not xold)) "x" "")
1043                       (if (if yflag
1044                               (if (eq yflag 0) yold
1045                                 (> (prefix-numeric-value yflag) 0))
1046                             (not yold)) "y" "")
1047                       (if (if zflag
1048                               (if (eq zflag 0) zold
1049                                 (> (prefix-numeric-value zflag) 0))
1050                             (not zold)) "z" "")))
1051     (calc-graph-set-command "logscale" (if (not (equal str "")) str))))
1052
1053 (defun calc-graph-line-style (style)
1054   (interactive "P")
1055   (calc-graph-set-styles (and style (prefix-numeric-value style)) t))
1056
1057 (defun calc-graph-point-style (style)
1058   (interactive "P")
1059   (calc-graph-set-styles t (and style (prefix-numeric-value style))))
1060
1061 (defun calc-graph-set-styles (lines points)
1062   (calc-graph-init)
1063   (save-excursion
1064     (set-buffer calc-gnuplot-input)
1065     (or (calc-graph-find-plot nil nil)
1066         (error "No data points have been set!"))
1067     (let ((base (point))
1068           (mode nil) (lstyle nil) (pstyle nil)
1069           start end lenbl penbl)
1070       (re-search-forward "[,\n]")
1071       (forward-char -1)
1072       (setq end (point) start end)
1073       (goto-char base)
1074       (if (looking-at "[^,\n]*[^,\n \t]\\([ \t]+with\\)")
1075           (progn
1076             (setq start (match-beginning 1))
1077             (goto-char (match-end 0))
1078             (if (looking-at "[ \t]+\\([a-z]+\\)")
1079                 (setq mode (buffer-substring (match-beginning 1)
1080                                              (match-end 1))))
1081             (if (looking-at "[ \ta-z]+\\([0-9]+\\)")
1082                 (setq lstyle (string-to-number
1083                               (buffer-substring (match-beginning 1)
1084                                                 (match-end 1)))))
1085             (if (looking-at "[ \ta-z]+[0-9]+[ \t]+\\([0-9]+\\)")
1086                 (setq pstyle (string-to-number
1087                               (buffer-substring (match-beginning 1)
1088                                                 (match-end 1)))))))
1089       (setq lenbl (or (equal mode "lines") (equal mode "linespoints"))
1090             penbl (or (equal mode "points") (equal mode "linespoints")))
1091       (if lines
1092           (or (eq lines t)
1093               (setq lstyle lines
1094                     lenbl (>= lines 0)))
1095         (setq lenbl (not lenbl)))
1096       (if points
1097           (or (eq points t)
1098               (setq pstyle points
1099                     penbl (>= points 0)))
1100         (setq penbl (not penbl)))
1101       (delete-region start end)
1102       (goto-char start)
1103       (insert " with "
1104               (if lenbl
1105                   (if penbl "linespoints" "lines")
1106                 (if penbl "points" "dots")))
1107       (if (and pstyle (> pstyle 0))
1108           (insert " " (if (and lstyle (> lstyle 0)) (int-to-string lstyle) "1")
1109                   " " (int-to-string pstyle))
1110         (if (and lstyle (> lstyle 0))
1111             (insert " " (int-to-string lstyle))))))
1112   (calc-graph-view-commands))
1113
1114 (defun calc-graph-zero-x (flag)
1115   (interactive "P")
1116   (calc-graph-set-command "noxzeroaxis"
1117                           (and (if flag
1118                                    (<= (prefix-numeric-value flag) 0)
1119                                  (not (calc-graph-find-command "noxzeroaxis")))
1120                                " ")))
1121
1122 (defun calc-graph-zero-y (flag)
1123   (interactive "P")
1124   (calc-graph-set-command "noyzeroaxis"
1125                           (and (if flag
1126                                    (<= (prefix-numeric-value flag) 0)
1127                                  (not (calc-graph-find-command "noyzeroaxis")))
1128                                " ")))
1129
1130 (defun calc-graph-name (name)
1131   (interactive "sTitle for current curve: ")
1132   (calc-graph-init)
1133   (save-excursion
1134     (set-buffer calc-gnuplot-input)
1135     (or (calc-graph-find-plot nil nil)
1136         (error "No data points have been set!"))
1137     (let ((base (point))
1138           start
1139           end)
1140       (re-search-forward "[,\n]\\|[ \t]+with")
1141       (setq end (match-beginning 0))
1142       (goto-char base)
1143       (if (looking-at "[^,\n]*[^,\n \t]\\([ \t]+title\\)")
1144           (progn
1145             (goto-char (match-beginning 1))
1146             (delete-region (point) end))
1147         (goto-char end))
1148       (insert " title " (prin1-to-string name))))
1149   (calc-graph-view-commands))
1150
1151 (defun calc-graph-hide (flag)
1152   (interactive "P")
1153   (calc-graph-init)
1154   (and (calc-graph-find-plot nil nil)
1155        (progn
1156          (or (looking-at "{")
1157              (error "Can't hide this curve (wrong format)"))
1158          (forward-char 1)
1159          (if (looking-at "*")
1160              (if (or (null flag) (<= (prefix-numeric-value flag) 0))
1161                  (delete-char 1))
1162            (if (or (null flag) (> (prefix-numeric-value flag) 0))
1163                (insert "*"))))))
1164
1165 (defun calc-graph-header (title)
1166   (interactive "sTitle for entire graph: ")
1167   (calc-graph-set-command "title" (if (not (equal title ""))
1168                                       (prin1-to-string title))))
1169
1170 (defun calc-graph-border (flag)
1171   (interactive "P")
1172   (calc-graph-set-command "noborder"
1173                           (and (if flag
1174                                    (<= (prefix-numeric-value flag) 0)
1175                                  (not (calc-graph-find-command "noborder")))
1176                                " ")))
1177
1178 (defun calc-graph-grid (flag)
1179   (interactive "P")
1180   (calc-graph-set-command "grid" (and (if flag
1181                                           (> (prefix-numeric-value flag) 0)
1182                                         (not (calc-graph-find-command "grid")))
1183                                       " ")))
1184
1185 (defun calc-graph-key (flag)
1186   (interactive "P")
1187   (calc-graph-set-command "key" (and (if flag
1188                                          (> (prefix-numeric-value flag) 0)
1189                                        (not (calc-graph-find-command "key")))
1190                                      " ")))
1191
1192 (defun calc-graph-num-points (res flag)
1193   (interactive "sNumber of data points: \nP")
1194   (if flag
1195       (if (> (prefix-numeric-value flag) 0)
1196           (if (equal res "")
1197               (message "Default resolution is %d"
1198                        calc-graph-default-resolution)
1199             (setq calc-graph-default-resolution (string-to-number res)))
1200         (if (equal res "")
1201             (message "Default 3D resolution is %d"
1202                      calc-graph-default-resolution-3d)
1203           (setq calc-graph-default-resolution-3d (string-to-number res))))
1204     (calc-graph-set-command "samples" (if (not (equal res "")) res))))
1205
1206 (defun calc-graph-device (name flag)
1207   (interactive "sDevice name: \nP")
1208   (if (equal name "?")
1209       (progn
1210         (calc-gnuplot-command "set terminal")
1211         (calc-graph-view-trail))
1212     (if flag
1213         (if (> (prefix-numeric-value flag) 0)
1214             (if (equal name "")
1215                 (message "Default GNUPLOT device is \"%s\""
1216                          calc-gnuplot-default-device)
1217               (setq calc-gnuplot-default-device name))
1218           (if (equal name "")
1219               (message "GNUPLOT device for Print command is \"%s\""
1220                        calc-gnuplot-print-device)
1221             (setq calc-gnuplot-print-device name)))
1222       (calc-graph-set-command "terminal" (if (not (equal name ""))
1223                                              name)))))
1224
1225 (defun calc-graph-output (name flag)
1226   (interactive "FOutput file name: \np")
1227   (cond ((string-match "\\<[aA][uU][tT][oO]$" name)
1228          (setq name "auto"))
1229         ((string-match "\\<[tT][tT][yY]$" name)
1230          (setq name "tty"))
1231         ((string-match "\\<[sS][tT][dD][oO][uU][tT]$" name)
1232          (setq name "STDOUT"))
1233         ((equal (file-name-nondirectory name) "")
1234          (setq name ""))
1235         (t (setq name (expand-file-name name))))
1236   (if flag
1237       (if (> (prefix-numeric-value flag) 0)
1238           (if (equal name "")
1239               (message "Default GNUPLOT output file is \"%s\""
1240                        calc-gnuplot-default-output)
1241             (setq calc-gnuplot-default-output name))
1242         (if (equal name "")
1243             (message "GNUPLOT output file for Print command is \"%s\""
1244                      calc-gnuplot-print-output)
1245           (setq calc-gnuplot-print-output name)))
1246     (calc-graph-set-command "output" (if (not (equal name ""))
1247                                          (prin1-to-string name)))))
1248
1249 (defun calc-graph-display (name)
1250   (interactive "sX display name: ")
1251   (if (equal name "")
1252       (message "Current X display is \"%s\""
1253                (or calc-gnuplot-display "<none>"))
1254     (setq calc-gnuplot-display name)
1255     (if (calc-gnuplot-alive)
1256         (calc-gnuplot-command "exit"))))
1257
1258 (defun calc-graph-geometry (name)
1259   (interactive "sX geometry spec (or \"default\"): ")
1260   (if (equal name "")
1261       (message "Current X geometry is \"%s\""
1262                (or calc-gnuplot-geometry "default"))
1263     (setq calc-gnuplot-geometry (and (not (equal name "default")) name))
1264     (if (calc-gnuplot-alive)
1265         (calc-gnuplot-command "exit"))))
1266
1267 (defun calc-graph-find-command (cmd)
1268   (calc-graph-init)
1269   (save-excursion
1270     (set-buffer calc-gnuplot-input)
1271     (goto-char (point-min))
1272     (if (re-search-forward (concat "^set[ \t]+" cmd "[ \t]*\\(.*\\)$") nil t)
1273         (buffer-substring (match-beginning 1) (match-end 1)))))
1274
1275 (defun calc-graph-set-command (cmd &rest args)
1276   (calc-graph-init)
1277   (save-excursion
1278     (set-buffer calc-gnuplot-input)
1279     (goto-char (point-min))
1280     (if (re-search-forward (concat "^set[ \t]+" cmd "[ \t\n]") nil t)
1281         (progn
1282           (forward-char -1)
1283           (end-of-line)
1284           (let ((end (point)))
1285             (beginning-of-line)
1286             (delete-region (point) (1+ end))))
1287       (if (calc-graph-find-plot t t)
1288           (if (eq (preceding-char) ?\n)
1289               (forward-char -1))
1290         (goto-char (1- (point-max)))))
1291     (if (and args (car args))
1292         (progn
1293           (or (bolp)
1294               (insert "\n"))
1295           (insert "set " (mapconcat 'identity (cons cmd args) " ") "\n"))))
1296   (calc-graph-view-commands))
1297
1298 (defun calc-graph-command (cmd)
1299   (interactive "sGNUPLOT command: ")
1300   (calc-wrapper
1301    (calc-graph-init)
1302    (calc-graph-view-trail)
1303    (calc-gnuplot-command cmd)
1304    (accept-process-output)
1305    (calc-graph-view-trail)))
1306
1307 (defun calc-graph-kill (&optional no-view)
1308   (interactive)
1309   (calc-graph-delete-temps)
1310   (if (calc-gnuplot-alive)
1311       (calc-wrapper
1312        (or no-view (calc-graph-view-trail))
1313        (let ((calc-graph-no-wait t))
1314          (calc-gnuplot-command "exit"))
1315        (sit-for 1)
1316        (if (process-status calc-gnuplot-process)
1317            (delete-process calc-gnuplot-process))
1318        (setq calc-gnuplot-process nil))))
1319
1320 (defun calc-graph-quit ()
1321   (interactive)
1322   (if (get-buffer-window calc-gnuplot-input)
1323       (calc-graph-view-commands t))
1324   (if (get-buffer-window calc-gnuplot-buffer)
1325       (calc-graph-view-trail t))
1326   (calc-graph-kill t))
1327
1328 (defun calc-graph-view-commands (&optional no-need)
1329   (interactive "p")
1330   (or calc-graph-no-auto-view (calc-graph-init-buffers))
1331   (calc-graph-view calc-gnuplot-input calc-gnuplot-buffer (null no-need)))
1332
1333 (defun calc-graph-view-trail (&optional no-need)
1334   (interactive "p")
1335   (or calc-graph-no-auto-view (calc-graph-init-buffers))
1336   (calc-graph-view calc-gnuplot-buffer calc-gnuplot-input (null no-need)))
1337
1338 (defun calc-graph-view (buf other-buf need)
1339   (let (win)
1340     (or calc-graph-no-auto-view
1341         (if (setq win (get-buffer-window buf))
1342             (or need
1343                 (and (eq buf calc-gnuplot-buffer)
1344                      (save-excursion
1345                        (set-buffer buf)
1346                        (not (pos-visible-in-window-p (point-max) win))))
1347                 (progn
1348                   (bury-buffer buf)
1349                   (bury-buffer other-buf)
1350                   (let ((curwin (selected-window)))
1351                     (select-window win)
1352                     (switch-to-buffer nil)
1353                     (select-window curwin))))
1354           (if (setq win (get-buffer-window other-buf))
1355               (set-window-buffer win buf)
1356             (if (eq major-mode 'calc-mode)
1357                 (if (or need
1358                         (< (window-height) (1- (frame-height))))
1359                     (display-buffer buf))
1360               (switch-to-buffer buf)))))
1361     (save-excursion
1362       (set-buffer buf)
1363       (if (and (eq buf calc-gnuplot-buffer)
1364                (setq win (get-buffer-window buf))
1365                (not (pos-visible-in-window-p (point-max) win)))
1366           (progn
1367             (goto-char (point-max))
1368             (vertical-motion (- 6 (window-height win)))
1369             (set-window-start win (point))
1370             (goto-char (point-max)))))
1371     (or calc-graph-no-auto-view (sit-for 0))))
1372
1373 (defun calc-gnuplot-check-for-errors ()
1374   (if (save-excursion
1375         (prog2
1376          (progn
1377            (set-buffer calc-gnuplot-buffer)
1378            (goto-char calc-gnuplot-last-error-pos))
1379          (re-search-forward "^[ \t]+\\^$" nil t)
1380          (goto-char (point-max))
1381          (setq calc-gnuplot-last-error-pos (point-max))))
1382       (calc-graph-view-trail)))
1383
1384 (defun calc-gnuplot-command (&rest args)
1385   (calc-graph-init)
1386   (let ((cmd (concat (mapconcat 'identity args " ") "\n")))
1387     (accept-process-output)
1388     (save-excursion
1389       (set-buffer calc-gnuplot-buffer)
1390       (calc-gnuplot-check-for-errors)
1391       (goto-char (point-max))
1392       (setq calc-gnuplot-trail-mark (point))
1393       (or (>= calc-gnuplot-version 3)
1394           (insert cmd))
1395       (set-marker (process-mark calc-gnuplot-process) (point))
1396       (process-send-string calc-gnuplot-process cmd)
1397       (if (get-buffer-window calc-gnuplot-buffer)
1398           (calc-graph-view-trail))
1399       (accept-process-output (and (not calc-graph-no-wait)
1400                                   calc-gnuplot-process))
1401       (calc-gnuplot-check-for-errors)
1402       (if (get-buffer-window calc-gnuplot-buffer)
1403           (calc-graph-view-trail)))))
1404
1405 (defun calc-graph-init-buffers ()
1406   (or (and calc-gnuplot-buffer
1407            (buffer-name calc-gnuplot-buffer))
1408       (setq calc-gnuplot-buffer (get-buffer-create "*Gnuplot Trail*")))
1409   (or (and calc-gnuplot-input
1410            (buffer-name calc-gnuplot-input))
1411       (setq calc-gnuplot-input (get-buffer-create "*Gnuplot Commands*"))))
1412
1413 (defun calc-graph-init ()
1414   (or (calc-gnuplot-alive)
1415       (let ((process-connection-type t)
1416             origin)
1417         (if calc-gnuplot-process
1418             (progn
1419               (delete-process calc-gnuplot-process)
1420               (setq calc-gnuplot-process nil)))
1421         (calc-graph-init-buffers)
1422         (save-excursion
1423           (set-buffer calc-gnuplot-buffer)
1424           (insert "\nStarting gnuplot...\n")
1425           (setq origin (point)))
1426         (setq calc-graph-last-device nil)
1427         (setq calc-graph-last-output nil)
1428         (condition-case err
1429             (let ((args (append (and calc-gnuplot-display
1430                                      (not (equal calc-gnuplot-display
1431                                                  (getenv "DISPLAY")))
1432                                      (list "-display"
1433                                            calc-gnuplot-display))
1434                                 (and calc-gnuplot-geometry
1435                                      (list "-geometry"
1436                                            calc-gnuplot-geometry)))))
1437               (setq calc-gnuplot-process
1438                     (apply 'start-process
1439                            "gnuplot"
1440                            calc-gnuplot-buffer
1441                            calc-gnuplot-name
1442                            args))
1443               (set-process-query-on-exit-flag calc-gnuplot-process nil))
1444           (file-error
1445            (error "Sorry, can't find \"%s\" on your system"
1446                   calc-gnuplot-name)))
1447         (save-excursion
1448           (set-buffer calc-gnuplot-buffer)
1449           (while (and (not (save-excursion
1450                              (goto-char origin)
1451                              (search-forward "gnuplot> " nil t)))
1452                       (memq (process-status calc-gnuplot-process) '(run stop)))
1453             (accept-process-output calc-gnuplot-process))
1454           (or (memq (process-status calc-gnuplot-process) '(run stop))
1455               (error "Unable to start GNUPLOT process"))
1456           (if (save-excursion
1457                 (goto-char origin)
1458                 (re-search-forward
1459                  "G N U P L O T.*\n.*version \\([0-9]+\\)\\." nil t))
1460               (setq calc-gnuplot-version (string-to-number (buffer-substring
1461                                                          (match-beginning 1)
1462                                                          (match-end 1))))
1463             (setq calc-gnuplot-version 1))
1464           (goto-char (point-max)))))
1465   (save-excursion
1466     (set-buffer calc-gnuplot-input)
1467     (if (= (buffer-size) 0)
1468         (insert "# Commands for running gnuplot\n\n\n")
1469       (or calc-graph-no-auto-view
1470           (eq (char-after (1- (point-max))) ?\n)
1471           (progn
1472             (goto-char (point-max))
1473             (insert "\n"))))))
1474
1475 (provide 'calc-graph)
1476
1477 ;;; arch-tag: e4b06a52-c386-4d54-a2bb-7c0a0ef533c2
1478 ;;; calc-graph.el ends here