BUG: Force observations to NumPy arrays in Stochastic.observe.
BUG: Do not sum over all axes in utils.linalg.chol_logdet

Given an array of Cholesky factors, the function now computes the log
determinant for each factor without summing them together. This caused
bugs for the lower bound evaluation in GaussianArrayARD.