Good news for anyone who does a lot of matrix manipulation in their Python codes: PEP 465 was accepted for inclusion into Python 3.5!

The upshot is that an expression like:

  >>> S = (H.dot(beta) - r).T.dot(inv(H.dot(V).dot(H.T))).dot(H.dot(beta) - r)

turns into something much more reasonable and obvious:

 >>> S = (H @ beta - r).T @ inv(H @ V @ H.T) @ (H @ beta - r)

This may be a tasty enough treat to start winning scientists over to Python 3.



    4           7