#include "mathops.h" #ifdef TEST #include <stdlib.h> int main(void) { unsigned u; for(u=0; u<65536; u++) { unsigned s = u*u; unsigned root = ff_sqrt(s); unsigned root_m1 = ff_sqrt(s-1); if (s && root != u) { fprintf(stderr, "ff_sqrt failed at %u with %u\n", s, root); return 1; } if (u && root_m1 != u - 1) { fprintf(stderr, "ff_sqrt failed at %u with %u\n", s, root); return 1; } } return 0; } #endif /* TEST */