----------------------------------------------------------- perldl> help minv minv Compute inverse of a general square matrix using LU factorization. Supports inplace and threading. Uses getrf and getri or cgetrf and cgetri from Lapack and return "inverse, info" in array context. PDL(inv) = minv(PDL) my $a = random(10,10); my $inv = minv($a); Docs from /usr/lib/perl5/PDL/LinearAlgebra.pm perldl> help msolve msolve Solve linear system of equations using LU decomposition. A * X = B Returns X in scalar context else X, LU, pivot vector and info. B is overwrited by X if its inplace flag is set. Supports threading. Uses gesv or cgesv from Lapack. (PDL(X), (PDL(LU), PDL(pivot), PDL(info))) = msolve(PDL(A), PDL(B)) my $a = random(5,5); my $b = random(10,5); my $X = msolve($a, $b); Docs from /usr/lib/perl5/PDL/LinearAlgebra.pm