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
On Sat, Nov 17, 2012 at 9:20 AM, Liam Healy lhealy@common-lisp.net wrote:
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.
Stuff deleted ...
So, if anyone has thoughts on FLAME or any other library, I'll be happy to contribute my ideas if it will help.
I would like to see links to additional special functions, such as http://www.netlib.org/amos/readme for Bessel functions of complex order. This one is written in FORTRAN-77. Several years ago, I did make links to it, but it was an ugly and painful effort using another library that does not function anymore. The next time the need arises, I will try again, and post the results to encourage wider usage.
Mirko
On Sun, Nov 25, 2012 at 11:05 AM, Mirko Vukovic mirko.vukovic@gmail.comwrote:
On Sat, Nov 17, 2012 at 9:20 AM, Liam Healy lhealy@common-lisp.netwrote:
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.
Stuff deleted ...
So, if anyone has thoughts on FLAME or any other library, I'll be happy to contribute my ideas if it will help.
I would like to see links to additional special functions, such as http://www.netlib.org/amos/readme for Bessel functions of complex order. This one is written in FORTRAN-77. Several years ago, I did make links to it, but it was an ugly and painful effort using another library that does not function anymore. The next time the need arises, I will try again, and post the results to encourage wider usage.
Mirko
I would also like to see interpolation on irregular grids using kridging
or interpolation of Delauny triangulated grids.
One of these days, when the needs really strikes, I will try to link CL to the triangle library (https://www.cs.cmu.edu/~quake/triangle.html -- I tested the C code on linux, and there is documentation for linking it as a library.)
Mirko
-
Liam Healy -
Mirko Vukovic