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 will have to do.
On Tue, Apr 9, 2013 at 1:53 AM, Michal D. <[email protected]>wrote:
> 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 which two of the three you're referring: TCL, Haskell, J.
> ;-)
>
> Cheers,
>
> Mike
>
>
> On Mon, Apr 8, 2013 at 10:15 AM, Johann Hibschman <[email protected]
> >wrote:
>
> > On Mon, Apr 8, 2013 at 12:45 AM, Michal D. <[email protected]>
> > 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.
> > ----------------------------------------------------------------------
> > For information about J forums see http://www.jsoftware.com/forums.htm
> >
> ----------------------------------------------------------------------
> For information about J forums see http://www.jsoftware.com/forums.htm
>
--
Devon McCormick, CFA
^me^ at acm.
org is my
preferred e-mail
----------------------------------------------------------------------
For information about J forums see http://www.jsoftware.com/forums.htm
