#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 */
int abs(int);
int /*<<< orphan*/ fprintf(int /*<<< orphan*/, char *, int);
int /*<<< orphan*/ printf(char *, int, int);
int sqrt(int);
int /*<<< orphan*/ stderr;
int main() {
  int i, cons[1024], grad[1024];
  for (i = 0; i < 512; i++) {
    float A0 = i + 512, x0 = (sqrt(A0) + sqrt(A0 + 1)) * (1 << 6);
    cons[i] = A0 / x0 * (1 << 25) + x0 * (1 << 11);
    grad[i] = (1 << 24) / x0;
  }
  for (i = 512; i < 1024; i++) {
    float A0 = i, x0 = (sqrt(A0) + sqrt(A0 + 1)) * (1 << 6) * sqrt(2);
    cons[i] = A0 / x0 * (1 << 25) + x0 * (1 << 10);
    grad[i] = (1 << 24) / x0;
  }
  int A0, A1;
  int diffMax = 0;
  for (A0 = 0; A0 < (1 << 10); A0++) {
    for (A1 = 0; A1 < (1 << 14); A1++) {
      union {
        int i;
        float f;
      } u;
      u.i = (63 << 24) | (A0 << 14) | A1;
      u.f = sqrt(u.f);
      int M = (1 << 23) + (u.i & ((1 << 23) - 1));
      int fsqrtman;
      fsqrtman = cons[A0] + ((grad[A0] * A1) >> 13);
      int diff = abs(M - fsqrtman);
      if (diffMax < diff)
        diffMax = diff;
    }
  }
  fprintf(stderr, "diffMax = %d\n", diffMax);
  for (i = 0; i < 1024; i++) {
    int a, b;
    a = (cons[i] << 9) | (grad[i] >> 4);
    b = grad[i] & ((1 << 4) - 1);
    printf("%08X%X\n", a, b);
  }
  return 0;
}


x86 assembly (clang, -O0):
<main>:
push   %rbp
mov    %rsp,%rbp
sub    $0x2050,%rsp
movl   $0x0,-0x4(%rbp)
movl   $0x0,-0x8(%rbp)
cmpl   $0x200,-0x8(%rbp)
jge    122b <main+0xeb>
mov    -0x8(%rbp),%eax
add    $0x200,%eax
cvtsi2ss %eax,%xmm0
movss  %xmm0,-0x2014(%rbp)
ff
cvttss2si -0x2014(%rbp),%edi
ff
call   1060 <sqrt@plt>
mov    %eax,-0x2048(%rbp)
movss  0xe77(%rip),%xmm0
addss  -0x2014(%rbp),%xmm0
ff
cvttss2si %xmm0,%edi
call   1060 <sqrt@plt>
mov    %eax,%ecx
mov    -0x2048(%rbp),%eax
add    %ecx,%eax
shl    $0x6,%eax
cvtsi2ss %eax,%xmm0
movss  %xmm0,-0x2018(%rbp)
ff
movss  -0x2014(%rbp),%xmm0
ff
divss  -0x2018(%rbp),%xmm0
ff
movss  0xe39(%rip),%xmm1
mulss  -0x2018(%rbp),%xmm1
ff
movss  0xe1d(%rip),%xmm2
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
cvttss2si %xmm0,%ecx
movslq -0x8(%rbp),%rax
mov    %ecx,-0x1010(%rbp,%rax,4)
movss  0xdfa(%rip),%xmm0
divss  -0x2018(%rbp),%xmm0
ff
cvttss2si %xmm0,%ecx
movslq -0x8(%rbp),%rax
mov    %ecx,-0x2010(%rbp,%rax,4)
mov    -0x8(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x8(%rbp)
jmp    1159 <main+0x19>
movl   $0x200,-0x8(%rbp)
cmpl   $0x400,-0x8(%rbp)
jge    1318 <main+0x1d8>
cvtsi2ssl -0x8(%rbp),%xmm0
movss  %xmm0,-0x201c(%rbp)
ff
cvttss2si -0x201c(%rbp),%edi
ff
call   1060 <sqrt@plt>
mov    %eax,-0x2050(%rbp)
movss  0xda5(%rip),%xmm0
addss  -0x201c(%rbp),%xmm0
ff
cvttss2si %xmm0,%edi
call   1060 <sqrt@plt>
mov    %eax,%ecx
mov    -0x2050(%rbp),%eax
add    %ecx,%eax
shl    $0x6,%eax
mov    %eax,-0x204c(%rbp)
mov    $0x2,%edi
call   1060 <sqrt@plt>
mov    %eax,%ecx
mov    -0x204c(%rbp),%eax
imul   %ecx,%eax
cvtsi2ss %eax,%xmm0
movss  %xmm0,-0x2020(%rbp)
ff
movss  -0x201c(%rbp),%xmm0
ff
divss  -0x2020(%rbp),%xmm0
ff
movss  0xd44(%rip),%xmm1
mulss  -0x2020(%rbp),%xmm1
ff
movss  0xd30(%rip),%xmm2
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
cvttss2si %xmm0,%ecx
movslq -0x8(%rbp),%rax
mov    %ecx,-0x1010(%rbp,%rax,4)
movss  0xd0d(%rip),%xmm0
divss  -0x2020(%rbp),%xmm0
ff
cvttss2si %xmm0,%ecx
movslq -0x8(%rbp),%rax
mov    %ecx,-0x2010(%rbp,%rax,4)
mov    -0x8(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x8(%rbp)
jmp    1232 <main+0xf2>
movl   $0x0,-0x202c(%rbp)
movl   $0x0,-0x2024(%rbp)
cmpl   $0x400,-0x2024(%rbp)
jge    1434 <main+0x2f4>
movl   $0x0,-0x2028(%rbp)
cmpl   $0x4000,-0x2028(%rbp)
jge    141b <main+0x2db>
mov    -0x2024(%rbp),%eax
shl    $0xe,%eax
or     $0x3f000000,%eax
or     -0x2028(%rbp),%eax
mov    %eax,-0x2030(%rbp)
cvttss2si -0x2030(%rbp),%edi
ff
call   1060 <sqrt@plt>
cvtsi2ss %eax,%xmm0
movss  %xmm0,-0x2030(%rbp)
ff
mov    -0x2030(%rbp),%eax
and    $0x7fffff,%eax
add    $0x800000,%eax
mov    %eax,-0x2034(%rbp)
movslq -0x2024(%rbp),%rax
mov    -0x1010(%rbp,%rax,4),%eax
movslq -0x2024(%rbp),%rcx
mov    -0x2010(%rbp,%rcx,4),%ecx
imul   -0x2028(%rbp),%ecx
sar    $0xd,%ecx
add    %ecx,%eax
mov    %eax,-0x2038(%rbp)
mov    -0x2034(%rbp),%edi
sub    -0x2038(%rbp),%edi
call   1030 <abs@plt>
mov    %eax,-0x203c(%rbp)
mov    -0x202c(%rbp),%eax
cmp    -0x203c(%rbp),%eax
jge    1402 <main+0x2c2>
mov    -0x203c(%rbp),%eax
mov    %eax,-0x202c(%rbp)
jmp    1407 <main+0x2c7>
mov    -0x2028(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2028(%rbp)
jmp    1346 <main+0x206>
jmp    1420 <main+0x2e0>
mov    -0x2024(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2024(%rbp)
jmp    132c <main+0x1ec>
mov    0x2ba5(%rip),%rax
mov    (%rax),%edi
mov    -0x202c(%rbp),%edx
lea    0xbca(%rip),%rsi
call   1050 <fprintf@plt>
movl   $0x0,-0x8(%rbp)
cmpl   $0x400,-0x8(%rbp)
jge    14c1 <main+0x381>
movslq -0x8(%rbp),%rax
mov    -0x1010(%rbp,%rax,4),%eax
shl    $0x9,%eax
movslq -0x8(%rbp),%rcx
mov    -0x2010(%rbp,%rcx,4),%ecx
sar    $0x4,%ecx
or     %ecx,%eax
mov    %eax,-0x2040(%rbp)
movslq -0x8(%rbp),%rax
mov    -0x2010(%rbp,%rax,4),%eax
and    $0xf,%eax
mov    %eax,-0x2044(%rbp)
mov    -0x2040(%rbp),%esi
mov    -0x2044(%rbp),%edx
lea    0xb74(%rip),%rdi
call   1040 <printf@plt>
mov    -0x8(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x8(%rbp)
jmp    1456 <main+0x316>
xor    %eax,%eax
add    $0x2050,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<main>:
push   %rbp
push   %r15
push   %r14
push   %rbx
sub    $0x2018,%rsp
xor    %ebp,%ebp
nop
lea    0x200(%rbp),%eax
xorps  %xmm0,%xmm0
cvtsi2ss %eax,%xmm0
movss  %xmm0,0xc(%rsp)
cvttss2si %xmm0,%edi
call   1050 <sqrt@plt>
mov    %eax,%ebx
movss  0xc(%rsp),%xmm0
addss  0xe94(%rip),%xmm0
cvttss2si %xmm0,%edi
call   1050 <sqrt@plt>
movss  0xe8f(%rip),%xmm1
movss  0xe83(%rip),%xmm2
add    %ebx,%eax
shl    $0x6,%eax
xorps  %xmm0,%xmm0
cvtsi2ss %eax,%xmm0
movss  0xc(%rsp),%xmm3
divss  %xmm0,%xmm3
divss  %xmm0,%xmm1
mulss  0xe5d(%rip),%xmm0
mulss  %xmm2,%xmm3
addss  %xmm0,%xmm3
cvttss2si %xmm3,%eax
cvttss2si %xmm1,%ecx
mov    %eax,0x1010(%rsp,%rbp,4)
mov    %ecx,0x10(%rsp,%rbp,4)
add    $0x1,%rbp
cmp    $0x200,%rbp
jne    1140 <main+0x10>
mov    $0x200,%ebx
nopl   0x0(%rax,%rax,1)
xorps  %xmm0,%xmm0
cvtsi2ss %ebx,%xmm0
movss  %xmm0,0xc(%rsp)
cvttss2si %xmm0,%edi
call   1050 <sqrt@plt>
mov    %eax,%r14d
movss  0xc(%rsp),%xmm0
addss  0xdf9(%rip),%xmm0
cvttss2si %xmm0,%edi
call   1050 <sqrt@plt>
mov    %eax,%ebp
add    %r14d,%ebp
mov    $0x2,%edi
call   1050 <sqrt@plt>
imul   %ebp,%eax
shl    $0x6,%eax
xorps  %xmm0,%xmm0
cvtsi2ss %eax,%xmm0
movss  0xc(%rsp),%xmm2
divss  %xmm0,%xmm2
movss  0xdce(%rip),%xmm1
divss  %xmm0,%xmm1
mulss  0xdc6(%rip),%xmm0
mulss  0xdb6(%rip),%xmm2
addss  %xmm0,%xmm2
cvttss2si %xmm2,%eax
cvttss2si %xmm1,%ecx
mov    %eax,0x1010(%rsp,%rbx,4)
mov    %ecx,0x10(%rsp,%rbx,4)
add    $0x1,%rbx
cmp    $0x400,%rbx
jne    11e0 <main+0xb0>
xor    %r15d,%r15d
xor    %edx,%edx
nop
mov    %r15d,%r14d
shl    $0xe,%r14d
xor    %ebp,%ebp
nopl   0x0(%rax)
mov    %edx,%ebx
mov    %r14d,%eax
or     %ebp,%eax
or     $0x3f000000,%eax
movd   %eax,%xmm0
cvttss2si %xmm0,%edi
call   1050 <sqrt@plt>
xorps  %xmm0,%xmm0
cvtsi2ss %eax,%xmm0
movd   %xmm0,%eax
and    $0x7fffff,%eax
mov    0x10(%rsp,%r15,4),%ecx
imul   %ebp,%ecx
sar    $0xd,%ecx
or     $0x800000,%eax
add    0x1010(%rsp,%r15,4),%ecx
sub    %ecx,%eax
mov    %eax,%edx
neg    %edx
cmovs  %eax,%edx
cmp    %edx,%ebx
cmovge %ebx,%edx
add    $0x1,%ebp
cmp    $0x4000,%ebp
jne    1290 <main+0x160>
add    $0x1,%r15
cmp    $0x400,%r15
jne    1280 <main+0x150>
mov    0x2ce2(%rip),%rax
mov    (%rax),%edi
lea    0xd0d(%rip),%rsi
call   1040 <fprintf@plt>
xor    %ebp,%ebp
lea    0xd0d(%rip),%rbx
cs nopw 0x0(%rax,%rax,1)
nop
mov    0x1010(%rsp,%rbp,4),%eax
shl    $0x9,%eax
mov    0x10(%rsp,%rbp,4),%edx
mov    %edx,%esi
sar    $0x4,%esi
or     %eax,%esi
and    $0xf,%edx
mov    %rbx,%rdi
call   1030 <printf@plt>
add    $0x1,%rbp
cmp    $0x400,%rbp
jne    1320 <main+0x1f0>
xor    %eax,%eax
add    $0x2018,%rsp
pop    %rbx
pop    %r14
pop    %r15
pop    %rbp
ret

LLVM IR (-O0):
%union.anon = type { i32 }
@stderr = global i32 0
@.str = [14 x i8] c"diffMax = %d\0A\00"
@.str.1 = [8 x i8] c"%08X%X\0A\00"
define i32 @main() {
%1 = alloca i32
%2 = alloca i32
%3 = alloca [1024 x i32]
%4 = alloca [1024 x i32]
%5 = alloca float
%6 = alloca float
%7 = alloca float
%8 = alloca float
%9 = alloca i32
%10 = alloca i32
%11 = alloca i32
%12 = alloca %union.anon
%13 = alloca i32
%14 = alloca i32
%15 = alloca i32
%16 = alloca i32
%17 = alloca i32
store i32 0, i32* %1
store i32 0, i32* %2
br label %18
18: ; preds = %51, %0
%19 = load i32, i32* %2
%20 = icmp slt i32 %19, 512
br i1 %20, label %21, label %54
21: ; preds = %18
%22 = load i32, i32* %2
%23 = add nsw i32 %22, 512
%24 = sitofp i32 %23 to float
store float %24, float* %5
%25 = load float, float* %5
%26 = fptosi float %25 to i32
%27 = call i32 @sqrt(i32 %26)
%28 = load float, float* %5
%29 = fadd float %28, 1.000000e+00
%30 = fptosi float %29 to i32
%31 = call i32 @sqrt(i32 %30)
%32 = add nsw i32 %27, %31
%33 = mul nsw i32 %32, 64
%34 = sitofp i32 %33 to float
store float %34, float* %6
%35 = load float, float* %5
%36 = load float, float* %6
%37 = fdiv float %35, %36
%38 = load float, float* %6
%39 = fmul float %38, 2.048000e+03
%40 = call float @llvm.fmuladd.f32(float %37, float 0x4180000000000000, float %39)
%41 = fptosi float %40 to i32
%42 = load i32, i32* %2
%43 = sext i32 %42 to i64
%44 = getelementptr inbounds [1024 x i32], [1024 x i32]* %3, i64 0, i64 %43
store i32 %41, i32* %44
%45 = load float, float* %6
%46 = fdiv float 0x4170000000000000, %45
%47 = fptosi float %46 to i32
%48 = load i32, i32* %2
%49 = sext i32 %48 to i64
%50 = getelementptr inbounds [1024 x i32], [1024 x i32]* %4, i64 0, i64 %49
store i32 %47, i32* %50
br label %51
51: ; preds = %21
%52 = load i32, i32* %2
%53 = add nsw i32 %52, 1
store i32 %53, i32* %2
br label %18
54: ; preds = %18
store i32 512, i32* %2
br label %55
55: ; preds = %89, %54
%56 = load i32, i32* %2
%57 = icmp slt i32 %56, 1024
br i1 %57, label %58, label %92
58: ; preds = %55
%59 = load i32, i32* %2
%60 = sitofp i32 %59 to float
store float %60, float* %7
%61 = load float, float* %7
%62 = fptosi float %61 to i32
%63 = call i32 @sqrt(i32 %62)
%64 = load float, float* %7
%65 = fadd float %64, 1.000000e+00
%66 = fptosi float %65 to i32
%67 = call i32 @sqrt(i32 %66)
%68 = add nsw i32 %63, %67
%69 = mul nsw i32 %68, 64
%70 = call i32 @sqrt(i32 2)
%71 = mul nsw i32 %69, %70
%72 = sitofp i32 %71 to float
store float %72, float* %8
%73 = load float, float* %7
%74 = load float, float* %8
%75 = fdiv float %73, %74
%76 = load float, float* %8
%77 = fmul float %76, 1.024000e+03
%78 = call float @llvm.fmuladd.f32(float %75, float 0x4180000000000000, float %77)
%79 = fptosi float %78 to i32
%80 = load i32, i32* %2
%81 = sext i32 %80 to i64
%82 = getelementptr inbounds [1024 x i32], [1024 x i32]* %3, i64 0, i64 %81
store i32 %79, i32* %82
%83 = load float, float* %8
%84 = fdiv float 0x4170000000000000, %83
%85 = fptosi float %84 to i32
%86 = load i32, i32* %2
%87 = sext i32 %86 to i64
%88 = getelementptr inbounds [1024 x i32], [1024 x i32]* %4, i64 0, i64 %87
store i32 %85, i32* %88
br label %89
89: ; preds = %58
%90 = load i32, i32* %2
%91 = add nsw i32 %90, 1
store i32 %91, i32* %2
br label %55
92: ; preds = %55
store i32 0, i32* %11
store i32 0, i32* %9
br label %93
93: ; preds = %143, %92
%94 = load i32, i32* %9
%95 = icmp slt i32 %94, 1024
br i1 %95, label %96, label %146
96: ; preds = %93
store i32 0, i32* %10
br label %97
97: ; preds = %139, %96
%98 = load i32, i32* %10
%99 = icmp slt i32 %98, 16384
br i1 %99, label %100, label %142
100: ; preds = %97
%101 = load i32, i32* %9
%102 = shl i32 %101, 14
%103 = or i32 1056964608, %102
%104 = load i32, i32* %10
%105 = or i32 %103, %104
%106 = bitcast %union.anon* %12 to i32*
store i32 %105, i32* %106
%107 = bitcast %union.anon* %12 to float*
%108 = load float, float* %107
%109 = fptosi float %108 to i32
%110 = call i32 @sqrt(i32 %109)
%111 = sitofp i32 %110 to float
%112 = bitcast %union.anon* %12 to float*
store float %111, float* %112
%113 = bitcast %union.anon* %12 to i32*
%114 = load i32, i32* %113
%115 = and i32 %114, 8388607
%116 = add nsw i32 8388608, %115
store i32 %116, i32* %13
%117 = load i32, i32* %9
%118 = sext i32 %117 to i64
%119 = getelementptr inbounds [1024 x i32], [1024 x i32]* %3, i64 0, i64 %118
%120 = load i32, i32* %119
%121 = load i32, i32* %9
%122 = sext i32 %121 to i64
%123 = getelementptr inbounds [1024 x i32], [1024 x i32]* %4, i64 0, i64 %122
%124 = load i32, i32* %123
%125 = load i32, i32* %10
%126 = mul nsw i32 %124, %125
%127 = ashr i32 %126, 13
%128 = add nsw i32 %120, %127
store i32 %128, i32* %14
%129 = load i32, i32* %13
%130 = load i32, i32* %14
%131 = sub nsw i32 %129, %130
%132 = call i32 @abs(i32 %131)
store i32 %132, i32* %15
%133 = load i32, i32* %11
%134 = load i32, i32* %15
%135 = icmp slt i32 %133, %134
br i1 %135, label %136, label %138
136: ; preds = %100
%137 = load i32, i32* %15
store i32 %137, i32* %11
br label %138
138: ; preds = %136, %100
br label %139
139: ; preds = %138
%140 = load i32, i32* %10
%141 = add nsw i32 %140, 1
store i32 %141, i32* %10
br label %97
142: ; preds = %97
br label %143
143: ; preds = %142
%144 = load i32, i32* %9
%145 = add nsw i32 %144, 1
store i32 %145, i32* %9
br label %93
146: ; preds = %93
%147 = load i32, i32* @stderr
%148 = load i32, i32* %11
%149 = call i32 @fprintf(i32 %147, i8* getelementptr inbounds ([14 x i8], [14 x i8]* @.str, i64 0, i64 0), i32 %148)
store i32 0, i32* %2
br label %150
150: ; preds = %173, %146
%151 = load i32, i32* %2
%152 = icmp slt i32 %151, 1024
br i1 %152, label %153, label %176
153: ; preds = %150
%154 = load i32, i32* %2
%155 = sext i32 %154 to i64
%156 = getelementptr inbounds [1024 x i32], [1024 x i32]* %3, i64 0, i64 %155
%157 = load i32, i32* %156
%158 = shl i32 %157, 9
%159 = load i32, i32* %2
%160 = sext i32 %159 to i64
%161 = getelementptr inbounds [1024 x i32], [1024 x i32]* %4, i64 0, i64 %160
%162 = load i32, i32* %161
%163 = ashr i32 %162, 4
%164 = or i32 %158, %163
store i32 %164, i32* %16
%165 = load i32, i32* %2
%166 = sext i32 %165 to i64
%167 = getelementptr inbounds [1024 x i32], [1024 x i32]* %4, i64 0, i64 %166
%168 = load i32, i32* %167
%169 = and i32 %168, 15
store i32 %169, i32* %17
%170 = load i32, i32* %16
%171 = load i32, i32* %17
%172 = call i32 @printf(i8* getelementptr inbounds ([8 x i8], [8 x i8]* @.str.1, i64 0, i64 0), i32 %170, i32 %171)
br label %173
173: ; preds = %153
%174 = load i32, i32* %2
%175 = add nsw i32 %174, 1
store i32 %175, i32* %2
br label %150
176: ; preds = %150
ret i32 0
}

LLVM IR (-O2):
@stderr = global i32 0
@.str = [14 x i8] c"diffMax = %d\0A\00"
@.str.1 = [8 x i8] c"%08X%X\0A\00"
define i32 @main() {
%1 = alloca [1024 x i32]
%2 = alloca [1024 x i32]
%3 = bitcast [1024 x i32]* %1 to i8*
%4 = bitcast [1024 x i32]* %2 to i8*
br label %5
5: ; preds = %0, %5
%6 = phi i64 [ 0, %0 ], [ %26, %5 ]
%7 = trunc i64 %6 to i32
%8 = add i32 %7, 512
%9 = sitofp i32 %8 to float
%10 = fptosi float %9 to i32
%11 = tail call i32 @sqrt(i32 %10)
%12 = fadd float %9, 1.000000e+00
%13 = fptosi float %12 to i32
%14 = tail call i32 @sqrt(i32 %13)
%15 = add nsw i32 %14, %11
%16 = shl nsw i32 %15, 6
%17 = sitofp i32 %16 to float
%18 = fdiv float %9, %17
%19 = fmul float %17, 2.048000e+03
%20 = tail call float @llvm.fmuladd.f32(float %18, float 0x4180000000000000, float %19)
%21 = fptosi float %20 to i32
%22 = getelementptr inbounds [1024 x i32], [1024 x i32]* %1, i64 0, i64 %6
store i32 %21, i32* %22
%23 = fdiv float 0x4170000000000000, %17
%24 = fptosi float %23 to i32
%25 = getelementptr inbounds [1024 x i32], [1024 x i32]* %2, i64 0, i64 %6
store i32 %24, i32* %25
%26 = add nuw nsw i64 %6, 1
%27 = icmp eq i64 %26, 512
br i1 %27, label %28, label %5
28: ; preds = %5, %28
%29 = phi i64 [ %50, %28 ], [ 512, %5 ]
%30 = trunc i64 %29 to i32
%31 = sitofp i32 %30 to float
%32 = fptosi float %31 to i32
%33 = tail call i32 @sqrt(i32 %32)
%34 = fadd float %31, 1.000000e+00
%35 = fptosi float %34 to i32
%36 = tail call i32 @sqrt(i32 %35)
%37 = add nsw i32 %36, %33
%38 = shl nsw i32 %37, 6
%39 = tail call i32 @sqrt(i32 2)
%40 = mul nsw i32 %38, %39
%41 = sitofp i32 %40 to float
%42 = fdiv float %31, %41
%43 = fmul float %41, 1.024000e+03
%44 = tail call float @llvm.fmuladd.f32(float %42, float 0x4180000000000000, float %43)
%45 = fptosi float %44 to i32
%46 = getelementptr inbounds [1024 x i32], [1024 x i32]* %1, i64 0, i64 %29
store i32 %45, i32* %46
%47 = fdiv float 0x4170000000000000, %41
%48 = fptosi float %47 to i32
%49 = getelementptr inbounds [1024 x i32], [1024 x i32]* %2, i64 0, i64 %29
store i32 %48, i32* %49
%50 = add nuw nsw i64 %29, 1
%51 = icmp eq i64 %50, 1024
br i1 %51, label %52, label %28
52: ; preds = %28, %82
%53 = phi i64 [ %83, %82 ], [ 0, %28 ]
%54 = phi i32 [ %79, %82 ], [ 0, %28 ]
%55 = getelementptr inbounds [1024 x i32], [1024 x i32]* %1, i64 0, i64 %53
%56 = getelementptr inbounds [1024 x i32], [1024 x i32]* %2, i64 0, i64 %53
%57 = trunc i64 %53 to i32
%58 = shl i32 %57, 14
br label %59
59: ; preds = %52, %59
%60 = phi i32 [ 0, %52 ], [ %80, %59 ]
%61 = phi i32 [ %54, %52 ], [ %79, %59 ]
%62 = or i32 %58, %60
%63 = or i32 %62, 1056964608
%64 = bitcast i32 %63 to float
%65 = fptosi float %64 to i32
%66 = tail call i32 @sqrt(i32 %65)
%67 = sitofp i32 %66 to float
%68 = bitcast float %67 to i32
%69 = and i32 %68, 8388607
%70 = load i32, i32* %55
%71 = load i32, i32* %56
%72 = mul nsw i32 %71, %60
%73 = ashr i32 %72, 13
%74 = or i32 %69, 8388608
%75 = add i32 %70, %73
%76 = sub i32 %74, %75
%77 = tail call i32 @llvm.abs.i32(i32 %76, i1 true)
%78 = icmp slt i32 %61, %77
%79 = select i1 %78, i32 %77, i32 %61
%80 = add nuw nsw i32 %60, 1
%81 = icmp eq i32 %80, 16384
br i1 %81, label %82, label %59
82: ; preds = %59
%83 = add nuw nsw i64 %53, 1
%84 = icmp eq i64 %83, 1024
br i1 %84, label %85, label %52
85: ; preds = %82
%86 = load i32, i32* @stderr
%87 = tail call i32 @fprintf(i32 %86, i8* getelementptr inbounds ([14 x i8], [14 x i8]* @.str, i64 0, i64 0), i32 %79)
br label %88
88: ; preds = %85, %88
%89 = phi i64 [ 0, %85 ], [ %99, %88 ]
%90 = getelementptr inbounds [1024 x i32], [1024 x i32]* %1, i64 0, i64 %89
%91 = load i32, i32* %90
%92 = shl i32 %91, 9
%93 = getelementptr inbounds [1024 x i32], [1024 x i32]* %2, i64 0, i64 %89
%94 = load i32, i32* %93
%95 = ashr i32 %94, 4
%96 = or i32 %95, %92
%97 = and i32 %94, 15
%98 = tail call i32 @printf(i8* getelementptr inbounds ([8 x i8], [8 x i8]* @.str.1, i64 0, i64 0), i32 %96, i32 %97)
%99 = add nuw nsw i64 %89, 1
%100 = icmp eq i64 %99, 1024
br i1 %100, label %101, label %88
101: ; preds = %88
ret i32 0
}