PDL does have a binding to "the" solver library, LAPACK (linear algebra package): PDL::LinearAlgebra. The msolve function is almost certainly what you want. See the doc for PDL::MatrixOps/lu_backsub for a quick comparison of other linear algebra libraries in the PDLverse.