[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
RE: Octave 3.0
From: |
Andy Adler |
Subject: |
RE: Octave 3.0 |
Date: |
Fri, 28 Feb 2003 11:23:38 -0500 (EST) |
On Fri, 28 Feb 2003, Schloegl Alois wrote:
> I was referring to John's priority list. IMO, it is a good thing, if sparse
> arrays become part of core octave.
>
> Unfortunately, I've only access to some outdated versions of Octave (2.1.40),
> without a recent version of octave-forge. Hence, I can not say whether I miss
> a
> specific feature in your sparse package. Only, the following example should
> return the result below:
>
I get exactly your results:
octave:1> x = [-inf, -1, 0, 1, inf,nan];
octave:2> sx= sparse(x);
octave:3> full(x'*x)
ans =
Inf Inf NaN -Inf -Inf NaN
Inf 1 0 -1 -Inf NaN
NaN 0 0 0 NaN NaN
-Inf -1 0 1 Inf NaN
-Inf -Inf NaN Inf Inf NaN
NaN NaN NaN NaN NaN NaN
octave:4> full(sx'*sx)
ans =
Inf Inf 0 -Inf -Inf NaN
Inf 1 0 -1 -Inf NaN
0 0 0 0 0 0
-Inf -1 0 1 Inf NaN
-Inf -Inf 0 Inf Inf NaN
NaN NaN 0 NaN NaN NaN
andy