powered by
batch_normalization_differential(delta_y, mu, sigma_2, x, x_hat, y, gamma, beta, epsilon = exp(-12), with_BN = T)