compute_kmmc function

Compute KMMC statistics