[unitaryhack] Improve the calculation of multidimensional hermite polynomials about thewalrus HOT 10 OPEN

A second improvement would be to find a way to parallelize the loop in line

This might require a nontrivial amount of algorithmic thinking since this loop implements a recursion relation.

from thewalrus.

Finally, it would also be interesting to think about porting this part of the code to numba.

from thewalrus.

Hi @nquesada, I didn't get the second improvement you proposed. Can you elaborate on it? Thanks!

from thewalrus.

I just meant, to write an implementation of Hermite multidimensional in numba.

from thewalrus.

A second improvement would be to find a way to parallelize the loop in line

thewalrus/include/hermite_multidimensional.hpp

Line 114 in bc53689

for (ullint jj = 0; jj < Hdim-1; jj++) {

.
This might require a nontrivial amount of algorithmic thinking since this loop implements a recursion relation.

I was referring to this comment, can you elaborate on it? Thanks!

from thewalrus.

Hi @TripleR47, the idea here is to use numba to replace that line (a serial for loop) with a parallelized version

from thewalrus.

To add to @co9olguy : the challenge is that the for loop in line 114 is very long (dimensions Hdim) and implements a recursion relation, where the next element is calculated based on the values of previous ones. My comment is just to highlight that breaking a parallelzing a recursive loop might be non-trivial!

from thewalrus.

Would love to take a stab at this. Okay if I work on this?

from thewalrus.

Hi @e-eight, yes this issue is still open. Feel free to tackle it if interested :)

from thewalrus.

Alright. Will give it a shot.

from thewalrus.