move sqrt/rsqrt code to roots.rs