WIP deformation radius calculations

need to refine sparse diffmat function

need to test with irregular data intervals
add tests for sparse first derivative