void internal_pointop(im, res, size, lut, lutsize, origin, increment,
                      warnings) register double *im,
    *res, *lut;
register double origin, increment;
register int size, lutsize, warnings;
{
  register int i, index;
  register double pos;
  register int l_unwarned = warnings;
  register int r_unwarned = warnings;
  lutsize = lutsize - 2;
  if (increment > 0)
    for (i = 0; i < size; i++) {
      pos = (im[i] - origin) / increment;
      index = (int)pos;
      if (index < 0) {
        index = 0;
        if (l_unwarned) {
          printf("Warning: Extrapolating to left of lookup table...\n");
          l_unwarned = 0;
        }
      } else if (index > lutsize) {
        index = lutsize;
        if (r_unwarned) {
          printf("Warning: Extrapolating to right of lookup table...\n");
          r_unwarned = 0;
        }
      }
      res[i] = lut[index] + (lut[index + 1] - lut[index]) * (pos - index);
      if (isnan(res[i]))
        printf("**NAN: lut[%d]=%f lut[%d]=%f pos=%f index=%d\n", index,
               lut[index], index + 1, lut[index + 1], pos, index);
    }
  else
    for (i = 0; i < size; i++) {
      res[i] = *lut;
    }
}


x86 assembly (clang, -O0):
<internal_pointop>:
push   %rbp
mov    %rsp,%rbp
sub    $0x60,%rsp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %edx,-0x14(%rbp)
mov    %rcx,-0x20(%rbp)
mov    %r8d,-0x24(%rbp)
movsd  %xmm0,-0x30(%rbp)
movsd  %xmm1,-0x38(%rbp)
mov    %r9d,-0x3c(%rbp)
mov    -0x3c(%rbp),%eax
mov    %eax,-0x54(%rbp)
mov    -0x3c(%rbp),%eax
mov    %eax,-0x58(%rbp)
mov    -0x24(%rbp),%eax
sub    $0x2,%eax
mov    %eax,-0x24(%rbp)
movsd  -0x38(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jbe    12d9 <internal_pointop+0x1b9>
movl   $0x0,-0x40(%rbp)
mov    -0x40(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    12d4 <internal_pointop+0x1b4>
mov    -0x8(%rbp),%rax
movslq -0x40(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
subsd  -0x30(%rbp),%xmm0
divsd  -0x38(%rbp),%xmm0
movsd  %xmm0,-0x50(%rbp)
cvttsd2si -0x50(%rbp),%eax
mov    %eax,-0x44(%rbp)
cmpl   $0x0,-0x44(%rbp)
jge    11dc <internal_pointop+0xbc>
movl   $0x0,-0x44(%rbp)
cmpl   $0x0,-0x54(%rbp)
je     11d7 <internal_pointop+0xb7>
lea    0xe37(%rip),%rdi
mov    $0x0,%al
call   1030 <printf@plt>
movl   $0x0,-0x54(%rbp)
jmp    1217 <internal_pointop+0xf7>
mov    -0x44(%rbp),%eax
cmp    -0x24(%rbp),%eax
jle    1212 <internal_pointop+0xf2>
mov    -0x24(%rbp),%eax
mov    %eax,-0x44(%rbp)
cmpl   $0x0,-0x58(%rbp)
je     120d <internal_pointop+0xed>
lea    0xe34(%rip),%rdi
mov    $0x0,%al
call   1030 <printf@plt>
movl   $0x0,-0x58(%rbp)
jmp    1212 <internal_pointop+0xf2>
jmp    1217 <internal_pointop+0xf7>
mov    -0x20(%rbp),%rax
movslq -0x44(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
mov    -0x20(%rbp),%rax
mov    -0x44(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x20(%rbp),%rax
movslq -0x44(%rbp),%rcx
subsd  (%rax,%rcx,8),%xmm0
movsd  -0x50(%rbp),%xmm2
cvtsi2sdl -0x44(%rbp),%xmm3
subsd  %xmm3,%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x10(%rbp),%rax
movslq -0x40(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x10(%rbp),%rax
movslq -0x40(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    $0x1,%al
call   1040 <isnan@plt>
cmp    $0x0,%eax
je     12c1 <internal_pointop+0x1a1>
mov    -0x44(%rbp),%esi
mov    -0x20(%rbp),%rax
movslq -0x44(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x44(%rbp),%edx
add    $0x1,%edx
mov    -0x20(%rbp),%rax
mov    -0x44(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm1
movsd  -0x50(%rbp),%xmm2
mov    -0x44(%rbp),%ecx
lea    0xdad(%rip),%rdi
mov    $0x3,%al
call   1030 <printf@plt>
jmp    12c6 <internal_pointop+0x1a6>
mov    -0x40(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x40(%rbp)
jmp    1177 <internal_pointop+0x57>
jmp    1314 <internal_pointop+0x1f4>
movl   $0x0,-0x40(%rbp)
mov    -0x40(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    130f <internal_pointop+0x1ef>
mov    -0x20(%rbp),%rax
movsd  (%rax),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x40(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x40(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x40(%rbp)
jmp    12e0 <internal_pointop+0x1c0>
jmp    1314 <internal_pointop+0x1f4>
add    $0x60,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<internal_pointop>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x38,%rsp
movsd  %xmm0,0x28(%rsp)
mov    %r8d,0x10(%rsp)
mov    %rcx,%rbx
mov    %rsi,%r14
mov    %rdi,0x20(%rsp)
xorpd  %xmm0,%xmm0
movsd  %xmm1,0x30(%rsp)
ucomisd %xmm0,%xmm1
jbe    127e <internal_pointop+0x14e>
mov    %r9d,0xc(%rsp)
test   %edx,%edx
jle    13cd <internal_pointop+0x29d>
addl   $0xfffffffe,0x10(%rsp)
mov    %edx,%r15d
xor    %r12d,%r12d
mov    0xc(%rsp),%eax
mov    %eax,0x14(%rsp)
jmp    119d <internal_pointop+0x6d>
nopw   0x0(%rax,%rax,1)
add    $0x1,%r12
cmp    %r12,%r15
je     13cd <internal_pointop+0x29d>
mov    0x20(%rsp),%rax
movsd  (%rax,%r12,8),%xmm0
subsd  0x28(%rsp),%xmm0
divsd  0x30(%rsp),%xmm0
cvttsd2si %xmm0,%ebp
test   %ebp,%ebp
movsd  %xmm0,0x18(%rsp)
js     11f0 <internal_pointop+0xc0>
mov    0x10(%rsp),%eax
cmp    %ebp,%eax
jge    1213 <internal_pointop+0xe3>
lea    0xe5d(%rip),%rdi
mov    %eax,%ebp
cmpl   $0x0,0xc(%rsp)
movl   $0x0,0xc(%rsp)
jne    1208 <internal_pointop+0xd8>
jmp    1213 <internal_pointop+0xe3>
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
xor    %ebp,%ebp
lea    0xe68(%rip),%rdi
cmpl   $0x0,0x14(%rsp)
movl   $0x0,0x14(%rsp)
je     1213 <internal_pointop+0xe3>
call   1030 <puts@plt>
movsd  0x18(%rsp),%xmm0
movslq %ebp,%r13
movsd  (%rbx,%r13,8),%xmm1
movsd  0x8(%rbx,%r13,8),%xmm2
subsd  %xmm1,%xmm2
xorps  %xmm3,%xmm3
cvtsi2sd %r13d,%xmm3
subsd  %xmm3,%xmm0
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,(%r14,%r12,8)
mov    $0x1,%al
call   1050 <isnan@plt>
test   %eax,%eax
je     1190 <internal_pointop+0x60>
lea    0x1(%r13),%rdx
movsd  (%rbx,%r13,8),%xmm0
movsd  (%rbx,%rdx,8),%xmm1
lea    0xd9a(%rip),%rdi
mov    %r13d,%esi
movsd  0x18(%rsp),%xmm2
mov    %r13d,%ecx
mov    $0x3,%al
call   1040 <printf@plt>
jmp    1190 <internal_pointop+0x60>
test   %edx,%edx
jle    13cd <internal_pointop+0x29d>
mov    %edx,%eax
cmp    $0x4,%edx
jb     12a3 <internal_pointop+0x173>
lea    0x8(%rbx),%rcx
cmp    %r14,%rcx
jbe    1319 <internal_pointop+0x1e9>
lea    (%r14,%rax,8),%rcx
cmp    %rbx,%rcx
jbe    1319 <internal_pointop+0x1e9>
xor    %ecx,%ecx
mov    %rcx,%rdx
not    %rdx
add    %rax,%rdx
mov    %rax,%rsi
and    $0x3,%rsi
je     12d4 <internal_pointop+0x1a4>
nopw   0x0(%rax,%rax,1)
movsd  (%rbx),%xmm0
movsd  %xmm0,(%r14,%rcx,8)
add    $0x1,%rcx
add    $0xffffffffffffffff,%rsi
jne    12c0 <internal_pointop+0x190>
cmp    $0x3,%rdx
jb     13cd <internal_pointop+0x29d>
xchg   %ax,%ax
movsd  (%rbx),%xmm0
movsd  %xmm0,(%r14,%rcx,8)
movsd  (%rbx),%xmm0
movsd  %xmm0,0x8(%r14,%rcx,8)
movsd  (%rbx),%xmm0
movsd  %xmm0,0x10(%r14,%rcx,8)
movsd  (%rbx),%xmm0
movsd  %xmm0,0x18(%r14,%rcx,8)
add    $0x4,%rcx
cmp    %rcx,%rax
jne    12e0 <internal_pointop+0x1b0>
jmp    13cd <internal_pointop+0x29d>
mov    %eax,%ecx
and    $0xfffffffc,%ecx
lea    -0x4(%rcx),%rdi
mov    %rdi,%rsi
shr    $0x2,%rsi
add    $0x1,%rsi
mov    %esi,%edx
and    $0x3,%edx
cmp    $0xc,%rdi
jae    133c <internal_pointop+0x20c>
xor    %edi,%edi
jmp    1391 <internal_pointop+0x261>
and    $0xfffffffffffffffc,%rsi
movsd  (%rbx),%xmm0
unpcklpd %xmm0,%xmm0
xor    %edi,%edi
nopw   0x0(%rax,%rax,1)
movupd %xmm0,(%r14,%rdi,8)
movupd %xmm0,0x10(%r14,%rdi,8)
movupd %xmm0,0x20(%r14,%rdi,8)
movupd %xmm0,0x30(%r14,%rdi,8)
movupd %xmm0,0x40(%r14,%rdi,8)
movupd %xmm0,0x50(%r14,%rdi,8)
movupd %xmm0,0x60(%r14,%rdi,8)
movupd %xmm0,0x70(%r14,%rdi,8)
add    $0x10,%rdi
add    $0xfffffffffffffffc,%rsi
jne    1350 <internal_pointop+0x220>
test   %rdx,%rdx
je     13c4 <internal_pointop+0x294>
movsd  (%rbx),%xmm0
unpcklpd %xmm0,%xmm0
lea    (%r14,%rdi,8),%rsi
add    $0x10,%rsi
shl    $0x5,%rdx
xor    %edi,%edi
nopl   0x0(%rax)
movupd %xmm0,-0x10(%rsi,%rdi,1)
movupd %xmm0,(%rsi,%rdi,1)
add    $0x20,%rdi
cmp    %rdi,%rdx
jne    13b0 <internal_pointop+0x280>
cmp    %rax,%rcx
jne    12a5 <internal_pointop+0x175>
add    $0x38,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

LLVM IR (-O0):
@.str = [51 x i8] c"Warning: Extrapolating to left of lookup table...\0A\00"
@.str.1 = [52 x i8] c"Warning: Extrapolating to right of lookup table...\0A\00"
@.str.2 = [46 x i8] c"**NAN: lut[%d]=%f lut[%d]=%f pos=%f index=%d\0A\00"
define void @internal_pointop(double* %0, double* %1, i32 %2, double* %3, i32 %4, double %5, double %6, i32 %7) {
%9 = alloca double*
%10 = alloca double*
%11 = alloca i32
%12 = alloca double*
%13 = alloca i32
%14 = alloca double
%15 = alloca double
%16 = alloca i32
%17 = alloca i32
%18 = alloca i32
%19 = alloca double
%20 = alloca i32
%21 = alloca i32
store double* %0, double** %9
store double* %1, double** %10
store i32 %2, i32* %11
store double* %3, double** %12
store i32 %4, i32* %13
store double %5, double* %14
store double %6, double* %15
store i32 %7, i32* %16
%22 = load i32, i32* %16
store i32 %22, i32* %20
%23 = load i32, i32* %16
store i32 %23, i32* %21
%24 = load i32, i32* %13
%25 = sub nsw i32 %24, 2
store i32 %25, i32* %13
%26 = load double, double* %15
%27 = fcmp ogt double %26, 0.000000e+00
br i1 %27, label %28, label %122
28: ; preds = %8
store i32 0, i32* %17
br label %29
29: ; preds = %118, %28
%30 = load i32, i32* %17
%31 = load i32, i32* %11
%32 = icmp slt i32 %30, %31
br i1 %32, label %33, label %121
33: ; preds = %29
%34 = load double*, double** %9
%35 = load i32, i32* %17
%36 = sext i32 %35 to i64
%37 = getelementptr inbounds double, double* %34, i64 %36
%38 = load double, double* %37
%39 = load double, double* %14
%40 = fsub double %38, %39
%41 = load double, double* %15
%42 = fdiv double %40, %41
store double %42, double* %19
%43 = load double, double* %19
%44 = fptosi double %43 to i32
store i32 %44, i32* %18
%45 = load i32, i32* %18
%46 = icmp slt i32 %45, 0
br i1 %46, label %47, label %53
47: ; preds = %33
store i32 0, i32* %18
%48 = load i32, i32* %20
%49 = icmp ne i32 %48, 0
br i1 %49, label %50, label %52
50: ; preds = %47
%51 = call i32 (i8*, ...) @printf(i8* getelementptr inbounds ([51 x i8], [51 x i8]* @.str, i64 0, i64 0))
store i32 0, i32* %20
br label %52
52: ; preds = %50, %47
br label %65
53: ; preds = %33
%54 = load i32, i32* %18
%55 = load i32, i32* %13
%56 = icmp sgt i32 %54, %55
br i1 %56, label %57, label %64
57: ; preds = %53
%58 = load i32, i32* %13
store i32 %58, i32* %18
%59 = load i32, i32* %21
%60 = icmp ne i32 %59, 0
br i1 %60, label %61, label %63
61: ; preds = %57
%62 = call i32 (i8*, ...) @printf(i8* getelementptr inbounds ([52 x i8], [52 x i8]* @.str.1, i64 0, i64 0))
store i32 0, i32* %21
br label %63
63: ; preds = %61, %57
br label %64
64: ; preds = %63, %53
br label %65
65: ; preds = %64, %52
%66 = load double*, double** %12
%67 = load i32, i32* %18
%68 = sext i32 %67 to i64
%69 = getelementptr inbounds double, double* %66, i64 %68
%70 = load double, double* %69
%71 = load double*, double** %12
%72 = load i32, i32* %18
%73 = add nsw i32 %72, 1
%74 = sext i32 %73 to i64
%75 = getelementptr inbounds double, double* %71, i64 %74
%76 = load double, double* %75
%77 = load double*, double** %12
%78 = load i32, i32* %18
%79 = sext i32 %78 to i64
%80 = getelementptr inbounds double, double* %77, i64 %79
%81 = load double, double* %80
%82 = fsub double %76, %81
%83 = load double, double* %19
%84 = load i32, i32* %18
%85 = sitofp i32 %84 to double
%86 = fsub double %83, %85
%87 = call double @llvm.fmuladd.f64(double %82, double %86, double %70)
%88 = load double*, double** %10
%89 = load i32, i32* %17
%90 = sext i32 %89 to i64
%91 = getelementptr inbounds double, double* %88, i64 %90
store double %87, double* %91
%92 = load double*, double** %10
%93 = load i32, i32* %17
%94 = sext i32 %93 to i64
%95 = getelementptr inbounds double, double* %92, i64 %94
%96 = load double, double* %95
%97 = call i32 (double, ...) bitcast (i32 (...)* @isnan to i32 (double, ...)*)(double %96)
%98 = icmp ne i32 %97, 0
br i1 %98, label %99, label %117
99: ; preds = %65
%100 = load i32, i32* %18
%101 = load double*, double** %12
%102 = load i32, i32* %18
%103 = sext i32 %102 to i64
%104 = getelementptr inbounds double, double* %101, i64 %103
%105 = load double, double* %104
%106 = load i32, i32* %18
%107 = add nsw i32 %106, 1
%108 = load double*, double** %12
%109 = load i32, i32* %18
%110 = add nsw i32 %109, 1
%111 = sext i32 %110 to i64
%112 = getelementptr inbounds double, double* %108, i64 %111
%113 = load double, double* %112
%114 = load double, double* %19
%115 = load i32, i32* %18
%116 = call i32 (i8*, ...) @printf(i8* getelementptr inbounds ([46 x i8], [46 x i8]* @.str.2, i64 0, i64 0), i32 %100, double %105, i32 %107, double %113, double %114, i32 %115)
br label %117
117: ; preds = %99, %65
br label %118
118: ; preds = %117
%119 = load i32, i32* %17
%120 = add nsw i32 %119, 1
store i32 %120, i32* %17
br label %29
121: ; preds = %29
br label %138
122: ; preds = %8
store i32 0, i32* %17
br label %123
123: ; preds = %134, %122
%124 = load i32, i32* %17
%125 = load i32, i32* %11
%126 = icmp slt i32 %124, %125
br i1 %126, label %127, label %137
127: ; preds = %123
%128 = load double*, double** %12
%129 = load double, double* %128
%130 = load double*, double** %10
%131 = load i32, i32* %17
%132 = sext i32 %131 to i64
%133 = getelementptr inbounds double, double* %130, i64 %132
store double %129, double* %133
br label %134
134: ; preds = %127
%135 = load i32, i32* %17
%136 = add nsw i32 %135, 1
store i32 %136, i32* %17
br label %123
137: ; preds = %123
br label %138
138: ; preds = %137, %121
ret void
}

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