July 2010 - gsll-devel - mailman3.common-lisp.net

[Gsll-devel] patch for calculus/monte-carlo.lisp
by Sam Steingold 12 Aug '10

12 Aug '10

PI has type long-float, not double-float, so the following patch is necessary: ========================================================= diff --git a/calculus/monte-carlo.lisp b/calculus/monte-carlo.lisp index cb18bb1..32bae65 100644 --- a/calculus/monte-carlo.lisp +++ b/calculus/monte-carlo.lisp @@ -211,7 +211,8 @@ (defparameter *mc-lower* #m(0.0d0 0.0d0 0.0d0)) (defparameter *mc-upper* - (grid:make-foreign-array 'double-float :initial-contents (list pi pi pi))) + (grid:make-foreign-array 'double-float :initial-contents + (list #1=#.(float pi 0d0) #1# #1#))) (defun random-walk-plain-example (&optional (nsamples 500000)) (monte-carlo-integrate-plain 'mcrw *mc-lower* *mc-upper* nsamples)) ========================================================= alternatively, you may need to call COERCE in grid:make-foreign-array or something. -- Sam Steingold <http://sds.podval.org>

2 5

Re: [Gsll-devel] status.text & clisp
by Sam Steingold 12 Aug '10

12 Aug '10

Hi, status.text still needs updating: CLISP 2.49 on x86_64 (64-bit build): TOTAL: 1350 assertions passed, 21 failed, 13 execution errors. I don't know what "execution errors" are (probably things like FUNCALL: undefined function NIL which I don't know how to reproduce), but I would appreciate bug reports. Some assertion errors should probably not be there: Expected (0.9999999999999997d0 2.0d0 30.0d0) but saw (0.9999999999999999d0 2.0d0 30.0d0) (also, it would be nice to align the expected and actual output). thanks. Sam. On 8/19/09, Sam Steingold <sds(a)gnu.org> wrote: > Hi, > > http://common-lisp.net/project/gsll/status.text says: > > 64-bit CLISP: not tested; bug in CLISP causes crash on loading > > I do not have this problem with 64-bit clisp: > TOTAL: 1345 assertions passed, 26 failed, 13 execution errors. > (BTW, how do I debug those 13 errors? i.e., how do I list them, rerun > them, see the code which triggered them?) > > also, status.text says > > CHEBYSHEV: Should have signalled TYPE-ERROR but saw #<SIMPLE-ERROR #x20EB4F1E> > > - how do I reproduce that? > > FFI::FOREIGN-CALL-OUT: 3.1415926535897932385L0 cannot be converted to > the foreign type DOUBLE-FLOAT > > - pi is a long float, you have to convert it to a double explicitly, > like (float pi 0d0). > > thanks. -- Sam Steingold <http://sds.podval.org>

2 7

[Gsll-devel] GSLL failed with CCL Trunk
by Leo 07 Aug '10

07 Aug '10

Hello, I am trying to run the latest GSLL to make sure no crash in CCL. Unfortunately somehow I could not load GSLL at all. My GSLL repo is update-to-date. I'm getting an error (as attached below) and file ASDF-TMP-mathematical.dx64fsl is empty (i.e. 0B). Any ideas? Thank you. Leo -------------------------------- > Error: DIVISION-BY-ZERO detected > performing LOG on (1.0) > While executing: CCL::%FP-ERROR-FROM-STATUS, in process listener(1). > Type :GO to continue, :POP to abort, :R for a list of available restarts. > If continued: Skip loading init file. > Type :? for other options. 1 > :b (1A18BC8) : 0 (%FP-ERROR-FROM-STATUS 6404 LOG 1.0 NIL) 317 (1A18C08) : 1 (%SINGLE-FLOAT-LOG 1.0) 205 (1A18C28) : 2 (LOG 1 2) 373 (1A18C60) : 3 (INIT-STREAM-IOBLOCK #<BASIC-FILE-BINARY-OUTPUT-STREAM (NIL/:closed #x30200101128D> :INSIZE NIL :OUTSIZE 4096 :SHARE-BUFFERS-P NIL :ELEMENT-TYPE (UNSIGNED-BYTE 8) :DEVICE 6 :ADVANCE-FUNCTION NIL :LISTEN-FUNCTION CCL::FD-STREAM-LISTEN :EOFP-FUNCTION NIL :FORCE-OUTPUT-FUNCTION CCL::OUTPUT-FILE-FORCE-OUTPUT :CLOSE-FUNCTION CCL::FD-STREAM-CLOSE :ELEMENT-SHIFT NIL :INTERACTIVE NIL :SHARING :PRIVATE :CHARACTER-P NIL :ENCODING NIL :LINE-TERMINATION NIL :INPUT-TIMEOUT NIL :OUTPUT-TIMEOUT NIL :DEADLINE NIL) 2861 (1A18DD8) : 4 (MAKE-IOBLOCK-STREAM #<BUILT-IN-CLASS CCL::BASIC-FILE-BINARY-OUTPUT-STREAM> :INSIZE NIL :OUTSIZE 4096 :SHARE-BUFFERS-P NIL :INTERACTIVE NIL :DIRECTION :OUTPUT :ELEMENT-TYPE (INTEGER 0 255) :DIRECTION :OUTPUT :LISTEN-FUNCTION CCL::FD-STREAM-LISTEN :CLOSE-FUNCTION CCL::FD-STREAM-CLOSE :ADVANCE-FUNCTION NIL :FORCE-OUTPUT-FUNCTION CCL::OUTPUT-FILE-FORCE-OUTPUT :DEVICE 6 :ENCODING NIL :EXTERNAL-FORMAT :BINARY :SHARING :PRIVATE :LINE-TERMINATION NIL :CHARACTER-P NIL) 381 (1A18E10) : 5 (MAKE-FILE-STREAM "/Users/Leo/.cache/common-lisp/clozure-1.6-fasl94-darwin-x86-64/Users/Leo/.ccl/site/gsll/mathematical/ASDF-TMP-mathematical.dx64fsl" :OUTPUT (INTEGER 0 255) :SUPERSEDE :CREATE CCL::BASIC-FILE-STREAM :DEFAULT :PRIVATE T) 2709 (1A18F18) : 6 (OPEN "/Users/Leo/.cache/common-lisp/clozure-1.6-fasl94-darwin-x86-64/Users/Leo/.ccl/site/gsll/mathematical/ASDF-TMP-mathematical.dx64fsl" :DIRECTION :OUTPUT :ELEMENT-TYPE (UNSIGNED-BYTE 8) :IF-EXISTS :SUPERSEDE :IF-DOES-NOT-EXIST :CREATE :EXTERNAL-FORMAT :DEFAULT :CLASS FILE-STREAM :SHARING :PRIVATE :BASIC T) 797 (1A18FF0) : 7 (FASL-DUMP-FILE NIL NIL ((23 83) (47 "home:.ccl;site;gsll;mathematical;mathematical.lisp.newest") (70 #<# *#:#-# #>) (4 #) (70 #<# *#:#-# #>) ...) #<HASH-TABLE :TEST EQ size 55/567 #x30200101E61D> "/Users/Leo/.cache/common-lisp/clozure-1.6-fasl94-darwin-x86-64/Users/Leo/.ccl/site/gsll/mathematical/ASDF-TMP-mathematical.dx64fsl") 373 (1A190C0) : 8 (FASL-SCAN-FORMS-AND-DUMP-FILE ((23 83) (47 "home:.ccl;site;gsll;mathematical;mathematical.lisp.newest") (70 #<# *#:#-# #>) (4 #) (70 #<# *#:#-# #>) ...) "/Users/Leo/.cache/common-lisp/clozure-1.6-fasl94-darwin-x86-64/Users/Leo/.ccl/site/gsll/mathematical/ASDF-TMP-mathematical.dx64fsl" #<CCL::LEXICAL-ENVIRONMENT #x302000EFE39D>) 349 (1A19178) : 9 (%COMPILE-FILE "/Users/Leo/.ccl/site/gsll/mathematical/mathematical.lisp" "/Users/Leo/.cache/common-lisp/clozure-1.6-fasl94-darwin-x86-64/Users/Leo/.ccl/site/gsll/mathematical/ASDF-TMP-mathematical.dx64fsl" NIL NIL NIL T T NIL T :DEFER NIL #<BACKEND DARWINX8664 #x3020001BE85D> :DEFAULT NIL 0) 3101 (1A19448) : 10 (COMPILE-FILE #P"/Users/Leo/.ccl/site/gsll/mathematical/mathematical.lisp" :OUTPUT-FILE #P"/Users/Leo/.cache/common-lisp/clozure-1.6-fasl94-darwin-x86-64/Users/Leo/.ccl/site/gsll/mathematical/ASDF-TMP-mathematical.dx64fsl" :VERBOSE NIL :PRINT NIL :LOAD NIL :FEATURES NIL :TARGET :DARWINX8664 :SAVE-LOCAL-SYMBOLS T :SAVE-DOC-STRINGS T :SAVE-DEFINITIONS NIL :SAVE-SOURCE-LOCATIONS T :EXTERNAL-FORMAT :DEFAULT :FORCE NIL :COMPILE-FILE-ORIGINAL-TRUENAME NIL :COMPILE-FILE-ORIGINAL-BUFFER-OFFSET 0 :BREAK-ON-PROGRAM-ERRORS :DEFER) 1317 (1A19590) : 11 (COMPILE-FILE* #P"/Users/Leo/.ccl/site/gsll/mathematical/mathematical.lisp" :OUTPUT-FILE #P"/Users/Leo/.cache/common-lisp/clozure-1.6-fasl94-darwin-x86-64/Users/Leo/.ccl/site/gsll/mathematical/mathematical.dx64fsl") 309 (1A195E0) : 12 (FUNCALL #'#<#<STANDARD-METHOD ASDF:PERFORM (ASDF:COMPILE-OP ASDF:CL-SOURCE-FILE)>> #<COMPILE-OP NIL #x302000BC64BD> #<CL-SOURCE-FILE "gsll" "mathematical" "mathematical">) 357 (1A19660) : 13 (%%BEFORE-AND-AFTER-COMBINED-METHOD-DCODE (NIL #<STANDARD-METHOD ASDF:PERFORM #> . 3420901)) 885 (1A196D0) : 14 (%%STANDARD-COMBINED-METHOD-DCODE ((#<#>) (#<#>) #<STANDARD-METHOD ASDF:PERFORM #>) 3420901) 237 (1A19740) : 15 (%CALL-NEXT-METHOD (NIL #<STANDARD-METHOD ASDF::PERFORM-WITH-RESTARTS #> . 3420950)) 1021 (1A197D8) : 16 (FUNCALL #'#<#<STANDARD-METHOD ASDF::PERFORM-WITH-RESTARTS (ASDF:COMPILE-OP ASDF:CL-SOURCE-FILE)>> #<COMPILE-OP NIL #x302000BC64BD> #<CL-SOURCE-FILE "gsll" "mathematical" "mathematical">) 637 (1A19848) : 17 (%%STANDARD-COMBINED-METHOD-DCODE (#<STANDARD-METHOD ASDF::PERFORM-WITH-RESTARTS #> #<STANDARD-METHOD ASDF::PERFORM-WITH-RESTARTS #>) 3420950) 709 (1A198C8) : 18 (FUNCALL #'#<(:INTERNAL CCL::WITH-COMPILATION-UNIT-BODY (ASDF:OPERATE (T T)))>) 917 (1A19928) : 19 (CALL-WITH-COMPILATION-UNIT #<COMPILED-LEXICAL-CLOSURE (:INTERNAL CCL::WITH-COMPILATION-UNIT-BODY #) #x70F16F> :OVERRIDE NIL) 197 (1A19970) : 20 (FUNCALL #'#<#<STANDARD-METHOD ASDF:OPERATE (T T)>> ASDF:LOAD-OP :QDSOA) 813 (1A19A50) : 21 (CALL-CHECK-REGS ASDF:LOAD-SYSTEM :QDSOA) 229 (1A19A88) : 22 (FUNCALL #'#<(:INTERNAL CCL::WITH-COMPILATION-UNIT-BODY CCL::LOAD-FROM-STREAM)>) 757 (1A19B10) : 23 (CALL-WITH-COMPILATION-UNIT #<COMPILED-LEXICAL-CLOSURE (:INTERNAL CCL::WITH-COMPILATION-UNIT-BODY CCL::LOAD-FROM-STREAM) #x70F22F> :OVERRIDE NIL) 197 (1A19B58) : 24 (LOAD-FROM-STREAM #<BASIC-FILE-CHARACTER-INPUT-STREAM ("/Users/Leo/\\.ccl-init.lisp"/4 ISO-8859-1) #x30200053311D> NIL) 357 (1A19B98) : 25 (%LOAD #P"/Users/Leo/\\.ccl-init.lisp" NIL NIL NIL :DEFAULT) 3389 (1A19CE8) : 26 (LOAD "home:\\.ccl-init" :VERBOSE NIL :PRINT NIL :IF-DOES-NOT-EXIST NIL :EXTERNAL-FORMAT :DEFAULT) 989 (1A19D80) : 27 (STARTUP-CCL ("home:ccl-init" "home:\\.ccl-init")) 749 (1A19E00) : 28 (FUNCALL #'#<(:INTERNAL (TOPLEVEL-FUNCTION (CCL::LISP-DEVELOPMENT-SYSTEM T)))>) 77 (1A19E20) : 29 (FUNCALL #'#<(:INTERNAL CCL::MAKE-MCL-LISTENER-PROCESS)>) 645 (1A19EB8) : 30 (RUN-PROCESS-INITIAL-FORM #<TTY-LISTENER listener(1) [Active] #x30200052F48D> (#<COMPILED-LEXICAL-CLOSURE # #x30200052EEAF>)) 685 (1A19F48) : 31 (FUNCALL #'#<(:INTERNAL (CCL::%PROCESS-PRESET-INTERNAL (PROCESS)))> #<TTY-LISTENER listener(1) [Active] #x30200052F48D> (#<COMPILED-LEXICAL-CLOSURE # #x30200052EEAF>)) 389 (1A19F98) : 32 (FUNCALL #'#<(:INTERNAL CCL::THREAD-MAKE-STARTUP-FUNCTION)>) 301 1 > -- Any Emacs contains an ad hoc, informally-specified, bug-ridden, slow implementation of half of Common Lisp.

2 5

[Gsll-devel] Status of status.text
by Jason Nielsen 21 Jul '10

21 Jul '10

Dear Liam, I just pulled the latest version of gsll. Apparently you have busy as there are quite a few changes. That is great, keep up the good work! With all the recent e-mail on the test suite I'm wondering how up to date is the status.text file? When I run the test suite on an x86-64 box running Ubuntu 10.04 with SBCL 1.0.39 and a custom built gsl library (mostly some CFLAGS tweaks) I get the following two failures: EXPONENTIAL: (MULTIPLE-VALUE-LIST (LET ((RNG (MAKE-RANDOM-NUMBER-GENERATOR +MT19937+ 0))) (LOOP FOR I FROM 0 TO 10 COLLECT (SAMPLE RNG :EXPONENTIAL :MU 10.0d0)))) failed: Expected ((0.0025828444588394794d0 18.145581427987647d0 12.636598054339759d0 0.5424387252062355d0 14.624994234158105d0 7.236607929535993d0 0.4345362449683603d0 2.95303920904529d0 6.161052939065796d0 3.011686333539114d0 2.7451079819355364d0)) but saw ((82.61578216370394d0 1.77823538531874d0 3.3214653339021365d0 29.412645812312775d0 2.6351843583734817d0 6.635371453335953d0 31.57709150859036d0 13.637712930094189d0 7.766237034468241d0 13.468928238240315d0 14.268822088229632d0)) EXPONENTIAL: 11 assertions passed, 1 failed. and: LU: (MULTIPLE-VALUE-LIST (TEST-LU-SOLVE-DIM (CREATE-COMPLEX-MATRIX 7))) failed: Expected (#2m(#C(24.0717272023734d0 -9.84612797621247d0) #C(-269.338853034031d0 87.5455232472528d0) #C(2966.61356736296d0 -1026.24473923993d0) #C(-18207.3812124749d0 5673.8447304241d0) #C(55769.3879019068d0 -16154.0963210502d0) #C(-78894.1207561151d0 19505.3812987858d0) #C(39554.8551241728d0 -7765.93696255317d0))) but saw (#2m(#C(-81.70598196064235d0 -81.70599595009693d0) #C(3354.209927435102d0 3354.278264591346d0) #C(-31954.85554225817d0 -31956.353685805567d0) #C(119027.61090077131d0 119037.36678048137d0) #C(-203627.71695137964d0 -203653.92708809028d0) #C(160626.5466779804d0 160657.2640539496d0) #C(-47344.00882810219d0 -47357.00253607087d0))) LU: 11 assertions passed, 1 failed. TOTAL: 4014 assertions passed, 2 failed, 0 execution errors. a few quick tests of the LU decomp routine seem to suggest that it is working fine so I'm wondering if this is the kind of "numeric fuzz" you were alluding to in a previous post. With regards, Jason

6 11

[Gsll-devel] Unhandled exception while running (lisp-unit:run-tests)
by Leo 21 Jul '10

21 Jul '10

Hello all, I am running ccl trunk (checked out on 2010-06-13) with GSLL (gsl 1.4) on Linux: openSUSE 11.0 (X86-64) VERSION = 11.0 The machine has 4 CPUs. (asdf:load-system :gsll) (asdf:load-system :gsll-tests) (in-package :gsl) (lisp-unit:run-tests) and get this: ...... LU: (MULTIPLE-VALUE-LIST (TEST-LU-SOLVE-DIM (CREATE-COMPLEX-MATRIX 7))) failed: Expected (#<VECTOR-COMPLEX-DOUBLE-FLOAT #(#C(24.0717272023734D0 -9.84612797621247D0) #C(-269.338853034031D0 87.5455232472528D0) #C(2966.61356736296D0 -1026.24473923993D0) #C(-18207.3812124749D0 5673.8447304241D0) #C(55769.3879019068D0 -16154.0963210502D0) #C(-78894.1207561151D0 19505.3812987858D0) #C(39554.8551241728D0 -7765.93696255317D0))>) but saw (#<VECTOR-COMPLEX-DOUBLE-FLOAT #(#C(-81.70598196064235D0 -81.70599595009693D0) #C(3354.209927435102D0 3354.278264591346D0) #C(-31954.85554225817D0 -31956.353685805567D0) #C(119027.61090077131D0 119037.36678048137D0) #C(-203627.71695137964D0 -203653.92708809028D0) #C(160626.5466779804D0 160657.2640539496D0) #C(-47344.00882810219D0 -47357.00253607087D0))>) LU: 11 assertions passed, 1 failed. MATHEMATICAL: 3 assertions passed, 0 failed. MATHIEU: 75 assertions passed, 0 failed. MATRIX-ADD: 10 assertions passed, 0 failed. MATRIX-COPY: 12 assertions passed, 0 failed. MATRIX-COPY-TO-CL-AND-BACK: 12 assertions passed, 0 failed. MATRIX-DIV: 10 assertions passed, 0 failed. MATRIX-MAX: 10 assertions passed, 0 failed. MATRIX-MAX-INDEX: 10 assertions passed, 0 failed. MATRIX-MEAN: 10 assertions passed, 0 failed. MATRIX-MIN: 10 assertions passed, 0 failed. MATRIX-MIN-INDEX: 10 assertions passed, 0 failed. MATRIX-MINMAX: 10 assertions passed, 0 failed. MATRIX-MINMAX-INDEX: 10 assertions passed, 0 failed. MATRIX-MULT: 10 assertions passed, 0 failed. MATRIX-PRODUCT: 8 assertions passed, 0 failed. MATRIX-PRODUCT-HERMITIAN: 4 assertions passed, 0 failed. MATRIX-PRODUCT-NONSQUARE: 4 assertions passed, 0 failed. MATRIX-PRODUCT-SYMMETRIC: 4 assertions passed, 0 failed. MATRIX-PRODUCT-TRIANGULAR: 8 assertions passed, 0 failed. MATRIX-SET-ALL: 12 assertions passUnhandled exception 11 at 0x7f40eb6dd1c8, context->regs at #x7fff3fc87178 Exception occurred while executing foreign code received signal 11; faulting address: (nil) unexpected si_code value: 128 ? for help [18645] Clozure CL kernel debugger: B current thread: tcr = 0x7f40ec2146f0, native thread ID = 0x48d5, interrupts enabled (#x00007F40EB664DF0) #x00003000000767BC : #<Function %DO-FF-CALL #x000030000007676F> + 77 (#x00007F40EB664E00) #x0000300000077774 : #<Function %FF-CALL #x00003000000768BF> + 3765 (#x00007F40EB664E90) #x000030200101EA74 : #<Function (:INTERNAL (INITIALIZE-INSTANCE :AFTER (MATHIEU))) #x000030200101E9FF> + 117 (#x00007F40EB664EB8) #x00003000000E2A8C : #<Function DRAIN-TERMINATION-QUEUE #x00003000000E288F> + 509 (#x00007F40EB664F40) #x0000300000A73AAC : #<Function (:INTERNAL ADD-GC-HOOK) #x0000300000A73A4F> + 93 (#x00007F40EB664F60) #x000030000045AD84 : #<Function HOUSEKEEPING #x000030000045AD4F> + 53 (#x00007F40EB664F90) #x000030000045B57C : #<Function HOUSEKEEPING-LOOP #x000030000045B36F> + 525 [18645] Clozure CL kernel debugger: Leo

2 8

[Gsll-devel] grid: map-n-grids example
by Mirko Vukovic 20 Jul '10

20 Jul '10

Hello, I just spent a few minutes on figuring out how to use map-n-grids as there are no examples in the distribution or documentation. Here is one example, here is one that accomplishes the following (in pseudo-code): (map 'grid #'(lambda (x1 x2) (- x1 x2)) grid1 grid2) (map-n-grids :sources `((,*array-3-4-double-float* ,(affi *array-3-4-double-float*)) (,*array-3-4-double-float* ,(affi *array-3-4-double-float*))) :combination-function #'(lambda (i x1 x2) (print i) (- x1 x2))) In this example, I am not sure of the role of the first variable (I did not have more time to stare at map.list) Mirko

2 1

[Gsll-devel] New GSD array organization
by Liam Healy 19 Jul '10

19 Jul '10

GSLL's underlying system for generalized array, the "Grid Structured Data" (gsd) has been reorganized. There are some user-visible changes: 1. It is simpler to exchange arrays between foreign libraries, and lisp code. For example, a foreign array may be created in Lisp, passed to GSLL, then passed to another foreign library, then processed once again in Lisp. 2. Grid operations include common array operations for mathematical use, such as row, column, diagonal, transpose, slice and others; these are built on the more general map-grid which in turn is built on map-n-grids, so that more complicated array operations can be defined generically. 3. The function #'grid:copy will copy an array, possibly with a change of grid type ('array to 'foreign-array or vice versa), or other property. 4. Grids types possible at present are CL arrays and foreign arrays. 5. On SBCL, ECL, CCL, Allegro, and Lispworks, foreign-arrays are optionally built using static-vectors (http://gitorious.org/iolib/static-vectors) For foreign arrays made in Lisp (not by foreign libraries), the function cl-array is the same array in Lisp. Note that it is not a copy; the contents of memory are the same. 6. Changes 1. The grid and foreign-array systems both use the grid package; the c-array package doesn't exist any more. 2. All the array-related definitions of GSLL are now in the system foreign-array, including the individual classes such as grid:vector-double-float. To access elements, use grid:gref (gsl:maref has been removed). 3. Internally, foreign-array objects have a slot foreign-metadata. Each application is free to add properties to this list; GSLL adds 'gsll:mpointer and 'gsll:block-data. You should do a fresh pull of both GSLL and GSD. Liam

1 0