pub fn dpsifn( x: c_double, n: c_int, k: c_int, m: c_int, ans: *mut c_double, nz: *mut c_int, ierr: *mut c_int, )