Return to sqrt.c CVS log

Up to [MW Coherent from dump] / coherent / b / lib / libm

Annotation of coherent/b/lib/libm/sqrt.c, revision 1.1.1.1

1.1       root        1: /*
                      2:  * Square root function.
                      3:  */
                      4: #include <math.h>
                      5: 
                      6: #if    EMU87
                      7: #include "emumath.h"
                      8: #endif
                      9: 
                     10: double
                     11: sqrt(x)
                     12: double x;
                     13: {
                     14:        double s;
                     15:        int i;
                     16:        register int n;
                     17: 
                     18:        if (x < 0.0) {
                     19:                errno = EDOM;
                     20:                return (0.0);
                     21:        }
                     22:        if (x == 0.)
                     23:                return(0.);
                     24:        n = L2L2P;
                     25:        /* The assignment in the following line is to avoid compiler bug. */
                     26:        s = frexp(x, &i);
                     27:        s = ldexp(1.0, i/2);
                     28:        do {
                     29:                s = (s + x/s) / 2.0;
                     30:        } while (--n);
                     31:        return (s);
                     32: }

This archive runs on limited infrastructure. Preserving old code on modern bandwidth. Automated agents are requested to crawl responsibly.