Logspline density estimation routines