Update of /project/oct/cvsroot/oct In directory clnet:/tmp/cvs-serv20277 Modified Files: qd-const.lisp qd-io.lisp qd-methods.lisp Log Message: qd-io.lisp: o Get rid of the #q reader that returns an internal %qd-real. qd-methods.lisp: o QPHASE was returning the wrong kind of object for positive values. o Replace use of #q0 with the construction of a qd-real 0. qd-const.lisp: o Replace all uses of #q with #.(qd-from-string ...). --- /project/oct/cvsroot/oct/qd-const.lisp 2007/05/31 02:47:37 1.13 +++ /project/oct/cvsroot/oct/qd-const.lisp 2007/08/24 21:45:16 1.14 @@ -2201,69 +2201,69 @@ +qd-pi/4+ ;; Do we need to make these values more accurate? (The ;; reader has quite a bit of roundoff.) - #q0.46364760900080611621425623146121440202853705428612026381093308872018q0 - #q0.24497866312686415417208248121127581091414409838118406712737591466738q0 - #q0.12435499454676143503135484916387102557317019176980408991511411911572q0 - #q0.062418809995957348473979112985505113606273887797499194607527816898697q0 - #q0.031239833430268276253711744892490977032495663725400040255315586255793q0 - #q0.0156237286204768308028015212565703189111141398009054178814105073966645q0 - #q0.0078123410601011112964633918421992816212228117250147235574539022483893q0 - #q0.003906230131966971827628665311424387140357490115202856215213095149011q0 - #q0.00195312251647881868512148262507671393161074677723351033905753396043094q0 - #q9.7656218955931943040343019971729085163419701581008759004900725226773q-4 - #q4.8828121119489827546923962564484866619236113313500303710940335348752q-4 - #q2.4414062014936176401672294325965998621241779097061761180790046091019q-4 - #q1.22070311893670204239058646117956300930829409015787498451939837846645q-4 - #q6.1035156174208775021662569173829153785143536833346179337671134316588q-5 - #q3.0517578115526096861825953438536019750949675119437837531021156883611q-5 - #q1.5258789061315762107231935812697885137429238144575874846241186407446q-5 - #q7.6293945311019702633884823401050905863507439184680771577638306965336q-6 - #q3.8146972656064962829230756163729937228052573039688663101874392503939q-6 - #q1.9073486328101870353653693059172441687143421654501533666700577234671q-6 - #q9.53674316405960879420670689923112390019634124498790160133611802076q-7 - #q4.7683715820308885992758382144924707587049404378664196740053215887142q-7 - #q2.3841857910155798249094797721893269783096898769063155913766911372218q-7 - #q1.19209289550780685311368497137922112645967587664586735576738225215437q-7 - #q5.9604644775390554413921062141788874250030195782366297314294565710003q-8 - #q2.9802322387695303676740132767709503349043907067445107249258477840843q-8 - #q1.4901161193847655147092516595963247108248930025964720012170057805491q-8 - #q7.4505805969238279871365645744953921132066925545665870075947601416172q-9 - #q3.725290298461914045267070571811923583671948328737040524231998269239q-9 - #q1.8626451492309570290958838214764904345065282835738863513491050124951q-9 - #q9.3132257461547851535573547768456130389292649614929067394376854242196q-10 - #q4.6566128730773925777884193471057016297347863891561617421323492554414q-10 - #q2.32830643653869628902042741838821270371274293204981860525486662280605q-10 - #q1.16415321826934814452599092729852658796396457380014290026584979170883q-10 - #q5.8207660913467407226496761591231582349549156257795272423976206167147q-11 - #q2.9103830456733703613273032698903947793693632003639830495829934525029q-11 - #q1.4551915228366851806639597837362993474211703608936710732067270213307q-11 - #q7.2759576141834259033201841046703741842764629388821429640111752890838q-12 - #q3.6379788070917129516601402005837967730345578669779258118296083646486q-12 - #q1.81898940354585647583007611882297459662931973336029253714520765350336q-12 - #q9.094947017729282379150388117278718245786649666696631862264792881855q-13 - #q4.5474735088646411895751949990348397807233312083369623012466392138249q-13 - #q2.2737367544323205947875976170668549725904164010421166413578155299654q-13 - #q1.1368683772161602973937988232271068715738020501302644662229139921281q-13 - #q5.6843418860808014869689941345026335894672525626628305471702634435609q-14 - #q2.8421709430404007434844970695472041986834065703328538172835210852389q-14 - #q1.42108547152020037174224853506058802483542582129160672712566632799217q-14 - #q7.1054273576010018587112426756616725310442822766145084088962160950957q-15 - #q3.5527136788005009293556213378756778163805352845768135511116874239215q-15 - #q1.7763568394002504646778106689434441020475669105721016938889503158663q-15 - #q8.881784197001252323389053344724227002559458638215127117361184578544q-16 - #q4.440892098500626161694526672362989312819932329776890889670147968684q-16 - #q2.22044604925031308084726333618160413285249154122211136120876849284695q-16 - #q1.11022302462515654042363166809081575098156144265276392015109606150467q-16 - #q5.5511151231257827021181583404540958606019518033159549001888700768492q-17 - #q2.7755575615628913510590791702270500685127439754144943625236087596052q-17 - #q1.3877787807814456755295395851135253015328429969268117953154510949506q-17 - #q6.9388939039072283776476979255676268417598037461585147441443138686883q-18 - #q3.4694469519536141888238489627838134626418504682698143430180392335861q-18 - #q1.7347234759768070944119244813919067365411688085337267928772549041983q-18 - #q8.673617379884035472059622406959533689231148510667158491096568630248q-19 - #q4.336808689942017736029811203479766845431237313833394811387071078781q-19 - #q2.16840434497100886801490560173988342281757653922917435142338388484765q-19 - #q1.08420217248550443400745280086994171142153300490364679392792298560597q-19 + #.(qd-from-string "0.46364760900080611621425623146121440202853705428612026381093308872018q0") + #.(qd-from-string "0.24497866312686415417208248121127581091414409838118406712737591466738q0") + #.(qd-from-string "0.12435499454676143503135484916387102557317019176980408991511411911572q0") + #.(qd-from-string "0.062418809995957348473979112985505113606273887797499194607527816898697q0") + #.(qd-from-string "0.031239833430268276253711744892490977032495663725400040255315586255793q0") + #.(qd-from-string "0.0156237286204768308028015212565703189111141398009054178814105073966645q0") + #.(qd-from-string "0.0078123410601011112964633918421992816212228117250147235574539022483893q0") + #.(qd-from-string "0.003906230131966971827628665311424387140357490115202856215213095149011q0") + #.(qd-from-string "0.00195312251647881868512148262507671393161074677723351033905753396043094q0") + #.(qd-from-string "9.7656218955931943040343019971729085163419701581008759004900725226773q-4") + #.(qd-from-string "4.8828121119489827546923962564484866619236113313500303710940335348752q-4") + #.(qd-from-string "2.4414062014936176401672294325965998621241779097061761180790046091019q-4") + #.(qd-from-string "1.22070311893670204239058646117956300930829409015787498451939837846645q-4") + #.(qd-from-string "6.1035156174208775021662569173829153785143536833346179337671134316588q-5") + #.(qd-from-string "3.0517578115526096861825953438536019750949675119437837531021156883611q-5") + #.(qd-from-string "1.5258789061315762107231935812697885137429238144575874846241186407446q-5") + #.(qd-from-string "7.6293945311019702633884823401050905863507439184680771577638306965336q-6") + #.(qd-from-string "3.8146972656064962829230756163729937228052573039688663101874392503939q-6") + #.(qd-from-string "1.9073486328101870353653693059172441687143421654501533666700577234671q-6") + #.(qd-from-string "9.53674316405960879420670689923112390019634124498790160133611802076q-7") + #.(qd-from-string "4.7683715820308885992758382144924707587049404378664196740053215887142q-7") + #.(qd-from-string "2.3841857910155798249094797721893269783096898769063155913766911372218q-7") + #.(qd-from-string "1.19209289550780685311368497137922112645967587664586735576738225215437q-7") + #.(qd-from-string "5.9604644775390554413921062141788874250030195782366297314294565710003q-8") + #.(qd-from-string "2.9802322387695303676740132767709503349043907067445107249258477840843q-8") + #.(qd-from-string "1.4901161193847655147092516595963247108248930025964720012170057805491q-8") + #.(qd-from-string "7.4505805969238279871365645744953921132066925545665870075947601416172q-9") + #.(qd-from-string "3.725290298461914045267070571811923583671948328737040524231998269239q-9") + #.(qd-from-string "1.8626451492309570290958838214764904345065282835738863513491050124951q-9") + #.(qd-from-string "9.3132257461547851535573547768456130389292649614929067394376854242196q-10") + #.(qd-from-string "4.6566128730773925777884193471057016297347863891561617421323492554414q-10") + #.(qd-from-string "2.32830643653869628902042741838821270371274293204981860525486662280605q-10") + #.(qd-from-string "1.16415321826934814452599092729852658796396457380014290026584979170883q-10") + #.(qd-from-string "5.8207660913467407226496761591231582349549156257795272423976206167147q-11") + #.(qd-from-string "2.9103830456733703613273032698903947793693632003639830495829934525029q-11") + #.(qd-from-string "1.4551915228366851806639597837362993474211703608936710732067270213307q-11") + #.(qd-from-string "7.2759576141834259033201841046703741842764629388821429640111752890838q-12") + #.(qd-from-string "3.6379788070917129516601402005837967730345578669779258118296083646486q-12") + #.(qd-from-string "1.81898940354585647583007611882297459662931973336029253714520765350336q-12") + #.(qd-from-string "9.094947017729282379150388117278718245786649666696631862264792881855q-13") + #.(qd-from-string "4.5474735088646411895751949990348397807233312083369623012466392138249q-13") + #.(qd-from-string "2.2737367544323205947875976170668549725904164010421166413578155299654q-13") + #.(qd-from-string "1.1368683772161602973937988232271068715738020501302644662229139921281q-13") + #.(qd-from-string "5.6843418860808014869689941345026335894672525626628305471702634435609q-14") + #.(qd-from-string "2.8421709430404007434844970695472041986834065703328538172835210852389q-14") + #.(qd-from-string "1.42108547152020037174224853506058802483542582129160672712566632799217q-14") + #.(qd-from-string "7.1054273576010018587112426756616725310442822766145084088962160950957q-15") + #.(qd-from-string "3.5527136788005009293556213378756778163805352845768135511116874239215q-15") + #.(qd-from-string "1.7763568394002504646778106689434441020475669105721016938889503158663q-15") + #.(qd-from-string "8.881784197001252323389053344724227002559458638215127117361184578544q-16") + #.(qd-from-string "4.440892098500626161694526672362989312819932329776890889670147968684q-16") + #.(qd-from-string "2.22044604925031308084726333618160413285249154122211136120876849284695q-16") + #.(qd-from-string "1.11022302462515654042363166809081575098156144265276392015109606150467q-16") + #.(qd-from-string "5.5511151231257827021181583404540958606019518033159549001888700768492q-17") + #.(qd-from-string "2.7755575615628913510590791702270500685127439754144943625236087596052q-17") + #.(qd-from-string "1.3877787807814456755295395851135253015328429969268117953154510949506q-17") + #.(qd-from-string "6.9388939039072283776476979255676268417598037461585147441443138686883q-18") + #.(qd-from-string "3.4694469519536141888238489627838134626418504682698143430180392335861q-18") + #.(qd-from-string "1.7347234759768070944119244813919067365411688085337267928772549041983q-18") + #.(qd-from-string "8.673617379884035472059622406959533689231148510667158491096568630248q-19") + #.(qd-from-string "4.336808689942017736029811203479766845431237313833394811387071078781q-19") + #.(qd-from-string "2.16840434497100886801490560173988342281757653922917435142338388484765q-19") + #.(qd-from-string "1.08420217248550443400745280086994171142153300490364679392792298560597q-19") )) "Table of atan(2^(-k)) for k = 1 to 64. But the first three entries are 1") @@ -2553,7 +2553,7 @@ "Table of (2^(-k)) for k = -2 to 64. But the first three entries are 1") (defconstant +cordic-scale+ - #q0.065865828601599636584870082133151126045971796871364763285694473524426q0) + #.(qd-from-string "0.065865828601599636584870082133151126045971796871364763285694473524426q0"))  (defun dump-qd (qd) (flet ((dump-d (d) --- /project/oct/cvsroot/oct/qd-io.lisp 2007/06/13 16:34:52 1.11 +++ /project/oct/cvsroot/oct/qd-io.lisp 2007/08/24 21:45:16 1.12 @@ -442,5 +442,6 @@ (cl::with-input-from-string (s string) (read-qd s))) +#+nil (set-dispatch-macro-character #\# #\Q #'qd-reader) --- /project/oct/cvsroot/oct/qd-methods.lisp 2007/08/23 17:15:06 1.43 +++ /project/oct/cvsroot/oct/qd-methods.lisp 2007/08/24 21:45:16 1.44 @@ -761,7 +761,7 @@ (defmethod qphase ((x qd-real)) (if (minusp x) (- +pi+) - #q0)) + (make-instance 'qd-real :value (make-qd-d 0d0)))) (declaim (inline phase)) (defun phase (x)