sparsify.with.simmelian function

Extract Nick et al's (2013) Simmelian backbone