Hello, Finally my double-precision-Pd efforts resulted in code decent enough to be useful in practice. It's all documented on: http://www.katjaas.nl/doubleprecision/doubleprecision.html