INNER CODE UNIT · Python
efficiency
QianyeSu/Skyborn · src/skyborn/calc/GPI/core.py:520
efficiency = (sst_k - t0_arr) / t0_arr
valid_eff = efficiency > 0
valid_pi = pi_arr > 0
lneff = np.full(efficiency.shape, np.nan, dtype=float)
lneff[valid_eff] = np.log(efficiency[valid_eff])
lnpi = np.full(pi_arr.shape, np.nan, dtype=float)
lnpi[valid_pi] = 2.0 * np.log(pi_arr[valid_pi])
lndiseq = np.full(pi_arr.shape, np.nan, dtype=float)
valid = valid_eff & valid_pi
lndiseq[valid] = lnpi[valid] - lneff[valid] - lnCKCD
return (
float(lnpi) if np.ndim(lnpi) == 0 else lnpi,
float(lneff) if np.ndim(lneff) == 0 else lneff,
float(lndiseq) if np.ndim(lndiseq) == 0 else lndiseq,