h8Q#<> quad-qd.lisp by nil on RJF-E510 at 2007/1/11 14:27:52\ using 8.0 [Windows] (Oct 23, 2006 22:40).\ Optimization settings at wfasl time:\ ((safety 0) (space 0) (speed 3) (compilation-speed 0) (debug 2)) *)"8.0 [Windows] (Oct 23, 2006 22:40)Ag39744<*source-pathname* pathname quad-qdlisp unspecific  quad-qd.lisp  /UV$VFɋuÐ null 0+ @ quad-qd.lispLN<top-level-form <QAg39745@<qd internal-in-package / UV$F׋^W'ɋu null 0+ @<top-level-form <nQAg39746@<optimize speed safety space compilation-speed compiler-proclaim / UV$F^׋ɋuÐ null 0+ @top-level-form <LN<<QAg39747@<qd.dllload / UV$F^׋ɋuÐ null 0+ @<<Qqd legendre_pd and typep k quote fixnum ga<= !assert-1 !?abs = 1+  two-arg-*  two-arg-/ 1-  two-arg-- /  UVB]tE̋׋FWWWRP^׋EЉU3^W';tփ}uFV"3ɱɋuÃ}uE3V3ɱ]EڋU3۳É]ċ ;]E3؉]E3^&׋Ћ^*3W';?E̋^.׋ЋE̋^2W']̃3Ҳ^2W'^63ҲW'ЋE3ɱU]؋W';^UEډU]؋U3^2W'^:׉EE3ЋE^2W'EE3ЋE^2W'EEċ^:׉EE3ЋE^2W'ЋE^>W'Uċ^6W'@E]]؋]3ҲEċ]ċUE`E3؉]E3ЋE^2W'EE3ЋE^>W'ЋE̋^2W'EE3؉]E3ЋE^2W'^:׋ЋE^6W'QEډE3ЋE!Ë_W' kx0 2{N !i Dp3+ @c_qd_divlanguage c LNc_qd_sub!!LNc_qd_mul!!LNc_qd_add!!LN legendre_pdbf and typep k quote fixnum <= !assert-1 !into abs = 1+  two-arg-*  two-arg-/ progn structure-allocate-instance !aqd "M!.inv-structure-ref !!r ?"Ndetermine-foreign-address !< "!!!!"M!!!!r @"N!!!!"M!!!!r @"N!!< "!!!!"M!!!!r ?"N!!!!"M!!!!r @"N!!< "!!!!"M!!!!r ?"N!!< "!!!!"M!!!!r @"N!!< "!!!!"M!!!!r ?"N!!!!"M!!!!r @"N!!< "!!!!"M!!!!r @"N!!< "!!!!"M!!!!r @"N!!< "!!!!"M!!! !r ?"N!!<"!!!!"M!!!!!r @"N!!< "!!!!"M!!!"!r @"N!!< "!!!!"M!!!#!r @"N!!< "1- !!!!"M!!!%!r ?"N!!<"/' mUVtUtE׋FWWWRP^׋EЉU3^W';t֋E^׃}u333ɱɋuÃ}u 3Ҳ3ɱE3^׉E܋E^׉E3^׉E3^jEȋ]EԋEȋ us;u؏E̋E^"׋Ћ^&3W';E^*׋ЋE^.W'] 3Ҳ^.W'^23ҲW'ЋE3ɱ+Ӊ]ЋW';t]Љ]Ћ]ԋSF6@|Sv:xtp|I $tIL$pIL$śoL$ 7iy jA7x7YSu79y ta 17^>[UԋBl^B[vFhd`lI $dIL$`IL$śoL$ 7iy jA7h7YSu79y ta 17^B[VJRFN@vR\XTPXI $TIL$PIL$śoL$ 7iy jA7\7YSu79y ta 17^N[UBL^V[vZHD@LI $DIL$@IL$śoL$ 7iy jA7H7YSu79y ta 17^V[U؋B<^^[vb840 two-arg-- -DT! @ two-arg-* cos  .inv-aref_1d / lUV\E]IË׉]ԋ^W'׃WWhRP3ɱ;}t ]ԋS;tL]jEE[tW;uЋW+EȋE̋׃WWujuRP3ɱEEVEE^W'ЋF^"W'E3۳Ë Uu;E~,!ËW+>U]܋W';unEɋuËÉ]܋^&׉EЋV*E܋^.W'ЋF2^6W'U؋^6W'^:׋U uRP^>3Ҳ]upE܋]^]떋Ë_W' n02{Nn ar j nph p3+ @ improvedroots 1- aref_1d into  .inv-aref_1d improve_lg_root /UVTE~G]]jE̋]ẺU܋ u1;~BjEЋE܋^׉EċEЋ ]uP;E~dYjE븋Ӊ]؋W';t]jEԉ]؋E܋^׉EȋEԋ ]un;EsUċW';u EɋuËEUЋ^W'^׋U uRP^3Ҳ]EЋE=UȋW';u 3Ҳ]ugpaE̋]؋UEUԋ^W'ЋE܋^"W'U uRP^3Ҳ]u(p"EԋE Ë_W'rË_W'땋Ë_W'Ԑ an &optional iter0 2{N iter a n j ! k  p3+ @improve_lg_root  legendre_pd =  two-arg-/  two-arg-- newton iteration failed at ~s/p8UV4U^W'r3A׉E‰E3ҋ^W';tE3ҋ^W';t"EɋuËEU܋^W'ЋE^W'UF"W' nguess02{Nguess val deriv p3+ @improvedrootsbf 1- aref_1d  legendre_pdbf progn structure-allocate-instance quote aqd "M!.inv-structure-ref !!r ?"Ndetermine-foreign-address !<"! !< "/UVEjE‰E܋^׉E؋E ]u;E~EɋuËU؋W';tEU^W'ЋE܋^W'r3A׉EȉUċEU^W'@Eԋ]ȋ[UċRF@v"E]UEMI $MIL$MIL$śoL$ 7iy jA7]7YSu79y ta 17^[EU]^W'@v&EEMԍI $MIL$MIL$śoL$ 7iy jA7]7YSu79y ta 17EU^W'3Ҳ]up EE<Ë_W' an0(2{Na j n a1 val deriv a2 !!tt !p3+ @improve_lg_rootbf  legendre_pdbf progn structure-allocate-instance quote aqd "M!.inv-structure-ref !!r ?"Ndetermine-foreign-address !<"!!< "/XUVtU^W'r3A׋]EԋCE؋Eԋ@REЋF@vẺUȉEċMЍI $MȍIL$MčIL$śoL$ 7iy jA7]̋7YSu79y ta 17^[UBvE]EM؍I $MIL$MIL$śoL$ 7iy jA7]7YSu79y ta 17EɋuÐ nguess0( 2{N guess a1 val ! a2 tt ! !  p3+ @ legendre_pdwt into and typep k quote fixnum < !assert-1 !1- abs =  two-arg-*  two-arg-/ progn structure-allocate-instance !aqd "M!.inv-structure-ref !!r @"N!!!!"M!!!!r @"Ndetermine-foreign-address !< "!!!!"M!!!!r ?"N!!!!"M!!!!r @"N!!< "!!!!"M!!!!r ?"N!!< "!!!!"M!!!!r @"N!!< "!!!!"M!!!!r ?"N!!!!"M!!!!r @"N!!< "!!!!"M!!!!r @"N!!< "!!!!"M!!!!r @"N!!< "!!!!"M!!!!r ?"N!!<"!!!!"M!!!!r ?"N!!< "!!!!"M!!!!r @"N!!< "!!!!"M!!! !r @"N!!< "!!!!"M!!!!!r ?"N!!< "!!!!"M!!!"!r ?"N!!!!"M!!!#!r @"N!!< "!!!!"M!!!$!r @"N!!< "!!!!"M!!!%!r @"N!!< "!!!!"M!!!&!r @"N!!< "!!!!"M!!!'!r @"N!!< "!!!!"M!!!(!r ?"N!!!!"M!!!)!r @"N!!< "!!!!"M!!!*!r ?"N!!<"/42UVU܉EЋ^׉E]t׋FWWWRP^Ӌ^3W';tڋE܋^׉EuE3^׉Ej`3^u܏EjE̋EЋ^"׉EȋE̋ ] ;EE܋^&׋Ћ^*3W';: E^"׋ЋE^.W']Ѓ.3Ҳ^.W'^23ҲW'Eċ^6[UB^:[v>I $IL$IL$śoL$ 7iy jA77YSu79y ta 17^:[VBRFF@vJI $IL$IL$śoL$ 7iy jA77YSu79y ta 17^F[U܋RFN@vRI $IL$IL$śoL$ 7iy jA77YSu79y ta 17^N[U؋RFV@vZ|I $IL$|IL$śoL$ 7iy jA77YSu79y ta 17^V[]]SF^@xVbRvftplxI $pIL$lIL$śoL$ 7iy jA7t7YSu79y ta 17^b[UԋRFj@vnhd`\dI $`IL$\IL$śoL$ 7iy jA7h7YSu79y ta 17^j[VrRvvXTPMI $TIL$PIL$śoL$ 7iy jA7X7YSu79y ta 17^r[URFz@v~LHD@HI $DIL$@IL$śoL$ 7iy jA7L7YSu79y ta 17EċVz^.W'ЋE^.W'Ћ^2W'ɋuËUȋW';E̋^׋ЉU3^.W'^"׋@]܋S\@XTP\I $TIL$PIL$śoL$ 7iy jA7X7YSu79y ta 17[U؋BL[HD@LI $DIL$@IL$śoL$ 7iy jA7H7YSu79y ta 17[E]^"׋@]ԋ[R<8408I $4IL$0IL$śoL$ 7iy jA7<7YSu79y ta 17[R,($MI $(IL$$IL$śoL$ 7iy jA7,7YSu79y ta 17U^2W'u؏EԉE3Ҳ]5+E̋E]S@ R I $IL$IL$śoL$ 7iy jA77YSu79y ta 17[]]܋SE؉P@ I $IL$IL$śoL$ 7iy jA7 7YSu79y ta 17[UԋB[I $IL$IL$śoL$ 7iy jA77YSu79y ta 17[RMI $IL$IL$śoL$ 7iy jA77YSu79y ta 17[]]܋SCRI $IL$IL$śoL$ 7iy jA77YSu79y ta 17[R@I $IL$IL$śoL$ 7iy jA77YSu79y ta 17[RMI $IL$IL$śoL$ 7iy jA77YSu79y ta 17EUË_W' kx0(j2{Njkf x t1 t0 !ii Ӝi a1 !1!1!1!1a2 tt !1!2!3!1!2!3!2!3!1!2!3!1!2!3!1!2!3!2!3!1!2!3!1!2!3!2!3!1!2!3!1!2!3!1!2!3!1!2!3!1!2!3!1!2!3!2!3!1!2!3kp3+ @A  legendre_pdwt / UV$][[ЋC^W'ɋu q0 2{Nclosure np1 pp3+ @A /(UV4][[3Eu3Aɋu q0 2{Nclosure fun pp3+ @ int_gs_l1bf lexically-visible-function < !< setroots  improvedroots  two-arg-+ into improvedrootsbf array map 1- aref_1d progn structure-allocate-instance quote aqd "M! .inv-structure-ref !! r ?"Ndetermine-foreign-address ! < "!! < "/tUVUQtxYt^CEGE]^CEGE]QtxYE^׋U jRP^׉E؋E3Ҳ^"W'|3^&׉EjEċE؋U^*W'3Ҳ]p~EċE3Ҳu};~̍x]]F.Ӄ uRP^2׉E܍x]]F.Ӄ uRP^2׉EjEԋE^6׉E̋Eԋ ]u*;E~>3Ë_W'uW';v=ŰW';u EɋuË]SE܉UȋUԋ^:W'@EEЋUԋ^:W'@^>[vBEE]MI $MIL$MIL$śoL$ 7iy jA7]7YSu79y ta 17^>[UBvFE]EMȍI $MIL$MIL$śoL$ 7iy jA7]7YSu79y ta 173Ҳ]up EԋELË_W' funn0(2{Nsum n weights  abscissae i vals a1 !!a2 tt !!closure np1 |p!fun tpp3+ @internal  int_gs_l1bf LN!!LN!!< << <W@ c_qd_copylanguage c LN int_gs_l2bf  two-arg-+ into ?? two-arg-/ ? two-arg-- -DT! @ two-arg-* cos improve_lg_root improve_lg_rootbf determine-foreign-address quote <" legendre_pdwt progn structure-allocate-instance ! aqd "M!.inv-structure-ref !! r @"N! ! < "!!! !"M!!!! r ?"N! ! < "! ! <"/UVE܋‰E3Ҳ^W'E3^׉EVE^W'ЋF^"W'E3۳UÉ]̋ u ;~EɋuÉUԋW';tUԉUԋV&E̋^*W'ЋF.^2W'UЋ^2W'^6jEЋE^:W'E3Ҳ]EE3Ҳ;Eȋ^jEċЋEUȋ^>W'@]ȋSvBEEUMI $MIL$śoL$7iy jA7]7YSu79y ta 173Ҳ]u.p(.Ë_W'%W';.UË_W'EċE3Ҳ; ]SE؉UUȋ^FW'@EEȋu3A@^J[vNEE]MI $MIL$MIL$śoL$ 7iy jA7]7YSu79y ta 17^J[VRRvVE]UMI $MIL$MIL$śoL$ 7iy jA7]7YSu79y ta 17^R[UBvZE|x|I $xIL$śoL$7iy jA7]7YSu79y ta 173Ҳ]u-p'E̋]̋UlW';EË_W'ϐ funn0(2{Nsum n fun np1 Ӟnph j v k a1 !!!tt a2 !!!!!p3+ @ int_gs_l3bf into 1- aref_1d  two-arg-*  two-arg-+ /zUVLEU3^׉E܋EIuBE3 ]̋[];}u^3ۉE؋Ë Uu;EEɋuËʋU؉]ԋW';t؋EUԋ^W'u3AEЋEUԋ^W'ЋEЋ^W'ЋE܋^"W'E3Ҳ]upEԋ]ԋE`]뛋Ë_W' fun abscissaeweights02{Nfun  abscissae sum n i p3+ @A legendre_pdwt / UV$][[Ћ^W'ɋuÐ q0 2{Nclosure np1 pp3+ @ ab_and_wts lexically-visible-function <setroots  improvedroots  two-arg-+ improvedrootsbf array map /dUVT^CEGẺ]ЉE^׋U jRP^׉EE3Ҳ^W'EjE܋EU^"W'3Ҳ]uDp>E܋E3Ҳu:;~Ѝ]ĉ]ԍ]ڋF&Ӄ uRP^*׋ЋE3ɱɋuËË_W'븋W';t닐 n02{Nn  abscissae i closure np1 pp3+ @internal  ab_and_wts LN!!<<W@ setroots2 calc-array-size ? two-arg-+ ? two-arg-/ 1- > two-arg-- -DT! @ two-arg-* cos  .inv-aref_1d / |UVdEԺcW'E]IË׉]Ћ^W'׃WWhRP3ɱ;}t ]ЋS;tL]jEE[tW;uЋW+EċEȋ׃WWujuRP3ɱEEVEEԋ^W'ЋF^"W'E3۳Ë Uu;E~,!ËW+>U]܋W';unEɋuËÉ]܋^&׉E̋V*E܋^.W'ЋF2^6W'U؋^6W'^:׋U uRP^>3Ҳ]upE܋]^]떋Ë_W' n0 2{N halfn ar j nph n  p3+ @improvedrootsbf2 1- aref_1d  legendre_pdbf progn structure-allocate-instance quote aqd "M!.inv-structure-ref !!r ?"Ndetermine-foreign-address !<"! !< "/UVEjE‰EܺcW'^׉E؋E ]u;E~EɋuËU؋W';tEU^W'ЋE܋^W'r3A׉EȉUċEU^W'@Eԋ]ȋ[UċRF@v"E]UEMI $MIL$MIL$śoL$ 7iy jA7]7YSu79y ta 17^[EU]^W'@v&EEMԍI $MIL$MIL$śoL$ 7iy jA7]7YSu79y ta 17EU^W'3Ҳ]up EE<Ë_W' an0(2{Na j n a1 val deriv a2 !!tt !p3+ @improvedroots2 1- aref_1d into  .inv-aref_1d improve_lg_root /UVTE~U]]jE̋]ẺU܋ u?;~PjEЋEܺcW'^׉EċEЋ ]u^;E~rgjE몋Ӊ]؋W';t]jEԉ]؋EܺcW'^׉EȋEԋ ]un;EsUċW';u EɋuËEUЋ^W'^׋U uRP^3Ҳ]EЋE/UȋW';u 3Ҳ]ugpaE̋]؋UEUԋ^W'ЋE܋^"W'U uRP^3Ҳ]u(p"EԋE Ë_W'rË_W'땋Ë_W'Ԑ an &optional iter0 2{N iter a n -j ! k .2 p3+ @ ab_and_wts2  setroots2 improvedroots2  two-arg-+  two-arg-- calc-array-size improvedrootsbf2 1- into  .inv-aref_1d aref_1d  legendre_pdwt  two-arg-* / UVE^׋U jRP^׉E؋E3Ҳ^W'EE3Ҳ^W'E̋EcW'E܋]IË׉]^"W'׃WWhRP3ɱ;}t ]S;tL]jEE[tW;uЋW+EE׃WWujuRP3ɱEԋEԋ]EԋIuDË׉]^"W'׃WWhRP3ɱ;}t%]S;tgËW+.ËW+뱋]jEE[tW;uЋW+EE׃WWujuRP3ɱEЋEjEEЋE؋U^&W'3Ҳ]uBpW'U uRP^2uЏEE̋ ]uY+MȉMpOEЋUȋ^6W'U uRP^23Ҳ]u6p0EȋEdUȋW'EiUȋW'E렋Ë_W' n02{Nn np1 halfn a  abscissae weights nm1 i <7:C!DFGp3+ @ int_gs_l4bf  two-arg-+ into ?? two-arg-/  legendre_pdwt  two-arg-* ? two-arg-- -DT! @cos improve_lg_root improve_lg_rootbf determine-foreign-address quote <"progn structure-allocate-instance ! aqd "M! .inv-structure-ref !! r @"N! ! < "! ! ! !"M! !!! r ?"N! ! < "! ! < "/PUVE‰E3Ҳ^W'E3^׉EVE܋^W'ЋF^"W'EԋEܺcW'3۳EЋË UuU;E~fE3ҲW'3u3AEE3ҋ^&W'ЋE^*W'ЋE^W'ɋuËUЉ]̋W';t]̋V.É]̋^2W'ЋF6^*W'Uԋ^*W'^:jEЋE܋^>W'E3Ҳ]u'p!EE3Ҳu;#$EvË_W'ՋW';E^jEЋE܉Uȋ^BW'@]ȋSvFEEUMI $MIL$śoL$7iy jA7]7YSu79y ta 173Ҳ]UKEE3ҲF;O]SE؉UċUȋ^&W'@EEȋu3A@EUȸ^*W'u3A@^J[vNEE]MI $MIL$MIL$śoL$ 7iy jA7]7YSu79y ta 17^J[VRRvVE|xMI $|IL$xIL$śoL$ 7iy jA7]7YSu79y ta 17^R[UBvZtplMčI $pIL$lIL$śoL$ 7iy jA7t7YSu79y ta 173Ҳ]u6p06UË_W'W';EË_W'E̋]̋E funn0(2{Nsum fun n np1 nph halfn j fv a1 !!k v !!tt a2 !!!!! p3+ Ag40847@<top-level-form < macro-calls declaim  direct-calls store-doc-forms  load-scm-info update-xref-info load compiler-proclaim internal-in-package sets *source-pathname*  references !  legendre_pd ! logand_2op 1+ abs = +_2op 1-  two-arg-*  two-arg--  two-arg-/ >_2op values eq assert-1 <= fixnump not !logand oddp + 1+ loop-really-desetq * - twobytwo / setf > when loop eql case list typep unless assert  legendre_pdbf !!!!!!!!!qd_div qd_sub qd_mul qd_add  structure-ref aqd-q !!!into !!!!!!!!#!%!&! !!!" with-temps dsetv !(!)!*!+!,!'!-!.!/!0setroots !! .inv-aref_1d !!cos !! two-arg-+ !length cdr !null calc-array-size cons listp !! !!!" .inv-s-aref !$!#!'!(!)!*+ !%!&!- make-array  improvedroots !!8!!<aref_1d improve_lg_root !!!! !!!"!Earef !'!(!)!* !J!!&!%!$!!!error !! improvedrootsbf !!!3!2!6!7!I!1!! !! !!!"!:!K!(!)!* improve_lg_rootbf !!:!!3!2!6!7!1 legendre_pdwt !!evenp ! !!!"!:!9!(!)!*!#!%!&!'!-!.!/!0!!2!!!!!!3!4!6!7!!!!!< !!!8 int_gs_l1bf !!K!:!'!(!)! !!!"!*!%!F !!5!4!I!6!7!map !!!M!8!>!;!H<!funcall <!!O int_gs_l2bf !!5!4!T!O qd_copy_into !N!6!7!!J!!!=!!!8!>!!:! !!!"!$!#!'!(!)!*!&!%!F int_gs_l3bf !! !!!"!K!#!%!F!'!(!)!* !!!I!T!!>!!?!!8  ab_and_wts !!(!)! !!!"!*!%!F !!S!!!!M!>!;!H <!!O setroots2 !! !!!"!E!$!#!'!(!)!*!F!%!&!-!G!!!<!!!=!!!>!!?!@!!A!B!C!Dash improvedrootsbf2 !!!3!2!6!7!I!1!!Z! !! !!!"!:!K!(!)!* improvedroots2 !!8!!<!I!J!Z!! !! !!!"!E!K!'!(!)!*  ab_and_wts2 !!!- !K!Elet* !'!(!)! !!!"!*!-!G!$!%!F!!!8!!!-_2op !O!I!<!!!![!?!@!!A!B!C!D!Z!!>!Y!\ int_gs_l4bf !!!!4!5!T!O!V!N!6!7!!J!!!=!!Z!!8!>!!!!:! !!!"!$!#!'!(!)!*!&!%!F!/ UV$F^׋ɋuÐ null 0+ @<top-level-form <H*QAg40848@< int_gs_l4bf operator !! ab_and_wts2 !!!improvedroots2 !!!improvedrootsbf2 !!! setroots2 !!! ab_and_wts !!! int_gs_l3bf !!! int_gs_l2bf !!! int_gs_l1bf !! ! legendre_pdwt !! !improve_lg_rootbf !! !improvedrootsbf !! !improve_lg_root !! ! improvedroots !!!setroots !!! legendre_pdbf !!! legendre_pd !!! load-scm-info / UV$F^׋ɋuÐ null 0+ @<top-level-form <H*QAg40849@< legendre_pd  .de-args. M legendre_pdbf !setroots ! improvedroots !improve_lg_root !improvedrootsbf !improve_lg_rootbf ! legendre_pdwt ! int_gs_l1bf ! int_gs_l2bf ! int_gs_l3bf ! ab_and_wts ! setroots2 !improvedrootsbf2 !improvedroots2 ! ab_and_wts2 ! int_gs_l4bf !store-doc-forms / UV$F^׋ɋuÐ null 0+ @<top-level-form <H*Q