[Numpy-discussion] determinant of a scalar not handled

Skipper Seabold jsseabold at gmail.com
Tue Jul 27 08:51:08 EDT 2010


On Mon, Jul 26, 2010 at 10:05 PM, Alan G Isaac <aisaac at american.edu> wrote:
> But I am still confused about the use case.
> What is the scalar- (or 1d-array-) returning procedure
> invokedbefore taking the determinant?
>

Recently I ran into this trying to make the log-likelihood of a
multivariate and univariate autoregressive process use the same
function.  One has log(sigma_scalar) and one calls for
logdet(Sigma_matrix).  I also ran in to again yesterday working on the
Kalman filter, depending on the process being modeled and how the user
writes a function if the needed coefficient arrays depend on
parameters.  To be more general, I have to put in atleast_2d, even
though these checks are really in slogdet.

Skipper



More information about the NumPy-Discussion mailing list