#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 */
float exp(double);
float sqrt(double);
float Normal_nD(int ndim, float *mean, float *sigma, float *x) {
  int i;
  float p, srpi2;
  srpi2 = sqrt(6.283185307179586);
  p = 1.0;
  for (i = 0; i < ndim; i++) {
    p *= (exp(-0.5 * (x[i] - mean[i]) * (x[i] - mean[i]) /
              (sigma[i] * sigma[i])) /
          (srpi2 * sigma[i]));
  }
  return (p);
}


x86 assembly (clang, -O0):
<Normal_nD>:
push   %rbp
mov    %rsp,%rbp
sub    $0x30,%rsp
mov    %edi,-0x4(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %rcx,-0x20(%rbp)
movsd  0xec9(%rip),%xmm0
call   1040 <sqrt@plt>
movss  %xmm0,-0x2c(%rbp)
movss  0xeaf(%rip),%xmm0
movss  %xmm0,-0x28(%rbp)
movl   $0x0,-0x24(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    120a <Normal_nD+0xea>
mov    -0x20(%rbp),%rax
movslq -0x24(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x24(%rbp),%rcx
subss  (%rax,%rcx,4),%xmm0
cvtss2sd %xmm0,%xmm1
movsd  0xe81(%rip),%xmm0
mulsd  %xmm1,%xmm0
mov    -0x20(%rbp),%rax
movslq -0x24(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm1
mov    -0x10(%rbp),%rax
movslq -0x24(%rbp),%rcx
subss  (%rax,%rcx,4),%xmm1
cvtss2sd %xmm1,%xmm1
mulsd  %xmm1,%xmm0
mov    -0x18(%rbp),%rax
movslq -0x24(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm1
mov    -0x18(%rbp),%rax
movslq -0x24(%rbp),%rcx
mulss  (%rax,%rcx,4),%xmm1
cvtss2sd %xmm1,%xmm1
divsd  %xmm1,%xmm0
call   1030 <exp@plt>
movss  -0x2c(%rbp),%xmm1
mov    -0x18(%rbp),%rax
movslq -0x24(%rbp),%rcx
mulss  (%rax,%rcx,4),%xmm1
divss  %xmm1,%xmm0
mulss  -0x28(%rbp),%xmm0
movss  %xmm0,-0x28(%rbp)
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    115d <Normal_nD+0x3d>
movss  -0x28(%rbp),%xmm0
add    $0x30,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<Normal_nD>:
push   %rbp
push   %r15
push   %r14
push   %r12
push   %rbx
sub    $0x10,%rsp
mov    %rcx,%r14
mov    %rdx,%rbx
mov    %rsi,%r15
mov    %edi,%ebp
movsd  0xec1(%rip),%xmm0
call   1040 <sqrt@plt>
test   %ebp,%ebp
jle    11e2 <Normal_nD+0xc2>
movaps %xmm0,%xmm2
mov    %ebp,%r12d
movss  0xeb6(%rip),%xmm0
xor    %ebp,%ebp
movss  %xmm2,0xc(%rsp)
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
movss  %xmm0,0x8(%rsp)
movss  (%r14,%rbp,4),%xmm0
subss  (%r15,%rbp,4),%xmm0
xorps  %xmm1,%xmm1
cvtss2sd %xmm0,%xmm1
movaps %xmm1,%xmm0
mulsd  0xe74(%rip),%xmm0
mulsd  %xmm1,%xmm0
movss  (%rbx,%rbp,4),%xmm1
mulss  %xmm1,%xmm1
cvtss2sd %xmm1,%xmm1
divsd  %xmm1,%xmm0
call   1030 <exp@plt>
movss  0xc(%rsp),%xmm2
movss  (%rbx,%rbp,4),%xmm1
mulss  %xmm2,%xmm1
divss  %xmm1,%xmm0
movss  0x8(%rsp),%xmm1
mulss  %xmm0,%xmm1
movss  %xmm1,0x8(%rsp)
movss  0x8(%rsp),%xmm0
add    $0x1,%rbp
cmp    %rbp,%r12
jne    1170 <Normal_nD+0x50>
jmp    11ea <Normal_nD+0xca>
movss  0xe26(%rip),%xmm0
add    $0x10,%rsp
pop    %rbx
pop    %r12
pop    %r14
pop    %r15
pop    %rbp
ret

LLVM IR (-O0):
define float @Normal_nD(i32 %0, float* %1, float* %2, float* %3) {
%5 = alloca i32
%6 = alloca float*
%7 = alloca float*
%8 = alloca float*
%9 = alloca i32
%10 = alloca float
%11 = alloca float
store i32 %0, i32* %5
store float* %1, float** %6
store float* %2, float** %7
store float* %3, float** %8
%12 = call float @sqrt(double 0x401921FB54442D18)
store float %12, float* %11
store float 1.000000e+00, float* %10
store i32 0, i32* %9
br label %13
13: ; preds = %68, %4
%14 = load i32, i32* %9
%15 = load i32, i32* %5
%16 = icmp slt i32 %14, %15
br i1 %16, label %17, label %71
17: ; preds = %13
%18 = load float*, float** %8
%19 = load i32, i32* %9
%20 = sext i32 %19 to i64
%21 = getelementptr inbounds float, float* %18, i64 %20
%22 = load float, float* %21
%23 = load float*, float** %6
%24 = load i32, i32* %9
%25 = sext i32 %24 to i64
%26 = getelementptr inbounds float, float* %23, i64 %25
%27 = load float, float* %26
%28 = fsub float %22, %27
%29 = fpext float %28 to double
%30 = fmul double -5.000000e-01, %29
%31 = load float*, float** %8
%32 = load i32, i32* %9
%33 = sext i32 %32 to i64
%34 = getelementptr inbounds float, float* %31, i64 %33
%35 = load float, float* %34
%36 = load float*, float** %6
%37 = load i32, i32* %9
%38 = sext i32 %37 to i64
%39 = getelementptr inbounds float, float* %36, i64 %38
%40 = load float, float* %39
%41 = fsub float %35, %40
%42 = fpext float %41 to double
%43 = fmul double %30, %42
%44 = load float*, float** %7
%45 = load i32, i32* %9
%46 = sext i32 %45 to i64
%47 = getelementptr inbounds float, float* %44, i64 %46
%48 = load float, float* %47
%49 = load float*, float** %7
%50 = load i32, i32* %9
%51 = sext i32 %50 to i64
%52 = getelementptr inbounds float, float* %49, i64 %51
%53 = load float, float* %52
%54 = fmul float %48, %53
%55 = fpext float %54 to double
%56 = fdiv double %43, %55
%57 = call float @exp(double %56)
%58 = load float, float* %11
%59 = load float*, float** %7
%60 = load i32, i32* %9
%61 = sext i32 %60 to i64
%62 = getelementptr inbounds float, float* %59, i64 %61
%63 = load float, float* %62
%64 = fmul float %58, %63
%65 = fdiv float %57, %64
%66 = load float, float* %10
%67 = fmul float %66, %65
store float %67, float* %10
br label %68
68: ; preds = %17
%69 = load i32, i32* %9
%70 = add nsw i32 %69, 1
store i32 %70, i32* %9
br label %13
71: ; preds = %13
%72 = load float, float* %10
ret float %72
}

LLVM IR (-O2):
define float @Normal_nD(i32 %0, float* nocapture %1, float* nocapture %2, float* nocapture %3) {
%5 = tail call float @sqrt(double 0x401921FB54442D18)
%6 = icmp sgt i32 %0, 0
br i1 %6, label %7, label %32
7: ; preds = %4
%8 = zext i32 %0 to i64
br label %9
9: ; preds = %7, %9
%10 = phi i64 [ 0, %7 ], [ %30, %9 ]
%11 = phi float [ 1.000000e+00, %7 ], [ %29, %9 ]
%12 = getelementptr inbounds float, float* %3, i64 %10
%13 = load float, float* %12
%14 = getelementptr inbounds float, float* %1, i64 %10
%15 = load float, float* %14
%16 = fsub float %13, %15
%17 = fpext float %16 to double
%18 = fmul double %17, -5.000000e-01
%19 = fmul double %18, %17
%20 = getelementptr inbounds float, float* %2, i64 %10
%21 = load float, float* %20
%22 = fmul float %21, %21
%23 = fpext float %22 to double
%24 = fdiv double %19, %23
%25 = tail call float @exp(double %24)
%26 = load float, float* %20
%27 = fmul float %5, %26
%28 = fdiv float %25, %27
%29 = fmul float %11, %28
%30 = add nuw nsw i64 %10, 1
%31 = icmp eq i64 %30, %8
br i1 %31, label %32, label %9
32: ; preds = %9, %4
%33 = phi float [ 1.000000e+00, %4 ], [ %29, %9 ]
ret float %33
}