1 ;;; calc-units.el --- unit conversion functions for Calc
3 ;; Copyright (C) 1990-1993, 2001-2011 Free Software Foundation, Inc.
5 ;; Author: David Gillespie <daveg@synaptics.com>
6 ;; Maintainer: Jay Belanger <jay.p.belanger@gmail.com>
8 ;; This file is part of GNU Emacs.
10 ;; GNU Emacs is free software: you can redistribute it and/or modify
11 ;; it under the terms of the GNU General Public License as published by
12 ;; the Free Software Foundation, either version 3 of the License, or
13 ;; (at your option) any later version.
15 ;; GNU Emacs is distributed in the hope that it will be useful,
16 ;; but WITHOUT ANY WARRANTY; without even the implied warranty of
17 ;; MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
18 ;; GNU General Public License for more details.
20 ;; You should have received a copy of the GNU General Public License
21 ;; along with GNU Emacs. If not, see <http://www.gnu.org/licenses/>.
27 ;; This file is autoloaded from calc-ext.el.
36 ;;; Units table last updated 9-Jan-91 by Ulrich Mueller (ulm@vsnhd1.cern.ch)
37 ;;; with some additions by Przemek Klosowski (przemek@rrdstrad.nist.gov)
38 ;;; Updated April 2002 by Jochen Küpper
40 ;;; Updated August 2007, using
41 ;;; CODATA (http://physics.nist.gov/cuu/Constants/index.html)
42 ;;; NIST (http://physics.nist.gov/Pubs/SP811/appenB9.html)
43 ;;; ESUWM (Encyclopaedia of Scientific Units, Weights and
44 ;;; Measures, by François Cardarelli)
45 ;;; All conversions are exact unless otherwise noted.
47 (defvar math-standard-units
50 ( in
"254*10^(-2) cm" "Inch" nil
54 ( mi
"5280 ft" "Mile" )
55 ( au
"149597870691. m" "Astronomical Unit" nil
57 ;; (approx) NASA JPL (http://neo.jpl.nasa.gov/glossary/au.html)
58 ( lyr
"c yr" "Light Year" )
59 ( pc
"3.0856775854*10^16 m" "Parsec (**)" nil
60 "3.0856775854 10^16 m (*)") ;; (approx) ESUWM
61 ( nmi
"1852 m" "Nautical Mile" )
62 ( fath
"6 ft" "Fathom" )
63 ( fur
"660 ft" "Furlong")
64 ( mu
"1 um" "Micron" )
65 ( mil
"(1/1000) in" "Mil" )
66 ( point
"(1/72) in" "Point (PostScript convention)" )
67 ( Ang
"10^(-10) m" "Angstrom" )
68 ( mfi
"mi+ft+in" "Miles + feet + inches" )
70 ( texpt
"(100/7227) in" "Point (TeX convention) (**)" )
71 ( texpc
"12 texpt" "Pica (TeX convention) (**)" )
72 ( texbp
"point" "Big point (TeX convention) (**)" )
73 ( texdd
"(1238/1157) texpt" "Didot point (TeX convention) (**)" )
74 ( texcc
"12 texdd" "Cicero (TeX convention) (**)" )
75 ( texsp
"(1/65536) texpt" "Scaled TeX point (TeX convention) (**)" )
78 ( hect
"10000 m^2" "*Hectare" )
80 ( acre
"(1/640) mi^2" "Acre" )
81 ( b
"10^(-28) m^2" "Barn" )
84 ( L
"10^(-3) m^3" "*Liter" )
86 ( gal
"4 qt" "US Gallon" )
88 ( pt
"2 cup" "Pint (**)" )
89 ( cup
"8 ozfl" "Cup" )
90 ( ozfl
"2 tbsp" "Fluid Ounce" )
91 ( floz
"2 tbsp" "Fluid Ounce" )
92 ( tbsp
"3 tsp" "Tablespoon" )
93 ;; ESUWM defines a US gallon as 231 in^3.
94 ;; That gives the following exact value for tsp.
95 ( tsp
"492892159375*10^(-11) ml" "Teaspoon" nil
97 ( vol
"tsp+tbsp+ozfl+cup+pt+qt+gal" "Gallons + ... + teaspoons" nil
98 "tsp+tbsp+ozfl+cup+pt+qt+gal")
99 ( galC
"galUK" "Canadian Gallon" )
100 ( galUK
"454609*10^(-5) L" "UK Gallon" nil
106 ( min
"60 s" "Minute" )
107 ( hr
"60 min" "Hour" )
108 ( day
"24 hr" "Day" )
109 ( wk
"7 day" "Week" )
110 ( hms
"wk+day+hr+min+s" "Hours, minutes, seconds" )
111 ( yr
"36525*10^(-2) day" "Year (Julian)" nil
116 ( mph
"mi/hr" "*Miles per hour" )
117 ( kph
"km/hr" "Kilometers per hour" )
118 ( knot
"nmi/hr" "Knot" )
119 ( c
"299792458 m/s" "Speed of light" ) ;;; CODATA
122 ( ga
"980665*10^(-5) m/s^2" "*\"g\" acceleration" nil
123 "9.80665 m / s^2") ;; CODATA
127 ( lb
"16 oz" "Pound (mass)" )
128 ( oz
"28349523125*10^(-9) g" "Ounce (mass)" nil
129 "28.349523125 g") ;; ESUWM
130 ( ton
"2000 lb" "Ton" )
131 ( tpo
"ton+lb+oz" "Tons + pounds + ounces (mass)" )
132 ( t
"1000 kg" "Metric ton" )
133 ( tonUK
"10160469088*10^(-7) kg" "UK ton" nil
134 "1016.0469088 kg") ;; ESUWM
135 ( lbt
"12 ozt" "Troy pound" )
136 ( ozt
"311034768*10^(-7) g" "Troy ounce" nil
137 "31.10347680 g") ;; ESUWM, 1/12 exact value for lbt
138 ( ct
"(2/10) g" "Carat" nil
140 ( u
"1.660538782*10^(-27) kg" "Unified atomic mass" nil
141 "1.660538782 10^-27 kg (*)");;(approx) CODATA
144 ( N
"m kg/s^2" "*Newton" )
145 ( dyn
"10^(-5) N" "Dyne" )
146 ( gf
"ga g" "Gram (force)" )
147 ( lbf
"ga lb" "Pound (force)" )
148 ( kip
"1000 lbf" "Kilopound (force)" )
149 ( pdl
"138254954376*10^(-12) N" "Poundal" nil
150 "0.138254954376 N") ;; ESUWM
154 ( erg
"10^(-7) J" "Erg" )
155 ( cal
"41868*10^(-4) J" "International Table Calorie" nil
157 ( calth
"4184*10^(-3) J" "Thermochemical Calorie" nil
159 ( Cal
"1000 cal" "Large Calorie")
160 ( Btu
"105505585262*10^(-8) J" "International Table Btu" nil
161 "1055.05585262 J") ;; ESUWM
162 ( eV
"ech V" "Electron volt" )
163 ( ev
"eV" "Electron volt" )
164 ( therm
"105506000 J" "EEC therm" )
165 ( invcm
"h c/cm" "Energy in inverse centimeters" )
166 ( Kayser
"invcm" "Kayser (inverse centimeter energy)" )
167 ( men
"100/invcm" "Inverse energy in meters" )
168 ( Hzen
"h Hz" "Energy in Hertz")
169 ( Ken
"k K" "Energy in Kelvins")
170 ( Wh
"W hr" "Watt hour")
171 ( Ws
"W s" "Watt second")
175 ( hp
"550 ft lbf/s" "Horsepower") ;;ESUWM
176 ( hpm
"75 m kgf/s" "Metric Horsepower") ;;ESUWM
179 ( K nil
"*Degree Kelvin" K
)
180 ( dK
"K" "Degree Kelvin" K
)
181 ( degK
"K" "Degree Kelvin" K
)
182 ( dC
"K" "Degree Celsius" C
)
183 ( degC
"K" "Degree Celsius" C
)
184 ( dF
"(5/9) K" "Degree Fahrenheit" F
)
185 ( degF
"(5/9) K" "Degree Fahrenheit" F
)
188 ( Pa
"N/m^2" "*Pascal" )
189 ( bar
"10^5 Pa" "Bar" )
190 ( atm
"101325 Pa" "Standard atmosphere" ) ;; CODATA
191 ( Torr
"(1/760) atm" "Torr")
192 ( mHg
"1000 Torr" "Meter of mercury" )
193 ( inHg
"254*10^(-1) mmHg" "Inch of mercury" nil
195 ( inH2O
"2.490889*10^2 Pa" "Inch of water" nil
196 "2.490889 10^2 Pa (*)") ;;(approx) NIST
197 ( psi
"lbf/in^2" "Pounds per square inch" )
200 ( P
"(1/10) Pa s" "*Poise" )
201 ( St
"10^(-4) m^2/s" "Stokes" )
205 ( C
"A s" "Coulomb" )
206 ( Fdy
"ech Nav" "Faraday" )
207 ( e
"ech" "Elementary charge" )
208 ( ech
"1.602176487*10^(-19) C" "Elementary charge" nil
209 "1.602176487 10^-19 C (*)") ;;(approx) CODATA
214 ( S
"A/V" "Siemens" )
217 ( T
"Wb/m^2" "Tesla" )
218 ( Gs
"10^(-4) T" "Gauss" )
221 ;; Luminous intensity
222 ( cd nil
"*Candela" )
223 ( sb
"10000 cd/m^2" "Stilb" )
224 ( lm
"cd sr" "Lumen" )
225 ( lx
"lm/m^2" "Lux" )
226 ( ph
"10000 lx" "Phot" )
227 ( fc
"lm/ft^2" "Footcandle") ;; ESUWM
228 ( lam
"10000 lm/m^2" "Lambert" )
229 ( flam
"(1/pi) cd/ft^2" "Footlambert") ;; ESUWM
232 ( Bq
"1/s" "*Becquerel" )
233 ( Ci
"37*10^9 Bq" "Curie" ) ;; ESUWM
235 ( Sv
"Gy" "Sievert" )
236 ( R
"258*10^(-6) C/kg" "Roentgen" ) ;; NIST
237 ( rd
"(1/100) Gy" "Rad" )
240 ;; Amount of substance
244 ( rad nil
"*Radian" )
245 ( circ
"2 pi rad" "Full circle" )
246 ( rev
"circ" "Full revolution" )
247 ( deg
"circ/360" "Degree" )
248 ( arcmin
"deg/60" "Arc minute" )
249 ( arcsec
"arcmin/60" "Arc second" )
250 ( grad
"circ/400" "Grade" )
251 ( rpm
"rev/min" "Revolutions per minute" )
254 ( sr nil
"*Steradian" )
256 ;; Other physical quantities
257 ;; The values are from CODATA, and are approximate.
258 ( h
"6.62606896*10^(-34) J s" "*Planck's constant" nil
259 "6.62606896 10^-34 J s (*)")
260 ( hbar
"h / (2 pi)" "Planck's constant" ) ;; Exact
261 ( mu0
"4 pi 10^(-7) H/m" "Permeability of vacuum") ;; Exact
262 ( μ
0 "mu0" "Permeability of vacuum") ;; Exact
263 ( eps0
"1 / (mu0 c^2)" "Permittivity of vacuum" )
264 ( ε
0 "eps0" "Permittivity of vacuum" )
265 ( G
"6.67428*10^(-11) m^3/(kg s^2)" "Gravitational constant" nil
266 "6.67428 10^-11 m^3/(kg s^2) (*)")
267 ( Nav
"6.02214179*10^(23) / mol" "Avogadro's constant" nil
268 "6.02214179 10^23 / mol (*)")
269 ( me
"9.10938215*10^(-31) kg" "Electron rest mass" nil
270 "9.10938215 10^-31 kg (*)")
271 ( mp
"1.672621637*10^(-27) kg" "Proton rest mass" nil
272 "1.672621637 10^-27 kg (*)")
273 ( mn
"1.674927211*10^(-27) kg" "Neutron rest mass" nil
274 "1.674927211 10^-27 kg (*)")
275 ( mmu
"1.88353130*10^(-28) kg" "Muon rest mass" nil
276 "1.88353130 10^-28 kg (*)")
277 ( mμ
"mmu" "Muon rest mass" nil
278 "1.88353130 10^-28 kg (*)")
279 ( Ryd
"10973731.568527 /m" "Rydberg's constant" nil
280 "10973731.568527 /m (*)")
281 ( k
"1.3806504*10^(-23) J/K" "Boltzmann's constant" nil
282 "1.3806504 10^-23 J/K (*)")
283 ( alpha
"7.2973525376*10^(-3)" "Fine structure constant" nil
284 "7.2973525376 10^-3 (*)")
285 ( α
"alpha" "Fine structure constant" nil
286 "7.2973525376 10^-3 (*)")
287 ( muB
"927.400915*10^(-26) J/T" "Bohr magneton" nil
288 "927.400915 10^-26 J/T (*)")
289 ( muN
"5.05078324*10^(-27) J/T" "Nuclear magneton" nil
290 "5.05078324 10^-27 J/T (*)")
291 ( mue
"-928.476377*10^(-26) J/T" "Electron magnetic moment" nil
292 "-928.476377 10^-26 J/T (*)")
293 ( mup
"1.410606662*10^(-26) J/T" "Proton magnetic moment" nil
294 "1.410606662 10^-26 J/T (*)")
295 ( R0
"8.314472 J/(mol K)" "Molar gas constant" nil
296 "8.314472 J/(mol K) (*)")
297 ( V0
"22.710981*10^(-3) m^3/mol" "Standard volume of ideal gas" nil
298 "22.710981 10^-3 m^3/mol (*)")
301 ( dB
"(ln(10)/20) Np" "decibel")))
304 (defvar math-additional-units nil
305 "*Additional units table for user-defined units.
306 Must be formatted like `math-standard-units'.
307 If you change this, be sure to set `math-units-table' to nil to ensure
308 that the combined units table will be rebuilt.")
310 (defvar math-unit-prefixes
311 '( ( ?Y
(^
10 24) "Yotta" )
312 ( ?Z
(^
10 21) "Zetta" )
313 ( ?E
(^
10 18) "Exa" )
314 ( ?P
(^
10 15) "Peta" )
315 ( ?T
(^
10 12) "Tera" )
316 ( ?G
(^
10 9) "Giga" )
317 ( ?M
(^
10 6) "Mega" )
318 ( ?k
(^
10 3) "Kilo" )
319 ( ?K
(^
10 3) "Kilo" )
320 ( ?h
(^
10 2) "Hecto" )
321 ( ?H
(^
10 2) "Hecto" )
322 ( ?D
(^
10 1) "Deka" )
324 ( ?d
(^
10 -
1) "Deci" )
325 ( ?c
(^
10 -
2) "Centi" )
326 ( ?m
(^
10 -
3) "Milli" )
327 ( ?u
(^
10 -
6) "Micro" )
328 ( ?μ
(^
10 -
6) "Micro" )
329 ( ?n
(^
10 -
9) "Nano" )
330 ( ?p
(^
10 -
12) "Pico" )
331 ( ?f
(^
10 -
15) "Femto" )
332 ( ?a
(^
10 -
18) "Atto" )
333 ( ?z
(^
10 -
21) "zepto" )
334 ( ?y
(^
10 -
24) "yocto" )))
336 (defvar math-standard-units-systems
338 ( si
( ( g
'(/ (var kg var-kg
) 1000) ) ) )
339 ( mks
( ( g
'(/ (var kg var-kg
) 1000) ) ) )
340 ( cgs
( ( m
'(* (var cm var-cm
) 100 ) ) ) )))
342 (defvar math-units-table nil
343 "Internal units table.
344 Derived from `math-standard-units' and `math-additional-units'.
345 Entries are (SYMBOL EXPR DOC-STRING TEMP-TYPE BASE-UNITS).")
347 (defvar math-units-table-buffer-valid nil
)
351 (defun calc-base-units ()
354 (let ((calc-autorange-units nil
))
355 (calc-enter-result 1 "bsun" (math-simplify-units
356 (math-to-standard-units (calc-top-n 1)
359 (defun calc-quick-units ()
362 (let* ((num (- last-command-event ?
0))
363 (pos (if (= num
0) 10 num
))
364 (units (calc-var-value 'var-Units
))
365 (expr (calc-top-n 1)))
366 (unless (and (>= num
0) (<= num
9))
367 (error "Bad unit number"))
368 (unless (math-vectorp units
)
369 (error "No \"quick units\" are defined"))
370 (unless (< pos
(length units
))
371 (error "Unit number %d not defined" pos
))
372 (if (math-units-in-expr-p expr nil
)
373 (calc-enter-result 1 (format "cun%d" num
)
374 (math-convert-units expr
(nth pos units
)))
375 (calc-enter-result 1 (format "*un%d" num
)
377 (math-mul expr
(nth pos units
))))))))
379 (defun math-get-standard-units (expr)
380 "Return the standard units in EXPR."
383 (math-to-standard-units expr nil
))))
385 (defun math-get-units (expr)
386 "Return the units in EXPR."
388 (math-extract-units expr
)))
390 (defun math-make-unit-string (expr)
391 "Return EXPR in string form.
392 If EXPR is nil, return nil."
394 (let ((cexpr (math-compose-expr expr
0)))
395 (replace-regexp-in-string
399 (math-composition-to-string cexpr
))))))
401 (defvar math-default-units-table
402 (make-hash-table :test
'equal
)
403 "A table storing previously converted units.")
405 (defun math-get-default-units (expr)
406 "Get default units to use when converting the units in EXPR."
407 (let* ((units (math-get-units expr
))
408 (standard-units (math-get-standard-units expr
))
409 (default-units (gethash
411 math-default-units-table
)))
412 (if (equal units
(car default-units
))
413 (math-make-unit-string (cadr default-units
))
414 (math-make-unit-string (car default-units
)))))
416 (defun math-put-default-units (expr)
417 "Put the units in EXPR in the default units table."
418 (let* ((units (math-get-units expr
))
419 (standard-units (math-get-standard-units expr
))
420 (default-units (gethash
422 math-default-units-table
)))
425 (puthash standard-units
(list units
) math-default-units-table
))
426 ((not (equal units
(car default-units
)))
427 (puthash standard-units
428 (list units
(car default-units
))
429 math-default-units-table
)))))
432 (defun calc-convert-units (&optional old-units new-units
)
435 (let ((expr (calc-top-n 1))
440 (unless (math-units-in-expr-p expr t
)
441 (let ((uold (or old-units
443 (setq uoldname
(read-string "Old units: "))
444 (if (equal uoldname
"")
448 (if (string-match "\\` */" uoldname
)
449 (setq uoldname
(concat "1" uoldname
)))
450 (math-read-expr uoldname
))))))
451 (when (eq (car-safe uold
) 'error
)
452 (error "Bad format in units expression: %s" (nth 1 uold
)))
453 (setq expr
(math-mul expr uold
))))
455 (setq defunits
(math-get-default-units expr
))
459 (concat "Old units: "
471 (string= new-units
"")
473 (setq new-units defunits
)))
474 (when (string-match "\\` */" new-units
)
475 (setq new-units
(concat "1" new-units
)))
476 (setq units
(math-read-expr new-units
))
477 (when (eq (car-safe units
) 'error
)
478 (error "Bad format in units expression: %s" (nth 2 units
)))
479 (math-put-default-units units
)
480 (let ((unew (math-units-in-expr-p units t
))
481 (std (and (eq (car-safe units
) 'var
)
482 (assq (nth 1 units
) math-standard-units-systems
))))
484 (calc-enter-result 1 "cvun" (math-simplify-units
485 (math-to-standard-units expr
488 (error "No units specified"))
489 (calc-enter-result 1 "cvun"
492 (and uoldname
(not (equal uoldname
"1"))))))))))
494 (defun calc-autorange-units (arg)
497 (calc-change-mode 'calc-autorange-units arg nil t
)
498 (message (if calc-autorange-units
499 "Adjusting target unit prefix automatically"
500 "Using target units exactly"))))
502 (defun calc-convert-temperature (&optional old-units new-units
)
505 (let ((expr (calc-top-n 1))
510 (setq uold
(or old-units
511 (let ((units (math-single-units-in-expr-p expr
)))
514 (list 'var
(car units
)
515 (intern (concat "var-"
518 (error "Not a pure temperature expression"))
520 (setq uoldname
(read-string
521 "Old temperature units: ")))))))
522 (when (eq (car-safe uold
) 'error
)
523 (error "Bad format in units expression: %s" (nth 2 uold
)))
524 (or (math-units-in-expr-p expr nil
)
525 (setq expr
(math-mul expr uold
)))
526 (setq defunits
(math-get-default-units expr
))
527 (setq unew
(or new-units
531 (concat "Old temperature units: "
534 "New temperature units")
540 (setq unew
(math-read-expr (if (string= unew
"") defunits unew
)))
541 (when (eq (car-safe unew
) 'error
)
542 (error "Bad format in units expression: %s" (nth 2 unew
)))
543 (math-put-default-units unew
)
544 (let ((ntemp (calc-normalize
546 (math-convert-temperature expr uold unew
548 (if (Math-zerop ntemp
)
549 (setq ntemp
(list '* ntemp unew
)))
550 (let ((calc-simplify-mode 'none
))
551 (calc-enter-result 1 "cvtm" ntemp
))))))
553 (defun calc-remove-units ()
556 (calc-enter-result 1 "rmun" (math-simplify-units
557 (math-remove-units (calc-top-n 1))))))
559 (defun calc-extract-units ()
562 (calc-enter-result 1 "rmun" (math-simplify-units
563 (math-extract-units (calc-top-n 1))))))
565 ;; The variables calc-num-units and calc-den-units are local to
566 ;; calc-explain-units, but are used by calc-explain-units-rec,
567 ;; which is called by calc-explain-units.
568 (defvar calc-num-units
)
569 (defvar calc-den-units
)
571 (defun calc-explain-units ()
574 (let ((calc-num-units nil
)
575 (calc-den-units nil
))
576 (calc-explain-units-rec (calc-top-n 1) 1)
577 (and calc-den-units
(string-match "^[^(].* .*[^)]$" calc-den-units
)
578 (setq calc-den-units
(concat "(" calc-den-units
")")))
581 (message "%s per %s" calc-num-units calc-den-units
)
582 (message "%s" calc-num-units
))
584 (message "1 per %s" calc-den-units
)
585 (message "No units in expression"))))))
587 (defun calc-explain-units-rec (expr pow
)
588 (let ((u (math-check-unit-name expr
))
590 (if (and u
(not (math-zerop pow
)))
591 (let ((name (or (nth 2 u
) (symbol-name (car u
)))))
592 (if (eq (aref name
0) ?\
*)
593 (setq name
(substring name
1)))
594 (if (string-match "[^a-zA-Zα-ωΑ-Ω0-9']" name
)
595 (if (string-match "^[a-zA-Zα-ωΑ-Ω0-9' ()]*$" name
)
596 (while (setq pos
(string-match "[ ()]" name
))
597 (setq name
(concat (substring name
0 pos
)
598 (if (eq (aref name pos
) 32) "-" "")
599 (substring name
(1+ pos
)))))
600 (setq name
(concat "(" name
")"))))
601 (or (eq (nth 1 expr
) (car u
))
602 (setq name
(concat (nth 2 (assq (aref (symbol-name
605 (if (and (string-match "[^a-zA-Zα-ωΑ-Ω0-9']" name
)
606 (not (memq (car u
) '(mHg gf
))))
609 (cond ((or (math-equal-int pow
1)
610 (math-equal-int pow -
1)))
611 ((or (math-equal-int pow
2)
612 (math-equal-int pow -
2))
613 (if (equal (nth 4 u
) '((m .
1)))
614 (setq name
(concat "Square-" name
))
615 (setq name
(concat name
"-squared"))))
616 ((or (math-equal-int pow
3)
617 (math-equal-int pow -
3))
618 (if (equal (nth 4 u
) '((m .
1)))
619 (setq name
(concat "Cubic-" name
))
620 (setq name
(concat name
"-cubed"))))
622 (setq name
(concat name
"^"
623 (math-format-number (math-abs pow
))))))
625 (setq calc-num-units
(if calc-num-units
626 (concat calc-num-units
" " name
)
628 (setq calc-den-units
(if calc-den-units
629 (concat calc-den-units
" " name
)
631 (cond ((eq (car-safe expr
) '*)
632 (calc-explain-units-rec (nth 1 expr
) pow
)
633 (calc-explain-units-rec (nth 2 expr
) pow
))
634 ((eq (car-safe expr
) '/)
635 (calc-explain-units-rec (nth 1 expr
) pow
)
636 (calc-explain-units-rec (nth 2 expr
) (- pow
)))
637 ((memq (car-safe expr
) '(neg + -
))
638 (calc-explain-units-rec (nth 1 expr
) pow
))
639 ((and (eq (car-safe expr
) '^
)
640 (math-realp (nth 2 expr
)))
641 (calc-explain-units-rec (nth 1 expr
)
642 (math-mul pow
(nth 2 expr
))))))))
644 (defun calc-simplify-units ()
647 (calc-with-default-simplification
648 (calc-enter-result 1 "smun" (math-simplify-units (calc-top-n 1))))))
650 (defun calc-view-units-table (n)
652 (and n
(setq math-units-table-buffer-valid nil
))
653 (let ((win (get-buffer-window "*Units Table*")))
656 math-units-table-buffer-valid
)
658 (bury-buffer (window-buffer win
))
659 (let ((curwin (selected-window)))
661 (switch-to-buffer nil
)
662 (select-window curwin
)))
663 (math-build-units-table-buffer nil
))))
665 (defun calc-enter-units-table (n)
667 (and n
(setq math-units-table-buffer-valid nil
))
668 (math-build-units-table-buffer t
)
669 (message "%s" (substitute-command-keys "Type \\[calc] to return to the Calculator")))
671 (defun calc-define-unit (uname desc
&optional disp
)
672 (interactive "SDefine unit name: \nsDescription: \nP")
673 (if disp
(setq disp
(read-string "Display definition: ")))
675 (let ((form (calc-top-n 1))
676 (unit (assq uname math-additional-units
)))
678 (setq math-additional-units
679 (cons (setq unit
(list uname nil nil nil nil
))
680 math-additional-units
)
681 math-units-table nil
))
682 (setcar (cdr unit
) (and (not (and (eq (car-safe form
) 'var
)
683 (eq (nth 1 form
) uname
)))
684 (not (math-equal-int form
1))
685 (math-format-flat-expr form
0)))
686 (setcar (cdr (cdr unit
)) (and (not (equal desc
""))
689 (setcar (cdr (cdr (cdr (cdr unit
)))) disp
))))
690 (calc-invalidate-units-table))
692 (defun calc-undefine-unit (uname)
693 (interactive "SUndefine unit name: ")
695 (let ((unit (assq uname math-additional-units
)))
697 (if (assq uname math-standard-units
)
698 (error "\"%s\" is a predefined unit name" uname
)
699 (error "Unit name \"%s\" not found" uname
)))
700 (setq math-additional-units
(delq unit math-additional-units
)
701 math-units-table nil
)))
702 (calc-invalidate-units-table))
704 (defun calc-invalidate-units-table ()
705 (setq math-units-table nil
)
706 (let ((buf (get-buffer "*Units Table*")))
708 (with-current-buffer buf
710 (goto-char (point-min))
711 (if (looking-at "Calculator Units Table")
712 (let ((inhibit-read-only t
))
713 (insert "(Obsolete) "))))))))
715 (defun calc-get-unit-definition (uname)
716 (interactive "SGet definition for unit: ")
718 (math-build-units-table)
719 (let ((unit (assq uname math-units-table
)))
721 (error "Unit name \"%s\" not found" uname
))
722 (let ((msg (nth 2 unit
)))
724 (if (string-match "^\\*" msg
)
725 (setq msg
(substring msg
1)))
726 (setq msg
(symbol-name uname
)))
729 (calc-enter-result 0 "ugdf" (nth 1 unit
))
730 (message "Derived unit: %s" msg
))
731 (calc-enter-result 0 "ugdf" (list 'var uname
734 (symbol-name uname
)))))
735 (message "Base unit: %s" msg
))))))
737 (defun calc-permanent-units ()
741 (set-buffer (find-file-noselect (substitute-in-file-name
742 calc-settings-file
)))
743 (goto-char (point-min))
744 (if (and (search-forward ";;; Custom units stored by Calc" nil t
)
748 (search-forward "\n;;; End of custom units" nil t
)))
752 (delete-region pos
(point)))
753 (goto-char (point-max))
756 (insert ";;; Custom units stored by Calc on " (current-time-string) "\n")
757 (if math-additional-units
759 (insert "(setq math-additional-units '(\n")
760 (let ((list math-additional-units
))
762 (insert " (" (symbol-name (car (car list
))) " "
763 (if (nth 1 (car list
))
764 (if (stringp (nth 1 (car list
)))
765 (prin1-to-string (nth 1 (car list
)))
766 (prin1-to-string (math-format-flat-expr
767 (nth 1 (car list
)) 0)))
770 (prin1-to-string (nth 2 (car list
)))
772 (setq list
(cdr list
))))
774 (insert ";;; (no custom units defined)\n"))
775 (insert ";;; End of custom units\n")
779 ;; The variable math-cu-unit-list is local to math-build-units-table,
780 ;; but is used by math-compare-unit-names, which is called (indirectly)
781 ;; by math-build-units-table.
782 ;; math-cu-unit-list is also local to math-convert-units, but is used
783 ;; by math-convert-units-rec, which is called by math-convert-units.
784 (defvar math-cu-unit-list
)
786 (defun math-build-units-table ()
788 (let* ((combined-units (append math-additional-units
789 math-standard-units
))
790 (math-cu-unit-list (mapcar 'car combined-units
))
792 (message "Building units table...")
793 (setq math-units-table-buffer-valid nil
)
794 (setq tab
(mapcar (function
798 (if (stringp (nth 1 x
))
799 (let ((exp (math-read-plain-expr
801 (if (eq (car-safe exp
) 'error
)
802 (error "Format error in definition of %s in units table: %s"
809 (list (cons (car x
) 1)))
812 (let ((math-units-table tab
))
813 (mapc 'math-find-base-units tab
))
814 (message "Building units table...done")
815 (setq math-units-table tab
))))
817 ;; The variables math-fbu-base and math-fbu-entry are local to
818 ;; math-find-base-units, but are used by math-find-base-units-rec,
819 ;; which is called by math-find-base-units.
820 (defvar math-fbu-base
)
821 (defvar math-fbu-entry
)
823 (defun math-find-base-units (math-fbu-entry)
824 (if (eq (nth 4 math-fbu-entry
) 'boom
)
825 (error "Circular definition involving unit %s" (car math-fbu-entry
)))
826 (or (nth 4 math-fbu-entry
)
828 (setcar (nthcdr 4 math-fbu-entry
) 'boom
)
829 (math-find-base-units-rec (nth 1 math-fbu-entry
) 1)
831 (error "Dimensionless definition for unit %s" (car math-fbu-entry
)))
832 (while (eq (cdr (car math-fbu-base
)) 0)
833 (setq math-fbu-base
(cdr math-fbu-base
)))
834 (let ((b math-fbu-base
))
836 (if (eq (cdr (car (cdr b
))) 0)
837 (setcdr b
(cdr (cdr b
)))
839 (setq math-fbu-base
(sort math-fbu-base
'math-compare-unit-names
))
840 (setcar (nthcdr 4 math-fbu-entry
) math-fbu-base
)
843 (defun math-compare-unit-names (a b
)
844 (memq (car b
) (cdr (memq (car a
) math-cu-unit-list
))))
846 (defun math-find-base-units-rec (expr pow
)
847 (let ((u (math-check-unit-name expr
)))
849 (let ((ulist (math-find-base-units u
)))
851 (let ((p (* (cdr (car ulist
)) pow
))
852 (old (assq (car (car ulist
)) math-fbu-base
)))
854 (setcdr old
(+ (cdr old
) p
))
856 (cons (cons (car (car ulist
)) p
) math-fbu-base
))))
857 (setq ulist
(cdr ulist
)))))
858 ((math-scalarp expr
))
859 ((and (eq (car expr
) '^
)
860 (integerp (nth 2 expr
)))
861 (math-find-base-units-rec (nth 1 expr
) (* pow
(nth 2 expr
))))
863 (math-find-base-units-rec (nth 1 expr
) pow
)
864 (math-find-base-units-rec (nth 2 expr
) pow
))
866 (math-find-base-units-rec (nth 1 expr
) pow
)
867 (math-find-base-units-rec (nth 2 expr
) (- pow
)))
868 ((eq (car expr
) 'neg
)
869 (math-find-base-units-rec (nth 1 expr
) pow
))
871 (math-find-base-units-rec (nth 1 expr
) pow
))
872 ((eq (car expr
) 'var
)
873 (or (eq (nth 1 expr
) 'pi
)
874 (error "Unknown name %s in defining expression for unit %s"
875 (nth 1 expr
) (car math-fbu-entry
))))
876 ((equal expr
'(calcFunc-ln 10)))
877 (t (error "Malformed defining expression for unit %s" (car math-fbu-entry
))))))
880 (defun math-units-in-expr-p (expr sub-exprs
)
882 (if (eq (car expr
) 'var
)
883 (math-check-unit-name expr
)
885 (memq (car expr
) '(* / ^
)))
886 (or (math-units-in-expr-p (nth 1 expr
) sub-exprs
)
887 (math-units-in-expr-p (nth 2 expr
) sub-exprs
))))))
889 (defun math-only-units-in-expr-p (expr)
891 (if (eq (car expr
) 'var
)
892 (math-check-unit-name expr
)
893 (if (memq (car expr
) '(* /))
894 (and (math-only-units-in-expr-p (nth 1 expr
))
895 (math-only-units-in-expr-p (nth 2 expr
)))
896 (and (eq (car expr
) '^
)
897 (and (math-only-units-in-expr-p (nth 1 expr
))
898 (math-realp (nth 2 expr
))))))))
900 (defun math-single-units-in-expr-p (expr)
901 (cond ((math-scalarp expr
) nil
)
902 ((eq (car expr
) 'var
)
903 (math-check-unit-name expr
))
905 (let ((u1 (math-single-units-in-expr-p (nth 1 expr
)))
906 (u2 (math-single-units-in-expr-p (nth 2 expr
))))
907 (or (and u1 u2
'wrong
)
911 (if (math-units-in-expr-p (nth 2 expr
) nil
)
913 (math-single-units-in-expr-p (nth 1 expr
))))
916 (defun math-check-unit-name (v)
917 (and (eq (car-safe v
) 'var
)
918 (or (assq (nth 1 v
) (or math-units-table
(math-build-units-table)))
919 (let ((name (symbol-name (nth 1 v
))))
920 (and (> (length name
) 1)
921 (assq (aref name
0) math-unit-prefixes
)
922 (or (assq (intern (substring name
1)) math-units-table
)
923 (and (eq (aref name
0) ?M
)
925 (eq (aref name
1) ?e
)
926 (eq (aref name
2) ?g
)
927 (assq (intern (substring name
3))
928 math-units-table
))))))))
930 ;; The variable math-which-standard is local to math-to-standard-units,
931 ;; but is used by math-to-standard-rec, which is called by
932 ;; math-to-standard-units.
933 (defvar math-which-standard
)
935 (defun math-to-standard-units (expr math-which-standard
)
936 (math-to-standard-rec expr
))
938 (defun math-to-standard-rec (expr)
939 (if (eq (car-safe expr
) 'var
)
940 (let ((u (math-check-unit-name expr
))
945 (setq expr
(math-to-standard-rec (nth 1 u
)))
946 (let ((st (assq (car u
) math-which-standard
)))
948 (setq expr
(nth 1 st
))
949 (setq expr
(list 'var
(car u
)
950 (intern (concat "var-"
956 (nth 1 (assq (aref (symbol-name base
) 0)
963 (if (Math-primp expr
)
966 (mapcar 'math-to-standard-rec
(cdr expr
))))))
968 (defun math-apply-units (expr units ulist
&optional pure
)
969 (setq expr
(math-simplify-units expr
))
973 (or (math-numberp expr
)
974 (error "Incompatible units"))
976 (setq value
(math-div expr
(nth 1 (car ulist
)))
977 value
(math-floor (let ((calc-internal-prec
978 (1- calc-internal-prec
)))
979 (math-normalize value
)))
980 new
(math-add new
(math-mul value
(car (car ulist
))))
981 expr
(math-sub expr
(math-mul value
(nth 1 (car ulist
))))
983 (math-add new
(math-mul (math-div expr
(nth 1 (car ulist
)))
987 (math-simplify-units (list '* expr units
)))))
989 (defvar math-decompose-units-cache nil
)
990 (defun math-decompose-units (units)
991 (let ((u (math-check-unit-name units
)))
992 (and u
(eq (car-safe (nth 1 u
)) '+)
993 (setq units
(nth 1 u
))))
994 (setq units
(calcFunc-expand units
))
995 (and (eq (car-safe units
) '+)
996 (let ((entry (list units calc-internal-prec calc-prefer-frac
)))
997 (or (equal entry
(car math-decompose-units-cache
))
1001 (while (eq (car-safe utemp
) '+)
1002 (setq ulist
(cons (math-decompose-unit-part (nth 2 utemp
))
1004 utemp
(nth 1 utemp
)))
1005 (setq ulist
(cons (math-decompose-unit-part utemp
) ulist
)
1007 (while (setq utemp
(cdr utemp
))
1008 (unless (equal (nth 2 (car utemp
)) (nth 2 (car ulist
)))
1009 (error "Inconsistent units in sum")))
1010 (setq math-decompose-units-cache
1015 (not (Math-lessp (nth 1 x
)
1017 (cdr math-decompose-units-cache
))))
1019 (defun math-decompose-unit-part (unit)
1021 (math-is-multiple (math-simplify-units (math-to-standard-units
1025 ;; The variable math-fcu-u is local to math-find-compatible-unit,
1026 ;; but is used by math-find-compatible-rec which is called by
1027 ;; math-find-compatible-unit.
1030 (defun math-find-compatible-unit (expr unit
)
1031 (let ((math-fcu-u (math-check-unit-name unit
)))
1033 (math-find-compatible-unit-rec expr
1))))
1035 (defun math-find-compatible-unit-rec (expr pow
)
1036 (cond ((eq (car-safe expr
) '*)
1037 (or (math-find-compatible-unit-rec (nth 1 expr
) pow
)
1038 (math-find-compatible-unit-rec (nth 2 expr
) pow
)))
1039 ((eq (car-safe expr
) '/)
1040 (or (math-find-compatible-unit-rec (nth 1 expr
) pow
)
1041 (math-find-compatible-unit-rec (nth 2 expr
) (- pow
))))
1042 ((and (eq (car-safe expr
) '^
)
1043 (integerp (nth 2 expr
)))
1044 (math-find-compatible-unit-rec (nth 1 expr
) (* pow
(nth 2 expr
))))
1046 (let ((u2 (math-check-unit-name expr
)))
1047 (if (equal (nth 4 math-fcu-u
) (nth 4 u2
))
1048 (cons expr pow
))))))
1050 ;; The variables math-cu-new-units and math-cu-pure are local to
1051 ;; math-convert-units, but are used by math-convert-units-rec,
1052 ;; which is called by math-convert-units.
1053 (defvar math-cu-new-units
)
1054 (defvar math-cu-pure
)
1056 (defun math-convert-units (expr math-cu-new-units
&optional math-cu-pure
)
1057 (if (eq (car-safe math-cu-new-units
) 'var
)
1058 (let ((unew (assq (nth 1 math-cu-new-units
)
1059 (math-build-units-table))))
1060 (if (eq (car-safe (nth 1 unew
)) '+)
1061 (setq math-cu-new-units
(nth 1 unew
)))))
1062 (math-with-extra-prec 2
1063 (let ((compat (and (not math-cu-pure
)
1064 (math-find-compatible-unit expr math-cu-new-units
)))
1065 (math-cu-unit-list nil
)
1066 (math-combining-units nil
))
1068 (math-simplify-units
1069 (math-mul (math-mul (math-simplify-units
1070 (math-div expr
(math-pow (car compat
)
1072 (math-pow math-cu-new-units
(cdr compat
)))
1073 (math-simplify-units
1074 (math-to-standard-units
1075 (math-pow (math-div (car compat
) math-cu-new-units
)
1078 (when (setq math-cu-unit-list
(math-decompose-units math-cu-new-units
))
1079 (setq math-cu-new-units
(nth 2 (car math-cu-unit-list
))))
1080 (when (eq (car-safe expr
) '+)
1081 (setq expr
(math-simplify-units expr
)))
1082 (if (math-units-in-expr-p expr t
)
1083 (math-convert-units-rec expr
)
1084 (math-apply-units (math-to-standard-units
1085 (list '/ expr math-cu-new-units
) nil
)
1086 math-cu-new-units math-cu-unit-list math-cu-pure
))))))
1088 (defun math-convert-units-rec (expr)
1089 (if (math-units-in-expr-p expr nil
)
1090 (math-apply-units (math-to-standard-units
1091 (list '/ expr math-cu-new-units
) nil
)
1092 math-cu-new-units math-cu-unit-list math-cu-pure
)
1093 (if (Math-primp expr
)
1096 (mapcar 'math-convert-units-rec
(cdr expr
))))))
1098 (defun math-convert-temperature (expr old new
&optional pure
)
1099 (let* ((units (math-single-units-in-expr-p expr
))
1101 (if (or (null units
)
1102 (equal (nth 1 old
) (car units
)))
1103 (math-check-unit-name old
)
1104 (error "Inconsistent temperature units"))
1106 (unew (math-check-unit-name new
)))
1107 (unless (and (consp unew
) (nth 3 unew
))
1108 (error "Not a valid temperature unit"))
1109 (unless (and (consp uold
) (nth 3 uold
))
1110 (error "Not a pure temperature expression"))
1111 (let ((v (car uold
)))
1112 (setq expr
(list '/ expr
(list 'var v
1113 (intern (concat "var-"
1114 (symbol-name v
)))))))
1115 (or (eq (nth 3 uold
) (nth 3 unew
))
1116 (cond ((eq (nth 3 uold
) 'K
)
1117 (setq expr
(list '- expr
'(/ 27315 100)))
1118 (if (eq (nth 3 unew
) 'F
)
1119 (setq expr
(list '+ (list '* expr
'(/ 9 5)) 32))))
1120 ((eq (nth 3 uold
) 'C
)
1121 (if (eq (nth 3 unew
) 'F
)
1122 (setq expr
(list '+ (list '* expr
'(/ 9 5)) 32))
1123 (setq expr
(list '+ expr
'(/ 27315 100)))))
1125 (setq expr
(list '* (list '- expr
32) '(/ 5 9)))
1126 (if (eq (nth 3 unew
) 'K
)
1127 (setq expr
(list '+ expr
'(/ 27315 100)))))))
1130 (list '* expr new
))))
1134 (defun math-simplify-units (a)
1135 (let ((math-simplifying-units t
)
1136 (calc-matrix-mode 'scalar
))
1138 (defalias 'calcFunc-usimplify
'math-simplify-units
)
1140 ;; The function created by math-defsimplify uses the variable
1141 ;; math-simplify-expr, and so is used by functions in math-defsimplify
1142 (defvar math-simplify-expr
)
1144 (math-defsimplify (+ -
)
1145 (and math-simplifying-units
1146 (math-units-in-expr-p (nth 1 math-simplify-expr
) nil
)
1147 (let* ((units (math-extract-units (nth 1 math-simplify-expr
)))
1148 (ratio (math-simplify (math-to-standard-units
1149 (list '/ (nth 2 math-simplify-expr
) units
) nil
))))
1150 (if (math-units-in-expr-p ratio nil
)
1152 (calc-record-why "*Inconsistent units" math-simplify-expr
)
1154 (list '* (math-add (math-remove-units (nth 1 math-simplify-expr
))
1155 (if (eq (car math-simplify-expr
) '-
)
1156 (math-neg ratio
) ratio
))
1160 (math-simplify-units-prod))
1162 (defun math-simplify-units-prod ()
1163 (and math-simplifying-units
1164 calc-autorange-units
1165 (Math-realp (nth 1 math-simplify-expr
))
1166 (let* ((num (math-float (nth 1 math-simplify-expr
)))
1167 (xpon (calcFunc-xpon num
))
1168 (unitp (cdr (cdr math-simplify-expr
)))
1170 (pow (if (eq (car math-simplify-expr
) '*) 1 -
1))
1172 (and (eq (car-safe unit
) '*)
1173 (setq unitp
(cdr unit
)
1175 (and (eq (car-safe unit
) '^
)
1176 (integerp (nth 2 unit
))
1177 (setq pow
(* pow
(nth 2 unit
))
1180 (and (setq u
(math-check-unit-name unit
))
1183 (>= xpon
(if (eq (car u
) 'm
) 1 3)))
1186 (reverse math-unit-prefixes
)
1187 math-unit-prefixes
))
1190 (or (eq (car u
) (nth 1 unit
))
1195 math-unit-prefixes
))))))
1196 (setq xpon
(+ xpon uxpon
))
1198 (or (memq (car (car p
)) '(?d ?D ?h ?H
))
1199 (and (eq (car (car p
)) ?c
)
1200 (not (eq (car u
) 'm
)))
1201 (< xpon
(setq pxpon
(* (nth 2 (nth 1 (car p
)))
1204 (setq pname
(math-build-var-name
1205 (if (eq (car (car p
)) 0)
1207 (concat (char-to-string
1211 (and (/= (car (car p
)) 0)
1213 math-units-table
)))))
1217 (or (not (eq p pref
))
1218 (< xpon
(+ pxpon
(* (math-abs pow
) 3))))
1220 (setcar (cdr math-simplify-expr
)
1221 (let ((calc-prefer-frac nil
))
1222 (calcFunc-scf (nth 1 math-simplify-expr
)
1224 (setcar unitp pname
)
1225 math-simplify-expr
)))))))
1227 (defvar math-try-cancel-units
)
1230 (and math-simplifying-units
1231 (let ((np (cdr math-simplify-expr
))
1232 (math-try-cancel-units 0)
1234 (setq n
(if (eq (car-safe (nth 2 math-simplify-expr
)) '*)
1235 (cdr (nth 2 math-simplify-expr
))
1236 (nthcdr 2 math-simplify-expr
)))
1237 (if (math-realp (car n
))
1239 (setcar (cdr math-simplify-expr
) (math-mul (nth 1 math-simplify-expr
)
1240 (let ((calc-prefer-frac nil
))
1241 (math-div 1 (car n
)))))
1243 (while (eq (car-safe (setq n
(car np
))) '*)
1244 (math-simplify-units-divisor (cdr n
) (cdr (cdr math-simplify-expr
)))
1245 (setq np
(cdr (cdr n
))))
1246 (math-simplify-units-divisor np
(cdr (cdr math-simplify-expr
)))
1247 (if (eq math-try-cancel-units
0)
1248 (let* ((math-simplifying-units nil
)
1249 (base (math-simplify
1250 (math-to-standard-units math-simplify-expr nil
))))
1251 (if (Math-numberp base
)
1252 (setq math-simplify-expr base
))))
1253 (if (eq (car-safe math-simplify-expr
) '/)
1254 (math-simplify-units-prod))
1255 math-simplify-expr
)))
1257 (defun math-simplify-units-divisor (np dp
)
1260 (while (eq (car-safe (setq d
(car dp
))) '*)
1261 (when (setq temp
(math-simplify-units-quotient n
(nth 1 d
)))
1262 (setcar np
(setq n temp
))
1264 (setq dp
(cdr (cdr d
))))
1265 (when (setq temp
(math-simplify-units-quotient n d
))
1266 (setcar np
(setq n temp
))
1269 ;; Simplify, e.g., "in / cm" to "2.54" in a units expression.
1270 (defun math-simplify-units-quotient (n d
)
1273 (when (and (eq (car-safe n
) '^
)
1274 (integerp (nth 2 n
)))
1275 (setq pow1
(nth 2 n
) n
(nth 1 n
)))
1276 (when (and (eq (car-safe d
) '^
)
1277 (integerp (nth 2 d
)))
1278 (setq pow2
(nth 2 d
) d
(nth 1 d
)))
1279 (let ((un (math-check-unit-name n
))
1280 (ud (math-check-unit-name d
)))
1282 (if (and (equal (nth 4 un
) (nth 4 ud
))
1285 (math-to-standard-units (list '/ n d
) nil
)
1286 (list '^
(math-to-standard-units (list '/ n d
) nil
) pow1
))
1293 (and (eq (car (car un
)) (car (car ud1
)))
1294 (setq math-try-cancel-units
1295 (+ math-try-cancel-units
1296 (- (* (cdr (car un
)) pow1
)
1297 (* (cdr (car ud
)) pow2
)))))
1298 (setq ud1
(cdr ud1
)))
1303 (and math-simplifying-units
1304 (math-realp (nth 2 math-simplify-expr
))
1305 (if (memq (car-safe (nth 1 math-simplify-expr
)) '(* /))
1306 (list (car (nth 1 math-simplify-expr
))
1307 (list '^
(nth 1 (nth 1 math-simplify-expr
))
1308 (nth 2 math-simplify-expr
))
1309 (list '^
(nth 2 (nth 1 math-simplify-expr
))
1310 (nth 2 math-simplify-expr
)))
1311 (math-simplify-units-pow (nth 1 math-simplify-expr
)
1312 (nth 2 math-simplify-expr
)))))
1314 (math-defsimplify calcFunc-sqrt
1315 (and math-simplifying-units
1316 (if (memq (car-safe (nth 1 math-simplify-expr
)) '(* /))
1317 (list (car (nth 1 math-simplify-expr
))
1318 (list 'calcFunc-sqrt
(nth 1 (nth 1 math-simplify-expr
)))
1319 (list 'calcFunc-sqrt
(nth 2 (nth 1 math-simplify-expr
))))
1320 (math-simplify-units-pow (nth 1 math-simplify-expr
) '(frac 1 2)))))
1322 (math-defsimplify (calcFunc-floor
1332 (and math-simplifying-units
1333 (= (length math-simplify-expr
) 2)
1334 (if (math-only-units-in-expr-p (nth 1 math-simplify-expr
))
1335 (nth 1 math-simplify-expr
)
1336 (if (and (memq (car-safe (nth 1 math-simplify-expr
)) '(* /))
1337 (or (math-only-units-in-expr-p
1338 (nth 1 (nth 1 math-simplify-expr
)))
1339 (math-only-units-in-expr-p
1340 (nth 2 (nth 1 math-simplify-expr
)))))
1341 (list (car (nth 1 math-simplify-expr
))
1342 (cons (car math-simplify-expr
)
1343 (cons (nth 1 (nth 1 math-simplify-expr
))
1344 (cdr (cdr math-simplify-expr
))))
1345 (cons (car math-simplify-expr
)
1346 (cons (nth 2 (nth 1 math-simplify-expr
))
1347 (cdr (cdr math-simplify-expr
)))))))))
1349 (defun math-simplify-units-pow (a pow
)
1350 (if (and (eq (car-safe a
) '^
)
1351 (math-check-unit-name (nth 1 a
))
1352 (math-realp (nth 2 a
)))
1353 (list '^
(nth 1 a
) (math-mul pow
(nth 2 a
)))
1354 (let* ((u (math-check-unit-name a
))
1355 (pf (math-to-simple-fraction pow
))
1356 (d (and (eq (car-safe pf
) 'frac
) (nth 2 pf
))))
1358 (math-units-are-multiple u d
)
1359 (list '^
(math-to-standard-units a nil
) pow
)))))
1362 (defun math-units-are-multiple (u n
)
1364 (while (and u
(= (%
(cdr (car u
)) n
) 0))
1368 (math-defsimplify calcFunc-sin
1369 (and math-simplifying-units
1370 (math-units-in-expr-p (nth 1 math-simplify-expr
) nil
)
1371 (let ((rad (math-simplify-units
1373 (math-to-standard-units (nth 1 math-simplify-expr
) nil
))))
1374 (calc-angle-mode 'rad
))
1375 (and (eq (car-safe rad
) '*)
1376 (math-realp (nth 1 rad
))
1377 (eq (car-safe (nth 2 rad
)) 'var
)
1378 (eq (nth 1 (nth 2 rad
)) 'rad
)
1379 (list 'calcFunc-sin
(nth 1 rad
))))))
1381 (math-defsimplify calcFunc-cos
1382 (and math-simplifying-units
1383 (math-units-in-expr-p (nth 1 math-simplify-expr
) nil
)
1384 (let ((rad (math-simplify-units
1386 (math-to-standard-units (nth 1 math-simplify-expr
) nil
))))
1387 (calc-angle-mode 'rad
))
1388 (and (eq (car-safe rad
) '*)
1389 (math-realp (nth 1 rad
))
1390 (eq (car-safe (nth 2 rad
)) 'var
)
1391 (eq (nth 1 (nth 2 rad
)) 'rad
)
1392 (list 'calcFunc-cos
(nth 1 rad
))))))
1394 (math-defsimplify calcFunc-tan
1395 (and math-simplifying-units
1396 (math-units-in-expr-p (nth 1 math-simplify-expr
) nil
)
1397 (let ((rad (math-simplify-units
1399 (math-to-standard-units (nth 1 math-simplify-expr
) nil
))))
1400 (calc-angle-mode 'rad
))
1401 (and (eq (car-safe rad
) '*)
1402 (math-realp (nth 1 rad
))
1403 (eq (car-safe (nth 2 rad
)) 'var
)
1404 (eq (nth 1 (nth 2 rad
)) 'rad
)
1405 (list 'calcFunc-tan
(nth 1 rad
))))))
1407 (math-defsimplify calcFunc-sec
1408 (and math-simplifying-units
1409 (math-units-in-expr-p (nth 1 math-simplify-expr
) nil
)
1410 (let ((rad (math-simplify-units
1412 (math-to-standard-units (nth 1 math-simplify-expr
) nil
))))
1413 (calc-angle-mode 'rad
))
1414 (and (eq (car-safe rad
) '*)
1415 (math-realp (nth 1 rad
))
1416 (eq (car-safe (nth 2 rad
)) 'var
)
1417 (eq (nth 1 (nth 2 rad
)) 'rad
)
1418 (list 'calcFunc-sec
(nth 1 rad
))))))
1420 (math-defsimplify calcFunc-csc
1421 (and math-simplifying-units
1422 (math-units-in-expr-p (nth 1 math-simplify-expr
) nil
)
1423 (let ((rad (math-simplify-units
1425 (math-to-standard-units (nth 1 math-simplify-expr
) nil
))))
1426 (calc-angle-mode 'rad
))
1427 (and (eq (car-safe rad
) '*)
1428 (math-realp (nth 1 rad
))
1429 (eq (car-safe (nth 2 rad
)) 'var
)
1430 (eq (nth 1 (nth 2 rad
)) 'rad
)
1431 (list 'calcFunc-csc
(nth 1 rad
))))))
1433 (math-defsimplify calcFunc-cot
1434 (and math-simplifying-units
1435 (math-units-in-expr-p (nth 1 math-simplify-expr
) nil
)
1436 (let ((rad (math-simplify-units
1438 (math-to-standard-units (nth 1 math-simplify-expr
) nil
))))
1439 (calc-angle-mode 'rad
))
1440 (and (eq (car-safe rad
) '*)
1441 (math-realp (nth 1 rad
))
1442 (eq (car-safe (nth 2 rad
)) 'var
)
1443 (eq (nth 1 (nth 2 rad
)) 'rad
)
1444 (list 'calcFunc-cot
(nth 1 rad
))))))
1447 (defun math-remove-units (expr)
1448 (if (math-check-unit-name expr
)
1450 (if (Math-primp expr
)
1453 (mapcar 'math-remove-units
(cdr expr
))))))
1455 (defun math-extract-units (expr)
1456 (if (memq (car-safe expr
) '(* /))
1458 (mapcar 'math-extract-units
(cdr expr
)))
1459 (if (math-check-unit-name expr
) expr
1)))
1461 (defun math-build-units-table-buffer (enter-buffer)
1462 (if (not (and math-units-table math-units-table-buffer-valid
1463 (get-buffer "*Units Table*")))
1464 (let ((buf (get-buffer-create "*Units Table*"))
1465 (uptr (math-build-units-table))
1466 (calc-language (if (eq calc-language
'big
) nil calc-language
))
1467 (calc-float-format '(float 0))
1468 (calc-group-digits nil
)
1469 (calc-number-radix 10)
1470 (calc-twos-complement-mode nil
)
1471 (calc-point-char ".")
1475 (message "Formatting units table...")
1477 (let ((inhibit-read-only t
))
1479 (insert "Calculator Units Table:\n\n")
1480 (insert "(All definitions are exact unless marked with an asterisk (*).)\n\n")
1481 (insert "Unit Type Definition Description\n\n")
1485 (when (eq (car u
) 'm
)
1487 (setq shadowed
(and std
(assq (car u
) math-additional-units
)))
1490 (eq (aref name
0) ?\
*))
1491 (unless (eq uptr math-units-table
)
1493 (setq name
(substring name
1)))
1495 (and shadowed
(insert "("))
1496 (insert (symbol-name (car u
)))
1497 (and shadowed
(insert ")"))
1501 (insert (symbol-name (nth 3 u
))))
1507 (and shadowed
(insert "("))
1511 (insert (math-format-value (nth 1 u
) 80))
1512 (insert (symbol-name (car u
)))))
1513 (and shadowed
(insert ")"))
1519 (insert " (redefined above)")
1521 (insert " (base unit)")))
1523 (setq uptr
(cdr uptr
)))
1524 (insert "\n\nUnit Prefix Table:\n\n")
1525 (setq uptr math-unit-prefixes
)
1528 (insert " " (char-to-string (car u
)))
1529 (if (equal (nth 1 u
) (nth 1 (nth 1 uptr
)))
1530 (insert " " (char-to-string (car (car (setq uptr
(cdr uptr
)))))
1533 (insert "10^" (int-to-string (nth 2 (nth 1 u
))))
1535 (insert " " (nth 2 u
) "\n")
1536 (while (eq (car (car (setq uptr
(cdr uptr
)))) 0)))
1538 (insert "(**) When in TeX or LaTeX display mode, the TeX specific unit\n"
1539 "names will not use the `tex' prefix; the unit name for a\n"
1540 "TeX point will be `pt' instead of `texpt', for example.\n"
1541 "To avoid conflicts, the unit names for pint and parsec will\n"
1542 "be `pint' and `parsec' instead of `pt' and `pc'."))
1544 (message "Formatting units table...done"))
1545 (setq math-units-table-buffer-valid t
)
1546 (let ((oldbuf (current-buffer)))
1548 (goto-char (point-min))
1549 (set-buffer oldbuf
))
1552 (display-buffer buf
)))
1554 (pop-to-buffer (get-buffer "*Units Table*"))
1555 (display-buffer (get-buffer "*Units Table*")))))
1557 ;;; Logarithmic units functions
1559 (defvar math-logunits
'((var dB var-dB
)
1562 (defun math-extract-logunits (expr)
1563 (if (memq (car-safe expr
) '(* /))
1565 (mapcar 'math-extract-logunits
(cdr expr
)))
1566 (if (memq (car-safe expr
) '(^
))
1567 (list '^
(math-extract-logunits (nth 1 expr
)) (nth 2 expr
))
1568 (if (member expr math-logunits
) expr
1))))
1570 (defun math-logcombine (a b neg
)
1571 (let ((aunit (math-simplify (math-extract-logunits a
))))
1572 (if (not (eq (car-safe aunit
) 'var
))
1573 (calc-record-why "*Improper logarithmic unit" aunit
)
1574 (let* ((units (math-extract-units a
))
1575 (acoeff (math-simplify (math-remove-units a
)))
1576 (bcoeff (math-simplify (math-to-standard-units
1577 (list '/ b units
) nil
))))
1578 (if (math-units-in-expr-p bcoeff nil
)
1579 (calc-record-why "*Inconsistent units" nil
)
1581 (or (math-lessp acoeff bcoeff
)
1582 (math-equal acoeff bcoeff
)))
1583 (calc-record-why "*Improper coefficients" nil
)
1585 (if (equal aunit
'(var dB var-dB
))
1590 (math-pow 10 (math-div acoeff
10))
1591 (math-pow 10 (math-div bcoeff
10)))
1593 (math-pow 10 (math-div acoeff
10))
1594 (math-pow 10 (math-div bcoeff
10))))))
1598 (calcFunc-exp acoeff
)
1599 (calcFunc-exp bcoeff
))
1601 (calcFunc-exp acoeff
)
1602 (calcFunc-exp bcoeff
)))))
1605 (defun calcFunc-luplus (a b
)
1606 (math-logcombine a b nil
))
1608 (defun calcFunc-luminus (a b
)
1609 (math-logcombine a b t
))
1611 (defun calc-luplus (arg)
1614 (if (calc-is-inverse)
1615 (calc-binary-op "lu-" 'calcFunc-luminus arg
)
1616 (calc-binary-op "lu+" 'calcFunc-luplus arg
))))
1618 (defun calc-luminus (arg)
1621 (if (calc-is-inverse)
1622 (calc-binary-op "lu+" 'calcFunc-luplus arg
)
1623 (calc-binary-op "lu-" 'calcFunc-luminus arg
))))
1625 ;(defun calcFunc-lmul (a b)
1628 (defun math-logunit-level (val ref power
)
1629 (let ((lunit (math-simplify (math-extract-logunits val
))))
1630 (if (not (eq (car-safe lunit
) 'var
))
1631 (calc-record-why "*Improper logarithmic unit" lunit
)
1632 (if (not (eq 1 (math-simplify (math-extract-units (math-div val lunit
)))))
1633 (calc-record-why "*Inappropriate units" nil
)
1634 (let ((coeff (math-simplify (math-div val lunit
))))
1635 (if (equal lunit
'(var dB var-dB
))
1650 (defvar calc-default-field-reference-level
)
1651 (defvar calc-default-power-reference-level
)
1653 (defun calcFunc-fieldlevel (val &optional ref
)
1655 (setq ref
(math-read-expr calc-default-field-reference-level
)))
1656 (math-logunit-level val ref nil
))
1658 (defun calcFunc-powerlevel (val &optional ref
)
1660 (setq ref
(math-read-expr calc-default-power-reference-level
)))
1661 (math-logunit-level val ref t
))
1663 (defun calc-level (arg)
1666 (if (calc-is-hyperbolic)
1667 (if (calc-is-option)
1668 (calc-binary-op "plvl" 'calcFunc-powerlevel arg
)
1669 (calc-unary-op "plvl" 'calcFunc-powerlevel arg
))
1670 (if (calc-is-option)
1671 (calc-binary-op "flvl" 'calcFunc-fieldlevel arg
)
1672 (calc-unary-op "flvl" 'calcFunc-fieldlevel arg
)))))
1674 (provide 'calc-units
)
1680 ;;; calc-units.el ends here