diffpen()# liesel_ptm.diffpen(ncol, diff=2)[source]# A P-spline penalty matrix based on diff-order differences.