This commit starts introducing `massiv` in the codebase,
initially for simple functions like `termDivNan`. The main
goal is to extend the linear algebra toolkit up to the
point where we can implement `distributional` in terms of
`massive` and measure its performance.

The main idea is trying to refactor/improve the existing linear algebra
functions one function at the time, using reference implementations and
benchmarks along the way.