In PyKeops package, there is no available formula for Matrix-Matrix multiplication. Instead, they have implemented something similar to numpy.tensordot
PyKeops
numpy.tensordot