[SciPy-dev] Shape mismatch sparse/dense matrix vector multiplication

[Nils Wagner]

On Fri, 20 Oct 2006 11:54:09 -0500
  Robert Kern <robert.kern at gmail.com> wrote:
> Nils Wagner wrote:
>> Hi all,
>> 
>> I am confused by the the result of the attached script.
>> Is that a bug ?
> 
> What result did you get? What result did you expect?
> 
> -- 
> Robert Kern
> 
> "I have come to believe that the whole world is an 
>enigma, a harmless enigma
>  that is made terrible by our own mad attempt to 
>interpret it as though it had
>  an underlying truth."
>   -- Umberto Eco
> _________________________________,______________
> Scipy-dev mailing list
> Scipy-dev at scipy.org
> http://projects.scipy.org/mailman/listinfo/scipy-dev

What it makes worse is the fact that
the norm of $\| yd-ys\|$ is not zero but something
like 1.e-7. I have expected a value near machine precision
(1.e-16).

Any idea how to fix that ?

Nils