Commit 405df618a38d3b8ddaae691f865bbf068e931ac5

Authored by Raymond Toy
1 parent c388f817

Add Euler's constant and fix dependencies.

oct.asd:
o qd-const2 depends on qd-const.

qd-const.lisp:
o Add qd value for Euler's constant

qd-const2.lisp:
o Add qd-real value for Euler's constant.
... ... @@ -52,7 +52,7 @@
52 52 :depends-on ("qd" "qd-const"))
53 53 (:file "qd-class"
54 54 :depends-on ("qd-fun"))
55   - (:file "qd-const2" :depends-on ("qd-class"))
  55 + (:file "qd-const2" :depends-on ("qd-class" "qd-const"))
56 56 (:file "qd-methods"
57 57 :depends-on ("qd-class"))
58 58 (:file "qd-reader"
... ...
qd-const.lisp
... ... @@ -3752,6 +3752,15 @@
3752 3752 (%make-qd-d q0 q1 q2 q3))))
3753 3753 "A table of cos(k/1024) for k = 1 to 256")
3754 3754  
  3755 +(defconstant +qd-%gamma+
  3756 + (multiple-value-bind (q0 q1 q2 q3)
  3757 + (renorm-5 (scale-float (float 5199096506725913 1.0d0) -53)
  3758 + (scale-float (float -6416270704492459 1.0d0) -110)
  3759 + (scale-float (float -5430032177773649 1.0d0) -164)
  3760 + (scale-float (float 7163319935831933 1.0d0) -218)
  3761 + (scale-float (float -6507777405222912 1.0d0) -273))
  3762 + (%make-qd-d q0 q1 q2 q3))
  3763 + "Quad-double value of Euler's constant")
3755 3764  
3756 3765 #||
3757 3766  
... ...
qd-const2.lisp
... ... @@ -82,3 +82,7 @@
82 82 (defconstant +qd-real-one+
83 83 (make-instance 'qd-real :value (make-qd-d 1d0))
84 84 "QD-REAL representation of 1")
  85 +
  86 +(defconstant +%gamma+
  87 + (make-instance 'qd-real :value octi::+qd-%gamma+)
  88 + "Euler's constant")
85 89 \ No newline at end of file
... ...