libk/src/math/sqrt.c


1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36

#include <math.h>

float sqrtf(float num) {
    if (num < 0)
        return NAN;

    if (num < 2)
        return num;
 
    float y = num;
    float z = (y + (num / y)) / 2;
 
    while (fabsf(y - z) >= 0.00001) {
        y = z;
        z = (y + (num / y)) / 2;
    }
    return z;
}

double sqrt(double num) {
    if (num < 0)
        return NAN;

    if (num < 2)
        return num;
 
    double y = num;
    double z = (y + (num / y)) / 2;
 
    while (fabs(y - z) >= 0.00001) {
        y = z;
        z = (y + (num / y)) / 2;
    }
    return z;
}