August 2008 - gsll-devel - mailman3.common-lisp.net

[Gsll-devel] Thanks Liam!
by Mirko Vukovic 02 Sep '08

02 Sep '08

Used the mean&variance stuff yesterday. Thanks (again)! Mirko

2 2

[Gsll-devel] Fwd: [Help-gsl] Functions without documentation
by Liam Healy 12 Aug '08

12 Aug '08

I asked on the help-gsl mailing list, and the status of the matrix exponential is that, while present in GSL, it is "not supported" due to lack of testing. If you would like to take up the issue of how to get it supported, I suggest you send email to the help-gsl list. Liam ---------- Forwarded message ---------- From: Brian Gough <bjg(a)gnu.org> Date: Mon, Aug 11, 2008 at 12:21 PM Subject: Re: [Help-gsl] Functions without documentation To: Liam Healy <lnp(a)healy.washington.dc.us> Cc: Help GSL <Help-gsl(a)gnu.org> At Sun, 10 Aug 2008 16:59:20 -0400, Liam Healy wrote: > > There are three functions > gsl_linalg_matmult > gsl_linalg_matmult_mod > gsl_linalg_exponential_ss > in /usr/include/gsl/gsl_linalg.h and the library, but apparently don't > have any documentation. Are they in working order? Is there any > documentation for them? Should it be included with the rest of the > GSL documentation? Hello, Functions not documented in the manual are considered obsolete or unsupported. In particular gsl_linalg_exponential_ss is lacking any tests, so it is not supported. -- Brian Gough GNU Scientific Library - http://www.gnu.org/software/gsl/

2 1

[Gsll-devel] gsl_linalg_exponential_ss
by Tamas K Papp 11 Aug '08

11 Aug '08

Hi Liam, Is there a binding for gsl_linalg_exponential_ss (that calculates the exponential of a matrix?) I could not find one. Thanks, Tamas

2 5

[Gsll-devel] keeping m-variables in closures
by Mirko Vukovic 10 Aug '08

10 Aug '08

Hello folks (the following is not a showstopper, but I am just curious). I am trying to define a closure (still doing LU decompositions and matrix solutions), and would like to keep some of the variables around for further processing. Here is one example (defun parallel-plate-equilibrium (eps0 T0 eps1- eps1+ eps2 T2 &key (q1 0d0)) (letm ((dim 5) (mat (make-array (list dim dim) :element-type 'double-float :initial-contents (list (list 1d0 (1- eps0) 0d0 0d0 0d0) (list (1- eps1-) 1d0 0d0 0d0 (0- eps1-)) (list 0d0 0d0 1d0 (1- eps1+) (0- eps1+)) (list 0d0 0d0 (1- eps2) 1d0 0d0) (list (0- eps1-) 0d0 0d0 (0- eps1+) (+ eps1- eps1+))))) (mmat (matrix-double-float mat)) ;; (mmat0 mmat) (rhs (make-array dim :element-type 'double-float :initial-contents (list (st4 T0 eps0) 0d0 0d0 (st4 T2 eps2) q1))) (mrhs (vector-double-float rhs)) (per (permutation dim)) (res (vector-double-float dim)) (0-vec (vector-double-float (make-array dim :initial-element 0d0))) (x (vector-double-float dim))) (lu-decomp mmat per) (lu-solve mmat per mrhs x) (data x) (lambda (command) (case command ;; these two will fail because mmat, mrhs, mmat0 are nil ;; (:solve (lu-solve mmat per mrhs x)) ;; (:check (gemv :notrans 1d0 mmat0 x 0d0 0-vec)) (:mat mat) (:rhs rhs) (:temp (expt (/ (aref (data x) 4) +sigma+) 0.25)) (:fluxes (data x)))))) However, it seems that gsll's variables are discarded once the letm form is done (as per documentation, and also with macroexpansion). But the code for letm mentions some "freers", which gave me hope that I can specify the freer as NIL. But I could not figure that one out quite yet. So, is it possible to save some of gsll's variables for later use? Thanks, Mirko

2 6

Re: [Gsll-devel] cannot compile lu.lisp's invert-matrix example
by Mirko Vukovic 06 Aug '08

06 Aug '08

Got it. Indeed, it works fine. In my lisp code this message implies a conceptual or programming error. So, I always try eliminate it. Thus the knee-jerk reaction. I am also wondering: here do I ship the beer (or whatever) :-) This package is sooo useful Thanks, Mirko On Tue, Aug 5, 2008 at 9:21 PM, Liam Healy <liam(a)healy.washington.dc.us> wrote: > If that's all you got, nothing is wrong. Notes are not errors and not > even warnings, > they are just that, notes. Did you try to run it? Please see the SBCL > docs for > an explanation of notes. > > On Tue, Aug 5, 2008 at 1:29 PM, Mirko Vukovic <mirko.vukovic(a)gmail.com> > wrote: > > Hello Liam (& others) > > > > I am using SBCL 1.014 and GSL with a timestamp from march of this year. > > > > I am trying to use the lu.lisp library and cannot compile the > invert-matrix > > example. The error message is "deleting unreachable code" at the (letm. > > If I comment out the calls to lu-... and (data-inv) and do instead (data > > mmat), then it will compile. > > > > I am including the transcript at the end of this message. > > > > Any clue as to what could be wrong? > > > > Thanks, > > > > Mirko > > > > > > CL-USER> (in-package :gsl) > > #<PACKAGE "GSLL"> > > GSL> (defun invert-matrix (mat) > > "Invert the matrix." > > (letm ((mmat (matrix-double-float mat)) > > (dim (array-dimension mat 0)) > > (per (permutation dim)) > > (inv (matrix-double-float dim dim))) > > (lu-decomp mmat per) > > (lu-invert mmat per inv) > > (data inv))) > > ; in: LAMBDA NIL > > ; (GSLL:LETM > > ; ((GSLL::MMAT (GSLL:MATRIX-DOUBLE-FLOAT GSLL::MAT)) > > ; (GSLL::DIM (ARRAY-DIMENSION GSLL::MAT 0)) > > ; (GSLL::PER (GSLL:PERMUTATION GSLL::DIM)) > > ; (GSLL::INV (GSLL:MATRIX-DOUBLE-FLOAT GSLL::DIM GSLL::DIM))) > > ; (GSLL:LU-DECOMP GSLL::MMAT GSLL::PER) > > ; (GSLL:LU-INVERT GSLL::MMAT GSLL::PER GSLL::INV) (GSLL:DATA > GSLL::INV)) > > ; --> LET* UNWIND-PROTECT FLET BLOCK MULTIPLE-VALUE-BIND > > ; --> MULTIPLE-VALUE-CALL BLOCK SB-C::%WITHIN-CLEANUP RETURN-FROM PROGN > > ; --> UNLESS COND IF PROGN SETF LET* > > ; ==> > > ; GSLL::PER > > ; > > ; note: deleting unreachable code > > > > ; ==> > > ; GSLL::INV > > ; > > ; note: deleting unreachable code > > ; > > ; compilation unit finished > > ; printed 2 notes > > STYLE-WARNING: redefining INVERT-MATRIX in DEFUN > > INVERT-MATRIX > > > > _______________________________________________ > > Gsll-devel mailing list > > Gsll-devel(a)common-lisp.net > > http://common-lisp.net/cgi-bin/mailman/listinfo/gsll-devel > > > > >

1 0

[Gsll-devel] cannot compile lu.lisp's invert-matrix example
by Mirko Vukovic 05 Aug '08

05 Aug '08

Hello Liam (& others) I am using SBCL 1.014 and GSL with a timestamp from march of this year. I am trying to use the lu.lisp library and cannot compile the invert-matrix example. The error message is "deleting unreachable code" at the (letm. If I comment out the calls to lu-... and (data-inv) and do instead (data mmat), then it will compile. I am including the transcript at the end of this message. Any clue as to what could be wrong? Thanks, Mirko CL-USER> (in-package :gsl) #<PACKAGE "GSLL"> GSL> (defun invert-matrix (mat) "Invert the matrix." (letm ((mmat (matrix-double-float mat)) (dim (array-dimension mat 0)) (per (permutation dim)) (inv (matrix-double-float dim dim))) (lu-decomp mmat per) (lu-invert mmat per inv) (data inv))) ; in: LAMBDA NIL ; (GSLL:LETM ; ((GSLL::MMAT (GSLL:MATRIX-DOUBLE-FLOAT GSLL::MAT)) ; (GSLL::DIM (ARRAY-DIMENSION GSLL::MAT 0)) ; (GSLL::PER (GSLL:PERMUTATION GSLL::DIM)) ; (GSLL::INV (GSLL:MATRIX-DOUBLE-FLOAT GSLL::DIM GSLL::DIM))) ; (GSLL:LU-DECOMP GSLL::MMAT GSLL::PER) ; (GSLL:LU-INVERT GSLL::MMAT GSLL::PER GSLL::INV) (GSLL:DATA GSLL::INV)) ; --> LET* UNWIND-PROTECT FLET BLOCK MULTIPLE-VALUE-BIND ; --> MULTIPLE-VALUE-CALL BLOCK SB-C::%WITHIN-CLEANUP RETURN-FROM PROGN ; --> UNLESS COND IF PROGN SETF LET* ; ==> ; GSLL::PER ; ; note: deleting unreachable code ; ==> ; GSLL::INV ; ; note: deleting unreachable code ; ; compilation unit finished ; printed 2 notes STYLE-WARNING: redefining INVERT-MATRIX in DEFUN INVERT-MATRIX

1 0