#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 HUGE_VAL;
int /*<<< orphan*/ exp(double const);
double pow(int /*<<< orphan*/, double const);
double lcfit2_compute_weights(const size_t n, const double *lnl,
                              const double alpha, double *w) {
  double max_lnl = -HUGE_VAL;
  for (size_t i = 0; i < n; ++i) {
    if (lnl[i] > max_lnl) {
      max_lnl = lnl[i];
    }
  }
  for (size_t i = 0; i < n; ++i) {
    w[i] = pow(exp(lnl[i] - max_lnl), alpha);
  }
  return max_lnl;
}


x86 assembly (clang, -O0):
<lcfit2_compute_weights>:
push   %rbp
mov    %rsp,%rbp
sub    $0x40,%rsp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
movsd  %xmm0,-0x18(%rbp)
mov    %rdx,-0x20(%rbp)
mov    0x2ea8(%rip),%rax
movsd  (%rax),%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
movsd  %xmm0,-0x28(%rbp)
movq   $0x0,-0x30(%rbp)
mov    -0x30(%rbp),%rax
cmp    -0x8(%rbp),%rax
jae    11b6 <lcfit2_compute_weights+0x96>
mov    -0x10(%rbp),%rax
mov    -0x30(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
ucomisd -0x28(%rbp),%xmm0
jbe    11a0 <lcfit2_compute_weights+0x80>
mov    -0x10(%rbp),%rax
mov    -0x30(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x28(%rbp)
jmp    11a5 <lcfit2_compute_weights+0x85>
mov    -0x30(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x30(%rbp)
jmp    1168 <lcfit2_compute_weights+0x48>
movq   $0x0,-0x38(%rbp)
mov    -0x38(%rbp),%rax
cmp    -0x8(%rbp),%rax
jae    120d <lcfit2_compute_weights+0xed>
mov    -0x10(%rbp),%rax
mov    -0x38(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
subsd  -0x28(%rbp),%xmm0
call   1030 <exp@plt>
mov    %eax,%edi
movsd  -0x18(%rbp),%xmm0
call   1040 <pow@plt>
mov    -0x20(%rbp),%rax
mov    -0x38(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x38(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x38(%rbp)
jmp    11be <lcfit2_compute_weights+0x9e>
movsd  -0x28(%rbp),%xmm0
add    $0x40,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<lcfit2_compute_weights>:
push   %r15
push   %r14
push   %r12
push   %rbx
sub    $0x18,%rsp
movaps %xmm0,%xmm2
mov    0x2eb3(%rip),%rax
movsd  (%rax),%xmm0
xorpd  0xebf(%rip),%xmm0
test   %rdi,%rdi
je     1225 <lcfit2_compute_weights+0x105>
mov    %rdx,%r14
mov    %rsi,%r12
mov    %rdi,%r15
lea    -0x1(%rdi),%rcx
mov    %r15d,%eax
and    $0x3,%eax
cmp    $0x3,%rcx
jae    1167 <lcfit2_compute_weights+0x47>
xor    %ecx,%ecx
jmp    11a4 <lcfit2_compute_weights+0x84>
mov    %r15,%rdx
and    $0xfffffffffffffffc,%rdx
xor    %ecx,%ecx
movsd  (%r12,%rcx,8),%xmm1
maxsd  %xmm0,%xmm1
movsd  0x8(%r12,%rcx,8),%xmm0
maxsd  %xmm1,%xmm0
movsd  0x10(%r12,%rcx,8),%xmm1
maxsd  %xmm0,%xmm1
movsd  0x18(%r12,%rcx,8),%xmm0
maxsd  %xmm1,%xmm0
add    $0x4,%rcx
cmp    %rcx,%rdx
jne    1170 <lcfit2_compute_weights+0x50>
movapd %xmm0,%xmm1
test   %rax,%rax
je     11d6 <lcfit2_compute_weights+0xb6>
lea    (%r12,%rcx,8),%rcx
xor    %edx,%edx
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
movsd  (%rcx,%rdx,8),%xmm1
maxsd  %xmm0,%xmm1
add    $0x1,%rdx
movapd %xmm1,%xmm0
cmp    %rdx,%rax
jne    11c0 <lcfit2_compute_weights+0xa0>
test   %r15,%r15
je     1221 <lcfit2_compute_weights+0x101>
xor    %ebx,%ebx
movsd  %xmm2,0x8(%rsp)
movsd  %xmm1,0x10(%rsp)
nopl   0x0(%rax)
movsd  (%r12,%rbx,8),%xmm0
subsd  %xmm1,%xmm0
call   1030 <exp@plt>
mov    %eax,%edi
movsd  0x8(%rsp),%xmm0
call   1040 <pow@plt>
movsd  0x10(%rsp),%xmm1
movsd  %xmm0,(%r14,%rbx,8)
add    $0x1,%rbx
cmp    %rbx,%r15
jne    11f0 <lcfit2_compute_weights+0xd0>
movapd %xmm1,%xmm0
add    $0x18,%rsp
pop    %rbx
pop    %r12
pop    %r14
pop    %r15
ret

LLVM IR (-O0):
@HUGE_VAL = global double 0.000000e+00
define double @lcfit2_compute_weights(i64 %0, double* %1, double %2, double* %3) {
%5 = alloca i64
%6 = alloca double*
%7 = alloca double
%8 = alloca double*
%9 = alloca double
%10 = alloca i64
%11 = alloca i64
store i64 %0, i64* %5
store double* %1, double** %6
store double %2, double* %7
store double* %3, double** %8
%12 = load double, double* @HUGE_VAL
%13 = fneg double %12
store double %13, double* %9
store i64 0, i64* %10
br label %14
14: ; preds = %31, %4
%15 = load i64, i64* %10
%16 = load i64, i64* %5
%17 = icmp ult i64 %15, %16
br i1 %17, label %18, label %34
18: ; preds = %14
%19 = load double*, double** %6
%20 = load i64, i64* %10
%21 = getelementptr inbounds double, double* %19, i64 %20
%22 = load double, double* %21
%23 = load double, double* %9
%24 = fcmp ogt double %22, %23
br i1 %24, label %25, label %30
25: ; preds = %18
%26 = load double*, double** %6
%27 = load i64, i64* %10
%28 = getelementptr inbounds double, double* %26, i64 %27
%29 = load double, double* %28
store double %29, double* %9
br label %30
30: ; preds = %25, %18
br label %31
31: ; preds = %30
%32 = load i64, i64* %10
%33 = add i64 %32, 1
store i64 %33, i64* %10
br label %14
34: ; preds = %14
store i64 0, i64* %11
br label %35
35: ; preds = %52, %34
%36 = load i64, i64* %11
%37 = load i64, i64* %5
%38 = icmp ult i64 %36, %37
br i1 %38, label %39, label %55
39: ; preds = %35
%40 = load double*, double** %6
%41 = load i64, i64* %11
%42 = getelementptr inbounds double, double* %40, i64 %41
%43 = load double, double* %42
%44 = load double, double* %9
%45 = fsub double %43, %44
%46 = call i32 @exp(double %45)
%47 = load double, double* %7
%48 = call double @pow(i32 %46, double %47)
%49 = load double*, double** %8
%50 = load i64, i64* %11
%51 = getelementptr inbounds double, double* %49, i64 %50
store double %48, double* %51
br label %52
52: ; preds = %39
%53 = load i64, i64* %11
%54 = add i64 %53, 1
store i64 %54, i64* %11
br label %35
55: ; preds = %35
%56 = load double, double* %9
ret double %56
}

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