January 2013 - gsll-devel - mailman3.common-lisp.net

[GSLL-devel] Fixing up make-grid-data for clisp (was making copy-to working in clisp+grid: fix to element-type)
by Mirko Vukovic 26 Dec '13

26 Dec '13

This is a follow up on my posted fix to copy-to. For some reason this has triggered another error. I do not understand how that other error is triggered. The error and the fix are discussed at the bottom of the email On Fri, Jan 4, 2013 at 10:39 AM, Mirko Vukovic <mirko.vukovic(a)gmail.com>wrote: > The following works in sbcl+grid: > > (copy-to (make-array 2 :element-type 'double-float :initial-contents > '(1.0d0 2.0d0)) 'grid:foreign-array) > > But it does not work in clisp+grid. > > The root cause lies in the function grid/array.lisp/element-type. It uses > `(array-element-type grid)'. > > But this can present a problem. Quoting hyperspec: > > (Because of *array* <http://26_glo_a.htm#array> *upgrading*, this *type > specifier* <http://26_glo_t.htm#type_specifier> can in some cases denote > a *supertype* <http://26_glo_s.htm#supertype> of the *expressed array > element type* <http://26_glo_e.htm#expressed_array_element_type> of the * > array*.) > > In CLISP, array-element-type returns `T' when passed #(1d0 2d0) > > It returns T even when passed a simple array: > > (array-element-type (make-array 2 :element-type 'double-float > :initial-contents > '(1.0d0 2.0d0) > :adjustable nil > :fill-pointer nil > :displaced-to nil) ) > > The proposed fix is > > (defmethod element-type ((grid array)) > (type-of (row-major-aref grid 0)) > #+skip(array-element-type grid)) > > Now copy-to works in clisp as well. > > Mirko > > > After finishing this, I recompiled afresh, and got another error that I traced to make-grid-data and make-array In clisp, make-array will return an array of NIL's even if the element-type is specified as double float. In SBCL make-array will return an array filled with 1d0. I had to add some clisp specific code to make-grid-data in order for it to initialize properly. First, a helper function: (defun default-element (element-type) "Return a default element depending on element-type " (let ((a-list '((double-float . 1d0) (symbol . T)))) (let ((match (cdr (assoc element-type a-list)))) (assert match () "Default element type undefined for element-type ~a" element-type) match))) And second, a bit of set-up code in make-grid-data. I broke apart the (let* statement to provide an explicit initial element if it has not been specified already: (defmethod make-grid-data ((type (eql 'array)) dimensions rest-spec &key (initial-element nil initial-element-p) (initial-contents nil initial-contents-p)) (let ((element-type (top-spec-type rest-spec))) #+clisp (unless (or initial-element-p initial-contents-p) (setf initial-element (default-element element-type) initial-element-p t)) (let ((array ... UNCHANGED Summary: To fix to copy-to, I modified (defmethod element-type ((grid array)) For some unknown reason, that triggered an error in the array initialization routine. To fix that, I had to add code to make-grid-data, so that it fills double-float arrays with 0d0 instead of NIL's I will continue using this code to see if there are additional problems. BTW, I realize that CLISP is not a preferred platform for speedy computation, but that is what I use on Windows, until SBCL gets officially ported to it. And in addition, it is a nice cross-check on portability. Mirko

2 2

[GSLL-devel] making copy-to working in clisp+grid: fix to element-type
by Mirko Vukovic 26 Dec '13

26 Dec '13

The following works in sbcl+grid: (copy-to (make-array 2 :element-type 'double-float :initial-contents '(1.0d0 2.0d0)) 'grid:foreign-array) But it does not work in clisp+grid. The root cause lies in the function grid/array.lisp/element-type. It uses `(array-element-type grid)'. But this can present a problem. Quoting hyperspec: (Because of *array* <26_glo_a.htm#array> *upgrading*, this *type specifier*<26_glo_t.htm#type_specifier>can in some cases denote a *supertype* <26_glo_s.htm#supertype> of the *expressed array element type*<26_glo_e.htm#expressed_array_element_type>of the *array*.) In CLISP, array-element-type returns `T' when passed #(1d0 2d0) It returns T even when passed a simple array: (array-element-type (make-array 2 :element-type 'double-float :initial-contents '(1.0d0 2.0d0) :adjustable nil :fill-pointer nil :displaced-to nil) ) The proposed fix is (defmethod element-type ((grid array)) (type-of (row-major-aref grid 0)) #+skip(array-element-type grid)) Now copy-to works in clisp as well. Mirko

2 1

[GSLL-devel] added basic fft to my dataframes
by Mirko Vukovic 14 Jan '13

14 Jan '13

In preparation for experiments that start tomorrow, I added today real-array fft and power spectrum capabilities to my dataframe code. It uses GSLL's fft functions (search for file fft.lisp). The classes, generic functions, and methods are still a bit rough, and will be refined. You can read more about it: https://github.com/mirkov/data-table/blob/master/user/example1/README.org That same directory contains examples-w-fft.lisp that shows the code usage. I annotated the code. The png files show sunspot data and two power spectra. The example shows how I create one data-table with raw data, and then another one with slightly massaged data. I do the fft on the latter one. The example also shows that the current syntax is `functional' vs `declarative'. I explicitly build the second data table, instead of declaring its dependency on the first table. I think that the latter approach is preferable from a user's point of view. But I will refrain from trying to design and implement this declarative approach. I need more real-world experience before trying to do this. I may add windowing capability to the fft code later on. Mirko

1 0

[GSLL-devel] Error installing GSLL
by Patrick M. Groulx 02 Jan '13

02 Jan '13

Hi all, I have been trying to install GSLL for the past couple of days on Ubuntu and Debian, both running SBCL, and I keep encountering the same problem. I install Quicklisp without any difficulty but when I attempt to install GSLL I get the output that I have included below. It appears that the error is in static-vectors and the use of the component "CFFI-GROVEL-FILE." I'm not seeing where in the source files that I have the component defined anywhere, so I am at a loss as to what it does. I have included the output of (list-all-packages) at the very bottom in case I am missing something. I'm hoping that someone with a more in-depth understanding of the different components can help me out. Thanks, Patrick * (ql:quickload "gsll") ; Loading system definition from ; /usr/share/common-lisp/source/cl-cffi/cffi-grovel.asd into #<PACKAGE "ASDF2"> ; Registering #<SYSTEM "cffi-grovel"> ; Loading system definition from ; /usr/share/common-lisp/source/alexandria/alexandria.asd into ; #<PACKAGE "ASDF2"> ; Registering #<SYSTEM "alexandria"> ; Loading system definition from /usr/share/common-lisp/source/cl-cffi/cffi.asd ; into #<PACKAGE "ASDF2"> ; Registering #<SYSTEM "cffi"> ; Loading system definition from /usr/share/common-lisp/source/babel/babel.asd ; into #<PACKAGE "ASDF2"> ; Registering #<SYSTEM "babel"> ; Loading system definition from ; /usr/share/common-lisp/source/trivial-features/trivial-features.asd into ; #<PACKAGE "ASDF2"> ; Registering #<SYSTEM "trivial-features"> ; Loading system definition from ; /home/pgroulx/quicklisp/dists/quicklisp/software/asdf-system-connections-201 01006-darcs/asdf-system-connections.asd ; into #<PACKAGE "ASDF1"> ; Registering #<SYSTEM "asdf-system-connections"> To load "gsll": Load 1 ASDF system: gsll ; Loading "gsll" debugger invoked on a LOAD-SYSTEM-DEFINITION-ERROR in thread #<THREAD "initial thread" RUNNING {AAE5749}>: Error while trying to load definition for system antik from pathname /home/pgroulx/quicklisp/dists/quicklisp/software/antik-20120909-git/antik.as d: Error while trying to load definition for system static-vectors from pathname /home/pgroulx/quicklisp/dists/quicklisp/software/static-vectors-1.4/src/stat ic-vectors.asd: don't recognize component type CFFI-GROVEL-FILE Type HELP for debugger help, or (SB-EXT:QUIT) to exit from SBCL. restarts (invokable by number or by possibly-abbreviated name): 0: [REINITIALIZE-SOURCE-REGISTRY-AND-RETRY] Retry finding system static-vectors after reinitializing the source-registry. 1: Retry finding system antik after reinitializing the source-registry. 2: [ABORT ] Give up on "gsll" 3: Exit debugger, returning to top level. ((FLET #:LAMBDA2060) #<LOAD-SYSTEM-DEFINITION-ERROR {B6E6AD9}>) 0] 0 .......; ; compilation aborted because of fatal error: ; SB-INT:SIMPLE-READER-PACKAGE-ERROR at 4208 (line 102, column 38) on #<SB-SYS:FD-STREAM ; for "file /home/pgroulx/quicklisp/dists/quicklisp/software/antik-20120909-git/foreign- array/foreign-array.lisp" ; {C1CA541}>: ; package "STATIC-VECTORS" not found ; debugger invoked on a ASDF:COMPILE-ERROR in thread #<THREAD "initial thread" RUNNING {AAE5749}>: Error while invoking #<COMPILE-OP (:VERBOSE NIL) {B53E519}> on #<CL-SOURCE-FILE "antik" "foreign-array" "foreign-array"> Type HELP for debugger help, or (SB-EXT:QUIT) to exit from SBCL. restarts (invokable by number or by possibly-abbreviated name): 0: [RETRY ] Retry compiling #<CL-SOURCE-FILE "antik" "foreign-array" "foreign-array">. 1: [ACCEPT] Continue, treating compiling #<CL-SOURCE-FILE "antik" "foreign-array" "foreign-array"> as having been successful. 2: [ABORT ] Give up on "gsll" 3: Exit debugger, returning to top level. ((SB-PCL::FAST-METHOD ASDF:PERFORM (ASDF:COMPILE-OP ASDF:CL-SOURCE-FILE)) #<unavailable argument> #<unavailable argument> #<ASDF:COMPILE-OP (:VERBOSE NIL) {B53E519}> #<ASDF:CL-SOURCE-FILE "antik" "foreign-array" "foreign-array">) 0] 2 ; ; compilation unit aborted ; caught 2 fatal ERROR conditions ("gsll") (#<PACKAGE "AFFI"> #<PACKAGE "ANTIK-USER"> #<PACKAGE "GRID"> #<PACKAGE "ANTIK"> #<PACKAGE "ITERATE"> #<PACKAGE "SPLIT-SEQUENCE"> #<PACKAGE "FOREIGN-STRUCTURES-BY-VALUE"> #<PACKAGE "DRAKMA"> #<PACKAGE "PURI"> #<PACKAGE "CL-BASE64"> #<PACKAGE "CHUNGA"> #<PACKAGE "CL-PPCRE"> #<PACKAGE "USOCKET"> #<PACKAGE "CL+SSL"> #<PACKAGE "FLEXI-STREAMS"> #<PACKAGE "TRIVIAL-GRAY-STREAMS"> #<PACKAGE "BORDEAUX-THREADS"> #<PACKAGE "PURI-SYSTEM"> #<PACKAGE "CL-BASE64-SYSTEM"> #<PACKAGE "CL-PPCRE-ASD"> #<PACKAGE "TRIVIAL-GRAY-STREAMS-SYSTEM"> #<PACKAGE "FLEXI-STREAMS-SYSTEM"> #<PACKAGE "CL+SSL-SYSTEM"> #<PACKAGE "ASDF-SYSTEM-CONNECTIONS-SYSTEM"> #<PACKAGE "CFFI-GROVEL"> #<PACKAGE "CFFI-FEATURES"> #<PACKAGE "CFFI"> #<PACKAGE "CFFI-SYS"> #<PACKAGE "BABEL"> #<PACKAGE "BABEL-ENCODINGS"> #<PACKAGE "ALEXANDRIA.0.DEV"> #<PACKAGE "DRAKMA-ASD"> #<PACKAGE "TRIVIAL-GARBAGE"> #<PACKAGE "QL-SBCL"> #<PACKAGE "SB-BSD-SOCKETS"> #<PACKAGE "SB-BSD-SOCKETS-SYSTEM"> #<PACKAGE "SB-BSD-SOCKETS-INTERNAL"> #<PACKAGE "SB-POSIX-SYSTEM"> #<PACKAGE "SB-GROVEL-SYSTEM"> #<PACKAGE "SB-GROVEL"> #<PACKAGE "SB-POSIX"> #<PACKAGE "QL-MKCL"> #<PACKAGE "QL-ECL"> #<PACKAGE "QL-LISPWORKS"> #<PACKAGE "QL-SCL"> #<PACKAGE "QL-CMUCL"> #<PACKAGE "QL-CLISP"> #<PACKAGE "QL-CCL"> #<PACKAGE "QL-ABCL"> #<PACKAGE "QL-ALLEGRO"> #<PACKAGE "QUICKLISP-CLIENT"> #<PACKAGE "QL-DIST-USER"> #<PACKAGE "QL-DIST"> #<PACKAGE "QL-CDB"> #<PACKAGE "QL-GUNZIPPER"> #<PACKAGE "QL-MINITAR"> #<PACKAGE "QL-HTTP"> #<PACKAGE "QL-PROGRESS"> #<PACKAGE "QL-NETWORK"> #<PACKAGE "QL-IMPL-UTIL"> #<PACKAGE "QL-IMPL"> #<PACKAGE "QL-CONFIG"> #<PACKAGE "QL-SETUP"> #<PACKAGE "QL-UTIL"> #<PACKAGE "ASDF"> #<PACKAGE "COMMON-LISP-CONTROLLER"> #<PACKAGE "COMMON-LISP-USER"> #<PACKAGE "SB-LOOP"> #<PACKAGE "SB-EVAL"> #<PACKAGE "SB-WALKER"> #<PACKAGE "SB-UNIX"> #<PACKAGE "SB-SEQUENCE"> #<PACKAGE "SB-PROFILE"> #<PACKAGE "SB-PRETTY"> #<PACKAGE "SB-MOP"> #<PACKAGE "SB-THREAD"> #<PACKAGE "SB-GRAY"> #<PACKAGE "SB-FORMAT"> #<PACKAGE "SB-FASL"> #<PACKAGE "SB-DISASSEM"> #<PACKAGE "SB-DEBUG"> #<PACKAGE "SB-C"> #<PACKAGE "SB-BIGNUM"> #<PACKAGE "SB-ASSEM"> #<PACKAGE "SB-ALIEN"> #<PACKAGE "SB-PCL"> #<PACKAGE "SB-ALIEN-INTERNALS"> #<PACKAGE "SB-DI"> #<PACKAGE "KEYWORD"> #<PACKAGE "SB-IMPL"> #<PACKAGE "SB-SYS"> #<PACKAGE "SB-KERNEL"> #<PACKAGE "SB-VM"> #<PACKAGE "SB-INT"> #<PACKAGE "SB-EXT"> #<PACKAGE "COMMON-LISP">)

1 0