double sqrt(x)
double x;
{
  double extra;
  unsigned char *px;
  unsigned char *pextra;
  int i;
  if (x == 0.0)
    return 0.0;
  if (x < 0.0)
    return 0.0;
  px = &x;
  pextra = &extra;
  extra = 0.707;
  pextra[5] = (px[5] >> 1) ^ 64;
  i = 7;
  while (--i) {
    extra += x / extra;
    --pextra[5];
  }
  return extra;
}


x86 assembly (clang, -O0):
<sqrt>:
push   %rbp
mov    %rsp,%rbp
movsd  %xmm0,-0x10(%rbp)
movsd  -0x10(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    112e <sqrt+0x2e>
jp     112e <sqrt+0x2e>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x8(%rbp)
jmp    11c6 <sqrt+0xc6>
xorps  %xmm0,%xmm0
ucomisd -0x10(%rbp),%xmm0
jbe    1149 <sqrt+0x49>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x8(%rbp)
jmp    11c6 <sqrt+0xc6>
lea    -0x10(%rbp),%rax
mov    %rax,-0x20(%rbp)
lea    -0x18(%rbp),%rax
mov    %rax,-0x28(%rbp)
movsd  0xe9f(%rip),%xmm0
movsd  %xmm0,-0x18(%rbp)
mov    -0x20(%rbp),%rax
movzbl 0x5(%rax),%eax
sar    $0x1,%eax
xor    $0x40,%eax
mov    %al,%cl
mov    -0x28(%rbp),%rax
mov    %cl,0x5(%rax)
movl   $0x7,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
add    $0xffffffff,%eax
mov    %eax,-0x2c(%rbp)
cmp    $0x0,%eax
je     11bc <sqrt+0xbc>
movsd  -0x10(%rbp),%xmm0
divsd  -0x18(%rbp),%xmm0
addsd  -0x18(%rbp),%xmm0
movsd  %xmm0,-0x18(%rbp)
mov    -0x28(%rbp),%rax
mov    0x5(%rax),%cl
add    $0xff,%cl
mov    %cl,0x5(%rax)
jmp    1184 <sqrt+0x84>
movsd  -0x18(%rbp),%xmm0
movsd  %xmm0,-0x8(%rbp)
movsd  -0x8(%rbp),%xmm0
pop    %rbp
ret

x86 assembly (clang, -O2):
<sqrt>:
xorpd  %xmm1,%xmm1
ucomisd %xmm0,%xmm1
jae    1217 <sqrt+0x117>
movq   %xmm0,%rax
shr    %rax
movabs $0x7f0000000000,%rcx
and    %rax,%rcx
movabs $0x3fe640be76c8b439,%rax
40 e6 3f
xor    %rcx,%rax
movq   %rax,%xmm1
movapd %xmm0,%xmm2
divsd  %xmm1,%xmm2
addsd  %xmm1,%xmm2
movq   %xmm2,%rdx
movabs $0xff0000000000,%rax
lea    (%rdx,%rax,1),%rsi
and    %rax,%rsi
movabs $0xffff00ffffffffff,%rcx
and    %rcx,%rdx
or     %rsi,%rdx
movq   %rdx,%xmm1
movapd %xmm0,%xmm2
divsd  %xmm1,%xmm2
addsd  %xmm1,%xmm2
movq   %xmm2,%rdx
lea    (%rdx,%rax,1),%rsi
and    %rax,%rsi
and    %rcx,%rdx
or     %rsi,%rdx
movq   %rdx,%xmm1
movapd %xmm0,%xmm2
divsd  %xmm1,%xmm2
addsd  %xmm1,%xmm2
movq   %xmm2,%rdx
lea    (%rdx,%rax,1),%rsi
and    %rax,%rsi
and    %rcx,%rdx
or     %rsi,%rdx
movq   %rdx,%xmm1
movapd %xmm0,%xmm2
divsd  %xmm1,%xmm2
addsd  %xmm1,%xmm2
movq   %xmm2,%rdx
lea    (%rdx,%rax,1),%rsi
and    %rax,%rsi
and    %rcx,%rdx
or     %rsi,%rdx
movq   %rdx,%xmm1
movapd %xmm0,%xmm2
divsd  %xmm1,%xmm2
addsd  %xmm1,%xmm2
movq   %xmm2,%rdx
lea    (%rdx,%rax,1),%rsi
and    %rax,%rsi
and    %rcx,%rdx
or     %rsi,%rdx
movq   %rdx,%xmm1
divsd  %xmm1,%xmm0
addsd  %xmm1,%xmm0
movq   %xmm0,%rdx
lea    (%rdx,%rax,1),%rsi
and    %rax,%rsi
and    %rcx,%rdx
or     %rsi,%rdx
movq   %rdx,%xmm1
movapd %xmm1,%xmm0
ret

LLVM IR (-O0):
define double @sqrt(double %0) {
%2 = alloca double
%3 = alloca double
%4 = alloca double
%5 = alloca i8*
%6 = alloca i8*
%7 = alloca i32
store double %0, double* %3
%8 = load double, double* %3
%9 = fcmp oeq double %8, 0.000000e+00
br i1 %9, label %10, label %11
10: ; preds = %1
store double 0.000000e+00, double* %2
br label %43
11: ; preds = %1
%12 = load double, double* %3
%13 = fcmp olt double %12, 0.000000e+00
br i1 %13, label %14, label %15
14: ; preds = %11
store double 0.000000e+00, double* %2
br label %43
15: ; preds = %11
%16 = bitcast double* %3 to i8*
store i8* %16, i8** %5
%17 = bitcast double* %4 to i8*
store i8* %17, i8** %6
store double 7.070000e-01, double* %4
%18 = load i8*, i8** %5
%19 = getelementptr inbounds i8, i8* %18, i64 5
%20 = load i8, i8* %19
%21 = zext i8 %20 to i32
%22 = ashr i32 %21, 1
%23 = xor i32 %22, 64
%24 = trunc i32 %23 to i8
%25 = load i8*, i8** %6
%26 = getelementptr inbounds i8, i8* %25, i64 5
store i8 %24, i8* %26
store i32 7, i32* %7
br label %27
27: ; preds = %31, %15
%28 = load i32, i32* %7
%29 = add nsw i32 %28, -1
store i32 %29, i32* %7
%30 = icmp ne i32 %29, 0
br i1 %30, label %31, label %41
31: ; preds = %27
%32 = load double, double* %3
%33 = load double, double* %4
%34 = fdiv double %32, %33
%35 = load double, double* %4
%36 = fadd double %35, %34
store double %36, double* %4
%37 = load i8*, i8** %6
%38 = getelementptr inbounds i8, i8* %37, i64 5
%39 = load i8, i8* %38
%40 = add i8 %39, -1
store i8 %40, i8* %38
br label %27
41: ; preds = %27
%42 = load double, double* %4
store double %42, double* %2
br label %43
43: ; preds = %41, %14, %10
%44 = load double, double* %2
ret double %44
}

LLVM IR (-O2):
define double @sqrt(double %0) {
%2 = fcmp ugt double %0, 0.000000e+00
br i1 %2, label %3, label %57
3: ; preds = %1
%4 = bitcast double %0 to i64
%5 = lshr i64 %4, 1
%6 = and i64 %5, 139637976727552
%7 = xor i64 %6, 4604438855813739577
%8 = bitcast i64 %7 to double
%9 = fdiv double %0, %8
%10 = fadd double %9, %8
%11 = bitcast double %10 to i64
%12 = add i64 %11, 280375465082880
%13 = and i64 %12, 280375465082880
%14 = and i64 %11, -280375465082881
%15 = or i64 %13, %14
%16 = bitcast i64 %15 to double
%17 = fdiv double %0, %16
%18 = fadd double %17, %16
%19 = bitcast double %18 to i64
%20 = add i64 %19, 280375465082880
%21 = and i64 %20, 280375465082880
%22 = and i64 %19, -280375465082881
%23 = or i64 %21, %22
%24 = bitcast i64 %23 to double
%25 = fdiv double %0, %24
%26 = fadd double %25, %24
%27 = bitcast double %26 to i64
%28 = add i64 %27, 280375465082880
%29 = and i64 %28, 280375465082880
%30 = and i64 %27, -280375465082881
%31 = or i64 %29, %30
%32 = bitcast i64 %31 to double
%33 = fdiv double %0, %32
%34 = fadd double %33, %32
%35 = bitcast double %34 to i64
%36 = add i64 %35, 280375465082880
%37 = and i64 %36, 280375465082880
%38 = and i64 %35, -280375465082881
%39 = or i64 %37, %38
%40 = bitcast i64 %39 to double
%41 = fdiv double %0, %40
%42 = fadd double %41, %40
%43 = bitcast double %42 to i64
%44 = add i64 %43, 280375465082880
%45 = and i64 %44, 280375465082880
%46 = and i64 %43, -280375465082881
%47 = or i64 %45, %46
%48 = bitcast i64 %47 to double
%49 = fdiv double %0, %48
%50 = fadd double %49, %48
%51 = bitcast double %50 to i64
%52 = add i64 %51, 280375465082880
%53 = and i64 %52, 280375465082880
%54 = and i64 %51, -280375465082881
%55 = or i64 %53, %54
%56 = bitcast i64 %55 to double
br label %57
57: ; preds = %3, %1
%58 = phi double [ 0.000000e+00, %1 ], [ %56, %3 ]
ret double %58
}