emsig / empymod

Full 3D electromagnetic modeller for 1D VTI media
https://empymod.emsig.xyz
Apache License 2.0
85 stars 22 forks source link

BUG: Issues with `numba>=0.54.0` #138

Closed prisae closed 2 years ago

prisae commented 3 years ago

There seems to be a problem with numba>=0.54.0 and empymod.

I have not yet figured out what the issue is. For now, if you run into problems with empymod, please try to downgrade to numba<0.54.0 and see if the issue persists.

prisae commented 2 years ago

MWE

from empymod import dipole
dipole([0, 0, 1], [1, 0, -1], 0, [1, 2], 1)

Works with numba<0.54, fails with a ZeroDivisionError with numba>=0.54.

prisae commented 2 years ago

It seems to come from fastmath=True for empymod.kernel.wavenumber, which seem to have changed how it handles division by zero or similar.