Add JIT-compiled functions to support multi-index multiplication

Multi-index multiplication is now supported by JIT-compiled functions for performance, namely the functions to compute cross-product-and-sum and to find unique indices from a lexicographically sorted indices. These are shown to accelerate the computation.

Additionally several tests are included specifically for these new functionalities.

This MR should resolve Issue #145 (closed).