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


x86 assembly (clang, -O0):
<fvec_isplat_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_isplat_add+0x5b>
mov    -0x8(%rbp),%rax
mov    -0x18(%rbp),%rcx
movslq -0x24(%rbp),%rdx
movslq (%rcx,%rdx,4),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x20(%rbp),%rax
movslq -0x24(%rbp),%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_isplat_add+0x1a>
pop    %rbp
ret

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

LLVM IR (-O0):
define void @fvec_isplat_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** %7
%17 = load i32, i32* %9
%18 = sext i32 %17 to i64
%19 = getelementptr inbounds i32, i32* %16, i64 %18
%20 = load i32, i32* %19
%21 = sext i32 %20 to i64
%22 = getelementptr inbounds float, float* %15, i64 %21
%23 = load float, float* %22
%24 = load float*, float** %8
%25 = load i32, i32* %9
%26 = sext i32 %25 to i64
%27 = getelementptr inbounds float, float* %24, i64 %26
%28 = load float, float* %27
%29 = fadd float %28, %23
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_isplat_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 i32, i32* %2, i64 %13
%16 = load i32, i32* %15
%17 = sext i32 %16 to i64
%18 = getelementptr inbounds float, float* %0, i64 %17
%19 = load float, float* %18
%20 = getelementptr inbounds float, float* %3, i64 %13
%21 = load float, float* %20
%22 = fadd float %19, %21
store float %22, float* %20
%23 = or i64 %13, 1
%24 = getelementptr inbounds i32, i32* %2, i64 %23
%25 = load i32, i32* %24
%26 = sext i32 %25 to i64
%27 = getelementptr inbounds float, float* %0, i64 %26
%28 = load float, float* %27
%29 = getelementptr inbounds float, float* %3, i64 %23
%30 = load float, float* %29
%31 = fadd float %28, %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 i32, i32* %2, i64 %36
%40 = load i32, i32* %39
%41 = sext i32 %40 to i64
%42 = getelementptr inbounds float, float* %0, i64 %41
%43 = load float, float* %42
%44 = getelementptr inbounds float, float* %3, i64 %36
%45 = load float, float* %44
%46 = fadd float %43, %45
store float %46, float* %44
br label %47
47: ; preds = %38, %35, %4
ret void
}