void arcsum(unsigned long iin[], unsigned long iout[], unsigned long ja,
            int nwk, unsigned long nrad, unsigned long nc) {
  int j, karry = 0;
  unsigned long jtmp;
  for (j = nwk; j > nc; j--) {
    jtmp = ja;
    ja /= nrad;
    iout[j] = iin[j] + (jtmp - ja * nrad) + karry;
    if (iout[j] >= nrad) {
      iout[j] -= nrad;
      karry = 1;
    } else
      karry = 0;
  }
  iout[nc] = iin[nc] + ja + karry;
}


x86 assembly (clang, -O0):
<arcsum>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %ecx,-0x1c(%rbp)
mov    %r8,-0x28(%rbp)
mov    %r9,-0x30(%rbp)
movl   $0x0,-0x38(%rbp)
mov    -0x1c(%rbp),%eax
mov    %eax,-0x34(%rbp)
movslq -0x34(%rbp),%rax
cmp    -0x30(%rbp),%rax
jbe    11d4 <arcsum+0xd4>
mov    -0x18(%rbp),%rax
mov    %rax,-0x40(%rbp)
mov    -0x28(%rbp),%rcx
mov    -0x18(%rbp),%rax
xor    %edx,%edx
div    %rcx
mov    %rax,-0x18(%rbp)
mov    -0x8(%rbp),%rax
movslq -0x34(%rbp),%rcx
mov    (%rax,%rcx,8),%rdx
mov    -0x40(%rbp),%rax
mov    -0x18(%rbp),%rcx
imul   -0x28(%rbp),%rcx
sub    %rcx,%rax
add    %rax,%rdx
movslq -0x38(%rbp),%rax
add    %rax,%rdx
mov    -0x10(%rbp),%rax
movslq -0x34(%rbp),%rcx
mov    %rdx,(%rax,%rcx,8)
mov    -0x10(%rbp),%rax
movslq -0x34(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
cmp    -0x28(%rbp),%rax
jb     11ba <arcsum+0xba>
mov    -0x28(%rbp),%rsi
mov    -0x10(%rbp),%rax
movslq -0x34(%rbp),%rcx
mov    (%rax,%rcx,8),%rdx
sub    %rsi,%rdx
mov    %rdx,(%rax,%rcx,8)
movl   $0x1,-0x38(%rbp)
jmp    11c1 <arcsum+0xc1>
movl   $0x0,-0x38(%rbp)
jmp    11c6 <arcsum+0xc6>
mov    -0x34(%rbp),%eax
add    $0xffffffff,%eax
mov    %eax,-0x34(%rbp)
jmp    1128 <arcsum+0x28>
mov    -0x8(%rbp),%rax
mov    -0x30(%rbp),%rcx
mov    (%rax,%rcx,8),%rdx
add    -0x18(%rbp),%rdx
movslq -0x38(%rbp),%rax
add    %rax,%rdx
mov    -0x10(%rbp),%rax
mov    -0x30(%rbp),%rcx
mov    %rdx,(%rax,%rcx,8)
pop    %rbp
ret

x86 assembly (clang, -O2):
<arcsum>:
mov    %rdx,%r10
movslq %ecx,%r11
xor    %ecx,%ecx
cmp    %r9,%r11
ja     115d <arcsum+0x5d>
mov    %r10,%rax
add    %rax,%rcx
add    (%rdi,%r9,8),%rcx
mov    %rcx,(%rsi,%r9,8)
ret
nopl   0x0(%rax)
mov    %r10,%rax
xor    %edx,%edx
div    %r8
mov    %rax,%rdx
add    %rcx,%r10
add    (%rdi,%r11,8),%r10
imul   %r8,%rdx
sub    %rdx,%r10
xor    %ecx,%ecx
cmp    %r8,%r10
setae  %cl
mov    $0x0,%edx
cmovae %r8,%rdx
sub    %rdx,%r10
mov    %r10,(%rsi,%r11,8)
add    $0xffffffffffffffff,%r11
mov    %rax,%r10
cmp    %r9,%r11
jbe    1110 <arcsum+0x10>
mov    %r10,%rax
or     %r8,%rax
shr    $0x20,%rax
jne    1120 <arcsum+0x20>
mov    %r10d,%eax
xor    %edx,%edx
div    %r8d
jmp    1128 <arcsum+0x28>

LLVM IR (-O0):
define void @arcsum(i64* %0, i64* %1, i64 %2, i32 %3, i64 %4, i64 %5) {
%7 = alloca i64*
%8 = alloca i64*
%9 = alloca i64
%10 = alloca i32
%11 = alloca i64
%12 = alloca i64
%13 = alloca i32
%14 = alloca i32
%15 = alloca i64
store i64* %0, i64** %7
store i64* %1, i64** %8
store i64 %2, i64* %9
store i32 %3, i32* %10
store i64 %4, i64* %11
store i64 %5, i64* %12
store i32 0, i32* %14
%16 = load i32, i32* %10
store i32 %16, i32* %13
br label %17
17: ; preds = %62, %6
%18 = load i32, i32* %13
%19 = sext i32 %18 to i64
%20 = load i64, i64* %12
%21 = icmp ugt i64 %19, %20
br i1 %21, label %22, label %65
22: ; preds = %17
%23 = load i64, i64* %9
store i64 %23, i64* %15
%24 = load i64, i64* %11
%25 = load i64, i64* %9
%26 = udiv i64 %25, %24
store i64 %26, i64* %9
%27 = load i64*, i64** %7
%28 = load i32, i32* %13
%29 = sext i32 %28 to i64
%30 = getelementptr inbounds i64, i64* %27, i64 %29
%31 = load i64, i64* %30
%32 = load i64, i64* %15
%33 = load i64, i64* %9
%34 = load i64, i64* %11
%35 = mul i64 %33, %34
%36 = sub i64 %32, %35
%37 = add i64 %31, %36
%38 = load i32, i32* %14
%39 = sext i32 %38 to i64
%40 = add i64 %37, %39
%41 = load i64*, i64** %8
%42 = load i32, i32* %13
%43 = sext i32 %42 to i64
%44 = getelementptr inbounds i64, i64* %41, i64 %43
store i64 %40, i64* %44
%45 = load i64*, i64** %8
%46 = load i32, i32* %13
%47 = sext i32 %46 to i64
%48 = getelementptr inbounds i64, i64* %45, i64 %47
%49 = load i64, i64* %48
%50 = load i64, i64* %11
%51 = icmp uge i64 %49, %50
br i1 %51, label %52, label %60
52: ; preds = %22
%53 = load i64, i64* %11
%54 = load i64*, i64** %8
%55 = load i32, i32* %13
%56 = sext i32 %55 to i64
%57 = getelementptr inbounds i64, i64* %54, i64 %56
%58 = load i64, i64* %57
%59 = sub i64 %58, %53
store i64 %59, i64* %57
store i32 1, i32* %14
br label %61
60: ; preds = %22
store i32 0, i32* %14
br label %61
61: ; preds = %60, %52
br label %62
62: ; preds = %61
%63 = load i32, i32* %13
%64 = add nsw i32 %63, -1
store i32 %64, i32* %13
br label %17
65: ; preds = %17
%66 = load i64*, i64** %7
%67 = load i64, i64* %12
%68 = getelementptr inbounds i64, i64* %66, i64 %67
%69 = load i64, i64* %68
%70 = load i64, i64* %9
%71 = add i64 %69, %70
%72 = load i32, i32* %14
%73 = sext i32 %72 to i64
%74 = add i64 %71, %73
%75 = load i64*, i64** %8
%76 = load i64, i64* %12
%77 = getelementptr inbounds i64, i64* %75, i64 %76
store i64 %74, i64* %77
ret void
}

LLVM IR (-O2):
define void @arcsum(i64* nocapture %0, i64* nocapture writeonly %1, i64 %2, i32 %3, i64 %4, i64 %5) {
%7 = sext i32 %3 to i64
%8 = icmp ugt i64 %7, %5
br i1 %8, label %9, label %30
9: ; preds = %6, %9
%10 = phi i64 [ %26, %9 ], [ %7, %6 ]
%11 = phi i64 [ %13, %9 ], [ %2, %6 ]
%12 = phi i64 [ %25, %9 ], [ 0, %6 ]
%13 = udiv i64 %11, %4
%14 = getelementptr inbounds i64, i64* %0, i64 %10
%15 = load i64, i64* %14
%16 = mul i64 %13, %4
%17 = add i64 %11, %12
%18 = add i64 %17, %15
%19 = sub i64 %18, %16
%20 = getelementptr inbounds i64, i64* %1, i64 %10
%21 = icmp ult i64 %19, %4
%22 = select i1 %21, i64 0, i64 %4
%23 = sub i64 %19, %22
%24 = xor i1 %21, true
%25 = zext i1 %24 to i64
store i64 %23, i64* %20
%26 = add i64 %10, -1
%27 = icmp ugt i64 %26, %5
br i1 %27, label %9, label %28
28: ; preds = %9
%29 = zext i1 %24 to i64
br label %30
30: ; preds = %28, %6
%31 = phi i64 [ 0, %6 ], [ %29, %28 ]
%32 = phi i64 [ %2, %6 ], [ %13, %28 ]
%33 = getelementptr inbounds i64, i64* %0, i64 %5
%34 = load i64, i64* %33
%35 = add i64 %32, %31
%36 = add i64 %35, %34
%37 = getelementptr inbounds i64, i64* %1, i64 %5
store i64 %36, i64* %37
ret void
}