It looks like the determinant conjunction . gives functions that are
explicitly rank 2:
-/ .* b. 0
2 _ _
This means we can't use a matrix of rank-1 polynomials as input without
redefining the conjunction ourselves.
Marshall
On Sun, Jun 10, 2012 at 10:05 PM, Raul Miller <rauldmil...@gmail.com> wrote:
> I am looking for a concise expression to find eigenvalues in J (for
> small matrices)
>
> lapack is not universally available, so I am ignoring it.
>
> Here's the best I have been able to find (loosely based on an "At Play
> With J" article)
>
> ev=: 1 {:: [: p.@> [: -&.>/ .(+//.@(*/)&.>) ] <@,"0 -@=@i.@#
>
> Is there some way to eliminate the boxing and unboxing here?
>
> Thanks,
>
> --
> Raul
> ----------------------------------------------------------------------
> For information about J forums see http://www.jsoftware.com/forums.htm
>
----------------------------------------------------------------------
For information about J forums see http://www.jsoftware.com/forums.htm
