float modulus(float x, float y) {
  float z = x;
  int n;
  if (z < 0) {
    n = (int)((-z) / y) + 1;
    z += n * y;
  }
  n = (int)(z / y);
  z -= n * y;
  return z;
}


x86 assembly (clang, -O0):
<modulus>:
push   %rbp
mov    %rsp,%rbp
movss  %xmm0,-0x4(%rbp)
movss  %xmm1,-0x8(%rbp)
movss  -0x4(%rbp),%xmm0
movss  %xmm0,-0xc(%rbp)
xorps  %xmm0,%xmm0
ucomiss -0xc(%rbp),%xmm0
jbe    1162 <modulus+0x62>
movss  -0xc(%rbp),%xmm0
movd   %xmm0,%eax
xor    $0x80000000,%eax
movd   %eax,%xmm0
divss  -0x8(%rbp),%xmm0
cvttss2si %xmm0,%eax
add    $0x1,%eax
mov    %eax,-0x10(%rbp)
cvtsi2ssl -0x10(%rbp),%xmm0
movss  -0x8(%rbp),%xmm2
movss  -0xc(%rbp),%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,-0xc(%rbp)
movss  -0xc(%rbp),%xmm0
divss  -0x8(%rbp),%xmm0
cvttss2si %xmm0,%eax
mov    %eax,-0x10(%rbp)
cvtsi2ssl -0x10(%rbp),%xmm0
movss  -0x8(%rbp),%xmm2
movss  -0xc(%rbp),%xmm1
movd   %xmm0,%eax
xor    $0x80000000,%eax
movd   %eax,%xmm0
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,-0xc(%rbp)
movss  -0xc(%rbp),%xmm0
pop    %rbp
ret

x86 assembly (clang, -O2):
<modulus>:
xorps  %xmm2,%xmm2
ucomiss %xmm0,%xmm2
jbe    112c <modulus+0x2c>
movaps 0xef1(%rip),%xmm2
xorps  %xmm0,%xmm2
divss  %xmm1,%xmm2
cvttss2si %xmm2,%eax
add    $0x1,%eax
xorps  %xmm2,%xmm2
cvtsi2ss %eax,%xmm2
mulss  %xmm1,%xmm2
addss  %xmm2,%xmm0
movaps %xmm0,%xmm2
divss  %xmm1,%xmm2
cvttps2dq %xmm2,%xmm2
cvtdq2ps %xmm2,%xmm2
mulss  %xmm1,%xmm2
subss  %xmm2,%xmm0
ret

LLVM IR (-O0):
define float @modulus(float %0, float %1) {
%3 = alloca float
%4 = alloca float
%5 = alloca float
%6 = alloca i32
store float %0, float* %3
store float %1, float* %4
%7 = load float, float* %3
store float %7, float* %5
%8 = load float, float* %5
%9 = fcmp olt float %8, 0.000000e+00
br i1 %9, label %10, label %22
10: ; preds = %2
%11 = load float, float* %5
%12 = fneg float %11
%13 = load float, float* %4
%14 = fdiv float %12, %13
%15 = fptosi float %14 to i32
%16 = add nsw i32 %15, 1
store i32 %16, i32* %6
%17 = load i32, i32* %6
%18 = sitofp i32 %17 to float
%19 = load float, float* %4
%20 = load float, float* %5
%21 = call float @llvm.fmuladd.f32(float %18, float %19, float %20)
store float %21, float* %5
br label %22
22: ; preds = %10, %2
%23 = load float, float* %5
%24 = load float, float* %4
%25 = fdiv float %23, %24
%26 = fptosi float %25 to i32
store i32 %26, i32* %6
%27 = load i32, i32* %6
%28 = sitofp i32 %27 to float
%29 = load float, float* %4
%30 = load float, float* %5
%31 = fneg float %28
%32 = call float @llvm.fmuladd.f32(float %31, float %29, float %30)
store float %32, float* %5
%33 = load float, float* %5
ret float %33
}

LLVM IR (-O2):
define float @modulus(float %0, float %1) {
%3 = fcmp olt float %0, 0.000000e+00
br i1 %3, label %4, label %11
4: ; preds = %2
%5 = fneg float %0
%6 = fdiv float %5, %1
%7 = fptosi float %6 to i32
%8 = add nsw i32 %7, 1
%9 = sitofp i32 %8 to float
%10 = tail call float @llvm.fmuladd.f32(float %9, float %1, float %0)
br label %11
11: ; preds = %4, %2
%12 = phi float [ %10, %4 ], [ %0, %2 ]
%13 = fdiv float %12, %1
%14 = fptosi float %13 to i32
%15 = sitofp i32 %14 to float
%16 = fneg float %15
%17 = tail call float @llvm.fmuladd.f32(float %16, float %1, float %12)
ret float %17
}