Hello!
Just cloned git repository. In both SBCL 1.0.23 and ClozureCL 1.2-r11583M on Linux PPC, with GSL 1.9-3:
CL-USER> (gsll:invert-matrix #m((1 0 2) (-1 3 1) (1 1 0))) #<MATRIX-DOUBLE-FLOAT #2A((0.11111111111111116D0 -0.2222222222222222D0 0.6666666666666666D0) (-0.11111111111111109D0 0.22222222222222224D0 0.3333333333333333D0) (0.4444444444444444D0 0.1111111111111111D0 -0.3333333333333333D0))> CL-USER> (gsll:matrix-product * #m(3 2 1)) #<VECTOR-DOUBLE-FLOAT #(0.5555555555555557D0 0.44444444444444453D0 1.222222222222222D0)> CL-USER> (gsll:inverse-matrix-product #m((1 0 2) (-1 3 1) (1 1 0)) #m(3 2 1)) #<VECTOR-DOUBLE-FLOAT #(1D+-0 #| not-a-number |# -1D++0 1D++0)>
Do I misunderstand something?
Sincerely yours, Dmitri
PS and to add insult to injury ;) nonsymmetric-generalized.lisp fails to compile, complaining about invalid function call: ("gsl_eigen_gen" "gsl_eigen_gen_QZ")
On Wed, Feb 25, 2009 at 7:25 AM, Dmitri Hrapof hrapof@common-lisp.ruwrote:
Hello!
Just cloned git repository. In both SBCL 1.0.23 and ClozureCL 1.2-r11583M on Linux PPC, with GSL 1.9-3:
CL-USER> (gsll:invert-matrix #m((1 0 2) (-1 3 1) (1 1 0))) #<MATRIX-DOUBLE-FLOAT #2A((0.11111111111111116D0 -0.2222222222222222D0 0.6666666666666666D0) (-0.11111111111111109D0 0.22222222222222224D0 0.3333333333333333D0) (0.4444444444444444D0 0.1111111111111111D0 -0.3333333333333333D0))> CL-USER> (gsll:matrix-product * #m(3 2 1)) #<VECTOR-DOUBLE-FLOAT #(0.5555555555555557D0 0.44444444444444453D0 1.222222222222222D0)> CL-USER> (gsll:inverse-matrix-product #m((1 0 2) (-1 3 1) (1 1 0)) #m(3 2 1)) #<VECTOR-DOUBLE-FLOAT #(1D+-0 #| not-a-number |# -1D++0 1D++0)>
Do I misunderstand something?
Sincerely yours, Dmitri
My reading of gsl documentation ( http://www.gnu.org/software/gsl/manual/html_node/Level-3-GSL-BLAS-Interface....) is that the matrix is supposed to be triangular (upper or lower).
PS and to add insult to injury ;) nonsymmetric-generalized.lisp fails to compile, complaining about invalid function call: ("gsl_eigen_gen" "gsl_eigen_gen_QZ")
O
Gsll-devel mailing list Gsll-devel@common-lisp.net http://common-lisp.net/cgi-bin/mailman/listinfo/gsll-devel
Mirko Vukovic пишет:
My reading of gsl documentation (http://www.gnu.org/software/gsl/manual/html_node/Level-3-GSL-BLAS-Interface....) is that the matrix is supposed to be triangular (upper or lower).
Hvala!
On Wed, Feb 25, 2009 at 7:25 AM, Dmitri Hrapof hrapof@common-lisp.ru wrote:
Hello!
Just cloned git repository. In both SBCL 1.0.23 and ClozureCL 1.2-r11583M on Linux PPC, with GSL 1.9-3:
CL-USER> (gsll:invert-matrix #m((1 0 2) (-1 3 1) (1 1 0))) #<MATRIX-DOUBLE-FLOAT #2A((0.11111111111111116D0 -0.2222222222222222D0 0.6666666666666666D0) (-0.11111111111111109D0 0.22222222222222224D0 0.3333333333333333D0) (0.4444444444444444D0 0.1111111111111111D0 -0.3333333333333333D0))> CL-USER> (gsll:matrix-product * #m(3 2 1)) #<VECTOR-DOUBLE-FLOAT #(0.5555555555555557D0 0.44444444444444453D0 1.222222222222222D0)> CL-USER> (gsll:inverse-matrix-product #m((1 0 2) (-1 3 1) (1 1 0)) #m(3 2 1)) #<VECTOR-DOUBLE-FLOAT #(1D+-0 #| not-a-number |# -1D++0 1D++0)>
Do I misunderstand something?
No, I get a similar result. I don't know why. The macroexpansion to the foreign function call looks right. You don't by any chance have a C program calling GSL/BLAS that does this calculation, do you? That would help isolate whether the problem is in GSL or GSLL.
Sincerely yours, Dmitri
PS and to add insult to injury ;) nonsymmetric-generalized.lisp fails to compile, complaining about invalid function call: ("gsl_eigen_gen" "gsl_eigen_gen_QZ")
I've fixed and pushed this. There was a bug in how it signaled obsolete-gsl-version (generalized eigensystem problems are not supported in GSL 1.9 so if you try to call those functions it should give you an error).
Thanks for the report.
Liam
Liam Healy пишет:
You don't by any chance have a C program calling GSL/BLAS that does this calculation, do you?
Unfortunately, no.
That would help isolate whether the problem is in GSL or GSLL.
It seems Mirko is right, otherwise uplo woudn't have been supplied a default value.
Thanks for the report.
You are welcome; thank you for the fix!
On Wed, Feb 25, 2009 at 7:25 AM, Dmitri Hrapof hrapof@common-lisp.ruwrote:
Hello!
Just cloned git repository. In both SBCL 1.0.23 and ClozureCL 1.2-r11583M on Linux PPC, with GSL 1.9-3:
stuff deleted ...
Sincerely yours, Dmitri
PS and to add insult to injury ;) nonsymmetric-generalized.lisp fails to compile, complaining about invalid function call: ("gsl_eigen_gen" "gsl_eigen_gen_QZ")
I responded too quickly on my previous email:
On my sbcl 1.0.24 this file compiled fine. Maybe something is amiss with your gsl? But I would have no clue as to how to debug that.
Mirko
On Wed, Feb 25, 2009 at 10:05 AM, Mirko Vukovic mirko.vukovic@gmail.com wrote:
PS and to add insult to injury ;) nonsymmetric-generalized.lisp fails to compile, complaining about invalid function call: ("gsl_eigen_gen" "gsl_eigen_gen_QZ")
I responded too quickly on my previous email:
On my sbcl 1.0.24 this file compiled fine. Maybe something is amiss with your gsl? But I would have no clue as to how to debug that.
No, it's your version of GSL, which is probably >=1.10. The bug was not triggered unless you had an old GSL version (in his case, 1.9). But it was a genuine bug; it records the name of the GSL function so it can print it out when the error is signaled on function call, but didn't take into account the possibility that there might be a list of functions, as in this case. It's been fixed.
Liam
-
Dmitri Hrapof -
Liam Healy -
Mirko Vukovic