Differences in Home vs CAS

[parisse]

The line t := (TRN(r)*r) / (TRN(r)*A*r); should be rewritten as t := (TRN(r)*r) * (TRN(r)*A*r)^(-1)
Try debug(conjugate_gradient(...)), this will spot the line where the error occurs.