float fmodf_c(float x, float y) {
  int n;
  union {
    float f;
    int i;
  } yinv;
  float a;
  yinv.f = y;
  n = 0x3F800000 - (yinv.i & 0x7F800000);
  yinv.i = yinv.i + n;
  yinv.f = 1.41176471f - 0.47058824f * yinv.f;
  yinv.i = yinv.i + n;
  a = 2.0 - yinv.f * y;
  yinv.f = yinv.f * a;
  a = 2.0 - yinv.f * y;
  yinv.f = yinv.f * a;
  a = 2.0 - yinv.f * y;
  yinv.f = yinv.f * a;
  a = 2.0 - yinv.f * y;
  yinv.f = yinv.f * a;
  n = (int)(x * yinv.f);
  x = x - ((float)n) * y;
  return x;
}


x86 assembly (clang, -O0):
<fmodf_c>:
push   %rbp
mov    %rsp,%rbp
movss  %xmm0,-0x4(%rbp)
movss  %xmm1,-0x8(%rbp)
movss  -0x8(%rbp),%xmm0
movss  %xmm0,-0x10(%rbp)
mov    -0x10(%rbp),%ecx
and    $0x7f800000,%ecx
mov    $0x3f800000,%eax
sub    %ecx,%eax
mov    %eax,-0xc(%rbp)
mov    -0x10(%rbp),%eax
add    -0xc(%rbp),%eax
mov    %eax,-0x10(%rbp)
movss  -0x10(%rbp),%xmm0
movss  0xec7(%rip),%xmm1
mulss  %xmm1,%xmm0
movss  0xebf(%rip),%xmm1
addss  %xmm1,%xmm0
movss  %xmm0,-0x10(%rbp)
mov    -0x10(%rbp),%eax
add    -0xc(%rbp),%eax
mov    %eax,-0x10(%rbp)
movss  -0x10(%rbp),%xmm0
mulss  -0x8(%rbp),%xmm0
cvtss2sd %xmm0,%xmm1
movsd  0xe8b(%rip),%xmm0
subsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,-0x14(%rbp)
movss  -0x10(%rbp),%xmm0
mulss  -0x14(%rbp),%xmm0
movss  %xmm0,-0x10(%rbp)
movss  -0x10(%rbp),%xmm0
mulss  -0x8(%rbp),%xmm0
cvtss2sd %xmm0,%xmm1
movsd  0xe59(%rip),%xmm0
subsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,-0x14(%rbp)
movss  -0x10(%rbp),%xmm0
mulss  -0x14(%rbp),%xmm0
movss  %xmm0,-0x10(%rbp)
movss  -0x10(%rbp),%xmm0
mulss  -0x8(%rbp),%xmm0
cvtss2sd %xmm0,%xmm1
movsd  0xe27(%rip),%xmm0
subsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,-0x14(%rbp)
movss  -0x10(%rbp),%xmm0
mulss  -0x14(%rbp),%xmm0
movss  %xmm0,-0x10(%rbp)
movss  -0x10(%rbp),%xmm0
mulss  -0x8(%rbp),%xmm0
cvtss2sd %xmm0,%xmm1
movsd  0xdf5(%rip),%xmm0
subsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,-0x14(%rbp)
movss  -0x10(%rbp),%xmm0
mulss  -0x14(%rbp),%xmm0
movss  %xmm0,-0x10(%rbp)
movss  -0x4(%rbp),%xmm0
mulss  -0x10(%rbp),%xmm0
cvttss2si %xmm0,%eax
mov    %eax,-0xc(%rbp)
movss  -0x4(%rbp),%xmm1
cvtsi2ssl -0xc(%rbp),%xmm0
movss  -0x8(%rbp),%xmm2
movd   %xmm0,%eax
xor    $0x80000000,%eax
movd   %eax,%xmm0
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,-0x4(%rbp)
movss  -0x4(%rbp),%xmm0
pop    %rbp
ret

x86 assembly (clang, -O2):
<fmodf_c>:
movd   %xmm1,%eax
mov    %eax,%ecx
and    $0x7f800000,%ecx
mov    $0x3f800000,%edx
sub    %ecx,%edx
add    %edx,%eax
movd   %eax,%xmm2
mulss  0xedf(%rip),%xmm2
addss  0xedb(%rip),%xmm2
movd   %xmm2,%eax
add    %edx,%eax
movd   %eax,%xmm2
movdqa %xmm2,%xmm3
mulss  %xmm1,%xmm3
movss  0xec5(%rip),%xmm4
movaps %xmm4,%xmm5
subss  %xmm3,%xmm5
mulss  %xmm2,%xmm5
movaps %xmm5,%xmm2
mulss  %xmm1,%xmm2
movaps %xmm4,%xmm3
subss  %xmm2,%xmm3
mulss  %xmm5,%xmm3
movaps %xmm3,%xmm2
mulss  %xmm1,%xmm2
movaps %xmm4,%xmm5
subss  %xmm2,%xmm5
mulss  %xmm3,%xmm5
movaps %xmm5,%xmm2
mulss  %xmm1,%xmm2
subss  %xmm2,%xmm4
mulss  %xmm5,%xmm4
mulss  %xmm0,%xmm4
cvttps2dq %xmm4,%xmm2
cvtdq2ps %xmm2,%xmm2
mulss  %xmm1,%xmm2
subss  %xmm2,%xmm0
ret

LLVM IR (-O0):
%union.anon = type { float }
define float @fmodf_c(float %0, float %1) {
%3 = alloca float
%4 = alloca float
%5 = alloca i32
%6 = alloca %union.anon
%7 = alloca float
store float %0, float* %3
store float %1, float* %4
%8 = load float, float* %4
%9 = bitcast %union.anon* %6 to float*
store float %8, float* %9
%10 = bitcast %union.anon* %6 to i32*
%11 = load i32, i32* %10
%12 = and i32 %11, 2139095040
%13 = sub nsw i32 1065353216, %12
store i32 %13, i32* %5
%14 = bitcast %union.anon* %6 to i32*
%15 = load i32, i32* %14
%16 = load i32, i32* %5
%17 = add nsw i32 %15, %16
%18 = bitcast %union.anon* %6 to i32*
store i32 %17, i32* %18
%19 = bitcast %union.anon* %6 to float*
%20 = load float, float* %19
%21 = call float @llvm.fmuladd.f32(float 0xBFDE1E1E20000000, float %20, float 0x3FF69696A0000000)
%22 = bitcast %union.anon* %6 to float*
store float %21, float* %22
%23 = bitcast %union.anon* %6 to i32*
%24 = load i32, i32* %23
%25 = load i32, i32* %5
%26 = add nsw i32 %24, %25
%27 = bitcast %union.anon* %6 to i32*
store i32 %26, i32* %27
%28 = bitcast %union.anon* %6 to float*
%29 = load float, float* %28
%30 = load float, float* %4
%31 = fmul float %29, %30
%32 = fpext float %31 to double
%33 = fsub double 2.000000e+00, %32
%34 = fptrunc double %33 to float
store float %34, float* %7
%35 = bitcast %union.anon* %6 to float*
%36 = load float, float* %35
%37 = load float, float* %7
%38 = fmul float %36, %37
%39 = bitcast %union.anon* %6 to float*
store float %38, float* %39
%40 = bitcast %union.anon* %6 to float*
%41 = load float, float* %40
%42 = load float, float* %4
%43 = fmul float %41, %42
%44 = fpext float %43 to double
%45 = fsub double 2.000000e+00, %44
%46 = fptrunc double %45 to float
store float %46, float* %7
%47 = bitcast %union.anon* %6 to float*
%48 = load float, float* %47
%49 = load float, float* %7
%50 = fmul float %48, %49
%51 = bitcast %union.anon* %6 to float*
store float %50, float* %51
%52 = bitcast %union.anon* %6 to float*
%53 = load float, float* %52
%54 = load float, float* %4
%55 = fmul float %53, %54
%56 = fpext float %55 to double
%57 = fsub double 2.000000e+00, %56
%58 = fptrunc double %57 to float
store float %58, float* %7
%59 = bitcast %union.anon* %6 to float*
%60 = load float, float* %59
%61 = load float, float* %7
%62 = fmul float %60, %61
%63 = bitcast %union.anon* %6 to float*
store float %62, float* %63
%64 = bitcast %union.anon* %6 to float*
%65 = load float, float* %64
%66 = load float, float* %4
%67 = fmul float %65, %66
%68 = fpext float %67 to double
%69 = fsub double 2.000000e+00, %68
%70 = fptrunc double %69 to float
store float %70, float* %7
%71 = bitcast %union.anon* %6 to float*
%72 = load float, float* %71
%73 = load float, float* %7
%74 = fmul float %72, %73
%75 = bitcast %union.anon* %6 to float*
store float %74, float* %75
%76 = load float, float* %3
%77 = bitcast %union.anon* %6 to float*
%78 = load float, float* %77
%79 = fmul float %76, %78
%80 = fptosi float %79 to i32
store i32 %80, i32* %5
%81 = load float, float* %3
%82 = load i32, i32* %5
%83 = sitofp i32 %82 to float
%84 = load float, float* %4
%85 = fneg float %83
%86 = call float @llvm.fmuladd.f32(float %85, float %84, float %81)
store float %86, float* %3
%87 = load float, float* %3
ret float %87
}

LLVM IR (-O2):
define float @fmodf_c(float %0, float %1) {
%3 = bitcast float %1 to i32
%4 = and i32 %3, 2139095040
%5 = sub nsw i32 1065353216, %4
%6 = add nsw i32 %5, %3
%7 = bitcast i32 %6 to float
%8 = tail call float @llvm.fmuladd.f32(float %7, float 0xBFDE1E1E20000000, float 0x3FF69696A0000000)
%9 = bitcast float %8 to i32
%10 = add nsw i32 %5, %9
%11 = bitcast i32 %10 to float
%12 = fmul float %11, %1
%13 = fsub float 2.000000e+00, %12
%14 = fmul float %13, %11
%15 = fmul float %14, %1
%16 = fsub float 2.000000e+00, %15
%17 = fmul float %14, %16
%18 = fmul float %17, %1
%19 = fsub float 2.000000e+00, %18
%20 = fmul float %17, %19
%21 = fmul float %20, %1
%22 = fsub float 2.000000e+00, %21
%23 = fmul float %20, %22
%24 = fmul float %23, %0
%25 = fptosi float %24 to i32
%26 = sitofp i32 %25 to float
%27 = fneg float %26
%28 = tail call float @llvm.fmuladd.f32(float %27, float %1, float %0)
ret float %28
}