it gives wrong numbers in calculating Weierstrass zeta function

[Zou-Bo]

I am doing the calculation that follows the equation (13) and (14) in https://doi.org/10.1063/1.5042618
I use tau = 1+sqrt(3)*im. It turns out that it gives incorrect results when tau varies slightly.

[chakravala]

Looks like it might be an issue with a multiplication by im somewhere? The values seem to be swapped by a multiplication of the imaginary unit