[SciPy-user] Interaction of sparse and dense matrices in scipy

Nils Wagner nwagner at mecha.uni-stuttgart.de
Tue Feb 21 06:04:20 CST 2006


Ed Schofield wrote:
>Nils Wagner wrote:
>
>  
>>Hi Ed,
>>
>>Thank you for your valuable comments.
>>The multiplication a*b, where a is dense a b is sparse, is certainly of
>>importance.
>>For example Block Krylov methods make heavily use of such products.
>>Anyway, it would be nice if scipy can handle such products :-)
>> 
>>
>>    
>It does work now where 'a' is an *array*, but not when 'a' is a (dense)
>matrix.  For now, just use a dense array :)
>
>The reason is that the matrix object overrides the * operator and calls
>dot(a, b).  I'll think about this problem when I have time and make some
>proposals on scipy-dev and/or numpy-discussion for how to solve it...
>
>-- Ed
>
>_______________________________________________
>SciPy-user mailing list
>SciPy-user at scipy.net
>http://www.scipy.net/mailman/listinfo/scipy-user
>  
Hi Ed,

May I send you my (sparse/dense) codes for solving large Sylvester
equations.
The results of the dense version are o.k. while the sparse version
doesn't work as expected.
I hope you can drop some hints ;-) .

Cheers,

Nils
 



More information about the SciPy-user mailing list