Support integral square root with unsigned fixed point output in llvm-libc

[PiJoules]

sqrtuk and co work really well but we don't have sqrtulk at the moment, so none of the fixed point sqrt functions will work out of the box for values larger than 2^16. Even though it's not part of ISO 18037, it would be nice to have an integral sqrt as an option.