#define NULL ((void *)0)
typedef unsigned long size_t; // Customize by platform.
typedef long intptr_t;
typedef unsigned long uintptr_t;
typedef long scalar_t__; // Either arithmetic or pointer type.
/* By default, we understand bool (as a convenience). */
typedef int bool;
#define false 0
#define true 1

/* Forward declarations */

/* Type definitions */

/* Variables and functions */
scalar_t__ exp(double);
double log(double);
double sumLogProbs(double *logprobs, int sz) {
  double max = 0.0;
  double p = 0.0;
  int i;
  for (i = 0; i < sz; ++i)
    if (i == 0 || logprobs[i] > max)
      max = logprobs[i];
  if (max == log(0))
    return max;
  for (i = 0; i < sz; ++i)
    p += exp(logprobs[i] - max);
  return max + log(p);
}


x86 assembly (clang, -O0):
<sumLogProbs>:
push   %rbp
mov    %rsp,%rbp
sub    $0x40,%rsp
mov    %rdi,-0x10(%rbp)
mov    %esi,-0x14(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x20(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x28(%rbp)
movl   $0x0,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    1199 <sumLogProbs+0x79>
cmpl   $0x0,-0x2c(%rbp)
je     1174 <sumLogProbs+0x54>
mov    -0x10(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
ucomisd -0x20(%rbp),%xmm0
jbe    1186 <sumLogProbs+0x66>
mov    -0x10(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x20(%rbp)
jmp    118b <sumLogProbs+0x6b>
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    1146 <sumLogProbs+0x26>
movsd  -0x20(%rbp),%xmm0
movsd  %xmm0,-0x38(%rbp)
xorps  %xmm0,%xmm0
call   1040 <log@plt>
movaps %xmm0,%xmm1
movsd  -0x38(%rbp),%xmm0
ucomisd %xmm1,%xmm0
jne    11d2 <sumLogProbs+0xb2>
jp     11d2 <sumLogProbs+0xb2>
movsd  -0x20(%rbp),%xmm0
movsd  %xmm0,-0x8(%rbp)
jmp    123e <sumLogProbs+0x11e>
movl   $0x0,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    1219 <sumLogProbs+0xf9>
mov    -0x10(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
subsd  -0x20(%rbp),%xmm0
call   1030 <exp@plt>
cvtsi2sd %rax,%xmm0
addsd  -0x28(%rbp),%xmm0
movsd  %xmm0,-0x28(%rbp)
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    11d9 <sumLogProbs+0xb9>
movsd  -0x20(%rbp),%xmm0
movsd  %xmm0,-0x40(%rbp)
movsd  -0x28(%rbp),%xmm0
call   1040 <log@plt>
movaps %xmm0,%xmm1
movsd  -0x40(%rbp),%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x8(%rbp)
movsd  -0x8(%rbp),%xmm0
add    $0x40,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<sumLogProbs>:
push   %rbp
push   %r15
push   %r14
push   %rbx
sub    $0x18,%rsp
test   %esi,%esi
jle    116a <sumLogProbs+0x4a>
mov    %esi,%r14d
mov    %rdi,%rbx
mov    %esi,%r15d
movsd  (%rdi),%xmm1
cmp    $0x1,%esi
jne    1190 <sumLogProbs+0x70>
xorpd  %xmm0,%xmm0
movsd  %xmm1,0x8(%rsp)
call   1040 <log@plt>
movsd  0x8(%rsp),%xmm1
ucomisd %xmm0,%xmm1
jne    123c <sumLogProbs+0x11c>
jnp    12a6 <sumLogProbs+0x186>
jmp    123c <sumLogProbs+0x11c>
xorpd  %xmm0,%xmm0
call   1040 <log@plt>
xorpd  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
xorpd  %xmm0,%xmm0
jne    1291 <sumLogProbs+0x171>
jnp    12a6 <sumLogProbs+0x186>
jmp    1291 <sumLogProbs+0x171>
lea    -0x1(%r15),%rdx
lea    -0x2(%r15),%rcx
mov    %edx,%eax
and    $0x3,%eax
cmp    $0x3,%rcx
jae    11aa <sumLogProbs+0x8a>
mov    $0x1,%ecx
jmp    11e5 <sumLogProbs+0xc5>
and    $0xfffffffffffffffc,%rdx
xor    %ecx,%ecx
movsd  0x8(%rbx,%rcx,8),%xmm0
maxsd  %xmm1,%xmm0
movsd  0x10(%rbx,%rcx,8),%xmm1
maxsd  %xmm0,%xmm1
movsd  0x18(%rbx,%rcx,8),%xmm0
maxsd  %xmm1,%xmm0
movsd  0x20(%rbx,%rcx,8),%xmm1
maxsd  %xmm0,%xmm1
add    $0x4,%rcx
cmp    %rcx,%rdx
jne    11b0 <sumLogProbs+0x90>
add    $0x1,%rcx
movapd %xmm1,%xmm0
test   %rax,%rax
je     1216 <sumLogProbs+0xf6>
lea    (%rbx,%rcx,8),%rcx
xor    %edx,%edx
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
movsd  (%rcx,%rdx,8),%xmm1
maxsd  %xmm0,%xmm1
add    $0x1,%rdx
movapd %xmm1,%xmm0
cmp    %rdx,%rax
jne    1200 <sumLogProbs+0xe0>
movsd  %xmm1,0x8(%rsp)
xorpd  %xmm0,%xmm0
call   1040 <log@plt>
movsd  0x8(%rsp),%xmm1
ucomisd %xmm0,%xmm1
jne    1233 <sumLogProbs+0x113>
jnp    12a6 <sumLogProbs+0x186>
test   %r14d,%r14d
xorpd  %xmm0,%xmm0
jle    1291 <sumLogProbs+0x171>
xorpd  %xmm0,%xmm0
xor    %ebp,%ebp
movsd  %xmm1,0x8(%rsp)
nopl   0x0(%rax,%rax,1)
movsd  %xmm0,0x10(%rsp)
movsd  (%rbx,%rbp,8),%xmm0
subsd  %xmm1,%xmm0
call   1030 <exp@plt>
movsd  0x8(%rsp),%xmm1
xorps  %xmm0,%xmm0
cvtsi2sd %rax,%xmm0
movsd  0x10(%rsp),%xmm2
addsd  %xmm0,%xmm2
movsd  %xmm2,0x10(%rsp)
movsd  0x10(%rsp),%xmm0
add    $0x1,%rbp
cmp    %rbp,%r15
jne    1250 <sumLogProbs+0x130>
movsd  %xmm1,0x8(%rsp)
call   1040 <log@plt>
movsd  0x8(%rsp),%xmm1
addsd  %xmm0,%xmm1
movapd %xmm1,%xmm0
add    $0x18,%rsp
pop    %rbx
pop    %r14
pop    %r15
pop    %rbp
ret

LLVM IR (-O0):
define double @sumLogProbs(double* %0, i32 %1) {
%3 = alloca double
%4 = alloca double*
%5 = alloca i32
%6 = alloca double
%7 = alloca double
%8 = alloca i32
store double* %0, double** %4
store i32 %1, i32* %5
store double 0.000000e+00, double* %6
store double 0.000000e+00, double* %7
store i32 0, i32* %8
br label %9
9: ; preds = %31, %2
%10 = load i32, i32* %8
%11 = load i32, i32* %5
%12 = icmp slt i32 %10, %11
br i1 %12, label %13, label %34
13: ; preds = %9
%14 = load i32, i32* %8
%15 = icmp eq i32 %14, 0
br i1 %15, label %24, label %16
16: ; preds = %13
%17 = load double*, double** %4
%18 = load i32, i32* %8
%19 = sext i32 %18 to i64
%20 = getelementptr inbounds double, double* %17, i64 %19
%21 = load double, double* %20
%22 = load double, double* %6
%23 = fcmp ogt double %21, %22
br i1 %23, label %24, label %30
24: ; preds = %16, %13
%25 = load double*, double** %4
%26 = load i32, i32* %8
%27 = sext i32 %26 to i64
%28 = getelementptr inbounds double, double* %25, i64 %27
%29 = load double, double* %28
store double %29, double* %6
br label %30
30: ; preds = %24, %16
br label %31
31: ; preds = %30
%32 = load i32, i32* %8
%33 = add nsw i32 %32, 1
store i32 %33, i32* %8
br label %9
34: ; preds = %9
%35 = load double, double* %6
%36 = call double @log(double 0.000000e+00)
%37 = fcmp oeq double %35, %36
br i1 %37, label %38, label %40
38: ; preds = %34
%39 = load double, double* %6
store double %39, double* %3
br label %65
40: ; preds = %34
store i32 0, i32* %8
br label %41
41: ; preds = %57, %40
%42 = load i32, i32* %8
%43 = load i32, i32* %5
%44 = icmp slt i32 %42, %43
br i1 %44, label %45, label %60
45: ; preds = %41
%46 = load double*, double** %4
%47 = load i32, i32* %8
%48 = sext i32 %47 to i64
%49 = getelementptr inbounds double, double* %46, i64 %48
%50 = load double, double* %49
%51 = load double, double* %6
%52 = fsub double %50, %51
%53 = call i64 @exp(double %52)
%54 = sitofp i64 %53 to double
%55 = load double, double* %7
%56 = fadd double %55, %54
store double %56, double* %7
br label %57
57: ; preds = %45
%58 = load i32, i32* %8
%59 = add nsw i32 %58, 1
store i32 %59, i32* %8
br label %41
60: ; preds = %41
%61 = load double, double* %6
%62 = load double, double* %7
%63 = call double @log(double %62)
%64 = fadd double %61, %63
store double %64, double* %3
br label %65
65: ; preds = %60, %38
%66 = load double, double* %3
ret double %66
}

LLVM IR (-O2):
define double @sumLogProbs(double* nocapture %0, i32 %1) {
%3 = icmp sgt i32 %1, 0
br i1 %3, label %4, label %64
4: ; preds = %2
%5 = zext i32 %1 to i64
%6 = load double, double* %0
%7 = icmp eq i32 %1, 1
br i1 %7, label %61, label %8
8: ; preds = %4
%9 = add nsw i64 %5, -1
%10 = add nsw i64 %5, -2
%11 = and i64 %9, 3
%12 = icmp ult i64 %10, 3
br i1 %12, label %41, label %13
13: ; preds = %8
%14 = and i64 %9, -4
br label %15
15: ; preds = %15, %13
%16 = phi i64 [ 1, %13 ], [ %38, %15 ]
%17 = phi double [ %6, %13 ], [ %37, %15 ]
%18 = phi i64 [ 0, %13 ], [ %39, %15 ]
%19 = getelementptr inbounds double, double* %0, i64 %16
%20 = load double, double* %19
%21 = fcmp ogt double %20, %17
%22 = select i1 %21, double %20, double %17
%23 = add nuw nsw i64 %16, 1
%24 = getelementptr inbounds double, double* %0, i64 %23
%25 = load double, double* %24
%26 = fcmp ogt double %25, %22
%27 = select i1 %26, double %25, double %22
%28 = add nuw nsw i64 %16, 2
%29 = getelementptr inbounds double, double* %0, i64 %28
%30 = load double, double* %29
%31 = fcmp ogt double %30, %27
%32 = select i1 %31, double %30, double %27
%33 = add nuw nsw i64 %16, 3
%34 = getelementptr inbounds double, double* %0, i64 %33
%35 = load double, double* %34
%36 = fcmp ogt double %35, %32
%37 = select i1 %36, double %35, double %32
%38 = add nuw nsw i64 %16, 4
%39 = add i64 %18, 4
%40 = icmp eq i64 %39, %14
br i1 %40, label %41, label %15
41: ; preds = %15, %8
%42 = phi double [ undef, %8 ], [ %37, %15 ]
%43 = phi i64 [ 1, %8 ], [ %38, %15 ]
%44 = phi double [ %6, %8 ], [ %37, %15 ]
%45 = icmp eq i64 %11, 0
br i1 %45, label %57, label %46
46: ; preds = %41, %46
%47 = phi i64 [ %54, %46 ], [ %43, %41 ]
%48 = phi double [ %53, %46 ], [ %44, %41 ]
%49 = phi i64 [ %55, %46 ], [ 0, %41 ]
%50 = getelementptr inbounds double, double* %0, i64 %47
%51 = load double, double* %50
%52 = fcmp ogt double %51, %48
%53 = select i1 %52, double %51, double %48
%54 = add nuw nsw i64 %47, 1
%55 = add i64 %49, 1
%56 = icmp eq i64 %55, %11
br i1 %56, label %57, label %46
57: ; preds = %46, %41
%58 = phi double [ %42, %41 ], [ %53, %46 ]
%59 = tail call double @log(double 0.000000e+00)
%60 = fcmp oeq double %58, %59
br i1 %60, label %87, label %67
61: ; preds = %4
%62 = tail call double @log(double 0.000000e+00)
%63 = fcmp oeq double %6, %62
br i1 %63, label %87, label %68
64: ; preds = %2
%65 = tail call double @log(double 0.000000e+00)
%66 = fcmp oeq double %65, 0.000000e+00
br i1 %66, label %87, label %82
67: ; preds = %57
br i1 %3, label %68, label %82
68: ; preds = %61, %67
%69 = phi double [ %58, %67 ], [ %6, %61 ]
%70 = zext i32 %1 to i64
br label %71
71: ; preds = %68, %71
%72 = phi i64 [ 0, %68 ], [ %80, %71 ]
%73 = phi double [ 0.000000e+00, %68 ], [ %79, %71 ]
%74 = getelementptr inbounds double, double* %0, i64 %72
%75 = load double, double* %74
%76 = fsub double %75, %69
%77 = tail call i64 @exp(double %76)
%78 = sitofp i64 %77 to double
%79 = fadd double %73, %78
%80 = add nuw nsw i64 %72, 1
%81 = icmp eq i64 %80, %70
br i1 %81, label %82, label %71
82: ; preds = %71, %64, %67
%83 = phi double [ %58, %67 ], [ 0.000000e+00, %64 ], [ %69, %71 ]
%84 = phi double [ 0.000000e+00, %67 ], [ 0.000000e+00, %64 ], [ %79, %71 ]
%85 = tail call double @log(double %84)
%86 = fadd double %83, %85
br label %87
87: ; preds = %61, %64, %57, %82
%88 = phi double [ %86, %82 ], [ %58, %57 ], [ 0.000000e+00, %64 ], [ %6, %61 ]
ret double %88
}