#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 computegradient(double *img, int w, int h, double *gx, double *gy) {
  int i, j, k;
  double glength;
  for (i = 1; i < h - 1; i++) {
    for (j = 1; j < w - 1; j++) {
      k = i * w + j;
      if ((img[k] > 0.0) && (img[k] < 1.0)) {
        gx[k] = -img[k - w - 1] - 1.4142136 * img[k - 1] - img[k + w - 1] +
                img[k - w + 1] + 1.4142136 * img[k + 1] + img[k + w + 1];
        gy[k] = -img[k - w - 1] - 1.4142136 * img[k - w] - img[k - w + 1] +
                img[k + w - 1] + 1.4142136 * img[k + w] + img[k + w + 1];
        glength = gx[k] * gx[k] + gy[k] * gy[k];
        if (glength > 0.0) {
          glength = sqrt(glength);
          gx[k] = gx[k] / glength;
          gy[k] = gy[k] / glength;
        }
      }
    }
  }
}


x86 assembly (clang, -O0):
<computegradient>:
push   %rbp
mov    %rsp,%rbp
sub    $0x40,%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    13df <computegradient+0x2cf>
movl   $0x1,-0x28(%rbp)
mov    -0x28(%rbp),%eax
mov    -0xc(%rbp),%ecx
sub    $0x1,%ecx
cmp    %ecx,%eax
jge    13cc <computegradient+0x2bc>
mov    -0x24(%rbp),%eax
imul   -0xc(%rbp),%eax
add    -0x28(%rbp),%eax
mov    %eax,-0x2c(%rbp)
mov    -0x8(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jbe    13b9 <computegradient+0x2a9>
mov    -0x8(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movsd  0xe6f(%rip),%xmm0
ucomisd (%rax,%rcx,8),%xmm0
jbe    13b9 <computegradient+0x2a9>
mov    -0x8(%rbp),%rax
mov    -0x2c(%rbp),%ecx
sub    -0xc(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mov    -0x8(%rbp),%rax
mov    -0x2c(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm1
movsd  0xe2e(%rip),%xmm2
mulsd  %xmm2,%xmm1
addsd  %xmm0,%xmm1
mov    -0x8(%rbp),%rax
mov    -0x2c(%rbp),%ecx
add    -0xc(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
subsd  (%rax,%rcx,8),%xmm1
mov    -0x8(%rbp),%rax
mov    -0x2c(%rbp),%ecx
sub    -0xc(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
addsd  (%rax,%rcx,8),%xmm1
mov    -0x8(%rbp),%rax
mov    -0x2c(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  0xdda(%rip),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x8(%rbp),%rax
mov    -0x2c(%rbp),%ecx
add    -0xc(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
addsd  (%rax,%rcx,8),%xmm0
mov    -0x18(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x8(%rbp),%rax
mov    -0x2c(%rbp),%ecx
sub    -0xc(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mov    -0x8(%rbp),%rax
mov    -0x2c(%rbp),%ecx
sub    -0xc(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm1
movsd  0xd72(%rip),%xmm2
mulsd  %xmm2,%xmm1
addsd  %xmm0,%xmm1
mov    -0x8(%rbp),%rax
mov    -0x2c(%rbp),%ecx
sub    -0xc(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
subsd  (%rax,%rcx,8),%xmm1
mov    -0x8(%rbp),%rax
mov    -0x2c(%rbp),%ecx
add    -0xc(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
addsd  (%rax,%rcx,8),%xmm1
mov    -0x8(%rbp),%rax
mov    -0x2c(%rbp),%ecx
add    -0xc(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  0xd1e(%rip),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x8(%rbp),%rax
mov    -0x2c(%rbp),%ecx
add    -0xc(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
addsd  (%rax,%rcx,8),%xmm0
mov    -0x20(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x18(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x18(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm2
mov    -0x20(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
mov    -0x20(%rbp),%rax
movslq -0x2c(%rbp),%rcx
mulsd  (%rax,%rcx,8),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x38(%rbp)
movsd  -0x38(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jbe    13b4 <computegradient+0x2a4>
movsd  -0x38(%rbp),%xmm0
call   1030 <sqrt@plt>
movsd  %xmm0,-0x38(%rbp)
mov    -0x18(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
divsd  -0x38(%rbp),%xmm0
mov    -0x18(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x20(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
divsd  -0x38(%rbp),%xmm0
mov    -0x20(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
jmp    13b9 <computegradient+0x2a9>
jmp    13be <computegradient+0x2ae>
mov    -0x28(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x28(%rbp)
jmp    1149 <computegradient+0x39>
jmp    13d1 <computegradient+0x2c1>
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    1131 <computegradient+0x21>
add    $0x40,%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     12f0 <computegradient+0x1e0>
mov    %rdi,%rbx
mov    (%rsp),%rax
add    $0xffffffff,%eax
mov    %rax,(%rsp)
mov    0x8(%rsp),%rdx
lea    -0x1(%rdx),%eax
movslq %edx,%rdx
lea    0x10(%rdi,%rdx,8),%r12
lea    0x0(,%rdx,8),%rsi
mov    %rsi,0x10(%rsp)
lea    (%rcx,%rdx,8),%rbp
add    $0x8,%rbp
lea    (%r8,%rdx,8),%r14
add    $0x8,%r14
shl    $0x4,%rdx
lea    (%rdx,%rdi,1),%r15
add    $0x10,%r15
add    $0x10,%rbx
lea    -0x8(,%rax,8),%rcx
ff
mov    $0x1,%edx
xorpd  %xmm2,%xmm2
movsd  0xe6c(%rip),%xmm3
movsd  0xe6c(%rip),%xmm4
movsd  0xe6c(%rip),%xmm5
mov    %rcx,0x20(%rsp)
jmp    11d7 <computegradient+0xc7>
nopl   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     12f0 <computegradient+0x1e0>
mov    %rdx,0x18(%rsp)
cmpl   $0x3,0x8(%rsp)
jl     11b0 <computegradient+0xa0>
xor    %r13d,%r13d
jmp    123d <computegradient+0x12d>
call   1030 <sqrt@plt>
movsd  0xe1b(%rip),%xmm5
movsd  0xe0b(%rip),%xmm4
movsd  0xdfb(%rip),%xmm3
xorpd  %xmm2,%xmm2
mov    0x20(%rsp),%rcx
movsd  0x0(%rbp,%r13,1),%xmm1
divsd  %xmm0,%xmm1
movsd  %xmm1,0x0(%rbp,%r13,1)
movsd  (%r14,%r13,1),%xmm1
divsd  %xmm0,%xmm1
movsd  %xmm1,(%r14,%r13,1)
add    $0x8,%r13
cmp    %r13,%rcx
je     11b0 <computegradient+0xa0>
movsd  -0x8(%r12,%r13,1),%xmm0
ucomisd %xmm2,%xmm0
jbe    1230 <computegradient+0x120>
ucomisd %xmm0,%xmm3
jbe    1230 <computegradient+0x120>
movsd  -0x10(%r12,%r13,1),%xmm0
mulsd  %xmm4,%xmm0
subsd  -0x10(%rbx,%r13,1),%xmm0
subsd  -0x10(%r15,%r13,1),%xmm0
addsd  (%rbx,%r13,1),%xmm0
movsd  (%r12,%r13,1),%xmm1
mulsd  %xmm5,%xmm1
addsd  %xmm0,%xmm1
addsd  (%r15,%r13,1),%xmm1
movsd  %xmm1,0x0(%rbp,%r13,1)
movsd  -0x8(%rbx,%r13,1),%xmm0
mulsd  %xmm4,%xmm0
subsd  -0x10(%rbx,%r13,1),%xmm0
subsd  (%rbx,%r13,1),%xmm0
addsd  -0x10(%r15,%r13,1),%xmm0
movsd  -0x8(%r15,%r13,1),%xmm1
mulsd  %xmm5,%xmm1
addsd  %xmm0,%xmm1
addsd  (%r15,%r13,1),%xmm1
movsd  %xmm1,(%r14,%r13,1)
movsd  0x0(%rbp,%r13,1),%xmm0
mulsd  %xmm1,%xmm1
mulsd  %xmm0,%xmm0
addsd  %xmm1,%xmm0
ucomisd %xmm2,%xmm0
jbe    1230 <computegradient+0x120>
jb     11e8 <computegradient+0xd8>
sqrtsd %xmm0,%xmm0
jmp    120e <computegradient+0xfe>
add    $0x28,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

LLVM IR (-O2):
define void @computegradient(double* nocapture %0, i32 %1, i32 %2, double* nocapture %3, double* 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 double, double* %0, i64 %20
%22 = load double, double* %21
%23 = fcmp ogt double %22, 0.000000e+00
%24 = fcmp olt double %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 double, double* %0, i64 %28
%30 = load double, double* %29
%31 = fneg double %30
%32 = add nsw i64 %20, -1
%33 = getelementptr inbounds double, double* %0, i64 %32
%34 = load double, double* %33
%35 = tail call double @llvm.fmuladd.f64(double %34, double 0xBFF6A09E7098EF50, double %31)
%36 = add nsw i64 %20, %11
%37 = add nsw i64 %36, -1
%38 = getelementptr inbounds double, double* %0, i64 %37
%39 = load double, double* %38
%40 = fsub double %35, %39
%41 = add nsw i64 %27, 1
%42 = getelementptr inbounds double, double* %0, i64 %41
%43 = load double, double* %42
%44 = fadd double %40, %43
%45 = add nsw i64 %20, 1
%46 = getelementptr inbounds double, double* %0, i64 %45
%47 = load double, double* %46
%48 = tail call double @llvm.fmuladd.f64(double %47, double 0x3FF6A09E7098EF50, double %44)
%49 = add nsw i64 %36, 1
%50 = getelementptr inbounds double, double* %0, i64 %49
%51 = load double, double* %50
%52 = fadd double %48, %51
%53 = getelementptr inbounds double, double* %3, i64 %20
store double %52, double* %53
%54 = load double, double* %29
%55 = fneg double %54
%56 = getelementptr inbounds double, double* %0, i64 %27
%57 = load double, double* %56
%58 = tail call double @llvm.fmuladd.f64(double %57, double 0xBFF6A09E7098EF50, double %55)
%59 = load double, double* %42
%60 = fsub double %58, %59
%61 = load double, double* %38
%62 = fadd double %60, %61
%63 = getelementptr inbounds double, double* %0, i64 %36
%64 = load double, double* %63
%65 = tail call double @llvm.fmuladd.f64(double %64, double 0x3FF6A09E7098EF50, double %62)
%66 = load double, double* %50
%67 = fadd double %65, %66
%68 = getelementptr inbounds double, double* %4, i64 %20
store double %67, double* %68
%69 = load double, double* %53
%70 = fmul double %67, %67
%71 = tail call double @llvm.fmuladd.f64(double %69, double %69, double %70)
%72 = fcmp ogt double %71, 0.000000e+00
br i1 %72, label %73, label %79
73: ; preds = %26
%74 = tail call double @sqrt(double %71)
%75 = load double, double* %53
%76 = fdiv double %75, %74
store double %76, double* %53
%77 = load double, double* %68
%78 = fdiv double %77, %74
store double %78, double* %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
}