share/contrib/floatproperties.lisp

   1 ;; Expose some properties of double floating point numbers to Maxima.
   2 ;; Note: floatbits is one plus the number of bits in the fractional part.
   3
   4 (defun $float_eps ()
   5   flonum-epsilon)
   6
   7 (defmvar $largest_float most-positive-flonum)
   8 (setf (get '$largest_float 'assign) 'neverset)
   9
  10 (defmvar $least_positive_float least-positive-flonum)
  11 (setf (get '$least_positive_float 'assign) 'neverset)
  12
  13 (defun $float_bits ()
  14   (float-digits 0.0))
  15
  16 (defun $bigfloat_eps ()
  17   (let ((r ($bfloat (div 1 (expt 2 fpprec)))))
  18     (list (first r) (incf (second r)) (third r))))