Use modified square root routine

[wbl]

This routine returns either sqrt(x) or sqrt(h*x), where h is a known constant. The additional cost is an extra iteration of the loop in Tonelli-Shanks.