#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 log10(double);
double calculate_rms(double *samples[7], int ch, int n) {
  double rms_all, rms_left, rms_right;
  int i;
  if (ch == 1) {
    rms_all = 0;
    for (i = 0; i < n; i++)
      rms_all += (samples[0][i] * samples[0][i]);
    rms_all /= n;
  } else {
    rms_left = rms_right = 0;
    for (i = 0; i < n; i++) {
      rms_left += (samples[0][i] * samples[0][i]);
      rms_right += (samples[1][i] * samples[1][i]);
    }
    rms_all = (rms_left + rms_right) / (2.0 * n);
  }
  rms_all = 10.0 * log10(rms_all + 1e-10);
  return rms_all;
}


x86 assembly (clang, -O0):
<calculate_rms>:
push   %rbp
mov    %rsp,%rbp
sub    $0x30,%rsp
mov    %rdi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
mov    %edx,-0x10(%rbp)
cmpl   $0x1,-0xc(%rbp)
jne    119f <calculate_rms+0x8f>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x18(%rbp)
movl   $0x0,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0x10(%rbp),%eax
jge    1187 <calculate_rms+0x77>
mov    -0x8(%rbp),%rax
mov    (%rax),%rax
movslq -0x2c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x8(%rbp),%rax
mov    (%rax),%rax
movslq -0x2c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm2
movsd  -0x18(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x18(%rbp)
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    113b <calculate_rms+0x2b>
cvtsi2sdl -0x10(%rbp),%xmm1
movsd  -0x18(%rbp),%xmm0
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x18(%rbp)
jmp    125a <calculate_rms+0x14a>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x28(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x20(%rbp)
movl   $0x0,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0x10(%rbp),%eax
jge    1236 <calculate_rms+0x126>
mov    -0x8(%rbp),%rax
mov    (%rax),%rax
movslq -0x2c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x8(%rbp),%rax
mov    (%rax),%rax
movslq -0x2c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm2
movsd  -0x20(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x20(%rbp)
mov    -0x8(%rbp),%rax
mov    0x8(%rax),%rax
movslq -0x2c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x8(%rbp),%rax
mov    0x8(%rax),%rax
movslq -0x2c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm2
movsd  -0x28(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x28(%rbp)
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    11b6 <calculate_rms+0xa6>
movsd  -0x20(%rbp),%xmm0
addsd  -0x28(%rbp),%xmm0
cvtsi2sdl -0x10(%rbp),%xmm2
movsd  0xdb3(%rip),%xmm1
mulsd  %xmm2,%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x18(%rbp)
movsd  0xdae(%rip),%xmm0
addsd  -0x18(%rbp),%xmm0
call   1030 <log10@plt>
movaps %xmm0,%xmm1
movsd  0xd91(%rip),%xmm0
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0x18(%rbp)
movsd  -0x18(%rbp),%xmm0
add    $0x30,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<calculate_rms>:
push   %rax
cmp    $0x1,%esi
jne    1139 <calculate_rms+0x29>
test   %edx,%edx
jle    1166 <calculate_rms+0x56>
mov    (%rdi),%rcx
mov    %edx,%esi
lea    -0x1(%rsi),%rdi
mov    %esi,%eax
and    $0x3,%eax
cmp    $0x3,%rdi
jae    117c <calculate_rms+0x6c>
xorpd  %xmm0,%xmm0
xor    %edi,%edi
jmp    11d0 <calculate_rms+0xc0>
test   %edx,%edx
jle    1173 <calculate_rms+0x63>
mov    (%rdi),%rcx
mov    0x8(%rdi),%rsi
mov    %edx,%eax
lea    -0x1(%rax),%rdi
mov    %eax,%r8d
and    $0x3,%r8d
cmp    $0x3,%rdi
jae    1202 <calculate_rms+0xf2>
xorpd  %xmm1,%xmm1
xor    %edi,%edi
jmp    126b <calculate_rms+0x15b>
xorpd  %xmm0,%xmm0
cvtsi2sd %edx,%xmm1
jmp    12b2 <calculate_rms+0x1a2>
xorpd  %xmm0,%xmm0
jmp    12a7 <calculate_rms+0x197>
and    $0xfffffffc,%esi
xorpd  %xmm0,%xmm0
xor    %edi,%edi
cs nopw 0x0(%rax,%rax,1)
nop
movsd  (%rcx,%rdi,8),%xmm1
movsd  0x8(%rcx,%rdi,8),%xmm2
mulsd  %xmm1,%xmm1
addsd  %xmm0,%xmm1
mulsd  %xmm2,%xmm2
addsd  %xmm1,%xmm2
movsd  0x10(%rcx,%rdi,8),%xmm1
mulsd  %xmm1,%xmm1
addsd  %xmm2,%xmm1
movsd  0x18(%rcx,%rdi,8),%xmm0
mulsd  %xmm0,%xmm0
addsd  %xmm1,%xmm0
add    $0x4,%rdi
cmp    %rdi,%rsi
jne    1190 <calculate_rms+0x80>
test   %rax,%rax
je     11f6 <calculate_rms+0xe6>
lea    (%rcx,%rdi,8),%rcx
xor    %esi,%esi
nopl   0x0(%rax,%rax,1)
movsd  (%rcx,%rsi,8),%xmm1
mulsd  %xmm1,%xmm1
addsd  %xmm1,%xmm0
add    $0x1,%rsi
cmp    %rsi,%rax
jne    11e0 <calculate_rms+0xd0>
xorps  %xmm1,%xmm1
cvtsi2sd %edx,%xmm1
jmp    12b2 <calculate_rms+0x1a2>
and    $0xfffffffc,%eax
xorpd  %xmm0,%xmm0
xor    %edi,%edi
nopl   0x0(%rax,%rax,1)
movsd  (%rsi,%rdi,8),%xmm2
movsd  0x8(%rsi,%rdi,8),%xmm3
movsd  0x10(%rsi,%rdi,8),%xmm4
movsd  0x18(%rsi,%rdi,8),%xmm1
movhpd (%rcx,%rdi,8),%xmm2
mulpd  %xmm2,%xmm2
addpd  %xmm0,%xmm2
movhpd 0x8(%rcx,%rdi,8),%xmm3
mulpd  %xmm3,%xmm3
addpd  %xmm2,%xmm3
movhpd 0x10(%rcx,%rdi,8),%xmm4
mulpd  %xmm4,%xmm4
addpd  %xmm3,%xmm4
movhpd 0x18(%rcx,%rdi,8),%xmm1
mulpd  %xmm1,%xmm1
addpd  %xmm4,%xmm1
add    $0x4,%rdi
movapd %xmm1,%xmm0
cmp    %rdi,%rax
jne    1210 <calculate_rms+0x100>
test   %r8,%r8
je     129b <calculate_rms+0x18b>
lea    (%rsi,%rdi,8),%rax
lea    (%rcx,%rdi,8),%rcx
xor    %esi,%esi
nopw   0x0(%rax,%rax,1)
movsd  (%rax,%rsi,8),%xmm0
movhpd (%rcx,%rsi,8),%xmm0
mulpd  %xmm0,%xmm0
addpd  %xmm0,%xmm1
add    $0x1,%rsi
cmp    %rsi,%r8
jne    1280 <calculate_rms+0x170>
movapd %xmm1,%xmm0
unpckhpd %xmm1,%xmm0
addsd  %xmm1,%xmm0
xorps  %xmm1,%xmm1
cvtsi2sd %edx,%xmm1
addsd  %xmm1,%xmm1
divsd  %xmm1,%xmm0
addsd  0xd42(%rip),%xmm0
call   1030 <log10@plt>
mulsd  0xd3d(%rip),%xmm0
pop    %rax
ret

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

LLVM IR (-O2):
define double @calculate_rms(double** nocapture %0, i32 %1, i32 %2) {
%4 = icmp eq i32 %1, 1
%5 = icmp sgt i32 %2, 0
br i1 %4, label %17, label %6
6: ; preds = %3
br i1 %5, label %7, label %129
7: ; preds = %6
%8 = load double*, double** %0
%9 = getelementptr inbounds double*, double** %0, i64 1
%10 = load double*, double** %9
%11 = zext i32 %2 to i64
%12 = add nsw i64 %11, -1
%13 = and i64 %11, 3
%14 = icmp ult i64 %12, 3
br i1 %14, label %105, label %15
15: ; preds = %7
%16 = and i64 %11, 4294967292
br label %67
17: ; preds = %3
br i1 %5, label %18, label %63
18: ; preds = %17
%19 = load double*, double** %0
%20 = zext i32 %2 to i64
%21 = add nsw i64 %20, -1
%22 = and i64 %20, 3
%23 = icmp ult i64 %21, 3
br i1 %23, label %48, label %24
24: ; preds = %18
%25 = and i64 %20, 4294967292
br label %26
26: ; preds = %26, %24
%27 = phi i64 [ 0, %24 ], [ %45, %26 ]
%28 = phi double [ 0.000000e+00, %24 ], [ %44, %26 ]
%29 = phi i64 [ 0, %24 ], [ %46, %26 ]
%30 = getelementptr inbounds double, double* %19, i64 %27
%31 = load double, double* %30
%32 = tail call double @llvm.fmuladd.f64(double %31, double %31, double %28)
%33 = or i64 %27, 1
%34 = getelementptr inbounds double, double* %19, i64 %33
%35 = load double, double* %34
%36 = tail call double @llvm.fmuladd.f64(double %35, double %35, double %32)
%37 = or i64 %27, 2
%38 = getelementptr inbounds double, double* %19, i64 %37
%39 = load double, double* %38
%40 = tail call double @llvm.fmuladd.f64(double %39, double %39, double %36)
%41 = or i64 %27, 3
%42 = getelementptr inbounds double, double* %19, i64 %41
%43 = load double, double* %42
%44 = tail call double @llvm.fmuladd.f64(double %43, double %43, double %40)
%45 = add nuw nsw i64 %27, 4
%46 = add i64 %29, 4
%47 = icmp eq i64 %46, %25
br i1 %47, label %48, label %26
48: ; preds = %26, %18
%49 = phi double [ undef, %18 ], [ %44, %26 ]
%50 = phi i64 [ 0, %18 ], [ %45, %26 ]
%51 = phi double [ 0.000000e+00, %18 ], [ %44, %26 ]
%52 = icmp eq i64 %22, 0
br i1 %52, label %63, label %53
53: ; preds = %48, %53
%54 = phi i64 [ %60, %53 ], [ %50, %48 ]
%55 = phi double [ %59, %53 ], [ %51, %48 ]
%56 = phi i64 [ %61, %53 ], [ 0, %48 ]
%57 = getelementptr inbounds double, double* %19, i64 %54
%58 = load double, double* %57
%59 = tail call double @llvm.fmuladd.f64(double %58, double %58, double %55)
%60 = add nuw nsw i64 %54, 1
%61 = add i64 %56, 1
%62 = icmp eq i64 %61, %22
br i1 %62, label %63, label %53
63: ; preds = %48, %53, %17
%64 = phi double [ 0.000000e+00, %17 ], [ %49, %48 ], [ %59, %53 ]
%65 = sitofp i32 %2 to double
%66 = fdiv double %64, %65
br label %134
67: ; preds = %67, %15
%68 = phi i64 [ 0, %15 ], [ %102, %67 ]
%69 = phi <2 x double> [ zeroinitializer, %15 ], [ %101, %67 ]
%70 = phi i64 [ 0, %15 ], [ %103, %67 ]
%71 = getelementptr inbounds double, double* %8, i64 %68
%72 = load double, double* %71
%73 = getelementptr inbounds double, double* %10, i64 %68
%74 = load double, double* %73
%75 = insertelement <2 x double> poison, double %74, i64 0
%76 = insertelement <2 x double> %75, double %72, i64 1
%77 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %76, <2 x double> %76, <2 x double> %69)
%78 = or i64 %68, 1
%79 = getelementptr inbounds double, double* %8, i64 %78
%80 = load double, double* %79
%81 = getelementptr inbounds double, double* %10, i64 %78
%82 = load double, double* %81
%83 = insertelement <2 x double> poison, double %82, i64 0
%84 = insertelement <2 x double> %83, double %80, i64 1
%85 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %84, <2 x double> %84, <2 x double> %77)
%86 = or i64 %68, 2
%87 = getelementptr inbounds double, double* %8, i64 %86
%88 = load double, double* %87
%89 = getelementptr inbounds double, double* %10, i64 %86
%90 = load double, double* %89
%91 = insertelement <2 x double> poison, double %90, i64 0
%92 = insertelement <2 x double> %91, double %88, i64 1
%93 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %92, <2 x double> %92, <2 x double> %85)
%94 = or i64 %68, 3
%95 = getelementptr inbounds double, double* %8, i64 %94
%96 = load double, double* %95
%97 = getelementptr inbounds double, double* %10, i64 %94
%98 = load double, double* %97
%99 = insertelement <2 x double> poison, double %98, i64 0
%100 = insertelement <2 x double> %99, double %96, i64 1
%101 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %100, <2 x double> %100, <2 x double> %93)
%102 = add nuw nsw i64 %68, 4
%103 = add i64 %70, 4
%104 = icmp eq i64 %103, %16
br i1 %104, label %105, label %67
105: ; preds = %67, %7
%106 = phi <2 x double> [ undef, %7 ], [ %101, %67 ]
%107 = phi i64 [ 0, %7 ], [ %102, %67 ]
%108 = phi <2 x double> [ zeroinitializer, %7 ], [ %101, %67 ]
%109 = icmp eq i64 %13, 0
br i1 %109, label %124, label %110
110: ; preds = %105, %110
%111 = phi i64 [ %121, %110 ], [ %107, %105 ]
%112 = phi <2 x double> [ %120, %110 ], [ %108, %105 ]
%113 = phi i64 [ %122, %110 ], [ 0, %105 ]
%114 = getelementptr inbounds double, double* %8, i64 %111
%115 = load double, double* %114
%116 = getelementptr inbounds double, double* %10, i64 %111
%117 = load double, double* %116
%118 = insertelement <2 x double> poison, double %117, i64 0
%119 = insertelement <2 x double> %118, double %115, i64 1
%120 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %119, <2 x double> %119, <2 x double> %112)
%121 = add nuw nsw i64 %111, 1
%122 = add i64 %113, 1
%123 = icmp eq i64 %122, %13
br i1 %123, label %124, label %110
124: ; preds = %110, %105
%125 = phi <2 x double> [ %106, %105 ], [ %120, %110 ]
%126 = shufflevector <2 x double> %125, <2 x double> poison, <2 x i32> <i32 1, i32 undef>
%127 = fadd <2 x double> %126, %125
%128 = extractelement <2 x double> %127, i64 0
br label %129
129: ; preds = %124, %6
%130 = phi double [ %128, %124 ], [ 0.000000e+00, %6 ]
%131 = sitofp i32 %2 to double
%132 = fmul double %131, 2.000000e+00
%133 = fdiv double %130, %132
br label %134
134: ; preds = %129, %63
%135 = phi double [ %66, %63 ], [ %133, %129 ]
%136 = fadd double %135, 1.000000e-10
%137 = tail call double @log10(double %136)
%138 = fmul double %137, 1.000000e+01
ret double %138
}