double _sqt(x)
double x;
{
  extern double _fef();
  int exponent;
  double val;
  if (x <= 0) {
    if (x < 0)
      error(3);
    return 0;
  }
  val = _fef(x, &exponent);
  if (exponent & 1) {
    exponent--;
    val *= 2;
  }
  val = ldexp(val + 1.0, exponent / 2 - 1);
  for (exponent = 5 - 1; exponent >= 0; exponent--) {
    val = (val + x / val) / 2.0;
  }
  return val;
}


x86 assembly (clang, -O0):
<_sqt>:
push   %rbp
mov    %rsp,%rbp
sub    $0x20,%rsp
movsd  %xmm0,-0x10(%rbp)
xorps  %xmm0,%xmm0
ucomisd -0x10(%rbp),%xmm0
jb     1172 <_sqt+0x42>
xorps  %xmm0,%xmm0
ucomisd -0x10(%rbp),%xmm0
jbe    1165 <_sqt+0x35>
mov    $0x3,%edi
mov    $0x0,%al
call   1050 <error@plt>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x8(%rbp)
jmp    1225 <_sqt+0xf5>
movsd  -0x10(%rbp),%xmm0
lea    -0x14(%rbp),%rdi
mov    $0x1,%al
call   1040 <_fef@plt>
movsd  %xmm0,-0x20(%rbp)
mov    -0x14(%rbp),%eax
and    $0x1,%eax
cmp    $0x0,%eax
je     11b1 <_sqt+0x81>
mov    -0x14(%rbp),%eax
add    $0xffffffff,%eax
mov    %eax,-0x14(%rbp)
movsd  0xe59(%rip),%xmm0
mulsd  -0x20(%rbp),%xmm0
movsd  %xmm0,-0x20(%rbp)
movsd  0xe4f(%rip),%xmm0
addsd  -0x20(%rbp),%xmm0
mov    -0x14(%rbp),%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
mov    %eax,%edi
sub    $0x1,%edi
call   1030 <ldexp@plt>
movsd  %xmm0,-0x20(%rbp)
movl   $0x4,-0x14(%rbp)
cmpl   $0x0,-0x14(%rbp)
jl     121b <_sqt+0xeb>
movsd  -0x20(%rbp),%xmm0
movsd  -0x10(%rbp),%xmm1
divsd  -0x20(%rbp),%xmm1
addsd  %xmm1,%xmm0
movsd  0xdfc(%rip),%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x20(%rbp)
mov    -0x14(%rbp),%eax
add    $0xffffffff,%eax
mov    %eax,-0x14(%rbp)
jmp    11df <_sqt+0xaf>
movsd  -0x20(%rbp),%xmm0
movsd  %xmm0,-0x8(%rbp)
movsd  -0x8(%rbp),%xmm0
add    $0x20,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<_sqt>:
sub    $0x18,%rsp
xorpd  %xmm3,%xmm3
ucomisd %xmm0,%xmm3
jae    11dc <_sqt+0xac>
lea    0xc(%rsp),%rdi
movsd  %xmm0,0x10(%rsp)
mov    $0x1,%al
call   1040 <_fef@plt>
mov    0xc(%rsp),%eax
test   $0x1,%al
je     1167 <_sqt+0x37>
add    $0xffffffff,%eax
mov    %eax,0xc(%rsp)
addsd  %xmm0,%xmm0
addsd  0xe91(%rip),%xmm0
mov    %eax,%edi
shr    $0x1f,%edi
add    %eax,%edi
sar    %edi
add    $0xffffffff,%edi
call   1030 <ldexp@plt>
movsd  0x10(%rsp),%xmm3
movapd %xmm3,%xmm1
divsd  %xmm0,%xmm1
addsd  %xmm0,%xmm1
movsd  0xe6e(%rip),%xmm0
mulsd  %xmm0,%xmm1
movapd %xmm3,%xmm2
divsd  %xmm1,%xmm2
addsd  %xmm1,%xmm2
mulsd  %xmm0,%xmm2
movapd %xmm3,%xmm1
divsd  %xmm2,%xmm1
addsd  %xmm2,%xmm1
mulsd  %xmm0,%xmm1
movapd %xmm3,%xmm2
divsd  %xmm1,%xmm2
addsd  %xmm1,%xmm2
mulsd  %xmm0,%xmm2
divsd  %xmm2,%xmm3
addsd  %xmm2,%xmm3
mulsd  %xmm0,%xmm3
jmp    11ee <_sqt+0xbe>
jbe    11ee <_sqt+0xbe>
mov    $0x3,%edi
xor    %eax,%eax
call   1050 <error@plt>
xorpd  %xmm3,%xmm3
movapd %xmm3,%xmm0
add    $0x18,%rsp
ret

LLVM IR (-O0):
define double @_sqt(double %0) {
%2 = alloca double
%3 = alloca double
%4 = alloca i32
%5 = alloca double
store double %0, double* %3
%6 = load double, double* %3
%7 = fcmp ole double %6, 0.000000e+00
br i1 %7, label %8, label %14
8: ; preds = %1
%9 = load double, double* %3
%10 = fcmp olt double %9, 0.000000e+00
br i1 %10, label %11, label %13
11: ; preds = %8
%12 = call i32 (i32, ...) bitcast (i32 (...)* @error to i32 (i32, ...)*)(i32 3)
br label %13
13: ; preds = %11, %8
store double 0.000000e+00, double* %2
br label %47
14: ; preds = %1
%15 = load double, double* %3
%16 = call double (double, i32*, ...) bitcast (double (...)* @_fef to double (double, i32*, ...)*)(double %15, i32* %4)
store double %16, double* %5
%17 = load i32, i32* %4
%18 = and i32 %17, 1
%19 = icmp ne i32 %18, 0
br i1 %19, label %20, label %25
20: ; preds = %14
%21 = load i32, i32* %4
%22 = add nsw i32 %21, -1
store i32 %22, i32* %4
%23 = load double, double* %5
%24 = fmul double %23, 2.000000e+00
store double %24, double* %5
br label %25
25: ; preds = %20, %14
%26 = load double, double* %5
%27 = fadd double %26, 1.000000e+00
%28 = load i32, i32* %4
%29 = sdiv i32 %28, 2
%30 = sub nsw i32 %29, 1
%31 = call double @ldexp(double %27, i32 %30)
store double %31, double* %5
store i32 4, i32* %4
br label %32
32: ; preds = %42, %25
%33 = load i32, i32* %4
%34 = icmp sge i32 %33, 0
br i1 %34, label %35, label %45
35: ; preds = %32
%36 = load double, double* %5
%37 = load double, double* %3
%38 = load double, double* %5
%39 = fdiv double %37, %38
%40 = fadd double %36, %39
%41 = fdiv double %40, 2.000000e+00
store double %41, double* %5
br label %42
42: ; preds = %35
%43 = load i32, i32* %4
%44 = add nsw i32 %43, -1
store i32 %44, i32* %4
br label %32
45: ; preds = %32
%46 = load double, double* %5
store double %46, double* %2
br label %47
47: ; preds = %45, %13
%48 = load double, double* %2
ret double %48
}

LLVM IR (-O2):
define double @_sqt(double %0) {
%2 = alloca i32
%3 = bitcast i32* %2 to i8*
%4 = fcmp ugt double %0, 0.000000e+00
br i1 %4, label %9, label %5
5: ; preds = %1
%6 = fcmp olt double %0, 0.000000e+00
br i1 %6, label %7, label %39
7: ; preds = %5
%8 = tail call i32 (i32, ...) bitcast (i32 (...)* @error to i32 (i32, ...)*)(i32 3)
br label %39
9: ; preds = %1
%10 = call double (double, i32*, ...) bitcast (double (...)* @_fef to double (double, i32*, ...)*)(double %0, i32* %2)
%11 = load i32, i32* %2
%12 = and i32 %11, 1
%13 = icmp eq i32 %12, 0
br i1 %13, label %17, label %14
14: ; preds = %9
%15 = add nsw i32 %11, -1
store i32 %15, i32* %2
%16 = fmul double %10, 2.000000e+00
br label %17
17: ; preds = %14, %9
%18 = phi i32 [ %15, %14 ], [ %11, %9 ]
%19 = phi double [ %16, %14 ], [ %10, %9 ]
%20 = fadd double %19, 1.000000e+00
%21 = sdiv i32 %18, 2
%22 = add nsw i32 %21, -1
%23 = call double @ldexp(double %20, i32 %22)
%24 = fdiv double %0, %23
%25 = fadd double %23, %24
%26 = fmul double %25, 5.000000e-01
%27 = fdiv double %0, %26
%28 = fadd double %26, %27
%29 = fmul double %28, 5.000000e-01
%30 = fdiv double %0, %29
%31 = fadd double %29, %30
%32 = fmul double %31, 5.000000e-01
%33 = fdiv double %0, %32
%34 = fadd double %32, %33
%35 = fmul double %34, 5.000000e-01
%36 = fdiv double %0, %35
%37 = fadd double %35, %36
%38 = fmul double %37, 5.000000e-01
br label %39
39: ; preds = %17, %5, %7
%40 = phi double [ 0.000000e+00, %7 ], [ 0.000000e+00, %5 ], [ %38, %17 ]
ret double %40
}