double modulo(double a, double b) {
  double sign = a / (a > 0 ? a : -a);
  a = (a > 0 ? a : -a);
  b = (b > 0 ? b : -b);
  while (a >= b)
    a -= b;
  return a * sign;
}


x86 assembly (clang, -O0):
<modulo>:
push   %rbp
mov    %rsp,%rbp
movsd  %xmm0,-0x8(%rbp)
movsd  %xmm1,-0x10(%rbp)
movsd  -0x8(%rbp),%xmm0
movsd  %xmm0,-0x20(%rbp)
movsd  -0x8(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jbe    1139 <modulo+0x39>
movsd  -0x8(%rbp),%xmm0
movsd  %xmm0,-0x28(%rbp)
jmp    115a <modulo+0x5a>
movsd  -0x8(%rbp),%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
movsd  %xmm0,-0x28(%rbp)
movsd  -0x20(%rbp),%xmm0
movsd  -0x28(%rbp),%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x18(%rbp)
movsd  -0x8(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jbe    118e <modulo+0x8e>
movsd  -0x8(%rbp),%xmm0
movsd  %xmm0,-0x30(%rbp)
jmp    11af <modulo+0xaf>
movsd  -0x8(%rbp),%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
movsd  %xmm0,-0x30(%rbp)
movsd  -0x30(%rbp),%xmm0
movsd  %xmm0,-0x8(%rbp)
movsd  -0x10(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jbe    11da <modulo+0xda>
movsd  -0x10(%rbp),%xmm0
movsd  %xmm0,-0x38(%rbp)
jmp    11fb <modulo+0xfb>
movsd  -0x10(%rbp),%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
movsd  %xmm0,-0x38(%rbp)
movsd  -0x38(%rbp),%xmm0
movsd  %xmm0,-0x10(%rbp)
movsd  -0x8(%rbp),%xmm0
ucomisd -0x10(%rbp),%xmm0
jb     122d <modulo+0x12d>
movsd  -0x10(%rbp),%xmm1
movsd  -0x8(%rbp),%xmm0
subsd  %xmm1,%xmm0
movsd  %xmm0,-0x8(%rbp)
jmp    1205 <modulo+0x105>
movsd  -0x8(%rbp),%xmm0
mulsd  -0x18(%rbp),%xmm0
pop    %rbp
ret

x86 assembly (clang, -O2):
<modulo>:
movapd 0xef8(%rip),%xmm4
movapd %xmm0,%xmm5
xorpd  %xmm4,%xmm5
xorpd  %xmm3,%xmm3
xorpd  %xmm2,%xmm2
cmpltsd %xmm0,%xmm2
movapd %xmm2,%xmm6
andnpd %xmm5,%xmm6
andpd  %xmm0,%xmm2
orpd   %xmm6,%xmm2
xorpd  %xmm1,%xmm4
cmpltsd %xmm1,%xmm3
movapd %xmm3,%xmm5
andnpd %xmm4,%xmm5
andpd  %xmm1,%xmm3
orpd   %xmm5,%xmm3
movapd %xmm2,%xmm1
ucomisd %xmm3,%xmm2
jb     115a <modulo+0x5a>
subsd  %xmm3,%xmm1
ucomisd %xmm3,%xmm1
jae    1150 <modulo+0x50>
divsd  %xmm2,%xmm0
mulsd  %xmm1,%xmm0
ret

LLVM IR (-O0):
define double @modulo(double %0, double %1) {
%3 = alloca double
%4 = alloca double
%5 = alloca double
store double %0, double* %3
store double %1, double* %4
%6 = load double, double* %3
%7 = load double, double* %3
%8 = fcmp ogt double %7, 0.000000e+00
br i1 %8, label %9, label %11
9: ; preds = %2
%10 = load double, double* %3
br label %14
11: ; preds = %2
%12 = load double, double* %3
%13 = fneg double %12
br label %14
14: ; preds = %11, %9
%15 = phi double [ %10, %9 ], [ %13, %11 ]
%16 = fdiv double %6, %15
store double %16, double* %5
%17 = load double, double* %3
%18 = fcmp ogt double %17, 0.000000e+00
br i1 %18, label %19, label %21
19: ; preds = %14
%20 = load double, double* %3
br label %24
21: ; preds = %14
%22 = load double, double* %3
%23 = fneg double %22
br label %24
24: ; preds = %21, %19
%25 = phi double [ %20, %19 ], [ %23, %21 ]
store double %25, double* %3
%26 = load double, double* %4
%27 = fcmp ogt double %26, 0.000000e+00
br i1 %27, label %28, label %30
28: ; preds = %24
%29 = load double, double* %4
br label %33
30: ; preds = %24
%31 = load double, double* %4
%32 = fneg double %31
br label %33
33: ; preds = %30, %28
%34 = phi double [ %29, %28 ], [ %32, %30 ]
store double %34, double* %4
br label %35
35: ; preds = %39, %33
%36 = load double, double* %3
%37 = load double, double* %4
%38 = fcmp oge double %36, %37
br i1 %38, label %39, label %43
39: ; preds = %35
%40 = load double, double* %4
%41 = load double, double* %3
%42 = fsub double %41, %40
store double %42, double* %3
br label %35
43: ; preds = %35
%44 = load double, double* %3
%45 = load double, double* %5
%46 = fmul double %44, %45
ret double %46
}

LLVM IR (-O2):
define double @modulo(double %0, double %1) {
%3 = fcmp ogt double %0, 0.000000e+00
%4 = fneg double %0
%5 = select i1 %3, double %0, double %4
%6 = fcmp ogt double %1, 0.000000e+00
%7 = fneg double %1
%8 = select i1 %6, double %1, double %7
%9 = fcmp ult double %5, %8
br i1 %9, label %14, label %10
10: ; preds = %2, %10
%11 = phi double [ %12, %10 ], [ %5, %2 ]
%12 = fsub double %11, %8
%13 = fcmp ult double %12, %8
br i1 %13, label %14, label %10
14: ; preds = %10, %2
%15 = phi double [ %5, %2 ], [ %12, %10 ]
%16 = fdiv double %0, %5
%17 = fmul double %16, %15
ret double %17
}