Given raw (pre-sigmoid) logits and binary target (0s and 1s), return the mean binary cross-entropy loss: for each element, `-(target * log(p) + (1 - target)

• log(1 - p))wherep = sigmoid(logit), averaged over all elements. Compute it in a numerically stable way directly from the logits rather than from sigmoid(logits)`.