Started to port the code from the sharpness paper (PLoS CB 2013)