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