void fvec_splat_add(const float *a, int n, const int *assign, float *accu) {
  int i;
  for (i = 0; i < n; i++) {
    accu[assign[i]] += a[i];
  }
}


x86 assembly (clang, -O0):
<fvec_splat_add>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %rcx,-0x20(%rbp)
movl   $0x0,-0x24(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    115b <fvec_splat_add+0x5b>
mov    -0x8(%rbp),%rax
movslq -0x24(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x20(%rbp),%rax
mov    -0x18(%rbp),%rcx
movslq -0x24(%rbp),%rdx
movslq (%rcx,%rdx,4),%rcx
addss  (%rax,%rcx,4),%xmm0
movss  %xmm0,(%rax,%rcx,4)
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    111a <fvec_splat_add+0x1a>
pop    %rbp
ret

x86 assembly (clang, -O2):
<fvec_splat_add>:
test   %esi,%esi
jle    116a <fvec_splat_add+0x6a>
mov    %esi,%r8d
cmp    $0x1,%esi
jne    1110 <fvec_splat_add+0x10>
xor    %esi,%esi
jmp    1151 <fvec_splat_add+0x51>
mov    %r8d,%r9d
and    $0xfffffffe,%r9d
xor    %esi,%esi
nopl   0x0(%rax)
movss  (%rdi,%rsi,4),%xmm0
movslq (%rdx,%rsi,4),%rax
addss  (%rcx,%rax,4),%xmm0
movss  %xmm0,(%rcx,%rax,4)
movss  0x4(%rdi,%rsi,4),%xmm0
movslq 0x4(%rdx,%rsi,4),%rax
addss  (%rcx,%rax,4),%xmm0
movss  %xmm0,(%rcx,%rax,4)
add    $0x2,%rsi
cmp    %rsi,%r9
jne    1120 <fvec_splat_add+0x20>
test   $0x1,%r8b
je     116a <fvec_splat_add+0x6a>
movss  (%rdi,%rsi,4),%xmm0
movslq (%rdx,%rsi,4),%rax
addss  (%rcx,%rax,4),%xmm0
movss  %xmm0,(%rcx,%rax,4)
ret

LLVM IR (-O0):
define void @fvec_splat_add(float* %0, i32 %1, i32* %2, float* %3) {
%5 = alloca float*
%6 = alloca i32
%7 = alloca i32*
%8 = alloca float*
%9 = alloca i32
store float* %0, float** %5
store i32 %1, i32* %6
store i32* %2, i32** %7
store float* %3, float** %8
store i32 0, i32* %9
br label %10
10: ; preds = %30, %4
%11 = load i32, i32* %9
%12 = load i32, i32* %6
%13 = icmp slt i32 %11, %12
br i1 %13, label %14, label %33
14: ; preds = %10
%15 = load float*, float** %5
%16 = load i32, i32* %9
%17 = sext i32 %16 to i64
%18 = getelementptr inbounds float, float* %15, i64 %17
%19 = load float, float* %18
%20 = load float*, float** %8
%21 = load i32*, i32** %7
%22 = load i32, i32* %9
%23 = sext i32 %22 to i64
%24 = getelementptr inbounds i32, i32* %21, i64 %23
%25 = load i32, i32* %24
%26 = sext i32 %25 to i64
%27 = getelementptr inbounds float, float* %20, i64 %26
%28 = load float, float* %27
%29 = fadd float %28, %19
store float %29, float* %27
br label %30
30: ; preds = %14
%31 = load i32, i32* %9
%32 = add nsw i32 %31, 1
store i32 %32, i32* %9
br label %10
33: ; preds = %10
ret void
}

LLVM IR (-O2):
define void @fvec_splat_add(float* nocapture %0, i32 %1, i32* nocapture %2, float* nocapture %3) {
%5 = icmp sgt i32 %1, 0
br i1 %5, label %6, label %47
6: ; preds = %4
%7 = zext i32 %1 to i64
%8 = and i64 %7, 1
%9 = icmp eq i32 %1, 1
br i1 %9, label %35, label %10
10: ; preds = %6
%11 = and i64 %7, 4294967294
br label %12
12: ; preds = %12, %10
%13 = phi i64 [ 0, %10 ], [ %32, %12 ]
%14 = phi i64 [ 0, %10 ], [ %33, %12 ]
%15 = getelementptr inbounds float, float* %0, i64 %13
%16 = load float, float* %15
%17 = getelementptr inbounds i32, i32* %2, i64 %13
%18 = load i32, i32* %17
%19 = sext i32 %18 to i64
%20 = getelementptr inbounds float, float* %3, i64 %19
%21 = load float, float* %20
%22 = fadd float %16, %21
store float %22, float* %20
%23 = or i64 %13, 1
%24 = getelementptr inbounds float, float* %0, i64 %23
%25 = load float, float* %24
%26 = getelementptr inbounds i32, i32* %2, i64 %23
%27 = load i32, i32* %26
%28 = sext i32 %27 to i64
%29 = getelementptr inbounds float, float* %3, i64 %28
%30 = load float, float* %29
%31 = fadd float %25, %30
store float %31, float* %29
%32 = add nuw nsw i64 %13, 2
%33 = add i64 %14, 2
%34 = icmp eq i64 %33, %11
br i1 %34, label %35, label %12
35: ; preds = %12, %6
%36 = phi i64 [ 0, %6 ], [ %32, %12 ]
%37 = icmp eq i64 %8, 0
br i1 %37, label %47, label %38
38: ; preds = %35
%39 = getelementptr inbounds float, float* %0, i64 %36
%40 = load float, float* %39
%41 = getelementptr inbounds i32, i32* %2, i64 %36
%42 = load i32, i32* %41
%43 = sext i32 %42 to i64
%44 = getelementptr inbounds float, float* %3, i64 %43
%45 = load float, float* %44
%46 = fadd float %40, %45
store float %46, float* %44
br label %47
47: ; preds = %38, %35, %4
ret void
}