import empymod
import numpy as np
empymod.dipole(
src=[0, 0, 0], rec=[5000, 0, 0], depth=[], res=1,
freqtime=np.logspace(-2, 2, 201),
signal=0, ft='fftlog', ftarg=[5, [-2, 1]]
)
---------------------------------------------------------------------------
ValueError Traceback (most recent call last)
<ipython-input-1-88dbe35cac4d> in <module>
3
4 empymod.dipole([0, 0, 0], [5000, 0, 0], [], 1, np.logspace(-2, 2, 201),
----> 5 signal=0, ft='fftlog', ftarg=[5, [-2, 1]])
~/anaconda3/lib/python3.7/site-packages/empymod/model.py in dipole(src, rec, depth, res, freqtime, signal, ab, aniso, epermH, epermV, mpermH, mpermV, xdirect, ht, htarg, ft, ftarg, opt, loop, verb)
947 # Do f->t transform if required
948 if signal is not None:
--> 949 EM, conv = tem(EM, off, freq, time, signal, ft, ftarg)
950
951 # In case of QWE/QUAD, print Warning if not converged
~/anaconda3/lib/python3.7/site-packages/empymod/model.py in tem(fEM, off, freq, time, signal, ft, ftarg, conv)
1587 tEM = np.zeros((time.size, off.size))
1588 for i in range(off.size):
-> 1589 out = getattr(transform, ft)(fEM[:, i]*fact, time, freq, ftarg)
1590 tEM[:, i] += out[0]
1591 conv *= out[1]
~/anaconda3/lib/python3.7/site-packages/empymod/transform.py in fftlog(fEM, time, freq, ftarg)
723 ar = a[2*m-1]
724 ai = a[2*m]
--> 725 a[2*m-1] = ar*argcos[:-1] - ai*argsin[:-1]
726 a[2*m] = ar*argsin[:-1] + ai*argcos[:-1]
727 # problematical last element, for even n
ValueError: operands could not be broadcast together with shapes (17,) (16,)