float sum_AB_lf(float *AB, float phi[], int k, int nsample, int nparam) {
  float sum;
  int i;
  sum = 0.;
  for (i = 0; i < nparam; i++)
    sum += phi[i] * AB[k * nparam + i];
  return (sum);
}


x86 assembly (clang, -O0):
<sum_AB_lf>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %edx,-0x14(%rbp)
mov    %ecx,-0x18(%rbp)
mov    %r8d,-0x1c(%rbp)
xorps  %xmm0,%xmm0
movss  %xmm0,-0x20(%rbp)
movl   $0x0,-0x24(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0x1c(%rbp),%eax
jge    1174 <sum_AB_lf+0x74>
mov    -0x10(%rbp),%rax
movslq -0x24(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x14(%rbp),%ecx
imul   -0x1c(%rbp),%ecx
add    -0x24(%rbp),%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm2
movss  -0x20(%rbp),%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,-0x20(%rbp)
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    1125 <sum_AB_lf+0x25>
movss  -0x20(%rbp),%xmm0
pop    %rbp
ret

x86 assembly (clang, -O2):
<sum_AB_lf>:
test   %r8d,%r8d
jle    1127 <sum_AB_lf+0x27>
imul   %r8d,%edx
movslq %edx,%r9
mov    %r8d,%eax
lea    -0x1(%rax),%rcx
mov    %eax,%r8d
and    $0x3,%r8d
cmp    $0x3,%rcx
jae    112b <sum_AB_lf+0x2b>
xorps  %xmm0,%xmm0
xor    %edx,%edx
jmp    1187 <sum_AB_lf+0x87>
xorps  %xmm0,%xmm0
ret
and    $0xfffffffc,%eax
lea    (%rdi,%r9,4),%rcx
add    $0xc,%rcx
xorps  %xmm0,%xmm0
xor    %edx,%edx
nopl   0x0(%rax,%rax,1)
movss  (%rsi,%rdx,4),%xmm1
movss  0x4(%rsi,%rdx,4),%xmm2
mulss  -0xc(%rcx,%rdx,4),%xmm1
mulss  -0x8(%rcx,%rdx,4),%xmm2
addss  %xmm0,%xmm1
movss  0x8(%rsi,%rdx,4),%xmm3
mulss  -0x4(%rcx,%rdx,4),%xmm3
addss  %xmm1,%xmm2
movss  0xc(%rsi,%rdx,4),%xmm0
mulss  (%rcx,%rdx,4),%xmm0
addss  %xmm2,%xmm3
addss  %xmm3,%xmm0
add    $0x4,%rdx
cmp    %rdx,%rax
jne    1140 <sum_AB_lf+0x40>
test   %r8,%r8
je     11b7 <sum_AB_lf+0xb7>
add    %rdx,%r9
lea    (%rdi,%r9,4),%rax
lea    (%rsi,%rdx,4),%rcx
xor    %edx,%edx
nopl   0x0(%rax)
movss  (%rcx,%rdx,4),%xmm1
mulss  (%rax,%rdx,4),%xmm1
addss  %xmm1,%xmm0
add    $0x1,%rdx
cmp    %rdx,%r8
jne    11a0 <sum_AB_lf+0xa0>
ret

LLVM IR (-O0):
define float @sum_AB_lf(float* %0, float* %1, i32 %2, i32 %3, i32 %4) {
%6 = alloca float*
%7 = alloca float*
%8 = alloca i32
%9 = alloca i32
%10 = alloca i32
%11 = alloca float
%12 = alloca i32
store float* %0, float** %6
store float* %1, float** %7
store i32 %2, i32* %8
store i32 %3, i32* %9
store i32 %4, i32* %10
store float 0.000000e+00, float* %11
store i32 0, i32* %12
br label %13
13: ; preds = %34, %5
%14 = load i32, i32* %12
%15 = load i32, i32* %10
%16 = icmp slt i32 %14, %15
br i1 %16, label %17, label %37
17: ; preds = %13
%18 = load float*, float** %7
%19 = load i32, i32* %12
%20 = sext i32 %19 to i64
%21 = getelementptr inbounds float, float* %18, i64 %20
%22 = load float, float* %21
%23 = load float*, float** %6
%24 = load i32, i32* %8
%25 = load i32, i32* %10
%26 = mul nsw i32 %24, %25
%27 = load i32, i32* %12
%28 = add nsw i32 %26, %27
%29 = sext i32 %28 to i64
%30 = getelementptr inbounds float, float* %23, i64 %29
%31 = load float, float* %30
%32 = load float, float* %11
%33 = call float @llvm.fmuladd.f32(float %22, float %31, float %32)
store float %33, float* %11
br label %34
34: ; preds = %17
%35 = load i32, i32* %12
%36 = add nsw i32 %35, 1
store i32 %36, i32* %12
br label %13
37: ; preds = %13
%38 = load float, float* %11
ret float %38
}

LLVM IR (-O2):
define float @sum_AB_lf(float* nocapture %0, float* nocapture %1, i32 %2, i32 %3, i32 %4) {
%6 = icmp sgt i32 %4, 0
br i1 %6, label %7, label %68
7: ; preds = %5
%8 = mul nsw i32 %4, %2
%9 = sext i32 %8 to i64
%10 = zext i32 %4 to i64
%11 = add nsw i64 %10, -1
%12 = and i64 %10, 3
%13 = icmp ult i64 %11, 3
br i1 %13, label %50, label %14
14: ; preds = %7
%15 = and i64 %10, 4294967292
br label %16
16: ; preds = %16, %14
%17 = phi i64 [ 0, %14 ], [ %47, %16 ]
%18 = phi float [ 0.000000e+00, %14 ], [ %46, %16 ]
%19 = phi i64 [ 0, %14 ], [ %48, %16 ]
%20 = getelementptr inbounds float, float* %1, i64 %17
%21 = load float, float* %20
%22 = add nsw i64 %17, %9
%23 = getelementptr inbounds float, float* %0, i64 %22
%24 = load float, float* %23
%25 = tail call float @llvm.fmuladd.f32(float %21, float %24, float %18)
%26 = or i64 %17, 1
%27 = getelementptr inbounds float, float* %1, i64 %26
%28 = load float, float* %27
%29 = add nsw i64 %26, %9
%30 = getelementptr inbounds float, float* %0, i64 %29
%31 = load float, float* %30
%32 = tail call float @llvm.fmuladd.f32(float %28, float %31, float %25)
%33 = or i64 %17, 2
%34 = getelementptr inbounds float, float* %1, i64 %33
%35 = load float, float* %34
%36 = add nsw i64 %33, %9
%37 = getelementptr inbounds float, float* %0, i64 %36
%38 = load float, float* %37
%39 = tail call float @llvm.fmuladd.f32(float %35, float %38, float %32)
%40 = or i64 %17, 3
%41 = getelementptr inbounds float, float* %1, i64 %40
%42 = load float, float* %41
%43 = add nsw i64 %40, %9
%44 = getelementptr inbounds float, float* %0, i64 %43
%45 = load float, float* %44
%46 = tail call float @llvm.fmuladd.f32(float %42, float %45, float %39)
%47 = add nuw nsw i64 %17, 4
%48 = add i64 %19, 4
%49 = icmp eq i64 %48, %15
br i1 %49, label %50, label %16
50: ; preds = %16, %7
%51 = phi float [ undef, %7 ], [ %46, %16 ]
%52 = phi i64 [ 0, %7 ], [ %47, %16 ]
%53 = phi float [ 0.000000e+00, %7 ], [ %46, %16 ]
%54 = icmp eq i64 %12, 0
br i1 %54, label %68, label %55
55: ; preds = %50, %55
%56 = phi i64 [ %65, %55 ], [ %52, %50 ]
%57 = phi float [ %64, %55 ], [ %53, %50 ]
%58 = phi i64 [ %66, %55 ], [ 0, %50 ]
%59 = getelementptr inbounds float, float* %1, i64 %56
%60 = load float, float* %59
%61 = add nsw i64 %56, %9
%62 = getelementptr inbounds float, float* %0, i64 %61
%63 = load float, float* %62
%64 = tail call float @llvm.fmuladd.f32(float %60, float %63, float %57)
%65 = add nuw nsw i64 %56, 1
%66 = add i64 %58, 1
%67 = icmp eq i64 %66, %12
br i1 %67, label %68, label %55
68: ; preds = %50, %55, %5
%69 = phi float [ 0.000000e+00, %5 ], [ %51, %50 ], [ %64, %55 ]
ret float %69
}