float sin(rad)
float rad;
{
  float app;
  float diff;
  int inc = 1;
  while (rad > 2 * 3.14159265358979323846)
    rad -= 2 * 3.14159265358979323846;
  while (rad < -2 * 3.14159265358979323846)
    rad += 2 * 3.14159265358979323846;
  app = diff = rad;
  diff = (diff * (-(rad * rad))) / ((2.0 * inc) * (2.0 * inc + 1.0));
  app = app + diff;
  inc++;
  while (fabs(diff) >= 0.00001) {
    diff = (diff * (-(rad * rad))) / ((2.0 * inc) * (2.0 * inc + 1.0));
    app = app + diff;
    inc++;
  }
  return (app);
}


x86 assembly (clang, -O0):
<sin>:
push   %rbp
mov    %rsp,%rbp
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,-0x4(%rbp)
movl   $0x1,-0x10(%rbp)
movss  -0x4(%rbp),%xmm0
cvtss2sd %xmm0,%xmm0
movsd  0xedb(%rip),%xmm1
ucomisd %xmm1,%xmm0
jbe    1152 <sin+0x52>
movss  -0x4(%rbp),%xmm0
cvtss2sd %xmm0,%xmm0
movsd  0xec0(%rip),%xmm1
subsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,-0x4(%rbp)
jmp    1114 <sin+0x14>
jmp    1157 <sin+0x57>
movss  -0x4(%rbp),%xmm0
cvtss2sd %xmm0,%xmm1
movsd  0xea0(%rip),%xmm0
ucomisd %xmm1,%xmm0
jbe    1195 <sin+0x95>
movss  -0x4(%rbp),%xmm0
cvtss2sd %xmm0,%xmm0
movsd  0xe7d(%rip),%xmm1
addsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,-0x4(%rbp)
jmp    1157 <sin+0x57>
movss  -0x4(%rbp),%xmm0
movss  %xmm0,-0xc(%rbp)
movss  %xmm0,-0x8(%rbp)
movss  -0xc(%rbp),%xmm0
movss  -0x4(%rbp),%xmm1
mulss  -0x4(%rbp),%xmm1
movd   %xmm1,%eax
xor    $0x80000000,%eax
movd   %eax,%xmm1
mulss  %xmm1,%xmm0
cvtss2sd %xmm0,%xmm0
cvtsi2sdl -0x10(%rbp),%xmm2
movsd  0xe43(%rip),%xmm1
mulsd  %xmm2,%xmm1
cvtsi2sdl -0x10(%rbp),%xmm2
addsd  %xmm2,%xmm2
movsd  0xe26(%rip),%xmm3
addsd  %xmm3,%xmm2
mulsd  %xmm2,%xmm1
divsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,-0xc(%rbp)
movss  -0x8(%rbp),%xmm0
addss  -0xc(%rbp),%xmm0
movss  %xmm0,-0x8(%rbp)
mov    -0x10(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x10(%rbp)
movss  -0xc(%rbp),%xmm0
cvtss2sd %xmm0,%xmm0
movaps 0xe09(%rip),%xmm1
pand   %xmm1,%xmm0
movsd  0xded(%rip),%xmm1
ucomisd %xmm1,%xmm0
jb     12b5 <sin+0x1b5>
movss  -0xc(%rbp),%xmm0
movss  -0x4(%rbp),%xmm1
mulss  -0x4(%rbp),%xmm1
movd   %xmm1,%eax
xor    $0x80000000,%eax
movd   %eax,%xmm1
mulss  %xmm1,%xmm0
cvtss2sd %xmm0,%xmm0
cvtsi2sdl -0x10(%rbp),%xmm2
movsd  0xdaa(%rip),%xmm1
mulsd  %xmm2,%xmm1
cvtsi2sdl -0x10(%rbp),%xmm2
addsd  %xmm2,%xmm2
movsd  0xd8d(%rip),%xmm3
addsd  %xmm3,%xmm2
mulsd  %xmm2,%xmm1
divsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,-0xc(%rbp)
movss  -0x8(%rbp),%xmm0
addss  -0xc(%rbp),%xmm0
movss  %xmm0,-0x8(%rbp)
mov    -0x10(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x10(%rbp)
jmp    1217 <sin+0x117>
movss  -0x8(%rbp),%xmm0
pop    %rbp
ret

x86 assembly (clang, -O2):
<sin>:
cvtsd2ss %xmm0,%xmm0
cvtss2sd %xmm0,%xmm1
ucomisd 0xef0(%rip),%xmm1
jbe    1148 <sin+0x48>
movsd  0xeee(%rip),%xmm2
movsd  0xede(%rip),%xmm3
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
addsd  %xmm2,%xmm1
xorps  %xmm0,%xmm0
cvtsd2ss %xmm1,%xmm0
xorps  %xmm1,%xmm1
cvtss2sd %xmm0,%xmm1
ucomisd %xmm3,%xmm1
ja     1130 <sin+0x30>
movsd  0xeb8(%rip),%xmm2
ucomisd %xmm1,%xmm2
jbe    1178 <sin+0x78>
movsd  0xea2(%rip),%xmm3
xchg   %ax,%ax
addsd  %xmm3,%xmm1
xorps  %xmm0,%xmm0
cvtsd2ss %xmm1,%xmm0
xorps  %xmm1,%xmm1
cvtss2sd %xmm0,%xmm1
ucomisd %xmm1,%xmm2
ja     1160 <sin+0x60>
movaps 0xea1(%rip),%xmm1
xorps  %xmm0,%xmm1
mulss  %xmm0,%xmm1
movaps %xmm0,%xmm2
mulss  %xmm1,%xmm2
divss  0xeab(%rip),%xmm2
addss  %xmm2,%xmm0
movaps 0xe90(%rip),%xmm3
andps  %xmm2,%xmm3
cvtss2sd %xmm3,%xmm3
ucomisd 0xe61(%rip),%xmm3
jb     120e <sin+0x10e>
mov    $0x2,%eax
movsd  0xe5a(%rip),%xmm3
movaps 0xe6b(%rip),%xmm4
movsd  0xe43(%rip),%xmm5
nopl   (%rax)
mulss  %xmm1,%xmm2
cvtss2sd %xmm2,%xmm2
xorps  %xmm6,%xmm6
cvtsi2sd %eax,%xmm6
addsd  %xmm6,%xmm6
movapd %xmm6,%xmm7
addsd  %xmm3,%xmm7
mulsd  %xmm6,%xmm7
divsd  %xmm7,%xmm2
cvtsd2ss %xmm2,%xmm2
addss  %xmm2,%xmm0
add    $0x1,%eax
movaps %xmm2,%xmm6
andps  %xmm4,%xmm6
cvtss2sd %xmm6,%xmm6
ucomisd %xmm5,%xmm6
jae    11d0 <sin+0xd0>
ret

LLVM IR (-O0):
define float @sin(double %0) {
%2 = alloca float
%3 = alloca float
%4 = alloca float
%5 = alloca i32
%6 = fptrunc double %0 to float
store float %6, float* %2
store i32 1, i32* %5
br label %7
7: ; preds = %11, %1
%8 = load float, float* %2
%9 = fpext float %8 to double
%10 = fcmp ogt double %9, 0x401921FB54442D18
br i1 %10, label %11, label %16
11: ; preds = %7
%12 = load float, float* %2
%13 = fpext float %12 to double
%14 = fsub double %13, 0x401921FB54442D18
%15 = fptrunc double %14 to float
store float %15, float* %2
br label %7
16: ; preds = %7
br label %17
17: ; preds = %21, %16
%18 = load float, float* %2
%19 = fpext float %18 to double
%20 = fcmp olt double %19, 0xC01921FB54442D18
br i1 %20, label %21, label %26
21: ; preds = %17
%22 = load float, float* %2
%23 = fpext float %22 to double
%24 = fadd double %23, 0x401921FB54442D18
%25 = fptrunc double %24 to float
store float %25, float* %2
br label %17
26: ; preds = %17
%27 = load float, float* %2
store float %27, float* %4
store float %27, float* %3
%28 = load float, float* %4
%29 = load float, float* %2
%30 = load float, float* %2
%31 = fmul float %29, %30
%32 = fneg float %31
%33 = fmul float %28, %32
%34 = fpext float %33 to double
%35 = load i32, i32* %5
%36 = sitofp i32 %35 to double
%37 = fmul double 2.000000e+00, %36
%38 = load i32, i32* %5
%39 = sitofp i32 %38 to double
%40 = call double @llvm.fmuladd.f64(double 2.000000e+00, double %39, double 1.000000e+00)
%41 = fmul double %37, %40
%42 = fdiv double %34, %41
%43 = fptrunc double %42 to float
store float %43, float* %4
%44 = load float, float* %3
%45 = load float, float* %4
%46 = fadd float %44, %45
store float %46, float* %3
%47 = load i32, i32* %5
%48 = add nsw i32 %47, 1
store i32 %48, i32* %5
br label %49
49: ; preds = %54, %26
%50 = load float, float* %4
%51 = fpext float %50 to double
%52 = call double @llvm.fabs.f64(double %51)
%53 = fcmp oge double %52, 1.000000e-05
br i1 %53, label %54, label %76
54: ; preds = %49
%55 = load float, float* %4
%56 = load float, float* %2
%57 = load float, float* %2
%58 = fmul float %56, %57
%59 = fneg float %58
%60 = fmul float %55, %59
%61 = fpext float %60 to double
%62 = load i32, i32* %5
%63 = sitofp i32 %62 to double
%64 = fmul double 2.000000e+00, %63
%65 = load i32, i32* %5
%66 = sitofp i32 %65 to double
%67 = call double @llvm.fmuladd.f64(double 2.000000e+00, double %66, double 1.000000e+00)
%68 = fmul double %64, %67
%69 = fdiv double %61, %68
%70 = fptrunc double %69 to float
store float %70, float* %4
%71 = load float, float* %3
%72 = load float, float* %4
%73 = fadd float %71, %72
store float %73, float* %3
%74 = load i32, i32* %5
%75 = add nsw i32 %74, 1
store i32 %75, i32* %5
br label %49
76: ; preds = %49
%77 = load float, float* %3
ret float %77
}

LLVM IR (-O2):
define float @sin(double %0) {
%2 = fptrunc double %0 to float
%3 = fpext float %2 to double
%4 = fcmp ogt double %3, 0x401921FB54442D18
br i1 %4, label %9, label %5
5: ; preds = %9, %1
%6 = phi double [ %3, %1 ], [ %13, %9 ]
%7 = phi float [ %2, %1 ], [ %12, %9 ]
%8 = fcmp olt double %6, 0xC01921FB54442D18
br i1 %8, label %15, label %21
9: ; preds = %1, %9
%10 = phi double [ %13, %9 ], [ %3, %1 ]
%11 = fadd double %10, 0xC01921FB54442D18
%12 = fptrunc double %11 to float
%13 = fpext float %12 to double
%14 = fcmp ogt double %13, 0x401921FB54442D18
br i1 %14, label %9, label %5
15: ; preds = %5, %15
%16 = phi double [ %19, %15 ], [ %6, %5 ]
%17 = fadd double %16, 0x401921FB54442D18
%18 = fptrunc double %17 to float
%19 = fpext float %18 to double
%20 = fcmp olt double %19, 0xC01921FB54442D18
br i1 %20, label %15, label %21
21: ; preds = %15, %5
%22 = phi float [ %7, %5 ], [ %18, %15 ]
%23 = fneg float %22
%24 = fmul float %22, %23
%25 = fmul float %22, %24
%26 = fdiv float %25, 6.000000e+00
%27 = fadd float %22, %26
%28 = tail call float @llvm.fabs.f32(float %26)
%29 = fpext float %28 to double
%30 = fcmp ult double %29, 1.000000e-05
br i1 %30, label %48, label %31
31: ; preds = %21, %31
%32 = phi i32 [ %44, %31 ], [ 2, %21 ]
%33 = phi float [ %42, %31 ], [ %26, %21 ]
%34 = phi float [ %43, %31 ], [ %27, %21 ]
%35 = fmul float %24, %33
%36 = fpext float %35 to double
%37 = sitofp i32 %32 to double
%38 = fmul double %37, 2.000000e+00
%39 = tail call double @llvm.fmuladd.f64(double %37, double 2.000000e+00, double 1.000000e+00)
%40 = fmul double %38, %39
%41 = fdiv double %36, %40
%42 = fptrunc double %41 to float
%43 = fadd float %34, %42
%44 = add nuw nsw i32 %32, 1
%45 = tail call float @llvm.fabs.f32(float %42)
%46 = fpext float %45 to double
%47 = fcmp ult double %46, 1.000000e-05
br i1 %47, label %48, label %31
48: ; preds = %31, %21
%49 = phi float [ %27, %21 ], [ %43, %31 ]
ret float %49
}