#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__ ceil(double);
int /*<<< orphan*/ do_vect(double **, unsigned long);
double exp(double);
int gidz_cube(int, int, int, unsigned int);
int /*<<< orphan*/ printf(char *, unsigned int, unsigned long);
void create_gaussian(double **phi, unsigned long *nvox, unsigned *ext,
                     double sigmap, double sigma_factor) {
  int exth, indx, indy, indz;
  double dsqu;
  double mscale;
  unsigned long count;
  double bvalue, cvalue;
  exth = (int)ceil(sigma_factor * sigmap);
  *ext = 2 * exth - 1;
  *nvox = *ext * *ext * *ext;
  printf("lib_vwk> Generating Gaussian kernel with %d^3 = %lu voxels.\n", *ext,
         *nvox);
  do_vect(phi, *nvox);
  bvalue = -1.0 / (2.0 * sigmap * sigmap);
  cvalue = sigma_factor * sigma_factor * sigmap * sigmap;
  mscale = 0.0;
  for (indz = 0; indz < (int)*ext; indz++)
    for (indy = 0; indy < (int)*ext; indy++)
      for (indx = 0; indx < (int)*ext; indx++) {
        dsqu = (indx - exth + 1) * (indx - exth + 1) +
               (indy - exth + 1) * (indy - exth + 1) +
               (indz - exth + 1) * (indz - exth + 1);
        if (dsqu <= cvalue)
          *(*phi + gidz_cube(indz, indy, indx, *ext)) = exp(dsqu * bvalue);
        mscale += *(*phi + gidz_cube(indz, indy, indx, *ext));
      }
  for (count = 0; count < *nvox; count++)
    *(*phi + count) /= mscale;
}


x86 assembly (clang, -O0):
<create_gaussian>:
push   %rbp
mov    %rsp,%rbp
sub    $0x80,%rsp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
movsd  %xmm0,-0x20(%rbp)
movsd  %xmm1,-0x28(%rbp)
movsd  -0x28(%rbp),%xmm0
mulsd  -0x20(%rbp),%xmm0
call   1030 <ceil@plt>
mov    %eax,-0x2c(%rbp)
mov    -0x2c(%rbp),%ecx
shl    $0x1,%ecx
sub    $0x1,%ecx
mov    -0x18(%rbp),%rax
mov    %ecx,(%rax)
mov    -0x18(%rbp),%rax
mov    (%rax),%eax
mov    -0x18(%rbp),%rcx
imul   (%rcx),%eax
mov    -0x18(%rbp),%rcx
imul   (%rcx),%eax
mov    %eax,%eax
mov    %eax,%ecx
mov    -0x10(%rbp),%rax
mov    %rcx,(%rax)
mov    -0x18(%rbp),%rax
mov    (%rax),%esi
mov    -0x10(%rbp),%rax
mov    (%rax),%rdx
lea    0xe4b(%rip),%rdi
call   1040 <printf@plt>
mov    -0x8(%rbp),%rdi
mov    -0x10(%rbp),%rax
mov    (%rax),%rsi
call   1060 <do_vect@plt>
movsd  0xe26(%rip),%xmm1
mulsd  -0x20(%rbp),%xmm1
mulsd  -0x20(%rbp),%xmm1
movsd  0xe0c(%rip),%xmm0
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x58(%rbp)
movsd  -0x28(%rbp),%xmm0
mulsd  -0x28(%rbp),%xmm0
mulsd  -0x20(%rbp),%xmm0
mulsd  -0x20(%rbp),%xmm0
movsd  %xmm0,-0x60(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x48(%rbp)
movl   $0x0,-0x38(%rbp)
mov    -0x38(%rbp),%eax
mov    -0x18(%rbp),%rcx
cmp    (%rcx),%eax
jge    1371 <create_gaussian+0x221>
movl   $0x0,-0x34(%rbp)
mov    -0x34(%rbp),%eax
mov    -0x18(%rbp),%rcx
cmp    (%rcx),%eax
jge    135e <create_gaussian+0x20e>
movl   $0x0,-0x30(%rbp)
mov    -0x30(%rbp),%eax
mov    -0x18(%rbp),%rcx
cmp    (%rcx),%eax
jge    134b <create_gaussian+0x1fb>
mov    -0x30(%rbp),%eax
sub    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    -0x30(%rbp),%ecx
sub    -0x2c(%rbp),%ecx
add    $0x1,%ecx
imul   %ecx,%eax
mov    -0x34(%rbp),%ecx
sub    -0x2c(%rbp),%ecx
add    $0x1,%ecx
mov    -0x34(%rbp),%edx
sub    -0x2c(%rbp),%edx
add    $0x1,%edx
imul   %edx,%ecx
add    %ecx,%eax
mov    -0x38(%rbp),%ecx
sub    -0x2c(%rbp),%ecx
add    $0x1,%ecx
mov    -0x38(%rbp),%edx
sub    -0x2c(%rbp),%edx
add    $0x1,%edx
imul   %edx,%ecx
add    %ecx,%eax
cvtsi2sd %eax,%xmm0
movsd  %xmm0,-0x40(%rbp)
movsd  -0x40(%rbp),%xmm1
movsd  -0x60(%rbp),%xmm0
ucomisd %xmm1,%xmm0
jb     1306 <create_gaussian+0x1b6>
movsd  -0x40(%rbp),%xmm0
mulsd  -0x58(%rbp),%xmm0
call   1050 <exp@plt>
movsd  %xmm0,-0x70(%rbp)
mov    -0x8(%rbp),%rax
mov    (%rax),%rax
mov    %rax,-0x68(%rbp)
mov    -0x38(%rbp),%edi
mov    -0x34(%rbp),%esi
mov    -0x30(%rbp),%edx
mov    -0x18(%rbp),%rax
mov    (%rax),%ecx
call   1070 <gidz_cube@plt>
movsd  -0x70(%rbp),%xmm0
mov    %eax,%ecx
mov    -0x68(%rbp),%rax
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x8(%rbp),%rax
mov    (%rax),%rax
mov    %rax,-0x78(%rbp)
mov    -0x38(%rbp),%edi
mov    -0x34(%rbp),%esi
mov    -0x30(%rbp),%edx
mov    -0x18(%rbp),%rax
mov    (%rax),%ecx
call   1070 <gidz_cube@plt>
mov    %eax,%ecx
mov    -0x78(%rbp),%rax
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
addsd  -0x48(%rbp),%xmm0
movsd  %xmm0,-0x48(%rbp)
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
jmp    1251 <create_gaussian+0x101>
jmp    1350 <create_gaussian+0x200>
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    123b <create_gaussian+0xeb>
jmp    1363 <create_gaussian+0x213>
mov    -0x38(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x38(%rbp)
jmp    1225 <create_gaussian+0xd5>
movq   $0x0,-0x50(%rbp)
mov    -0x50(%rbp),%rax
mov    -0x10(%rbp),%rcx
cmp    (%rcx),%rax
jae    13b9 <create_gaussian+0x269>
movsd  -0x48(%rbp),%xmm1
mov    -0x8(%rbp),%rax
mov    (%rax),%rax
mov    -0x50(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
divsd  %xmm1,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x50(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x50(%rbp)
jmp    1379 <create_gaussian+0x229>
add    $0x80,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<create_gaussian>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x58,%rsp
mov    %rdx,%r13
mov    %rsi,%rbx
mov    %rdi,%r15
movsd  %xmm0,0x20(%rsp)
movsd  %xmm1,0x8(%rsp)
mulsd  %xmm1,%xmm0
call   1030 <ceil@plt>
mov    %rax,%rbp
lea    (%rax,%rax,1),%esi
add    $0xffffffff,%esi
mov    %esi,0x0(%r13)
mov    %esi,%edx
imul   %edx,%edx
imul   %esi,%edx
mov    %rdx,(%rbx)
lea    0xe6d(%rip),%rdi
call   1040 <printf@plt>
mov    %rbx,0x38(%rsp)
mov    (%rbx),%rsi
mov    %r15,%rdi
call   1060 <do_vect@plt>
mov    0x0(%r13),%ecx
test   %ecx,%ecx
jle    1315 <create_gaussian+0x1c5>
movsd  0x20(%rsp),%xmm2
movapd %xmm2,%xmm0
addsd  %xmm2,%xmm0
mulsd  %xmm2,%xmm0
movsd  0xe2a(%rip),%xmm1
divsd  %xmm0,%xmm1
movsd  %xmm1,0x40(%rsp)
movsd  0x8(%rsp),%xmm1
mulsd  %xmm1,%xmm1
mulsd  %xmm2,%xmm1
mulsd  %xmm2,%xmm1
mov    $0x1,%eax
sub    %ebp,%eax
mov    %rax,0x50(%rsp)
movl   $0x0,0x4(%rsp)
xorpd  %xmm3,%xmm3
mov    %rbp,0x18(%rsp)
movsd  %xmm1,0x8(%rsp)
jmp    1238 <create_gaussian+0xe8>
nopw   0x0(%rax,%rax,1)
mov    0x4(%rsp),%eax
add    $0x1,%eax
mov    %eax,0x4(%rsp)
cmp    %ecx,%eax
mov    0x18(%rsp),%rbp
jge    1319 <create_gaussian+0x1c9>
test   %ecx,%ecx
jle    1220 <create_gaussian+0xd0>
mov    0x4(%rsp),%eax
sub    %ebp,%eax
add    $0x1,%eax
imul   %eax,%eax
mov    %eax,0x14(%rsp)
xor    %ebp,%ebp
jmp    1257 <create_gaussian+0x107>
add    $0x1,%ebp
cmp    %ecx,%ebp
jge    1220 <create_gaussian+0xd0>
test   %ecx,%ecx
jle    1250 <create_gaussian+0x100>
mov    %ebp,%r12d
sub    0x18(%rsp),%r12d
add    $0x1,%r12d
imul   %r12d,%r12d
add    0x14(%rsp),%r12d
xor    %r14d,%r14d
jmp    12b1 <create_gaussian+0x161>
cs nopw 0x0(%rax,%rax,1)
nop
mov    (%r15),%rbx
mov    0x4(%rsp),%edi
mov    %ebp,%esi
mov    %r14d,%edx
call   1070 <gidz_cube@plt>
cltq
movapd 0x20(%rsp),%xmm3
addsd  (%rbx,%rax,8),%xmm3
add    $0x1,%r14d
mov    0x0(%r13),%ecx
cmp    %ecx,%r14d
movsd  0x8(%rsp),%xmm1
jge    1250 <create_gaussian+0x100>
movapd %xmm3,0x20(%rsp)
mov    0x50(%rsp),%rax
add    %r14d,%eax
imul   %eax,%eax
add    %r12d,%eax
xorps  %xmm0,%xmm0
cvtsi2sd %eax,%xmm0
ucomisd %xmm0,%xmm1
jb     1280 <create_gaussian+0x130>
mulsd  0x40(%rsp),%xmm0
call   1050 <exp@plt>
movsd  %xmm0,0x48(%rsp)
mov    %r15,%rbx
mov    (%r15),%r15
mov    0x0(%r13),%ecx
mov    0x4(%rsp),%edi
mov    %ebp,%esi
mov    %r14d,%edx
call   1070 <gidz_cube@plt>
cltq
movsd  0x48(%rsp),%xmm0
movsd  %xmm0,(%r15,%rax,8)
mov    %rbx,%r15
mov    0x0(%r13),%ecx
jmp    1280 <create_gaussian+0x130>
xorpd  %xmm3,%xmm3
mov    0x38(%rsp),%rax
mov    (%rax),%rax
test   %rax,%rax
je     13cf <create_gaussian+0x27f>
mov    (%r15),%rcx
cmp    $0x1,%rax
jne    133a <create_gaussian+0x1ea>
xor    %edx,%edx
jmp    13b8 <create_gaussian+0x268>
mov    %rax,%rdx
and    $0xfffffffffffffffe,%rdx
movapd %xmm3,%xmm0
unpcklpd %xmm3,%xmm0
lea    -0x2(%rdx),%rsi
mov    %rsi,%rbp
shr    %rbp
add    $0x1,%rbp
test   %rsi,%rsi
je     13a0 <create_gaussian+0x250>
mov    %rbp,%rsi
and    $0xfffffffffffffffe,%rsi
xor    %edi,%edi
cs nopw 0x0(%rax,%rax,1)
nop
movupd (%rcx,%rdi,8),%xmm1
movupd 0x10(%rcx,%rdi,8),%xmm2
divpd  %xmm0,%xmm1
movupd %xmm1,(%rcx,%rdi,8)
divpd  %xmm0,%xmm2
movupd %xmm2,0x10(%rcx,%rdi,8)
add    $0x4,%rdi
add    $0xfffffffffffffffe,%rsi
jne    1370 <create_gaussian+0x220>
test   $0x1,%bpl
jne    13a8 <create_gaussian+0x258>
jmp    13ca <create_gaussian+0x27a>
xor    %edi,%edi
test   $0x1,%bpl
je     13ca <create_gaussian+0x27a>
movupd (%rcx,%rdi,8),%xmm1
divpd  %xmm0,%xmm1
movupd %xmm1,(%rcx,%rdi,8)
jmp    13ca <create_gaussian+0x27a>
movsd  (%rcx,%rdx,8),%xmm0
divsd  %xmm3,%xmm0
movsd  %xmm0,(%rcx,%rdx,8)
add    $0x1,%rdx
cmp    %rdx,%rax
jne    13b8 <create_gaussian+0x268>
add    $0x58,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

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