I don't want to fix `log(exposure)` in this PR. it's not backwards compatible and needs better design (for example merge `_offset_all = np.exp(exposure) + offset` to avoid the repeated sum).