RE: Octave 3.0

[Andy Adler]

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