Tag Info

Accepted

Complexity of matrix inversion in numpy

(This is getting too long for comments...) I'll assume you actually need to compute an inverse in your algorithm.1 First, it is important to note that these alternative algorithms are not actually ...
Accepted

Alternatives to numpy.einsum

Pyhton 3.5 will introduce a new operator @, which was proposed by NumPy devs to be the matrix multiplication operator. You may want to read PEP 465, but the syntax ...
• 255
Accepted

• 8,219

Calculate determinant of unitary matrices based on SVD implementation

For part 1: To my knowledge, the answer is no. For part 2: This question had several good answers, all of which were negative (there isn't really a faster way). I don't believe there is meaningful new ...
• 1,357
Accepted

Stochastic SIR using SDEint python package

This model is implemented using Julia's DifferentialEquations.jl in this tutorial. Here's a version of that code: ...
Accepted

Performing a random walk on a lattice that traps the particles

This can be solved much more efficiently if you recast the problem as finding the mean first passage time (MFPT) of a Markov chain for a given starting state. You can easily represent the random walk ...
• 3,954