181 karma · joined May 15, 2013
They started behind the curve, but they're moving quickly, relatively at least
There's a couple of cases with float64 vs. complex128 matrices, but I have been annoyed with those silent changes in Matlab where the answer is wrong but the code continues anyway.
The performance story is complex. Typically we're the same speed on small matrices (and using Go is faster if you include the cgo overhead). We currently have significant speed penalties on large matrices (300x300 or so), but Kunde21 is working on assembly kernels for the BLAS functions to close that gap
And yea, mat mul is not N^3 theoretically, but most implementations are. I've heard that some (mkl maybe) are 2.8, but haven't had someone point code to me. My personal attempts at implementing Strassen were slower than a tuned N^3 implementation, at least for matrices that fit into memory.
It's actually the second highest eigenvalue. The highest eigenvalue is always 1 for stochastic matrices.
> Go doesn't have a nice REPL like IPython, and it never will.
In particular, there is a nice illustration of "inertia doesn't matter" at 15 minutes.
Actually, No Free Lunch guarantees that such cases exist.