Commit 8f64641d authored by Roland Haas's avatar Roland Haas
Browse files

POWER: make work with numpy older than 1.17

parent 2051efc7
......@@ -520,8 +520,8 @@ def POWER(sim_path, radii, modes):
for i in range(1, amp_extrapolation_order+1):
A_amp = np.column_stack((A_amp, np.power(radii, -1*i)))
b_phase = np.empty_like(radii, shape=(len(radii), len(t)))
b_amp = np.empty_like(radii, shape=(len(radii), len(t)))
b_phase = np.empty(dtype=radii.dtype, shape=(len(radii), len(t)))
b_amp = np.empty(dtype=radii.dtype, shape=(len(radii), len(t)))
for j in range(len(radii)):
b_phase[j] = phase[j][:, 1]
b_amp[j] = amp[j][:, 1]
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment