smZ function

Generate Sparse Matrix Z for penalized spline smoothing