1 ;;; calc.el --- the GNU Emacs calculator
3 ;; Copyright (C) 1990, 1991, 1992, 1993, 2001, 2002, 2003, 2004, 2005,
4 ;; 2006, 2007, 2008, 2009 Free Software Foundation, Inc.
6 ;; Author: David Gillespie <daveg@synaptics.com>
7 ;; Maintainer: Jay Belanger <jay.p.belanger@gmail.com>
8 ;; Keywords: convenience, extensions
10 ;; This file is part of GNU Emacs.
12 ;; GNU Emacs is free software: you can redistribute it and/or modify
13 ;; it under the terms of the GNU General Public License as published by
14 ;; the Free Software Foundation, either version 3 of the License, or
15 ;; (at your option) any later version.
17 ;; GNU Emacs is distributed in the hope that it will be useful,
18 ;; but WITHOUT ANY WARRANTY; without even the implied warranty of
19 ;; MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
20 ;; GNU General Public License for more details.
22 ;; You should have received a copy of the GNU General Public License
23 ;; along with GNU Emacs. If not, see <http://www.gnu.org/licenses/>.
27 ;; Calc is split into many files. This file is the main entry point.
28 ;; This file includes autoload commands for various other basic Calc
29 ;; facilities. The more advanced features are based in calc-ext, which
30 ;; in turn contains autoloads for the rest of the Calc files. This
31 ;; odd set of interactions is designed to make Calc's loading time
32 ;; be as short as possible when only simple calculations are needed.
34 ;; Original author's address:
35 ;; Dave Gillespie, daveg@synaptics.com, uunet!synaptx!daveg.
36 ;; Synaptics, Inc., 2698 Orchard Parkway, San Jose, CA 95134.
38 ;; The old address daveg@csvax.cs.caltech.edu will continue to
39 ;; work for the foreseeable future.
41 ;; Bug reports and suggestions are always welcome! (Type M-x
42 ;; report-calc-bug to send them).
44 ;; All functions, macros, and Lisp variables defined here begin with one
45 ;; of the prefixes "math", "Math", or "calc", with the exceptions of
46 ;; "full-calc", "full-calc-keypad", "another-calc", "quick-calc",
47 ;; "report-calc-bug", and "defmath". User-accessible variables begin
52 ;; Fix rewrite mechanism to do less gratuitous rearrangement of terms.
53 ;; Implement a pattern-based "refers" predicate.
55 ;; Make it possible to Undo a selection command.
56 ;; Figure out how to allow selecting rows of matrices.
57 ;; If cursor was in selection before, move it after j n, j p, j L, etc.
58 ;; Consider reimplementing calc-delete-selection using rewrites.
60 ;; Implement line-breaking in non-flat compositions (is this desirable?).
61 ;; Implement matrix formatting with multi-line components.
63 ;; Have "Z R" define a user command based on a set of rewrite rules.
64 ;; Support "incf" and "decf" in defmath definitions.
65 ;; Have defmath generate calls to calc-binary-op or calc-unary-op.
66 ;; Make some way to define algebraic functions using keyboard macros.
68 ;; Allow calc-word-size=0 => Common Lisp-style signed bitwise arithmetic.
69 ;; Consider digamma function (and thus arb. prec. Euler's gamma constant).
70 ;; May as well make continued-fractions stuff available to the user.
72 ;; How about matrix eigenvalues, SVD, pseudo-inverse, etc.?
73 ;; Should cache matrix inverses as well as decompositions.
74 ;; If dividing by a non-square matrix, use least-squares automatically.
75 ;; Consider supporting matrix exponentials.
77 ;; Have ninteg detect and work around singularities at the endpoints.
78 ;; Use an adaptive subdivision algorithm for ninteg.
79 ;; Provide nsum and nprod to go along with ninteg.
81 ;; Handle TeX-mode parsing of \matrix{ ... } where ... contains braces.
82 ;; Support AmS-TeX's \{d,t,}frac, \{d,t,}binom notations.
83 ;; Format and parse sums and products in Eqn and Math modes.
85 ;; Get math-read-big-expr to read sums, products, etc.
86 ;; Change calc-grab-region to use math-read-big-expr.
87 ;; Have a way to define functions using := in Embedded Mode.
89 ;; Support polar plotting with GNUPLOT.
90 ;; Make a calc-graph-histogram function.
92 ;; Replace hokey formulas for complex functions with formulas designed
93 ;; to minimize roundoff while maintaining the proper branch cuts.
94 ;; Test accuracy of advanced math functions over whole complex plane.
95 ;; Extend Bessel functions to provide arbitrary precision.
96 ;; Extend advanced math functions to handle error forms and intervals.
97 ;; Provide a better implementation for math-sin-cos-raw.
98 ;; Provide a better implementation for math-hypot.
99 ;; Provide a better implementation for math-make-frac.
100 ;; Provide a better implementation for calcFunc-prfac.
101 ;; Provide a better implementation for calcFunc-factor.
103 ;; Provide more examples in the tutorial section of the manual.
104 ;; Cover in the tutorial: simplification modes, declarations,
105 ;; bitwise stuff, selections, matrix mapping, financial functions.
106 ;; Provide more Lisp programming examples in the manual.
107 ;; Finish the Internals section of the manual (and bring it up to date).
109 ;; Tim suggests adding spreadsheet-like features.
110 ;; Implement language modes for Gnuplot, Lisp, Ada, APL, ...?
112 ;; For atan series, if x > tan(pi/12) (about 0.268) reduce using the identity
113 ;; atan(x) = atan((x * sqrt(3) - 1) / (sqrt(3) + x)) + pi/6.
115 ;; A better integration algorithm:
116 ;; Use breadth-first instead of depth-first search, as follows:
117 ;; The integral cache allows unfinished integrals in symbolic notation
118 ;; on the righthand side. An entry with no unfinished integrals on the
119 ;; RHS is "complete"; references to it elsewhere are replaced by the
120 ;; integrated value. More than one cache entry for the same integral
121 ;; may exist, though if one becomes complete, the others may be deleted.
122 ;; The integrator works by using every applicable rule (such as
123 ;; substitution, parts, linearity, etc.) to generate possible righthand
124 ;; sides, all of which are entered into the cache. Now, as long as the
125 ;; target integral is not complete (and the time limit has not run out)
126 ;; choose an incomplete integral from the cache and, for every integral
127 ;; appearing in its RHS's, add those integrals to the cache using the
128 ;; same substitition, parts, etc. rules. The cache should be organized
129 ;; as a priority queue, choosing the "simplest" incomplete integral at
130 ;; each step, or choosing randomly among equally simple integrals.
131 ;; Simplicity equals small size, and few steps removed from the original
132 ;; target integral. Note that when the integrator finishes, incomplete
133 ;; integrals can be left in the cache, so the algorithm can start where
134 ;; it left off if another similar integral is later requested.
135 ;; Breadth-first search would avoid the nagging problem of, e.g., whether
136 ;; to use parts or substitution first, and which decomposition is best.
137 ;; All are tried, and any path that diverges will quickly be put on the
138 ;; back burner by the priority queue.
139 ;; Note: Probably a good idea to call math-simplify-extended before
140 ;; measuring a formula's simplicity.
146 ;; Declare functions which are defined elsewhere.
147 (declare-function calc-set-language
"calc-lang" (lang &optional option no-refresh
))
148 (declare-function calc-edit-finish
"calc-yank" (&optional keep
))
149 (declare-function calc-edit-cancel
"calc-yank" ())
150 (declare-function calc-do-quick-calc
"calc-aent" ())
151 (declare-function calc-do-calc-eval
"calc-aent" (str separator args
))
152 (declare-function calc-do-keypad
"calc-keypd" (&optional full-display interactive
))
153 (declare-function calcFunc-unixtime
"calc-forms" (date &optional zone
))
154 (declare-function math-parse-date
"calc-forms" (math-pd-str))
155 (declare-function math-lessp
"calc-ext" (a b
))
156 (declare-function calc-embedded-finish-command
"calc-embed" ())
157 (declare-function calc-embedded-select-buffer
"calc-embed" ())
158 (declare-function calc-embedded-mode-line-change
"calc-embed" ())
159 (declare-function calc-push-list-in-macro
"calc-prog" (vals m sels
))
160 (declare-function calc-replace-selections
"calc-sel" (n vals m
))
161 (declare-function calc-record-list
"calc-misc" (vals &optional prefix
))
162 (declare-function calc-normalize-fancy
"calc-ext" (val))
163 (declare-function calc-do-handle-whys
"calc-misc" ())
164 (declare-function calc-top-selected
"calc-sel" (&optional n m
))
165 (declare-function calc-sel-error
"calc-sel" ())
166 (declare-function calc-pop-stack-in-macro
"calc-prog" (n mm
))
167 (declare-function calc-embedded-stack-change
"calc-embed" ())
168 (declare-function calc-refresh-evaltos
"calc-ext" (&optional which-var
))
169 (declare-function calc-do-refresh
"calc-misc" ())
170 (declare-function calc-binary-op-fancy
"calc-ext" (name func arg ident unary
))
171 (declare-function calc-unary-op-fancy
"calc-ext" (name func arg
))
172 (declare-function calc-delete-selection
"calc-sel" (n))
173 (declare-function calc-alg-digit-entry
"calc-aent" ())
174 (declare-function calc-alg-entry
"calc-aent" (&optional initial prompt
))
175 (declare-function calc-dots
"calc-incom" ())
176 (declare-function calc-temp-minibuffer-message
"calc-misc" (m))
177 (declare-function math-read-radix-digit
"calc-misc" (dig))
178 (declare-function calc-digit-dots
"calc-incom" ())
179 (declare-function math-normalize-fancy
"calc-ext" (a))
180 (declare-function math-normalize-nonstandard
"calc-ext" ())
181 (declare-function math-recompile-eval-rules
"calc-alg" ())
182 (declare-function math-apply-rewrites
"calc-rewr" (expr rules
&optional heads math-apply-rw-ruleset
))
183 (declare-function calc-record-why
"calc-misc" (&rest stuff
))
184 (declare-function math-dimension-error
"calc-vec" ())
185 (declare-function calc-incomplete-error
"calc-incom" (a))
186 (declare-function math-float-fancy
"calc-arith" (a))
187 (declare-function math-neg-fancy
"calc-arith" (a))
188 (declare-function math-zerop
"calc-misc" (a))
189 (declare-function calc-add-fractions
"calc-frac" (a b
))
190 (declare-function math-add-objects-fancy
"calc-arith" (a b
))
191 (declare-function math-add-symb-fancy
"calc-arith" (a b
))
192 (declare-function math-mul-zero
"calc-arith" (a b
))
193 (declare-function calc-mul-fractions
"calc-frac" (a b
))
194 (declare-function math-mul-objects-fancy
"calc-arith" (a b
))
195 (declare-function math-mul-symb-fancy
"calc-arith" (a b
))
196 (declare-function math-reject-arg
"calc-misc" (&optional a p option
))
197 (declare-function math-div-by-zero
"calc-arith" (a b
))
198 (declare-function math-div-zero
"calc-arith" (a b
))
199 (declare-function math-make-frac
"calc-frac" (num den
))
200 (declare-function calc-div-fractions
"calc-frac" (a b
))
201 (declare-function math-div-objects-fancy
"calc-arith" (a b
))
202 (declare-function math-div-symb-fancy
"calc-arith" (a b
))
203 (declare-function math-compose-expr
"calccomp" (a prec
))
204 (declare-function math-comp-width
"calccomp" (c))
205 (declare-function math-composition-to-string
"calccomp" (c &optional width
))
206 (declare-function math-stack-value-offset-fancy
"calccomp" ())
207 (declare-function math-format-flat-expr-fancy
"calc-ext" (a prec
))
208 (declare-function math-adjust-fraction
"calc-ext" (a))
209 (declare-function math-format-binary
"calc-bin" (a))
210 (declare-function math-format-radix
"calc-bin" (a))
211 (declare-function math-group-float
"calc-ext" (str))
212 (declare-function math-mod
"calc-misc" (a b
))
213 (declare-function math-format-number-fancy
"calc-ext" (a prec
))
214 (declare-function math-format-bignum-fancy
"calc-ext" (a))
215 (declare-function math-read-number-fancy
"calc-ext" (s))
216 (declare-function calc-do-grab-region
"calc-yank" (top bot arg
))
217 (declare-function calc-do-grab-rectangle
"calc-yank" (top bot arg
&optional reduce
))
218 (declare-function calc-do-embedded
"calc-embed" (calc-embed-arg end obeg oend
))
219 (declare-function calc-do-embedded-activate
"calc-embed" (calc-embed-arg cbuf
))
220 (declare-function math-do-defmath
"calc-prog" (func args body
))
221 (declare-function calc-load-everything
"calc-ext" ())
228 :group
'applications
)
231 (defcustom calc-settings-file
232 (convert-standard-filename "~/.calc.el")
233 "File in which to record permanent settings."
237 (defcustom calc-language-alist
238 '((latex-mode . latex
)
240 (plain-tex-mode . tex
)
243 (pascal-mode . pascal
)
246 (fortran-mode . fortran
)
248 (texinfo-mode . calc-normal-language
))
249 "Alist of major modes with appropriate Calc languages."
251 :type
'(alist :key-type
(symbol :tag
"Major mode")
252 :value-type
(symbol :tag
"Calc language")))
254 (defcustom calc-embedded-announce-formula
255 "%Embed\n\\(% .*\n\\)*"
256 "A regular expression which is sure to be followed by a calc-embedded formula."
260 (defcustom calc-embedded-announce-formula-alist
261 '((c++-mode .
"//Embed\n\\(// .*\n\\)*")
262 (c-mode .
"/\\*Embed\\*/\n\\(/\\* .*\\*/\n\\)*")
263 (f90-mode .
"!Embed\n\\(! .*\n\\)*")
264 (fortran-mode .
"C Embed\n\\(C .*\n\\)*")
265 (html-helper-mode .
"<!-- Embed -->\n\\(<!-- .* -->\n\\)*")
266 (html-mode .
"<!-- Embed -->\n\\(<!-- .* -->\n\\)*")
267 (nroff-mode .
"\\\\\"Embed\n\\(\\\\\" .*\n\\)*")
268 (pascal-mode .
"{Embed}\n\\({.*}\n\\)*")
269 (sgml-mode .
"<!-- Embed -->\n\\(<!-- .* -->\n\\)*")
270 (xml-mode .
"<!-- Embed -->\n\\(<!-- .* -->\n\\)*")
271 (texinfo-mode .
"@c Embed\n\\(@c .*\n\\)*"))
272 "Alist of major modes with appropriate values for `calc-embedded-announce-formula'."
274 :type
'(alist :key-type
(symbol :tag
"Major mode")
275 :value-type
(regexp :tag
"Regexp to announce formula")))
277 (defcustom calc-embedded-open-formula
278 "\\`\\|^\n\\|\\$\\$?\\|\\\\\\[\\|^\\\\begin[^{].*\n\\|^\\\\begin{.*[^x]}.*\n\\|^@.*\n\\|^\\.EQ.*\n\\|\\\\(\\|^%\n\\|^\\.\\\\\"\n"
279 "A regular expression for the opening delimiter of a formula used by calc-embedded."
283 (defcustom calc-embedded-close-formula
284 "\\'\\|\n$\\|\\$\\$?\\|\\\\]\\|^\\\\end[^{].*\n\\|^\\\\end{.*[^x]}.*\n\\|^@.*\n\\|^\\.EN.*\n\\|\\\\)\\|\n%\n\\|^\\.\\\\\"\n"
285 "A regular expression for the closing delimiter of a formula used by calc-embedded."
289 (defcustom calc-embedded-open-close-formula-alist
291 "Alist of major modes with pairs of formula delimiters used by calc-embedded."
293 :type
'(alist :key-type
(symbol :tag
"Major mode")
294 :value-type
(list (regexp :tag
"Opening formula delimiter")
295 (regexp :tag
"Closing formula delimiter"))))
297 (defcustom calc-embedded-word-regexp
298 "[-+]?[0-9]+\\(\\.[0-9]+\\)?\\([eE][-+]?[0-9]+\\)?"
299 "A regular expression determining a word for calc-embedded-word."
303 (defcustom calc-embedded-word-regexp-alist
305 "Alist of major modes with word regexps used by calc-embedded-word."
307 :type
'(alist :key-type
(symbol :tag
"Major mode")
308 :value-type
(regexp :tag
"Regexp for word")))
310 (defcustom calc-embedded-open-plain
312 "A string which is the opening delimiter for a \"plain\" formula.
313 If calc-show-plain mode is enabled, this is inserted at the front of
318 (defcustom calc-embedded-close-plain
320 "A string which is the closing delimiter for a \"plain\" formula.
321 See calc-embedded-open-plain."
325 (defcustom calc-embedded-open-close-plain-alist
326 '((c++-mode
"// %% " " %%\n")
327 (c-mode "/* %% " " %% */\n")
328 (f90-mode "! %% " " %%\n")
329 (fortran-mode "C %% " " %%\n")
330 (html-helper-mode "<!-- %% " " %% -->\n")
331 (html-mode "<!-- %% " " %% -->\n")
332 (nroff-mode "\\\" %% " " %%\n")
333 (pascal-mode "{%% " " %%}\n")
334 (sgml-mode "<!-- %% " " %% -->\n")
335 (xml-mode "<!-- %% " " %% -->\n")
336 (texinfo-mode "@c %% " " %%\n"))
337 "Alist of major modes with pairs of delimiters for \"plain\" formulas."
339 :type
'(alist :key-type
(symbol :tag
"Major mode")
340 :value-type
(list (string :tag
"Opening \"plain\" delimiter")
341 (string :tag
"Closing \"plain\" delimiter"))))
343 (defcustom calc-embedded-open-new-formula
345 "A string which is inserted at front of formula by calc-embedded-new-formula."
349 (defcustom calc-embedded-close-new-formula
351 "A string which is inserted at end of formula by calc-embedded-new-formula."
355 (defcustom calc-embedded-open-close-new-formula-alist
357 "Alist of major modes with pairs of new formula delimiters used by calc-embedded."
359 :type
'(alist :key-type
(symbol :tag
"Major mode")
360 :value-type
(list (string :tag
"Opening new formula delimiter")
361 (string :tag
"Closing new formula delimiter"))))
363 (defcustom calc-embedded-open-mode
365 "A string which should precede calc-embedded mode annotations.
366 This is not required to be present for user-written mode annotations."
370 (defcustom calc-embedded-close-mode
372 "A string which should follow calc-embedded mode annotations.
373 This is not required to be present for user-written mode annotations."
377 (defcustom calc-embedded-open-close-mode-alist
378 '((c++-mode
"// " "\n")
379 (c-mode "/* " " */\n")
381 (fortran-mode "C " "\n")
382 (html-helper-mode "<!-- " " -->\n")
383 (html-mode "<!-- " " -->\n")
384 (nroff-mode "\\\" " "\n")
385 (pascal-mode "{ " " }\n")
386 (sgml-mode "<!-- " " -->\n")
387 (xml-mode "<!-- " " -->\n")
388 (texinfo-mode "@c " "\n"))
389 "Alist of major modes with pairs of strings to delimit annotations."
391 :type
'(alist :key-type
(symbol :tag
"Major mode")
392 :value-type
(list (string :tag
"Opening annotation delimiter")
393 (string :tag
"Closing annotation delimiter"))))
395 (defcustom calc-gnuplot-name
396 (if (eq system-type
'windows-nt
) "pgnuplot" "gnuplot")
397 "Name of GNUPLOT program, for calc-graph features."
401 (defcustom calc-gnuplot-plot-command
403 "Name of command for displaying GNUPLOT output; %s = file name to print."
405 :type
'(choice (string) (sexp)))
407 (defcustom calc-gnuplot-print-command
409 "Name of command for printing GNUPLOT output; %s = file name to print."
411 :type
'(choice (string) (sexp)))
413 (defcustom calc-multiplication-has-precedence
415 "If non-nil, multiplication has precedence over division
420 (defcustom calc-undo-length
422 "The number of undo steps that will be preserved when Calc is quit."
426 (defvar calc-bug-address
"jay.p.belanger@gmail.com"
427 "Address of the maintainer of Calc, for use by `report-calc-bug'.")
429 (defvar calc-scan-for-dels t
430 "If t, scan keymaps to find all DEL-like keys.
431 if nil, only DEL itself is mapped to calc-pop.")
433 (defvar calc-stack
'((top-of-stack 1 nil
))
435 Entries are 3-lists: Formula, Height (in lines), Selection (or nil).")
437 (defvar calc-stack-top
1
438 "Index into `calc-stack' of \"top\" of stack.
439 This is 1 unless `calc-truncate-stack' has been used.")
441 (defvar calc-display-sci-high
0
442 "Floating-point numbers with this positive exponent or higher above the
443 current precision are displayed in scientific notation in calc-mode.")
445 (defvar calc-display-sci-low -
3
446 "Floating-point numbers with this negative exponent or lower are displayed
447 scientific notation in calc-mode.")
449 (defvar calc-other-modes nil
450 "List of used-defined strings to append to Calculator mode line.")
452 (defvar calc-Y-help-msgs nil
453 "List of strings for Y prefix help.")
455 (defvar calc-loaded-settings-file nil
456 "t if `calc-settings-file' has been loaded yet.")
459 (defvar calc-mode-var-list
'()
460 "List of variables used in customizing GNU Calc.")
462 (defmacro defcalcmodevar
(var defval
&optional doc
)
463 "Declare VAR as a Calc variable, with default value DEFVAL
465 The variable VAR will be added to `calc-mode-var-list'."
467 (defvar ,var
,defval
,doc
)
468 (add-to-list 'calc-mode-var-list
(list (quote ,var
) ,defval
))))
470 (defun calc-mode-var-list-restore-default-values ()
471 "Restore the default values of the variables in `calc-mode-var-list'."
472 (mapcar (function (lambda (v) (set (car v
) (nth 1 v
))))
475 (defun calc-mode-var-list-restore-saved-values ()
476 "Restore the user-saved values of the variables in `calc-mode-var-list'."
477 (let ((newvarlist '()))
480 (file (substitute-in-file-name calc-settings-file
)))
482 (file-regular-p file
)
483 (set-buffer (find-file-noselect file
))
484 (goto-char (point-min))
485 (search-forward ";;; Mode settings stored by Calc" nil t
)
489 (search-forward "\n;;; End of mode settings" nil t
)))
491 (calc-mode-var-list-restore-default-values)
492 (eval-region pos
(point))
493 (let ((varlist calc-mode-var-list
))
495 (let ((var (car varlist
)))
497 (cons (list (car var
) (symbol-value (car var
)))
499 (setq varlist
(cdr varlist
)))))))
501 (mapcar (function (lambda (v) (set (car v
) (nth 1 v
))))
503 (calc-mode-var-list-restore-default-values))))
505 (defcalcmodevar calc-always-load-extensions nil
506 "If non-nil, load the calc-ext module automatically when Calc is loaded.")
508 (defcalcmodevar calc-line-numbering t
509 "If non-nil, display line numbers in Calculator stack.")
511 (defcalcmodevar calc-line-breaking t
512 "If non-nil, break long values across multiple lines in Calculator stack.")
514 (defcalcmodevar calc-display-just nil
515 "If nil, stack display is left-justified.
516 If `right', stack display is right-justified.
517 If `center', stack display is centered.")
519 (defcalcmodevar calc-display-origin nil
520 "Horizontal origin of displayed stack entries.
521 In left-justified mode, this is effectively indentation. (Default 0).
522 In right-justified mode, this is effectively window width.
523 In centered mode, center of stack entry is placed here.")
525 (defcalcmodevar calc-number-radix
10
526 "Radix for entry and display of numbers in calc-mode, 2-36.")
528 (defcalcmodevar calc-leading-zeros nil
529 "If non-nil, leading zeros are provided to pad integers to calc-word-size.")
531 (defcalcmodevar calc-group-digits nil
532 "If non-nil, group digits in large displayed integers by inserting spaces.
533 If an integer, group that many digits at a time.
534 If t, use 4 for binary and hex, 3 otherwise.")
536 (defcalcmodevar calc-group-char
","
537 "The character (in the form of a string) to be used for grouping digits.
538 This is used only when calc-group-digits mode is on.")
540 (defcalcmodevar calc-point-char
"."
541 "The character (in the form of a string) to be used as a decimal point.")
543 (defcalcmodevar calc-frac-format
'(":" nil
)
544 "Format of displayed fractions; a string of one or two of \":\" or \"/\".")
546 (defcalcmodevar calc-prefer-frac nil
547 "If non-nil, prefer fractional over floating-point results.")
549 (defcalcmodevar calc-hms-format
"%s@ %s' %s\""
550 "Format of displayed hours-minutes-seconds angles, a format string.
551 String must contain three %s marks for hours, minutes, seconds respectively.")
553 (defcalcmodevar calc-date-format
'((H ":" mm C SS pp
" ")
554 Www
" " Mmm
" " D
", " YYYY
)
555 "Format of displayed date forms.")
557 (defcalcmodevar calc-float-format
'(float 0)
558 "Format to use for display of floating-point numbers in calc-mode.
559 Must be a list of one of the following forms:
560 (float 0) Floating point format, display full precision.
561 (float N) N > 0: Floating point format, at most N significant figures.
562 (float -N) -N < 0: Floating point format, calc-internal-prec - N figs.
563 (fix N) N >= 0: Fixed point format, N places after decimal point.
564 (sci 0) Scientific notation, full precision.
565 (sci N) N > 0: Scientific notation, N significant figures.
566 (sci -N) -N < 0: Scientific notation, calc-internal-prec - N figs.
567 (eng 0) Engineering notation, full precision.
568 (eng N) N > 0: Engineering notation, N significant figures.
569 (eng -N) -N < 0: Engineering notation, calc-internal-prec - N figs.")
571 (defcalcmodevar calc-full-float-format
'(float 0)
572 "Format to use when full precision must be displayed.")
574 (defcalcmodevar calc-complex-format nil
575 "Format to use for display of complex numbers in calc-mode. Must be one of:
580 (defcalcmodevar calc-complex-mode
'cplx
581 "Preferred form, either `cplx' or `polar', for complex numbers.")
583 (defcalcmodevar calc-infinite-mode nil
584 "If nil, 1 / 0 is left unsimplified.
585 If 0, 1 / 0 is changed to inf (zeros are considered positive).
586 Otherwise, 1 / 0 is changed to uinf (undirected infinity).")
588 (defcalcmodevar calc-display-strings nil
589 "If non-nil, display vectors of byte-sized integers as strings.")
591 (defcalcmodevar calc-matrix-just
'center
592 "If nil, vector elements are left-justified.
593 If `right', vector elements are right-justified.
594 If `center', vector elements are centered.")
596 (defcalcmodevar calc-break-vectors nil
597 "If non-nil, display vectors one element per line.")
599 (defcalcmodevar calc-full-vectors t
600 "If non-nil, display long vectors in full. If nil, use abbreviated form.")
602 (defcalcmodevar calc-full-trail-vectors t
603 "If non-nil, display long vectors in full in the trail.")
605 (defcalcmodevar calc-vector-commas
","
606 "If non-nil, separate elements of displayed vectors with this string.")
608 (defcalcmodevar calc-vector-brackets
"[]"
609 "If non-nil, surround displayed vectors with these characters.")
611 (defcalcmodevar calc-matrix-brackets
'(R O
)
612 "A list of code-letter symbols that control \"big\" matrix display.
613 If `R' is present, display inner brackets for matrices.
614 If `O' is present, display outer brackets for matrices (above/below).
615 If `C' is present, display outer brackets for matrices (centered).")
617 (defcalcmodevar calc-language nil
618 "Language or format for entry and display of stack values. Must be one of:
619 nil Use standard Calc notation.
620 flat Use standard Calc notation, one-line format.
621 big Display formulas in 2-d notation (enter w/std notation).
622 unform Use unformatted display: add(a, mul(b,c)).
623 c Use C language notation.
624 pascal Use Pascal language notation.
625 fortran Use Fortran language notation.
626 tex Use TeX notation.
627 latex Use LaTeX notation.
628 eqn Use eqn notation.
629 yacas Use Yacas notation.
630 maxima Use Maxima notation.
631 giac Use Giac notation.
632 math Use Mathematica(tm) notation.
633 maple Use Maple notation.")
635 (defcalcmodevar calc-language-option nil
636 "Numeric prefix argument for the command that set `calc-language'.")
638 (defcalcmodevar calc-left-label
""
639 "Label to display at left of formula.")
641 (defcalcmodevar calc-right-label
""
642 "Label to display at right of formula.")
644 (defcalcmodevar calc-word-size
32
645 "Minimum number of bits per word, if any, for binary operations in calc-mode.")
647 (defcalcmodevar calc-previous-modulo nil
648 "Most recently used value of M in a modulo form.")
650 (defcalcmodevar calc-simplify-mode nil
651 "Type of simplification applied to results.
652 If `none', results are not simplified when pushed on the stack.
653 If `num', functions are simplified only when args are constant.
654 If nil, only fast simplifications are applied.
655 If `binary', `math-clip' is applied if appropriate.
656 If `alg', `math-simplify' is applied.
657 If `ext', `math-simplify-extended' is applied.
658 If `units', `math-simplify-units' is applied.")
660 (defcalcmodevar calc-auto-recompute t
661 "If non-nil, recompute evalto's automatically when necessary.")
663 (defcalcmodevar calc-display-raw nil
664 "If non-nil, display shows unformatted Lisp exprs. (For debugging)")
666 (defcalcmodevar calc-internal-prec
12
667 "Number of digits of internal precision for calc-mode calculations.")
669 (defcalcmodevar calc-angle-mode
'deg
670 "If deg, angles are in degrees; if rad, angles are in radians.
671 If hms, angles are in degrees-minutes-seconds.")
673 (defcalcmodevar calc-algebraic-mode nil
674 "If non-nil, numeric entry accepts whole algebraic expressions.
675 If nil, algebraic expressions must be preceded by \"'\".")
677 (defcalcmodevar calc-incomplete-algebraic-mode nil
678 "Like calc-algebraic-mode except only affects ( and [ keys.")
680 (defcalcmodevar calc-symbolic-mode nil
681 "If non-nil, inexact numeric computations like sqrt(2) are postponed.
682 If nil, computations on numbers always yield numbers where possible.")
684 (defcalcmodevar calc-matrix-mode nil
685 "If `matrix', variables are assumed to be matrix-valued.
686 If a number, variables are assumed to be NxN matrices.
687 If `sqmatrix', variables are assumed to be square matrices of an unspecified size.
688 If `scalar', variables are assumed to be scalar-valued.
689 If nil, symbolic math routines make no assumptions about variables.")
691 (defcalcmodevar calc-shift-prefix nil
692 "If non-nil, shifted letter keys are prefix keys rather than normal meanings.")
694 (defcalcmodevar calc-window-height
7
695 "Initial height of Calculator window.")
697 (defcalcmodevar calc-display-trail t
698 "If non-nil, M-x calc creates a window to display Calculator trail.")
700 (defcalcmodevar calc-show-selections t
701 "If non-nil, selected sub-formulas are shown by obscuring rest of formula.
702 If nil, selected sub-formulas are highlighted by obscuring the sub-formulas.")
704 (defcalcmodevar calc-use-selections t
705 "If non-nil, commands operate only on selected portions of formulas.
706 If nil, selections displayed but ignored.")
708 (defcalcmodevar calc-assoc-selections t
709 "If non-nil, selection hides deep structure of associative formulas.")
711 (defcalcmodevar calc-display-working-message
'lots
712 "If non-nil, display \"Working...\" for potentially slow Calculator commands.")
714 (defcalcmodevar calc-auto-why
'maybe
715 "If non-nil, automatically execute a \"why\" command to explain odd results.")
717 (defcalcmodevar calc-timing nil
718 "If non-nil, display timing information on each slow command.")
720 (defcalcmodevar calc-mode-save-mode
'local
)
722 (defcalcmodevar calc-standard-date-formats
724 "<H:mm:SSpp >Www Mmm D, YYYY"
725 "D Mmm YYYY<, h:mm:SS>"
726 "Www Mmm BD< hh:mm:ss> YYYY"
734 (defcalcmodevar calc-autorange-units nil
735 "If non-nil, automatically set unit prefixes to keep units in a reasonable range.")
737 (defcalcmodevar calc-was-keypad-mode nil
738 "Non-nil if Calc was last invoked in keypad mode.")
740 (defcalcmodevar calc-full-mode nil
741 "Non-nil if Calc was last invoked in full-screen mode.")
743 (defcalcmodevar calc-user-parse-tables nil
744 "Alist of languages with user-defined parse rules.")
746 (defcalcmodevar calc-gnuplot-default-device
"default"
747 "The default device name for GNUPLOT plotting.")
749 (defcalcmodevar calc-gnuplot-default-output
"STDOUT"
750 "The default output file for GNUPLOT plotting.")
752 (defcalcmodevar calc-gnuplot-print-device
"postscript"
753 "The default device name for GNUPLOT printing.")
755 (defcalcmodevar calc-gnuplot-print-output
"auto"
756 "The default output for GNUPLOT printing.")
758 (defcalcmodevar calc-gnuplot-geometry nil
759 "The default geometry for the GNUPLOT window.")
761 (defcalcmodevar calc-graph-default-resolution
15
762 "The default number of data points when plotting curves.")
764 (defcalcmodevar calc-graph-default-resolution-3d
5
765 "The default number of x- and y- data points when plotting surfaces.")
767 (defcalcmodevar calc-invocation-macro nil
768 "A user defined macro for starting Calc.
769 Used by `calc-user-invocation'.")
771 (defcalcmodevar calc-show-banner t
772 "*If non-nil, show a friendly greeting above the stack.")
774 (defconst calc-local-var-list
'(calc-stack
778 calc-always-load-extensions
789 calc-incomplete-algebraic-mode
805 calc-standard-date-formats
807 calc-full-float-format
811 calc-full-trail-vectors
825 calc-assoc-selections
829 (defvar calc-mode-hook nil
830 "Hook run when entering calc-mode.")
832 (defvar calc-trail-mode-hook nil
833 "Hook run when entering calc-trail-mode.")
835 (defvar calc-start-hook nil
836 "Hook run when calc is started.")
838 (defvar calc-end-hook nil
839 "Hook run when calc is quit.")
841 (defvar calc-load-hook nil
842 "Hook run when calc.el is loaded.")
844 (defvar calc-window-hook nil
845 "Hook called to create the Calc window.")
847 (defvar calc-trail-window-hook nil
848 "Hook called to create the Calc trail window.")
850 (defvar calc-embedded-new-buffer-hook nil
851 "Hook run when starting embedded mode in a new buffer.")
853 (defvar calc-embedded-new-formula-hook nil
854 "Hook run when starting embedded mode in a new formula.")
856 (defvar calc-embedded-mode-hook nil
857 "Hook run when starting embedded mode.")
859 ;; Set up the autoloading linkage.
860 (let ((name (and (fboundp 'calc-dispatch
)
861 (eq (car-safe (symbol-function 'calc-dispatch
)) 'autoload
)
862 (nth 1 (symbol-function 'calc-dispatch
))))
865 ;; If Calc files exist on the load-path, we're all set.
866 (while (and p
(not (file-exists-p
867 (expand-file-name "calc-misc.elc" (car p
)))))
871 ;; If Calc is autoloaded using a path name, look there for Calc files.
872 ;; This works for both relative ("calc/calc.elc") and absolute paths.
873 (and name
(file-name-directory name
)
875 (name2 (concat (file-name-directory name
)
877 (while (and p2
(not (file-exists-p
878 (expand-file-name name2
(car p2
)))))
881 (setq load-path
(nconc load-path
886 name
(car p2
))))))))))))
888 ;; The following modes use specially-formatted data.
889 (put 'calc-mode
'mode-class
'special
)
890 (put 'calc-trail-mode
'mode-class
'special
)
892 ;; Define "inexact-result" as an e-lisp error symbol.
893 (put 'inexact-result
'error-conditions
'(error inexact-result calc-error
))
894 (put 'inexact-result
'error-message
"Calc internal error (inexact-result)")
896 ;; Define "math-overflow" and "math-underflow" as e-lisp error symbols.
897 (put 'math-overflow
'error-conditions
'(error math-overflow calc-error
))
898 (put 'math-overflow
'error-message
"Floating-point overflow occurred")
899 (put 'math-underflow
'error-conditions
'(error math-underflow calc-error
))
900 (put 'math-underflow
'error-message
"Floating-point underflow occurred")
902 (defvar calc-trail-pointer nil
903 "The \"current\" entry in trail buffer.")
904 (defvar calc-trail-overlay nil
905 "The value of overlay-arrow-string.")
906 (defvar calc-undo-list nil
907 "The list of previous operations for undo.")
908 (defvar calc-redo-list nil
909 "The list of recent undo operations.")
910 (defvar calc-main-buffer nil
911 "A pointer to Calculator buffer.")
912 (defvar calc-buffer-list nil
913 "A list of all Calc buffers.")
914 (defvar calc-trail-buffer nil
915 "A pointer to Calc Trail buffer.")
917 "Explanations of most recent errors.")
918 (defvar calc-next-why nil
)
919 (defvar calc-inverse-flag nil
920 "If non-nil, next operation is Inverse.")
921 (defvar calc-hyperbolic-flag nil
922 "If non-nil, next operation is Hyperbolic.")
923 (defvar calc-keep-args-flag nil
924 "If non-nil, next operation should not remove its arguments from stack.")
925 (defvar calc-function-open
"("
926 "Open-parenthesis string for function call notation.")
927 (defvar calc-function-close
")"
928 "Close-parenthesis string for function call notation.")
929 (defvar calc-language-output-filter nil
930 "Function through which to pass strings after formatting.")
931 (defvar calc-language-input-filter nil
932 "Function through which to pass strings before parsing.")
933 (defvar calc-radix-formatter nil
934 "Formatting function used for non-decimal numbers.")
935 (defvar calc-lang-slash-idiv nil
936 "A list of languages in which / might represent integer division.")
937 (defvar calc-lang-allow-underscores nil
938 "A list of languages which allow underscores in variable names.")
939 (defvar calc-lang-allow-percentsigns nil
940 "A list of languages which allow percent signs in variable names.")
941 (defvar calc-lang-c-type-hex nil
942 "Languages in which octal and hex numbers are written with leading 0 and 0x,")
943 (defvar calc-lang-brackets-are-subscripts nil
944 "Languages in which subscripts are indicated by brackets.")
945 (defvar calc-lang-parens-are-subscripts nil
946 "Languages in which subscripts are indicated by parentheses.")
948 (defvar calc-last-kill nil
949 "The last number killed in calc-mode.")
950 (defvar calc-dollar-values nil
951 "Values to be used for '$'.")
952 (defvar calc-dollar-used nil
953 "The highest order of '$' that occurred.")
954 (defvar calc-hashes-used nil
955 "The highest order of '#' that occurred.")
956 (defvar calc-quick-prev-results nil
957 "Previous results from Quick Calc.")
958 (defvar calc-said-hello nil
959 "Non-nil if the welcomd message has been displayed.")
960 (defvar calc-executing-macro nil
961 "Non-nil if a keyboard macro is executing from the \"K\" key.")
962 (defvar calc-any-selections nil
963 "Non-nil if there are selections present.")
964 (defvar calc-help-phase
0
965 "The number of consecutive \"?\" keystrokes.")
966 (defvar calc-full-help-flag nil
967 "Non-nil if `calc-full-help' is being executed.")
968 (defvar calc-refresh-count
0
969 "The number of `calc-refresh' calls.")
970 (defvar calc-display-dirty nil
971 "Non-nil if the stack display might not reflect the latest mode settings.")
972 (defvar calc-prepared-composition nil
)
973 (defvar calc-selection-cache-default-entry nil
)
974 (defvar calc-embedded-info nil
975 "If non-nil, a vector consisting of information for embedded mode.")
976 (defvar calc-embedded-active nil
977 "Alist of buffers with sorted lists of calc-embedded-infos.")
978 (defvar calc-standalone-flag nil
979 "Non-nil if Emacs started with standalone Calc.")
980 (defvar var-EvalRules nil
981 "User defined rules that Calc will apply automatically.")
982 (defvar math-eval-rules-cache-tag t
)
983 (defvar math-radix-explicit-format t
)
984 (defvar math-expr-function-mapping nil
985 "Alist of language specific functions with Calc functions.")
986 (defvar math-expr-variable-mapping nil
987 "Alist of language specific variables with Calc variables.")
988 (defvar math-read-expr-quotes nil
)
989 (defvar math-working-step nil
)
990 (defvar math-working-step-2 nil
)
991 (defvar var-i
'(special-const (math-imaginary 1)))
992 (defvar var-pi
'(special-const (math-pi)))
993 (defvar var-e
'(special-const (math-e)))
994 (defvar var-phi
'(special-const (math-phi)))
995 (defvar var-gamma
'(special-const (math-gamma-const)))
996 (defvar var-Modes
'(special-const (math-get-modes-vec)))
998 (mapc (lambda (v) (or (boundp v
) (set v nil
)))
1001 (defvar calc-mode-map
1002 (let ((map (make-keymap)))
1003 (suppress-keymap map t
)
1004 (define-key map
"+" 'calc-plus
)
1005 (define-key map
"-" 'calc-minus
)
1006 (define-key map
"*" 'calc-times
)
1007 (define-key map
"/" 'calc-divide
)
1008 (define-key map
"%" 'calc-mod
)
1009 (define-key map
"&" 'calc-inv
)
1010 (define-key map
"^" 'calc-power
)
1011 (define-key map
"\M-%" 'calc-percent
)
1012 (define-key map
"e" 'calcDigit-start
)
1013 (define-key map
"i" 'calc-info
)
1014 (define-key map
"n" 'calc-change-sign
)
1015 (define-key map
"q" 'calc-quit
)
1016 (define-key map
"Y" 'nil
)
1017 (define-key map
"Y?" 'calc-shift-Y-prefix-help
)
1018 (define-key map
"?" 'calc-help
)
1019 (define-key map
" " 'calc-enter
)
1020 (define-key map
"'" 'calc-algebraic-entry
)
1021 (define-key map
"$" 'calc-auto-algebraic-entry
)
1022 (define-key map
"\"" 'calc-auto-algebraic-entry
)
1023 (define-key map
"\t" 'calc-roll-down
)
1024 (define-key map
"\M-\t" 'calc-roll-up
)
1025 (define-key map
"\C-x\C-t" 'calc-transpose-lines
)
1026 (define-key map
"\C-m" 'calc-enter
)
1027 (define-key map
"\M-\C-m" 'calc-last-args-stub
)
1028 (define-key map
"\C-j" 'calc-over
)
1029 (define-key map
"\C-y" 'calc-yank
)
1030 (define-key map
[mouse-2
] 'calc-yank
)
1032 (mapc (lambda (x) (define-key map
(char-to-string x
) 'undefined
))
1034 (mapc (lambda (x) (define-key map
(char-to-string x
) 'calc-missing-key
))
1035 (concat "ABCDEFGHIJKLMNPQRSTUVXZabcdfghjkmoprstuvwxyz"
1036 ":\\|!()[]<>{},;=~`\C-k\C-w\C-_"))
1037 (define-key map
"\M-w" 'calc-missing-key
)
1038 (define-key map
"\M-k" 'calc-missing-key
)
1039 (define-key map
"\M-\C-w" 'calc-missing-key
)
1040 (mapc (lambda (x) (define-key map
(char-to-string x
) 'calcDigit-start
))
1043 "The key map for Calc.")
1047 (defvar calc-digit-map
1048 (let ((map (make-keymap)))
1049 (if (featurep 'xemacs
)
1050 (map-keymap (function
1052 (define-key map keys
1053 (if (eq bind
'undefined
)
1054 'undefined
'calcDigit-nondigit
))))
1056 (let ((cmap (nth 1 calc-mode-map
))
1061 (if (eq (aref cmap i
) 'undefined
)
1062 'undefined
'calcDigit-nondigit
))
1064 (mapc (lambda (x) (define-key map
(char-to-string x
) 'calcDigit-key
))
1065 "_0123456789.e+-:n#@oh'\"mspM")
1066 (mapc (lambda (x) (define-key map
(char-to-string x
) 'calcDigit-letter
))
1067 "abcdfgijklqrtuvwxyzABCDEFGHIJKLNOPQRSTUVWXYZ")
1068 (define-key map
"'" 'calcDigit-algebraic
)
1069 (define-key map
"`" 'calcDigit-edit
)
1070 (define-key map
"\C-g" 'abort-recursive-edit
)
1072 "The key map for entering Calc digits.")
1077 (define-key calc-digit-map x
'calcDigit-backspace
)
1078 (define-key calc-mode-map x
'calc-pop
)
1079 (define-key calc-mode-map
1081 (if (featurep 'xemacs
)
1082 (if (= (length x
) 1)
1083 (vector (if (consp (aref x
0))
1084 (cons 'meta
(aref x
0))
1085 (list 'meta
(aref x
0))))
1091 (if calc-scan-for-dels
1092 (append (where-is-internal 'delete-backward-char global-map
)
1093 (where-is-internal 'backward-delete-char global-map
)
1094 (where-is-internal 'backward-delete-char-untabify global-map
)
1098 (defvar calc-dispatch-map
1099 (let ((map (make-keymap)))
1101 (let* ((x-chr (car x
))
1102 (x-str (char-to-string x-chr
))
1104 (define-key map x-str x-def
)
1105 (when (string-match "[a-z]" x-str
)
1106 ;; Map upper case char to same definition.
1107 (define-key map
(upcase x-str
) x-def
)
1108 (unless (string-match "[gmv]" x-str
)
1109 ;; Map control prefixed char to same definition.
1110 (define-key map
(vector (list 'control x-chr
)) x-def
)))
1111 (define-key map
(format "\e%c" x-chr
) x-def
)))
1112 '( ( ?a . calc-embedded-activate
)
1113 ( ?b . calc-big-or-small
)
1115 ( ?d . calc-embedded-duplicate
)
1116 ( ?e . calc-embedded
)
1117 ( ?f . calc-embedded-new-formula
)
1118 ( ?g . calc-grab-region
)
1119 ( ?h . calc-dispatch-help
)
1121 ( ?j . calc-embedded-select
)
1122 ( ?k . calc-keypad
)
1123 ( ?l . calc-load-everything
)
1124 ( ?m . read-kbd-macro
)
1125 ( ?n . calc-embedded-next
)
1126 ( ?o . calc-other-window
)
1127 ( ?p . calc-embedded-previous
)
1129 ( ?r . calc-grab-rectangle
)
1130 ( ?s . calc-info-summary
)
1131 ( ?t . calc-tutorial
)
1132 ( ?u . calc-embedded-update-formula
)
1133 ( ?w . calc-embedded-word
)
1135 ( ?y . calc-copy-to-buffer
)
1136 ( ?z . calc-user-invocation
)
1137 ( ?
\' . calc-embedded-new-formula
)
1138 ( ?\
` . calc-embedded-edit
)
1139 ( ?
: . calc-grab-sum-down
)
1140 ( ?_ . calc-grab-sum-across
)
1142 ( ?? . calc-dispatch-help
)
1143 ( ?
# . calc-same-interface
)
1144 ( ?
& . calc-same-interface
)
1145 ( ?
\\ . calc-same-interface
)
1146 ( ?
= . calc-same-interface
)
1147 ( ?
* . calc-same-interface
)
1148 ( ?
/ . calc-same-interface
)
1149 ( ?
+ . calc-same-interface
)
1150 ( ?- . calc-same-interface
) ))
1152 "The key map for starting Calc.")
1155 ;;;; (Autoloads here)
1156 (load "calc-loaddefs.el" nil t
)
1158 ;;;###autoload (define-key ctl-x-map "*" 'calc-dispatch)
1161 (defun calc-dispatch (&optional arg
)
1162 "Invoke the GNU Emacs Calculator. See `calc-dispatch-help' for details."
1164 ; (sit-for echo-keystrokes)
1165 (condition-case err
; look for other keys bound to calc-dispatch
1166 (let ((keys (this-command-keys)))
1167 (unless (or (not (stringp keys
))
1168 (string-match "\\`\C-u\\|\\`\e[-0-9#]\\|`[\M--\M-0-\M-9]" keys
)
1169 (eq (lookup-key calc-dispatch-map keys
) 'calc-same-interface
))
1170 (when (and (string-match "\\`[\C-@-\C-_]" keys
)
1172 (lookup-key calc-dispatch-map
(substring keys
0 1))))
1173 (define-key calc-dispatch-map
(substring keys
0 1) nil
))
1174 (define-key calc-dispatch-map keys
'calc-same-interface
)))
1176 (calc-do-dispatch arg
))
1178 (defvar calc-dispatch-help nil
)
1179 (defun calc-do-dispatch (arg)
1180 "Start the Calculator."
1181 (let ((key (calc-read-key-sequence
1182 (if calc-dispatch-help
1183 "Calc options: Calc, Keypad, Quick, Embed; eXit; Info, Tutorial; Grab; ?=more"
1184 (format "%s (Type ? for a list of Calc options)"
1185 (key-description (this-command-keys))))
1186 calc-dispatch-map
)))
1187 (setq key
(lookup-key calc-dispatch-map key
))
1191 (or (commandp key
) (require 'calc-ext
))
1192 (call-interactively key
))
1195 (defun calc-read-key-sequence (prompt map
)
1196 "Read keys, with prompt PROMPT and keymap MAP."
1197 (let ((prompt2 (format "%s " (key-description (this-command-keys))))
1198 (glob (current-global-map))
1199 (loc (current-local-map)))
1200 (or (input-pending-p) (message "%s" prompt
))
1201 (let ((key (calc-read-key t
)))
1202 (calc-unread-command (cdr key
))
1205 (use-global-map map
)
1207 (read-key-sequence nil
))
1208 (use-global-map glob
)
1209 (use-local-map loc
)))))
1211 (defvar calc-alg-map
) ; Defined in calc-ext.el
1214 (defvar calc-embedded-modes
) ; Defined in calc-embed.el
1215 (defvar calc-override-minor-modes
) ; Defined in calc-embed.el
1216 (defun calc-kill-stack-buffer ()
1217 "Check to see if user wants to kill the Calc stack buffer.
1218 This will look for buffers using the Calc buffer for embedded mode,
1219 and inform the user if there are any.
1220 If the user wants to kill the Calc buffer, this will remove
1221 embedded information from the appropriate buffers and tidy up
1223 (let ((cb (current-buffer))
1227 (cea calc-embedded-active
))
1228 ;; Get a list of all buffers using this buffer for
1231 (when (and (eq cb
(aref (nth 1 (car cea
)) 1))
1232 (buffer-name (car (car cea
))))
1233 (setq info-list
(cons (car cea
) info-list
)))
1234 (setq cea
(cdr cea
)))
1235 ;; Eventually, prompt user with a list of buffers using embedded mode.
1239 (concat "This Calc stack is being used for embedded mode. Kill anyway?")))
1241 (with-current-buffer (car (car info-list
))
1242 (when calc-embedded-info
1243 (setq calc-embedded-info nil
1244 mode-line-buffer-identification
(car calc-embedded-modes
)
1245 truncate-lines
(nth 2 calc-embedded-modes
)
1246 buffer-read-only nil
)
1247 (use-local-map (nth 1 calc-embedded-modes
))
1248 (setq minor-mode-overriding-map-alist
1249 (remq calc-override-minor-modes minor-mode-overriding-map-alist
))
1250 (let ((str mode-line-buffer-identification
))
1251 (setq mode-line-buffer-identification str
))
1252 (set-buffer-modified-p (buffer-modified-p))))
1253 (setq calc-embedded-active
1254 (delete (car info-list
) calc-embedded-active
))
1255 (setq info-list
(cdr info-list
))))
1258 (setq calc-buffer-list
(delete cb calc-buffer-list
))
1259 (with-current-buffer calc-trail-buffer
1260 (if (eq cb calc-main-buffer
)
1261 ;; If there are other Calc stacks, make another one
1262 ;; the calc-main-buffer ...
1263 (if calc-buffer-list
1264 (setq calc-main-buffer
(car calc-buffer-list
))
1265 ;; ... otherwise kill the trail and its windows.
1266 (let ((wl (get-buffer-window-list calc-trail-buffer
)))
1268 (delete-window (car wl
))
1269 (setq wl
(cdr wl
))))
1270 (kill-buffer calc-trail-buffer
)
1271 (setq calc-trail-buffer nil
))))
1275 "Calculator major mode.
1277 This is an RPN calculator featuring arbitrary-precision integer, rational,
1278 floating-point, complex, matrix, and symbolic arithmetic.
1280 RPN calculation: 2 RET 3 + produces 5.
1281 Algebraic style: ' 2+3 RET produces 5.
1283 Basic operators are +, -, *, /, ^, & (reciprocal), % (modulo), n (change-sign).
1285 Press ? repeatedly for more complete help. Press `h i' to read the
1286 Calc manual on-line, `h s' to read the summary, or `h t' for the tutorial.
1288 Notations: 3.14e6 3.14 * 10^6
1289 _23 negative number -23 (or type `23 n')
1290 17:3 the fraction 17/3
1291 5:2:3 the fraction 5 and 2/3
1292 16#12C the integer 12C base 16 = 300 base 10
1293 8#177:100 the fraction 177:100 base 8 = 127:64 base 10
1294 (2, 4) complex number 2 + 4i
1295 (2; 4) polar complex number (r; theta)
1296 [1, 2, 3] vector ([[1, 2], [3, 4]] is a matrix)
1297 [1 .. 4) semi-open interval, 1 <= x < 4
1298 2 +/- 3 (p key) number with mean 2, standard deviation 3
1299 2 mod 3 (M key) number 2 computed modulo 3
1300 <1 jan 91> Date form (enter using ' key)
1307 (lambda (v) (set-default v
(symbol-value v
)))) calc-local-var-list
)
1308 (kill-all-local-variables)
1309 (use-local-map (if (eq calc-algebraic-mode
'total
)
1310 (progn (require 'calc-ext
) calc-alg-map
) calc-mode-map
))
1311 (mapc (function (lambda (v) (make-local-variable v
))) calc-local-var-list
)
1312 (make-local-variable 'overlay-arrow-position
)
1313 (make-local-variable 'overlay-arrow-string
)
1314 (add-hook 'change-major-mode-hook
'font-lock-defontify nil t
)
1315 (add-hook 'kill-buffer-query-functions
1316 'calc-kill-stack-buffer
1318 (setq truncate-lines t
)
1319 (setq buffer-read-only t
)
1320 (setq major-mode
'calc-mode
)
1321 (setq mode-name
"Calculator")
1322 (setq calc-stack-top
(length (or (memq (assq 'top-of-stack calc-stack
)
1324 (setq calc-stack
(list (list 'top-of-stack
1326 (setq calc-stack-top
(- (length calc-stack
) calc-stack-top -
1))
1327 (or calc-loaded-settings-file
1328 (null calc-settings-file
)
1329 (equal calc-settings-file user-init-file
)
1331 (setq calc-loaded-settings-file t
)
1332 (load (file-name-sans-extension calc-settings-file
) t
))) ; t = missing-ok
1333 (let ((p command-line-args
))
1335 (and (equal (car p
) "-f")
1336 (string-match "calc" (nth 1 p
))
1337 (string-match "full" (nth 1 p
))
1338 (setq calc-standalone-flag t
))
1340 (require 'calc-menu
)
1341 (run-mode-hooks 'calc-mode-hook
)
1343 (calc-set-mode-line)
1344 (calc-check-defines)
1345 (if calc-buffer-list
(setq calc-stack
(copy-sequence calc-stack
)))
1346 (add-to-list 'calc-buffer-list
(current-buffer) t
))
1348 (defvar calc-check-defines
'calc-check-defines
) ; suitable for run-hooks
1349 (defun calc-check-defines ()
1350 (if (symbol-plist 'calc-define
)
1351 (let ((plist (copy-sequence (symbol-plist 'calc-define
))))
1352 (while (and plist
(null (nth 1 plist
)))
1353 (setq plist
(cdr (cdr plist
))))
1357 (require 'calc-macs
)
1358 (set-buffer "*Calculator*")
1360 (put 'calc-define
(car plist
) nil
)
1361 (eval (nth 1 plist
))
1362 (setq plist
(cdr (cdr plist
))))
1363 ;; See if this has added any more calc-define properties.
1364 (calc-check-defines))
1365 (setplist 'calc-define nil
)))))
1367 (defun calc-trail-mode (&optional buf
)
1369 This mode is used by the *Calc Trail* buffer, which records all results
1370 obtained by the GNU Emacs Calculator.
1372 Calculator commands beginning with the `t' key are used to manipulate
1375 This buffer uses the same key map as the *Calculator* buffer; calculator
1376 commands given here will actually operate on the *Calculator* stack."
1379 (use-local-map calc-mode-map
)
1380 (setq major-mode
'calc-trail-mode
)
1381 (setq mode-name
"Calc Trail")
1382 (setq truncate-lines t
)
1383 (setq buffer-read-only t
)
1384 (make-local-variable 'overlay-arrow-position
)
1385 (make-local-variable 'overlay-arrow-string
)
1387 (set (make-local-variable 'calc-main-buffer
) buf
))
1388 (when (= (buffer-size) 0)
1389 (let ((buffer-read-only nil
))
1390 (insert (propertize (concat "Emacs Calculator Trail\n")
1391 'font-lock-face
'italic
))))
1392 (run-mode-hooks 'calc-trail-mode-hook
))
1394 (defun calc-create-buffer ()
1395 "Create and initialize a buffer for the Calculator."
1396 (set-buffer (get-buffer-create "*Calculator*"))
1397 (or (eq major-mode
'calc-mode
)
1399 (setq max-lisp-eval-depth
(max max-lisp-eval-depth
1000))
1400 (when calc-always-load-extensions
1401 (require 'calc-ext
))
1404 (calc-set-language calc-language calc-language-option t
)))
1407 (defun calc (&optional arg full-display interactive
)
1408 "The Emacs Calculator. Full documentation is listed under \"calc-mode\"."
1409 (interactive "P\ni\np")
1413 (if (= (prefix-numeric-value arg
) -
1)
1414 (calc-grab-region (region-beginning) (region-end) nil
)
1415 (when (= (prefix-numeric-value arg
) -
2)
1417 (when (get-buffer-window "*Calc Keypad*")
1419 (set-buffer (window-buffer (selected-window))))
1420 (if (eq major-mode
'calc-mode
)
1422 (let ((oldbuf (current-buffer)))
1423 (calc-create-buffer)
1424 (setq calc-was-keypad-mode nil
)
1425 (if (or (eq full-display t
)
1426 (and (null full-display
) calc-full-mode
))
1427 (switch-to-buffer (current-buffer) t
)
1428 (if (get-buffer-window (current-buffer))
1429 (select-window (get-buffer-window (current-buffer)))
1430 (if calc-window-hook
1431 (run-hooks 'calc-window-hook
)
1432 (let ((w (get-largest-window)))
1433 (if (and pop-up-windows
1434 (> (window-height w
)
1435 (+ window-min-height calc-window-height
2)))
1437 (setq w
(split-window w
1438 (- (window-height w
)
1439 calc-window-height
2)
1441 (set-window-buffer w
(current-buffer))
1443 (pop-to-buffer (current-buffer)))))))
1445 (set-buffer (calc-trail-buffer))
1446 (and calc-display-trail
1447 (= (window-width) (frame-width))
1448 (calc-trail-display 1 t
)))
1449 (message "Welcome to the GNU Emacs Calculator! Press `?' or `h' for help, `q' to quit")
1450 (run-hooks 'calc-start-hook
)
1451 (and (windowp full-display
)
1452 (window-point full-display
)
1453 (select-window full-display
))
1454 (calc-check-defines)
1455 (when (and calc-said-hello interactive
)
1458 (setq calc-said-hello t
)))))
1461 (defun full-calc (&optional interactive
)
1462 "Invoke the Calculator and give it a full-sized window."
1464 (calc nil t interactive
))
1466 (defun calc-same-interface (arg)
1467 "Invoke the Calculator using the most recent interface (calc or calc-keypad)."
1469 (if (and (equal (buffer-name) "*Gnuplot Trail*")
1470 (> (recursion-depth) 0))
1471 (exit-recursive-edit)
1472 (if (eq major-mode
'calc-edit-mode
)
1473 (calc-edit-finish arg
)
1474 (if calc-was-keypad-mode
1476 (calc arg calc-full-mode t
)))))
1478 (defun calc-quit (&optional non-fatal interactive
)
1479 "Quit the Calculator in an appropriate manner."
1480 (interactive "i\np")
1481 (and calc-standalone-flag
(not non-fatal
)
1482 (save-buffers-kill-emacs nil
))
1483 (if (and (equal (buffer-name) "*Gnuplot Trail*")
1484 (> (recursion-depth) 0))
1485 (exit-recursive-edit))
1486 (if (eq major-mode
'calc-edit-mode
)
1488 (if (and interactive
1490 (eq (current-buffer) (aref calc-embedded-info
0)))
1492 (unless (eq major-mode
'calc-mode
)
1493 (calc-create-buffer))
1494 (run-hooks 'calc-end-hook
)
1495 (if (integerp calc-undo-length
)
1497 ((= calc-undo-length
0)
1498 (setq calc-undo-list nil calc-redo-list nil
))
1499 ((> calc-undo-length
0)
1500 (let ((tail (nthcdr (1- calc-undo-length
) calc-undo-list
)))
1501 (if tail
(setcdr tail nil
)))
1502 (setq calc-redo-list nil
))))
1503 (mapc (function (lambda (v) (set-default v
(symbol-value v
))))
1504 calc-local-var-list
)
1505 (let ((buf (current-buffer))
1506 (win (get-buffer-window (current-buffer)))
1507 (kbuf (get-buffer "*Calc Keypad*")))
1508 (delete-windows-on (calc-trail-buffer))
1509 ;; The next few lines will set `calc-window-height' so that the
1510 ;; next time Calc is called, the window will be the same size
1511 ;; as the current window.
1513 (< (window-height win
) (1- (frame-height)))
1514 (window-full-width-p win
) ; avoid calc-keypad
1515 (not (get-buffer-window "*Calc Keypad*")))
1516 (setq calc-window-height
(- (window-height win
) 2)))
1518 (delete-windows-on buf
)
1519 (and kbuf
(delete-windows-on kbuf
)))
1521 (bury-buffer calc-trail-buffer
)
1522 (and kbuf
(bury-buffer kbuf
))))))
1525 (defun quick-calc ()
1526 "Do a quick calculation in the minibuffer without invoking full Calculator."
1528 (calc-do-quick-calc))
1531 (defun calc-eval (str &optional separator
&rest args
)
1532 "Do a quick calculation and return the result as a string.
1533 Return value will either be the formatted result in string form,
1534 or a list containing a character position and an error message in string form."
1535 (calc-do-calc-eval str separator args
))
1538 (defun calc-keypad (&optional interactive
)
1539 "Invoke the Calculator in \"visual keypad\" mode.
1540 This is most useful in the X window system.
1541 In this mode, click on the Calc \"buttons\" using the left mouse button.
1542 Or, position the cursor manually and do M-x calc-keypad-press."
1545 (calc-do-keypad calc-full-mode interactive
))
1548 (defun full-calc-keypad (&optional interactive
)
1549 "Invoke the Calculator in full-screen \"visual keypad\" mode.
1550 See calc-keypad for details."
1553 (calc-do-keypad t interactive
))
1556 (defvar calc-aborted-prefix nil
)
1557 (defvar calc-start-time nil
)
1558 (defvar calc-command-flags
)
1559 (defvar calc-final-point-line
)
1560 (defvar calc-final-point-column
)
1561 ;;; Note that modifications to this function may break calc-pass-errors.
1562 (defun calc-do (do-body &optional do-slow
)
1563 (calc-check-defines)
1564 (let* ((calc-command-flags nil
)
1565 (calc-start-time (and calc-timing
(not calc-start-time
)
1567 (current-time-string)))
1568 (gc-cons-threshold (max gc-cons-threshold
1569 (if calc-timing
2000000 100000)))
1570 calc-final-point-line calc-final-point-column
)
1571 (setq calc-aborted-prefix
"")
1575 (if calc-embedded-info
1576 (calc-embedded-select-buffer)
1577 (calc-select-buffer))
1578 (and (eq calc-algebraic-mode
'total
)
1580 (use-local-map calc-alg-map
))
1581 (when (and do-slow calc-display-working-message
)
1582 (message "Working...")
1583 (calc-set-command-flag 'clear-message
))
1585 (setq calc-aborted-prefix nil
)
1586 (when (memq 'renum-stack calc-command-flags
)
1587 (calc-renumber-stack))
1588 (when (memq 'clear-message calc-command-flags
)
1591 (if (and (eq (car err
) 'error
)
1592 (stringp (nth 1 err
))
1593 (string-match "max-specpdl-size\\|max-lisp-eval-depth"
1595 (error "Computation got stuck or ran too long. Type `M' to increase the limit")
1596 (setq calc-aborted-prefix nil
)
1597 (signal (car err
) (cdr err
)))))
1598 (when calc-aborted-prefix
1599 (calc-record "<Aborted>" calc-aborted-prefix
))
1600 (and calc-start-time
1601 (let* ((calc-internal-prec 12)
1602 (calc-date-format nil
)
1603 (end-time (current-time-string))
1604 (time (if (equal calc-start-time end-time
)
1607 (calcFunc-unixtime (math-parse-date end-time
) 0)
1608 (calcFunc-unixtime (math-parse-date calc-start-time
)
1610 (if (math-lessp 1 time
)
1611 (calc-record time
"(t)"))))
1612 (or (memq 'no-align calc-command-flags
)
1613 (eq major-mode
'calc-trail-mode
)
1614 (calc-align-stack-window))
1615 (and (memq 'position-point calc-command-flags
)
1616 (if (eq major-mode
'calc-mode
)
1618 (goto-char (point-min))
1619 (forward-line (1- calc-final-point-line
))
1620 (move-to-column calc-final-point-column
))
1621 (save-current-buffer
1622 (calc-select-buffer)
1623 (goto-char (point-min))
1624 (forward-line (1- calc-final-point-line
))
1625 (move-to-column calc-final-point-column
))))
1626 (unless (memq 'keep-flags calc-command-flags
)
1628 (calc-select-buffer)
1629 (setq calc-inverse-flag nil
1630 calc-hyperbolic-flag nil
1631 calc-keep-args-flag nil
)))
1632 (when (memq 'do-edit calc-command-flags
)
1633 (switch-to-buffer (get-buffer-create "*Calc Edit*")))
1634 (calc-set-mode-line)
1635 (when calc-embedded-info
1636 (calc-embedded-finish-command))))
1637 (identity nil
)) ; allow a GC after timing is done
1640 (defun calc-set-command-flag (f)
1641 (unless (memq f calc-command-flags
)
1642 (setq calc-command-flags
(cons f calc-command-flags
))))
1644 (defun calc-select-buffer ()
1645 (or (eq major-mode
'calc-mode
)
1646 (if calc-main-buffer
1647 (set-buffer calc-main-buffer
)
1648 (let ((buf (get-buffer "*Calculator*")))
1651 (error "Calculator buffer not available"))))))
1653 (defun calc-cursor-stack-index (&optional index
)
1654 (goto-char (point-max))
1655 (forward-line (- (calc-substack-height (or index
1)))))
1657 (defun calc-stack-size ()
1658 (- (length calc-stack
) calc-stack-top
))
1660 (defun calc-substack-height (n)
1663 (setq n
(+ n calc-stack-top
))
1664 (while (and (> n
0) stack
)
1665 (setq sum
(+ sum
(nth 1 (car stack
)))
1670 (defun calc-set-mode-line ()
1672 (calc-select-buffer)
1673 (let* ((fmt (car calc-float-format
))
1674 (figs (nth 1 calc-float-format
))
1676 (format "Calc%s%s: %d %s %-14s"
1677 (if (and calc-embedded-info
1678 (eq (aref calc-embedded-info
1) (current-buffer)))
1680 (if (and (> (length (buffer-name)) 12)
1681 (equal (substring (buffer-name) 0 12)
1683 (substring (buffer-name) 12)
1686 (capitalize (symbol-name calc-angle-mode
))
1689 ;; Input-related modes
1690 (if (eq calc-algebraic-mode
'total
) "Alg* "
1691 (if calc-algebraic-mode
"Alg "
1692 (if calc-incomplete-algebraic-mode
"Alg[( " "")))
1694 ;; Computational modes
1695 (if calc-symbolic-mode
"Symb " "")
1696 (cond ((eq calc-matrix-mode
'matrix
) "Matrix ")
1697 ((integerp calc-matrix-mode
)
1698 (format "Matrix%d " calc-matrix-mode
))
1699 ((eq calc-matrix-mode
'sqmatrix
) "SqMatrix ")
1700 ((eq calc-matrix-mode
'scalar
) "Scalar ")
1702 (if (eq calc-complex-mode
'polar
) "Polar " "")
1703 (if calc-prefer-frac
"Frac " "")
1704 (cond ((null calc-infinite-mode
) "")
1705 ((eq calc-infinite-mode
1) "+Inf ")
1707 (cond ((eq calc-simplify-mode
'none
) "NoSimp ")
1708 ((eq calc-simplify-mode
'num
) "NumSimp ")
1709 ((eq calc-simplify-mode
'binary
)
1710 (format "BinSimp%d " calc-word-size
))
1711 ((eq calc-simplify-mode
'alg
) "AlgSimp ")
1712 ((eq calc-simplify-mode
'ext
) "ExtSimp ")
1713 ((eq calc-simplify-mode
'units
) "UnitSimp ")
1717 (cond ((= calc-number-radix
10) "")
1718 ((= calc-number-radix
2) "Bin ")
1719 ((= calc-number-radix
8) "Oct ")
1720 ((= calc-number-radix
16) "Hex ")
1721 (t (format "Radix%d " calc-number-radix
)))
1722 (if calc-leading-zeros
"Zero " "")
1723 (cond ((null calc-language
) "")
1724 ((get calc-language
'math-lang-name
)
1725 (concat (get calc-language
'math-lang-name
) " "))
1727 (capitalize (symbol-name calc-language
))
1729 (cond ((eq fmt
'float
)
1730 (if (zerop figs
) "" (format "Norm%d " figs
)))
1731 ((eq fmt
'fix
) (format "Fix%d " figs
))
1733 (if (zerop figs
) "Sci " (format "Sci%d " figs
)))
1735 (if (zerop figs
) "Eng " (format "Eng%d " figs
))))
1736 (cond ((not calc-display-just
)
1737 (if calc-display-origin
1738 (format "Left%d " calc-display-origin
) ""))
1739 ((eq calc-display-just
'right
)
1740 (if calc-display-origin
1741 (format "Right%d " calc-display-origin
)
1744 (if calc-display-origin
1745 (format "Center%d " calc-display-origin
)
1747 (cond ((integerp calc-line-breaking
)
1748 (format "Wid%d " calc-line-breaking
))
1749 (calc-line-breaking "")
1752 ;; Miscellaneous other modes/indicators
1753 (if calc-assoc-selections
"" "Break ")
1754 (cond ((eq calc-mode-save-mode
'save
) "Save ")
1755 ((not calc-embedded-info
) "")
1756 ((eq calc-mode-save-mode
'local
) "Local ")
1757 ((eq calc-mode-save-mode
'edit
) "LocEdit ")
1758 ((eq calc-mode-save-mode
'perm
) "LocPerm ")
1759 ((eq calc-mode-save-mode
'global
) "Global ")
1761 (if calc-auto-recompute
"" "Manual ")
1762 (if (and (fboundp 'calc-gnuplot-alive
)
1763 (calc-gnuplot-alive)) "Graph " "")
1764 (if (and calc-embedded-info
1765 (> (calc-stack-size) 0)
1766 (calc-top 1 'sel
)) "Sel " "")
1767 (if calc-display-dirty
"Dirty " "")
1768 (if calc-inverse-flag
"Inv " "")
1769 (if calc-hyperbolic-flag
"Hyp " "")
1770 (if calc-keep-args-flag
"Keep " "")
1771 (if (/= calc-stack-top
1) "Narrow " "")
1772 (apply 'concat calc-other-modes
)))))
1773 (if (equal new-mode-string mode-line-buffer-identification
)
1775 (setq mode-line-buffer-identification new-mode-string
)
1776 (set-buffer-modified-p (buffer-modified-p))
1777 (and calc-embedded-info
(calc-embedded-mode-line-change))))))
1779 (defun calc-align-stack-window ()
1780 (if (eq major-mode
'calc-mode
)
1782 (let ((win (get-buffer-window (current-buffer))))
1785 (calc-cursor-stack-index 0)
1786 (vertical-motion (- 2 (window-height win
)))
1787 (set-window-start win
(point)))))
1788 (calc-cursor-stack-index 0)
1789 (if (looking-at " *\\.$")
1790 (goto-char (1- (match-end 0)))))
1792 (calc-select-buffer)
1793 (calc-align-stack-window))))
1795 (defun calc-check-stack (n)
1796 (if (> n
(calc-stack-size))
1797 (error "Too few elements on stack"))
1799 (error "Invalid argument")))
1801 (defun calc-push-list (vals &optional m sels
)
1803 (if calc-executing-macro
1804 (calc-push-list-in-macro vals m sels
)
1806 (calc-select-buffer)
1807 (let* ((val (car vals
))
1808 (entry (list val
1 (car sels
)))
1809 (mm (+ (or m
1) calc-stack-top
)))
1810 (calc-cursor-stack-index (1- (or m
1)))
1812 (setcdr (nthcdr (- mm
2) calc-stack
)
1813 (cons entry
(nthcdr (1- mm
) calc-stack
)))
1814 (setq calc-stack
(cons entry calc-stack
)))
1815 (let ((buffer-read-only nil
))
1816 (insert (math-format-stack-value entry
) "\n"))
1817 (calc-record-undo (list 'push mm
))
1818 (calc-set-command-flag 'renum-stack
))))
1819 (setq vals
(cdr vals
)
1822 (defun calc-pop-push-list (n vals
&optional m sels
)
1823 (if (and calc-any-selections
(null sels
))
1824 (calc-replace-selections n vals m
)
1825 (calc-pop-stack n m sels
)
1826 (calc-push-list vals m sels
)))
1828 (defun calc-pop-push-record-list (n prefix vals
&optional m sels
)
1829 (or (and (consp vals
)
1830 (or (integerp (car vals
))
1831 (consp (car vals
))))
1832 (and vals
(setq vals
(list vals
)
1833 sels
(and sels
(list sels
)))))
1834 (calc-check-stack (+ n
(or m
1) -
1))
1837 (calc-record-list vals prefix
)
1838 (calc-record (car vals
) prefix
)))
1839 (calc-pop-push-list n vals m sels
))
1841 (defun calc-enter-result (n prefix vals
&optional m
)
1842 (setq calc-aborted-prefix prefix
)
1843 (if (and (consp vals
)
1844 (or (integerp (car vals
))
1845 (consp (car vals
))))
1846 (setq vals
(mapcar 'calc-normalize vals
))
1847 (setq vals
(calc-normalize vals
)))
1848 (or (and (consp vals
)
1849 (or (integerp (car vals
))
1850 (consp (car vals
))))
1851 (setq vals
(list vals
)))
1852 (if (equal vals
'((nil)))
1854 (calc-pop-push-record-list n prefix vals m
)
1857 (defun calc-normalize (val)
1858 (if (memq calc-simplify-mode
'(nil none num
))
1859 (math-normalize val
)
1861 (calc-normalize-fancy val
)))
1863 (defun calc-handle-whys ()
1865 (calc-do-handle-whys)))
1868 (defun calc-pop-stack (&optional n m sel-ok
) ; pop N objs at level M of stack.
1871 (or calc-keep-args-flag
1872 (let ((mm (+ m calc-stack-top
)))
1873 (if (and calc-any-selections
(not sel-ok
)
1874 (calc-top-selected n m
))
1876 (if calc-executing-macro
1877 (calc-pop-stack-in-macro n mm
)
1878 (calc-record-undo (list 'pop mm
(calc-top-list n m
'full
)))
1880 (calc-select-buffer)
1881 (let ((buffer-read-only nil
))
1884 (calc-cursor-stack-index (1- m
))
1885 (let ((bot (point)))
1886 (calc-cursor-stack-index (+ n m -
1))
1887 (delete-region (point) bot
))
1888 (setcdr (nthcdr (- mm
2) calc-stack
)
1889 (nthcdr (+ n mm -
1) calc-stack
)))
1890 (calc-cursor-stack-index n
)
1891 (setq calc-stack
(nthcdr n calc-stack
))
1892 (delete-region (point) (point-max))))
1893 (calc-set-command-flag 'renum-stack
))))))
1896 (defun calc-get-stack-element (x)
1897 (cond ((eq sel-mode
'entry
)
1901 ((or (null (nth 2 x
))
1903 (not calc-use-selections
))
1909 ;; Get the Nth element of the stack (N=1 is the top element).
1910 (defun calc-top (&optional n sel-mode
)
1912 (calc-check-stack n
)
1913 (calc-get-stack-element (nth (+ n calc-stack-top -
1) calc-stack
)))
1915 (defun calc-top-n (&optional n sel-mode
) ; in case precision has changed
1916 (math-check-complete (calc-normalize (calc-top n sel-mode
))))
1918 (defun calc-top-list (&optional n m sel-mode
)
1921 (calc-check-stack (+ n m -
1))
1923 (let ((top (copy-sequence (nthcdr (+ m calc-stack-top -
1)
1925 (setcdr (nthcdr (1- n
) top
) nil
)
1926 (nreverse (mapcar 'calc-get-stack-element top
)))))
1928 (defun calc-top-list-n (&optional n m sel-mode
)
1929 (mapcar 'math-check-complete
1930 (mapcar 'calc-normalize
(calc-top-list n m sel-mode
))))
1933 (defun calc-renumber-stack ()
1934 (if calc-line-numbering
1936 (calc-cursor-stack-index 0)
1938 (buffer-read-only nil
)
1939 (stack (nthcdr calc-stack-top calc-stack
)))
1940 (if (re-search-forward "^[0-9]+[:*]" nil t
)
1943 (while (re-search-forward "^[0-9]+[:*]" nil t
)
1944 (let ((buffer-read-only nil
))
1948 (calc-cursor-stack-index 0)))
1949 (while (re-search-backward "^[0-9]+[:*]" nil t
)
1952 (insert (format "%03d%s" (% lnum
1000)
1953 (if (and (nth 2 (car stack
))
1954 calc-use-selections
) "*" ":")))
1955 (let ((prefix (int-to-string lnum
)))
1956 (insert prefix
(if (and (nth 2 (car stack
))
1957 calc-use-selections
) "*" ":")
1958 (make-string (- 3 (length prefix
)) 32))))
1960 (setq lnum
(1+ lnum
)
1961 stack
(cdr stack
))))))
1962 (and calc-embedded-info
(calc-embedded-stack-change)))
1964 (defvar calc-any-evaltos nil
)
1965 (defun calc-refresh (&optional align
)
1967 (and (eq major-mode
'calc-mode
)
1968 (not calc-executing-macro
)
1969 (let* ((buffer-read-only nil
)
1970 (save-point (point))
1971 (save-mark (condition-case err
(mark) (error nil
)))
1972 (save-aligned (looking-at "\\.$"))
1974 (calc-any-evaltos nil
))
1975 (setq calc-any-selections nil
)
1977 (when calc-show-banner
1978 (insert (propertize "--- Emacs Calculator Mode ---\n"
1979 'font-lock-face
'italic
)))
1981 (goto-char (point-min))
1982 (when calc-show-banner
1984 (insert (math-format-stack-value (car thing
)) "\n")
1985 (setq thing
(cdr thing
)))
1986 (calc-renumber-stack)
1987 (if calc-display-dirty
1988 (calc-wrapper (setq calc-display-dirty nil
)))
1989 (and calc-any-evaltos calc-auto-recompute
1990 (calc-wrapper (calc-refresh-evaltos)))
1991 (if (or align save-aligned
)
1992 (calc-align-stack-window)
1993 (goto-char save-point
))
1994 (if save-mark
(set-mark save-mark
))))
1995 (and calc-embedded-info
(not (eq major-mode
'calc-mode
))
1997 (set-buffer (aref calc-embedded-info
1))
1998 (calc-refresh align
)))
1999 (setq calc-refresh-count
(1+ calc-refresh-count
)))
2001 ;;;; The Calc Trail buffer.
2003 (defun calc-check-trail-aligned ()
2005 (let ((win (get-buffer-window (current-buffer))))
2007 (pos-visible-in-window-p (1- (point-max)) win
)))))
2009 (defun calc-trail-buffer ()
2010 (and (or (null calc-trail-buffer
)
2011 (null (buffer-name calc-trail-buffer
)))
2013 (setq calc-trail-buffer
(get-buffer-create "*Calc Trail*"))
2014 (let ((buf (or (and (not (eq major-mode
'calc-mode
))
2015 (get-buffer "*Calculator*"))
2017 (set-buffer calc-trail-buffer
)
2018 (or (eq major-mode
'calc-trail-mode
)
2019 (calc-trail-mode buf
)))))
2020 (or (and calc-trail-pointer
2021 (eq (marker-buffer calc-trail-pointer
) calc-trail-buffer
))
2023 (set-buffer calc-trail-buffer
)
2024 (goto-char (point-min))
2026 (setq calc-trail-pointer
(point-marker))))
2029 (defun calc-record (val &optional prefix
)
2030 (setq calc-aborted-prefix nil
)
2031 (or calc-executing-macro
2032 (let* ((mainbuf (current-buffer))
2033 (buf (calc-trail-buffer))
2034 (calc-display-raw nil
)
2035 (calc-can-abbrev-vectors t
)
2039 (math-showing-full-precision
2040 (math-format-flat-expr val
0)))
2044 (let ((aligned (calc-check-trail-aligned))
2045 (buffer-read-only nil
))
2046 (goto-char (point-max))
2047 (cond ((null prefix
) (insert " "))
2048 ((and (> (length prefix
) 4)
2049 (string-match " " prefix
4))
2050 (insert (substring prefix
0 4) " "))
2051 (t (insert (format "%4s " prefix
))))
2053 (let ((win (get-buffer-window buf
)))
2054 (if (and aligned win
(not (memq 'hold-trail calc-command-flags
)))
2056 (goto-char (1- (point-max))))))))
2060 (defun calc-trail-display (flag &optional no-refresh interactive
)
2061 (interactive "P\ni\np")
2062 (let ((win (get-buffer-window (calc-trail-buffer))))
2063 (if (setq calc-display-trail
2064 (not (if flag
(memq flag
'(nil 0)) win
)))
2067 (if calc-trail-window-hook
2068 (run-hooks 'calc-trail-window-hook
)
2069 (let ((w (split-window nil
(/ (* (window-width) 2) 3) t
)))
2070 (set-window-buffer w calc-trail-buffer
)))
2072 (setq overlay-arrow-string calc-trail-overlay
2073 overlay-arrow-position calc-trail-pointer
)
2085 (calc-refresh))))))))
2088 (defun calc-trail-here ()
2090 (if (eq major-mode
'calc-trail-mode
)
2097 (if (or (bobp) (eobp))
2098 (setq overlay-arrow-position nil
) ; trail is empty
2099 (set-marker calc-trail-pointer
(point) (current-buffer))
2100 (setq calc-trail-overlay
(concat (buffer-substring (point)
2103 overlay-arrow-string calc-trail-overlay
2104 overlay-arrow-position calc-trail-pointer
)
2106 (let ((win (get-buffer-window (current-buffer))))
2109 (forward-line (/ (window-height win
) 2))
2110 (forward-line (- 1 (window-height win
)))
2111 (set-window-start win
(point))
2112 (set-window-point win
(+ calc-trail-pointer
4))
2113 (set-buffer calc-main-buffer
)
2114 (setq overlay-arrow-string calc-trail-overlay
2115 overlay-arrow-position calc-trail-pointer
))))))
2116 (error "Not in Calc Trail buffer")))
2123 (defun calc-record-undo (rec)
2124 (or calc-executing-macro
2125 (if (memq 'undo calc-command-flags
)
2126 (setq calc-undo-list
(cons (cons rec
(car calc-undo-list
))
2127 (cdr calc-undo-list
)))
2128 (setq calc-undo-list
(cons (list rec
) calc-undo-list
)
2130 (calc-set-command-flag 'undo
))))
2135 ;;; Arithmetic commands.
2137 (defun calc-binary-op (name func arg
&optional ident unary func2
)
2138 (setq calc-aborted-prefix name
)
2140 (calc-enter-result 2 name
(cons (or func2 func
)
2141 (mapcar 'math-check-complete
2142 (calc-top-list 2))))
2144 (calc-binary-op-fancy name func arg ident unary
)))
2146 (defun calc-unary-op (name func arg
&optional func2
)
2147 (setq calc-aborted-prefix name
)
2149 (calc-enter-result 1 name
(list (or func2 func
)
2150 (math-check-complete (calc-top 1))))
2152 (calc-unary-op-fancy name func arg
)))
2155 (defun calc-plus (arg)
2158 (calc-binary-op "+" 'calcFunc-add arg
0 nil
'+)))
2160 (defun calc-minus (arg)
2163 (calc-binary-op "-" 'calcFunc-sub arg
0 'neg
'-
)))
2165 (defun calc-times (arg)
2168 (calc-binary-op "*" 'calcFunc-mul arg
1 nil
'*)))
2170 (defun calc-divide (arg)
2173 (calc-binary-op "/" 'calcFunc-div arg
0 'calcFunc-inv
'/)))
2175 (defun calc-left-divide (arg)
2178 (calc-binary-op "ldiv" 'calcFunc-ldiv arg
0 nil nil
)))
2180 (defun calc-change-sign (arg)
2183 (calc-unary-op "chs" 'neg arg
)))
2187 ;;; Stack management commands.
2189 (defun calc-enter (n)
2193 (calc-push-list (calc-top-list 1 (- n
))))
2195 (calc-push-list (calc-top-list (calc-stack-size))))
2197 (calc-push-list (calc-top-list n
))))))
2203 (let* ((nn (prefix-numeric-value n
))
2204 (top (and (null n
) (calc-top 1))))
2205 (cond ((and (null n
)
2206 (eq (car-safe top
) 'incomplete
)
2207 (> (length top
) (if (eq (nth 1 top
) 'intv
) 3 2)))
2208 (calc-pop-push-list 1 (let ((tt (copy-sequence top
)))
2209 (setcdr (nthcdr (- (length tt
) 2) tt
) nil
)
2212 (if (and calc-any-selections
2213 (calc-top-selected 1 (- nn
)))
2214 (calc-delete-selection (- nn
))
2215 (calc-pop-stack 1 (- nn
) t
)))
2217 (calc-pop-stack (calc-stack-size) 1 t
))
2219 (if (and calc-any-selections
2221 (calc-top-selected 1 1))
2222 (calc-delete-selection 1)
2223 (calc-pop-stack nn
)))))))
2228 ;;;; Reading a number using the minibuffer.
2229 (defvar calc-buffer
)
2230 (defvar calc-prev-char
)
2231 (defvar calc-prev-prev-char
)
2232 (defvar calc-digit-value
)
2233 (defun calcDigit-start ()
2236 (if (or calc-algebraic-mode
2237 (and (> calc-number-radix
14) (eq last-command-event ?e
)))
2238 (calc-alg-digit-entry)
2239 (calc-unread-command)
2240 (setq calc-aborted-prefix nil
)
2241 (let* ((calc-digit-value nil
)
2242 (calc-prev-char nil
)
2243 (calc-prev-prev-char nil
)
2244 (calc-buffer (current-buffer))
2245 (buf (if (featurep 'xemacs
)
2247 (catch 'execute-kbd-macro
2249 (read-from-minibuffer
2250 "Calc: " "" calc-digit-map
)))
2251 (error "XEmacs requires RET after %s"
2252 "digit entry in kbd macro"))
2253 (let ((old-esc (lookup-key global-map
"\e")))
2256 (define-key global-map
"\e" nil
)
2257 (read-from-minibuffer "Calc: " "" calc-digit-map
))
2258 (define-key global-map
"\e" old-esc
))))))
2259 (or calc-digit-value
(setq calc-digit-value
(math-read-number buf
)))
2260 (if (stringp calc-digit-value
)
2261 (calc-alg-entry calc-digit-value
)
2262 (if calc-digit-value
2263 (calc-push-list (list (calc-record (calc-normalize
2264 calc-digit-value
))))))
2265 (if (eq calc-prev-char
'dots
)
2270 (defsubst calc-minibuffer-size
()
2271 (- (point-max) (minibuffer-prompt-end)))
2273 (defun calcDigit-nondigit ()
2275 ;; Exercise for the reader: Figure out why this is a good precaution!
2276 (or (boundp 'calc-buffer
)
2277 (use-local-map minibuffer-local-map
))
2278 (let ((str (minibuffer-contents)))
2279 (setq calc-digit-value
(save-excursion
2280 (set-buffer calc-buffer
)
2281 (math-read-number str
))))
2282 (if (and (null calc-digit-value
) (> (calc-minibuffer-size) 0))
2285 (calc-temp-minibuffer-message " [Bad format]"))
2286 (or (memq last-command-event
'(32 13))
2287 (progn (setq prefix-arg current-prefix-arg
)
2288 (calc-unread-command (if (and (eq last-command-event
27)
2289 (>= last-input-event
128))
2295 (defun calc-minibuffer-contains (rex)
2297 (goto-char (minibuffer-prompt-end))
2300 (defun calcDigit-key ()
2302 (goto-char (point-max))
2303 (if (or (and (memq last-command-event
'(?
+ ?-
))
2305 (/= (preceding-char) ?e
))
2306 (and (memq last-command-event
'(?m ?s
))
2307 (not (calc-minibuffer-contains "[-+]?[0-9]+\\.?0*[@oh].*"))
2308 (not (calc-minibuffer-contains "[-+]?\\(1[1-9]\\|[2-9][0-9]\\)#.*"))))
2309 (calcDigit-nondigit)
2310 (if (calc-minibuffer-contains "\\([-+]?\\|.* \\)\\'")
2311 (cond ((memq last-command-event
'(?. ?
@)) (insert "0"))
2312 ((and (memq last-command-event
'(?o ?h ?m
))
2313 (not (calc-minibuffer-contains ".*#.*"))) (insert "0"))
2314 ((memq last-command-event
'(?
: ?e
)) (insert "1"))
2315 ((eq last-command-event ?
#)
2316 (insert (int-to-string calc-number-radix
)))))
2317 (if (and (calc-minibuffer-contains "\\([-+]?[0-9]+#\\|[^:]*:\\)\\'")
2318 (eq last-command-event ?
:))
2320 (if (and (calc-minibuffer-contains "[-+]?[0-9]+#\\'")
2321 (eq last-command-event ?.
))
2323 (if (and (calc-minibuffer-contains "[-+]?0*\\([2-9]\\|1[0-4]\\)#\\'")
2324 (eq last-command-event ?e
))
2326 (if (or (and (memq last-command-event
'(?h ?o ?m ?s ?p
))
2327 (calc-minibuffer-contains ".*#.*"))
2328 (and (eq last-command-event ?e
)
2329 (calc-minibuffer-contains "[-+]?\\(1[5-9]\\|[2-9][0-9]\\)#.*"))
2330 (and (eq last-command-event ?n
)
2331 (calc-minibuffer-contains "[-+]?\\(2[4-9]\\|[3-9][0-9]\\)#.*")))
2332 (setq last-command-event
(upcase last-command-event
)))
2334 ((memq last-command-event
'(?_ ?n
))
2335 (goto-char (minibuffer-prompt-end))
2336 (if (and (search-forward " +/- " nil t
)
2337 (not (search-forward "e" nil t
)))
2339 (and (not (calc-minibuffer-contains "[-+]?\\(1[5-9]\\|[2-9][0-9]\\)#.*"))
2340 (search-forward "e" nil t
))
2341 (if (looking-at "+")
2343 (if (looking-at "-")
2346 (goto-char (point-max)))
2347 ((eq last-command-event ?p
)
2348 (if (or (calc-minibuffer-contains ".*\\+/-.*")
2349 (calc-minibuffer-contains ".*mod.*")
2350 (calc-minibuffer-contains ".*#.*")
2351 (calc-minibuffer-contains ".*[-+e:]\\'"))
2353 (if (not (calc-minibuffer-contains ".* \\'"))
2356 ((and (eq last-command-event ?M
)
2357 (not (calc-minibuffer-contains
2358 "[-+]?\\(2[3-9]\\|[3-9][0-9]\\)#.*")))
2359 (if (or (calc-minibuffer-contains ".*\\+/-.*")
2360 (calc-minibuffer-contains ".*mod *[^ ]+")
2361 (calc-minibuffer-contains ".*[-+e:]\\'"))
2363 (if (calc-minibuffer-contains ".*mod \\'")
2364 (if calc-previous-modulo
2365 (insert (math-format-flat-expr calc-previous-modulo
0))
2367 (if (not (calc-minibuffer-contains ".* \\'"))
2371 (insert (char-to-string last-command-event
))
2372 (if (or (and (calc-minibuffer-contains "[-+]?\\(.*\\+/- *\\|.*mod *\\)?\\([0-9][0-9]?\\)#[0-9a-zA-Z]*\\(:[0-9a-zA-Z]*\\(:[0-9a-zA-Z]*\\)?\\|.[0-9a-zA-Z]*\\(e[-+]?[0-9]*\\)?\\)?\\'")
2373 (let ((radix (string-to-number
2375 (match-beginning 2) (match-end 2)))))
2378 (or (memq last-command-event
'(?
# ?
: ?. ?e ?
+ ?-
))
2379 (let ((dig (math-read-radix-digit
2380 (upcase last-command-event
))))
2383 (calc-minibuffer-contains
2384 "[-+]?\\(.*\\+/- *\\|.*mod *\\)?\\([0-9]+\\.?0*[@oh] *\\)?\\([0-9]+\\.?0*['m] *\\)?[0-9]*\\(\\.?[0-9]*\\(e[-+]?[0-3]?[0-9]?[0-9]?[0-9]?[0-9]?[0-9]?[0-9]?\\)?\\|[0-9]:\\([0-9]+:\\)?[0-9]*\\)?[\"s]?\\'"))
2385 (if (and (memq last-command-event
'(?
@ ?o ?h ?
\' ?m
))
2386 (string-match " " calc-hms-format
))
2388 (if (and (eq this-command last-command
)
2389 (eq last-command-event ?.
))
2393 (delete-backward-char 1)
2395 (calc-temp-minibuffer-message " [Bad format]"))))))
2396 (setq calc-prev-prev-char calc-prev-char
2397 calc-prev-char last-command-event
))
2400 (defun calcDigit-backspace ()
2402 (goto-char (point-max))
2403 (cond ((calc-minibuffer-contains ".* \\+/- \\'")
2404 (backward-delete-char 5))
2405 ((calc-minibuffer-contains ".* mod \\'")
2406 (backward-delete-char 5))
2407 ((calc-minibuffer-contains ".* \\'")
2408 (backward-delete-char 2))
2409 ((eq last-command
'calcDigit-start
)
2411 (t (backward-delete-char 1)))
2412 (if (= (calc-minibuffer-size) 0)
2414 (setq last-command-event
13)
2415 (calcDigit-nondigit))))
2420 (defconst math-bignum-digit-length
2421 (truncate (/ (log10 (/ most-positive-fixnum
2)) 2))
2422 "The length of a \"digit\" in Calc bignums.
2423 If a big integer is of the form (bigpos N0 N1 ...), this is the
2424 length of the allowable Emacs integers N0, N1,...
2425 The value of 2*10^(2*MATH-BIGNUM-DIGIT-LENGTH) must be less than the
2426 largest Emacs integer.")
2428 (defconst math-bignum-digit-size
2429 (expt 10 math-bignum-digit-length
)
2430 "An upper bound for the size of the \"digit\"s in Calc bignums.")
2432 (defconst math-small-integer-size
2433 (expt math-bignum-digit-size
2)
2434 "An upper bound for the size of \"small integer\"s in Calc.")
2437 ;;;; Arithmetic routines.
2439 ;;; An object as manipulated by one of these routines may take any of the
2440 ;;; following forms:
2442 ;;; integer An integer. For normalized numbers, this format
2443 ;;; is used only for
2444 ;;; negative math-small-integer-size + 1 to
2445 ;;; math-small-integer-size - 1
2447 ;;; (bigpos N0 N1 N2 ...) A big positive integer,
2448 ;;; N0 + N1*math-bignum-digit-size
2449 ;;; + N2*(math-bignum-digit-size)^2 ...
2450 ;;; (bigneg N0 N1 N2 ...) A big negative integer,
2451 ;;; - N0 - N1*math-bignum-digit-size ...
2452 ;;; Each digit N is in the range
2453 ;;; 0 ... math-bignum-digit-size -1.
2454 ;;; Normalized, always at least three N present,
2455 ;;; and the most significant N is nonzero.
2457 ;;; (frac NUM DEN) A fraction. NUM and DEN are small or big integers.
2458 ;;; Normalized, DEN > 1.
2460 ;;; (float NUM EXP) A floating-point number, NUM * 10^EXP;
2461 ;;; NUM is a small or big integer, EXP is a small int.
2462 ;;; Normalized, NUM is not a multiple of 10, and
2463 ;;; abs(NUM) < 10^calc-internal-prec.
2464 ;;; Normalized zero is stored as (float 0 0).
2466 ;;; (cplx REAL IMAG) A complex number; REAL and IMAG are any of above.
2467 ;;; Normalized, IMAG is nonzero.
2469 ;;; (polar R THETA) Polar complex number. Normalized, R > 0 and THETA
2470 ;;; is neither zero nor 180 degrees (pi radians).
2472 ;;; (vec A B C ...) Vector of objects A, B, C, ... A matrix is a
2473 ;;; vector of vectors.
2475 ;;; (hms H M S) Angle in hours-minutes-seconds form. All three
2476 ;;; components have the same sign; H and M must be
2477 ;;; numerically integers; M and S are expected to
2478 ;;; lie in the range [0,60).
2480 ;;; (date N) A date or date/time object. N is an integer to
2481 ;;; store a date only, or a fraction or float to
2482 ;;; store a date and time.
2484 ;;; (sdev X SIGMA) Error form, X +/- SIGMA. When normalized,
2485 ;;; SIGMA > 0. X is any complex number and SIGMA
2486 ;;; is real numbers; or these may be symbolic
2487 ;;; expressions where SIGMA is assumed real.
2489 ;;; (intv MASK LO HI) Interval form. MASK is 0=(), 1=(], 2=[), or 3=[].
2490 ;;; LO and HI are any real numbers, or symbolic
2491 ;;; expressions which are assumed real, and LO < HI.
2492 ;;; For [LO..HI], if LO = HI normalization produces LO,
2493 ;;; and if LO > HI normalization produces [LO..LO).
2494 ;;; For other intervals, if LO > HI normalization
2495 ;;; sets HI equal to LO.
2497 ;;; (mod N M) Number modulo M. When normalized, 0 <= N < M.
2498 ;;; N and M are real numbers.
2500 ;;; (var V S) Symbolic variable. V is a Lisp symbol which
2501 ;;; represents the variable's visible name. S is
2502 ;;; the symbol which actually stores the variable's
2503 ;;; value: (var pi var-pi).
2505 ;;; In general, combining rational numbers in a calculation always produces
2506 ;;; a rational result, but if either argument is a float, result is a float.
2508 ;;; In the following comments, [x y z] means result is x, args must be y, z,
2509 ;;; respectively, where the code letters are:
2511 ;;; O Normalized object (vector or number)
2512 ;;; V Normalized vector
2513 ;;; N Normalized number of any type
2514 ;;; N Normalized complex number
2515 ;;; R Normalized real number (float or rational)
2516 ;;; F Normalized floating-point number
2517 ;;; T Normalized rational number
2518 ;;; I Normalized integer
2519 ;;; B Normalized big integer
2520 ;;; S Normalized small integer
2521 ;;; D Digit (small integer, 0..999)
2522 ;;; L Normalized bignum digit list (without "bigpos" or "bigneg" symbol)
2523 ;;; or normalized vector element list (without "vec")
2524 ;;; P Predicate (truth value)
2525 ;;; X Any Lisp object
2528 ;;; Lower-case letters signify possibly un-normalized values.
2529 ;;; "L.D" means a cons of an L and a D.
2530 ;;; [N N; n n] means result will be normalized if argument is.
2531 ;;; Also, [Public] marks routines intended to be called from outside.
2532 ;;; [This notation has been neglected in many recent routines.]
2534 (defvar math-eval-rules-cache
)
2535 (defvar math-eval-rules-cache-other
)
2536 ;;; Reduce an object to canonical (normalized) form. [O o; Z Z] [Public]
2538 (defvar math-normalize-a
)
2539 (defun math-normalize (math-normalize-a)
2541 ((not (consp math-normalize-a
))
2542 (if (integerp math-normalize-a
)
2543 (if (or (>= math-normalize-a math-small-integer-size
)
2544 (<= math-normalize-a
(- math-small-integer-size
)))
2545 (math-bignum math-normalize-a
)
2548 ((eq (car math-normalize-a
) 'bigpos
)
2549 (if (eq (nth (1- (length math-normalize-a
)) math-normalize-a
) 0)
2550 (let* ((last (setq math-normalize-a
2551 (copy-sequence math-normalize-a
))) (digs math-normalize-a
))
2552 (while (setq digs
(cdr digs
))
2553 (or (eq (car digs
) 0) (setq last digs
)))
2555 (if (cdr (cdr (cdr math-normalize-a
)))
2558 ((cdr (cdr math-normalize-a
)) (+ (nth 1 math-normalize-a
)
2559 (* (nth 2 math-normalize-a
)
2560 math-bignum-digit-size
)))
2561 ((cdr math-normalize-a
) (nth 1 math-normalize-a
))
2563 ((eq (car math-normalize-a
) 'bigneg
)
2564 (if (eq (nth (1- (length math-normalize-a
)) math-normalize-a
) 0)
2565 (let* ((last (setq math-normalize-a
(copy-sequence math-normalize-a
)))
2566 (digs math-normalize-a
))
2567 (while (setq digs
(cdr digs
))
2568 (or (eq (car digs
) 0) (setq last digs
)))
2570 (if (cdr (cdr (cdr math-normalize-a
)))
2573 ((cdr (cdr math-normalize-a
)) (- (+ (nth 1 math-normalize-a
)
2574 (* (nth 2 math-normalize-a
)
2575 math-bignum-digit-size
))))
2576 ((cdr math-normalize-a
) (- (nth 1 math-normalize-a
)))
2578 ((eq (car math-normalize-a
) 'float
)
2579 (math-make-float (math-normalize (nth 1 math-normalize-a
))
2580 (nth 2 math-normalize-a
)))
2581 ((or (memq (car math-normalize-a
)
2582 '(frac cplx polar hms date mod sdev intv vec var quote
2583 special-const calcFunc-if calcFunc-lambda
2584 calcFunc-quote calcFunc-condition
2586 (integerp (car math-normalize-a
))
2587 (and (consp (car math-normalize-a
))
2588 (not (eq (car (car math-normalize-a
)) 'lambda
))))
2590 (math-normalize-fancy math-normalize-a
))
2592 (or (and calc-simplify-mode
2594 (math-normalize-nonstandard))
2595 (let ((args (mapcar 'math-normalize
(cdr math-normalize-a
))))
2596 (or (condition-case err
2598 (assq (car math-normalize-a
) '( ( + . math-add
)
2605 ( | . math-concat
) ))))
2606 (or (and var-EvalRules
2608 (or (eq var-EvalRules math-eval-rules-cache-tag
)
2611 (math-recompile-eval-rules)))
2612 (and (or math-eval-rules-cache-other
2613 (assq (car math-normalize-a
)
2614 math-eval-rules-cache
))
2615 (math-apply-rewrites
2616 (cons (car math-normalize-a
) args
)
2617 (cdr math-eval-rules-cache
)
2618 nil math-eval-rules-cache
))))
2620 (apply (cdr func
) args
)
2621 (and (or (consp (car math-normalize-a
))
2622 (fboundp (car math-normalize-a
))
2623 (and (not (featurep 'calc-ext
))
2625 (fboundp (car math-normalize-a
))))
2626 (apply (car math-normalize-a
) args
)))))
2627 (wrong-number-of-arguments
2628 (calc-record-why "*Wrong number of arguments"
2629 (cons (car math-normalize-a
) args
))
2631 (wrong-type-argument
2633 (calc-record-why "Wrong type of argument"
2634 (cons (car math-normalize-a
) args
)))
2637 (calc-record-why "*Argument out of range"
2638 (cons (car math-normalize-a
) args
))
2641 (calc-record-why "No exact representation for result"
2642 (cons (car math-normalize-a
) args
))
2645 (calc-record-why "*Floating-point overflow occurred"
2646 (cons (car math-normalize-a
) args
))
2649 (calc-record-why "*Floating-point underflow occurred"
2650 (cons (car math-normalize-a
) args
))
2653 (if (eq (nth 1 err
) 'var-EvalRules
)
2655 (setq var-EvalRules nil
)
2656 (math-normalize (cons (car math-normalize-a
) args
)))
2657 (calc-record-why "*Variable is void" (nth 1 err
)))))
2658 (if (consp (car math-normalize-a
))
2659 (math-dimension-error)
2660 (cons (car math-normalize-a
) args
))))))))
2664 ;;; True if A is a floating-point real or complex number. [P x] [Public]
2665 (defun math-floatp (a)
2666 (cond ((eq (car-safe a
) 'float
) t
)
2667 ((memq (car-safe a
) '(cplx polar mod sdev intv
))
2668 (or (math-floatp (nth 1 a
))
2669 (math-floatp (nth 2 a
))
2670 (and (eq (car a
) 'intv
) (math-floatp (nth 3 a
)))))
2671 ((eq (car-safe a
) 'date
)
2672 (math-floatp (nth 1 a
)))))
2676 ;;; Verify that A is a complete object and return A. [x x] [Public]
2677 (defun math-check-complete (a)
2678 (cond ((integerp a
) a
)
2679 ((eq (car-safe a
) 'incomplete
)
2680 (calc-incomplete-error a
))
2682 (t (error "Invalid data object encountered"))))
2686 ;;; Coerce integer A to be a bignum. [B S]
2687 (defun math-bignum (a)
2689 (cons 'bigpos
(math-bignum-big a
))
2690 (cons 'bigneg
(math-bignum-big (- a
)))))
2692 (defun math-bignum-big (a) ; [L s]
2695 (cons (% a math-bignum-digit-size
)
2696 (math-bignum-big (/ a math-bignum-digit-size
)))))
2699 ;;; Build a normalized floating-point number. [F I S]
2700 (defun math-make-float (mant exp
)
2703 (let* ((ldiff (- calc-internal-prec
(math-numdigs mant
))))
2705 (setq mant
(math-scale-rounding mant ldiff
)
2706 exp
(- exp ldiff
))))
2708 (let ((digs (cdr mant
)))
2709 (if (= (%
(car digs
) 10) 0)
2711 (while (= (car digs
) 0)
2712 (setq digs
(cdr digs
)
2713 exp
(+ exp math-bignum-digit-length
)))
2714 (while (= (%
(car digs
) 10) 0)
2715 (setq digs
(math-div10-bignum digs
)
2717 (setq mant
(math-normalize (cons (car mant
) digs
))))))
2718 (while (= (% mant
10) 0)
2719 (setq mant
(/ mant
10)
2721 (if (and (<= exp -
4000000)
2722 (<= (+ exp
(math-numdigs mant
) -
1) -
4000000))
2723 (signal 'math-underflow nil
)
2724 (if (and (>= exp
3000000)
2725 (>= (+ exp
(math-numdigs mant
) -
1) 4000000))
2726 (signal 'math-overflow nil
)
2727 (list 'float mant exp
)))))
2729 (defun math-div10-bignum (a) ; [l l]
2731 (cons (+ (/ (car a
) 10) (* (%
(nth 1 a
) 10)
2732 (expt 10 (1- math-bignum-digit-length
))))
2733 (math-div10-bignum (cdr a
)))
2734 (list (/ (car a
) 10))))
2736 ;;; Coerce A to be a float. [F N; V V] [Public]
2737 (defun math-float (a)
2738 (cond ((Math-integerp a
) (math-make-float a
0))
2739 ((eq (car a
) 'frac
) (math-div (math-float (nth 1 a
)) (nth 2 a
)))
2740 ((eq (car a
) 'float
) a
)
2741 ((memq (car a
) '(cplx polar vec hms date sdev mod
))
2742 (cons (car a
) (mapcar 'math-float
(cdr a
))))
2743 (t (math-float-fancy a
))))
2747 (cond ((not (consp a
)) (- a
))
2748 ((eq (car a
) 'bigpos
) (cons 'bigneg
(cdr a
)))
2749 ((eq (car a
) 'bigneg
) (cons 'bigpos
(cdr a
)))
2750 ((memq (car a
) '(frac float
))
2751 (list (car a
) (Math-integer-neg (nth 1 a
)) (nth 2 a
)))
2752 ((memq (car a
) '(cplx vec hms date calcFunc-idn
))
2753 (cons (car a
) (mapcar 'math-neg
(cdr a
))))
2754 (t (math-neg-fancy a
))))
2757 ;;; Compute the number of decimal digits in integer A. [S I]
2758 (defun math-numdigs (a)
2761 (let* ((len (1- (length a
)))
2763 (+ (* (1- len
) math-bignum-digit-length
) (math-numdigs top
)))
2765 (cond ((>= a
100) (+ (math-numdigs (/ a
1000)) 3))
2771 (t (math-numdigs (- a
))))))
2773 ;;; Multiply (with truncation toward 0) the integer A by 10^N. [I i S]
2774 (defun math-scale-int (a n
)
2776 ((> n
0) (math-scale-left a n
))
2777 (t (math-normalize (math-scale-right a
(- n
))))))
2779 (defun math-scale-left (a n
) ; [I I S]
2783 (cons (car a
) (math-scale-left-bignum (cdr a
) n
))
2784 (if (>= n math-bignum-digit-length
)
2785 (if (or (>= a math-bignum-digit-size
)
2786 (<= a
(- math-bignum-digit-size
)))
2787 (math-scale-left (math-bignum a
) n
)
2788 (math-scale-left (* a math-bignum-digit-size
)
2789 (- n math-bignum-digit-length
)))
2790 (let ((sz (expt 10 (- (* 2 math-bignum-digit-length
) n
))))
2791 (if (or (>= a sz
) (<= a
(- sz
)))
2792 (math-scale-left (math-bignum a
) n
)
2793 (* a
(expt 10 n
))))))))
2795 (defun math-scale-left-bignum (a n
)
2796 (if (>= n math-bignum-digit-length
)
2797 (while (>= (setq a
(cons 0 a
)
2798 n
(- n math-bignum-digit-length
))
2799 math-bignum-digit-length
)))
2801 (math-mul-bignum-digit a
(expt 10 n
) 0)
2804 (defun math-scale-right (a n
) ; [i i S]
2808 (cons (car a
) (math-scale-right-bignum (cdr a
) n
))
2812 (- (math-scale-right (- a
) n
)))
2813 (if (>= n math-bignum-digit-length
)
2814 (while (and (> (setq a
(/ a math-bignum-digit-size
)) 0)
2815 (>= (setq n
(- n math-bignum-digit-length
))
2816 math-bignum-digit-length
))))
2821 (defun math-scale-right-bignum (a n
) ; [L L S; l l S]
2822 (if (>= n math-bignum-digit-length
)
2823 (setq a
(nthcdr (/ n math-bignum-digit-length
) a
)
2824 n
(% n math-bignum-digit-length
)))
2826 (cdr (math-mul-bignum-digit a
(expt 10 (- math-bignum-digit-length n
)) 0))
2829 ;;; Multiply (with rounding) the integer A by 10^N. [I i S]
2830 (defun math-scale-rounding (a n
)
2832 (math-scale-left a n
))
2836 (let ((val (if (< n
(- math-bignum-digit-length
))
2837 (math-scale-right-bignum
2839 (- (- math-bignum-digit-length
) n
))
2841 (math-mul-bignum-digit
2843 (expt 10 (+ math-bignum-digit-length n
)) 0)
2844 (cdr a
))))) ; n = -math-bignum-digit-length
2845 (if (and val
(>= (car val
) (/ math-bignum-digit-size
2)))
2847 (if (eq (car (cdr val
)) (1- math-bignum-digit-size
))
2848 (math-add-bignum (cdr val
) '(1))
2849 (cons (1+ (car (cdr val
))) (cdr (cdr val
))))
2854 (- (math-scale-rounding (- a
) n
))
2857 (/ (+ (math-scale-right a
(- -
1 n
)) 5) 10))))))
2860 ;;; Compute the sum of A and B. [O O O] [Public]
2861 (defun math-add (a b
)
2863 (and (not (or (consp a
) (consp b
)))
2866 (if (or (<= a
(- math-small-integer-size
)) (>= a math-small-integer-size
))
2869 (and (Math-zerop a
) (not (eq (car-safe a
) 'mod
))
2870 (if (and (math-floatp a
) (Math-ratp b
)) (math-float b
) b
))
2871 (and (Math-zerop b
) (not (eq (car-safe b
) 'mod
))
2872 (if (and (math-floatp b
) (Math-ratp a
)) (math-float a
) a
))
2873 (and (Math-objvecp a
) (Math-objvecp b
)
2875 (and (Math-integerp a
) (Math-integerp b
)
2877 (or (consp a
) (setq a
(math-bignum a
)))
2878 (or (consp b
) (setq b
(math-bignum b
)))
2879 (if (eq (car a
) 'bigneg
)
2880 (if (eq (car b
) 'bigneg
)
2881 (cons 'bigneg
(math-add-bignum (cdr a
) (cdr b
)))
2883 (let ((diff (math-sub-bignum (cdr b
) (cdr a
))))
2885 (cons 'bigneg
(math-sub-bignum (cdr a
) (cdr b
)))
2886 (cons 'bigpos diff
)))))
2887 (if (eq (car b
) 'bigneg
)
2889 (let ((diff (math-sub-bignum (cdr a
) (cdr b
))))
2891 (cons 'bigneg
(math-sub-bignum (cdr b
) (cdr a
)))
2892 (cons 'bigpos diff
))))
2893 (cons 'bigpos
(math-add-bignum (cdr a
) (cdr b
)))))))
2894 (and (Math-ratp a
) (Math-ratp b
)
2896 (calc-add-fractions a b
))
2897 (and (Math-realp a
) (Math-realp b
)
2899 (or (and (consp a
) (eq (car a
) 'float
))
2900 (setq a
(math-float a
)))
2901 (or (and (consp b
) (eq (car b
) 'float
))
2902 (setq b
(math-float b
)))
2903 (math-add-float a b
)))
2904 (and (require 'calc-ext
)
2905 (math-add-objects-fancy a b
))))
2906 (and (require 'calc-ext
)
2907 (math-add-symb-fancy a b
))))
2909 (defun math-add-bignum (a b
) ; [L L L; l l l]
2912 (let* ((a (copy-sequence a
)) (aa a
) (carry nil
) sum
)
2915 (if (< (setq sum
(+ (car aa
) (car b
)))
2916 (1- math-bignum-digit-size
))
2918 (setcar aa
(1+ sum
))
2920 (setcar aa
(- sum
(1- math-bignum-digit-size
))))
2921 (if (< (setq sum
(+ (car aa
) (car b
))) math-bignum-digit-size
)
2923 (setcar aa
(- sum math-bignum-digit-size
))
2929 (nconc a
(math-add-bignum b
'(1)))
2930 (while (eq (car aa
) (1- math-bignum-digit-size
))
2935 (setcar aa
(1+ (car aa
)))
2944 (defun math-sub-bignum (a b
) ; [l l l]
2947 (let* ((a (copy-sequence a
)) (aa a
) (borrow nil
) sum diff
)
2950 (if (>= (setq diff
(- (car aa
) (car b
))) 1)
2952 (setcar aa
(1- diff
))
2954 (setcar aa
(+ diff
(1- math-bignum-digit-size
))))
2955 (if (>= (setq diff
(- (car aa
) (car b
))) 0)
2957 (setcar aa
(+ diff math-bignum-digit-size
))
2963 (while (eq (car aa
) 0)
2964 (setcar aa
(1- math-bignum-digit-size
))
2968 (setcar aa
(1- (car aa
)))
2971 (while (eq (car b
) 0)
2976 (while (eq (car b
) 0)
2982 (defun math-add-float (a b
) ; [F F F]
2983 (let ((ediff (- (nth 2 a
) (nth 2 b
))))
2985 (if (>= ediff
(+ calc-internal-prec calc-internal-prec
))
2987 (math-make-float (math-add (nth 1 b
)
2990 (math-scale-left (nth 1 a
) ediff
)))
2992 (if (>= (setq ediff
(- ediff
))
2993 (+ calc-internal-prec calc-internal-prec
))
2995 (math-make-float (math-add (nth 1 a
)
2996 (math-scale-left (nth 1 b
) ediff
))
2999 ;;; Compute the difference of A and B. [O O O] [Public]
3000 (defun math-sub (a b
)
3001 (if (or (consp a
) (consp b
))
3002 (math-add a
(math-neg b
))
3004 (if (or (<= a
(- math-small-integer-size
)) (>= a math-small-integer-size
))
3008 (defun math-sub-float (a b
) ; [F F F]
3009 (let ((ediff (- (nth 2 a
) (nth 2 b
))))
3011 (if (>= ediff
(+ calc-internal-prec calc-internal-prec
))
3013 (math-make-float (math-add (Math-integer-neg (nth 1 b
))
3016 (math-scale-left (nth 1 a
) ediff
)))
3018 (if (>= (setq ediff
(- ediff
))
3019 (+ calc-internal-prec calc-internal-prec
))
3021 (math-make-float (math-add (nth 1 a
)
3023 (math-scale-left (nth 1 b
) ediff
)))
3027 ;;; Compute the product of A and B. [O O O] [Public]
3028 (defun math-mul (a b
)
3030 (and (not (consp a
)) (not (consp b
))
3031 (< a math-bignum-digit-size
) (> a
(- math-bignum-digit-size
))
3032 (< b math-bignum-digit-size
) (> b
(- math-bignum-digit-size
))
3034 (and (Math-zerop a
) (not (eq (car-safe b
) 'mod
))
3035 (if (Math-scalarp b
)
3036 (if (and (math-floatp b
) (Math-ratp a
)) (math-float a
) a
)
3038 (math-mul-zero a b
)))
3039 (and (Math-zerop b
) (not (eq (car-safe a
) 'mod
))
3040 (if (Math-scalarp a
)
3041 (if (and (math-floatp a
) (Math-ratp b
)) (math-float b
) b
)
3043 (math-mul-zero b a
)))
3044 (and (Math-objvecp a
) (Math-objvecp b
)
3046 (and (Math-integerp a
) (Math-integerp b
)
3048 (or (consp a
) (setq a
(math-bignum a
)))
3049 (or (consp b
) (setq b
(math-bignum b
)))
3051 (cons (if (eq (car a
) (car b
)) 'bigpos
'bigneg
)
3054 (math-mul-bignum (cdr a
) (cdr b
))
3055 (math-mul-bignum-digit (cdr a
) (nth 1 b
) 0))
3056 (math-mul-bignum-digit (cdr b
) (nth 1 a
) 0))))))
3057 (and (Math-ratp a
) (Math-ratp b
)
3059 (calc-mul-fractions a b
))
3060 (and (Math-realp a
) (Math-realp b
)
3062 (or (and (consp a
) (eq (car a
) 'float
))
3063 (setq a
(math-float a
)))
3064 (or (and (consp b
) (eq (car b
) 'float
))
3065 (setq b
(math-float b
)))
3066 (math-make-float (math-mul (nth 1 a
) (nth 1 b
))
3067 (+ (nth 2 a
) (nth 2 b
)))))
3068 (and (require 'calc-ext
)
3069 (math-mul-objects-fancy a b
))))
3070 (and (require 'calc-ext
)
3071 (math-mul-symb-fancy a b
))))
3073 (defun math-infinitep (a &optional undir
)
3074 (while (and (consp a
) (memq (car a
) '(* / neg
)))
3075 (if (or (not (eq (car a
) '*)) (math-infinitep (nth 1 a
)))
3077 (setq a
(nth 2 a
))))
3080 (memq (nth 2 a
) '(var-inf var-uinf var-nan
))
3081 (if (and undir
(eq (nth 2 a
) 'var-inf
))
3082 '(var uinf var-uinf
)
3085 ;;; Multiply digit lists A and B. [L L L; l l l]
3086 (defun math-mul-bignum (a b
)
3088 (let* ((sum (if (<= (car b
) 1)
3092 (math-mul-bignum-digit a
(car b
) 0)))
3093 (sump sum
) c d aa ss prod
)
3094 (while (setq b
(cdr b
))
3095 (setq ss
(setq sump
(or (cdr sump
) (setcdr sump
(list 0))))
3100 (setcar ss
(%
(setq prod
(+ (+ (car ss
) (* (car aa
) d
))
3101 c
)) math-bignum-digit-size
))
3103 (setq c
(/ prod math-bignum-digit-size
)
3104 ss
(or (cdr ss
) (setcdr ss
(list 0)))))
3105 (if (>= prod math-bignum-digit-size
)
3107 (setcar (cdr ss
) (+ (/ prod math-bignum-digit-size
) (car (cdr ss
))))
3108 (setcdr ss
(list (/ prod math-bignum-digit-size
))))))
3111 ;;; Multiply digit list A by digit D. [L L D D; l l D D]
3112 (defun math-mul-bignum-digit (a d c
)
3116 (let* ((a (copy-sequence a
)) (aa a
) prod
)
3119 (%
(setq prod
(+ (* (car aa
) d
) c
))
3120 math-bignum-digit-size
))
3123 c
(/ prod math-bignum-digit-size
)))
3124 (if (>= prod math-bignum-digit-size
)
3125 (setcdr aa
(list (/ prod math-bignum-digit-size
))))
3131 ;;; Compute the integer (quotient . remainder) of A and B, which may be
3132 ;;; small or big integers. Type and consistency of truncation is undefined
3133 ;;; if A or B is negative. B must be nonzero. [I.I I I] [Public]
3134 (defun math-idivmod (a b
)
3136 (math-reject-arg a
"*Division by zero"))
3137 (if (or (consp a
) (consp b
))
3138 (if (and (natnump b
) (< b math-bignum-digit-size
))
3139 (let ((res (math-div-bignum-digit (cdr a
) b
)))
3141 (math-normalize (cons (car a
) (car res
)))
3143 (or (consp a
) (setq a
(math-bignum a
)))
3144 (or (consp b
) (setq b
(math-bignum b
)))
3145 (let ((res (math-div-bignum (cdr a
) (cdr b
))))
3147 (math-normalize (cons (if (eq (car a
) (car b
)) 'bigpos
'bigneg
)
3149 (math-normalize (cons (car a
) (cdr res
))))))
3150 (cons (/ a b
) (% a b
))))
3152 (defun math-quotient (a b
) ; [I I I] [Public]
3153 (if (and (not (consp a
)) (not (consp b
)))
3155 (math-reject-arg a
"*Division by zero")
3157 (if (and (natnump b
) (< b math-bignum-digit-size
))
3159 (math-reject-arg a
"*Division by zero")
3160 (math-normalize (cons (car a
)
3161 (car (math-div-bignum-digit (cdr a
) b
)))))
3162 (or (consp a
) (setq a
(math-bignum a
)))
3163 (or (consp b
) (setq b
(math-bignum b
)))
3164 (let* ((alen (1- (length a
)))
3165 (blen (1- (length b
)))
3166 (d (/ math-bignum-digit-size
(1+ (nth (1- blen
) (cdr b
)))))
3167 (res (math-div-bignum-big (math-mul-bignum-digit (cdr a
) d
0)
3168 (math-mul-bignum-digit (cdr b
) d
0)
3170 (math-normalize (cons (if (eq (car a
) (car b
)) 'bigpos
'bigneg
)
3174 ;;; Divide a bignum digit list by another. [l.l l L]
3175 ;;; The following division algorithm is borrowed from Knuth vol. II, sec. 4.3.1
3176 (defun math-div-bignum (a b
)
3178 (let* ((alen (length a
))
3180 (d (/ math-bignum-digit-size
(1+ (nth (1- blen
) b
))))
3181 (res (math-div-bignum-big (math-mul-bignum-digit a d
0)
3182 (math-mul-bignum-digit b d
0)
3187 (car (math-div-bignum-digit (cdr res
) d
)))))
3188 (let ((res (math-div-bignum-digit a
(car b
))))
3189 (cons (car res
) (list (cdr res
))))))
3191 ;;; Divide a bignum digit list by a digit. [l.D l D]
3192 (defun math-div-bignum-digit (a b
)
3194 (let* ((res (math-div-bignum-digit (cdr a
) b
))
3195 (num (+ (* (cdr res
) math-bignum-digit-size
) (car a
))))
3197 (cons (/ num b
) (car res
))
3201 (defun math-div-bignum-big (a b alen blen
) ; [l.l l L]
3204 (let* ((res (math-div-bignum-big (cdr a
) b
(1- alen
) blen
))
3205 (num (cons (car a
) (cdr res
)))
3206 (res2 (math-div-bignum-part num b blen
)))
3208 (cons (car res2
) (car res
))
3211 (defun math-div-bignum-part (a b blen
) ; a < b*math-bignum-digit-size [D.l l L]
3212 (let* ((num (+ (* (or (nth blen a
) 0) math-bignum-digit-size
)
3213 (or (nth (1- blen
) a
) 0)))
3214 (den (nth (1- blen
) b
))
3215 (guess (min (/ num den
) (1- math-bignum-digit-size
))))
3216 (math-div-bignum-try a b
(math-mul-bignum-digit b guess
0) guess
)))
3218 (defun math-div-bignum-try (a b c guess
) ; [D.l l l D]
3219 (let ((rem (math-sub-bignum a c
)))
3221 (math-div-bignum-try a b
(math-sub-bignum c b
) (1- guess
))
3225 ;;; Compute the quotient of A and B. [O O N] [Public]
3226 (defun math-div (a b
)
3230 (math-div-by-zero a b
))
3231 (and (Math-zerop a
) (not (eq (car-safe b
) 'mod
))
3232 (if (Math-scalarp b
)
3233 (if (and (math-floatp b
) (Math-ratp a
)) (math-float a
) a
)
3235 (math-div-zero a b
)))
3236 (and (Math-objvecp a
) (Math-objvecp b
)
3238 (and (Math-integerp a
) (Math-integerp b
)
3239 (let ((q (math-idivmod a b
)))
3242 (if calc-prefer-frac
3245 (math-make-frac a b
))
3246 (math-div-float (math-make-float a
0)
3247 (math-make-float b
0))))))
3248 (and (Math-ratp a
) (Math-ratp b
)
3250 (calc-div-fractions a b
))
3251 (and (Math-realp a
) (Math-realp b
)
3253 (or (and (consp a
) (eq (car a
) 'float
))
3254 (setq a
(math-float a
)))
3255 (or (and (consp b
) (eq (car b
) 'float
))
3256 (setq b
(math-float b
)))
3257 (math-div-float a b
)))
3258 (and (require 'calc-ext
)
3259 (math-div-objects-fancy a b
))))
3260 (and (require 'calc-ext
)
3261 (math-div-symb-fancy a b
))))
3263 (defun math-div-float (a b
) ; [F F F]
3264 (let ((ldiff (max (- (1+ calc-internal-prec
)
3265 (- (math-numdigs (nth 1 a
)) (math-numdigs (nth 1 b
))))
3267 (math-make-float (math-quotient (math-scale-int (nth 1 a
) ldiff
) (nth 1 b
))
3268 (- (- (nth 2 a
) (nth 2 b
)) ldiff
))))
3273 (defvar calc-selection-cache-entry
)
3274 ;;; Format the number A as a string. [X N; X Z] [Public]
3275 (defun math-format-stack-value (entry)
3276 (setq calc-selection-cache-entry calc-selection-cache-default-entry
)
3277 (let* ((a (car entry
))
3278 (math-comp-selected (nth 2 entry
))
3279 (c (cond ((null a
) "<nil>")
3280 ((eq calc-display-raw t
) (format "%s" a
))
3282 ((eq a
'top-of-stack
) (propertize "." 'font-lock-face
'bold
))
3283 (calc-prepared-composition
3284 calc-prepared-composition
)
3285 ((and (Math-scalarp a
)
3286 (memq calc-language
'(nil flat unform
))
3287 (null math-comp-selected
))
3288 (math-format-number a
))
3289 (t (require 'calc-ext
)
3290 (math-compose-expr a
0))))
3291 (off (math-stack-value-offset c
))
3293 (and math-comp-selected
(setq calc-any-selections t
))
3297 (setq c
(math-comp-concat (make-string off ?\s
) c
)))
3298 (or (equal calc-left-label
"")
3299 (setq c
(math-comp-concat (if (eq a
'top-of-stack
)
3300 (make-string (length calc-left-label
) ?\s
)
3303 (when calc-line-numbering
3304 (setq c
(math-comp-concat (if (eq calc-language
'big
)
3305 (if math-comp-selected
3310 (unless (or (equal calc-right-label
"")
3311 (eq a
'top-of-stack
))
3313 (setq c
(list 'horiz c
3314 (make-string (max (- w
(math-comp-width c
)
3315 (length calc-right-label
)) 0) ?\s
)
3318 (setq s
(if (stringp c
)
3319 (if calc-display-raw
3322 (math-composition-to-string c w
)))
3323 (when calc-language-output-filter
3324 (setq s
(funcall calc-language-output-filter s
)))
3325 (if (eq calc-language
'big
)
3326 (setq s
(concat s
"\n"))
3327 (when calc-line-numbering
3328 (setq s
(concat "1:" (substring s
2)))))
3329 (setcar (cdr entry
) (calc-count-lines s
))
3332 ;; The variables math-svo-c, math-svo-wid and math-svo-off are local
3333 ;; to math-stack-value-offset, but are used by math-stack-value-offset-fancy
3336 (defun math-stack-value-offset (math-svo-c)
3337 (let* ((num (if calc-line-numbering
4 0))
3338 (math-svo-wid (calc-window-width))
3340 (if calc-display-just
3343 (math-stack-value-offset-fancy))
3344 (setq math-svo-off
(or calc-display-origin
0))
3345 (when (integerp calc-line-breaking
)
3346 (setq math-svo-wid calc-line-breaking
)))
3347 (cons (max (- math-svo-off
(length calc-left-label
)) 0)
3348 (+ math-svo-wid num
))))
3350 (defun calc-count-lines (s)
3353 (while (setq pos
(string-match "\n" s pos
))
3358 (defun math-format-value (a &optional w
)
3359 (if (and (Math-scalarp a
)
3360 (memq calc-language
'(nil flat unform
)))
3361 (math-format-number a
)
3363 (let ((calc-line-breaking nil
))
3364 (math-composition-to-string (math-compose-expr a
0) w
))))
3366 (defun calc-window-width ()
3367 (if calc-embedded-info
3368 (let ((win (get-buffer-window (aref calc-embedded-info
0))))
3369 (1- (if win
(window-width win
) (frame-width))))
3370 (- (window-width (get-buffer-window (current-buffer)))
3371 (if calc-line-numbering
5 1))))
3373 (defun math-comp-concat (c1 c2
)
3374 (if (and (stringp c1
) (stringp c2
))
3376 (list 'horiz c1 c2
)))
3380 ;;; Format an expression as a one-line string suitable for re-reading.
3382 (defun math-format-flat-expr (a prec
)
3384 ((or (not (or (consp a
) (integerp a
)))
3385 (eq calc-display-raw t
))
3386 (let ((print-escape-newlines t
))
3387 (concat "'" (prin1-to-string a
))))
3389 (let ((calc-group-digits nil
)
3390 (calc-point-char ".")
3391 (calc-frac-format (if (> (length (car calc-frac-format
)) 1)
3392 '("::" nil
) '(":" nil
)))
3393 (calc-complex-format nil
)
3394 (calc-hms-format "%s@ %s' %s\"")
3395 (calc-language nil
))
3396 (math-format-number a
)))
3399 (math-format-flat-expr-fancy a prec
))))
3403 ;;; Format a number as a string.
3404 (defun math-format-number (a &optional prec
) ; [X N] [Public]
3406 ((eq calc-display-raw t
) (format "%s" a
))
3407 ((and (nth 1 calc-frac-format
) (Math-integerp a
))
3409 (math-format-number (math-adjust-fraction a
)))
3411 (if (not (or calc-group-digits calc-leading-zeros
))
3412 (if (= calc-number-radix
10)
3415 (concat "-" (math-format-number (- a
)))
3417 (if math-radix-explicit-format
3418 (if calc-radix-formatter
3419 (funcall calc-radix-formatter
3421 (if (= calc-number-radix
2)
3422 (math-format-binary a
)
3423 (math-format-radix a
)))
3424 (format "%d#%s" calc-number-radix
3425 (if (= calc-number-radix
2)
3426 (math-format-binary a
)
3427 (math-format-radix a
))))
3428 (math-format-radix a
))))
3429 (math-format-number (math-bignum a
))))
3431 ((not (consp a
)) (prin1-to-string a
))
3432 ((eq (car a
) 'bigpos
) (math-format-bignum (cdr a
)))
3433 ((eq (car a
) 'bigneg
) (concat "-" (math-format-bignum (cdr a
))))
3434 ((and (eq (car a
) 'float
) (= calc-number-radix
10))
3435 (if (Math-integer-negp (nth 1 a
))
3436 (concat "-" (math-format-number (math-neg a
)))
3437 (let ((mant (nth 1 a
))
3439 (fmt (car calc-float-format
))
3440 (figs (nth 1 calc-float-format
))
3441 (point calc-point-char
)
3443 (if (and (eq fmt
'fix
)
3444 (or (and (< figs
0) (setq figs
(- figs
)))
3445 (> (+ exp
(math-numdigs mant
)) (- figs
))))
3447 (setq mant
(math-scale-rounding mant
(+ exp figs
))
3448 str
(if (integerp mant
)
3449 (int-to-string mant
)
3450 (math-format-bignum-decimal (cdr mant
))))
3451 (if (<= (length str
) figs
)
3452 (setq str
(concat (make-string (1+ (- figs
(length str
))) ?
0)
3455 (setq str
(concat (substring str
0 (- figs
)) point
3456 (substring str
(- figs
))))
3457 (setq str
(concat str point
)))
3458 (when calc-group-digits
3460 (setq str
(math-group-float str
))))
3462 (setq figs
(+ calc-internal-prec figs
)))
3464 (let ((adj (- figs
(math-numdigs mant
))))
3466 (setq mant
(math-scale-rounding mant adj
)
3468 (setq str
(if (integerp mant
)
3469 (int-to-string mant
)
3470 (math-format-bignum-decimal (cdr mant
))))
3471 (let* ((len (length str
))
3473 (if (and (eq fmt
'float
)
3474 (<= dpos
(+ calc-internal-prec calc-display-sci-high
))
3475 (>= dpos
(+ calc-display-sci-low
2)))
3479 (setq str
(concat "0" point str
)))
3480 ((and (<= exp
0) (> dpos
0))
3481 (setq str
(concat (substring str
0 dpos
) point
3482 (substring str dpos
))))
3484 (setq str
(concat str
(make-string exp ?
0) point
)))
3486 (setq str
(concat "0" point
3487 (make-string (- dpos
) ?
0) str
))))
3488 (when calc-group-digits
3490 (setq str
(math-group-float str
))))
3491 (let* ((eadj (+ exp len
))
3492 (scale (if (eq fmt
'eng
)
3493 (1+ (math-mod (+ eadj
300002) 3))
3495 (if (> scale
(length str
))
3496 (setq str
(concat str
(make-string (- scale
(length str
))
3498 (if (< scale
(length str
))
3499 (setq str
(concat (substring str
0 scale
) point
3500 (substring str scale
))))
3501 (when calc-group-digits
3503 (setq str
(math-group-float str
)))
3504 (setq str
(format (if (memq calc-language
'(math maple
))
3505 (if (and prec
(> prec
191))
3506 "(%s*10.^%d)" "%s*10.^%d")
3508 str
(- eadj scale
)))))))
3512 (math-format-number-fancy a prec
))))
3514 (defun math-format-bignum (a) ; [X L]
3515 (if (and (= calc-number-radix
10)
3516 (not calc-leading-zeros
)
3517 (not calc-group-digits
))
3518 (math-format-bignum-decimal a
)
3520 (math-format-bignum-fancy a
)))
3522 (defun math-format-bignum-decimal (a) ; [X L]
3525 (while (cdr (cdr a
))
3529 (number-to-string (* 2 math-bignum-digit-length
))
3531 (+ (* (nth 1 a
) math-bignum-digit-size
) (car a
))) s
)
3533 (concat (int-to-string
3534 (+ (* (or (nth 1 a
) 0) math-bignum-digit-size
) (car a
))) s
))
3539 ;;; Parse a simple number in string form. [N X] [Public]
3540 (defun math-read-number (s &optional decimal
)
3541 "Convert the string S into a Calc number."
3545 ;; Integers (most common case)
3546 ((string-match "\\` *\\([0-9]+\\) *\\'" s
)
3547 (let ((digs (math-match-substring s
1)))
3548 (if (and (memq calc-language calc-lang-c-type-hex
)
3550 (eq (aref digs
0) ?
0)
3552 (math-read-number (concat "8#" digs
))
3553 (if (<= (length digs
) (* 2 math-bignum-digit-length
))
3554 (string-to-number digs
)
3555 (cons 'bigpos
(math-read-bignum digs
))))))
3557 ;; Clean up the string if necessary
3558 ((string-match "\\`\\(.*\\)[ \t\n]+\\([^\001]*\\)\\'" s
)
3559 (math-read-number (concat (math-match-substring s
1)
3560 (math-match-substring s
2))))
3562 ;; Plus and minus signs
3563 ((string-match "^[-_+]\\(.*\\)$" s
)
3564 (let ((val (math-read-number (math-match-substring s
1))))
3565 (and val
(if (eq (aref s
0) ?
+) val
(math-neg val
)))))
3567 ;; Forms that require extensions module
3568 ((string-match "[^-+0-9eE.]" s
)
3570 (math-read-number-fancy s
))
3573 ((string-match "^\\([0-9]*\\)\\.\\([0-9]*\\)$" s
)
3574 (let ((int (math-match-substring s
1))
3575 (frac (math-match-substring s
2)))
3576 (let ((ilen (length int
))
3577 (flen (length frac
)))
3578 (let ((int (if (> ilen
0) (math-read-number int t
) 0))
3579 (frac (if (> flen
0) (math-read-number frac t
) 0)))
3580 (and int frac
(or (> ilen
0) (> flen
0))
3582 (math-add (math-scale-int int flen
) frac
)
3586 ((string-match "^\\(.*\\)[eE]\\([-+]?[0-9]+\\)$" s
)
3587 (let ((mant (math-match-substring s
1))
3588 (exp (math-match-substring s
2)))
3589 (let ((mant (if (> (length mant
) 0) (math-read-number mant t
) 1))
3590 (exp (if (<= (length exp
) (if (memq (aref exp
0) '(?
+ ?-
)) 8 7))
3591 (string-to-number exp
))))
3592 (and mant exp
(Math-realp mant
) (> exp -
4000000) (< exp
4000000)
3593 (let ((mant (math-float mant
)))
3594 (list 'float
(nth 1 mant
) (+ (nth 2 mant
) exp
)))))))
3599 ;;; Parse a very simple number, keeping all digits.
3600 (defun math-read-number-simple (s)
3601 "Convert the string S into a Calc number.
3602 S is assumed to be a simple number (integer or float without an exponent)
3603 and all digits are kept, regardless of Calc's current precision."
3606 ((string-match "^[0-9]+$" s
)
3607 (if (string-match "^\\(0+\\)" s
)
3608 (setq s
(substring s
(match-end 0))))
3609 (if (<= (length s
) (* 2 math-bignum-digit-length
))
3610 (string-to-number s
)
3611 (cons 'bigpos
(math-read-bignum s
))))
3613 ((string-match "^-[0-9]+$" s
)
3614 (if (<= (length s
) (1+ (* 2 math-bignum-digit-length
)))
3615 (string-to-number s
)
3616 (cons 'bigneg
(math-read-bignum (substring s
1)))))
3618 ((string-match "^\\(-?[0-9]*\\)\\.\\([0-9]*\\)$" s
)
3619 (let ((int (math-match-substring s
1))
3620 (frac (math-match-substring s
2)))
3621 (list 'float
(math-read-number-simple (concat int frac
))
3622 (- (length frac
)))))
3626 (defun math-match-substring (s n
)
3627 (if (match-beginning n
)
3628 (substring s
(match-beginning n
) (match-end n
))
3631 (defun math-read-bignum (s) ; [l X]
3632 (if (> (length s
) math-bignum-digit-length
)
3633 (cons (string-to-number (substring s
(- math-bignum-digit-length
)))
3634 (math-read-bignum (substring s
0 (- math-bignum-digit-length
))))
3635 (list (string-to-number s
))))
3637 (defconst math-standard-opers
3638 '( ( "_" calcFunc-subscr
1200 1201 )
3639 ( "%" calcFunc-percent
1100 -
1 )
3640 ( "u!" calcFunc-lnot -
1 1000 )
3641 ( "mod" mod
400 400 185 )
3642 ( "+/-" sdev
300 300 185 )
3643 ( "!!" calcFunc-dfact
210 -
1 )
3644 ( "!" calcFunc-fact
210 -
1 )
3647 ( "u+" ident -
1 197 )
3651 ( "\\" calcFunc-idiv
190 191 )
3655 ( "<" calcFunc-lt
160 161 )
3656 ( ">" calcFunc-gt
160 161 )
3657 ( "<=" calcFunc-leq
160 161 )
3658 ( ">=" calcFunc-geq
160 161 )
3659 ( "=" calcFunc-eq
160 161 )
3660 ( "==" calcFunc-eq
160 161 )
3661 ( "!=" calcFunc-neq
160 161 )
3662 ( "&&" calcFunc-land
110 111 )
3663 ( "||" calcFunc-lor
100 101 )
3664 ( "?" (math-read-if) 91 90 )
3665 ( "!!!" calcFunc-pnot -
1 85 )
3666 ( "&&&" calcFunc-pand
80 81 )
3667 ( "|||" calcFunc-por
75 76 )
3668 ( ":=" calcFunc-assign
51 50 )
3669 ( "::" calcFunc-condition
45 46 )
3670 ( "=>" calcFunc-evalto
40 41 )
3671 ( "=>" calcFunc-evalto
40 -
1 )))
3673 (defun math-standard-ops ()
3674 (if calc-multiplication-has-precedence
3679 math-standard-opers
))
3684 math-standard-opers
))))
3686 (defvar math-expr-opers
(math-standard-ops))
3688 (defun math-standard-ops-p ()
3689 (let ((meo (caar math-expr-opers
)))
3691 (string= meo
"*"))))
3693 (defun math-expr-ops ()
3694 (if (math-standard-ops-p)
3699 (defun calc-grab-region (top bot arg
)
3700 "Parse the region as a vector of numbers and push it on the Calculator stack."
3701 (interactive "r\nP")
3703 (calc-do-grab-region top bot arg
))
3706 (defun calc-grab-rectangle (top bot arg
)
3707 "Parse a rectangle as a matrix of numbers and push it on the Calculator stack."
3708 (interactive "r\nP")
3710 (calc-do-grab-rectangle top bot arg
))
3712 (defun calc-grab-sum-down (top bot arg
)
3713 "Parse a rectangle as a matrix of numbers and sum its columns."
3714 (interactive "r\nP")
3716 (calc-do-grab-rectangle top bot arg
'calcFunc-reduced
))
3718 (defun calc-grab-sum-across (top bot arg
)
3719 "Parse a rectangle as a matrix of numbers and sum its rows."
3720 (interactive "r\nP")
3722 (calc-do-grab-rectangle top bot arg
'calcFunc-reducea
))
3726 (defun calc-embedded (arg &optional end obeg oend
)
3727 "Start Calc Embedded mode on the formula surrounding point."
3730 (calc-do-embedded arg end obeg oend
))
3733 (defun calc-embedded-activate (&optional arg cbuf
)
3734 "Scan the current editing buffer for all embedded := and => formulas.
3735 Also looks for the equivalent TeX words, \\gets and \\evalto."
3737 (calc-do-embedded-activate arg cbuf
))
3739 (defun calc-user-invocation ()
3741 (unless calc-invocation-macro
3742 (error "Use `Z I' inside Calc to define a `C-x * Z' keyboard macro"))
3743 (execute-kbd-macro calc-invocation-macro nil
))
3745 ;;; User-programmability.
3748 (defmacro defmath
(func args
&rest body
) ; [Public]
3749 "Define Calc function.
3751 Like `defun' except that code in the body of the definition can
3752 make use of the full range of Calc data types and the usual
3753 arithmetic operations are converted to their Calc equivalents.
3755 The prefix `calcFunc-' is added to the specified name to get the
3756 actual Lisp function name.
3758 See Info node `(calc)Defining Functions'."
3759 (declare (doc-string 3))
3761 (math-do-defmath func args body
))
3763 ;;; Functions needed for Lucid Emacs support.
3765 (defun calc-read-key (&optional optkey
)
3766 (cond ((featurep 'xemacs
)
3767 (let ((event (next-command-event)))
3768 (let ((key (event-to-character event t t
)))
3769 (or key optkey
(error "Expected a plain keystroke"))
3772 (let ((key (read-event)))
3775 (defun calc-unread-command (&optional input
)
3776 (if (featurep 'xemacs
)
3777 (setq unread-command-event
3778 (if (integerp input
) (character-to-event input
)
3779 (or input last-command-event
)))
3780 (push (or input last-command-event
) unread-command-events
)))
3782 (defun calc-clear-unread-commands ()
3783 (if (featurep 'xemacs
)
3784 (setq unread-command-event nil
)
3785 (setq unread-command-events nil
)))
3787 (when calc-always-load-extensions
3789 (calc-load-everything))
3792 (run-hooks 'calc-load-hook
)
3796 ;; arch-tag: 0c3b170c-4ce6-4eaf-8d9b-5834d1fe938f
3797 ;;; calc.el ends here