unsigned long int schrage(unsigned long int a, unsigned long int b,
                          unsigned long int m) {
  unsigned long int q, t;
  if (a == 0UL)
    return 0UL;
  q = m / a;
  t = 2 * m - (m % a) * (b / q);
  if (t >= m)
    t -= m;
  t += a * (b % q);
  return (t >= m) ? (t - m) : t;
}


x86 assembly (clang, -O0):
<schrage>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x10(%rbp)
mov    %rsi,-0x18(%rbp)
mov    %rdx,-0x20(%rbp)
cmpq   $0x0,-0x10(%rbp)
jne    1128 <schrage+0x28>
movq   $0x0,-0x8(%rbp)
jmp    11e6 <schrage+0xe6>
mov    -0x20(%rbp),%rax
xor    %ecx,%ecx
mov    %ecx,%edx
divq   -0x10(%rbp)
mov    %rax,-0x28(%rbp)
mov    -0x20(%rbp),%rax
shl    $0x1,%rax
mov    %rax,-0x38(%rbp)
mov    -0x20(%rbp),%rax
xor    %ecx,%ecx
mov    %ecx,%edx
divq   -0x10(%rbp)
mov    %rdx,-0x40(%rbp)
mov    -0x18(%rbp),%rax
xor    %ecx,%ecx
mov    %ecx,%edx
divq   -0x28(%rbp)
mov    -0x40(%rbp),%rdx
mov    %rax,%rcx
mov    -0x38(%rbp),%rax
imul   %rcx,%rdx
sub    %rdx,%rax
mov    %rax,-0x30(%rbp)
mov    -0x30(%rbp),%rax
cmp    -0x20(%rbp),%rax
jb     1193 <schrage+0x93>
mov    -0x20(%rbp),%rcx
mov    -0x30(%rbp),%rax
sub    %rcx,%rax
mov    %rax,-0x30(%rbp)
mov    -0x10(%rbp),%rax
mov    %rax,-0x48(%rbp)
mov    -0x18(%rbp),%rax
xor    %ecx,%ecx
mov    %ecx,%edx
divq   -0x28(%rbp)
mov    -0x48(%rbp),%rax
imul   %rdx,%rax
add    -0x30(%rbp),%rax
mov    %rax,-0x30(%rbp)
mov    -0x30(%rbp),%rax
cmp    -0x20(%rbp),%rax
jb     11d6 <schrage+0xd6>
mov    -0x30(%rbp),%rax
sub    -0x20(%rbp),%rax
mov    %rax,-0x50(%rbp)
jmp    11de <schrage+0xde>
mov    -0x30(%rbp),%rax
mov    %rax,-0x50(%rbp)
mov    -0x50(%rbp),%rax
mov    %rax,-0x8(%rbp)
mov    -0x8(%rbp),%rax
pop    %rbp
ret

x86 assembly (clang, -O2):
<schrage>:
test   %rdi,%rdi
je     1124 <schrage+0x24>
mov    %rdx,%r10
mov    %rdx,%rax
or     %rdi,%rax
shr    $0x20,%rax
je     112a <schrage+0x2a>
mov    %r10,%rax
xor    %edx,%edx
div    %rdi
mov    %rdx,%r8
mov    %rax,%r9
jmp    1137 <schrage+0x37>
xor    %edx,%edx
mov    %rdx,%rax
ret
mov    %r10d,%eax
xor    %edx,%edx
div    %edi
mov    %edx,%r8d
mov    %eax,%r9d
lea    (%r10,%r10,1),%rcx
mov    %rsi,%rax
or     %r9,%rax
shr    $0x20,%rax
je     1151 <schrage+0x51>
mov    %rsi,%rax
xor    %edx,%edx
div    %r9
jmp    1158 <schrage+0x58>
mov    %esi,%eax
xor    %edx,%edx
div    %r9d
imul   %r8,%rax
sub    %rax,%rcx
xor    %eax,%eax
cmp    %r10,%rcx
mov    %r10,%rsi
cmovb  %rax,%rsi
sub    %rsi,%rcx
imul   %rdi,%rdx
add    %rcx,%rdx
cmp    %r10,%rdx
cmovb  %rax,%r10
sub    %r10,%rdx
mov    %rdx,%rax
ret

LLVM IR (-O0):
define i64 @schrage(i64 %0, i64 %1, i64 %2) {
%4 = alloca i64
%5 = alloca i64
%6 = alloca i64
%7 = alloca i64
%8 = alloca i64
%9 = alloca i64
store i64 %0, i64* %5
store i64 %1, i64* %6
store i64 %2, i64* %7
%10 = load i64, i64* %5
%11 = icmp eq i64 %10, 0
br i1 %11, label %12, label %13
12: ; preds = %3
store i64 0, i64* %4
br label %53
13: ; preds = %3
%14 = load i64, i64* %7
%15 = load i64, i64* %5
%16 = udiv i64 %14, %15
store i64 %16, i64* %8
%17 = load i64, i64* %7
%18 = mul i64 2, %17
%19 = load i64, i64* %7
%20 = load i64, i64* %5
%21 = urem i64 %19, %20
%22 = load i64, i64* %6
%23 = load i64, i64* %8
%24 = udiv i64 %22, %23
%25 = mul i64 %21, %24
%26 = sub i64 %18, %25
store i64 %26, i64* %9
%27 = load i64, i64* %9
%28 = load i64, i64* %7
%29 = icmp uge i64 %27, %28
br i1 %29, label %30, label %34
30: ; preds = %13
%31 = load i64, i64* %7
%32 = load i64, i64* %9
%33 = sub i64 %32, %31
store i64 %33, i64* %9
br label %34
34: ; preds = %30, %13
%35 = load i64, i64* %5
%36 = load i64, i64* %6
%37 = load i64, i64* %8
%38 = urem i64 %36, %37
%39 = mul i64 %35, %38
%40 = load i64, i64* %9
%41 = add i64 %40, %39
store i64 %41, i64* %9
%42 = load i64, i64* %9
%43 = load i64, i64* %7
%44 = icmp uge i64 %42, %43
br i1 %44, label %45, label %49
45: ; preds = %34
%46 = load i64, i64* %9
%47 = load i64, i64* %7
%48 = sub i64 %46, %47
br label %51
49: ; preds = %34
%50 = load i64, i64* %9
br label %51
51: ; preds = %49, %45
%52 = phi i64 [ %48, %45 ], [ %50, %49 ]
store i64 %52, i64* %4
br label %53
53: ; preds = %51, %12
%54 = load i64, i64* %4
ret i64 %54
}

LLVM IR (-O2):
define i64 @schrage(i64 %0, i64 %1, i64 %2) {
%4 = icmp eq i64 %0, 0
br i1 %4, label %21, label %5
5: ; preds = %3
%6 = udiv i64 %2, %0
%7 = shl i64 %2, 1
%8 = urem i64 %2, %0
%9 = udiv i64 %1, %6
%10 = mul i64 %9, %8
%11 = sub i64 %7, %10
%12 = icmp ult i64 %11, %2
%13 = select i1 %12, i64 0, i64 %2
%14 = sub i64 %11, %13
%15 = urem i64 %1, %6
%16 = mul i64 %15, %0
%17 = add i64 %14, %16
%18 = icmp ult i64 %17, %2
%19 = select i1 %18, i64 0, i64 %2
%20 = sub i64 %17, %19
br label %21
21: ; preds = %3, %5
%22 = phi i64 [ %20, %5 ], [ 0, %3 ]
ret i64 %22
}