INNER CODE UNIT · Python

lnCKCD

QianyeSu/Skyborn · src/skyborn/calc/GPI/core.py:518

    lnCKCD = float(np.log(CKCD))

    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,

View source record →