[Numpy-discussion] how to use lu decomposition in lapack library?