W durkin davies benchmarks 2

This is now ready. It supersedes and closes #170 .
What is the reason for continuing the other computation? To force it to this specific solution? Is this done in davies et al too?