double sin(double x) {
  double Result = x, Fac = 1.0, Xn = x, Precious = x;
  int n = 1, sign = 1;
  while (Precious > 1e-6) {
    n = n + 1;
    Fac = Fac * n * (n + 1);
    n = n + 1;
    Xn *= x * x;
    sign = -sign;
    Precious = Xn / Fac;
    Result = sign > 0 ? Result + Precious : Result - Precious;
  }
  return Result;
}


x86 assembly (clang, -O0):
<sin>:
push   %rbp
mov    %rsp,%rbp
movsd  %xmm0,-0x8(%rbp)
movsd  -0x8(%rbp),%xmm0
movsd  %xmm0,-0x10(%rbp)
movsd  0xee5(%rip),%xmm0
movsd  %xmm0,-0x18(%rbp)
movsd  -0x8(%rbp),%xmm0
movsd  %xmm0,-0x20(%rbp)
movsd  -0x8(%rbp),%xmm0
movsd  %xmm0,-0x28(%rbp)
movl   $0x1,-0x2c(%rbp)
movl   $0x1,-0x30(%rbp)
movsd  -0x28(%rbp),%xmm0
movsd  0xeb9(%rip),%xmm1
ucomisd %xmm1,%xmm0
jbe    11f3 <sin+0xf3>
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
movsd  -0x18(%rbp),%xmm0
cvtsi2sdl -0x2c(%rbp),%xmm1
mulsd  %xmm1,%xmm0
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
cvtsi2sd %eax,%xmm1
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0x18(%rbp)
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
movsd  -0x8(%rbp),%xmm0
mulsd  -0x8(%rbp),%xmm0
mulsd  -0x20(%rbp),%xmm0
movsd  %xmm0,-0x20(%rbp)
xor    %eax,%eax
sub    -0x30(%rbp),%eax
mov    %eax,-0x30(%rbp)
movsd  -0x20(%rbp),%xmm0
divsd  -0x18(%rbp),%xmm0
movsd  %xmm0,-0x28(%rbp)
cmpl   $0x0,-0x30(%rbp)
jle    11d5 <sin+0xd5>
movsd  -0x10(%rbp),%xmm0
addsd  -0x28(%rbp),%xmm0
movsd  %xmm0,-0x38(%rbp)
jmp    11e4 <sin+0xe4>
movsd  -0x10(%rbp),%xmm0
subsd  -0x28(%rbp),%xmm0
movsd  %xmm0,-0x38(%rbp)
movsd  -0x38(%rbp),%xmm0
movsd  %xmm0,-0x10(%rbp)
jmp    1142 <sin+0x42>
movsd  -0x10(%rbp),%xmm0
pop    %rbp
ret

x86 assembly (clang, -O2):
<sin>:
ucomisd 0xef8(%rip),%xmm0
jbe    119a <sin+0x9a>
movapd %xmm0,%xmm2
mulsd  %xmm0,%xmm2
movsd  0xeea(%rip),%xmm3
mov    $0xffffffff,%eax
mov    $0x1,%esi
mov    $0x2,%ecx
movapd 0xedb(%rip),%xmm4
movsd  0xec3(%rip),%xmm5
movapd %xmm0,%xmm1
jmp    1161 <sin+0x61>
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
addsd  %xmm7,%xmm1
add    $0x2,%ecx
ucomisd %xmm5,%xmm6
mov    %eax,%esi
mov    %edx,%eax
jbe    119b <sin+0x9b>
mov    %esi,%edx
xorps  %xmm6,%xmm6
cvtsi2sd %ecx,%xmm6
mulsd  %xmm3,%xmm6
lea    0x1(%rcx),%esi
xorps  %xmm3,%xmm3
cvtsi2sd %esi,%xmm3
mulsd  %xmm2,%xmm0
mulsd  %xmm6,%xmm3
movapd %xmm0,%xmm6
divsd  %xmm3,%xmm6
movapd %xmm6,%xmm7
test   %eax,%eax
jg     1150 <sin+0x50>
movapd %xmm6,%xmm7
xorpd  %xmm4,%xmm7
jmp    1150 <sin+0x50>
ret
movapd %xmm1,%xmm0
ret

LLVM IR (-O0):
define double @sin(double %0) {
%2 = alloca double
%3 = alloca double
%4 = alloca double
%5 = alloca double
%6 = alloca double
%7 = alloca i32
%8 = alloca i32
store double %0, double* %2
%9 = load double, double* %2
store double %9, double* %3
store double 1.000000e+00, double* %4
%10 = load double, double* %2
store double %10, double* %5
%11 = load double, double* %2
store double %11, double* %6
store i32 1, i32* %7
store i32 1, i32* %8
br label %12
12: ; preds = %48, %1
%13 = load double, double* %6
%14 = fcmp ogt double %13, 0x3EB0C6F7A0B5ED8D
br i1 %14, label %15, label %50
15: ; preds = %12
%16 = load i32, i32* %7
%17 = add nsw i32 %16, 1
store i32 %17, i32* %7
%18 = load double, double* %4
%19 = load i32, i32* %7
%20 = sitofp i32 %19 to double
%21 = fmul double %18, %20
%22 = load i32, i32* %7
%23 = add nsw i32 %22, 1
%24 = sitofp i32 %23 to double
%25 = fmul double %21, %24
store double %25, double* %4
%26 = load i32, i32* %7
%27 = add nsw i32 %26, 1
store i32 %27, i32* %7
%28 = load double, double* %2
%29 = load double, double* %2
%30 = fmul double %28, %29
%31 = load double, double* %5
%32 = fmul double %31, %30
store double %32, double* %5
%33 = load i32, i32* %8
%34 = sub nsw i32 0, %33
store i32 %34, i32* %8
%35 = load double, double* %5
%36 = load double, double* %4
%37 = fdiv double %35, %36
store double %37, double* %6
%38 = load i32, i32* %8
%39 = icmp sgt i32 %38, 0
br i1 %39, label %40, label %44
40: ; preds = %15
%41 = load double, double* %3
%42 = load double, double* %6
%43 = fadd double %41, %42
br label %48
44: ; preds = %15
%45 = load double, double* %3
%46 = load double, double* %6
%47 = fsub double %45, %46
br label %48
48: ; preds = %44, %40
%49 = phi double [ %43, %40 ], [ %47, %44 ]
store double %49, double* %3
br label %12
50: ; preds = %12
%51 = load double, double* %3
ret double %51
}

LLVM IR (-O2):
define double @sin(double %0) {
%2 = fcmp ogt double %0, 0x3EB0C6F7A0B5ED8D
br i1 %2, label %3, label %25
3: ; preds = %1
%4 = fmul double %0, %0
br label %5
5: ; preds = %3, %5
%6 = phi i32 [ 1, %3 ], [ %7, %5 ]
%7 = phi i32 [ -1, %3 ], [ %6, %5 ]
%8 = phi i32 [ 1, %3 ], [ %15, %5 ]
%9 = phi double [ %0, %3 ], [ %18, %5 ]
%10 = phi double [ 1.000000e+00, %3 ], [ %17, %5 ]
%11 = phi double [ %0, %3 ], [ %23, %5 ]
%12 = add nuw nsw i32 %8, 1
%13 = sitofp i32 %12 to double
%14 = fmul double %10, %13
%15 = add nuw nsw i32 %8, 2
%16 = sitofp i32 %15 to double
%17 = fmul double %14, %16
%18 = fmul double %4, %9
%19 = fdiv double %18, %17
%20 = icmp sgt i32 %7, 0
%21 = fneg double %19
%22 = select i1 %20, double %19, double %21
%23 = fadd double %11, %22
%24 = fcmp ogt double %19, 0x3EB0C6F7A0B5ED8D
br i1 %24, label %5, label %25
25: ; preds = %5, %1
%26 = phi double [ %0, %1 ], [ %23, %5 ]
ret double %26
}