[SciPy-dev] MCMC, Kalman Filtering, AI for SciPy?

[Travis Oliphant]

Fernando Perez wrote:

> eric jones schrieb:
>
>> Of course, we should get some testing feedback on the main platforms 
>> with the current CVS before doing this.  Besides that, does anyone 
>> have any specific (small) things to address before we do so?  If so, 
>> what is the timeframe?
>
>
> I have 2 things:
>
> 1. I'd like to have the implementation of Frobenius norms fixed.  The 
> current code, in scipy/linalg/basic.py does this:
>
>
>        elif ord in ['fro','f']:
>             X = scipy_base.mat(x)
>             return sqrt(sum(diag(X.H * X)))


good idea.


>
> 2. The matrixmultiply != bug I recently reported is _very_ serious for 
> people working with large matrices.  I also reported it in the numpy 
> list, but nobody replied.  Since scipy can trivially work around it 
> with the one-liner I showed, I think this really should be done.

I think you will see that I released a new version of Numeric that fixed 
this problem.   See 23.5

-Travis


>
> Code for #1 is here (but test it further), and I posted the one-liner 
> for #2 yesterday, so these two (IMHO important) fixes should take no 
> time for a committer.

We could also fix it in SciPy too.