#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 sqrtf(float);
void computegradient(const float *img, int w, int h, float *gx, float *gy) {
  int i, j, k, p, q;
  float glength, phi, phiscaled, ascaled, errsign, pfrac, qfrac, err0, err1,
      err;
  for (i = 1; i < h - 1; i++) {
    for (j = 1; j < w - 1; j++) {
      k = i * w + j;
      if ((img[k] > 0) && (img[k] < 1)) {
        gx[k] = -img[k - w - 1] - 1.4142136f * img[k - 1] - img[k + w - 1] +
                img[k - w + 1] + 1.4142136f * img[k + 1] + img[k + w + 1];
        gy[k] = -img[k - w - 1] - 1.4142136f * img[k - w] - img[k - w + 1] +
                img[k + w - 1] + 1.4142136f * img[k + w] + img[k + w + 1];
        glength = gx[k] * gx[k] + gy[k] * gy[k];
        if (glength > 0) {
          glength = sqrtf(glength);
          gx[k] = gx[k] / glength;
          gy[k] = gy[k] / glength;
        }
      }
    }
  }
}


x86 assembly (clang, -O0):
<computegradient>:
push   %rbp
mov    %rsp,%rbp
sub    $0x60,%rsp
mov    %rdi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
mov    %edx,-0x10(%rbp)
mov    %rcx,-0x18(%rbp)
mov    %r8,-0x20(%rbp)
movl   $0x1,-0x24(%rbp)
mov    -0x24(%rbp),%eax
mov    -0x10(%rbp),%ecx
sub    $0x1,%ecx
cmp    %ecx,%eax
jge    13c8 <computegradient+0x2b8>
movl   $0x1,-0x28(%rbp)
mov    -0x28(%rbp),%eax
mov    -0xc(%rbp),%ecx
sub    $0x1,%ecx
cmp    %ecx,%eax
jge    13b5 <computegradient+0x2a5>
mov    -0x24(%rbp),%eax
imul   -0xc(%rbp),%eax
add    -0x28(%rbp),%eax
mov    %eax,-0x2c(%rbp)
mov    -0x8(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
xorps  %xmm1,%xmm1
ucomiss %xmm1,%xmm0
jbe    13a2 <computegradient+0x292>
mov    -0x8(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movss  0xe70(%rip),%xmm0
ucomiss (%rax,%rcx,4),%xmm0
jbe    13a2 <computegradient+0x292>
mov    -0x8(%rbp),%rax
mov    -0x2c(%rbp),%ecx
sub    -0xc(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm0
movd   %xmm0,%eax
xor    $0x80000000,%eax
movd   %eax,%xmm0
mov    -0x8(%rbp),%rax
mov    -0x2c(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm1
movss  0xe32(%rip),%xmm2
mulss  %xmm2,%xmm1
addss  %xmm0,%xmm1
mov    -0x8(%rbp),%rax
mov    -0x2c(%rbp),%ecx
add    -0xc(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
subss  (%rax,%rcx,4),%xmm1
mov    -0x8(%rbp),%rax
mov    -0x2c(%rbp),%ecx
sub    -0xc(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
addss  (%rax,%rcx,4),%xmm1
mov    -0x8(%rbp),%rax
mov    -0x2c(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm0
movss  0xde2(%rip),%xmm2
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
mov    -0x8(%rbp),%rax
mov    -0x2c(%rbp),%ecx
add    -0xc(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
addss  (%rax,%rcx,4),%xmm0
mov    -0x18(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x8(%rbp),%rax
mov    -0x2c(%rbp),%ecx
sub    -0xc(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm0
movd   %xmm0,%eax
xor    $0x80000000,%eax
movd   %eax,%xmm0
mov    -0x8(%rbp),%rax
mov    -0x2c(%rbp),%ecx
sub    -0xc(%rbp),%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm1
movss  0xd80(%rip),%xmm2
mulss  %xmm2,%xmm1
addss  %xmm0,%xmm1
mov    -0x8(%rbp),%rax
mov    -0x2c(%rbp),%ecx
sub    -0xc(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
subss  (%rax,%rcx,4),%xmm1
mov    -0x8(%rbp),%rax
mov    -0x2c(%rbp),%ecx
add    -0xc(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
addss  (%rax,%rcx,4),%xmm1
mov    -0x8(%rbp),%rax
mov    -0x2c(%rbp),%ecx
add    -0xc(%rbp),%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm0
movss  0xd30(%rip),%xmm2
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
mov    -0x8(%rbp),%rax
mov    -0x2c(%rbp),%ecx
add    -0xc(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
addss  (%rax,%rcx,4),%xmm0
mov    -0x20(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x18(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x18(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm2
mov    -0x20(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm1
mov    -0x20(%rbp),%rax
movslq -0x2c(%rbp),%rcx
mulss  (%rax,%rcx,4),%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,-0x38(%rbp)
movss  -0x38(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomiss %xmm1,%xmm0
jbe    139d <computegradient+0x28d>
movss  -0x38(%rbp),%xmm0
call   1030 <sqrtf@plt>
movss  %xmm0,-0x38(%rbp)
mov    -0x18(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
divss  -0x38(%rbp),%xmm0
mov    -0x18(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x20(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
divss  -0x38(%rbp),%xmm0
mov    -0x20(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
jmp    13a2 <computegradient+0x292>
jmp    13a7 <computegradient+0x297>
mov    -0x28(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x28(%rbp)
jmp    1149 <computegradient+0x39>
jmp    13ba <computegradient+0x2aa>
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    1131 <computegradient+0x21>
add    $0x60,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<computegradient>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x28,%rsp
mov    %rsi,0x8(%rsp)
mov    %rdx,(%rsp)
cmp    $0x3,%edx
jl     12ec <computegradient+0x1dc>
mov    %rdi,%rbx
mov    (%rsp),%rax
add    $0xffffffff,%eax
mov    %rax,(%rsp)
mov    0x8(%rsp),%rdx
lea    -0x1(%rdx),%eax
movslq %edx,%rdx
lea    0x8(%rdi,%rdx,4),%r12
lea    0x0(,%rdx,4),%rsi
mov    %rsi,0x10(%rsp)
lea    (%rcx,%rdx,4),%rbp
add    $0x4,%rbp
lea    (%r8,%rdx,4),%r14
add    $0x4,%r14
lea    (%rdi,%rdx,8),%r15
add    $0x8,%r15
add    $0x8,%rbx
lea    -0x4(,%rax,4),%rcx
ff
mov    $0x1,%edx
xorps  %xmm2,%xmm2
movss  0xe71(%rip),%xmm3
movss  0xe6d(%rip),%xmm4
movss  0xe69(%rip),%xmm5
mov    %rcx,0x20(%rsp)
jmp    11d7 <computegradient+0xc7>
cs nopw 0x0(%rax,%rax,1)
mov    0x18(%rsp),%rdx
add    $0x1,%rdx
mov    0x10(%rsp),%rax
add    %rax,%r12
add    %rax,%rbp
add    %rax,%r14
add    %rax,%rbx
add    %rax,%r15
cmp    (%rsp),%rdx
je     12ec <computegradient+0x1dc>
mov    %rdx,0x18(%rsp)
cmpl   $0x3,0x8(%rsp)
jl     11b0 <computegradient+0xa0>
xor    %r13d,%r13d
jmp    123c <computegradient+0x12c>
call   1030 <sqrtf@plt>
movss  0xe13(%rip),%xmm5
movss  0xe07(%rip),%xmm4
movss  0xdfb(%rip),%xmm3
xorps  %xmm2,%xmm2
mov    0x20(%rsp),%rcx
movss  0x0(%rbp,%r13,1),%xmm1
divss  %xmm0,%xmm1
movss  %xmm1,0x0(%rbp,%r13,1)
movss  (%r14,%r13,1),%xmm1
divss  %xmm0,%xmm1
movss  %xmm1,(%r14,%r13,1)
add    $0x4,%r13
cmp    %r13,%rcx
je     11b0 <computegradient+0xa0>
movss  -0x4(%r12,%r13,1),%xmm0
ucomiss %xmm2,%xmm0
jbe    122f <computegradient+0x11f>
ucomiss %xmm0,%xmm3
jbe    122f <computegradient+0x11f>
movss  -0x8(%r12,%r13,1),%xmm0
mulss  %xmm4,%xmm0
subss  -0x8(%rbx,%r13,1),%xmm0
subss  -0x8(%r15,%r13,1),%xmm0
addss  (%rbx,%r13,1),%xmm0
movss  (%r12,%r13,1),%xmm1
mulss  %xmm5,%xmm1
addss  %xmm0,%xmm1
addss  (%r15,%r13,1),%xmm1
movss  %xmm1,0x0(%rbp,%r13,1)
movss  -0x4(%rbx,%r13,1),%xmm0
mulss  %xmm4,%xmm0
subss  -0x8(%rbx,%r13,1),%xmm0
subss  (%rbx,%r13,1),%xmm0
addss  -0x8(%r15,%r13,1),%xmm0
movss  -0x4(%r15,%r13,1),%xmm1
mulss  %xmm5,%xmm1
addss  %xmm0,%xmm1
addss  (%r15,%r13,1),%xmm1
movss  %xmm1,(%r14,%r13,1)
movss  0x0(%rbp,%r13,1),%xmm0
mulss  %xmm1,%xmm1
mulss  %xmm0,%xmm0
addss  %xmm1,%xmm0
ucomiss %xmm2,%xmm0
jbe    122f <computegradient+0x11f>
jb     11e8 <computegradient+0xd8>
sqrtss %xmm0,%xmm0
jmp    120d <computegradient+0xfd>
add    $0x28,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

LLVM IR (-O0):
define void @computegradient(float* %0, i32 %1, i32 %2, float* %3, float* %4) {
%6 = alloca float*
%7 = alloca i32
%8 = alloca i32
%9 = alloca float*
%10 = alloca float*
%11 = alloca i32
%12 = alloca i32
%13 = alloca i32
%14 = alloca i32
%15 = alloca i32
%16 = alloca float
%17 = alloca float
%18 = alloca float
%19 = alloca float
%20 = alloca float
%21 = alloca float
%22 = alloca float
%23 = alloca float
%24 = alloca float
%25 = alloca float
store float* %0, float** %6
store i32 %1, i32* %7
store i32 %2, i32* %8
store float* %3, float** %9
store float* %4, float** %10
store i32 1, i32* %11
br label %26
26: ; preds = %222, %5
%27 = load i32, i32* %11
%28 = load i32, i32* %8
%29 = sub nsw i32 %28, 1
%30 = icmp slt i32 %27, %29
br i1 %30, label %31, label %225
31: ; preds = %26
store i32 1, i32* %12
br label %32
32: ; preds = %218, %31
%33 = load i32, i32* %12
%34 = load i32, i32* %7
%35 = sub nsw i32 %34, 1
%36 = icmp slt i32 %33, %35
br i1 %36, label %37, label %221
37: ; preds = %32
%38 = load i32, i32* %11
%39 = load i32, i32* %7
%40 = mul nsw i32 %38, %39
%41 = load i32, i32* %12
%42 = add nsw i32 %40, %41
store i32 %42, i32* %13
%43 = load float*, float** %6
%44 = load i32, i32* %13
%45 = sext i32 %44 to i64
%46 = getelementptr inbounds float, float* %43, i64 %45
%47 = load float, float* %46
%48 = fcmp ogt float %47, 0.000000e+00
br i1 %48, label %49, label %217
49: ; preds = %37
%50 = load float*, float** %6
%51 = load i32, i32* %13
%52 = sext i32 %51 to i64
%53 = getelementptr inbounds float, float* %50, i64 %52
%54 = load float, float* %53
%55 = fcmp olt float %54, 1.000000e+00
br i1 %55, label %56, label %217
56: ; preds = %49
%57 = load float*, float** %6
%58 = load i32, i32* %13
%59 = load i32, i32* %7
%60 = sub nsw i32 %58, %59
%61 = sub nsw i32 %60, 1
%62 = sext i32 %61 to i64
%63 = getelementptr inbounds float, float* %57, i64 %62
%64 = load float, float* %63
%65 = fneg float %64
%66 = load float*, float** %6
%67 = load i32, i32* %13
%68 = sub nsw i32 %67, 1
%69 = sext i32 %68 to i64
%70 = getelementptr inbounds float, float* %66, i64 %69
%71 = load float, float* %70
%72 = call float @llvm.fmuladd.f32(float 0xBFF6A09E80000000, float %71, float %65)
%73 = load float*, float** %6
%74 = load i32, i32* %13
%75 = load i32, i32* %7
%76 = add nsw i32 %74, %75
%77 = sub nsw i32 %76, 1
%78 = sext i32 %77 to i64
%79 = getelementptr inbounds float, float* %73, i64 %78
%80 = load float, float* %79
%81 = fsub float %72, %80
%82 = load float*, float** %6
%83 = load i32, i32* %13
%84 = load i32, i32* %7
%85 = sub nsw i32 %83, %84
%86 = add nsw i32 %85, 1
%87 = sext i32 %86 to i64
%88 = getelementptr inbounds float, float* %82, i64 %87
%89 = load float, float* %88
%90 = fadd float %81, %89
%91 = load float*, float** %6
%92 = load i32, i32* %13
%93 = add nsw i32 %92, 1
%94 = sext i32 %93 to i64
%95 = getelementptr inbounds float, float* %91, i64 %94
%96 = load float, float* %95
%97 = call float @llvm.fmuladd.f32(float 0x3FF6A09E80000000, float %96, float %90)
%98 = load float*, float** %6
%99 = load i32, i32* %13
%100 = load i32, i32* %7
%101 = add nsw i32 %99, %100
%102 = add nsw i32 %101, 1
%103 = sext i32 %102 to i64
%104 = getelementptr inbounds float, float* %98, i64 %103
%105 = load float, float* %104
%106 = fadd float %97, %105
%107 = load float*, float** %9
%108 = load i32, i32* %13
%109 = sext i32 %108 to i64
%110 = getelementptr inbounds float, float* %107, i64 %109
store float %106, float* %110
%111 = load float*, float** %6
%112 = load i32, i32* %13
%113 = load i32, i32* %7
%114 = sub nsw i32 %112, %113
%115 = sub nsw i32 %114, 1
%116 = sext i32 %115 to i64
%117 = getelementptr inbounds float, float* %111, i64 %116
%118 = load float, float* %117
%119 = fneg float %118
%120 = load float*, float** %6
%121 = load i32, i32* %13
%122 = load i32, i32* %7
%123 = sub nsw i32 %121, %122
%124 = sext i32 %123 to i64
%125 = getelementptr inbounds float, float* %120, i64 %124
%126 = load float, float* %125
%127 = call float @llvm.fmuladd.f32(float 0xBFF6A09E80000000, float %126, float %119)
%128 = load float*, float** %6
%129 = load i32, i32* %13
%130 = load i32, i32* %7
%131 = sub nsw i32 %129, %130
%132 = add nsw i32 %131, 1
%133 = sext i32 %132 to i64
%134 = getelementptr inbounds float, float* %128, i64 %133
%135 = load float, float* %134
%136 = fsub float %127, %135
%137 = load float*, float** %6
%138 = load i32, i32* %13
%139 = load i32, i32* %7
%140 = add nsw i32 %138, %139
%141 = sub nsw i32 %140, 1
%142 = sext i32 %141 to i64
%143 = getelementptr inbounds float, float* %137, i64 %142
%144 = load float, float* %143
%145 = fadd float %136, %144
%146 = load float*, float** %6
%147 = load i32, i32* %13
%148 = load i32, i32* %7
%149 = add nsw i32 %147, %148
%150 = sext i32 %149 to i64
%151 = getelementptr inbounds float, float* %146, i64 %150
%152 = load float, float* %151
%153 = call float @llvm.fmuladd.f32(float 0x3FF6A09E80000000, float %152, float %145)
%154 = load float*, float** %6
%155 = load i32, i32* %13
%156 = load i32, i32* %7
%157 = add nsw i32 %155, %156
%158 = add nsw i32 %157, 1
%159 = sext i32 %158 to i64
%160 = getelementptr inbounds float, float* %154, i64 %159
%161 = load float, float* %160
%162 = fadd float %153, %161
%163 = load float*, float** %10
%164 = load i32, i32* %13
%165 = sext i32 %164 to i64
%166 = getelementptr inbounds float, float* %163, i64 %165
store float %162, float* %166
%167 = load float*, float** %9
%168 = load i32, i32* %13
%169 = sext i32 %168 to i64
%170 = getelementptr inbounds float, float* %167, i64 %169
%171 = load float, float* %170
%172 = load float*, float** %9
%173 = load i32, i32* %13
%174 = sext i32 %173 to i64
%175 = getelementptr inbounds float, float* %172, i64 %174
%176 = load float, float* %175
%177 = load float*, float** %10
%178 = load i32, i32* %13
%179 = sext i32 %178 to i64
%180 = getelementptr inbounds float, float* %177, i64 %179
%181 = load float, float* %180
%182 = load float*, float** %10
%183 = load i32, i32* %13
%184 = sext i32 %183 to i64
%185 = getelementptr inbounds float, float* %182, i64 %184
%186 = load float, float* %185
%187 = fmul float %181, %186
%188 = call float @llvm.fmuladd.f32(float %171, float %176, float %187)
store float %188, float* %16
%189 = load float, float* %16
%190 = fcmp ogt float %189, 0.000000e+00
br i1 %190, label %191, label %216
191: ; preds = %56
%192 = load float, float* %16
%193 = call float @sqrtf(float %192)
store float %193, float* %16
%194 = load float*, float** %9
%195 = load i32, i32* %13
%196 = sext i32 %195 to i64
%197 = getelementptr inbounds float, float* %194, i64 %196
%198 = load float, float* %197
%199 = load float, float* %16
%200 = fdiv float %198, %199
%201 = load float*, float** %9
%202 = load i32, i32* %13
%203 = sext i32 %202 to i64
%204 = getelementptr inbounds float, float* %201, i64 %203
store float %200, float* %204
%205 = load float*, float** %10
%206 = load i32, i32* %13
%207 = sext i32 %206 to i64
%208 = getelementptr inbounds float, float* %205, i64 %207
%209 = load float, float* %208
%210 = load float, float* %16
%211 = fdiv float %209, %210
%212 = load float*, float** %10
%213 = load i32, i32* %13
%214 = sext i32 %213 to i64
%215 = getelementptr inbounds float, float* %212, i64 %214
store float %211, float* %215
br label %216
216: ; preds = %191, %56
br label %217
217: ; preds = %216, %49, %37
br label %218
218: ; preds = %217
%219 = load i32, i32* %12
%220 = add nsw i32 %219, 1
store i32 %220, i32* %12
br label %32
221: ; preds = %32
br label %222
222: ; preds = %221
%223 = load i32, i32* %11
%224 = add nsw i32 %223, 1
store i32 %224, i32* %11
br label %26
225: ; preds = %26
ret void
}

LLVM IR (-O2):
define void @computegradient(float* nocapture %0, i32 %1, i32 %2, float* nocapture %3, float* nocapture %4) {
%6 = icmp sgt i32 %2, 2
br i1 %6, label %7, label %85
7: ; preds = %5
%8 = add nsw i32 %2, -1
%9 = add i32 %1, -1
%10 = icmp sgt i32 %1, 2
%11 = sext i32 %1 to i64
%12 = zext i32 %8 to i64
%13 = zext i32 %9 to i64
br label %14
14: ; preds = %7, %82
%15 = phi i64 [ 1, %7 ], [ %83, %82 ]
br i1 %10, label %16, label %82
16: ; preds = %14
%17 = mul nsw i64 %15, %11
br label %18
18: ; preds = %16, %79
%19 = phi i64 [ 1, %16 ], [ %80, %79 ]
%20 = add nsw i64 %19, %17
%21 = getelementptr inbounds float, float* %0, i64 %20
%22 = load float, float* %21
%23 = fcmp ogt float %22, 0.000000e+00
%24 = fcmp olt float %22, 1.000000e+00
%25 = and i1 %23, %24
br i1 %25, label %26, label %79
26: ; preds = %18
%27 = sub nsw i64 %20, %11
%28 = add nsw i64 %27, -1
%29 = getelementptr inbounds float, float* %0, i64 %28
%30 = load float, float* %29
%31 = fneg float %30
%32 = add nsw i64 %20, -1
%33 = getelementptr inbounds float, float* %0, i64 %32
%34 = load float, float* %33
%35 = tail call float @llvm.fmuladd.f32(float %34, float 0xBFF6A09E80000000, float %31)
%36 = add nsw i64 %20, %11
%37 = add nsw i64 %36, -1
%38 = getelementptr inbounds float, float* %0, i64 %37
%39 = load float, float* %38
%40 = fsub float %35, %39
%41 = add nsw i64 %27, 1
%42 = getelementptr inbounds float, float* %0, i64 %41
%43 = load float, float* %42
%44 = fadd float %40, %43
%45 = add nsw i64 %20, 1
%46 = getelementptr inbounds float, float* %0, i64 %45
%47 = load float, float* %46
%48 = tail call float @llvm.fmuladd.f32(float %47, float 0x3FF6A09E80000000, float %44)
%49 = add nsw i64 %36, 1
%50 = getelementptr inbounds float, float* %0, i64 %49
%51 = load float, float* %50
%52 = fadd float %48, %51
%53 = getelementptr inbounds float, float* %3, i64 %20
store float %52, float* %53
%54 = load float, float* %29
%55 = fneg float %54
%56 = getelementptr inbounds float, float* %0, i64 %27
%57 = load float, float* %56
%58 = tail call float @llvm.fmuladd.f32(float %57, float 0xBFF6A09E80000000, float %55)
%59 = load float, float* %42
%60 = fsub float %58, %59
%61 = load float, float* %38
%62 = fadd float %60, %61
%63 = getelementptr inbounds float, float* %0, i64 %36
%64 = load float, float* %63
%65 = tail call float @llvm.fmuladd.f32(float %64, float 0x3FF6A09E80000000, float %62)
%66 = load float, float* %50
%67 = fadd float %65, %66
%68 = getelementptr inbounds float, float* %4, i64 %20
store float %67, float* %68
%69 = load float, float* %53
%70 = fmul float %67, %67
%71 = tail call float @llvm.fmuladd.f32(float %69, float %69, float %70)
%72 = fcmp ogt float %71, 0.000000e+00
br i1 %72, label %73, label %79
73: ; preds = %26
%74 = tail call float @sqrtf(float %71)
%75 = load float, float* %53
%76 = fdiv float %75, %74
store float %76, float* %53
%77 = load float, float* %68
%78 = fdiv float %77, %74
store float %78, float* %68
br label %79
79: ; preds = %18, %73, %26
%80 = add nuw nsw i64 %19, 1
%81 = icmp eq i64 %80, %13
br i1 %81, label %82, label %18
82: ; preds = %79, %14
%83 = add nuw nsw i64 %15, 1
%84 = icmp eq i64 %83, %12
br i1 %84, label %85, label %14
85: ; preds = %82, %5
ret void
}