double wrap(double x, const double width, const double half_width) {
  if (width > 0.0) {
    while (x >= half_width)
      x -= width;
    while (x < -half_width)
      x += width;
  }
  return x;
}


x86 assembly (clang, -O0):
<wrap>:
push   %rbp
mov    %rsp,%rbp
movsd  %xmm0,-0x8(%rbp)
movsd  %xmm1,-0x10(%rbp)
movsd  %xmm2,-0x18(%rbp)
movsd  -0x10(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jbe    119b <wrap+0x9b>
jmp    112a <wrap+0x2a>
movsd  -0x8(%rbp),%xmm0
ucomisd -0x18(%rbp),%xmm0
jb     1152 <wrap+0x52>
movsd  -0x10(%rbp),%xmm1
movsd  -0x8(%rbp),%xmm0
subsd  %xmm1,%xmm0
movsd  %xmm0,-0x8(%rbp)
jmp    112a <wrap+0x2a>
jmp    1157 <wrap+0x57>
movsd  -0x8(%rbp),%xmm1
movsd  -0x18(%rbp),%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
ucomisd %xmm1,%xmm0
jbe    1196 <wrap+0x96>
movsd  -0x10(%rbp),%xmm0
addsd  -0x8(%rbp),%xmm0
movsd  %xmm0,-0x8(%rbp)
jmp    1157 <wrap+0x57>
jmp    119b <wrap+0x9b>
movsd  -0x8(%rbp),%xmm0
pop    %rbp
ret

x86 assembly (clang, -O2):
<wrap>:
xorpd  %xmm3,%xmm3
ucomisd %xmm3,%xmm1
jbe    113a <wrap+0x3a>
nopw   0x0(%rax,%rax,1)
ucomisd %xmm2,%xmm0
jb     1120 <wrap+0x20>
subsd  %xmm1,%xmm0
ucomisd %xmm2,%xmm0
jae    1116 <wrap+0x16>
xorpd  0xed8(%rip),%xmm2
ucomisd %xmm0,%xmm2
jbe    113a <wrap+0x3a>
xchg   %ax,%ax
addsd  %xmm1,%xmm0
ucomisd %xmm0,%xmm2
ja     1130 <wrap+0x30>
ret

LLVM IR (-O0):
define double @wrap(double %0, double %1, double %2) {
%4 = alloca double
%5 = alloca double
%6 = alloca double
store double %0, double* %4
store double %1, double* %5
store double %2, double* %6
%7 = load double, double* %5
%8 = fcmp ogt double %7, 0.000000e+00
br i1 %8, label %9, label %29
9: ; preds = %3
br label %10
10: ; preds = %14, %9
%11 = load double, double* %4
%12 = load double, double* %6
%13 = fcmp oge double %11, %12
br i1 %13, label %14, label %18
14: ; preds = %10
%15 = load double, double* %5
%16 = load double, double* %4
%17 = fsub double %16, %15
store double %17, double* %4
br label %10
18: ; preds = %10
br label %19
19: ; preds = %24, %18
%20 = load double, double* %4
%21 = load double, double* %6
%22 = fneg double %21
%23 = fcmp olt double %20, %22
br i1 %23, label %24, label %28
24: ; preds = %19
%25 = load double, double* %5
%26 = load double, double* %4
%27 = fadd double %26, %25
store double %27, double* %4
br label %19
28: ; preds = %19
br label %29
29: ; preds = %28, %3
%30 = load double, double* %4
ret double %30
}

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