As discussed, here are some suggested speedups for the algorithms in this repository.
All operators are now represented as numPy arrays. This should be backwards-compatible with any code that uses the old version, as long as that code uses only indexing operations on Python lists.
Most of the speedups occur because of numPy's more efficient implementations of array methods. Matrices are internally converted into sparse formats when this helps efficiency of matrix operations.
As discussed, here are some suggested speedups for the algorithms in this repository.
All operators are now represented as numPy arrays. This should be backwards-compatible with any code that uses the old version, as long as that code uses only indexing operations on Python lists.
Most of the speedups occur because of numPy's more efficient implementations of array methods. Matrices are internally converted into sparse formats when this helps efficiency of matrix operations.