#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 constantNormalVelocityLevelSetEvaluate(const int nPoints, const int nSpace,
                                            double b, const double *x,
                                            const double *u,
                                            const double *gradu, double *m,
                                            double *dm, double *f, double *df,
                                            double *H, double *dH) {
  double normgradu;
  int k, id;
  for (k = 0; k < nPoints; k++) {
    m[k] = u[k];
    dm[k] = 1.0;
    normgradu = 0.0;
    for (id = 0; id < nSpace; id++) {
      f[k * nSpace + id] = 0.0;
      df[k * nSpace + id] = 0.0;
      normgradu += gradu[k * nSpace + id] * gradu[k * nSpace + id];
    }
    normgradu = sqrt(normgradu);
    H[k] = b * normgradu;
    for (id = 0; id < nSpace; id++) {
      dH[k * nSpace + id] = gradu[k * nSpace + id] / (normgradu + 1.0e-8);
    }
  }
}


x86 assembly (clang, -O0):
<constantNormalVelocityLevelSetEvaluate>:
push   %rbp
mov    %rsp,%rbp
sub    $0x40,%rsp
mov    0x30(%rbp),%rax
mov    0x28(%rbp),%rax
mov    0x20(%rbp),%rax
mov    0x18(%rbp),%rax
mov    0x10(%rbp),%rax
mov    %edi,-0x4(%rbp)
mov    %esi,-0x8(%rbp)
movsd  %xmm0,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %rcx,-0x20(%rbp)
mov    %r8,-0x28(%rbp)
mov    %r9,-0x30(%rbp)
movl   $0x0,-0x3c(%rbp)
mov    -0x3c(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    12b9 <constantNormalVelocityLevelSetEvaluate+0x1a9>
mov    -0x20(%rbp),%rax
movslq -0x3c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x30(%rbp),%rax
movslq -0x3c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    0x10(%rbp),%rax
movslq -0x3c(%rbp),%rcx
movsd  0xe7c(%rip),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x38(%rbp)
movl   $0x0,-0x40(%rbp)
mov    -0x40(%rbp),%eax
cmp    -0x8(%rbp),%eax
jge    1222 <constantNormalVelocityLevelSetEvaluate+0x112>
mov    0x18(%rbp),%rax
mov    -0x3c(%rbp),%ecx
imul   -0x8(%rbp),%ecx
add    -0x40(%rbp),%ecx
movslq %ecx,%rcx
xorps  %xmm0,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    0x20(%rbp),%rax
mov    -0x3c(%rbp),%ecx
imul   -0x8(%rbp),%ecx
add    -0x40(%rbp),%ecx
movslq %ecx,%rcx
xorps  %xmm0,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x28(%rbp),%rax
mov    -0x3c(%rbp),%ecx
imul   -0x8(%rbp),%ecx
add    -0x40(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x28(%rbp),%rax
mov    -0x3c(%rbp),%ecx
imul   -0x8(%rbp),%ecx
add    -0x40(%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    -0x40(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x40(%rbp)
jmp    1198 <constantNormalVelocityLevelSetEvaluate+0x88>
movsd  -0x38(%rbp),%xmm0
call   1030 <sqrt@plt>
movsd  %xmm0,-0x38(%rbp)
movsd  -0x10(%rbp),%xmm0
mulsd  -0x38(%rbp),%xmm0
mov    0x28(%rbp),%rax
movslq -0x3c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
movl   $0x0,-0x40(%rbp)
mov    -0x40(%rbp),%eax
cmp    -0x8(%rbp),%eax
jge    12a6 <constantNormalVelocityLevelSetEvaluate+0x196>
mov    -0x28(%rbp),%rax
mov    -0x3c(%rbp),%ecx
imul   -0x8(%rbp),%ecx
add    -0x40(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  0xd8f(%rip),%xmm1
addsd  -0x38(%rbp),%xmm1
divsd  %xmm1,%xmm0
mov    0x30(%rbp),%rax
mov    -0x3c(%rbp),%ecx
imul   -0x8(%rbp),%ecx
add    -0x40(%rbp),%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x40(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x40(%rbp)
jmp    124f <constantNormalVelocityLevelSetEvaluate+0x13f>
jmp    12ab <constantNormalVelocityLevelSetEvaluate+0x19b>
mov    -0x3c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x3c(%rbp)
jmp    114e <constantNormalVelocityLevelSetEvaluate+0x3e>
add    $0x40,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<constantNormalVelocityLevelSetEvaluate>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x88,%rsp
mov    %r9,0x80(%rsp)
mov    %r8,0x8(%rsp)
test   %edi,%edi
jle    1532 <constantNormalVelocityLevelSetEvaluate+0x422>
movapd %xmm0,%xmm3
mov    0xe0(%rsp),%rax
mov    0xd8(%rsp),%r9
mov    0xd0(%rsp),%r10
mov    0xc8(%rsp),%r14
movslq %esi,%rbx
mov    %edi,%r8d
xorpd  %xmm4,%xmm4
mulsd  %xmm4,%xmm0
movsd  %xmm0,0x68(%rsp)
mov    %ebx,%edi
mov    %edi,%edx
and    $0xfffffffe,%edx
mov    %rdx,0x10(%rsp)
add    $0xfffffffffffffffe,%rdx
mov    %rdx,0x48(%rsp)
mov    %rdx,%rbp
shr    %rbp
add    $0x1,%rbp
mov    %rbp,%rdx
mov    %rbp,0x40(%rsp)
and    $0xfffffffffffffffe,%rbp
mov    %rbp,0x38(%rsp)
mov    %rdi,%rbp
mov    %rdi,(%rsp)
neg    %rdi
mov    %rdi,0x60(%rsp)
mov    0x8(%rsp),%rdx
lea    0x8(%rdx),%r15
mov    %rbx,0x70(%rsp)
lea    0x0(,%rbx,8),%rdi
mov    %rdi,0x78(%rsp)
lea    0x8(%r10),%r12
lea    0x8(%r14),%rbp
lea    0x10(%rax),%rdi
lea    0x10(%rdx),%r14
lea    0x8(%rax),%r13
xor    %ebx,%ebx
movabs $0x3ff0000000000000,%r11
mov    %r8,0x18(%rsp)
mov    %rcx,0x58(%rsp)
movsd  %xmm3,0x50(%rsp)
mov    %esi,0x2c(%rsp)
jmp    1230 <constantNormalVelocityLevelSetEvaluate+0x120>
xchg   %ax,%ax
movsd  0x68(%rsp),%xmm0
movsd  %xmm0,(%r9,%rbx,8)
add    $0x1,%rbx
mov    0x78(%rsp),%rax
add    %rax,%r15
add    %rax,%r12
add    %rax,%rbp
add    %rax,%rdi
add    %rax,%r14
add    %rax,%r13
cmp    %r8,%rbx
je     1532 <constantNormalVelocityLevelSetEvaluate+0x422>
movsd  (%rcx,%rbx,8),%xmm0
mov    0x80(%rsp),%rax
movsd  %xmm0,(%rax,%rbx,8)
mov    0xc0(%rsp),%rax
mov    %r11,(%rax,%rbx,8)
test   %esi,%esi
jle    1200 <constantNormalVelocityLevelSetEvaluate+0xf0>
cmp    $0x1,%esi
mov    %rdi,0x30(%rsp)
xorpd  %xmm0,%xmm0
jne    1270 <constantNormalVelocityLevelSetEvaluate+0x160>
xor    %eax,%eax
mov    (%rsp),%rdx
jmp    12c9 <constantNormalVelocityLevelSetEvaluate+0x1b9>
nopl   0x0(%rax,%rax,1)
xor    %eax,%eax
mov    (%rsp),%rdx
mov    0x10(%rsp),%rdi
nopl   0x0(%rax,%rax,1)
movq   $0x0,-0x8(%rbp,%rax,8)
movq   $0x0,-0x8(%r12,%rax,8)
movsd  -0x8(%r15,%rax,8),%xmm1
mulsd  %xmm1,%xmm1
addsd  %xmm0,%xmm1
movq   $0x0,0x0(%rbp,%rax,8)
movq   $0x0,(%r12,%rax,8)
movsd  (%r15,%rax,8),%xmm0
mulsd  %xmm0,%xmm0
addsd  %xmm1,%xmm0
add    $0x2,%rax
cmp    %rax,%rdi
jne    1280 <constantNormalVelocityLevelSetEvaluate+0x170>
mov    %rbx,%rdi
imul   0x70(%rsp),%rdi
test   $0x1,%dl
je     130c <constantNormalVelocityLevelSetEvaluate+0x1fc>
add    %rdi,%rax
mov    0xc8(%rsp),%rdx
movq   $0x0,(%rdx,%rax,8)
mov    0xd0(%rsp),%rdx
movq   $0x0,(%rdx,%rax,8)
mov    0x8(%rsp),%rdx
movsd  (%rdx,%rax,8),%xmm1
mulsd  %xmm1,%xmm1
addsd  %xmm1,%xmm0
ucomisd %xmm4,%xmm0
mov    %rdi,0x20(%rsp)
jb     1320 <constantNormalVelocityLevelSetEvaluate+0x210>
sqrtsd %xmm0,%xmm0
jmp    134f <constantNormalVelocityLevelSetEvaluate+0x23f>
nopl   (%rax)
call   1030 <sqrt@plt>
xorpd  %xmm4,%xmm4
mov    0x18(%rsp),%r8
mov    0xd8(%rsp),%r9
movabs $0x3ff0000000000000,%r11
mov    0x2c(%rsp),%esi
movsd  0x50(%rsp),%xmm3
mov    0x58(%rsp),%rcx
mov    0x30(%rsp),%rdi
movapd %xmm0,%xmm1
mulsd  %xmm3,%xmm1
movsd  %xmm1,(%r9,%rbx,8)
test   %esi,%esi
jle    120c <constantNormalVelocityLevelSetEvaluate+0xfc>
addsd  0xc8e(%rip),%xmm0
cmp    $0x1,%esi
jne    1382 <constantNormalVelocityLevelSetEvaluate+0x272>
xor    %edx,%edx
mov    (%rsp),%r10
jmp    1490 <constantNormalVelocityLevelSetEvaluate+0x380>
mov    %rdi,%r11
mov    (%rsp),%rax
mov    0x20(%rsp),%rdi
add    %rdi,%rax
mov    0xe0(%rsp),%rdx
lea    (%rdx,%rdi,8),%r8
mov    0x8(%rsp),%rdx
lea    (%rdx,%rax,8),%rdx
cmp    %rdx,%r8
jae    13e2 <constantNormalVelocityLevelSetEvaluate+0x2d2>
mov    0xe0(%rsp),%rdx
lea    (%rdx,%rax,8),%rax
mov    0x8(%rsp),%rdx
lea    (%rdx,%rdi,8),%rdx
cmp    %rax,%rdx
jae    13e2 <constantNormalVelocityLevelSetEvaluate+0x2d2>
xor    %edx,%edx
mov    0x18(%rsp),%r8
mov    (%rsp),%r10
mov    %r11,%rdi
movabs $0x3ff0000000000000,%r11
jmp    1490 <constantNormalVelocityLevelSetEvaluate+0x380>
movapd %xmm0,%xmm1
unpcklpd %xmm0,%xmm1
cmpq   $0x0,0x48(%rsp)
mov    0x18(%rsp),%r8
mov    %r11,%rdi
je     1511 <constantNormalVelocityLevelSetEvaluate+0x401>
mov    0x38(%rsp),%rdx
xor    %eax,%eax
mov    0x10(%rsp),%r10
movabs $0x3ff0000000000000,%r11
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
movupd -0x10(%r14,%rax,8),%xmm2
divpd  %xmm1,%xmm2
movupd %xmm2,-0x10(%rdi,%rax,8)
movupd (%r14,%rax,8),%xmm2
divpd  %xmm1,%xmm2
movupd %xmm2,(%rdi,%rax,8)
add    $0x4,%rax
add    $0xfffffffffffffffe,%rdx
jne    1420 <constantNormalVelocityLevelSetEvaluate+0x310>
testb  $0x1,0x40(%rsp)
je     1471 <constantNormalVelocityLevelSetEvaluate+0x361>
add    0x20(%rsp),%rax
mov    0x8(%rsp),%rdx
movupd (%rdx,%rax,8),%xmm2
divpd  %xmm1,%xmm2
mov    0xe0(%rsp),%rdx
movupd %xmm2,(%rdx,%rax,8)
mov    %r10,%rdx
mov    %r10,%rax
mov    (%rsp),%r10
cmp    %r10,%rax
je     120c <constantNormalVelocityLevelSetEvaluate+0xfc>
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
mov    %rdx,%rax
test   $0x1,%r10b
je     14c5 <constantNormalVelocityLevelSetEvaluate+0x3b5>
mov    0x20(%rsp),%r10
add    %rdx,%r10
mov    0x8(%rsp),%rax
movsd  (%rax,%r10,8),%xmm1
divsd  %xmm0,%xmm1
mov    0xe0(%rsp),%rax
movsd  %xmm1,(%rax,%r10,8)
mov    %rdx,%rax
or     $0x1,%rax
not    %rdx
cmp    0x60(%rsp),%rdx
mov    (%rsp),%rdx
je     120c <constantNormalVelocityLevelSetEvaluate+0xfc>
nopw   0x0(%rax,%rax,1)
movsd  -0x8(%r15,%rax,8),%xmm1
divsd  %xmm0,%xmm1
movsd  %xmm1,-0x8(%r13,%rax,8)
movsd  (%r15,%rax,8),%xmm1
divsd  %xmm0,%xmm1
movsd  %xmm1,0x0(%r13,%rax,8)
add    $0x2,%rax
cmp    %rax,%rdx
jne    14e0 <constantNormalVelocityLevelSetEvaluate+0x3d0>
jmp    120c <constantNormalVelocityLevelSetEvaluate+0xfc>
xor    %eax,%eax
mov    0x10(%rsp),%r10
movabs $0x3ff0000000000000,%r11
testb  $0x1,0x40(%rsp)
jne    1451 <constantNormalVelocityLevelSetEvaluate+0x341>
jmp    1471 <constantNormalVelocityLevelSetEvaluate+0x361>
add    $0x88,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

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