#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 sqrt(double);
void levelSetCurvatureCoefficientsEvaluate(int nPoints, int nSpace,
                                           double *grad_phi, double *u,
                                           double *f, double *r, double *dr) {
  int i, I;
  double norm_grad_phi;
  for (i = 0; i < nPoints; i++) {
    r[i] = u[i];
    dr[i] = 1.0;
    norm_grad_phi = 0.0;
    for (I = 0; I < nSpace; I++)
      norm_grad_phi += grad_phi[i * nSpace + I] * grad_phi[i * nSpace + I];
    norm_grad_phi = sqrt(norm_grad_phi);
    if (norm_grad_phi > 0.0) {
      for (I = 0; I < nSpace; I++)
        f[i * nSpace + I] = grad_phi[i * nSpace + I] / norm_grad_phi;
    } else
      f[i * nSpace + I] = 0.0;
  }
}


x86 assembly (clang, -O0):
<levelSetCurvatureCoefficientsEvaluate>:
push   %rbp
mov    %rsp,%rbp
sub    $0x40,%rsp
mov    0x10(%rbp),%rax
mov    %edi,-0x4(%rbp)
mov    %esi,-0x8(%rbp)
mov    %rdx,-0x10(%rbp)
mov    %rcx,-0x18(%rbp)
mov    %r8,-0x20(%rbp)
mov    %r9,-0x28(%rbp)
movl   $0x0,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    127f <levelSetCurvatureCoefficientsEvaluate+0x16f>
mov    -0x18(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x28(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    0x10(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movsd  0xe91(%rip),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x38(%rbp)
movl   $0x0,-0x30(%rbp)
mov    -0x30(%rbp),%eax
cmp    -0x8(%rbp),%eax
jge    11db <levelSetCurvatureCoefficientsEvaluate+0xcb>
mov    -0x10(%rbp),%rax
mov    -0x2c(%rbp),%ecx
imul   -0x8(%rbp),%ecx
add    -0x30(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x2c(%rbp),%ecx
imul   -0x8(%rbp),%ecx
add    -0x30(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm2
movsd  -0x38(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x38(%rbp)
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
jmp    1183 <levelSetCurvatureCoefficientsEvaluate+0x73>
movsd  -0x38(%rbp),%xmm0
call   1030 <sqrt@plt>
movsd  %xmm0,-0x38(%rbp)
movsd  -0x38(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jbe    1253 <levelSetCurvatureCoefficientsEvaluate+0x143>
movl   $0x0,-0x30(%rbp)
mov    -0x30(%rbp),%eax
cmp    -0x8(%rbp),%eax
jge    124e <levelSetCurvatureCoefficientsEvaluate+0x13e>
mov    -0x10(%rbp),%rax
mov    -0x2c(%rbp),%ecx
imul   -0x8(%rbp),%ecx
add    -0x30(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
divsd  -0x38(%rbp),%xmm0
mov    -0x20(%rbp),%rax
mov    -0x2c(%rbp),%ecx
imul   -0x8(%rbp),%ecx
add    -0x30(%rbp),%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
jmp    1203 <levelSetCurvatureCoefficientsEvaluate+0xf3>
jmp    126c <levelSetCurvatureCoefficientsEvaluate+0x15c>
mov    -0x20(%rbp),%rax
mov    -0x2c(%rbp),%ecx
imul   -0x8(%rbp),%ecx
add    -0x30(%rbp),%ecx
movslq %ecx,%rcx
xorps  %xmm0,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
jmp    1271 <levelSetCurvatureCoefficientsEvaluate+0x161>
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    1139 <levelSetCurvatureCoefficientsEvaluate+0x29>
add    $0x40,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<levelSetCurvatureCoefficientsEvaluate>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x98,%rsp
mov    %r9,0x90(%rsp)
mov    %rcx,0x88(%rsp)
mov    %rdx,0x10(%rsp)
mov    %rsi,0x28(%rsp)
test   %edi,%edi
jle    14b1 <levelSetCurvatureCoefficientsEvaluate+0x3a1>
mov    0xd0(%rsp),%r10
mov    0x28(%rsp),%rax
lea    -0x1(%rax),%edx
add    $0x1,%rdx
mov    %rdx,0x78(%rsp)
movslq %eax,%r9
mov    %edi,%ecx
mov    %r9d,%eax
mov    %rax,0x8(%rsp)
lea    -0x1(%rax),%rdx
mov    %rdx,0x70(%rsp)
mov    %eax,%edx
and    $0xfffffffe,%edx
mov    %rdx,0x58(%rsp)
add    $0xfffffffffffffffe,%rdx
mov    %rdx,0x50(%rsp)
mov    %rdx,%rbp
shr    %rbp
add    $0x1,%rbp
mov    %eax,%r12d
and    $0x3,%r12d
mov    %eax,%r14d
and    $0xfffffffc,%r14d
mov    %rbp,%rdx
mov    %rbp,0x48(%rsp)
and    $0xfffffffffffffffe,%rbp
mov    %rbp,0x40(%rsp)
neg    %rax
mov    %rax,0x60(%rsp)
mov    0x10(%rsp),%rbp
lea    0x18(%rbp),%r13
lea    0x0(,%r9,8),%rax
mov    %rax,0x80(%rsp)
xorpd  %xmm3,%xmm3
lea    0x10(%r8),%rdi
lea    0x8(%r8),%r15
xor    %ebx,%ebx
mov    %r9,0x20(%rsp)
mov    %rcx,0x18(%rsp)
mov    %r8,0x68(%rsp)
jmp    1225 <levelSetCurvatureCoefficientsEvaluate+0x115>
xchg   %ax,%ax
xor    %eax,%eax
cltq
add    %rsi,%rax
movq   $0x0,(%r8,%rax,8)
add    $0x1,%rbx
mov    0x80(%rsp),%rax
add    %rax,%r13
add    %rax,%rbp
mov    0x38(%rsp),%rdi
add    %rax,%rdi
add    %rax,%r15
cmp    %rcx,%rbx
je     14b1 <levelSetCurvatureCoefficientsEvaluate+0x3a1>
mov    %rbx,%rsi
imul   %r9,%rsi
mov    0x88(%rsp),%rax
movsd  (%rax,%rbx,8),%xmm0
mov    0x90(%rsp),%rax
movsd  %xmm0,(%rax,%rbx,8)
movabs $0x3ff0000000000000,%rax
mov    %rax,(%r10,%rbx,8)
test   %r9d,%r9d
mov    %rdi,0x38(%rsp)
jle    11f0 <levelSetCurvatureCoefficientsEvaluate+0xe0>
cmpq   $0x3,0x70(%rsp)
jae    1270 <levelSetCurvatureCoefficientsEvaluate+0x160>
xorpd  %xmm0,%xmm0
xor    %eax,%eax
jmp    12c5 <levelSetCurvatureCoefficientsEvaluate+0x1b5>
xchg   %ax,%ax
xorpd  %xmm0,%xmm0
xor    %eax,%eax
cs nopw 0x0(%rax,%rax,1)
movsd  -0x18(%r13,%rax,8),%xmm1
movsd  -0x10(%r13,%rax,8),%xmm2
mulsd  %xmm1,%xmm1
addsd  %xmm0,%xmm1
mulsd  %xmm2,%xmm2
addsd  %xmm1,%xmm2
movsd  -0x8(%r13,%rax,8),%xmm1
mulsd  %xmm1,%xmm1
addsd  %xmm2,%xmm1
movsd  0x0(%r13,%rax,8),%xmm0
mulsd  %xmm0,%xmm0
addsd  %xmm1,%xmm0
add    $0x4,%rax
cmp    %rax,%r14
jne    1280 <levelSetCurvatureCoefficientsEvaluate+0x170>
test   %r12,%r12
je     12f6 <levelSetCurvatureCoefficientsEvaluate+0x1e6>
lea    0x0(,%rax,8),%rax
add    %rbp,%rax
xor    %edx,%edx
nopw   0x0(%rax,%rax,1)
movsd  (%rax,%rdx,8),%xmm1
mulsd  %xmm1,%xmm1
addsd  %xmm1,%xmm0
add    $0x1,%rdx
cmp    %rdx,%r12
jne    12e0 <levelSetCurvatureCoefficientsEvaluate+0x1d0>
ucomisd %xmm3,%xmm0
mov    %rsi,0x30(%rsp)
jb     1310 <levelSetCurvatureCoefficientsEvaluate+0x200>
sqrtsd %xmm0,%xmm0
jmp    1335 <levelSetCurvatureCoefficientsEvaluate+0x225>
nopw   0x0(%rax,%rax,1)
call   1030 <sqrt@plt>
mov    0x30(%rsp),%rsi
xorpd  %xmm3,%xmm3
mov    0x18(%rsp),%rcx
mov    0x20(%rsp),%r9
mov    0xd0(%rsp),%r10
mov    0x68(%rsp),%r8
ucomisd %xmm3,%xmm0
mov    0x78(%rsp),%rax
jbe    11f2 <levelSetCurvatureCoefficientsEvaluate+0xe2>
cmpl   $0x0,0x28(%rsp)
jle    11ff <levelSetCurvatureCoefficientsEvaluate+0xef>
cmpl   $0x1,0x28(%rsp)
jne    1362 <levelSetCurvatureCoefficientsEvaluate+0x252>
xor    %edx,%edx
mov    0x8(%rsp),%r11
jmp    142e <levelSetCurvatureCoefficientsEvaluate+0x31e>
mov    0x8(%rsp),%rax
mov    0x30(%rsp),%rcx
add    %rcx,%rax
lea    (%r8,%rcx,8),%rsi
mov    0x10(%rsp),%r11
lea    (%r11,%rax,8),%rdx
cmp    %rdx,%rsi
jae    13a4 <levelSetCurvatureCoefficientsEvaluate+0x294>
lea    (%r8,%rax,8),%rax
lea    (%r11,%rcx,8),%rdx
cmp    %rax,%rdx
jae    13a4 <levelSetCurvatureCoefficientsEvaluate+0x294>
xor    %edx,%edx
mov    0x20(%rsp),%r9
mov    0x18(%rsp),%rcx
mov    0x8(%rsp),%r11
jmp    142e <levelSetCurvatureCoefficientsEvaluate+0x31e>
movapd %xmm0,%xmm1
unpcklpd %xmm0,%xmm1
cmpq   $0x0,0x50(%rsp)
mov    0x20(%rsp),%r9
je     149f <levelSetCurvatureCoefficientsEvaluate+0x38f>
mov    0x40(%rsp),%rdx
xor    %eax,%eax
mov    0x38(%rsp),%rdi
nopl   0x0(%rax)
movupd 0x0(%rbp,%rax,8),%xmm2
divpd  %xmm1,%xmm2
movupd %xmm2,-0x10(%rdi,%rax,8)
movupd 0x10(%rbp,%rax,8),%xmm2
divpd  %xmm1,%xmm2
movupd %xmm2,(%rdi,%rax,8)
add    $0x4,%rax
add    $0xfffffffffffffffe,%rdx
jne    13d0 <levelSetCurvatureCoefficientsEvaluate+0x2c0>
testb  $0x1,0x48(%rsp)
je     1413 <levelSetCurvatureCoefficientsEvaluate+0x303>
add    %rcx,%rax
movupd (%r11,%rax,8),%xmm2
divpd  %xmm1,%xmm2
movupd %xmm2,(%r8,%rax,8)
mov    0x58(%rsp),%rax
mov    %rax,%rdx
mov    0x8(%rsp),%r11
cmp    %r11,%rax
mov    0x18(%rsp),%rcx
je     11ff <levelSetCurvatureCoefficientsEvaluate+0xef>
mov    %rdx,%rax
test   $0x1,%r11b
je     145a <levelSetCurvatureCoefficientsEvaluate+0x34a>
mov    0x30(%rsp),%rsi
add    %rdx,%rsi
mov    0x10(%rsp),%rax
movsd  (%rax,%rsi,8),%xmm1
divsd  %xmm0,%xmm1
movsd  %xmm1,(%r8,%rsi,8)
mov    %rdx,%rax
or     $0x1,%rax
not    %rdx
cmp    0x60(%rsp),%rdx
mov    0x8(%rsp),%rdx
je     11ff <levelSetCurvatureCoefficientsEvaluate+0xef>
nopl   (%rax)
movsd  0x0(%rbp,%rax,8),%xmm1
divsd  %xmm0,%xmm1
movsd  %xmm1,-0x8(%r15,%rax,8)
movsd  0x8(%rbp,%rax,8),%xmm1
divsd  %xmm0,%xmm1
movsd  %xmm1,(%r15,%rax,8)
add    $0x2,%rax
cmp    %rax,%rdx
jne    1470 <levelSetCurvatureCoefficientsEvaluate+0x360>
jmp    11ff <levelSetCurvatureCoefficientsEvaluate+0xef>
xor    %eax,%eax
testb  $0x1,0x48(%rsp)
jne    1400 <levelSetCurvatureCoefficientsEvaluate+0x2f0>
jmp    1413 <levelSetCurvatureCoefficientsEvaluate+0x303>
add    $0x98,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

LLVM IR (-O2):
define void @levelSetCurvatureCoefficientsEvaluate(i32 %0, i32 %1, double* nocapture %2, double* nocapture %3, double* nocapture writeonly %4, double* nocapture writeonly %5, double* nocapture writeonly %6) {
%8 = icmp sgt i32 %0, 0
br i1 %8, label %9, label %176
9: ; preds = %7
%10 = icmp sgt i32 %1, 0
%11 = add i32 %1, -1
%12 = zext i32 %11 to i64
%13 = add nuw nsw i64 %12, 1
%14 = sext i32 %1 to i64
%15 = zext i32 %0 to i64
%16 = zext i32 %1 to i64
%17 = zext i32 %1 to i64
%18 = add nsw i64 %16, -1
%19 = and i64 %16, 4294967294
%20 = add nsw i64 %19, -2
%21 = lshr exact i64 %20, 1
%22 = add nuw i64 %21, 1
%23 = and i64 %16, 3
%24 = icmp ult i64 %18, 3
%25 = and i64 %16, 4294967292
%26 = icmp eq i64 %23, 0
%27 = icmp eq i32 %1, 1
%28 = and i64 %16, 4294967294
%29 = and i64 %22, 1
%30 = icmp eq i64 %20, 0
%31 = and i64 %22, -2
%32 = icmp eq i64 %29, 0
%33 = icmp eq i64 %28, %16
%34 = and i64 %16, 1
%35 = icmp eq i64 %34, 0
%36 = sub nsw i64 0, %16
br label %37
37: ; preds = %9, %173
%38 = phi i64 [ 0, %9 ], [ %174, %173 ]
%39 = mul i64 %38, %14
%40 = getelementptr double, double* %4, i64 %39
%41 = add i64 %39, %16
%42 = getelementptr double, double* %4, i64 %41
%43 = getelementptr double, double* %2, i64 %39
%44 = getelementptr double, double* %2, i64 %41
%45 = getelementptr inbounds double, double* %3, i64 %38
%46 = load double, double* %45
%47 = getelementptr inbounds double, double* %5, i64 %38
store double %46, double* %47
%48 = getelementptr inbounds double, double* %6, i64 %38
store double 1.000000e+00, double* %48
br i1 %10, label %49, label %166
49: ; preds = %37
%50 = mul nsw i64 %38, %14
br i1 %24, label %77, label %51
51: ; preds = %49, %51
%52 = phi i64 [ %74, %51 ], [ 0, %49 ]
%53 = phi double [ %73, %51 ], [ 0.000000e+00, %49 ]
%54 = phi i64 [ %75, %51 ], [ 0, %49 ]
%55 = add nsw i64 %52, %50
%56 = getelementptr inbounds double, double* %2, i64 %55
%57 = load double, double* %56
%58 = tail call double @llvm.fmuladd.f64(double %57, double %57, double %53)
%59 = or i64 %52, 1
%60 = add nsw i64 %59, %50
%61 = getelementptr inbounds double, double* %2, i64 %60
%62 = load double, double* %61
%63 = tail call double @llvm.fmuladd.f64(double %62, double %62, double %58)
%64 = or i64 %52, 2
%65 = add nsw i64 %64, %50
%66 = getelementptr inbounds double, double* %2, i64 %65
%67 = load double, double* %66
%68 = tail call double @llvm.fmuladd.f64(double %67, double %67, double %63)
%69 = or i64 %52, 3
%70 = add nsw i64 %69, %50
%71 = getelementptr inbounds double, double* %2, i64 %70
%72 = load double, double* %71
%73 = tail call double @llvm.fmuladd.f64(double %72, double %72, double %68)
%74 = add nuw nsw i64 %52, 4
%75 = add i64 %54, 4
%76 = icmp eq i64 %75, %25
br i1 %76, label %77, label %51
77: ; preds = %51, %49
%78 = phi double [ undef, %49 ], [ %73, %51 ]
%79 = phi i64 [ 0, %49 ], [ %74, %51 ]
%80 = phi double [ 0.000000e+00, %49 ], [ %73, %51 ]
br i1 %26, label %92, label %81
81: ; preds = %77, %81
%82 = phi i64 [ %89, %81 ], [ %79, %77 ]
%83 = phi double [ %88, %81 ], [ %80, %77 ]
%84 = phi i64 [ %90, %81 ], [ 0, %77 ]
%85 = add nsw i64 %82, %50
%86 = getelementptr inbounds double, double* %2, i64 %85
%87 = load double, double* %86
%88 = tail call double @llvm.fmuladd.f64(double %87, double %87, double %83)
%89 = add nuw nsw i64 %82, 1
%90 = add i64 %84, 1
%91 = icmp eq i64 %90, %23
br i1 %91, label %92, label %81
92: ; preds = %81, %77
%93 = phi double [ %78, %77 ], [ %88, %81 ]
%94 = tail call double @sqrt(double %93)
%95 = fcmp ogt double %94, 0.000000e+00
br i1 %95, label %96, label %166
96: ; preds = %92
br i1 %10, label %97, label %173
97: ; preds = %96
%98 = mul nsw i64 %38, %14
br i1 %27, label %138, label %99
99: ; preds = %97
%100 = icmp ult double* %40, %44
%101 = icmp ult double* %43, %42
%102 = and i1 %100, %101
br i1 %102, label %138, label %103
103: ; preds = %99
%104 = insertelement <2 x double> poison, double %94, i64 0
%105 = shufflevector <2 x double> %104, <2 x double> poison, <2 x i32> zeroinitializer
br i1 %30, label %127, label %106
106: ; preds = %103, %106
%107 = phi i64 [ %124, %106 ], [ 0, %103 ]
%108 = phi i64 [ %125, %106 ], [ 0, %103 ]
%109 = add nsw i64 %107, %98
%110 = getelementptr inbounds double, double* %2, i64 %109
%111 = bitcast double* %110 to <2 x double>*
%112 = load <2 x double>, <2 x double>* %111
%113 = fdiv <2 x double> %112, %105
%114 = getelementptr inbounds double, double* %4, i64 %109
%115 = bitcast double* %114 to <2 x double>*
store <2 x double> %113, <2 x double>* %115
%116 = or i64 %107, 2
%117 = add nsw i64 %116, %98
%118 = getelementptr inbounds double, double* %2, i64 %117
%119 = bitcast double* %118 to <2 x double>*
%120 = load <2 x double>, <2 x double>* %119
%121 = fdiv <2 x double> %120, %105
%122 = getelementptr inbounds double, double* %4, i64 %117
%123 = bitcast double* %122 to <2 x double>*
store <2 x double> %121, <2 x double>* %123
%124 = add nuw i64 %107, 4
%125 = add nuw i64 %108, 2
%126 = icmp eq i64 %125, %31
br i1 %126, label %127, label %106
127: ; preds = %106, %103
%128 = phi i64 [ 0, %103 ], [ %124, %106 ]
br i1 %32, label %137, label %129
129: ; preds = %127
%130 = add nsw i64 %128, %98
%131 = getelementptr inbounds double, double* %2, i64 %130
%132 = bitcast double* %131 to <2 x double>*
%133 = load <2 x double>, <2 x double>* %132
%134 = fdiv <2 x double> %133, %105
%135 = getelementptr inbounds double, double* %4, i64 %130
%136 = bitcast double* %135 to <2 x double>*
store <2 x double> %134, <2 x double>* %136
br label %137
137: ; preds = %127, %129
br i1 %33, label %173, label %138
138: ; preds = %99, %97, %137
%139 = phi i64 [ 0, %99 ], [ 0, %97 ], [ %28, %137 ]
%140 = xor i64 %139, -1
br i1 %35, label %148, label %141
141: ; preds = %138
%142 = add nsw i64 %139, %98
%143 = getelementptr inbounds double, double* %2, i64 %142
%144 = load double, double* %143
%145 = fdiv double %144, %94
%146 = getelementptr inbounds double, double* %4, i64 %142
store double %145, double* %146
%147 = or i64 %139, 1
br label %148
148: ; preds = %141, %138
%149 = phi i64 [ %139, %138 ], [ %147, %141 ]
%150 = icmp eq i64 %140, %36
br i1 %150, label %173, label %151
151: ; preds = %148, %151
%152 = phi i64 [ %164, %151 ], [ %149, %148 ]
%153 = add nsw i64 %152, %98
%154 = getelementptr inbounds double, double* %2, i64 %153
%155 = load double, double* %154
%156 = fdiv double %155, %94
%157 = getelementptr inbounds double, double* %4, i64 %153
store double %156, double* %157
%158 = add nuw nsw i64 %152, 1
%159 = add nsw i64 %158, %98
%160 = getelementptr inbounds double, double* %2, i64 %159
%161 = load double, double* %160
%162 = fdiv double %161, %94
%163 = getelementptr inbounds double, double* %4, i64 %159
store double %162, double* %163
%164 = add nuw nsw i64 %152, 2
%165 = icmp eq i64 %164, %17
br i1 %165, label %173, label %151
166: ; preds = %37, %92
%167 = phi i64 [ %13, %92 ], [ 0, %37 ]
%168 = mul nsw i64 %38, %14
%169 = shl i64 %167, 32
%170 = ashr exact i64 %169, 32
%171 = add nsw i64 %170, %168
%172 = getelementptr inbounds double, double* %4, i64 %171
store double 0.000000e+00, double* %172
br label %173
173: ; preds = %148, %151, %137, %96, %166
%174 = add nuw nsw i64 %38, 1
%175 = icmp eq i64 %174, %15
br i1 %175, label %176, label %37
176: ; preds = %173, %7
ret void
}