#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__ check_malloc(int);
float gaussian_rand(float *, float *, long *);
int /*<<< orphan*/ printf(char *, double, ...);
int /*<<< orphan*/ scanf(char *, int *);
double sqrt(float);
main(int ac, char **av) {
  float *real, *imag, amp;
  float rmean, imean, rsig, isig;
  int i, n;
  long seed = 0;
  float fsigma = 1.0;
  float fmean = 0.0;
  scanf("%d", &n);
  real = (float *)check_malloc(n * sizeof(float));
  imag = (float *)check_malloc(n * sizeof(float));
  amp = 0.0;
  for (i = 0; i < n; i++) {
    real[i] = gaussian_rand(&fsigma, &fmean, &seed);
    imag[i] = gaussian_rand(&fsigma, &fmean, &seed);
    rmean = rmean + real[i];
    imean = imean + imag[i];
    amp = amp + (real[i] * real[i] + imag[i] * imag[i]);
  }
  rmean = rmean / (float)(n);
  imean = imean / (float)(n);
  rsig = 0.0;
  isig = 0.0;
  for (i = 0; i < n; i++) {
    rsig = rsig + (real[i] - rmean) * (real[i] - rmean);
    isig = isig + (imag[i] - imean) * (imag[i] - imean);
  }
  rsig = sqrt(rsig / (float)(n - 1));
  isig = sqrt(isig / (float)(n - 1));
  printf("rmean= %f, rsig= %f\n", rmean, rsig);
  printf("imean= %f, isig= %f\n", imean, isig);
  printf("avg. amp= %f\n", sqrt(amp / (float)(n)));
}


x86 assembly (clang, -O0):
<main>:
push   %rbp
mov    %rsp,%rbp
sub    $0x50,%rsp
movl   $0x0,-0x4(%rbp)
mov    %edi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
movq   $0x0,-0x48(%rbp)
movss  0xe8a(%rip),%xmm0
movss  %xmm0,-0x4c(%rbp)
xorps  %xmm0,%xmm0
movss  %xmm0,-0x50(%rbp)
lea    0xe7a(%rip),%rdi
lea    -0x3c(%rbp),%rsi
call   1060 <scanf@plt>
movslq -0x3c(%rbp),%rax
shl    $0x2,%rax
mov    %eax,%edi
call   1050 <check_malloc@plt>
mov    %rax,-0x18(%rbp)
movslq -0x3c(%rbp),%rax
shl    $0x2,%rax
mov    %eax,%edi
call   1050 <check_malloc@plt>
mov    %rax,-0x20(%rbp)
xorps  %xmm0,%xmm0
movss  %xmm0,-0x24(%rbp)
movl   $0x0,-0x38(%rbp)
mov    -0x38(%rbp),%eax
cmp    -0x3c(%rbp),%eax
jge    1296 <main+0x146>
lea    -0x4c(%rbp),%rdi
lea    -0x50(%rbp),%rsi
lea    -0x48(%rbp),%rdx
call   1040 <gaussian_rand@plt>
mov    -0x18(%rbp),%rax
movslq -0x38(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
lea    -0x4c(%rbp),%rdi
lea    -0x50(%rbp),%rsi
lea    -0x48(%rbp),%rdx
call   1040 <gaussian_rand@plt>
mov    -0x20(%rbp),%rax
movslq -0x38(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
movss  -0x28(%rbp),%xmm0
mov    -0x18(%rbp),%rax
movslq -0x38(%rbp),%rcx
addss  (%rax,%rcx,4),%xmm0
movss  %xmm0,-0x28(%rbp)
movss  -0x2c(%rbp),%xmm0
mov    -0x20(%rbp),%rax
movslq -0x38(%rbp),%rcx
addss  (%rax,%rcx,4),%xmm0
movss  %xmm0,-0x2c(%rbp)
movss  -0x24(%rbp),%xmm0
mov    -0x18(%rbp),%rax
movslq -0x38(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm1
mov    -0x18(%rbp),%rax
movslq -0x38(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm3
mov    -0x20(%rbp),%rax
movslq -0x38(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm2
mov    -0x20(%rbp),%rax
movslq -0x38(%rbp),%rcx
mulss  (%rax,%rcx,4),%xmm2
mulss  %xmm3,%xmm1
addss  %xmm2,%xmm1
addss  %xmm1,%xmm0
movss  %xmm0,-0x24(%rbp)
mov    -0x38(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x38(%rbp)
jmp    11c8 <main+0x78>
movss  -0x28(%rbp),%xmm0
cvtsi2ssl -0x3c(%rbp),%xmm1
divss  %xmm1,%xmm0
movss  %xmm0,-0x28(%rbp)
movss  -0x2c(%rbp),%xmm0
cvtsi2ssl -0x3c(%rbp),%xmm1
divss  %xmm1,%xmm0
movss  %xmm0,-0x2c(%rbp)
xorps  %xmm0,%xmm0
movss  %xmm0,-0x30(%rbp)
xorps  %xmm0,%xmm0
movss  %xmm0,-0x34(%rbp)
movl   $0x0,-0x38(%rbp)
mov    -0x38(%rbp),%eax
cmp    -0x3c(%rbp),%eax
jge    1359 <main+0x209>
movss  -0x30(%rbp),%xmm1
mov    -0x18(%rbp),%rax
movslq -0x38(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
subss  -0x28(%rbp),%xmm0
mov    -0x18(%rbp),%rax
movslq -0x38(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm2
subss  -0x28(%rbp),%xmm2
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,-0x30(%rbp)
movss  -0x34(%rbp),%xmm1
mov    -0x20(%rbp),%rax
movslq -0x38(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
subss  -0x2c(%rbp),%xmm0
mov    -0x20(%rbp),%rax
movslq -0x38(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm2
subss  -0x2c(%rbp),%xmm2
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,-0x34(%rbp)
mov    -0x38(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x38(%rbp)
jmp    12d3 <main+0x183>
movss  -0x30(%rbp),%xmm0
mov    -0x3c(%rbp),%eax
sub    $0x1,%eax
cvtsi2ss %eax,%xmm1
divss  %xmm1,%xmm0
call   1070 <sqrt@plt>
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,-0x30(%rbp)
movss  -0x34(%rbp),%xmm0
mov    -0x3c(%rbp),%eax
sub    $0x1,%eax
cvtsi2ss %eax,%xmm1
divss  %xmm1,%xmm0
call   1070 <sqrt@plt>
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,-0x34(%rbp)
movss  -0x28(%rbp),%xmm0
cvtss2sd %xmm0,%xmm0
movss  -0x30(%rbp),%xmm1
cvtss2sd %xmm1,%xmm1
lea    0xc53(%rip),%rdi
mov    $0x2,%al
call   1030 <printf@plt>
movss  -0x2c(%rbp),%xmm0
cvtss2sd %xmm0,%xmm0
movss  -0x34(%rbp),%xmm1
cvtss2sd %xmm1,%xmm1
lea    0xc48(%rip),%rdi
mov    $0x2,%al
call   1030 <printf@plt>
movss  -0x24(%rbp),%xmm0
cvtsi2ssl -0x3c(%rbp),%xmm1
divss  %xmm1,%xmm0
call   1070 <sqrt@plt>
lea    0xc3c(%rip),%rdi
mov    $0x1,%al
call   1030 <printf@plt>
mov    -0x4(%rbp),%eax
add    $0x50,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<main>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x48,%rsp
movq   $0x0,0x40(%rsp)
movl   $0x3f800000,0x1c(%rsp)
3f
movl   $0x0,0x18(%rsp)
lea    0xe92(%rip),%rdi
lea    0xc(%rsp),%rsi
call   1060 <scanf@plt>
mov    0xc(%rsp),%edi
shl    $0x2,%edi
call   1050 <check_malloc@plt>
mov    %rax,%rbx
mov    0xc(%rsp),%edi
shl    $0x2,%edi
call   1050 <check_malloc@plt>
mov    %rax,%r14
mov    0xc(%rsp),%eax
test   %eax,%eax
jle    1299 <main+0x149>
xorps  %xmm1,%xmm1
xor    %ebp,%ebp
lea    0x1c(%rsp),%r15
lea    0x18(%rsp),%r12
lea    0x40(%rsp),%r13
cs nopw 0x0(%rax,%rax,1)
movaps %xmm3,0x30(%rsp)
movss  %xmm1,0x10(%rsp)
mov    %r15,%rdi
mov    %r12,%rsi
mov    %r13,%rdx
call   1040 <gaussian_rand@plt>
movss  %xmm0,(%rbx,%rbp,4)
mov    %r15,%rdi
mov    %r12,%rsi
mov    %r13,%rdx
call   1040 <gaussian_rand@plt>
movaps 0x30(%rsp),%xmm3
movss  %xmm0,(%r14,%rbp,4)
movss  (%rbx,%rbp,4),%xmm1
movaps %xmm0,%xmm2
unpcklps %xmm1,%xmm2
addps  %xmm2,%xmm3
mulss  %xmm0,%xmm0
mulss  %xmm1,%xmm1
addss  %xmm0,%xmm1
movss  0x10(%rsp),%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,0x10(%rsp)
movss  0x10(%rsp),%xmm1
add    $0x1,%rbp
movslq 0xc(%rsp),%rax
cmp    %rax,%rbp
jl     11d0 <main+0x80>
xorps  %xmm0,%xmm0
cvtsi2ss %eax,%xmm0
shufps $0x0,%xmm0,%xmm0
divps  %xmm0,%xmm3
test   %eax,%eax
movaps %xmm3,0x30(%rsp)
jle    12bb <main+0x16b>
mov    %eax,%ecx
cmp    $0x1,%eax
jne    12c8 <main+0x178>
xorps  %xmm4,%xmm4
xor    %edx,%edx
test   $0x1,%cl
je     132c <main+0x1dc>
movss  (%rbx,%rdx,4),%xmm0
movss  (%r14,%rdx,4),%xmm1
unpcklps %xmm0,%xmm1
subps  %xmm3,%xmm1
mulps  %xmm1,%xmm1
addps  %xmm1,%xmm4
movaps %xmm4,0x20(%rsp)
movaps %xmm4,%xmm0
shufps $0x55,%xmm4,%xmm0
jmp    1331 <main+0x1e1>
movaps 0xd60(%rip),%xmm0
movaps %xmm0,0x30(%rsp)
xorps  %xmm0,%xmm0
movss  %xmm0,0x10(%rsp)
xorps  %xmm0,%xmm0
xorps  %xmm1,%xmm1
movaps %xmm1,0x20(%rsp)
jmp    1331 <main+0x1e1>
xorps  %xmm0,%xmm0
xorps  %xmm2,%xmm2
movaps %xmm2,0x20(%rsp)
jmp    1331 <main+0x1e1>
mov    %ecx,%esi
and    $0xfffffffe,%esi
xorps  %xmm4,%xmm4
xor    %edx,%edx
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
movss  (%rbx,%rdx,4),%xmm0
movss  0x4(%rbx,%rdx,4),%xmm1
movss  (%r14,%rdx,4),%xmm2
unpcklps %xmm0,%xmm2
subps  %xmm3,%xmm2
mulps  %xmm2,%xmm2
addps  %xmm4,%xmm2
movss  0x4(%r14,%rdx,4),%xmm0
unpcklps %xmm1,%xmm0
subps  %xmm3,%xmm0
mulps  %xmm0,%xmm0
movaps %xmm0,%xmm4
addps  %xmm2,%xmm4
add    $0x2,%rdx
cmp    %rdx,%rsi
jne    12e0 <main+0x190>
movaps %xmm4,%xmm0
shufps $0x55,%xmm4,%xmm0
test   $0x1,%cl
jne    1271 <main+0x121>
movaps %xmm4,0x20(%rsp)
add    $0xffffffff,%eax
xorps  %xmm1,%xmm1
cvtsi2ss %eax,%xmm1
divss  %xmm1,%xmm0
call   1070 <sqrt@plt>
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,0x14(%rsp)
mov    0xc(%rsp),%eax
add    $0xffffffff,%eax
xorps  %xmm1,%xmm1
cvtsi2ss %eax,%xmm1
movaps 0x20(%rsp),%xmm0
divss  %xmm1,%xmm0
call   1070 <sqrt@plt>
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,0x20(%rsp)
movaps 0x30(%rsp),%xmm0
shufps $0x55,%xmm0,%xmm0
cvtss2sd %xmm0,%xmm0
movss  0x14(%rsp),%xmm1
cvtss2sd %xmm1,%xmm1
lea    0xc81(%rip),%rdi
mov    $0x2,%al
call   1030 <printf@plt>
movaps 0x30(%rsp),%xmm0
cvtss2sd %xmm0,%xmm0
movss  0x20(%rsp),%xmm1
cvtss2sd %xmm1,%xmm1
lea    0xc75(%rip),%rdi
mov    $0x2,%al
call   1030 <printf@plt>
xorps  %xmm1,%xmm1
cvtsi2ssl 0xc(%rsp),%xmm1
movss  0x10(%rsp),%xmm0
divss  %xmm1,%xmm0
call   1070 <sqrt@plt>
lea    0xc64(%rip),%rdi
mov    $0x1,%al
call   1030 <printf@plt>
xor    %eax,%eax
add    $0x48,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

LLVM IR (-O2):
@.str = [3 x i8] c"%d\00"
@.str.1 = [21 x i8] c"rmean= %f, rsig= %f\0A\00"
@.str.2 = [21 x i8] c"imean= %f, isig= %f\0A\00"
@.str.3 = [14 x i8] c"avg. amp= %f\0A\00"
define i32 @main(i32 %0, i8** nocapture readnone %1) {
%3 = alloca i32
%4 = alloca i64
%5 = alloca float
%6 = alloca float
%7 = bitcast i32* %3 to i8*
%8 = bitcast i64* %4 to i8*
store i64 0, i64* %4
%9 = bitcast float* %5 to i8*
store float 1.000000e+00, float* %5
%10 = bitcast float* %6 to i8*
store float 0.000000e+00, float* %6
%11 = call i32 @scanf(i8* getelementptr inbounds ([3 x i8], [3 x i8]* @.str, i64 0, i64 0), i32* %3)
%12 = load i32, i32* %3
%13 = shl i32 %12, 2
%14 = call i64 @check_malloc(i32 %13)
%15 = inttoptr i64 %14 to float*
%16 = load i32, i32* %3
%17 = shl i32 %16, 2
%18 = call i64 @check_malloc(i32 %17)
%19 = inttoptr i64 %18 to float*
%20 = load i32, i32* %3
%21 = icmp sgt i32 %20, 0
br i1 %21, label %22, label %97
22: ; preds = %2, %22
%23 = phi i64 [ %37, %22 ], [ 0, %2 ]
%24 = phi float [ %36, %22 ], [ 0.000000e+00, %2 ]
%25 = phi <2 x float> [ %33, %22 ], [ undef, %2 ]
%26 = call float @gaussian_rand(float* %5, float* %6, i64* %4)
%27 = getelementptr inbounds float, float* %15, i64 %23
store float %26, float* %27
%28 = call float @gaussian_rand(float* %5, float* %6, i64* %4)
%29 = getelementptr inbounds float, float* %19, i64 %23
store float %28, float* %29
%30 = load float, float* %27
%31 = insertelement <2 x float> poison, float %28, i64 0
%32 = insertelement <2 x float> %31, float %30, i64 1
%33 = fadd <2 x float> %25, %32
%34 = fmul float %28, %28
%35 = call float @llvm.fmuladd.f32(float %30, float %30, float %34)
%36 = fadd float %24, %35
%37 = add nuw nsw i64 %23, 1
%38 = load i32, i32* %3
%39 = sext i32 %38 to i64
%40 = icmp slt i64 %37, %39
br i1 %40, label %22, label %41
41: ; preds = %22
%42 = sitofp i32 %38 to float
%43 = insertelement <2 x float> poison, float %42, i64 0
%44 = shufflevector <2 x float> %43, <2 x float> poison, <2 x i32> zeroinitializer
%45 = fdiv <2 x float> %33, %44
%46 = icmp sgt i32 %38, 0
br i1 %46, label %47, label %97
47: ; preds = %41
%48 = zext i32 %38 to i64
%49 = and i64 %48, 1
%50 = icmp eq i32 %38, 1
br i1 %50, label %80, label %51
51: ; preds = %47
%52 = and i64 %48, 4294967294
br label %53
53: ; preds = %53, %51
%54 = phi i64 [ 0, %51 ], [ %74, %53 ]
%55 = phi <2 x float> [ zeroinitializer, %51 ], [ %73, %53 ]
%56 = phi i64 [ 0, %51 ], [ %75, %53 ]
%57 = getelementptr inbounds float, float* %15, i64 %54
%58 = load float, float* %57
%59 = getelementptr inbounds float, float* %19, i64 %54
%60 = load float, float* %59
%61 = insertelement <2 x float> poison, float %60, i64 0
%62 = insertelement <2 x float> %61, float %58, i64 1
%63 = fsub <2 x float> %62, %45
%64 = call <2 x float> @llvm.fmuladd.v2f32(<2 x float> %63, <2 x float> %63, <2 x float> %55)
%65 = or i64 %54, 1
%66 = getelementptr inbounds float, float* %15, i64 %65
%67 = load float, float* %66
%68 = getelementptr inbounds float, float* %19, i64 %65
%69 = load float, float* %68
%70 = insertelement <2 x float> poison, float %69, i64 0
%71 = insertelement <2 x float> %70, float %67, i64 1
%72 = fsub <2 x float> %71, %45
%73 = call <2 x float> @llvm.fmuladd.v2f32(<2 x float> %72, <2 x float> %72, <2 x float> %64)
%74 = add nuw nsw i64 %54, 2
%75 = add i64 %56, 2
%76 = icmp eq i64 %75, %52
br i1 %76, label %77, label %53
77: ; preds = %53
%78 = extractelement <2 x float> %73, i64 1
%79 = extractelement <2 x float> %73, i64 0
br label %80
80: ; preds = %77, %47
%81 = phi float [ undef, %47 ], [ %79, %77 ]
%82 = phi float [ undef, %47 ], [ %78, %77 ]
%83 = phi i64 [ 0, %47 ], [ %74, %77 ]
%84 = phi <2 x float> [ zeroinitializer, %47 ], [ %73, %77 ]
%85 = icmp eq i64 %49, 0
br i1 %85, label %97, label %86
86: ; preds = %80
%87 = getelementptr inbounds float, float* %15, i64 %83
%88 = load float, float* %87
%89 = getelementptr inbounds float, float* %19, i64 %83
%90 = load float, float* %89
%91 = insertelement <2 x float> poison, float %90, i64 0
%92 = insertelement <2 x float> %91, float %88, i64 1
%93 = fsub <2 x float> %92, %45
%94 = call <2 x float> @llvm.fmuladd.v2f32(<2 x float> %93, <2 x float> %93, <2 x float> %84)
%95 = extractelement <2 x float> %94, i64 0
%96 = extractelement <2 x float> %94, i64 1
br label %97
97: ; preds = %86, %80, %2, %41
%98 = phi i32 [ %38, %41 ], [ %20, %2 ], [ %38, %80 ], [ %38, %86 ]
%99 = phi float [ %36, %41 ], [ 0.000000e+00, %2 ], [ %36, %80 ], [ %36, %86 ]
%100 = phi float [ 0.000000e+00, %41 ], [ 0.000000e+00, %2 ], [ %82, %80 ], [ %96, %86 ]
%101 = phi float [ 0.000000e+00, %41 ], [ 0.000000e+00, %2 ], [ %81, %80 ], [ %95, %86 ]
%102 = phi <2 x float> [ %45, %41 ], [ <float 0x7FF8000000000000, float 0x7FF8000000000000>, %2 ], [ %45, %80 ], [ %45, %86 ]
%103 = add nsw i32 %98, -1
%104 = sitofp i32 %103 to float
%105 = fdiv float %100, %104
%106 = call double @sqrt(float %105)
%107 = fptrunc double %106 to float
%108 = load i32, i32* %3
%109 = add nsw i32 %108, -1
%110 = sitofp i32 %109 to float
%111 = fdiv float %101, %110
%112 = call double @sqrt(float %111)
%113 = fptrunc double %112 to float
%114 = extractelement <2 x float> %102, i64 1
%115 = fpext float %114 to double
%116 = fpext float %107 to double
%117 = call i32 (i8*, double, ...) @printf(i8* getelementptr inbounds ([21 x i8], [21 x i8]* @.str.1, i64 0, i64 0), double %115, double %116)
%118 = extractelement <2 x float> %102, i64 0
%119 = fpext float %118 to double
%120 = fpext float %113 to double
%121 = call i32 (i8*, double, ...) @printf(i8* getelementptr inbounds ([21 x i8], [21 x i8]* @.str.2, i64 0, i64 0), double %119, double %120)
%122 = load i32, i32* %3
%123 = sitofp i32 %122 to float
%124 = fdiv float %99, %123
%125 = call double @sqrt(float %124)
%126 = call i32 (i8*, double, ...) @printf(i8* getelementptr inbounds ([14 x i8], [14 x i8]* @.str.3, i64 0, i64 0), double %125)
ret i32 0
}