<br><br><div class="gmail_quote">On Wed, Feb 25, 2009 at 7:25 AM, Dmitri Hrapof <span dir="ltr"><<a href="mailto:hrapof@common-lisp.ru">hrapof@common-lisp.ru</a>></span> wrote:<br><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
Hello!<br>
<br>
Just cloned git repository. In both SBCL 1.0.23 and ClozureCL<br>
1.2-r11583M on Linux PPC, with GSL 1.9-3:<br>
<br>
CL-USER> (gsll:invert-matrix #m((1 0 2) (-1 3 1) (1 1 0)))<br>
#<MATRIX-DOUBLE-FLOAT #2A((0.11111111111111116D0 -0.2222222222222222D0<br>
0.6666666666666666D0)<br>
(-0.11111111111111109D0 0.22222222222222224D0<br>
0.3333333333333333D0)<br>
(0.4444444444444444D0 0.1111111111111111D0<br>
-0.3333333333333333D0))><br>
CL-USER> (gsll:matrix-product * #m(3 2 1))<br>
#<VECTOR-DOUBLE-FLOAT #(0.5555555555555557D0 0.44444444444444453D0<br>
1.222222222222222D0)><br>
CL-USER> (gsll:inverse-matrix-product #m((1 0 2) (-1 3 1) (1 1 0)) #m(3<br>
2 1))<br>
#<VECTOR-DOUBLE-FLOAT #(1D+-0 #| not-a-number |# -1D++0 1D++0)><br>
<br>
Do I misunderstand something?<br>
<br>
Sincerely yours,<br>
Dmitri<br>
</blockquote><div><br>My reading of gsl documentation (<a href="http://www.gnu.org/software/gsl/manual/html_node/Level-3-GSL-BLAS-Interface.html">http://www.gnu.org/software/gsl/manual/html_node/Level-3-GSL-BLAS-Interface.html</a>) is that the matrix is supposed to be triangular (upper or lower).<br>
</div><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;"><br>
PS and to add insult to injury ;) nonsymmetric-generalized.lisp fails to<br>
compile, complaining about invalid function call: ("gsl_eigen_gen"<br>
"gsl_eigen_gen_QZ")</blockquote><div><br>O <br></div><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;"><br>
<br>
_______________________________________________<br>
Gsll-devel mailing list<br>
<a href="mailto:Gsll-devel@common-lisp.net">Gsll-devel@common-lisp.net</a><br>
<a href="http://common-lisp.net/cgi-bin/mailman/listinfo/gsll-devel" target="_blank">http://common-lisp.net/cgi-bin/mailman/listinfo/gsll-devel</a><br>
</blockquote></div><br>