Efficient filtering and linear algebra routines for multidimensional arrays

An efficient sampler for discrete random variables

Compute `real(dot(x, y))` efficiently.

Implementation of Takahashi''s inverse subset algorithm for sparse matrices

-