sf_tris sf_spline4_init (int nd /* data length */) /*< initialize a tridiagonal matrix for cubic splines >*/ { ... }