Re: [Jprogramming] Eigen values in J - example

2013-04-09 Thread Michal D.
Thanks a lot, Igor, for the example and all the work you put into the library. Both methods (1) and (2) chug along for a few minutes but end up spitting out a result that looks reasonable. Cheers, Mike On Tue, Apr 9, 2013 at 12:51 AM, Igor Zhuravlov wrote: > On 8 апреля 2013 22:53:13 Michal

Re: [Jprogramming] Eigen values in J

2013-04-09 Thread Devon McCormick
Hi - did you see something like this from Igor Z? load 'math/mt' geevlnn=: {.@ggevlnn_mt_@(,:(idmat_mt_@c_mt_)) geevlnn i.2 2 _0.561553 3.56155 geevlnn p: i.4 4 104.336 _4.54479 1.02295 _1.81419 I don't know what's a canonical matrix to use for eigenvalue demonstration, so these wil

Re: [Jprogramming] Eigen values in J - mt addon fixed

2013-04-09 Thread Devon McCormick
Thanks - I had a good idea what it was from the name and usage. Thanks for putting this package together. It's good to have a portable J-based solution. On Mon, Apr 8, 2013 at 11:27 PM, Igor Zhuravlov wrote: > On 8 апреля 2013 11:51:05 Devon McCormick wrote: > > I attempted to use one of the

Re: [Jprogramming] Eigen values in J - example

2013-04-09 Thread Igor Zhuravlov
On 8 апреля 2013 22:53:13 Michal D. wrote: > I haven't been able to sort through the mt code to decipher how to call the > routines with a single input matrix. It looks like quite a substantial > amount of important work. load'math/mt' NB. generate random symmetric 4x4-matrix to play with

Re: [Jprogramming] Eigen values in J

2013-04-08 Thread Michal D.
The good news is my matrix is symmetric, the bad news is that the jacobi routine never seems to terminate. I haven't been able to sort through the mt code to decipher how to call the routines with a single input matrix. It looks like quite a substantial amount of important work. I'm not sure to

Re: [Jprogramming] Eigen values in J - mt addon fixed

2013-04-08 Thread Igor Zhuravlov
On 8 апреля 2013 11:51:05 Devon McCormick wrote: > I attempted to use one of the math/mt routines (geev) but get a value error > on "trsmu1x" which does not appear to be defined anywhere. Thank you for reporting this bug. Now it is fixed. Please, update mt addon to version 0.9.1. -- WBR Igor -

Re: [Jprogramming] Eigen values in J

2013-04-08 Thread Johann Hibschman
On Mon, Apr 8, 2013 at 12:45 AM, Michal D. wrote: > Ps. There is nothing I like better than generating random input data using > J. What took pages of Haskell and TCL was boiled down to 7 clear lines of > J. WOW =) Now there are two languages rarely mentioned in the same breath. ---

Re: [Jprogramming] Eigen values in J

2013-04-08 Thread km
Did you begin by loading ~addons/math/mt/mt.ijs ? --Kip Murray Sent from my iPad On Apr 8, 2013, at 10:51 AM, Devon McCormick wrote: > I attempted to use one of the math/mt routines (geev) but get a value error > on "trsmu1x" which does not appear to be defined anywhere. > > There is a Jaco

Re: [Jprogramming] Eigen values in J

2013-04-08 Thread Devon McCormick
I attempted to use one of the math/mt routines (geev) but get a value error on "trsmu1x" which does not appear to be defined anywhere. There is a Jacobi solver in ~addons/math/misc/jacobi that appears to work for a symmetric argument. Also, there's the code from Donald McIntyre's article on Jaco

Re: [Jprogramming] Eigen values in J

2013-04-08 Thread Kip Murray
Igor Zhuravlof provides j routines that model LAPACK routines for eigenvalues and eigenvectors. See "matrix toolbox" ~addons/math/mt with contents summarized in ~addons/math/mt/mt.ijs I have not tried them but would expect them to run in j701JHS. --Kip Murray On 4/8/2013 12:45 AM, Michal D

Re: [Jprogramming] Eigen values in J

2013-04-07 Thread bill lam
the lapack package has been updated to support 64-bit linux. Вс, 07 апр 2013, Michal D. писал(а): > Hi All, > > I have been trying to get J to compute the eigen values & vectors of a > matrix. I've come across > http://www.jsoftware.com/jwiki/Essays/Eigenvaluesbut none of the > algorithms there

[Jprogramming] Eigen values in J

2013-04-07 Thread Michal D.
Hi All, I have been trying to get J to compute the eigen values & vectors of a matrix. I've come across http://www.jsoftware.com/jwiki/Essays/Eigenvaluesbut none of the algorithms there seem stable on my input matrix and the lapack package doesn't seem to work on 64-bit J. I laboriously exported