ofs | hex dump | ascii |
---|
0000 | 21 3c 61 72 63 68 3e 0a 23 31 2f 32 30 20 20 20 20 20 20 20 20 20 20 20 31 36 39 36 39 37 30 36 | !<arch>.#1/20...........16969706 |
0020 | 39 33 20 20 35 30 31 20 20 20 32 30 20 20 20 20 31 30 30 36 34 34 20 20 32 38 20 20 20 20 20 20 | 93..501...20....100644..28...... |
0040 | 20 20 60 0a 5f 5f 2e 53 59 4d 44 45 46 20 53 4f 52 54 45 44 00 00 00 00 00 00 00 00 00 00 00 00 | ..`.__.SYMDEF.SORTED............ |
0060 | 23 31 2f 32 38 20 20 20 20 20 20 20 20 20 20 20 31 36 39 36 39 30 38 39 39 31 20 20 35 30 31 20 | #1/28...........1696908991..501. |
0080 | 20 20 32 30 20 20 20 20 31 30 30 36 34 34 20 20 36 39 32 20 20 20 20 20 20 20 60 0a 61 70 70 72 | ..20....100644..692.......`.appr |
00a0 | 6f 78 2c 64 65 72 69 76 61 74 69 76 65 2e 6c 69 73 70 00 00 00 00 00 00 28 69 6e 2d 70 61 63 6b | ox,derivative.lisp......(in-pack |
00c0 | 61 67 65 20 3a 6c 69 7a 66 63 6d 2e 61 70 70 72 6f 78 29 0a 0a 28 64 65 66 75 6e 20 63 65 6e 74 | age.:lizfcm.approx)..(defun.cent |
00e0 | 72 61 6c 2d 64 65 72 69 76 61 74 69 76 65 2d 61 74 20 28 66 20 78 20 26 6f 70 74 69 6f 6e 61 6c | ral-derivative-at.(f.x.&optional |
0100 | 20 28 64 65 6c 74 61 20 30 2e 30 31 29 29 0a 20 20 28 6c 65 74 2a 20 28 28 78 32 20 28 2b 20 78 | .(delta.0.01))...(let*.((x2.(+.x |
0120 | 20 64 65 6c 74 61 29 29 0a 20 20 20 20 20 20 20 20 20 28 78 31 20 28 2d 20 78 20 64 65 6c 74 61 | .delta))..........(x1.(-.x.delta |
0140 | 29 29 0a 20 20 20 20 20 20 20 20 20 28 79 32 20 28 61 70 70 6c 79 20 66 20 28 6c 69 73 74 20 78 | ))..........(y2.(apply.f.(list.x |
0160 | 32 29 29 29 0a 20 20 20 20 20 20 20 20 20 28 79 31 20 28 61 70 70 6c 79 20 66 20 28 6c 69 73 74 | 2)))..........(y1.(apply.f.(list |
0180 | 20 78 31 29 29 29 29 0a 20 20 20 20 28 2f 20 28 2d 20 79 32 20 79 31 29 0a 20 20 20 20 20 20 20 | .x1)))).....(/.(-.y2.y1)........ |
01a0 | 28 2d 20 78 32 20 78 31 29 29 29 29 0a 0a 28 64 65 66 75 6e 20 66 6f 72 77 61 72 64 2d 64 65 72 | (-.x2.x1))))..(defun.forward-der |
01c0 | 69 76 61 74 69 76 65 2d 61 74 20 28 66 20 78 20 26 6f 70 74 69 6f 6e 61 6c 20 28 64 65 6c 74 61 | ivative-at.(f.x.&optional.(delta |
01e0 | 20 30 2e 30 31 29 29 0a 20 20 28 6c 65 74 2a 20 28 28 78 32 20 28 2b 20 78 20 64 65 6c 74 61 29 | .0.01))...(let*.((x2.(+.x.delta) |
0200 | 29 0a 20 20 20 20 20 20 20 20 20 28 78 31 20 78 29 0a 20 20 20 20 20 20 20 20 20 28 79 32 20 28 | )..........(x1.x)..........(y2.( |
0220 | 61 70 70 6c 79 20 66 20 28 6c 69 73 74 20 78 32 29 29 29 0a 20 20 20 20 20 20 20 20 20 28 79 31 | apply.f.(list.x2)))..........(y1 |
0240 | 20 28 61 70 70 6c 79 20 66 20 28 6c 69 73 74 20 78 31 29 29 29 29 0a 20 20 20 20 28 2f 20 28 2d | .(apply.f.(list.x1)))).....(/.(- |
0260 | 20 79 32 20 79 31 29 0a 20 20 20 20 20 20 20 28 2d 20 78 32 20 78 31 29 29 29 29 0a 0a 28 64 65 | .y2.y1)........(-.x2.x1))))..(de |
0280 | 66 75 6e 20 62 61 63 6b 77 61 72 64 2d 64 65 72 69 76 61 74 69 76 65 2d 61 74 20 28 66 20 78 20 | fun.backward-derivative-at.(f.x. |
02a0 | 26 6f 70 74 69 6f 6e 61 6c 20 28 64 65 6c 74 61 20 30 2e 30 31 29 29 0a 20 20 28 6c 65 74 2a 20 | &optional.(delta.0.01))...(let*. |
02c0 | 28 28 78 32 20 78 29 0a 20 20 20 20 20 20 20 20 20 28 78 31 20 28 2d 20 78 20 64 65 6c 74 61 29 | ((x2.x)..........(x1.(-.x.delta) |
02e0 | 29 0a 20 20 20 20 20 20 20 20 20 28 79 32 20 28 61 70 70 6c 79 20 66 20 28 6c 69 73 74 20 78 32 | )..........(y2.(apply.f.(list.x2 |
0300 | 29 29 29 0a 20 20 20 20 20 20 20 20 20 28 79 31 20 28 61 70 70 6c 79 20 66 20 28 6c 69 73 74 20 | )))..........(y1.(apply.f.(list. |
0320 | 78 31 29 29 29 29 0a 20 20 20 20 28 2f 20 28 2d 20 79 32 20 79 31 29 0a 20 20 20 20 20 20 20 28 | x1)))).....(/.(-.y2.y1)........( |
0340 | 2d 20 78 32 20 78 31 29 29 29 29 0a 0a 0a 0a 0a 23 31 2f 32 38 20 20 20 20 20 20 20 20 20 20 20 | -.x2.x1)))).....#1/28........... |
0360 | 31 36 39 36 39 30 38 31 32 33 20 20 35 30 31 20 20 20 32 30 20 20 20 20 31 30 30 36 34 34 20 20 | 1696908123..501...20....100644.. |
0380 | 33 32 34 20 20 20 20 20 20 20 60 0a 61 70 70 72 6f 78 2c 6d 61 63 65 70 73 2e 6c 69 73 70 00 00 | 324.......`.approx,maceps.lisp.. |
03a0 | 00 00 00 00 00 00 00 00 28 69 6e 2d 70 61 63 6b 61 67 65 20 3a 6c 69 7a 66 63 6d 2e 61 70 70 72 | ........(in-package.:lizfcm.appr |
03c0 | 6f 78 29 0a 0a 28 64 65 66 75 6e 20 63 6f 6d 70 75 74 65 2d 6d 61 63 65 70 73 20 28 66 20 61 20 | ox)..(defun.compute-maceps.(f.a. |
03e0 | 69 6e 69 74 29 0a 20 20 28 6c 65 74 20 28 28 68 20 69 6e 69 74 29 0a 20 20 20 20 20 20 20 20 28 | init)...(let.((h.init).........( |
0400 | 65 72 72 20 69 6e 69 74 29 29 0a 20 20 20 20 28 6c 6f 6f 70 20 63 6f 6c 6c 65 63 74 20 28 6c 69 | err.init)).....(loop.collect.(li |
0420 | 73 74 20 61 20 68 20 65 72 72 29 0a 20 20 20 20 20 20 20 20 20 20 64 6f 0a 20 20 20 20 20 20 20 | st.a.h.err)...........do........ |
0440 | 20 20 20 28 73 65 74 66 20 68 20 28 2f 20 68 20 32 29 0a 20 20 20 20 20 20 20 20 20 20 20 20 20 | ...(setf.h.(/.h.2).............. |
0460 | 20 20 20 65 72 72 20 28 61 62 73 20 28 2d 20 28 66 75 6e 63 61 6c 6c 20 66 20 28 2b 20 61 20 68 | ...err.(abs.(-.(funcall.f.(+.a.h |
0480 | 29 29 0a 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 28 | )).............................( |
04a0 | 66 75 6e 63 61 6c 6c 20 66 20 61 29 29 29 29 0a 20 20 20 20 20 20 20 20 20 20 77 68 69 6c 65 20 | funcall.f.a))))...........while. |
04c0 | 28 3e 20 65 72 72 20 30 29 29 29 29 0a 0a 0a 0a 23 31 2f 32 38 20 20 20 20 20 20 20 20 20 20 20 | (>.err.0))))....#1/28........... |
04e0 | 31 36 39 36 39 30 38 39 39 39 20 20 35 30 31 20 20 20 32 30 20 20 20 20 31 30 30 36 34 34 20 20 | 1696908999..501...20....100644.. |
0500 | 32 32 38 20 20 20 20 20 20 20 60 0a 61 70 70 72 6f 78 2c 70 61 63 6b 61 67 65 2e 6c 69 73 70 00 | 228.......`.approx,package.lisp. |
0520 | 00 00 00 00 00 00 00 00 28 69 6e 2d 70 61 63 6b 61 67 65 20 3a 63 6c 2d 75 73 65 72 29 0a 28 64 | ........(in-package.:cl-user).(d |
0540 | 65 66 70 61 63 6b 61 67 65 20 6c 69 7a 66 63 6d 2e 61 70 70 72 6f 78 0a 20 20 28 3a 75 73 65 20 | efpackage.lizfcm.approx...(:use. |
0560 | 3a 63 6c 29 0a 20 20 28 3a 65 78 70 6f 72 74 20 3a 63 65 6e 74 72 61 6c 2d 64 65 72 69 76 61 74 | :cl)...(:export.:central-derivat |
0580 | 69 76 65 2d 61 74 0a 20 20 20 20 20 20 20 20 20 20 20 3a 66 6f 72 77 61 72 64 2d 64 65 72 69 76 | ive-at............:forward-deriv |
05a0 | 61 74 69 76 65 2d 61 74 0a 20 20 20 20 20 20 20 20 20 20 20 3a 62 61 63 6b 77 61 72 64 2d 64 65 | ative-at............:backward-de |
05c0 | 72 69 76 61 74 69 76 65 2d 61 74 0a 20 20 20 20 20 20 20 20 20 20 20 3a 63 6f 6d 70 75 74 65 2d | rivative-at............:compute- |
05e0 | 6d 61 63 65 70 73 29 29 0a 0a 0a 0a 0a 0a 0a 0a 23 31 2f 32 30 20 20 20 20 20 20 20 20 20 20 20 | maceps))........#1/20........... |
0600 | 31 36 39 36 39 36 39 32 37 32 20 20 35 30 31 20 20 20 32 30 20 20 20 20 31 30 30 36 34 34 20 20 | 1696969272..501...20....100644.. |
0620 | 31 35 35 36 20 20 20 20 20 20 60 0a 6c 69 7a 66 63 6d 2e 61 73 64 00 00 00 00 00 00 00 00 00 00 | 1556......`.lizfcm.asd.......... |
0640 | 28 61 73 64 66 3a 64 65 66 73 79 73 74 65 6d 20 22 6c 69 7a 66 63 6d 22 0a 20 20 20 20 20 20 20 | (asdf:defsystem."lizfcm"........ |
0660 | 20 20 20 20 20 20 20 20 20 3a 76 65 72 73 69 6f 6e 20 22 30 2e 31 2e 30 22 0a 20 20 20 20 20 20 | .........:version."0.1.0"....... |
0680 | 20 20 20 20 20 20 20 20 20 20 3a 61 75 74 68 6f 72 20 22 45 6c 69 7a 61 62 65 74 68 20 48 75 6e | ..........:author."Elizabeth.Hun |
06a0 | 74 22 0a 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 3a 6c 69 63 65 6e 73 65 20 22 4d 49 54 | t".................:license."MIT |
06c0 | 22 0a 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 3a 63 6f 6d 70 6f 6e 65 6e 74 73 0a 20 20 | ".................:components... |
06e0 | 20 20 20 20 20 20 20 20 20 20 20 20 20 20 28 28 3a 66 69 6c 65 20 22 75 74 69 6c 73 2c 77 69 74 | ..............((:file."utils,wit |
0700 | 68 69 6e 2d 72 61 6e 67 65 22 20 3a 64 65 70 65 6e 64 73 2d 6f 6e 20 28 22 75 74 69 6c 73 2c 70 | hin-range".:depends-on.("utils,p |
0720 | 61 63 6b 61 67 65 22 29 29 0a 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 28 3a 66 69 6c | ackage"))..................(:fil |
0740 | 65 20 22 75 74 69 6c 73 2c 74 61 62 6c 65 22 20 3a 64 65 70 65 6e 64 73 2d 6f 6e 20 28 22 75 74 | e."utils,table".:depends-on.("ut |
0760 | 69 6c 73 2c 70 61 63 6b 61 67 65 22 29 29 0a 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 | ils,package")).................. |
0780 | 28 3a 66 69 6c 65 20 22 75 74 69 6c 73 2c 70 61 63 6b 61 67 65 22 29 0a 20 20 20 20 20 20 20 20 | (:file."utils,package")......... |
07a0 | 20 20 20 20 20 20 20 20 20 28 3a 66 69 6c 65 20 22 61 70 70 72 6f 78 2c 6d 61 63 65 70 73 22 20 | .........(:file."approx,maceps". |
07c0 | 3a 64 65 70 65 6e 64 73 2d 6f 6e 20 28 22 61 70 70 72 6f 78 2c 70 61 63 6b 61 67 65 22 29 29 0a | :depends-on.("approx,package")). |
07e0 | 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 28 3a 66 69 6c 65 20 22 61 70 70 72 6f 78 2c | .................(:file."approx, |
0800 | 64 65 72 69 76 61 74 69 76 65 22 20 3a 64 65 70 65 6e 64 73 2d 6f 6e 20 28 22 61 70 70 72 6f 78 | derivative".:depends-on.("approx |
0820 | 2c 70 61 63 6b 61 67 65 22 29 29 0a 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 28 3a 66 | ,package"))..................(:f |
0840 | 69 6c 65 20 22 61 70 70 72 6f 78 2c 70 61 63 6b 61 67 65 22 29 0a 20 20 20 20 20 20 20 20 20 20 | ile."approx,package")........... |
0860 | 20 20 20 20 20 20 20 28 3a 66 69 6c 65 20 22 76 65 63 74 6f 72 2c 6c 65 61 73 74 2d 73 71 75 61 | .......(:file."vector,least-squa |
0880 | 72 65 73 22 29 0a 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 28 3a 66 69 6c 65 20 22 76 | res")..................(:file."v |
08a0 | 65 63 74 6f 72 2c 64 69 73 74 61 6e 63 65 22 20 3a 64 65 70 65 6e 64 73 2d 6f 6e 20 28 22 76 65 | ector,distance".:depends-on.("ve |
08c0 | 63 74 6f 72 2c 6e 6f 72 6d 22 20 22 76 65 63 74 6f 72 2c 70 61 63 6b 61 67 65 22 29 29 0a 20 20 | ctor,norm"."vector,package"))... |
08e0 | 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 28 3a 66 69 6c 65 20 22 76 65 63 74 6f 72 2c 6e 6f | ...............(:file."vector,no |
0900 | 72 6d 22 20 3a 64 65 70 65 6e 64 73 2d 6f 6e 20 28 22 76 65 63 74 6f 72 2c 70 61 63 6b 61 67 65 | rm".:depends-on.("vector,package |
0920 | 22 29 29 0a 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 28 3a 66 69 6c 65 20 22 76 65 63 | "))..................(:file."vec |
0940 | 74 6f 72 2c 70 61 63 6b 61 67 65 22 29 29 29 0a 0a 0a 28 61 73 64 66 3a 64 65 66 73 79 73 74 65 | tor,package")))...(asdf:defsyste |
0960 | 6d 20 22 6c 69 7a 66 63 6d 2f 74 65 73 74 73 22 0a 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 | m."lizfcm/tests"................ |
0980 | 20 3a 61 75 74 68 6f 72 20 22 45 6c 69 7a 61 62 65 74 68 20 48 75 6e 74 22 0a 20 20 20 20 20 20 | .:author."Elizabeth.Hunt"....... |
09a0 | 20 20 20 20 20 20 20 20 20 20 3a 6c 69 63 65 6e 73 65 20 22 4d 49 54 22 0a 20 20 20 20 20 20 20 | ..........:license."MIT"........ |
09c0 | 20 20 20 20 20 20 20 20 20 3a 64 65 70 65 6e 64 73 2d 6f 6e 0a 20 20 20 20 20 20 20 20 20 20 20 | .........:depends-on............ |
09e0 | 20 20 20 20 20 28 3a 66 69 76 65 61 6d 0a 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 | .....(:fiveam................... |
0a00 | 3a 6c 69 7a 66 63 6d 29 0a 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 3a 63 6f 6d 70 6f 6e | :lizfcm).................:compon |
0a20 | 65 6e 74 73 0a 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 28 28 3a 66 69 6c 65 20 22 74 65 | ents.................((:file."te |
0a40 | 73 74 73 2c 74 61 62 6c 65 22 20 3a 64 65 70 65 6e 64 73 2d 6f 6e 20 28 22 74 65 73 74 73 2c 73 | sts,table".:depends-on.("tests,s |
0a60 | 75 69 74 65 22 29 29 0a 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 28 3a 66 69 6c 65 20 | uite"))..................(:file. |
0a80 | 22 74 65 73 74 73 2c 6d 61 63 65 70 73 22 20 3a 64 65 70 65 6e 64 73 2d 6f 6e 20 28 22 74 65 73 | "tests,maceps".:depends-on.("tes |
0aa0 | 74 73 2c 73 75 69 74 65 22 29 29 0a 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 28 3a 66 | ts,suite"))..................(:f |
0ac0 | 69 6c 65 20 22 74 65 73 74 73 2c 61 70 70 72 6f 78 22 20 3a 64 65 70 65 6e 64 73 2d 6f 6e 20 28 | ile."tests,approx".:depends-on.( |
0ae0 | 22 74 65 73 74 73 2c 73 75 69 74 65 22 29 29 0a 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 | "tests,suite"))................. |
0b00 | 20 28 3a 66 69 6c 65 20 22 74 65 73 74 73 2c 76 65 63 74 6f 72 22 20 3a 64 65 70 65 6e 64 73 2d | .(:file."tests,vector".:depends- |
0b20 | 6f 6e 20 28 22 74 65 73 74 73 2c 73 75 69 74 65 22 29 29 0a 20 20 20 20 20 20 20 20 20 20 20 20 | on.("tests,suite"))............. |
0b40 | 20 20 20 20 20 28 3a 66 69 6c 65 20 22 74 65 73 74 73 2c 73 75 69 74 65 22 29 29 0a 20 20 20 20 | .....(:file."tests,suite"))..... |
0b60 | 20 20 20 20 20 20 20 20 20 20 20 20 3a 70 65 72 66 6f 72 6d 0a 20 20 20 20 20 20 20 20 20 20 20 | ............:perform............ |
0b80 | 20 20 20 20 20 28 61 73 64 66 3a 74 65 73 74 2d 6f 70 20 28 6f 20 63 29 20 28 75 69 6f 70 3a 73 | .....(asdf:test-op.(o.c).(uiop:s |
0ba0 | 79 6d 62 6f 6c 2d 63 61 6c 6c 0a 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 | ymbol-call...................... |
0bc0 | 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 3a 66 69 76 65 61 6d 20 3a 72 75 6e 21 0a 20 | .................:fiveam.:run!.. |
0be0 | 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 | ................................ |
0c00 | 20 20 20 20 20 28 75 69 6f 70 3a 66 69 6e 64 2d 73 79 6d 62 6f 6c 2a 20 3a 6c 69 7a 66 63 6d 2d | .....(uiop:find-symbol*.:lizfcm- |
0c20 | 74 65 73 74 2d 73 75 69 74 65 20 3a 6c 69 7a 66 63 6d 2f 74 65 73 74 73 29 29 29 29 0a 0a 0a 0a | test-suite.:lizfcm/tests)))).... |
0c40 | 23 31 2f 32 38 20 20 20 20 20 20 20 20 20 20 20 31 36 39 36 39 30 38 39 38 36 20 20 35 30 31 20 | #1/28...........1696908986..501. |
0c60 | 20 20 32 30 20 20 20 20 31 30 30 36 34 34 20 20 31 33 39 36 20 20 20 20 20 20 60 0a 74 65 73 74 | ..20....100644..1396......`.test |
0c80 | 73 2c 61 70 70 72 6f 78 2e 6c 69 73 70 00 00 00 00 00 00 00 00 00 00 00 28 64 65 66 70 61 63 6b | s,approx.lisp...........(defpack |
0ca0 | 61 67 65 20 6c 69 7a 66 63 6d 2f 74 65 73 74 73 2e 61 70 70 72 6f 78 0a 20 20 28 3a 75 73 65 20 | age.lizfcm/tests.approx...(:use. |
0cc0 | 3a 63 6c 0a 20 20 20 20 20 20 20 20 3a 66 69 76 65 61 6d 0a 20 20 20 20 20 20 20 20 3a 6c 69 7a | :cl.........:fiveam.........:liz |
0ce0 | 66 63 6d 2e 61 70 70 72 6f 78 0a 20 20 20 20 20 20 20 20 3a 6c 69 7a 66 63 6d 2e 75 74 69 6c 73 | fcm.approx.........:lizfcm.utils |
0d00 | 0a 20 20 20 20 20 20 20 20 3a 6c 69 7a 66 63 6d 2f 74 65 73 74 73 29 0a 20 20 28 3a 65 78 70 6f | .........:lizfcm/tests)...(:expo |
0d20 | 72 74 20 3a 61 70 70 72 6f 78 2d 73 75 69 74 65 29 29 0a 28 69 6e 2d 70 61 63 6b 61 67 65 20 3a | rt.:approx-suite)).(in-package.: |
0d40 | 6c 69 7a 66 63 6d 2f 74 65 73 74 73 2e 61 70 70 72 6f 78 29 0a 0a 28 64 65 66 2d 73 75 69 74 65 | lizfcm/tests.approx)..(def-suite |
0d60 | 20 61 70 70 72 6f 78 2d 73 75 69 74 65 0a 20 20 20 20 20 20 20 20 20 20 20 3a 69 6e 20 6c 69 7a | .approx-suite............:in.liz |
0d80 | 66 63 6d 2d 74 65 73 74 2d 73 75 69 74 65 29 0a 28 69 6e 2d 73 75 69 74 65 20 61 70 70 72 6f 78 | fcm-test-suite).(in-suite.approx |
0da0 | 2d 73 75 69 74 65 29 0a 0a 28 74 65 73 74 20 63 65 6e 74 72 61 6c 2d 64 65 72 69 76 61 74 69 76 | -suite)..(test.central-derivativ |
0dc0 | 65 2d 61 74 0a 20 20 20 20 20 20 3a 64 65 73 63 72 69 70 74 69 6f 6e 20 22 64 65 72 69 76 61 74 | e-at.......:description."derivat |
0de0 | 69 76 65 20 61 74 20 69 73 20 77 69 74 68 69 6e 20 62 6f 75 6e 64 73 22 0a 20 20 20 20 20 20 28 | ive.at.is.within.bounds".......( |
0e00 | 6c 65 74 20 28 28 66 20 28 6c 61 6d 62 64 61 20 28 78 29 20 28 2a 20 78 20 78 29 29 29 0a 20 20 | let.((f.(lambda.(x).(*.x.x)))... |
0e20 | 20 20 20 20 20 20 20 20 20 20 28 78 20 32 29 0a 20 20 20 20 20 20 20 20 20 20 20 20 28 61 63 63 | ..........(x.2).............(acc |
0e40 | 65 70 74 65 64 2d 64 65 6c 74 61 20 30 2e 30 32 29 0a 20 20 20 20 20 20 20 20 20 20 20 20 28 66 | epted-delta.0.02).............(f |
0e60 | 2d 70 72 69 6d 65 2d 61 74 2d 78 20 34 29 0a 20 20 20 20 20 20 20 20 20 20 20 20 28 64 65 6c 74 | -prime-at-x.4).............(delt |
0e80 | 61 20 30 2e 30 31 29 29 0a 20 20 20 20 20 20 20 20 28 69 73 20 28 77 69 74 68 69 6e 2d 72 61 6e | a.0.01)).........(is.(within-ran |
0ea0 | 67 65 2d 70 0a 20 20 20 20 20 20 20 20 20 20 20 20 20 20 28 63 65 6e 74 72 61 6c 2d 64 65 72 69 | ge-p...............(central-deri |
0ec0 | 76 61 74 69 76 65 2d 61 74 20 66 20 78 20 64 65 6c 74 61 29 0a 20 20 20 20 20 20 20 20 20 20 20 | vative-at.f.x.delta)............ |
0ee0 | 20 20 20 66 2d 70 72 69 6d 65 2d 61 74 2d 78 0a 20 20 20 20 20 20 20 20 20 20 20 20 20 20 61 63 | ...f-prime-at-x...............ac |
0f00 | 63 65 70 74 65 64 2d 64 65 6c 74 61 29 29 29 29 0a 0a 28 74 65 73 74 20 66 77 64 2d 64 65 72 69 | cepted-delta))))..(test.fwd-deri |
0f20 | 76 61 74 69 76 65 2d 61 74 0a 20 20 20 20 20 20 3a 64 65 73 63 72 69 70 74 69 6f 6e 20 22 66 6f | vative-at.......:description."fo |
0f40 | 72 77 61 72 64 20 64 65 72 69 76 61 74 69 76 65 20 61 74 20 69 73 20 77 69 74 68 69 6e 20 62 6f | rward.derivative.at.is.within.bo |
0f60 | 75 6e 64 73 22 0a 20 20 20 20 20 20 28 6c 65 74 20 28 28 66 20 28 6c 61 6d 62 64 61 20 28 78 29 | unds".......(let.((f.(lambda.(x) |
0f80 | 20 28 2a 20 78 20 78 29 29 29 0a 20 20 20 20 20 20 20 20 20 20 20 20 28 78 20 32 29 0a 20 20 20 | .(*.x.x))).............(x.2).... |
0fa0 | 20 20 20 20 20 20 20 20 20 28 61 63 63 65 70 74 65 64 2d 64 65 6c 74 61 20 30 2e 30 32 29 0a 20 | .........(accepted-delta.0.02).. |
0fc0 | 20 20 20 20 20 20 20 20 20 20 20 28 66 2d 70 72 69 6d 65 2d 61 74 2d 78 20 34 29 0a 20 20 20 20 | ...........(f-prime-at-x.4)..... |
0fe0 | 20 20 20 20 20 20 20 20 28 64 65 6c 74 61 20 30 2e 30 31 29 29 0a 20 20 20 20 20 20 20 20 28 69 | ........(delta.0.01)).........(i |
1000 | 73 20 28 77 69 74 68 69 6e 2d 72 61 6e 67 65 2d 70 0a 20 20 20 20 20 20 20 20 20 20 20 20 20 20 | s.(within-range-p............... |
1020 | 28 66 6f 72 77 61 72 64 2d 64 65 72 69 76 61 74 69 76 65 2d 61 74 20 66 20 78 20 64 65 6c 74 61 | (forward-derivative-at.f.x.delta |
1040 | 29 0a 20 20 20 20 20 20 20 20 20 20 20 20 20 20 66 2d 70 72 69 6d 65 2d 61 74 2d 78 0a 20 20 20 | )...............f-prime-at-x.... |
1060 | 20 20 20 20 20 20 20 20 20 20 20 61 63 63 65 70 74 65 64 2d 64 65 6c 74 61 29 29 29 29 0a 0a 28 | ...........accepted-delta))))..( |
1080 | 74 65 73 74 20 62 77 64 2d 64 65 72 69 76 61 74 69 76 65 2d 61 74 0a 20 20 20 20 20 20 3a 64 65 | test.bwd-derivative-at.......:de |
10a0 | 73 63 72 69 70 74 69 6f 6e 20 22 62 61 63 6b 77 61 72 64 20 64 65 72 69 76 61 74 69 76 65 20 61 | scription."backward.derivative.a |
10c0 | 74 20 69 73 20 77 69 74 68 69 6e 20 62 6f 75 6e 64 73 22 0a 20 20 20 20 20 20 28 6c 65 74 20 28 | t.is.within.bounds".......(let.( |
10e0 | 28 66 20 28 6c 61 6d 62 64 61 20 28 78 29 20 28 2a 20 78 20 78 29 29 29 0a 20 20 20 20 20 20 20 | (f.(lambda.(x).(*.x.x)))........ |
1100 | 20 20 20 20 20 28 78 20 32 29 0a 20 20 20 20 20 20 20 20 20 20 20 20 28 61 63 63 65 70 74 65 64 | .....(x.2).............(accepted |
1120 | 2d 64 65 6c 74 61 20 30 2e 30 32 29 0a 20 20 20 20 20 20 20 20 20 20 20 20 28 66 2d 70 72 69 6d | -delta.0.02).............(f-prim |
1140 | 65 2d 61 74 2d 78 20 34 29 0a 20 20 20 20 20 20 20 20 20 20 20 20 28 64 65 6c 74 61 20 30 2e 30 | e-at-x.4).............(delta.0.0 |
1160 | 31 29 29 0a 20 20 20 20 20 20 20 20 28 69 73 20 28 77 69 74 68 69 6e 2d 72 61 6e 67 65 2d 70 0a | 1)).........(is.(within-range-p. |
1180 | 20 20 20 20 20 20 20 20 20 20 20 20 20 20 28 62 61 63 6b 77 61 72 64 2d 64 65 72 69 76 61 74 69 | ..............(backward-derivati |
11a0 | 76 65 2d 61 74 20 66 20 78 20 64 65 6c 74 61 29 0a 20 20 20 20 20 20 20 20 20 20 20 20 20 20 66 | ve-at.f.x.delta)...............f |
11c0 | 2d 70 72 69 6d 65 2d 61 74 2d 78 0a 20 20 20 20 20 20 20 20 20 20 20 20 20 20 61 63 63 65 70 74 | -prime-at-x...............accept |
11e0 | 65 64 2d 64 65 6c 74 61 29 29 29 29 0a 0a 0a 0a 23 31 2f 32 38 20 20 20 20 20 20 20 20 20 20 20 | ed-delta))))....#1/28........... |
1200 | 31 36 39 36 39 30 38 31 32 34 20 20 35 30 31 20 20 20 32 30 20 20 20 20 31 30 30 36 34 34 20 20 | 1696908124..501...20....100644.. |
1220 | 31 30 39 32 20 20 20 20 20 20 60 0a 74 65 73 74 73 2c 6d 61 63 65 70 73 2e 6c 69 73 70 00 00 00 | 1092......`.tests,maceps.lisp... |
1240 | 00 00 00 00 00 00 00 00 28 64 65 66 70 61 63 6b 61 67 65 20 6c 69 7a 66 63 6d 2f 74 65 73 74 73 | ........(defpackage.lizfcm/tests |
1260 | 2e 6d 61 63 65 70 73 0a 20 20 28 3a 75 73 65 20 3a 63 6c 0a 20 20 20 20 20 20 20 20 3a 66 69 76 | .maceps...(:use.:cl.........:fiv |
1280 | 65 61 6d 0a 20 20 20 20 20 20 20 20 3a 6c 69 7a 66 63 6d 2e 61 70 70 72 6f 78 0a 20 20 20 20 20 | eam.........:lizfcm.approx...... |
12a0 | 20 20 20 3a 6c 69 7a 66 63 6d 2e 75 74 69 6c 73 0a 20 20 20 20 20 20 20 20 3a 6c 69 7a 66 63 6d | ...:lizfcm.utils.........:lizfcm |
12c0 | 2f 74 65 73 74 73 29 0a 20 20 28 3a 65 78 70 6f 72 74 20 3a 61 70 70 72 6f 78 2d 73 75 69 74 65 | /tests)...(:export.:approx-suite |
12e0 | 29 29 0a 28 69 6e 2d 70 61 63 6b 61 67 65 20 3a 6c 69 7a 66 63 6d 2f 74 65 73 74 73 2e 6d 61 63 | )).(in-package.:lizfcm/tests.mac |
1300 | 65 70 73 29 0a 0a 28 64 65 66 2d 73 75 69 74 65 20 6d 61 63 65 70 73 2d 73 75 69 74 65 0a 20 20 | eps)..(def-suite.maceps-suite... |
1320 | 20 20 20 20 20 20 20 20 20 3a 69 6e 20 6c 69 7a 66 63 6d 2d 74 65 73 74 2d 73 75 69 74 65 29 0a | .........:in.lizfcm-test-suite). |
1340 | 28 69 6e 2d 73 75 69 74 65 20 6d 61 63 65 70 73 2d 73 75 69 74 65 29 0a 0a 28 74 65 73 74 20 6d | (in-suite.maceps-suite)..(test.m |
1360 | 61 63 65 70 73 0a 20 20 20 20 20 20 3a 64 65 73 63 72 69 70 74 69 6f 6e 20 22 64 6f 75 62 6c 65 | aceps.......:description."double |
1380 | 20 70 72 65 63 69 73 69 6f 6e 20 70 72 6f 76 69 64 65 73 20 70 72 65 63 69 73 69 6f 6e 20 61 62 | .precision.provides.precision.ab |
13a0 | 6f 75 74 20 28 6d 61 63 20 65 70 73 20 6f 66 20 73 69 6e 67 6c 65 20 70 72 65 63 69 73 69 6f 6e | out.(mac.eps.of.single.precision |
13c0 | 29 20 73 71 75 61 72 65 64 22 0a 20 20 20 20 20 20 28 6c 65 74 2a 20 28 28 6d 61 63 65 70 73 2d | ).squared".......(let*.((maceps- |
13e0 | 63 6f 6d 70 75 74 61 74 69 6f 6e 2d 64 6f 75 62 6c 65 20 28 63 6f 6d 70 75 74 65 2d 6d 61 63 65 | computation-double.(compute-mace |
1400 | 70 73 20 28 6c 61 6d 62 64 61 20 28 78 29 20 78 29 0a 20 20 20 20 20 20 20 20 20 20 20 20 20 20 | ps.(lambda.(x).x)............... |
1420 | 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 | ................................ |
1440 | 20 20 20 20 20 20 20 20 20 20 31 2e 30 64 30 0a 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 | ..........1.0d0................. |
1460 | 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 | ................................ |
1480 | 20 20 20 20 20 20 20 20 31 2e 30 64 30 29 29 0a 20 20 20 20 20 20 20 20 20 20 20 20 20 28 6d 61 | ........1.0d0))..............(ma |
14a0 | 63 65 70 73 2d 63 6f 6d 70 75 74 61 74 69 6f 6e 2d 73 69 6e 67 6c 65 20 28 63 6f 6d 70 75 74 65 | ceps-computation-single.(compute |
14c0 | 2d 6d 61 63 65 70 73 20 28 6c 61 6d 62 64 61 20 28 78 29 20 78 29 0a 20 20 20 20 20 20 20 20 20 | -maceps.(lambda.(x).x).......... |
14e0 | 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 | ................................ |
1500 | 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 31 2e 30 0a 20 20 20 20 20 20 20 20 20 20 20 20 20 | ...............1.0.............. |
1520 | 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 | ................................ |
1540 | 20 20 20 20 20 20 20 20 20 20 20 31 2e 30 29 29 0a 20 20 20 20 20 20 20 20 20 20 20 20 20 28 6c | ...........1.0))..............(l |
1560 | 61 73 74 2d 64 6f 75 62 6c 65 2d 68 20 28 63 61 64 61 72 20 28 6c 61 73 74 20 6d 61 63 65 70 73 | ast-double-h.(cadar.(last.maceps |
1580 | 2d 63 6f 6d 70 75 74 61 74 69 6f 6e 2d 64 6f 75 62 6c 65 29 29 29 0a 20 20 20 20 20 20 20 20 20 | -computation-double))).......... |
15a0 | 20 20 20 20 28 6c 61 73 74 2d 73 69 6e 67 6c 65 2d 68 20 28 63 61 64 61 72 20 28 6c 61 73 74 20 | ....(last-single-h.(cadar.(last. |
15c0 | 6d 61 63 65 70 73 2d 63 6f 6d 70 75 74 61 74 69 6f 6e 2d 73 69 6e 67 6c 65 29 29 29 29 0a 20 20 | maceps-computation-single))))... |
15e0 | 20 20 20 20 20 20 28 69 73 20 28 77 69 74 68 69 6e 2d 72 61 6e 67 65 2d 70 0a 20 20 20 20 20 20 | ......(is.(within-range-p....... |
1600 | 20 20 20 20 20 20 20 20 28 2d 20 6c 61 73 74 2d 64 6f 75 62 6c 65 2d 68 20 28 2a 20 6c 61 73 74 | ........(-.last-double-h.(*.last |
1620 | 2d 73 69 6e 67 6c 65 2d 68 20 6c 61 73 74 2d 73 69 6e 67 6c 65 2d 68 29 29 0a 20 20 20 20 20 20 | -single-h.last-single-h))....... |
1640 | 20 20 20 20 20 20 20 20 30 0a 20 20 20 20 20 20 20 20 20 20 20 20 20 20 6c 61 73 74 2d 73 69 6e | ........0...............last-sin |
1660 | 67 6c 65 2d 68 29 29 29 29 0a 0a 0a 0a 0a 0a 0a 23 31 2f 32 30 20 20 20 20 20 20 20 20 20 20 20 | gle-h)))).......#1/20........... |
1680 | 31 36 39 36 39 30 38 31 32 34 20 20 35 30 31 20 20 20 32 30 20 20 20 20 31 30 30 36 34 34 20 20 | 1696908124..501...20....100644.. |
16a0 | 32 36 30 20 20 20 20 20 20 20 60 0a 74 65 73 74 73 2c 73 75 69 74 65 2e 6c 69 73 70 00 00 00 00 | 260.......`.tests,suite.lisp.... |
16c0 | 28 69 6e 2d 70 61 63 6b 61 67 65 20 3a 63 6c 2d 75 73 65 72 29 0a 28 64 65 66 70 61 63 6b 61 67 | (in-package.:cl-user).(defpackag |
16e0 | 65 20 6c 69 7a 66 63 6d 2f 74 65 73 74 73 0a 20 20 28 3a 75 73 65 20 3a 63 6c 0a 20 20 20 20 20 | e.lizfcm/tests...(:use.:cl...... |
1700 | 20 20 20 3a 66 69 76 65 61 6d 29 0a 20 20 28 3a 65 78 70 6f 72 74 20 3a 72 75 6e 21 0a 20 20 20 | ...:fiveam)...(:export.:run!.... |
1720 | 20 20 20 20 20 20 20 20 3a 6c 69 7a 66 63 6d 2d 74 65 73 74 2d 73 75 69 74 65 29 29 0a 28 69 6e | ........:lizfcm-test-suite)).(in |
1740 | 2d 70 61 63 6b 61 67 65 20 3a 6c 69 7a 66 63 6d 2f 74 65 73 74 73 29 0a 0a 28 64 65 66 2d 73 75 | -package.:lizfcm/tests)..(def-su |
1760 | 69 74 65 20 6c 69 7a 66 63 6d 2d 74 65 73 74 2d 73 75 69 74 65 0a 20 20 20 20 3a 64 65 73 63 72 | ite.lizfcm-test-suite.....:descr |
1780 | 69 70 74 69 6f 6e 20 22 54 68 65 20 75 6c 74 69 6d 61 74 65 20 70 61 72 65 6e 74 20 74 65 73 74 | iption."The.ultimate.parent.test |
17a0 | 20 73 75 69 74 65 22 29 0a 0a 0a 0a 0a 0a 0a 0a 23 31 2f 32 30 20 20 20 20 20 20 20 20 20 20 20 | .suite")........#1/20........... |
17c0 | 31 36 39 36 39 30 38 31 32 34 20 20 35 30 31 20 20 20 32 30 20 20 20 20 31 30 30 36 34 34 20 20 | 1696908124..501...20....100644.. |
17e0 | 39 39 36 20 20 20 20 20 20 20 60 0a 74 65 73 74 73 2c 74 61 62 6c 65 2e 6c 69 73 70 00 00 00 00 | 996.......`.tests,table.lisp.... |
1800 | 28 64 65 66 70 61 63 6b 61 67 65 20 6c 69 7a 66 63 6d 2f 74 65 73 74 73 2e 74 61 62 6c 65 0a 20 | (defpackage.lizfcm/tests.table.. |
1820 | 20 28 3a 75 73 65 20 3a 63 6c 0a 20 20 20 20 20 20 20 20 3a 66 69 76 65 61 6d 0a 20 20 20 20 20 | .(:use.:cl.........:fiveam...... |
1840 | 20 20 20 3a 6c 69 7a 66 63 6d 2e 75 74 69 6c 73 0a 20 20 20 20 20 20 20 20 3a 6c 69 7a 66 63 6d | ...:lizfcm.utils.........:lizfcm |
1860 | 2f 74 65 73 74 73 29 0a 20 20 28 3a 65 78 70 6f 72 74 20 3a 61 70 70 72 6f 78 2d 73 75 69 74 65 | /tests)...(:export.:approx-suite |
1880 | 29 29 0a 28 69 6e 2d 70 61 63 6b 61 67 65 20 3a 6c 69 7a 66 63 6d 2f 74 65 73 74 73 2e 74 61 62 | )).(in-package.:lizfcm/tests.tab |
18a0 | 6c 65 29 0a 0a 28 64 65 66 2d 73 75 69 74 65 20 74 61 62 6c 65 2d 73 75 69 74 65 0a 20 20 20 20 | le)..(def-suite.table-suite..... |
18c0 | 20 20 20 20 20 20 20 3a 69 6e 20 6c 69 7a 66 63 6d 2d 74 65 73 74 2d 73 75 69 74 65 29 0a 28 69 | .......:in.lizfcm-test-suite).(i |
18e0 | 6e 2d 73 75 69 74 65 20 74 61 62 6c 65 2d 73 75 69 74 65 29 0a 0a 28 64 65 66 75 6e 20 66 69 62 | n-suite.table-suite)..(defun.fib |
1900 | 20 28 6e 29 0a 20 20 28 63 6f 6e 64 20 28 28 3c 20 6e 20 32 29 20 6e 29 0a 20 20 20 20 20 20 20 | .(n)...(cond.((<.n.2).n)........ |
1920 | 20 28 74 20 28 2b 20 28 66 69 62 20 28 2d 20 6e 20 31 29 29 20 28 66 69 62 20 28 2d 20 6e 20 32 | .(t.(+.(fib.(-.n.1)).(fib.(-.n.2 |
1940 | 29 29 29 29 29 29 0a 0a 28 74 65 73 74 20 74 61 62 6c 65 2d 6f 66 2d 66 69 62 2d 76 61 6c 73 0a | ))))))..(test.table-of-fib-vals. |
1960 | 20 20 20 20 20 20 3a 64 65 73 63 72 69 70 74 69 6f 6e 20 22 74 61 62 6c 65 20 67 65 6e 65 72 61 | ......:description."table.genera |
1980 | 74 65 73 20 63 6f 72 72 65 63 74 6c 79 22 0a 20 20 20 20 20 20 28 6c 65 74 2a 20 28 28 68 65 61 | tes.correctly".......(let*.((hea |
19a0 | 64 65 72 73 20 27 28 22 6e 22 20 22 66 69 62 28 6e 29 22 29 29 0a 20 20 20 20 20 20 20 20 20 20 | ders.'("n"."fib(n)"))........... |
19c0 | 20 20 20 28 6e 2d 76 61 6c 75 65 73 20 27 28 28 31 29 20 28 32 29 20 28 33 29 20 28 34 29 29 29 | ...(n-values.'((1).(2).(3).(4))) |
19e0 | 0a 20 20 20 20 20 20 20 20 20 20 20 20 20 28 65 78 70 65 63 74 65 64 20 60 28 28 22 6e 22 20 22 | ..............(expected.`(("n"." |
1a00 | 66 69 62 28 6e 29 22 29 0a 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 | fib(n)")........................ |
1a20 | 20 20 28 31 20 2c 28 66 69 62 20 31 29 29 0a 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 | ..(1.,(fib.1)).................. |
1a40 | 20 20 20 20 20 20 20 20 28 32 20 2c 28 66 69 62 20 32 29 29 0a 20 20 20 20 20 20 20 20 20 20 20 | ........(2.,(fib.2))............ |
1a60 | 20 20 20 20 20 20 20 20 20 20 20 20 20 20 28 33 20 2c 28 66 69 62 20 33 29 29 0a 20 20 20 20 20 | ..............(3.,(fib.3))...... |
1a80 | 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 28 34 20 2c 28 66 69 62 20 34 29 29 | ....................(4.,(fib.4)) |
1aa0 | 29 29 0a 20 20 20 20 20 20 20 20 20 20 20 20 20 28 74 61 62 6c 65 64 20 28 6c 69 7a 66 63 6d 2e | ))..............(tabled.(lizfcm. |
1ac0 | 75 74 69 6c 73 3a 74 61 62 6c 65 20 28 3a 68 65 61 64 65 72 73 20 68 65 61 64 65 72 73 0a 20 20 | utils:table.(:headers.headers... |
1ae0 | 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 | ................................ |
1b00 | 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 3a 64 6f 6d 61 69 6e 2d 6f 72 64 65 72 20 28 | .................:domain-order.( |
1b20 | 6e 29 0a 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 | n).............................. |
1b40 | 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 3a 64 6f 6d 61 69 6e 2d 76 61 | ......................:domain-va |
1b60 | 6c 75 65 73 20 6e 2d 76 61 6c 75 65 73 29 0a 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 | lues.n-values).................. |
1b80 | 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 28 66 69 62 20 6e 29 29 | ........................(fib.n)) |
1ba0 | 29 29 0a 20 20 20 20 20 20 20 20 28 69 73 20 28 65 71 75 61 6c 20 65 78 70 65 63 74 65 64 20 74 | )).........(is.(equal.expected.t |
1bc0 | 61 62 6c 65 64 29 29 29 29 0a 0a 0a 0a 0a 0a 0a 23 31 2f 32 38 20 20 20 20 20 20 20 20 20 20 20 | abled)))).......#1/28........... |
1be0 | 31 36 39 36 39 37 30 35 32 31 20 20 35 30 31 20 20 20 32 30 20 20 20 20 31 30 30 36 34 34 20 20 | 1696970521..501...20....100644.. |
1c00 | 31 33 35 36 20 20 20 20 20 20 60 0a 74 65 73 74 73 2c 76 65 63 74 6f 72 2e 6c 69 73 70 00 00 00 | 1356......`.tests,vector.lisp... |
1c20 | 00 00 00 00 00 00 00 00 28 64 65 66 70 61 63 6b 61 67 65 20 6c 69 7a 66 63 6d 2f 74 65 73 74 73 | ........(defpackage.lizfcm/tests |
1c40 | 2e 76 65 63 74 6f 72 0a 20 20 28 3a 75 73 65 20 3a 63 6c 0a 20 20 20 20 20 20 20 20 3a 66 69 76 | .vector...(:use.:cl.........:fiv |
1c60 | 65 61 6d 0a 20 20 20 20 20 20 20 20 3a 6c 69 7a 66 63 6d 2e 76 65 63 74 6f 72 0a 20 20 20 20 20 | eam.........:lizfcm.vector...... |
1c80 | 20 20 20 3a 6c 69 7a 66 63 6d 2e 75 74 69 6c 73 0a 20 20 20 20 20 20 20 20 3a 6c 69 7a 66 63 6d | ...:lizfcm.utils.........:lizfcm |
1ca0 | 2f 74 65 73 74 73 29 0a 20 20 28 3a 65 78 70 6f 72 74 20 3a 76 65 63 74 6f 72 2d 73 75 69 74 65 | /tests)...(:export.:vector-suite |
1cc0 | 29 29 0a 28 69 6e 2d 70 61 63 6b 61 67 65 20 3a 6c 69 7a 66 63 6d 2f 74 65 73 74 73 2e 76 65 63 | )).(in-package.:lizfcm/tests.vec |
1ce0 | 74 6f 72 29 0a 0a 28 64 65 66 2d 73 75 69 74 65 20 76 65 63 74 6f 72 2d 73 75 69 74 65 0a 20 20 | tor)..(def-suite.vector-suite... |
1d00 | 20 20 20 20 20 20 20 20 20 3a 69 6e 20 6c 69 7a 66 63 6d 2d 74 65 73 74 2d 73 75 69 74 65 29 0a | .........:in.lizfcm-test-suite). |
1d20 | 28 69 6e 2d 73 75 69 74 65 20 76 65 63 74 6f 72 2d 73 75 69 74 65 29 0a 0a 28 74 65 73 74 20 70 | (in-suite.vector-suite)..(test.p |
1d40 | 2d 6e 6f 72 6d 0a 20 20 20 20 20 20 3a 64 65 73 63 72 69 70 74 69 6f 6e 20 22 63 6f 6d 70 75 74 | -norm.......:description."comput |
1d60 | 65 73 20 70 2d 6e 6f 72 6d 22 0a 20 20 20 20 20 20 28 6c 65 74 20 28 28 76 20 27 28 31 20 31 29 | es.p-norm".......(let.((v.'(1.1) |
1d80 | 29 0a 20 20 20 20 20 20 20 20 20 20 20 20 28 6c 65 6e 67 74 68 20 28 73 71 72 74 20 32 29 29 0a | ).............(length.(sqrt.2)). |
1da0 | 20 20 20 20 20 20 20 20 20 20 20 20 28 32 2d 6e 6f 72 6d 20 28 70 2d 6e 6f 72 6d 20 32 29 29 29 | ............(2-norm.(p-norm.2))) |
1dc0 | 0a 20 20 20 20 20 20 20 20 28 69 73 20 28 77 69 74 68 69 6e 2d 72 61 6e 67 65 2d 70 20 28 66 75 | .........(is.(within-range-p.(fu |
1de0 | 6e 63 61 6c 6c 20 32 2d 6e 6f 72 6d 20 76 29 0a 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 | ncall.2-norm.v)................. |
1e00 | 20 20 20 20 20 20 20 20 20 20 20 20 6c 65 6e 67 74 68 0a 20 20 20 20 20 20 20 20 20 20 20 20 20 | ............length.............. |
1e20 | 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 30 2e 30 30 30 30 31 29 29 29 29 0a 0a 28 74 65 73 | ...............0.00001))))..(tes |
1e40 | 74 20 76 65 63 74 6f 72 2d 64 69 73 74 61 6e 63 65 0a 20 20 20 20 20 20 3a 64 65 73 63 72 69 70 | t.vector-distance.......:descrip |
1e60 | 74 69 6f 6e 20 22 63 6f 6d 70 75 74 65 73 20 64 69 73 74 61 6e 63 65 20 76 69 61 20 6e 6f 72 6d | tion."computes.distance.via.norm |
1e80 | 22 0a 20 20 20 20 20 20 28 6c 65 74 20 28 28 76 31 20 27 28 30 20 30 29 29 0a 20 20 20 20 20 20 | ".......(let.((v1.'(0.0))....... |
1ea0 | 20 20 20 20 20 20 28 76 32 20 27 28 31 20 31 29 29 0a 20 20 20 20 20 20 20 20 20 20 20 20 28 64 | ......(v2.'(1.1)).............(d |
1ec0 | 69 73 74 20 28 73 71 72 74 20 32 29 29 29 0a 20 20 20 20 20 20 20 20 28 69 73 20 28 77 69 74 68 | ist.(sqrt.2))).........(is.(with |
1ee0 | 69 6e 2d 72 61 6e 67 65 2d 70 20 28 64 69 73 74 61 6e 63 65 20 76 31 20 76 32 20 28 70 2d 6e 6f | in-range-p.(distance.v1.v2.(p-no |
1f00 | 72 6d 20 32 29 29 0a 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 | rm.2)).......................... |
1f20 | 20 20 20 64 69 73 74 0a 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 | ...dist......................... |
1f40 | 20 20 20 20 30 2e 30 30 30 30 31 29 29 29 29 0a 0a 28 74 65 73 74 20 6c 65 61 73 74 2d 73 71 75 | ....0.00001))))..(test.least-squ |
1f60 | 61 72 65 73 0a 20 20 20 20 20 20 3a 64 65 73 63 72 69 70 74 69 6f 6e 20 22 6c 65 61 73 74 20 73 | ares.......:description."least.s |
1f80 | 71 75 61 72 65 73 20 69 73 20 63 6f 72 72 65 63 74 20 65 6e 6f 75 67 68 22 0a 20 20 20 20 20 20 | quares.is.correct.enough"....... |
1fa0 | 28 6c 65 74 20 28 28 78 20 27 28 30 20 31 20 32 20 33 20 34 29 29 0a 20 20 20 20 20 20 20 20 20 | (let.((x.'(0.1.2.3.4)).......... |
1fc0 | 20 20 20 28 79 20 27 28 31 20 32 20 33 20 34 20 35 29 29 29 0a 20 20 20 20 20 20 20 20 28 64 65 | ...(y.'(1.2.3.4.5))).........(de |
1fe0 | 73 74 72 75 63 74 75 72 69 6e 67 2d 62 69 6e 64 20 28 6d 20 62 29 20 28 6c 69 7a 66 63 6d 2e 76 | structuring-bind.(m.b).(lizfcm.v |
2000 | 65 63 74 6f 72 3a 6c 65 61 73 74 2d 73 71 75 61 72 65 73 2d 72 65 67 20 78 20 79 29 0a 20 20 20 | ector:least-squares-reg.x.y).... |
2020 | 20 20 20 20 20 20 20 28 69 73 20 28 77 69 74 68 69 6e 2d 72 61 6e 67 65 2d 70 20 6d 20 31 20 30 | .......(is.(within-range-p.m.1.0 |
2040 | 2e 30 30 30 30 31 29 29 0a 20 20 20 20 20 20 20 20 20 20 28 69 73 20 28 77 69 74 68 69 6e 2d 72 | .00001))...........(is.(within-r |
2060 | 61 6e 67 65 2d 70 20 62 20 31 20 30 2e 30 30 30 30 31 29 29 29 29 0a 20 20 20 20 20 20 28 6c 65 | ange-p.b.1.0.00001)))).......(le |
2080 | 74 20 28 28 78 20 27 28 31 20 32 20 33 20 34 20 35 20 36 20 37 29 29 0a 20 20 20 20 20 20 20 20 | t.((x.'(1.2.3.4.5.6.7))......... |
20a0 | 20 20 20 20 28 79 20 27 28 30 2e 35 20 33 20 32 20 33 2e 35 20 35 20 36 20 37 2e 35 29 29 29 0a | ....(y.'(0.5.3.2.3.5.5.6.7.5))). |
20c0 | 20 20 20 20 20 20 20 20 28 64 65 73 74 72 75 63 74 75 72 69 6e 67 2d 62 69 6e 64 20 28 6d 20 62 | ........(destructuring-bind.(m.b |
20e0 | 29 20 28 6c 69 7a 66 63 6d 2e 76 65 63 74 6f 72 3a 6c 65 61 73 74 2d 73 71 75 61 72 65 73 2d 72 | ).(lizfcm.vector:least-squares-r |
2100 | 65 67 20 78 20 79 29 0a 20 20 20 20 20 20 20 20 20 20 28 69 73 20 28 77 69 74 68 69 6e 2d 72 61 | eg.x.y)...........(is.(within-ra |
2120 | 6e 67 65 2d 70 20 6d 20 31 20 30 2e 33 29 29 29 29 29 20 3b 3b 20 6a 75 73 74 20 61 20 67 75 65 | nge-p.m.1.0.3))))).;;.just.a.gue |
2140 | 73 74 69 6d 61 74 65 20 66 6f 72 20 62 65 73 74 20 66 69 74 0a 0a 0a 0a 23 31 2f 32 38 20 20 20 | stimate.for.best.fit....#1/28... |
2160 | 20 20 20 20 20 20 20 20 31 36 39 36 39 30 38 31 32 33 20 20 35 30 31 20 20 20 32 30 20 20 20 20 | ........1696908123..501...20.... |
2180 | 31 30 30 36 34 34 20 20 31 34 30 20 20 20 20 20 20 20 60 0a 75 74 69 6c 73 2c 70 61 63 6b 61 67 | 100644..140.......`.utils,packag |
21a0 | 65 2e 6c 69 73 70 00 00 00 00 00 00 00 00 00 00 28 69 6e 2d 70 61 63 6b 61 67 65 20 3a 63 6c 2d | e.lisp..........(in-package.:cl- |
21c0 | 75 73 65 72 29 0a 28 64 65 66 70 61 63 6b 61 67 65 20 6c 69 7a 66 63 6d 2e 75 74 69 6c 73 0a 20 | user).(defpackage.lizfcm.utils.. |
21e0 | 20 28 3a 75 73 65 20 3a 63 6c 29 0a 20 20 28 3a 65 78 70 6f 72 74 20 3a 77 69 74 68 69 6e 2d 72 | .(:use.:cl)...(:export.:within-r |
2200 | 61 6e 67 65 2d 70 0a 20 20 20 20 20 20 20 20 20 20 20 3a 74 61 62 6c 65 29 29 0a 0a 0a 0a 0a 0a | ange-p............:table))...... |
2220 | 23 31 2f 32 30 20 20 20 20 20 20 20 20 20 20 20 31 36 39 36 39 30 38 31 32 33 20 20 35 30 31 20 | #1/20...........1696908123..501. |
2240 | 20 20 32 30 20 20 20 20 31 30 30 36 34 34 20 20 33 36 34 20 20 20 20 20 20 20 60 0a 75 74 69 6c | ..20....100644..364.......`.util |
2260 | 73 2c 74 61 62 6c 65 2e 6c 69 73 70 00 00 00 00 28 69 6e 2d 70 61 63 6b 61 67 65 20 3a 6c 69 7a | s,table.lisp....(in-package.:liz |
2280 | 66 63 6d 2e 75 74 69 6c 73 29 0a 0a 28 64 65 66 6d 61 63 72 6f 20 74 61 62 6c 65 20 28 28 26 6b | fcm.utils)..(defmacro.table.((&k |
22a0 | 65 79 20 68 65 61 64 65 72 73 20 64 6f 6d 61 69 6e 2d 6f 72 64 65 72 20 64 6f 6d 61 69 6e 2d 76 | ey.headers.domain-order.domain-v |
22c0 | 61 6c 75 65 73 29 20 26 62 6f 64 79 20 62 6f 64 79 29 0a 20 20 60 28 63 6f 6e 73 0a 20 20 20 20 | alues).&body.body)...`(cons..... |
22e0 | 20 2c 68 65 61 64 65 72 73 0a 20 20 20 20 20 28 6d 61 70 63 61 72 20 28 6c 61 6d 62 64 61 20 28 | .,headers......(mapcar.(lambda.( |
2300 | 74 75 70 6c 65 29 0a 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 28 64 65 73 74 72 75 63 74 75 | tuple)................(destructu |
2320 | 72 69 6e 67 2d 62 69 6e 64 20 2c 64 6f 6d 61 69 6e 2d 6f 72 64 65 72 20 74 75 70 6c 65 0a 20 20 | ring-bind.,domain-order.tuple... |
2340 | 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 28 61 70 70 65 6e 64 20 74 75 70 6c 65 0a 20 20 20 | ...............(append.tuple.... |
2360 | 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 28 6c 69 73 74 0a 20 20 20 20 | ......................(list..... |
2380 | 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 2c 40 62 6f 64 79 29 29 29 | .......................,@body))) |
23a0 | 29 0a 20 20 20 20 20 20 20 20 20 20 20 20 20 2c 64 6f 6d 61 69 6e 2d 76 61 6c 75 65 73 29 29 29 | )..............,domain-values))) |
23c0 | 0a 0a 0a 0a 0a 0a 0a 0a 23 31 2f 32 38 20 20 20 20 20 20 20 20 20 20 20 31 36 39 36 39 30 38 31 | ........#1/28...........16969081 |
23e0 | 32 34 20 20 35 30 31 20 20 20 32 30 20 20 20 20 31 30 30 36 34 34 20 20 31 37 32 20 20 20 20 20 | 24..501...20....100644..172..... |
2400 | 20 20 60 0a 75 74 69 6c 73 2c 77 69 74 68 69 6e 2d 72 61 6e 67 65 2e 6c 69 73 70 00 00 00 00 00 | ..`.utils,within-range.lisp..... |
2420 | 28 69 6e 2d 70 61 63 6b 61 67 65 20 3a 6c 69 7a 66 63 6d 2e 75 74 69 6c 73 29 0a 0a 28 64 65 66 | (in-package.:lizfcm.utils)..(def |
2440 | 75 6e 20 77 69 74 68 69 6e 2d 72 61 6e 67 65 2d 70 20 28 78 20 74 72 75 65 2d 76 61 6c 75 65 20 | un.within-range-p.(x.true-value. |
2460 | 64 65 6c 74 61 29 0a 20 20 28 61 6e 64 20 28 3c 20 78 20 28 2b 20 74 72 75 65 2d 76 61 6c 75 65 | delta)...(and.(<.x.(+.true-value |
2480 | 20 64 65 6c 74 61 29 29 0a 20 20 20 20 20 20 20 28 3e 20 78 20 28 2d 20 74 72 75 65 2d 76 61 6c | .delta))........(>.x.(-.true-val |
24a0 | 75 65 20 64 65 6c 74 61 29 29 29 29 0a 0a 0a 0a 23 31 2f 32 38 20 20 20 20 20 20 20 20 20 20 20 | ue.delta))))....#1/28........... |
24c0 | 31 36 39 36 39 30 38 31 32 34 20 20 35 30 31 20 20 20 32 30 20 20 20 20 31 30 30 36 34 34 20 20 | 1696908124..501...20....100644.. |
24e0 | 31 37 32 20 20 20 20 20 20 20 60 0a 76 65 63 74 6f 72 2c 64 69 73 74 61 6e 63 65 2e 6c 69 73 70 | 172.......`.vector,distance.lisp |
2500 | 00 00 00 00 00 00 00 00 28 69 6e 2d 70 61 63 6b 61 67 65 20 3a 6c 69 7a 66 63 6d 2e 76 65 63 74 | ........(in-package.:lizfcm.vect |
2520 | 6f 72 29 0a 0a 28 64 65 66 75 6e 20 64 69 73 74 61 6e 63 65 20 28 76 31 20 76 32 20 6e 6f 72 6d | or)..(defun.distance.(v1.v2.norm |
2540 | 29 0a 20 20 28 6c 65 74 2a 20 28 28 64 20 28 6d 61 70 63 61 72 20 23 27 2d 20 76 31 20 76 32 29 | )...(let*.((d.(mapcar.#'-.v1.v2) |
2560 | 29 0a 20 20 20 20 20 20 20 20 20 28 6c 65 6e 67 74 68 20 28 66 75 6e 63 61 6c 6c 20 6e 6f 72 6d | )..........(length.(funcall.norm |
2580 | 20 64 29 29 29 0a 20 20 20 20 6c 65 6e 67 74 68 29 29 0a 0a 0a 0a 0a 0a 23 31 2f 32 30 20 20 20 | .d))).....length))......#1/20... |
25a0 | 20 20 20 20 20 20 20 20 31 36 39 36 39 30 38 31 32 34 20 20 35 30 31 20 20 20 32 30 20 20 20 20 | ........1696908124..501...20.... |
25c0 | 31 30 30 36 34 34 20 20 33 30 30 20 20 20 20 20 20 20 60 0a 76 65 63 74 6f 72 2c 6e 6f 72 6d 2e | 100644..300.......`.vector,norm. |
25e0 | 6c 69 73 70 00 00 00 00 28 69 6e 2d 70 61 63 6b 61 67 65 20 3a 6c 69 7a 66 63 6d 2e 76 65 63 74 | lisp....(in-package.:lizfcm.vect |
2600 | 6f 72 29 0a 0a 28 64 65 66 75 6e 20 70 2d 6e 6f 72 6d 20 28 70 29 0a 20 20 28 6c 61 6d 62 64 61 | or)..(defun.p-norm.(p)...(lambda |
2620 | 20 28 76 29 0a 20 20 20 20 28 65 78 70 74 0a 20 20 20 20 20 20 28 72 65 64 75 63 65 20 23 27 2b | .(v).....(expt.......(reduce.#'+ |
2640 | 0a 20 20 20 20 20 20 20 20 20 20 20 20 20 20 28 6d 61 70 63 61 72 20 28 6c 61 6d 62 64 61 20 28 | ...............(mapcar.(lambda.( |
2660 | 78 29 0a 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 28 61 62 73 0a | x).........................(abs. |
2680 | 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 28 65 78 70 74 20 | ..........................(expt. |
26a0 | 78 20 70 29 29 29 0a 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 76 29 29 | x.p))).......................v)) |
26c0 | 0a 20 20 20 20 20 20 28 2f 20 31 20 70 29 29 29 29 0a 0a 28 64 65 66 75 6e 20 6d 61 78 2d 6e 6f | .......(/.1.p))))..(defun.max-no |
26e0 | 72 6d 20 28 76 29 0a 20 20 28 72 65 64 75 63 65 20 23 27 6d 61 78 20 76 29 29 0a 0a 0a 0a 0a 0a | rm.(v)...(reduce.#'max.v))...... |
2700 | 23 31 2f 32 38 20 20 20 20 20 20 20 20 20 20 20 31 36 39 36 39 36 32 32 36 38 20 20 35 30 31 20 | #1/28...........1696962268..501. |
2720 | 20 20 32 30 20 20 20 20 31 30 30 36 34 34 20 20 31 37 32 20 20 20 20 20 20 20 60 0a 76 65 63 74 | ..20....100644..172.......`.vect |
2740 | 6f 72 2c 70 61 63 6b 61 67 65 2e 6c 69 73 70 00 00 00 00 00 00 00 00 00 28 69 6e 2d 70 61 63 6b | or,package.lisp.........(in-pack |
2760 | 61 67 65 20 3a 63 6c 2d 75 73 65 72 29 0a 28 64 65 66 70 61 63 6b 61 67 65 20 6c 69 7a 66 63 6d | age.:cl-user).(defpackage.lizfcm |
2780 | 2e 76 65 63 74 6f 72 0a 20 20 28 3a 75 73 65 20 3a 63 6c 29 0a 20 20 28 3a 65 78 70 6f 72 74 0a | .vector...(:use.:cl)...(:export. |
27a0 | 20 20 20 20 3a 70 2d 6e 6f 72 6d 0a 20 20 20 20 3a 6d 61 78 2d 6e 6f 72 6d 0a 20 20 20 20 3a 64 | ....:p-norm.....:max-norm.....:d |
27c0 | 69 73 74 61 6e 63 65 0a 20 20 20 20 3a 6c 65 61 73 74 2d 73 71 75 61 72 65 73 2d 72 65 67 29 29 | istance.....:least-squares-reg)) |
27e0 | 0a 0a 0a 0a 0a 0a 0a 0a 23 31 2f 32 30 20 20 20 20 20 20 20 20 20 20 20 31 36 39 36 39 37 30 36 | ........#1/20...........16969706 |
2800 | 33 39 20 20 35 30 31 20 20 20 32 30 20 20 20 20 31 30 30 36 34 34 20 20 31 37 33 32 20 20 20 20 | 39..501...20....100644..1732.... |
2820 | 20 20 60 0a 6d 61 69 6e 2e 6c 69 73 70 00 00 00 00 00 00 00 00 00 00 00 28 6c 6f 61 64 20 22 6c | ..`.main.lisp...........(load."l |
2840 | 69 7a 66 63 6d 2e 61 73 64 22 29 0a 28 71 6c 3a 71 75 69 63 6b 6c 6f 61 64 20 3a 6c 69 7a 66 63 | izfcm.asd").(ql:quickload.:lizfc |
2860 | 6d 29 0a 0a 3b 3b 20 74 68 69 73 20 69 73 20 61 20 63 6f 6c 6c 65 63 74 69 6f 6e 20 73 68 6f 77 | m)..;;.this.is.a.collection.show |
2880 | 63 61 73 69 6e 67 20 74 68 65 20 6c 69 62 72 61 72 79 20 64 65 76 65 6c 6f 70 65 64 20 66 6f 72 | casing.the.library.developed.for |
28a0 | 20 6d 61 74 68 34 36 31 30 2c 20 72 65 71 75 69 72 65 64 0a 3b 3b 20 66 72 6f 6d 20 74 68 65 20 | .math4610,.required.;;.from.the. |
28c0 | 53 68 61 72 65 64 20 4c 69 62 72 61 72 79 20 64 65 66 69 6e 69 74 69 6f 6e 0a 0a 28 64 65 66 75 | Shared.Library.definition..(defu |
28e0 | 6e 20 73 6d 61 63 65 70 73 20 28 29 0a 20 20 28 63 61 64 61 72 20 28 6c 61 73 74 20 28 6c 69 7a | n.smaceps.()...(cadar.(last.(liz |
2900 | 66 63 6d 2e 61 70 70 72 6f 78 3a 63 6f 6d 70 75 74 65 2d 6d 61 63 65 70 73 0a 20 20 20 20 20 20 | fcm.approx:compute-maceps....... |
2920 | 20 20 20 20 20 20 20 20 20 20 20 28 6c 61 6d 62 64 61 20 28 78 29 20 78 29 20 31 2e 30 20 31 2e | ...........(lambda.(x).x).1.0.1. |
2940 | 30 29 29 29 29 0a 0a 28 64 65 66 75 6e 20 64 6d 61 63 65 70 73 20 28 29 0a 20 20 28 63 61 64 61 | 0))))..(defun.dmaceps.()...(cada |
2960 | 72 20 28 6c 61 73 74 20 28 6c 69 7a 66 63 6d 2e 61 70 70 72 6f 78 3a 63 6f 6d 70 75 74 65 2d 6d | r.(last.(lizfcm.approx:compute-m |
2980 | 61 63 65 70 73 0a 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 28 6c 61 6d 62 64 61 20 28 | aceps..................(lambda.( |
29a0 | 78 29 20 78 29 20 31 2e 30 64 30 20 31 2e 30 64 30 29 29 29 29 0a 0a 28 64 65 66 75 6e 20 6c 32 | x).x).1.0d0.1.0d0))))..(defun.l2 |
29c0 | 2d 6e 6f 72 6d 20 28 76 29 0a 20 20 28 6c 65 74 20 28 28 32 2d 6e 6f 72 6d 20 28 6c 69 7a 66 63 | -norm.(v)...(let.((2-norm.(lizfc |
29e0 | 6d 2e 76 65 63 74 6f 72 3a 70 2d 6e 6f 72 6d 20 32 29 29 29 0a 20 20 20 20 28 66 75 6e 63 61 6c | m.vector:p-norm.2))).....(funcal |
2a00 | 6c 20 32 2d 6e 6f 72 6d 20 76 29 29 29 0a 0a 28 64 65 66 75 6e 20 6c 31 2d 6e 6f 72 6d 20 28 76 | l.2-norm.v)))..(defun.l1-norm.(v |
2a20 | 29 0a 20 20 28 6c 65 74 20 28 28 31 2d 6e 6f 72 6d 20 28 6c 69 7a 66 63 6d 2e 76 65 63 74 6f 72 | )...(let.((1-norm.(lizfcm.vector |
2a40 | 3a 70 2d 6e 6f 72 6d 20 31 29 29 29 0a 20 20 20 20 28 66 75 6e 63 61 6c 6c 20 31 2d 6e 6f 72 6d | :p-norm.1))).....(funcall.1-norm |
2a60 | 20 76 29 29 29 0a 0a 28 64 65 66 75 6e 20 6c 69 6e 66 2d 6e 6f 72 6d 20 28 76 29 0a 20 20 28 6c | .v)))..(defun.linf-norm.(v)...(l |
2a80 | 69 7a 66 63 6d 2e 76 65 63 74 6f 72 3a 6d 61 78 2d 6e 6f 72 6d 20 76 29 29 0a 0a 28 64 65 66 75 | izfcm.vector:max-norm.v))..(defu |
2aa0 | 6e 20 6c 32 2d 64 69 73 74 61 6e 63 65 20 28 76 31 20 76 32 29 0a 20 20 28 6c 65 74 20 28 28 32 | n.l2-distance.(v1.v2)...(let.((2 |
2ac0 | 2d 6e 6f 72 6d 20 28 6c 69 7a 66 63 6d 2e 76 65 63 74 6f 72 3a 70 2d 6e 6f 72 6d 20 32 29 29 29 | -norm.(lizfcm.vector:p-norm.2))) |
2ae0 | 0a 20 20 20 20 28 6c 69 7a 66 63 6d 2e 76 65 63 74 6f 72 3a 64 69 73 74 61 6e 63 65 20 76 31 20 | .....(lizfcm.vector:distance.v1. |
2b00 | 76 32 20 32 2d 6e 6f 72 6d 29 29 29 0a 0a 28 64 65 66 75 6e 20 6c 31 2d 64 69 73 74 61 6e 63 65 | v2.2-norm)))..(defun.l1-distance |
2b20 | 20 28 76 31 20 76 32 29 0a 20 20 28 6c 65 74 20 28 28 31 2d 6e 6f 72 6d 20 28 6c 69 7a 66 63 6d | .(v1.v2)...(let.((1-norm.(lizfcm |
2b40 | 2e 76 65 63 74 6f 72 3a 70 2d 6e 6f 72 6d 20 31 29 29 29 0a 20 20 20 20 28 6c 69 7a 66 63 6d 2e | .vector:p-norm.1))).....(lizfcm. |
2b60 | 76 65 63 74 6f 72 3a 64 69 73 74 61 6e 63 65 20 76 31 20 76 32 20 31 2d 6e 6f 72 6d 29 29 29 0a | vector:distance.v1.v2.1-norm))). |
2b80 | 0a 28 64 65 66 75 6e 20 6c 69 6e 66 2d 64 69 73 74 61 6e 63 65 20 28 76 31 20 76 32 29 0a 20 20 | .(defun.linf-distance.(v1.v2)... |
2ba0 | 28 6c 69 7a 66 63 6d 2e 76 65 63 74 6f 72 3a 64 69 73 74 61 6e 63 65 20 76 31 20 76 32 20 27 6c | (lizfcm.vector:distance.v1.v2.'l |
2bc0 | 69 7a 66 63 6d 2e 76 65 63 74 6f 72 3a 6d 61 78 2d 6e 6f 72 6d 29 29 0a 0a 28 64 65 66 75 6e 20 | izfcm.vector:max-norm))..(defun. |
2be0 | 66 20 28 78 29 0a 20 20 28 2f 20 28 2d 20 78 20 31 29 20 28 2b 20 78 20 31 29 29 29 0a 0a 28 64 | f.(x)...(/.(-.x.1).(+.x.1)))..(d |
2c00 | 65 66 75 6e 20 66 70 72 69 6d 65 20 28 78 29 0a 20 20 28 2f 20 32 20 28 65 78 70 74 20 28 2b 20 | efun.fprime.(x)...(/.2.(expt.(+. |
2c20 | 78 20 31 29 20 32 29 29 29 0a 0a 28 64 65 66 6d 61 63 72 6f 20 73 68 6f 77 63 61 73 65 20 28 73 | x.1).2)))..(defmacro.showcase.(s |
2c40 | 2d 65 78 70 72 29 0a 20 20 60 28 66 6f 72 6d 61 74 20 74 20 22 7e 61 20 3d 20 7e 61 7e 25 22 20 | -expr)...`(format.t."~a.=.~a~%". |
2c60 | 2c 28 66 6f 72 6d 61 74 20 6e 69 6c 20 22 7e 61 22 20 73 2d 65 78 70 72 29 20 2c 73 2d 65 78 70 | ,(format.nil."~a".s-expr).,s-exp |
2c80 | 72 29 29 0a 0a 28 64 65 66 75 6e 20 6d 61 69 6e 20 28 29 0a 20 20 28 73 68 6f 77 63 61 73 65 20 | r))..(defun.main.()...(showcase. |
2ca0 | 28 73 6d 61 63 65 70 73 29 29 0a 20 20 28 73 68 6f 77 63 61 73 65 20 28 64 6d 61 63 65 70 73 29 | (smaceps))...(showcase.(dmaceps) |
2cc0 | 29 0a 20 20 28 73 68 6f 77 63 61 73 65 20 28 6c 32 2d 6e 6f 72 6d 20 27 28 31 20 32 29 29 29 0a | )...(showcase.(l2-norm.'(1.2))). |
2ce0 | 20 20 28 73 68 6f 77 63 61 73 65 20 28 6c 31 2d 6e 6f 72 6d 20 27 28 31 20 32 29 29 29 0a 20 20 | ..(showcase.(l1-norm.'(1.2)))... |
2d00 | 28 73 68 6f 77 63 61 73 65 20 28 6c 69 6e 66 2d 6e 6f 72 6d 20 27 28 31 20 32 29 29 29 0a 20 20 | (showcase.(linf-norm.'(1.2)))... |
2d20 | 28 73 68 6f 77 63 61 73 65 20 28 6c 31 2d 64 69 73 74 61 6e 63 65 20 27 28 31 20 32 29 20 27 28 | (showcase.(l1-distance.'(1.2).'( |
2d40 | 2d 33 20 34 29 29 29 0a 20 20 28 73 68 6f 77 63 61 73 65 20 28 6c 32 2d 64 69 73 74 61 6e 63 65 | -3.4)))...(showcase.(l2-distance |
2d60 | 20 27 28 31 20 32 29 20 27 28 2d 33 20 34 29 29 29 0a 20 20 28 73 68 6f 77 63 61 73 65 20 28 6c | .'(1.2).'(-3.4)))...(showcase.(l |
2d80 | 69 6e 66 2d 64 69 73 74 61 6e 63 65 20 27 28 31 20 32 29 20 27 28 2d 33 20 34 29 29 29 0a 20 20 | inf-distance.'(1.2).'(-3.4)))... |
2da0 | 28 73 68 6f 77 63 61 73 65 20 28 6c 69 7a 66 63 6d 2e 76 65 63 74 6f 72 3a 6c 65 61 73 74 2d 73 | (showcase.(lizfcm.vector:least-s |
2dc0 | 71 75 61 72 65 73 2d 72 65 67 20 27 28 31 20 32 20 33 20 34 20 35 20 36 20 37 29 0a 20 20 20 20 | quares-reg.'(1.2.3.4.5.6.7)..... |
2de0 | 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 | ................................ |
2e00 | 20 20 20 20 20 20 20 20 20 27 28 30 2e 35 20 33 20 32 20 33 2e 35 20 35 20 36 20 37 2e 35 29 29 | .........'(0.5.3.2.3.5.5.6.7.5)) |
2e20 | 29 0a 20 20 28 73 68 6f 77 63 61 73 65 20 28 6c 69 7a 66 63 6d 2e 61 70 70 72 6f 78 3a 66 6f 72 | )...(showcase.(lizfcm.approx:for |
2e40 | 77 61 72 64 2d 64 65 72 69 76 61 74 69 76 65 2d 61 74 20 27 66 20 31 20 30 2e 30 30 30 30 31 29 | ward-derivative-at.'f.1.0.00001) |
2e60 | 29 0a 20 20 28 73 68 6f 77 63 61 73 65 20 28 6c 69 7a 66 63 6d 2e 61 70 70 72 6f 78 3a 63 65 6e | )...(showcase.(lizfcm.approx:cen |
2e80 | 74 72 61 6c 2d 64 65 72 69 76 61 74 69 76 65 2d 61 74 20 27 66 20 31 20 30 2e 30 30 30 30 31 29 | tral-derivative-at.'f.1.0.00001) |
2ea0 | 29 0a 20 20 28 73 68 6f 77 63 61 73 65 20 28 6c 69 7a 66 63 6d 2e 61 70 70 72 6f 78 3a 62 61 63 | )...(showcase.(lizfcm.approx:bac |
2ec0 | 6b 77 61 72 64 2d 64 65 72 69 76 61 74 69 76 65 2d 61 74 20 27 66 20 31 20 30 2e 30 30 30 30 31 | kward-derivative-at.'f.1.0.00001 |
2ee0 | 29 29 29 0a 0a 0a 0a 0a 23 31 2f 33 36 20 20 20 20 20 20 20 20 20 20 20 31 36 39 36 39 36 32 32 | ))).....#1/36...........16969622 |
2f00 | 37 37 20 20 35 30 31 20 20 20 32 30 20 20 20 20 31 30 30 36 34 34 20 20 36 38 20 20 20 20 20 20 | 77..501...20....100644..68...... |
2f20 | 20 20 60 0a 76 65 63 74 6f 72 2c 6c 65 61 73 74 2d 73 71 75 61 72 65 73 2d 72 65 67 2e 6c 69 73 | ..`.vector,least-squares-reg.lis |
2f40 | 70 00 00 00 00 00 00 00 28 69 6e 2d 70 61 63 6b 61 67 65 20 3a 6c 69 7a 66 63 6d 2e 76 65 63 74 | p.......(in-package.:lizfcm.vect |
2f60 | 6f 72 29 0a 0a 0a 0a 0a 23 31 2f 33 36 20 20 20 20 20 20 20 20 20 20 20 31 36 39 36 39 36 39 33 | or).....#1/36...........16969693 |
2f80 | 38 38 20 20 35 30 31 20 20 20 32 30 20 20 20 20 31 30 30 36 34 34 20 20 35 32 34 20 20 20 20 20 | 88..501...20....100644..524..... |
2fa0 | 20 20 60 0a 76 65 63 74 6f 72 2c 6c 65 61 73 74 2d 73 71 75 61 72 65 73 2e 6c 69 73 70 00 00 00 | ..`.vector,least-squares.lisp... |
2fc0 | 00 00 00 00 00 00 00 00 28 69 6e 2d 70 61 63 6b 61 67 65 20 3a 6c 69 7a 66 63 6d 2e 76 65 63 74 | ........(in-package.:lizfcm.vect |
2fe0 | 6f 72 29 0a 0a 28 64 65 66 75 6e 20 6c 65 61 73 74 2d 73 71 75 61 72 65 73 2d 72 65 67 20 28 78 | or)..(defun.least-squares-reg.(x |
3000 | 20 79 29 0a 20 20 28 6c 65 74 2a 20 28 28 6e 20 28 6c 65 6e 67 74 68 20 78 29 29 0a 20 20 20 20 | .y)...(let*.((n.(length.x))..... |
3020 | 20 20 20 20 20 28 73 75 6d 2d 79 20 28 72 65 64 75 63 65 20 23 27 2b 20 79 29 29 0a 20 20 20 20 | .....(sum-y.(reduce.#'+.y))..... |
3040 | 20 20 20 20 20 28 73 75 6d 2d 78 20 28 72 65 64 75 63 65 20 23 27 2b 20 78 29 29 0a 20 20 20 20 | .....(sum-x.(reduce.#'+.x))..... |
3060 | 20 20 20 20 20 28 73 75 6d 2d 78 79 20 28 72 65 64 75 63 65 20 23 27 2b 20 28 6d 61 70 63 61 72 | .....(sum-xy.(reduce.#'+.(mapcar |
3080 | 20 23 27 2a 20 78 20 79 29 29 29 0a 20 20 20 20 20 20 20 20 20 28 73 75 6d 2d 78 73 71 75 61 72 | .#'*.x.y)))..........(sum-xsquar |
30a0 | 65 64 20 28 72 65 64 75 63 65 20 23 27 2b 20 28 6d 61 70 63 61 72 20 23 27 2a 20 78 20 78 29 29 | ed.(reduce.#'+.(mapcar.#'*.x.x)) |
30c0 | 29 0a 20 20 20 20 20 20 20 20 20 28 62 20 28 2f 20 28 2d 20 28 2a 20 73 75 6d 2d 79 20 73 75 6d | )..........(b.(/.(-.(*.sum-y.sum |
30e0 | 2d 78 73 71 75 61 72 65 64 29 20 28 2a 20 73 75 6d 2d 78 20 73 75 6d 2d 78 79 29 29 0a 20 20 20 | -xsquared).(*.sum-x.sum-xy)).... |
3100 | 20 20 20 20 20 20 20 20 20 20 20 20 28 2d 20 28 2a 20 6e 20 73 75 6d 2d 78 73 71 75 61 72 65 64 | ............(-.(*.n.sum-xsquared |
3120 | 29 20 28 2a 20 73 75 6d 2d 78 20 73 75 6d 2d 78 29 29 29 29 0a 20 20 20 20 20 20 20 20 20 28 61 | ).(*.sum-x.sum-x))))..........(a |
3140 | 20 28 2f 20 28 2d 20 28 2a 20 6e 20 73 75 6d 2d 78 79 29 20 28 2a 20 73 75 6d 2d 78 20 73 75 6d | .(/.(-.(*.n.sum-xy).(*.sum-x.sum |
3160 | 2d 79 29 29 0a 20 20 20 20 20 20 20 20 20 20 20 20 20 20 20 28 2d 20 28 2a 20 6e 20 73 75 6d 2d | -y))................(-.(*.n.sum- |
3180 | 78 73 71 75 61 72 65 64 29 20 28 2a 20 73 75 6d 2d 78 20 73 75 6d 2d 78 29 29 29 29 29 0a 20 20 | xsquared).(*.sum-x.sum-x)))))... |
31a0 | 20 20 28 6c 69 73 74 20 61 20 62 29 29 29 0a 0a | ..(list.a.b))).. |