#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 */
int /*<<< orphan*/ analysis(double *, int, double const *, int, int,
                            int const *);
int log10(double);
int /*<<< orphan*/ synthesis(double *, int, double *, int, int, int const *);
double coding_gain_1d(const double _r[2 * (1 << (2)) * 2 * (1 << (2))],
                      const int *_f) {
  int j;
  int i;
  double r[(1 << (2))];
  double rgt[2 * (1 << (2)) * (1 << (2))];
  double grgt[(1 << (2)) * (1 << (2))];
  double cg;
  for (j = 0; j < 2 * (1 << (2)); j++) {
    analysis(&rgt[j * (1 << (2))], 1, &_r[2 * (1 << (2)) * j], 1, 1, _f);
  }
  analysis(grgt, (1 << (2)), rgt, (1 << (2)), (1 << (2)), _f);
  for (j = 0; j < (1 << (2)); j++) {
    for (i = 0; i < (1 << (2)); i++) {
      r[i] = i == j ? 1 : 0;
    }
    synthesis(&rgt[j], (1 << (2)), r, 1, 1, _f);
  }
  cg = 0;
  for (j = 0; j < (1 << (2)); j++) {
    double h;
    h = 0;
    for (i = 0; i < 2 * (1 << (2)); i++) {
      h += rgt[i * (1 << (2)) + j] * rgt[i * (1 << (2)) + j];
    }
    cg -= 10 * log10(grgt[j * (1 << (2)) + j] * h);
  }
  return cg / (1 << (2));
}


x86 assembly (clang, -O0):
<coding_gain_1d>:
push   %rbp
mov    %rsp,%rbp
sub    $0x1d0,%rsp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
movl   $0x0,-0x14(%rbp)
cmpl   $0x8,-0x14(%rbp)
jge    11a0 <coding_gain_1d+0x70>
mov    -0x14(%rbp),%eax
shl    $0x2,%eax
cltq
lea    -0x140(%rbp),%rdi
shl    $0x3,%rax
add    %rax,%rdi
mov    -0x8(%rbp),%rdx
mov    -0x14(%rbp),%eax
shl    $0x3,%eax
cltq
shl    $0x3,%rax
add    %rax,%rdx
mov    -0x10(%rbp),%r9
mov    $0x1,%r8d
mov    %r8d,%esi
mov    %r8d,%ecx
call   1050 <analysis@plt>
mov    -0x14(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x14(%rbp)
jmp    114a <coding_gain_1d+0x1a>
lea    -0x1c0(%rbp),%rdi
lea    -0x140(%rbp),%rdx
mov    -0x10(%rbp),%r9
mov    $0x4,%r8d
mov    %r8d,%esi
mov    %r8d,%ecx
call   1050 <analysis@plt>
movl   $0x0,-0x14(%rbp)
cmpl   $0x4,-0x14(%rbp)
jge    124e <coding_gain_1d+0x11e>
movl   $0x0,-0x18(%rbp)
cmpl   $0x4,-0x18(%rbp)
jge    1213 <coding_gain_1d+0xe3>
mov    -0x18(%rbp),%edx
mov    -0x14(%rbp),%esi
xor    %eax,%eax
mov    $0x1,%ecx
cmp    %esi,%edx
cmove  %ecx,%eax
cvtsi2sd %eax,%xmm0
movslq -0x18(%rbp),%rax
movsd  %xmm0,-0x40(%rbp,%rax,8)
mov    -0x18(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x18(%rbp)
jmp    11db <coding_gain_1d+0xab>
movslq -0x14(%rbp),%rax
lea    -0x140(%rbp),%rdi
shl    $0x3,%rax
add    %rax,%rdi
lea    -0x40(%rbp),%rdx
mov    -0x10(%rbp),%r9
mov    $0x4,%esi
mov    $0x1,%r8d
mov    %r8d,%ecx
call   1040 <synthesis@plt>
mov    -0x14(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x14(%rbp)
jmp    11ca <coding_gain_1d+0x9a>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x1c8(%rbp)
ff
movl   $0x0,-0x14(%rbp)
cmpl   $0x4,-0x14(%rbp)
jge    131e <coding_gain_1d+0x1ee>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x1d0(%rbp)
ff
movl   $0x0,-0x18(%rbp)
cmpl   $0x8,-0x18(%rbp)
jge    12d4 <coding_gain_1d+0x1a4>
mov    -0x18(%rbp),%eax
shl    $0x2,%eax
add    -0x14(%rbp),%eax
cltq
movsd  -0x140(%rbp,%rax,8),%xmm0
ff ff
mov    -0x18(%rbp),%eax
shl    $0x2,%eax
add    -0x14(%rbp),%eax
cltq
movsd  -0x140(%rbp,%rax,8),%xmm2
ff ff
movsd  -0x1d0(%rbp),%xmm1
ff
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x1d0(%rbp)
ff
mov    -0x18(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x18(%rbp)
jmp    127c <coding_gain_1d+0x14c>
mov    -0x14(%rbp),%eax
shl    $0x2,%eax
add    -0x14(%rbp),%eax
cltq
movsd  -0x1c0(%rbp,%rax,8),%xmm0
ff ff
mulsd  -0x1d0(%rbp),%xmm0
ff
call   1030 <log10@plt>
imul   $0xa,%eax,%eax
cvtsi2sd %eax,%xmm1
movsd  -0x1c8(%rbp),%xmm0
ff
subsd  %xmm1,%xmm0
movsd  %xmm0,-0x1c8(%rbp)
ff
mov    -0x14(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x14(%rbp)
jmp    1260 <coding_gain_1d+0x130>
movsd  -0x1c8(%rbp),%xmm0
ff
movsd  0xcd2(%rip),%xmm1
divsd  %xmm1,%xmm0
add    $0x1d0,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<coding_gain_1d>:
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x1b0,%rsp
mov    %rsi,%r14
mov    %rdi,%rbx
lea    0x30(%rsp),%r15
mov    %r15,%rdi
mov    $0x1,%esi
mov    %rbx,%rdx
mov    $0x1,%ecx
mov    $0x1,%r8d
mov    %r14,%r9
call   1050 <analysis@plt>
lea    0x50(%rsp),%rdi
lea    0x40(%rbx),%rdx
mov    $0x1,%esi
mov    $0x1,%ecx
mov    $0x1,%r8d
mov    %r14,%r9
call   1050 <analysis@plt>
lea    0x70(%rsp),%rdi
lea    0x80(%rbx),%rdx
mov    $0x1,%esi
mov    $0x1,%ecx
mov    $0x1,%r8d
mov    %r14,%r9
call   1050 <analysis@plt>
lea    0x90(%rsp),%rdi
lea    0xc0(%rbx),%rdx
mov    $0x1,%esi
mov    $0x1,%ecx
mov    $0x1,%r8d
mov    %r14,%r9
call   1050 <analysis@plt>
lea    0xb0(%rsp),%rdi
lea    0x100(%rbx),%rdx
mov    $0x1,%esi
mov    $0x1,%ecx
mov    $0x1,%r8d
mov    %r14,%r9
call   1050 <analysis@plt>
lea    0xd0(%rsp),%rdi
lea    0x140(%rbx),%rdx
mov    $0x1,%esi
mov    $0x1,%ecx
mov    $0x1,%r8d
mov    %r14,%r9
call   1050 <analysis@plt>
lea    0xf0(%rsp),%rdi
lea    0x180(%rbx),%rdx
mov    $0x1,%esi
mov    $0x1,%ecx
mov    $0x1,%r8d
mov    %r14,%r9
call   1050 <analysis@plt>
lea    0x110(%rsp),%rdi
add    $0x1c0,%rbx
mov    $0x1,%esi
mov    %rbx,%rdx
mov    $0x1,%ecx
mov    $0x1,%r8d
mov    %r14,%r9
call   1050 <analysis@plt>
lea    0x130(%rsp),%r13
mov    %r13,%rdi
mov    $0x4,%esi
mov    %r15,%rdx
mov    $0x4,%ecx
mov    $0x4,%r8d
mov    %r14,%r9
call   1050 <analysis@plt>
movabs $0x3ff0000000000000,%rbx
mov    %rbx,0x10(%rsp)
xorpd  %xmm0,%xmm0
movupd %xmm0,0x18(%rsp)
movq   $0x0,0x28(%rsp)
lea    0x10(%rsp),%r12
mov    %r15,%rdi
mov    $0x4,%esi
mov    %r12,%rdx
mov    $0x1,%ecx
mov    $0x1,%r8d
mov    %r14,%r9
call   1040 <synthesis@plt>
xorpd  %xmm0,%xmm0
movhpd 0xd15(%rip),%xmm0
movapd %xmm0,0x10(%rsp)
lea    0x38(%rsp),%rdi
xorpd  %xmm0,%xmm0
movapd %xmm0,0x20(%rsp)
mov    $0x4,%esi
mov    %r12,%rdx
mov    $0x1,%ecx
mov    $0x1,%r8d
mov    %r14,%r9
call   1040 <synthesis@plt>
xorpd  %xmm0,%xmm0
movapd %xmm0,0x10(%rsp)
mov    %rbx,0x20(%rsp)
movq   $0x0,0x28(%rsp)
lea    0x40(%rsp),%rdi
mov    $0x4,%esi
mov    %r12,%rdx
mov    $0x1,%ecx
mov    $0x1,%r8d
mov    %r14,%r9
call   1040 <synthesis@plt>
xorpd  %xmm0,%xmm0
movapd %xmm0,0x10(%rsp)
movq   $0x0,0x20(%rsp)
mov    %rbx,0x28(%rsp)
lea    0x48(%rsp),%rdi
mov    $0x4,%esi
mov    %r12,%rdx
mov    $0x1,%ecx
mov    $0x1,%r8d
mov    %r14,%r9
call   1040 <synthesis@plt>
mov    $0xfffffffffffffffc,%rbx
xorpd  %xmm0,%xmm0
cs nopw 0x0(%rax,%rax,1)
movsd  %xmm0,0x8(%rsp)
movsd  0x50(%rsp,%rbx,8),%xmm0
movsd  0x70(%rsp,%rbx,8),%xmm1
mulsd  %xmm0,%xmm0
addsd  0xc52(%rip),%xmm0
mulsd  %xmm1,%xmm1
addsd  %xmm0,%xmm1
movsd  0x90(%rsp,%rbx,8),%xmm0
mulsd  %xmm0,%xmm0
addsd  %xmm1,%xmm0
movsd  0xb0(%rsp,%rbx,8),%xmm1
mulsd  %xmm1,%xmm1
addsd  %xmm0,%xmm1
movsd  0xd0(%rsp,%rbx,8),%xmm0
mulsd  %xmm0,%xmm0
addsd  %xmm1,%xmm0
movsd  0xf0(%rsp,%rbx,8),%xmm1
mulsd  %xmm1,%xmm1
addsd  %xmm0,%xmm1
movsd  0x110(%rsp,%rbx,8),%xmm2
mulsd  %xmm2,%xmm2
addsd  %xmm1,%xmm2
movsd  0x130(%rsp,%rbx,8),%xmm0
mulsd  %xmm0,%xmm0
addsd  %xmm2,%xmm0
mulsd  0x0(%r13),%xmm0
call   1030 <log10@plt>
add    %eax,%eax
lea    (%rax,%rax,4),%eax
xorps  %xmm0,%xmm0
cvtsi2sd %eax,%xmm0
movsd  0x8(%rsp),%xmm1
subsd  %xmm0,%xmm1
movsd  %xmm1,0x8(%rsp)
movsd  0x8(%rsp),%xmm0
add    $0x28,%r13
inc    %rbx
jne    13a0 <coding_gain_1d+0x270>
mulsd  0xb9a(%rip),%xmm0
add    $0x1b0,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
ret

LLVM IR (-O0):
define double @coding_gain_1d(double* %0, i32* %1) {
%3 = alloca double*
%4 = alloca i32*
%5 = alloca i32
%6 = alloca i32
%7 = alloca [4 x double]
%8 = alloca [32 x double]
%9 = alloca [16 x double]
%10 = alloca double
%11 = alloca double
store double* %0, double** %3
store i32* %1, i32** %4
store i32 0, i32* %5
br label %12
12: ; preds = %27, %2
%13 = load i32, i32* %5
%14 = icmp slt i32 %13, 8
br i1 %14, label %15, label %30
15: ; preds = %12
%16 = load i32, i32* %5
%17 = mul nsw i32 %16, 4
%18 = sext i32 %17 to i64
%19 = getelementptr inbounds [32 x double], [32 x double]* %8, i64 0, i64 %18
%20 = load double*, double** %3
%21 = load i32, i32* %5
%22 = mul nsw i32 8, %21
%23 = sext i32 %22 to i64
%24 = getelementptr inbounds double, double* %20, i64 %23
%25 = load i32*, i32** %4
%26 = call i32 @analysis(double* %19, i32 1, double* %24, i32 1, i32 1, i32* %25)
br label %27
27: ; preds = %15
%28 = load i32, i32* %5
%29 = add nsw i32 %28, 1
store i32 %29, i32* %5
br label %12
30: ; preds = %12
%31 = getelementptr inbounds [16 x double], [16 x double]* %9, i64 0, i64 0
%32 = getelementptr inbounds [32 x double], [32 x double]* %8, i64 0, i64 0
%33 = load i32*, i32** %4
%34 = call i32 @analysis(double* %31, i32 4, double* %32, i32 4, i32 4, i32* %33)
store i32 0, i32* %5
br label %35
35: ; preds = %62, %30
%36 = load i32, i32* %5
%37 = icmp slt i32 %36, 4
br i1 %37, label %38, label %65
38: ; preds = %35
store i32 0, i32* %6
br label %39
39: ; preds = %52, %38
%40 = load i32, i32* %6
%41 = icmp slt i32 %40, 4
br i1 %41, label %42, label %55
42: ; preds = %39
%43 = load i32, i32* %6
%44 = load i32, i32* %5
%45 = icmp eq i32 %43, %44
%46 = zext i1 %45 to i64
%47 = select i1 %45, i32 1, i32 0
%48 = sitofp i32 %47 to double
%49 = load i32, i32* %6
%50 = sext i32 %49 to i64
%51 = getelementptr inbounds [4 x double], [4 x double]* %7, i64 0, i64 %50
store double %48, double* %51
br label %52
52: ; preds = %42
%53 = load i32, i32* %6
%54 = add nsw i32 %53, 1
store i32 %54, i32* %6
br label %39
55: ; preds = %39
%56 = load i32, i32* %5
%57 = sext i32 %56 to i64
%58 = getelementptr inbounds [32 x double], [32 x double]* %8, i64 0, i64 %57
%59 = getelementptr inbounds [4 x double], [4 x double]* %7, i64 0, i64 0
%60 = load i32*, i32** %4
%61 = call i32 @synthesis(double* %58, i32 4, double* %59, i32 1, i32 1, i32* %60)
br label %62
62: ; preds = %55
%63 = load i32, i32* %5
%64 = add nsw i32 %63, 1
store i32 %64, i32* %5
br label %35
65: ; preds = %35
store double 0.000000e+00, double* %10
store i32 0, i32* %5
br label %66
66: ; preds = %108, %65
%67 = load i32, i32* %5
%68 = icmp slt i32 %67, 4
br i1 %68, label %69, label %111
69: ; preds = %66
store double 0.000000e+00, double* %11
store i32 0, i32* %6
br label %70
70: ; preds = %90, %69
%71 = load i32, i32* %6
%72 = icmp slt i32 %71, 8
br i1 %72, label %73, label %93
73: ; preds = %70
%74 = load i32, i32* %6
%75 = mul nsw i32 %74, 4
%76 = load i32, i32* %5
%77 = add nsw i32 %75, %76
%78 = sext i32 %77 to i64
%79 = getelementptr inbounds [32 x double], [32 x double]* %8, i64 0, i64 %78
%80 = load double, double* %79
%81 = load i32, i32* %6
%82 = mul nsw i32 %81, 4
%83 = load i32, i32* %5
%84 = add nsw i32 %82, %83
%85 = sext i32 %84 to i64
%86 = getelementptr inbounds [32 x double], [32 x double]* %8, i64 0, i64 %85
%87 = load double, double* %86
%88 = load double, double* %11
%89 = call double @llvm.fmuladd.f64(double %80, double %87, double %88)
store double %89, double* %11
br label %90
90: ; preds = %73
%91 = load i32, i32* %6
%92 = add nsw i32 %91, 1
store i32 %92, i32* %6
br label %70
93: ; preds = %70
%94 = load i32, i32* %5
%95 = mul nsw i32 %94, 4
%96 = load i32, i32* %5
%97 = add nsw i32 %95, %96
%98 = sext i32 %97 to i64
%99 = getelementptr inbounds [16 x double], [16 x double]* %9, i64 0, i64 %98
%100 = load double, double* %99
%101 = load double, double* %11
%102 = fmul double %100, %101
%103 = call i32 @log10(double %102)
%104 = mul nsw i32 10, %103
%105 = sitofp i32 %104 to double
%106 = load double, double* %10
%107 = fsub double %106, %105
store double %107, double* %10
br label %108
108: ; preds = %93
%109 = load i32, i32* %5
%110 = add nsw i32 %109, 1
store i32 %110, i32* %5
br label %66
111: ; preds = %66
%112 = load double, double* %10
%113 = fdiv double %112, 4.000000e+00
ret double %113
}

LLVM IR (-O2):
define double @coding_gain_1d(double* %0, i32* %1) {
%3 = alloca [4 x double]
%4 = alloca [32 x double]
%5 = alloca [16 x double]
%6 = bitcast [4 x double]* %3 to i8*
%7 = bitcast [32 x double]* %4 to i8*
%8 = bitcast [16 x double]* %5 to i8*
%9 = getelementptr inbounds [32 x double], [32 x double]* %4, i64 0, i64 0
%10 = call i32 @analysis(double* %9, i32 1, double* %0, i32 1, i32 1, i32* %1)
%11 = getelementptr inbounds [32 x double], [32 x double]* %4, i64 0, i64 4
%12 = getelementptr inbounds double, double* %0, i64 8
%13 = call i32 @analysis(double* %11, i32 1, double* %12, i32 1, i32 1, i32* %1)
%14 = getelementptr inbounds [32 x double], [32 x double]* %4, i64 0, i64 8
%15 = getelementptr inbounds double, double* %0, i64 16
%16 = call i32 @analysis(double* %14, i32 1, double* %15, i32 1, i32 1, i32* %1)
%17 = getelementptr inbounds [32 x double], [32 x double]* %4, i64 0, i64 12
%18 = getelementptr inbounds double, double* %0, i64 24
%19 = call i32 @analysis(double* %17, i32 1, double* %18, i32 1, i32 1, i32* %1)
%20 = getelementptr inbounds [32 x double], [32 x double]* %4, i64 0, i64 16
%21 = getelementptr inbounds double, double* %0, i64 32
%22 = call i32 @analysis(double* %20, i32 1, double* %21, i32 1, i32 1, i32* %1)
%23 = getelementptr inbounds [32 x double], [32 x double]* %4, i64 0, i64 20
%24 = getelementptr inbounds double, double* %0, i64 40
%25 = call i32 @analysis(double* %23, i32 1, double* %24, i32 1, i32 1, i32* %1)
%26 = getelementptr inbounds [32 x double], [32 x double]* %4, i64 0, i64 24
%27 = getelementptr inbounds double, double* %0, i64 48
%28 = call i32 @analysis(double* %26, i32 1, double* %27, i32 1, i32 1, i32* %1)
%29 = getelementptr inbounds [32 x double], [32 x double]* %4, i64 0, i64 28
%30 = getelementptr inbounds double, double* %0, i64 56
%31 = call i32 @analysis(double* %29, i32 1, double* %30, i32 1, i32 1, i32* %1)
%32 = getelementptr inbounds [16 x double], [16 x double]* %5, i64 0, i64 0
%33 = call i32 @analysis(double* %32, i32 4, double* %9, i32 4, i32 4, i32* %1)
%34 = getelementptr inbounds [4 x double], [4 x double]* %3, i64 0, i64 0
store double 1.000000e+00, double* %34
%35 = getelementptr inbounds [4 x double], [4 x double]* %3, i64 0, i64 1
%36 = getelementptr inbounds [4 x double], [4 x double]* %3, i64 0, i64 2
%37 = getelementptr inbounds [4 x double], [4 x double]* %3, i64 0, i64 3
%38 = bitcast double* %35 to i8*
%39 = call i32 @synthesis(double* %9, i32 4, double* %34, i32 1, i32 1, i32* %1)
%40 = bitcast [4 x double]* %3 to <2 x double>*
store <2 x double> <double 0.000000e+00, double 1.000000e+00>, <2 x double>* %40
%41 = getelementptr inbounds [32 x double], [32 x double]* %4, i64 0, i64 1
%42 = bitcast double* %36 to i8*
%43 = call i32 @synthesis(double* %41, i32 4, double* %34, i32 1, i32 1, i32* %1)
%44 = bitcast [4 x double]* %3 to i8*
store double 1.000000e+00, double* %36
store double 0.000000e+00, double* %37
%45 = getelementptr inbounds [32 x double], [32 x double]* %4, i64 0, i64 2
%46 = call i32 @synthesis(double* %45, i32 4, double* %34, i32 1, i32 1, i32* %1)
%47 = bitcast [4 x double]* %3 to i8*
store double 1.000000e+00, double* %37
%48 = getelementptr inbounds [32 x double], [32 x double]* %4, i64 0, i64 3
%49 = call i32 @synthesis(double* %48, i32 4, double* %34, i32 1, i32 1, i32* %1)
br label %50
50: ; preds = %2, %50
%51 = phi i64 [ 0, %2 ], [ %92, %50 ]
%52 = phi double [ 0.000000e+00, %2 ], [ %91, %50 ]
%53 = getelementptr inbounds [32 x double], [32 x double]* %4, i64 0, i64 %51
%54 = load double, double* %53
%55 = call double @llvm.fmuladd.f64(double %54, double %54, double 0.000000e+00)
%56 = add nuw nsw i64 %51, 4
%57 = getelementptr inbounds [32 x double], [32 x double]* %4, i64 0, i64 %56
%58 = load double, double* %57
%59 = call double @llvm.fmuladd.f64(double %58, double %58, double %55)
%60 = add nuw nsw i64 %51, 8
%61 = getelementptr inbounds [32 x double], [32 x double]* %4, i64 0, i64 %60
%62 = load double, double* %61
%63 = call double @llvm.fmuladd.f64(double %62, double %62, double %59)
%64 = add nuw nsw i64 %51, 12
%65 = getelementptr inbounds [32 x double], [32 x double]* %4, i64 0, i64 %64
%66 = load double, double* %65
%67 = call double @llvm.fmuladd.f64(double %66, double %66, double %63)
%68 = add nuw nsw i64 %51, 16
%69 = getelementptr inbounds [32 x double], [32 x double]* %4, i64 0, i64 %68
%70 = load double, double* %69
%71 = call double @llvm.fmuladd.f64(double %70, double %70, double %67)
%72 = add nuw nsw i64 %51, 20
%73 = getelementptr inbounds [32 x double], [32 x double]* %4, i64 0, i64 %72
%74 = load double, double* %73
%75 = call double @llvm.fmuladd.f64(double %74, double %74, double %71)
%76 = add nuw nsw i64 %51, 24
%77 = getelementptr inbounds [32 x double], [32 x double]* %4, i64 0, i64 %76
%78 = load double, double* %77
%79 = call double @llvm.fmuladd.f64(double %78, double %78, double %75)
%80 = add nuw nsw i64 %51, 28
%81 = getelementptr inbounds [32 x double], [32 x double]* %4, i64 0, i64 %80
%82 = load double, double* %81
%83 = call double @llvm.fmuladd.f64(double %82, double %82, double %79)
%84 = mul nuw nsw i64 %51, 5
%85 = getelementptr inbounds [16 x double], [16 x double]* %5, i64 0, i64 %84
%86 = load double, double* %85
%87 = fmul double %83, %86
%88 = call i32 @log10(double %87)
%89 = mul nsw i32 %88, 10
%90 = sitofp i32 %89 to double
%91 = fsub double %52, %90
%92 = add nuw nsw i64 %51, 1
%93 = icmp eq i64 %92, 4
br i1 %93, label %94, label %50
94: ; preds = %50
%95 = fmul double %91, 2.500000e-01
ret double %95
}