Hello I find your code in this link,
I know it is not the latest version, but because I am unable to find the corresponding part in PDL-2.4.9, I have to report the bug in the 2.4.3 version as the link above pointing to.
One sentence in the link
should be replaced by
for matrix calculation, as c is the coefficient of the unknown variables x. I have verified it in my own code.