slatdf(3) on Linux
The text bf0aa1793d3b, as Debian 3.1 (Sarge) shipped it.
use the LU factorization of the n-by-n matrix Z computed by SGETC2 and computes a contribution to the reciprocal Dif-estimate by solving Z * x = b for x, and choosing the r.h.s