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

[GSLL-devel] grid, foreign-array, fsbv and libffi
by Sumant Oemrawsingh 06 Dec '10

06 Dec '10

Hi Liam, Recently, I tried getting foreign-array working on windows (with sbcl). I noticed that foreign-array uses fsbv (apparently) and fsbv uses libffi. I haven't dug around in the code yet, but I was wondering if you have a quick explanation why libffi is used, and why cffi is not sufficient. The reason I ask is because I'm wondering if it would be possible to hack it under windows so it doesn't need libffi (which apparently is quite difficult to get working there). Thanks, Sumant

3 4

Re: [GSLL-devel] Patch to fsbf which makes installing gsll on fedora work
by Joshua Holbrook 30 Nov '10

30 Nov '10

Sumant said: > I don't think this should be added to fsbv in this hardcoded fashion, but should > rather be determined at (or before) grovel time This sounds like a better approach to me as well, especially considering that the installed version number could (for all I know) change with any given software update. Someone else will have to figure it out, though, since I'm pretty clueless about these things. --Josh

1 0

[GSLL-devel] Patch to fsbf which makes installing gsll on fedora work
by Joshua Holbrook 29 Nov '10

29 Nov '10

I don't know if this is the right place or not, but at any rate, the patch to the fsvb repository that was discussed on irc a week (or so) ago is attached. Basically, I just added the location of ffi.h for Fedora 12 to libffi-unix.lisp . Chances are it's not ready-for-prime-time since I'm not very good with lisp yet and have only really been able to test it on my machines, but hopefully it's useful. Thanks, --Joshua Holbrook

2 1

[GSLL-devel] Bug or feature: grid-map and complex type destinations
by Mirko Vukovic 23 Nov '10

23 Nov '10

Hi, I am reporting this in case it is a gsd bug (it must be a bug somewhere: in my head, my setup, or gsd) Mapping the `sqrt' function over a vector of negative numbers as follows: grid> (map-grid :source #1m(-1d0 -2d0 -3d0) :element-function #'sqrt :destination-specification '((foreign-array 3) complex-double-float)) gives an error: Element type COMPLEX-DOUBLE-FLOAT not permitted. [Condition of type SIMPLE-ERROR] I traced it to a `make-instance' statement. I also tried to instantiate a vector-complex-double-float object grid> (make-instance 'vector-complex-double-float :dimensions '(3) :element-type 'complex-double-float) that resulted in the same error message. Mirko

2 5

[GSLL-devel] *row-separator* error with cmucl
by levent guel 23 Nov '10

23 Nov '10

i got the *row-separator* error with cmucl, upon changing the (export *row-separator*) thing to (export '*row-separator*) the error disappeared and gsll is loading now. (library gsd, file foreign-arrays.lisp, in package grid) greetings Levent Guel levgue(a)gmx.net

2 1

[Gsll-devel] Allegro CL 8.2 on 10.5?
by Vipul Periwal 19 Nov '10

19 Nov '10

Has anyone succeeded in building GSLL for Allegro CL 8.2 on 10.5? I'm asking because I've been fixing a bunch of things to get pieces built one by one and it occurs to me that there is probably something wrong with what I'm doing. The last spot was getting through the libgsl- unix.lisp compilation where the appropriate compiler flag (-I/opt/ local/include/) was not being passed to cc. And then the resulting output has lower-case +edom+ etc. gslerrorno's while conditions.lisp is expecting +EDOM+ etc. I'm using mlisp so this last would probably not be a problem with alisp, but I hope someone can tell me of a workaround, or tell me to give up. Thanks. Vipul Periwal

2 1

[GSLL-devel] fsbv compilation error on cygwin: libffi-unix.c:54: error: 'FFI_UNIX64' undeclared
by Mirko Vukovic 18 Nov '10

18 Nov '10

Hello, I am trying to install fsbv on cygwin (windows xp) I downloaded the libffi package for the cygffi-4.dll and the gcc-4-core for ffi.h. I had to modify the library linking in init.lisp as follows (untested on my redhat but works on cygwin): (cffi:define-foreign-library libffi (:cygwin "cygffi-4.dll") (:unix "libffi.so") (t (:default "libffi.so"))) #| (cffi:load-foreign-library #+darwin "libffi.dylib" #+(and (not darwin) unix) "libffi.so") |# (cffi:use-foreign-library libffi) However, the compilation breaks with the error message: /home/external-libraries/common-lisp/fsbv/libffi-unix.c: In function 'main': /home/external-libraries/common-lisp/fsbv/libffi-unix.c:54: error: 'FFI_UNIX64' undeclared (first use in this function). Any ideas? Mirko

1 0

[GSLL-devel] gsd/grid: cannot convert #(...doubles...) to #m(...doubles...)
by Mirko Vukovic 16 Nov '10

16 Nov '10

I can go from foreign arrays to cl-arrays: (grid:copy-to #m(1.0d0 2.0d0) 'grid::array) (although the use of the *internal* symbol `array' puzzles me) but I cannot do the converse: (grid:copy-to #(1.0d0 2.0d0) 'grid:foreign-array) There is no class named GRID::VECTOR-T. [Condition of type SIMPLE-ERROR] Restarts: 0: [RETRY] Retry SLIME REPL evaluation request. 1: [*ABORT] Return to SLIME's top level. 2: [TERMINATE-THREAD] Terminate this thread (#<THREAD "repl-thread" RUNNING {BD57319}>) Backtrace: 0: (SB-PCL::FIND-CLASS-FROM-CELL GRID::VECTOR-T NIL T) 1: ((SB-PCL::FAST-METHOD MAKE-INSTANCE (SYMBOL)) #<unavailable argument> #<unavailable argument> GRID::VECTOR-T)[:EXTERNAL] 2: ((SB-PCL::FAST-METHOD GRID::MAKE-GRID-DATA ((EQL 'GRID:FOREIGN-ARRAY) T T)) #<unavailable argument> #<unavailable argument> #<unavailable argument> (2) (T))[:EXTERNAL] 3: (GRID:MAP-N-GRIDS)[:EXTERNAL] 4: (SB-INT:SIMPLE-EVAL-IN-LEXENV (GRID:COPY-TO #(1.0d0 2.0d0) 'GRID:FOREIGN-ARRAY) #<NULL-LEXENV>) I get the same results with the function `copy-to' I get the same problem if I am in the package :grid. >From tracing the call tree, I think that the `T' vector-T stems from a call to `array-element-type' which under SBCL gives the following: GRID> (array-element-type #(1.0 2.0)) T Thanks, Mirko

2 2

[GSLL-devel] Incorrect example code in gsl/grid documentation?
by Mirko Vukovic 15 Nov '10

15 Nov '10

Hello, (I downloaded gsd/gsll a few days ago) In the documentation section 3.2.1 map-grid, the following example does not work: (map-grid :source 'index-fill-decadal :source-dims '(3 4)) Instead, it should probably include the `destination-specification': (map-grid :source #'index-fill-decadal :source-dims '(3 4) :destination-specification '((array 3 4) double-float)). The map-grid documentation is a bit unclear. It says that if `destination-specification' is not provided, map-grid defaults to `source'. But source is a function. Looking at map-n-grids code in map.lisp, I don't see a test clause that handles the case of unspecified `destination-specification' Mirko

1 0

Re: [GSLL-devel] [Gsll-devel] Efficient access to externally generated double-float arrays?
by Sebastian Sturm 14 Nov '10

14 Nov '10

Here is something similar to what I suggested, although still very unfinished. It's probably not the proper way to write a macro, but I guess other people on the list will know how to do it correctly. Efficient linearization can be added in the same manner; also, the macro should be modified s.t. it accepts several array specifications at once, i.e. with-foreign-array ((array-1 :double) (array-2 :int) (array-3 :double) ...), etc. Of course, if you can figure out a way to incorporate the optimizations into gref without such a clumsy workaround, I'd be all for it. best regards, Sebastian (defun mapcons (fn x) (if (atom x) x (funcall fn (let ((a (mapcons fn (car x))) (d (mapcons fn (cdr x)))) (if (and (eql a (car x)) (eql d (cdr x))) x (cons a d)))))) (defmacro with-fast-access-to-single-foreign-array ((array element-type) &body body) (alexandria:with-unique-names (array-fptr) `(let ((,array-fptr (grid::foreign-pointer ,array))) ,@(mapcons (lambda (expr) (if (and (consp expr) (eq (first expr) 'grid:gref*) (eq (second expr) array)) (list 'cffi:mem-aref array-fptr element-type (elt expr 2)) expr)) body)))) (defun macro-force-function (dim) "Given an integer dim, this constructs a function that, when supplied with a N-dimensional vector Z and some output vector (-> pointer?), yields the corresponding forces" (declare (fixnum dim)) (let ((temp-values (make-array 2 :element-type 'double-float :initial-element 0.0d0))) (lambda (zvector output) (with-fast-access-to-single-foreign-array (output :double) (with-fast-access-to-single-foreign-array (zvector :double) (do ((i 0 (1+ i))) ((= i dim)) (declare (fixnum i)) (setf (aref temp-values 0) 0.0d0) (do ((m 0 (1+ m))) ((> m i)) (declare (fixnum m)) (do ((n i (1+ n))) ((= n dim)) (declare (fixnum n)) (setf (aref temp-values 1) 0.0d0) (do ((k m (1+ k))) ((> k n)) (declare (fixnum k)) (incf (aref temp-values 1) (grid:gref* zvector k))) (incf (aref temp-values 0) (expt (aref temp-values 1) -2)))) (setf (grid:gref* output i) (- (grid:gref* zvector i) (aref temp-values 0)))))))))

2 6