November 2012 - antik-devel - mailman3.common-lisp.net

[Antik-devel] FLAME and possible next library ports using grids
by Liam Healy 08 Dec '12

08 Dec '12

On this and other mailing lists, there have been recent discussions of "numerical lisp" which in one aspect or another has been the vision of many of us, some for many years. My approach is to write libraries I think will form a piece of the overall solution, with the hope that others will join in and write compatible libraries that form another piece. So, I initially wrote GSLL and then split the grid part out into Antik (which has other features I adapted from my personal libraries). The grid functions can use enhancement, but they are usable now for developing library ports. I have thought about and had discussions with others about what would be a good candidate library for the next port. These thoughts have revolved mainly around an FFT library (like fftw) or a linear algebra library (like *lapack) because those are widely used and, since GSL already provides some of that capability, offers a good opportunity to compare features and make a unified interface. Recently, I came across the FLAME linear algebra library (http://z.cs.utexas.edu/wiki/flame.wiki/FrontPage) It is designed with an intriguing modularity concept. I think it would be a good library to port to grid/Antik. Unfortunately, I don't have the time to work on this, but if anyone (particularly someone who needs linear algebra capability like this) wants to take it up, I encourage that and will help as much as I can. A big challenge I foresee in any foreign library interface is the allocation of arrays. If it is a "closed box," it makes it hard to use e.g. static-vectors and therefore the arrays have to be copied back and forth between the foreign side and the CL side. GSL has the advantage that there are functions that allow building GSL arrays out of supplied C arrays, so regular CL arrays can be allocated with a C pointer passed to these functions to wrap them with GSL metadata. I'm not sure many other libraries are that flexible; many (e.g. C++ libraries) consider it a virtue to seal off the raw storage (C arrays) from the user, and force you to go through a foreign interface to do the allocation. Since there are almost no CLs that permit construction of CL arrays using foreign data, copying is the only option left. So, if anyone has thoughts on FLAME or any other library, I'll be happy to contribute my ideas if it will help. Liam

2 2

[Antik-devel] test, disregard: does the list daemon accept non-plain-text emails
by Mirko Vukovic 25 Nov '12

25 Nov '12

Test email in gmail's rich-format. Please disregard. I apologize for this test. I could not find information on web whether only plain-text emails were accepted or not. This email is sent from gmail using rich formatting. Mirko

1 1

[Antik-devel] nested vector structure
by Mirko Vukovic 17 Nov '12

17 Nov '12

I have created a bit of code for gsll that uses a vector of vectors. I use this for tabular data (https://github.com/mirkov/data-table) I use raw CL for that, even though the outer vector can store foreign arrays. That way, I can do interpolation and fitting. The code is still a bit raw, and the object model may need cleaning up. A nice thing is that I can query and reshuffle the table, much like the database features in chapter 27 of Seibel's Practical Common Lisp. I would like to make this more grid compatible, but I am just too tired to figure it out. So, is it possible to - create a vector of vectors in grid - can grid use aref to access individual elements, or am I stuck with a nested aref Thanks, Mirko

2 3

[Antik-devel] Running antik on clisp+cygwin
by Mirko Vukovic 04 Nov '12

04 Nov '12

I am summarizing the changes I had to make to various pieces of code to get antik to compile and run on clisp+cygwin. I had to make changes in two libraries: fsbv and antik. They are all rather minor. *FSBV* You will need to install cygwin's libffi4 package Two modifications in libffi-unix.lisp: In the library loading sections, point to cygwin's ``ffi.h'' like this. #+cygwin (include "ffi.h") We comment out a unix64 section in ``cenum abi'' definition like so: (cenum abi ((:default-abi "FFI_DEFAULT_ABI")) ((:sysv "FFI_SYSV")) #-cygwin((:unix64 "FFI_UNIX64"))) (Cygwin is not a 32-bit application) To make sure things are workinig, you will need to modify, compile and load `examples.lisp'. The library section of "examples.lisp" is modified to load the correct cygwin libraries like so: (cffi:load-foreign-library #+(and unix (not cygwin)) "libgslcblas.so" #+cygwin "cyggslcblas-0.dll") ;; note that we load cygblas explicitly. #+cygwin (cffi:load-foreign-library "/usr/lib/lapack/cygblas-0.dll") (cffi:load-foreign-library #+(and unix (not cygwin)) "libgsl.so" #+cygwin (probe-file "/usr/bin/cyggsl-0.dll")) Once you have modified fsbv like that, make sure it compiles. The go to examples.lisp, compile and load. Then execute the few examples listed at the beginning of the file *ANTIK* In init/package.lisp (defpackage antik) I added the following symbols to be shaddowed: #:second #:fifth #:tenth #:rem Now, since CL's `second' is shaddowed, I had to prepend it with `CL:' in a few places. This was necessary in three files in physical-quantities: angle-component, physical-quantities and units.lisp. Here are the instances: >grep "cl:second" *.lisp angle-component.lisp: (cl:/ (first tangent) (cl:second tangent))) angle-component.lisp: (eql (cl:second (tangent angle)) label)) angle-component.lisp: (setf (cl:second (tangent angle)) run) physical-quantities.lisp: (parse-integer (cl:second prod)) physical-quantities.lisp: (all-same (cl:second mat) :test #'equal) ; all same units physical-quantities.lisp: (cl:second mat)) units.lisp: (or (cl:second (print-name unit)) units.lisp: (pux (cl:second expr) expon) units.lisp: (pux (cl:second expr) (cl:- expon))) ; one argument, invert units.lisp: (pux (cl:second expr) (cl:* expon (third expr)))) units.lisp: (if style (cl:second (print-name unit-sexp)) units.lisp: (mkstr (if (eql (cl:second unit-sexp) 1) units.lisp: (make-pq-string (cl:second unit-sexp) style)) units.lisp: (mkstr (make-pq-string (cl:second unit-sexp) style) units.lisp: (make-pq-string (cl:second units) (nf-option style)) With this, antik should compile cleanly on clisp+cygwin. HTH, Mirko

1 0