|
1 | 1 | import sys |
2 | 2 | import os |
3 | | -from scipy.integrate import trapz, quad |
| 3 | +from scipy.integrate import trapezoid, quad |
4 | 4 | from scipy.stats import iqr, norm |
5 | 5 | from scipy.interpolate import interp2d, SmoothBivariateSpline, Rbf, RectBivariateSpline |
6 | 6 | from scipy.fftpack import fft, ifft |
@@ -1111,7 +1111,9 @@ def fluxdens_to_fluxsum(R, I, axisratio): |
1111 | 1111 | S = np.zeros(len(R)) |
1112 | 1112 | S[0] = I[0] * np.pi * axisratio[0] * (R[0] ** 2) |
1113 | 1113 | for i in range(1, len(R)): |
1114 | | - S[i] = trapz(2 * np.pi * I[: i + 1] * R[: i + 1] * axisratio[: i + 1], R[: i + 1]) + S[0] |
| 1114 | + S[i] = ( |
| 1115 | + trapezoid(2 * np.pi * I[: i + 1] * R[: i + 1] * axisratio[: i + 1], R[: i + 1]) + S[0] |
| 1116 | + ) |
1115 | 1117 | return S |
1116 | 1118 |
|
1117 | 1119 |
|
@@ -1223,7 +1225,7 @@ def Fmode_fluxdens_to_fluxsum(R, I, parameters, A=None): |
1223 | 1225 | S[0] = I[0] * Aq[0] |
1224 | 1226 | Adiff = np.array([Aq[0]] + list(Aq[1:] - Aq[:-1])) |
1225 | 1227 | for i in range(1, len(R)): |
1226 | | - S[i] = trapz(I[: i + 1] * Adiff[: i + 1], R[: i + 1]) + S[0] |
| 1228 | + S[i] = trapezoid(I[: i + 1] * Adiff[: i + 1], R[: i + 1]) + S[0] |
1227 | 1229 | return S |
1228 | 1230 |
|
1229 | 1231 |
|
@@ -1343,7 +1345,7 @@ def SBprof_to_COG(R, SB, parameters): |
1343 | 1345 | # R = arcsec_to_pc(np.array(R), D) |
1344 | 1346 | # # Integrate up to each radius in the profile |
1345 | 1347 | # for i in range(1,len(R)): |
1346 | | - # m[i] = abs_mag_to_app_mag(L_to_mag(trapz(2*np.pi*I[:i+1]*R[:i+1]*axisratio[:i+1],R[:i+1]) + \ |
| 1348 | + # m[i] = abs_mag_to_app_mag(L_to_mag(trapezoid(2*np.pi*I[:i+1]*R[:i+1]*axisratio[:i+1],R[:i+1]) + \ |
1347 | 1349 | # mag_to_L(app_mag_to_abs_mag(m[0], D), band), band),D) |
1348 | 1350 | # elif method == 1: |
1349 | 1351 | # # Compute the starting point assuming constant SB within first isophote |
|
0 commit comments