c_axi[0][0, 0] = 1.9
cr_axi = {0: np.zeros((1, 2 * l + 1))}
cr_axi[0][0, 0] = 1.9
b = pd.zeros(1, dtype=complex)
br = pdr.zeros(1)
lfc.set_positions(spos_ac)
lfc.add(b, c_axi)
lfcr.set_positions(spos_ac)
lfcr.add(br, cr_axi)
a = pd.ifft(b)
ar = pdr.ifft(br)
equal(abs(a-ar).max(), 0, 1e-14)
if l == 0:
a = a[:, ::-1].copy()
b0 = pd.fft(a)
br0 = pdr.fft(a.real)
lfc.integrate(b0, c_axi)
lfcr.integrate(br0, cr_axi)
assert abs(c_axi[0][0]-cr_axi[0][0]).max() < 1e-14
c_axiv = {0: np.zeros((1, 2 * l + 1, 3), complex)}
cr_axiv = {0: np.zeros((1, 2 * l + 1, 3))}
lfc.derivative(b0, c_axiv)
lfcr.derivative(br0, cr_axiv)
assert abs(c_axiv[0][0]-cr_axiv[0][0]).max() < 1e-14
c_axi[0][0, 0] = 1.9
cr_axi = {0: np.zeros((1, 2 * l + 1))}
cr_axi[0][0, 0] = 1.9
b = pd.zeros(1, dtype=complex)
br = pdr.zeros(1)
lfc.set_positions(spos_ac)
lfc.add(b, c_axi)
lfcr.set_positions(spos_ac)
lfcr.add(br, cr_axi)
a = pd.ifft(b)
ar = pdr.ifft(br)
equal(abs(a-ar).max(), 0, 1e-14)
if l == 0:
a = a[:, ::-1].copy()
b0 = pd.fft(a)
br0 = pdr.fft(a.real)
lfc.integrate(b0, c_axi)
lfcr.integrate(br0, cr_axi)
assert abs(c_axi[0][0]-cr_axi[0][0]).max() < 1e-14
c_axiv = {0: np.zeros((1, 2 * l + 1, 3), complex)}
cr_axiv = {0: np.zeros((1, 2 * l + 1, 3))}
lfc.derivative(b0, c_axiv)
lfcr.derivative(br0, cr_axiv)
assert abs(c_axiv[0][0]-cr_axiv[0][0]).max() < 1e-14
