replace invSqrt() with 1/sqrt()

* the original implementation used 1/sqrt() instead of the fast
 invSqrt() function
* invSqrt() causes the estimate to diverge under constant input
* this commit fixes #17