Closed gonced8 closed 3 months ago
It makes sense to keep the optimized version. The original version was just a literal translation from the Algorithm 8 pseudocode (for m = 1 to M′ do:
), but because we are using matrix multiplications in other loops of this algorihtm, so it makes sense to change this too.
The for loop along the multiple samples is unnecessary and can be simplified with just a matrix multiplication. There are possible other places in the code where the same simplification can be implemented (probably search for loops computing np.outer).