Convergence issue when using complex current density at an angle.

[Kaeryv]

Hi,
First, thank you for sharing this software !
I encounter an issue when trying to input a current density with nonzero Jx and Jy where Jx and Jy are out of phase:
Jx = b*1j and Jy = a for example. The objective is to produce a circular Right-Handed or Left-Handed polarization.

Maybe I'm doing something wrong I tried to modify the number of iterations and the accuracy (eps).

Here is the code for b:

depth = 137
b = np.zeros((172, 172, depth), dtype=np.complex128)
b = b.copy(), b.copy(), b.copy()

ax = 1j / np.sqrt(2)
ay = 1 / np.sqrt(2)
# The plane at z=49 is a source
b[0][:,:,41+8] = -1j * omega * ax
b[1][:,:,41+8] = -1j * omega * ay

PMLs of 40 blocks in each direction are used around the whole structure.

I will link more information soon.
Have a nice day