float lms_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 (lms_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):
<lms_sin>:
push   %rbp
mov    %rsp,%rbp
sub    $0x10,%rsp
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,-0x4(%rbp)
movl   $0x1,-0x10(%rbp)
movss  -0x4(%rbp),%xmm0
cvtss2sd %xmm0,%xmm0
movsd  0xec7(%rip),%xmm1
ucomisd %xmm1,%xmm0
jbe    1166 <lms_sin+0x56>
movss  -0x4(%rbp),%xmm0
cvtss2sd %xmm0,%xmm0
movsd  0xeac(%rip),%xmm1
subsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,-0x4(%rbp)
jmp    1128 <lms_sin+0x18>
jmp    116b <lms_sin+0x5b>
movss  -0x4(%rbp),%xmm0
cvtss2sd %xmm0,%xmm1
movsd  0xe8c(%rip),%xmm0
ucomisd %xmm1,%xmm0
jbe    11a9 <lms_sin+0x99>
movss  -0x4(%rbp),%xmm0
cvtss2sd %xmm0,%xmm0
movsd  0xe69(%rip),%xmm1
addsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,-0x4(%rbp)
jmp    116b <lms_sin+0x5b>
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  0xe2f(%rip),%xmm1
mulsd  %xmm2,%xmm1
cvtsi2sdl -0x10(%rbp),%xmm2
addsd  %xmm2,%xmm2
movsd  0xe12(%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
mov    $0x1,%al
call   1030 <lms_fabs@plt>
cvtsi2sd %eax,%xmm0
movsd  0xdd9(%rip),%xmm1
ucomisd %xmm1,%xmm0
jb     12c9 <lms_sin+0x1b9>
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  0xd96(%rip),%xmm1
mulsd  %xmm2,%xmm1
cvtsi2sdl -0x10(%rbp),%xmm2
addsd  %xmm2,%xmm2
movsd  0xd79(%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    122b <lms_sin+0x11b>
movss  -0x8(%rbp),%xmm0
add    $0x10,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<lms_sin>:
push   %rbx
sub    $0x30,%rsp
cvtsd2ss %xmm0,%xmm3
xorps  %xmm0,%xmm0
cvtss2sd %xmm3,%xmm0
ucomisd 0xed8(%rip),%xmm0
jbe    1158 <lms_sin+0x48>
movsd  0xed6(%rip),%xmm1
movsd  0xec6(%rip),%xmm2
nopw   0x0(%rax,%rax,1)
addsd  %xmm1,%xmm0
xorps  %xmm3,%xmm3
cvtsd2ss %xmm0,%xmm3
xorps  %xmm0,%xmm0
cvtss2sd %xmm3,%xmm0
ucomisd %xmm2,%xmm0
ja     1140 <lms_sin+0x30>
movsd  0xea8(%rip),%xmm1
ucomisd %xmm0,%xmm1
jbe    1188 <lms_sin+0x78>
movsd  0xe92(%rip),%xmm2
xchg   %ax,%ax
addsd  %xmm2,%xmm0
xorps  %xmm3,%xmm3
cvtsd2ss %xmm0,%xmm3
xorps  %xmm0,%xmm0
cvtss2sd %xmm3,%xmm0
ucomisd %xmm0,%xmm1
ja     1170 <lms_sin+0x60>
movaps 0xe91(%rip),%xmm0
xorps  %xmm3,%xmm0
mulss  %xmm3,%xmm0
movaps %xmm3,%xmm1
movaps %xmm0,0x20(%rsp)
mulss  %xmm0,%xmm1
divss  0xe86(%rip),%xmm1
addss  %xmm1,%xmm3
movaps %xmm3,(%rsp)
movaps %xmm1,0x10(%rsp)
xorps  %xmm0,%xmm0
cvtss2sd %xmm1,%xmm0
mov    $0x1,%al
call   1030 <lms_fabs@plt>
test   %eax,%eax
jle    1240 <lms_sin+0x130>
mov    $0x2,%ebx
movaps (%rsp),%xmm0
movaps 0x10(%rsp),%xmm1
nopw   0x0(%rax,%rax,1)
mulss  0x20(%rsp),%xmm1
xorps  %xmm3,%xmm3
cvtss2sd %xmm1,%xmm3
xorps  %xmm1,%xmm1
cvtsi2sd %ebx,%xmm1
addsd  %xmm1,%xmm1
movapd %xmm1,%xmm2
addsd  0xe0c(%rip),%xmm2
mulsd  %xmm1,%xmm2
divsd  %xmm2,%xmm3
xorps  %xmm1,%xmm1
cvtsd2ss %xmm3,%xmm1
movaps %xmm1,0x10(%rsp)
addss  %xmm1,%xmm0
movaps %xmm0,(%rsp)
add    $0x1,%ebx
xorps  %xmm0,%xmm0
cvtss2sd %xmm1,%xmm0
mov    $0x1,%al
call   1030 <lms_fabs@plt>
movaps 0x10(%rsp),%xmm1
movaps (%rsp),%xmm0
test   %eax,%eax
jg     11e0 <lms_sin+0xd0>
jmp    1244 <lms_sin+0x134>
movaps (%rsp),%xmm0
add    $0x30,%rsp
pop    %rbx
ret

LLVM IR (-O0):
define float @lms_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 = %55, %26
%50 = load float, float* %4
%51 = fpext float %50 to double
%52 = call i32 (double, ...) bitcast (i32 (...)* @lms_fabs to i32 (double, ...)*)(double %51)
%53 = sitofp i32 %52 to double
%54 = fcmp oge double %53, 1.000000e-05
br i1 %54, label %55, label %77
55: ; preds = %49
%56 = load float, float* %4
%57 = load float, float* %2
%58 = load float, float* %2
%59 = fmul float %57, %58
%60 = fneg float %59
%61 = fmul float %56, %60
%62 = fpext float %61 to double
%63 = load i32, i32* %5
%64 = sitofp i32 %63 to double
%65 = fmul double 2.000000e+00, %64
%66 = load i32, i32* %5
%67 = sitofp i32 %66 to double
%68 = call double @llvm.fmuladd.f64(double 2.000000e+00, double %67, double 1.000000e+00)
%69 = fmul double %65, %68
%70 = fdiv double %62, %69
%71 = fptrunc double %70 to float
store float %71, float* %4
%72 = load float, float* %3
%73 = load float, float* %4
%74 = fadd float %72, %73
store float %74, float* %3
%75 = load i32, i32* %5
%76 = add nsw i32 %75, 1
store i32 %76, i32* %5
br label %49
77: ; preds = %49
%78 = load float, float* %3
ret float %78
}

LLVM IR (-O2):
define float @lms_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 = fpext float %26 to double
%29 = tail call i32 (double, ...) bitcast (i32 (...)* @lms_fabs to i32 (double, ...)*)(double %28)
%30 = icmp sgt i32 %29, 0
br i1 %30, label %31, label %48
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 = fpext float %42 to double
%46 = tail call i32 (double, ...) bitcast (i32 (...)* @lms_fabs to i32 (double, ...)*)(double %45)
%47 = icmp sgt i32 %46, 0
br i1 %47, label %31, label %48
48: ; preds = %31, %21
%49 = phi float [ %27, %21 ], [ %43, %31 ]
ret float %49
}