[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
Re: Different numerical results between sparse and full matrices
From: |
Jordi Gutiérrez Hermoso |
Subject: |
Re: Different numerical results between sparse and full matrices |
Date: |
Fri, 17 Feb 2012 08:02:05 -0500 |
On 17 February 2012 07:41, Michael Goffioul <address@hidden> wrote:
> octave.exe:28> sparse(x)./sparse(0)
> ans = Compressed Column Sparse (rows = 1, cols = 1, nnz = 1 [100%])
[snip]
> In other words, complex division by 0 in sparse matrices gives NaN
> instead of Inf. Does anybody know where I should look at to find where
> this could come from?
They're in liboctave/Sparse-op-defs.h... which seems to be one giant,
endless macro. I'm now itching to rewrite it without macros, if you
don't beat me to it.
- Jordi G. H.