February 2013 - antik-devel - mailman3.common-lisp.net

[Antik-devel] Bug in antik/.../sample/low-discrepancy-sequence.lisp ?
by Mirko Vukovic 17 Apr '13

17 Apr '13

Hello, this is on clisp 2.48 on cygwin. The function low-discrepancy-sample could not compile because iterate could not digest loop ***** Iterate, in clause (COLLECT (APPLY MAKE-FUNCTION (APPEND FIXED-PARAMETERS (MAPCAN (LAMBDA (X SEQ) (LIST (FIRST X) (VAL-FROM-RANGE (GREF V SEQ) (SECOND X) (THIRD X)))) SWEEP-PARAMETERS (LOOP FOR I FROM 0 BELOW (LENGTH SWEEP-PARAMETERS) COLLECT I))))): Iterate does not know how to handle the special form (MACROLET ((LOOP-FINISH NIL (SYSTEM::LOOP-FINISH-ERROR))) (BLOCK NIL (LET ((I 0) (#:LIMIT-273656 (LENGTH SWEEP-PARAMETERS))) (PROGN (LET ((#:ACCULIST-VAR-273657 NIL)) (MACROLET ((LOOP-FINISH NIL '(GO SYSTEM::END-LOOP))) (TAGBODY SYSTEM::BEGIN-LOOP (WHEN (COMMON-LISP:>= I #:LIMIT-273656) (LOOP-FINISH)) (PROGN (SETQ #:ACCULIST-VAR-273657 (CONS I #:ACCULIST-VAR-273657))) (PSETQ I (COMMON-LISP:+ I 1)) (GO SYSTEM::BEGIN-LOOP) SYSTEM::END-LOOP (MACROLET ((LOOP-FINISH NIL (SYSTEM::LOOP-FINISH-WARN) '(GO SYSTEM::END-LOOP))) (RETURN-FROM NIL (SYSTEM::LIST-NREVERSE #:ACCULIST-VAR-273657)))))))))) It will not be walked, which means that Iterate clauses inside it will not be seen. ***** If I replaced the loop form with the much more verbose (let (res) (dotimes (i (length sweep-parameters) (nreverse res)) (push i res))) the code compiled cleanly. The same error can be seen with a much simpler: (iter:iter (iter:for v iter:in '(1 2 3 4)) (iter:collect (loop for i in v))) Another alternative is to use iter instead of loop (iter:iter (iter:for i from 0 below (length sweep-parameters)) (iter:collect i)) It compiles cleanly, but I did not test the function accuracy. Mirko

2 1

Re: [Antik-devel] Matlisp
by Liam Healy 01 Mar '13

01 Mar '13

Hi Akshay, Thank you for your work and your email. I think a lot of people have had very similar thoughts to yours, and have started down the road meeting the needs of users they've identified. There are a lot of good and partially overlapping ideas within this loosely organized group of people. There was a similar discussion started in the context of common lisp stats (which I'm not a member of) a few months ago. Where I was participating, I cc'ed the antik-devel mailing list; see for example http://lists.common-lisp.net/pipermail/antik-devel/2012q4/000019.html and the archives in general http://lists.common-lisp.net/pipermail/antik-devel/2012q4/thread.html with the subject "numerical lisp". My feeling is that we need to break the desired functions into relatively narrow, compatible libraries. I started Antik as a split out of GSLL to handle "grids" (generalized arrays) to start down that road. It needs to be split up further, because it now includes other things such as calculations with physical dimensions. I've started that further split, but haven't done any work on it recently. I am interested in the work you've done; it's been a while since I looked at Matlisp. At the time, it did not meet my needs, but it seems you have done significant work on it and has a lot of potential to be part of the compatible family of libraries as I've described. I think a good first step is to identify the different areas (for example, generalized arrays as I've indicated, or statistics functions) and canvas the current state of software for features and capabilities. I've added Tamas Papp and Mirko Vukovic onto this email because they've been very active in thinking about these things. Also I've added the Antik and GSLL mailing lists for good measure. Liam On Wed, Feb 27, 2013 at 1:04 AM, Akshay Srinivasan <akshaysrinivasan(a)gmail.com> wrote: > Hello All, > I'm sort of the active developer of Matlisp these days (mostly whilst procrastinating > from doing more "academic" work). I know all of you have created your independent system for > scientific computing on Common Lisp, and it is really a pity that all of them (including the > old Matlisp) are not particularly usable when compared with something like SciPy or Matlab even; > more so since Lisp is so much nicer a language. > I think the version of Matlisp I'm developing has a lot of promise, and is taking a very > good shape (atleast I'd like to think so). I'm personally trying to use it for my work in Control > theory and such; but it's still not something which is particularly usable yet. The improvements over > the old version of Matlisp are: > - Use of a hacky static object system, so that functions can be generated for any user-defined "field". > There is some code demonstrating the use of Maxima to work with "symbolic-tensors", there isn't > code for integer matrices yet, but getting this to work is fairly trivial. > - The tensor (multi-array) structures can now be sliced inplace very easily like SciPy. > - Automatic fortran callback generation (see dlsode.lisp). There are plans to extend this to C, so > that we can make use of GSL. > - Lisp implementations of BLAS functions, and support for seamlessly switching between the Lisp and > Fortran versions of these routines. The lisp version of GEMM (and others) is nearly as fast as the > F77 reference BLAS one, on SBCL. > - The FFI has moved to CFFI, if you haven't noticed that yet; builds on CCL and CMUCL. > The basic structure of BLAS is out there, but the LAPACK wrappers are missing. However I do think that > the general structure of the code is in place. > > I just wanted to throw out the idea into the open, about the feasibility of merging all your > projects into Matlisp (and calling it by some other name, if that is an issue) ? Would any of you > be willing to take the time to go through my code and merge yours in ? I'm writing to you now > before I take time and polish things, because then there'd be little room to make architectural > changes later on. The code is here if you'd like to give it a try: > git://matlisp.git.sourceforge.net/gitroot/matlisp/matlisp tensor > or if you just want to browse the code: > https://github.com/enupten/matlisp > You may want to look at the (ugly) code in tests/ folder to get an idea about the code behaviour. > You'd also want to work in #:matlisp package (despite what the herald says). > > Getting the maxima bits to work is kind of tricky, but you don't really need it for numerical routines > to work. > > Akshay

2 1

[Antik-devel] Physics-constants library?
by Mirko Vukovic 26 Feb '13

26 Feb '13

I did not find a library with physics constants on cliki or quicklisp. There is one available from within GSLL, but loading GSLL just for physics constants is an overkill. Are there any other ones? I have one of my own (somewhat incomplete) that works for me, and I can clean it up and push it to quicklisp if nothing meaningful is out there. Mirko

2 2