pymc.math.matrix_inverse

pymc.math.matrix_inverse = Blockwise(core_op=MatrixInverse,signature=(m,m)->(m,m))

Generalizes a core Op to work with batched dimensions.

TODO: Dispatch JAX (should be easy with the vectorize macro) TODO: Dispatch Numba TODO: C implementation? TODO: Fuse Blockwise?