#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 */
double _logsumexp_pair(double, double);
double _logsumexp_sort_kahan_inplace(double *, int);
double log(double);
void _tram_update_lagrangian_mult(double *log_lagrangian_mult,
                                  double *biased_conf_energies,
                                  int *count_matrices, int *state_counts,
                                  int n_therm_states, int n_conf_states,
                                  double *scratch_M,
                                  double *new_log_lagrangian_mult) {
  int i, j, K, o;
  int Ki, Kj, KM, KMM;
  int CK, CKij;
  double divisor;
  for (K = 0; K < n_therm_states; ++K) {
    KM = K * n_conf_states;
    KMM = KM * n_conf_states;
    for (i = 0; i < n_conf_states; ++i) {
      Ki = KM + i;
      if (0 == state_counts[Ki]) {
        new_log_lagrangian_mult[Ki] = -(1.0 / 0.0);
        continue;
      }
      o = 0;
      for (j = 0; j < n_conf_states; ++j) {
        CKij = count_matrices[KMM + i * n_conf_states + j];
        if (i == j) {
          scratch_M[o++] = (0 == CKij) ? -(1.0 / 0.0) : log(0 + (double)CKij);
          continue;
        }
        CK = CKij + count_matrices[KMM + j * n_conf_states + i];
        if (0 == CK)
          continue;
        Kj = KM + j;
        divisor = _logsumexp_pair(
            log_lagrangian_mult[Kj] - biased_conf_energies[Ki] -
                log_lagrangian_mult[Ki] + biased_conf_energies[Kj],
            0.0);
        scratch_M[o++] = log((double)CK) - divisor;
      }
      new_log_lagrangian_mult[Ki] = _logsumexp_sort_kahan_inplace(scratch_M, o);
    }
  }
}


x86 assembly (clang, -O0):
<_tram_update_lagrangian_mult>:
push   %rbp
mov    %rsp,%rbp
sub    $0x60,%rsp
mov    0x18(%rbp),%rax
mov    0x10(%rbp),%rax
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %rcx,-0x20(%rbp)
mov    %r8d,-0x24(%rbp)
mov    %r9d,-0x28(%rbp)
movl   $0x0,-0x34(%rbp)
mov    -0x34(%rbp),%eax
cmp    -0x24(%rbp),%eax
jge    1345 <_tram_update_lagrangian_mult+0x215>
mov    -0x34(%rbp),%eax
imul   -0x28(%rbp),%eax
mov    %eax,-0x44(%rbp)
mov    -0x44(%rbp),%eax
imul   -0x28(%rbp),%eax
mov    %eax,-0x48(%rbp)
movl   $0x0,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0x28(%rbp),%eax
jge    1332 <_tram_update_lagrangian_mult+0x202>
mov    -0x44(%rbp),%eax
add    -0x2c(%rbp),%eax
mov    %eax,-0x3c(%rbp)
mov    -0x20(%rbp),%rcx
movslq -0x3c(%rbp),%rdx
xor    %eax,%eax
cmp    (%rcx,%rdx,4),%eax
jne    11c8 <_tram_update_lagrangian_mult+0x98>
mov    0x18(%rbp),%rax
movslq -0x3c(%rbp),%rcx
movsd  0xe42(%rip),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
jmp    1324 <_tram_update_lagrangian_mult+0x1f4>
movl   $0x0,-0x38(%rbp)
movl   $0x0,-0x30(%rbp)
mov    -0x30(%rbp),%eax
cmp    -0x28(%rbp),%eax
jge    130b <_tram_update_lagrangian_mult+0x1db>
mov    -0x18(%rbp),%rax
mov    -0x48(%rbp),%ecx
mov    -0x2c(%rbp),%edx
imul   -0x28(%rbp),%edx
add    %edx,%ecx
add    -0x30(%rbp),%ecx
movslq %ecx,%rcx
mov    (%rax,%rcx,4),%eax
mov    %eax,-0x50(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0x30(%rbp),%eax
jne    125e <_tram_update_lagrangian_mult+0x12e>
xor    %eax,%eax
cmp    -0x50(%rbp),%eax
jne    1227 <_tram_update_lagrangian_mult+0xf7>
movsd  0xde3(%rip),%xmm0
movsd  %xmm0,-0x60(%rbp)
jmp    123d <_tram_update_lagrangian_mult+0x10d>
cvtsi2sdl -0x50(%rbp),%xmm1
xorps  %xmm0,%xmm0
addsd  %xmm1,%xmm0
call   1050 <log@plt>
movsd  %xmm0,-0x60(%rbp)
movsd  -0x60(%rbp),%xmm0
mov    0x10(%rbp),%rax
mov    -0x38(%rbp),%ecx
mov    %ecx,%edx
add    $0x1,%edx
mov    %edx,-0x38(%rbp)
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
jmp    12fd <_tram_update_lagrangian_mult+0x1cd>
mov    -0x50(%rbp),%eax
mov    -0x18(%rbp),%rcx
mov    -0x48(%rbp),%edx
mov    -0x30(%rbp),%esi
imul   -0x28(%rbp),%esi
add    %esi,%edx
add    -0x2c(%rbp),%edx
movslq %edx,%rdx
add    (%rcx,%rdx,4),%eax
mov    %eax,-0x4c(%rbp)
xor    %eax,%eax
cmp    -0x4c(%rbp),%eax
jne    128d <_tram_update_lagrangian_mult+0x15d>
jmp    12fd <_tram_update_lagrangian_mult+0x1cd>
mov    -0x44(%rbp),%eax
add    -0x30(%rbp),%eax
mov    %eax,-0x40(%rbp)
mov    -0x8(%rbp),%rax
movslq -0x40(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x3c(%rbp),%rcx
subsd  (%rax,%rcx,8),%xmm0
mov    -0x8(%rbp),%rax
movslq -0x3c(%rbp),%rcx
subsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x40(%rbp),%rcx
addsd  (%rax,%rcx,8),%xmm0
xorps  %xmm1,%xmm1
call   1040 <_logsumexp_pair@plt>
movsd  %xmm0,-0x58(%rbp)
cvtsi2sdl -0x4c(%rbp),%xmm0
call   1050 <log@plt>
subsd  -0x58(%rbp),%xmm0
mov    0x10(%rbp),%rax
mov    -0x38(%rbp),%ecx
mov    %ecx,%edx
add    $0x1,%edx
mov    %edx,-0x38(%rbp)
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
jmp    11d6 <_tram_update_lagrangian_mult+0xa6>
mov    0x10(%rbp),%rdi
mov    -0x38(%rbp),%esi
call   1030 <_logsumexp_sort_kahan_inplace@plt>
mov    0x18(%rbp),%rax
movslq -0x3c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    1186 <_tram_update_lagrangian_mult+0x56>
jmp    1337 <_tram_update_lagrangian_mult+0x207>
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    115f <_tram_update_lagrangian_mult+0x2f>
add    $0x60,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<_tram_update_lagrangian_mult>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0xa8,%rsp
mov    %r9d,0xc(%rsp)
mov    %rcx,0x90(%rsp)
mov    %rdx,0x58(%rsp)
mov    %rsi,0x38(%rsp)
mov    %rdi,0x30(%rsp)
test   %r8d,%r8d
jle    13a7 <_tram_update_lagrangian_mult+0x277>
movslq 0xc(%rsp),%rax
mov    %r8d,%ecx
mov    %rcx,0x80(%rsp)
mov    %eax,%ecx
mov    %rcx,0x60(%rsp)
mov    %eax,%ecx
imul   %eax,%ecx
mov    %ecx,0x28(%rsp)
mov    %rax,%rcx
imul   %rax,%rcx
shl    $0x2,%rcx
mov    %rcx,0x78(%rsp)
lea    0x0(,%rax,4),%rcx
mov    %rcx,0x88(%rsp)
mov    %rax,0x68(%rsp)
lea    0x0(,%rax,8),%rax
mov    %rax,0x70(%rsp)
xor    %eax,%eax
mov    %rax,0x20(%rsp)
mov    0x38(%rsp),%rax
mov    %rax,0x48(%rsp)
mov    0x30(%rsp),%rax
mov    %rax,0x40(%rsp)
mov    0x58(%rsp),%rax
mov    %rax,0x18(%rsp)
xor    %eax,%eax
mov    %rax,0x10(%rsp)
jmp    123b <_tram_update_lagrangian_mult+0x10b>
cs nopw 0x0(%rax,%rax,1)
mov    0x10(%rsp),%rcx
add    $0x1,%rcx
mov    0x20(%rsp),%rax
add    0x28(%rsp),%eax
mov    %rax,0x20(%rsp)
mov    0x18(%rsp),%rax
add    0x78(%rsp),%rax
mov    %rax,0x18(%rsp)
mov    0x70(%rsp),%rax
add    %rax,0x40(%rsp)
add    %rax,0x48(%rsp)
mov    %rcx,%rax
mov    %rcx,0x10(%rsp)
cmp    0x80(%rsp),%rcx
je     13a7 <_tram_update_lagrangian_mult+0x277>
cmpl   $0x0,0xc(%rsp)
jle    11f0 <_tram_update_lagrangian_mult+0xc0>
mov    0x10(%rsp),%rax
imul   0x68(%rsp),%rax
mov    %rax,0x98(%rsp)
mov    0x18(%rsp),%r13
mov    0x20(%rsp),%rax
mov    %eax,%r12d
xor    %ebp,%ebp
jmp    12b2 <_tram_update_lagrangian_mult+0x182>
cs nopw 0x0(%rax,%rax,1)
mov    0xe0(%rsp),%rdi
mov    %r14d,%esi
call   1030 <_logsumexp_sort_kahan_inplace@plt>
mov    0x2c(%rsp),%r12d
mov    0xe8(%rsp),%rax
mov    0x50(%rsp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
add    $0x1,%rbp
add    $0x1,%r12d
add    0x88(%rsp),%r13
cmp    0x60(%rsp),%rbp
je     11f0 <_tram_update_lagrangian_mult+0xc0>
mov    0x98(%rsp),%rax
lea    (%rax,%rbp,1),%rcx
mov    0x90(%rsp),%rax
mov    %rcx,0x50(%rsp)
cmpl   $0x0,(%rax,%rcx,4)
movsd  0xd29(%rip),%xmm0
je     1285 <_tram_update_lagrangian_mult+0x155>
xor    %ebx,%ebx
mov    %r12d,0x2c(%rsp)
xor    %r14d,%r14d
jmp    1318 <_tram_update_lagrangian_mult+0x1e8>
cs nopw 0x0(%rax,%rax,1)
nop
movslq %r14d,%rax
add    $0x1,%r14d
mov    0xe0(%rsp),%rcx
movsd  %xmm0,(%rcx,%rax,8)
add    $0x1,%rbx
add    0xc(%rsp),%r12d
cmp    %rbx,0x60(%rsp)
je     1270 <_tram_update_lagrangian_mult+0x140>
mov    0x0(%r13,%rbx,4),%r15d
cmp    %rbx,%rbp
jne    1340 <_tram_update_lagrangian_mult+0x210>
movsd  0xcd6(%rip),%xmm0
test   %r15d,%r15d
je     12f0 <_tram_update_lagrangian_mult+0x1c0>
xorps  %xmm0,%xmm0
cvtsi2sd %r15d,%xmm0
call   1050 <log@plt>
jmp    12f0 <_tram_update_lagrangian_mult+0x1c0>
xchg   %ax,%ax
movslq %r12d,%rax
mov    0x58(%rsp),%rcx
add    (%rcx,%rax,4),%r15d
je     1304 <_tram_update_lagrangian_mult+0x1d4>
mov    0x40(%rsp),%rax
movsd  (%rax,%rbx,8),%xmm0
mov    0x38(%rsp),%rax
mov    0x50(%rsp),%rcx
subsd  (%rax,%rcx,8),%xmm0
mov    0x30(%rsp),%rax
subsd  (%rax,%rcx,8),%xmm0
mov    0x48(%rsp),%rax
addsd  (%rax,%rbx,8),%xmm0
xorps  %xmm1,%xmm1
call   1040 <_logsumexp_pair@plt>
movsd  %xmm0,0xa0(%rsp)
xorps  %xmm0,%xmm0
cvtsi2sd %r15d,%xmm0
call   1050 <log@plt>
subsd  0xa0(%rsp),%xmm0
jmp    12f0 <_tram_update_lagrangian_mult+0x1c0>
add    $0xa8,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

LLVM IR (-O0):
define void @_tram_update_lagrangian_mult(double* %0, double* %1, i32* %2, i32* %3, i32 %4, i32 %5, double* %6, double* %7) {
%9 = alloca double*
%10 = alloca double*
%11 = alloca i32*
%12 = alloca i32*
%13 = alloca i32
%14 = alloca i32
%15 = alloca double*
%16 = alloca double*
%17 = alloca i32
%18 = alloca i32
%19 = alloca i32
%20 = alloca i32
%21 = alloca i32
%22 = alloca i32
%23 = alloca i32
%24 = alloca i32
%25 = alloca i32
%26 = alloca i32
%27 = alloca double
store double* %0, double** %9
store double* %1, double** %10
store i32* %2, i32** %11
store i32* %3, i32** %12
store i32 %4, i32* %13
store i32 %5, i32* %14
store double* %6, double** %15
store double* %7, double** %16
store i32 0, i32* %19
br label %28
28: ; preds = %164, %8
%29 = load i32, i32* %19
%30 = load i32, i32* %13
%31 = icmp slt i32 %29, %30
br i1 %31, label %32, label %167
32: ; preds = %28
%33 = load i32, i32* %19
%34 = load i32, i32* %14
%35 = mul nsw i32 %33, %34
store i32 %35, i32* %23
%36 = load i32, i32* %23
%37 = load i32, i32* %14
%38 = mul nsw i32 %36, %37
store i32 %38, i32* %24
store i32 0, i32* %17
br label %39
39: ; preds = %160, %32
%40 = load i32, i32* %17
%41 = load i32, i32* %14
%42 = icmp slt i32 %40, %41
br i1 %42, label %43, label %163
43: ; preds = %39
%44 = load i32, i32* %23
%45 = load i32, i32* %17
%46 = add nsw i32 %44, %45
store i32 %46, i32* %21
%47 = load i32*, i32** %12
%48 = load i32, i32* %21
%49 = sext i32 %48 to i64
%50 = getelementptr inbounds i32, i32* %47, i64 %49
%51 = load i32, i32* %50
%52 = icmp eq i32 0, %51
br i1 %52, label %53, label %58
53: ; preds = %43
%54 = load double*, double** %16
%55 = load i32, i32* %21
%56 = sext i32 %55 to i64
%57 = getelementptr inbounds double, double* %54, i64 %56
store double 0xFFF0000000000000, double* %57
br label %160
58: ; preds = %43
store i32 0, i32* %20
store i32 0, i32* %18
br label %59
59: ; preds = %149, %58
%60 = load i32, i32* %18
%61 = load i32, i32* %14
%62 = icmp slt i32 %60, %61
br i1 %62, label %63, label %152
63: ; preds = %59
%64 = load i32*, i32** %11
%65 = load i32, i32* %24
%66 = load i32, i32* %17
%67 = load i32, i32* %14
%68 = mul nsw i32 %66, %67
%69 = add nsw i32 %65, %68
%70 = load i32, i32* %18
%71 = add nsw i32 %69, %70
%72 = sext i32 %71 to i64
%73 = getelementptr inbounds i32, i32* %64, i64 %72
%74 = load i32, i32* %73
store i32 %74, i32* %26
%75 = load i32, i32* %17
%76 = load i32, i32* %18
%77 = icmp eq i32 %75, %76
br i1 %77, label %78, label %94
78: ; preds = %63
%79 = load i32, i32* %26
%80 = icmp eq i32 0, %79
br i1 %80, label %81, label %82
81: ; preds = %78
br label %87
82: ; preds = %78
%83 = load i32, i32* %26
%84 = sitofp i32 %83 to double
%85 = fadd double 0.000000e+00, %84
%86 = call double @log(double %85)
br label %87
87: ; preds = %82, %81
%88 = phi double [ 0xFFF0000000000000, %81 ], [ %86, %82 ]
%89 = load double*, double** %15
%90 = load i32, i32* %20
%91 = add nsw i32 %90, 1
store i32 %91, i32* %20
%92 = sext i32 %90 to i64
%93 = getelementptr inbounds double, double* %89, i64 %92
store double %88, double* %93
br label %149
94: ; preds = %63
%95 = load i32, i32* %26
%96 = load i32*, i32** %11
%97 = load i32, i32* %24
%98 = load i32, i32* %18
%99 = load i32, i32* %14
%100 = mul nsw i32 %98, %99
%101 = add nsw i32 %97, %100
%102 = load i32, i32* %17
%103 = add nsw i32 %101, %102
%104 = sext i32 %103 to i64
%105 = getelementptr inbounds i32, i32* %96, i64 %104
%106 = load i32, i32* %105
%107 = add nsw i32 %95, %106
store i32 %107, i32* %25
%108 = load i32, i32* %25
%109 = icmp eq i32 0, %108
br i1 %109, label %110, label %111
110: ; preds = %94
br label %149
111: ; preds = %94
%112 = load i32, i32* %23
%113 = load i32, i32* %18
%114 = add nsw i32 %112, %113
store i32 %114, i32* %22
%115 = load double*, double** %9
%116 = load i32, i32* %22
%117 = sext i32 %116 to i64
%118 = getelementptr inbounds double, double* %115, i64 %117
%119 = load double, double* %118
%120 = load double*, double** %10
%121 = load i32, i32* %21
%122 = sext i32 %121 to i64
%123 = getelementptr inbounds double, double* %120, i64 %122
%124 = load double, double* %123
%125 = fsub double %119, %124
%126 = load double*, double** %9
%127 = load i32, i32* %21
%128 = sext i32 %127 to i64
%129 = getelementptr inbounds double, double* %126, i64 %128
%130 = load double, double* %129
%131 = fsub double %125, %130
%132 = load double*, double** %10
%133 = load i32, i32* %22
%134 = sext i32 %133 to i64
%135 = getelementptr inbounds double, double* %132, i64 %134
%136 = load double, double* %135
%137 = fadd double %131, %136
%138 = call double @_logsumexp_pair(double %137, double 0.000000e+00)
store double %138, double* %27
%139 = load i32, i32* %25
%140 = sitofp i32 %139 to double
%141 = call double @log(double %140)
%142 = load double, double* %27
%143 = fsub double %141, %142
%144 = load double*, double** %15
%145 = load i32, i32* %20
%146 = add nsw i32 %145, 1
store i32 %146, i32* %20
%147 = sext i32 %145 to i64
%148 = getelementptr inbounds double, double* %144, i64 %147
store double %143, double* %148
br label %149
149: ; preds = %111, %110, %87
%150 = load i32, i32* %18
%151 = add nsw i32 %150, 1
store i32 %151, i32* %18
br label %59
152: ; preds = %59
%153 = load double*, double** %15
%154 = load i32, i32* %20
%155 = call double @_logsumexp_sort_kahan_inplace(double* %153, i32 %154)
%156 = load double*, double** %16
%157 = load i32, i32* %21
%158 = sext i32 %157 to i64
%159 = getelementptr inbounds double, double* %156, i64 %158
store double %155, double* %159
br label %160
160: ; preds = %152, %53
%161 = load i32, i32* %17
%162 = add nsw i32 %161, 1
store i32 %162, i32* %17
br label %39
163: ; preds = %39
br label %164
164: ; preds = %163
%165 = load i32, i32* %19
%166 = add nsw i32 %165, 1
store i32 %166, i32* %19
br label %28
167: ; preds = %28
ret void
}

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