void warp_image(image, warped, u, v, nx, ny) float *image, *warped, *u, *v;
int nx, ny;
{
  float linear_interp(), new_i, new_j, remainder_i, remainder_j;
  int index, i, j, floor_i, floor_j, ceil_i, ceil_j;
  for (i = 0; i < ny; i++) {
    for (j = 0; j < nx; j++) {
      index = (i * nx) + j;
      new_i = ((float)i) - v[index];
      new_j = ((float)j) - u[index];
      if (new_j < 0 || new_j > (nx - 1) || new_i < 0 || new_i > (ny - 1))
        warped[index] = 0.0;
      else {
        floor_i = ((int)floor(new_i)) % ny;
        ceil_i = ((int)ceil(new_i)) % ny;
        remainder_i = new_i - ((float)floor_i);
        floor_j = ((int)floor(new_j)) % nx;
        ceil_j = ((int)ceil(new_j)) % nx;
        remainder_j = new_j - ((float)floor_j);
        warped[index] = linear_interp(
            linear_interp(image[(floor_i * nx) + floor_j],
                          image[(ceil_i * nx) + floor_j], remainder_i),
            linear_interp(image[(floor_i * nx) + ceil_j],
                          image[(ceil_i * nx) + ceil_j], remainder_i),
            remainder_j);
      }
    }
  }
}


x86 assembly (clang, -O0):
<warp_image>:
push   %rbp
mov    %rsp,%rbp
sub    $0x60,%rsp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %rcx,-0x20(%rbp)
mov    %r8d,-0x24(%rbp)
mov    %r9d,-0x28(%rbp)
movl   $0x0,-0x40(%rbp)
mov    -0x40(%rbp),%eax
cmp    -0x28(%rbp),%eax
jge    137a <warp_image+0x24a>
movl   $0x0,-0x44(%rbp)
mov    -0x44(%rbp),%eax
cmp    -0x24(%rbp),%eax
jge    1367 <warp_image+0x237>
mov    -0x40(%rbp),%eax
imul   -0x24(%rbp),%eax
add    -0x44(%rbp),%eax
mov    %eax,-0x3c(%rbp)
cvtsi2ssl -0x40(%rbp),%xmm0
mov    -0x20(%rbp),%rax
movslq -0x3c(%rbp),%rcx
subss  (%rax,%rcx,4),%xmm0
movss  %xmm0,-0x2c(%rbp)
cvtsi2ssl -0x44(%rbp),%xmm0
mov    -0x18(%rbp),%rax
movslq -0x3c(%rbp),%rcx
subss  (%rax,%rcx,4),%xmm0
movss  %xmm0,-0x30(%rbp)
xorps  %xmm0,%xmm0
ucomiss -0x30(%rbp),%xmm0
ja     11fb <warp_image+0xcb>
movss  -0x30(%rbp),%xmm0
mov    -0x24(%rbp),%eax
sub    $0x1,%eax
cvtsi2ss %eax,%xmm1
ucomiss %xmm1,%xmm0
ja     11fb <warp_image+0xcb>
xorps  %xmm0,%xmm0
ucomiss -0x2c(%rbp),%xmm0
ja     11fb <warp_image+0xcb>
movss  -0x2c(%rbp),%xmm0
mov    -0x28(%rbp),%eax
sub    $0x1,%eax
cvtsi2ss %eax,%xmm1
ucomiss %xmm1,%xmm0
jbe    1210 <warp_image+0xe0>
mov    -0x10(%rbp),%rax
movslq -0x3c(%rbp),%rcx
xorps  %xmm0,%xmm0
movss  %xmm0,(%rax,%rcx,4)
jmp    1354 <warp_image+0x224>
movss  -0x2c(%rbp),%xmm0
cvtss2sd %xmm0,%xmm0
call   1050 <floor@plt>
cvttsd2si %xmm0,%eax
cltd
idivl  -0x28(%rbp)
mov    %edx,-0x48(%rbp)
movss  -0x2c(%rbp),%xmm0
cvtss2sd %xmm0,%xmm0
call   1030 <ceil@plt>
cvttsd2si %xmm0,%eax
cltd
idivl  -0x28(%rbp)
mov    %edx,-0x50(%rbp)
movss  -0x2c(%rbp),%xmm0
cvtsi2ssl -0x48(%rbp),%xmm1
subss  %xmm1,%xmm0
movss  %xmm0,-0x34(%rbp)
movss  -0x30(%rbp),%xmm0
cvtss2sd %xmm0,%xmm0
call   1050 <floor@plt>
cvttsd2si %xmm0,%eax
cltd
idivl  -0x24(%rbp)
mov    %edx,-0x4c(%rbp)
movss  -0x30(%rbp),%xmm0
cvtss2sd %xmm0,%xmm0
call   1030 <ceil@plt>
cvttsd2si %xmm0,%eax
cltd
idivl  -0x24(%rbp)
mov    %edx,-0x54(%rbp)
movss  -0x30(%rbp),%xmm0
cvtsi2ssl -0x4c(%rbp),%xmm1
subss  %xmm1,%xmm0
movss  %xmm0,-0x38(%rbp)
mov    -0x8(%rbp),%rax
mov    -0x48(%rbp),%ecx
imul   -0x24(%rbp),%ecx
add    -0x4c(%rbp),%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm0
cvtss2sd %xmm0,%xmm0
mov    -0x8(%rbp),%rax
mov    -0x50(%rbp),%ecx
imul   -0x24(%rbp),%ecx
add    -0x4c(%rbp),%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm1
cvtss2sd %xmm1,%xmm1
movss  -0x34(%rbp),%xmm2
cvtss2sd %xmm2,%xmm2
mov    $0x3,%al
call   1040 <linear_interp@plt>
cvtss2sd %xmm0,%xmm0
movsd  %xmm0,-0x60(%rbp)
mov    -0x8(%rbp),%rax
mov    -0x48(%rbp),%ecx
imul   -0x24(%rbp),%ecx
add    -0x54(%rbp),%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm0
cvtss2sd %xmm0,%xmm0
mov    -0x8(%rbp),%rax
mov    -0x50(%rbp),%ecx
imul   -0x24(%rbp),%ecx
add    -0x54(%rbp),%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm1
cvtss2sd %xmm1,%xmm1
movss  -0x34(%rbp),%xmm2
cvtss2sd %xmm2,%xmm2
mov    $0x3,%al
call   1040 <linear_interp@plt>
movaps %xmm0,%xmm1
movsd  -0x60(%rbp),%xmm0
cvtss2sd %xmm1,%xmm1
movss  -0x38(%rbp),%xmm2
cvtss2sd %xmm2,%xmm2
mov    $0x3,%al
call   1040 <linear_interp@plt>
mov    -0x10(%rbp),%rax
movslq -0x3c(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
jmp    1359 <warp_image+0x229>
mov    -0x44(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x44(%rbp)
jmp    116a <warp_image+0x3a>
jmp    136c <warp_image+0x23c>
mov    -0x40(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x40(%rbp)
jmp    1157 <warp_image+0x27>
add    $0x60,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<warp_image>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x78,%rsp
mov    %r8,0x20(%rsp)
mov    %rdi,0x58(%rsp)
mov    %r9,0x38(%rsp)
test   %r9d,%r9d
jle    13a5 <warp_image+0x275>
mov    0x20(%rsp),%rdi
mov    %edi,%ebx
lea    -0x1(%rdi),%eax
cvtsi2ss %eax,%xmm1
mov    0x38(%rsp),%rbp
lea    -0x1(%rbp),%eax
cvtsi2ss %eax,%xmm0
movss  %xmm0,0x30(%rsp)
movslq %edi,%rax
mov    %ebp,%edi
mov    %rdi,0x40(%rsp)
shl    $0x2,%rax
mov    %rax,0x48(%rsp)
xor    %eax,%eax
mov    %rax,0x18(%rsp)
mov    %rbx,0x50(%rsp)
movss  %xmm1,0x2c(%rsp)
jmp    11ca <warp_image+0x9a>
nopl   (%rax)
mov    0x18(%rsp),%rdi
add    $0x1,%rdi
mov    0x48(%rsp),%rax
add    %rax,%rsi
add    %rax,%rcx
add    %rax,%rdx
mov    %rdi,%rax
mov    %rdi,0x18(%rsp)
cmp    0x40(%rsp),%rdi
je     13a5 <warp_image+0x275>
cmpl   $0x0,0x20(%rsp)
jle    11a0 <warp_image+0x70>
xorps  %xmm0,%xmm0
cvtsi2ssl 0x18(%rsp),%xmm0
movss  %xmm0,0x34(%rsp)
xor    %r13d,%r13d
mov    %rcx,0x70(%rsp)
mov    %rdx,0x68(%rsp)
mov    %rsi,0x60(%rsp)
jmp    120f <warp_image+0xdf>
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
movss  %xmm0,(%rsi,%r13,4)
add    $0x1,%r13
cmp    %r13,%rbx
je     11a0 <warp_image+0x70>
xorps  %xmm3,%xmm3
cvtsi2ss %r13d,%xmm3
subss  (%rdx,%r13,4),%xmm3
xorps  %xmm0,%xmm0
ucomiss %xmm3,%xmm0
ja     1200 <warp_image+0xd0>
ucomiss %xmm1,%xmm3
ja     1200 <warp_image+0xd0>
movss  0x34(%rsp),%xmm4
subss  (%rcx,%r13,4),%xmm4
xorps  %xmm2,%xmm2
ucomiss %xmm4,%xmm2
ja     1200 <warp_image+0xd0>
ucomiss 0x30(%rsp),%xmm4
ja     1200 <warp_image+0xd0>
xorps  %xmm0,%xmm0
cvtss2sd %xmm4,%xmm0
movsd  %xmm0,0x8(%rsp)
movss  %xmm3,0x4(%rsp)
movss  %xmm4,0x10(%rsp)
call   1050 <floor@plt>
cvttsd2si %xmm0,%eax
cltd
mov    0x38(%rsp),%rbp
idiv   %ebp
mov    %edx,%r14d
movsd  0x8(%rsp),%xmm0
call   1030 <ceil@plt>
cvttsd2si %xmm0,%eax
xorps  %xmm0,%xmm0
cvtsi2ss %r14d,%xmm0
cltd
idiv   %ebp
mov    %edx,%r15d
movss  0x10(%rsp),%xmm1
subss  %xmm0,%xmm1
movss  %xmm1,0x10(%rsp)
movss  0x4(%rsp),%xmm0
cvtss2sd %xmm0,%xmm0
movsd  %xmm0,0x8(%rsp)
call   1050 <floor@plt>
cvttsd2si %xmm0,%eax
cltd
idiv   %ebx
mov    %edx,%ebp
movsd  0x8(%rsp),%xmm0
call   1030 <ceil@plt>
cvttsd2si %xmm0,%eax
cltd
mov    0x20(%rsp),%rcx
idiv   %ecx
mov    %edx,%r12d
xorps  %xmm0,%xmm0
cvtsi2ss %ebp,%xmm0
movss  0x4(%rsp),%xmm1
subss  %xmm0,%xmm1
movss  %xmm1,0x4(%rsp)
imul   %ecx,%r14d
lea    (%r14,%rbp,1),%eax
cltq
mov    0x58(%rsp),%rbx
movss  (%rbx,%rax,4),%xmm0
cvtss2sd %xmm0,%xmm0
imul   %ecx,%r15d
add    %r15d,%ebp
movslq %ebp,%rax
movss  (%rbx,%rax,4),%xmm1
cvtss2sd %xmm1,%xmm1
movss  0x10(%rsp),%xmm2
cvtss2sd %xmm2,%xmm2
movsd  %xmm2,0x8(%rsp)
mov    $0x3,%al
call   1040 <linear_interp@plt>
cvtss2sd %xmm0,%xmm0
movsd  %xmm0,0x10(%rsp)
add    %r12d,%r14d
movslq %r14d,%rax
movss  (%rbx,%rax,4),%xmm0
cvtss2sd %xmm0,%xmm0
add    %r12d,%r15d
movslq %r15d,%rax
movss  (%rbx,%rax,4),%xmm1
mov    0x50(%rsp),%rbx
cvtss2sd %xmm1,%xmm1
movsd  0x8(%rsp),%xmm2
mov    $0x3,%al
call   1040 <linear_interp@plt>
xorps  %xmm1,%xmm1
cvtss2sd %xmm0,%xmm1
movss  0x4(%rsp),%xmm0
xorps  %xmm2,%xmm2
cvtss2sd %xmm0,%xmm2
movsd  0x10(%rsp),%xmm0
mov    $0x3,%al
call   1040 <linear_interp@plt>
movss  0x2c(%rsp),%xmm1
mov    0x60(%rsp),%rsi
mov    0x68(%rsp),%rdx
mov    0x70(%rsp),%rcx
jmp    1200 <warp_image+0xd0>
add    $0x78,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

LLVM IR (-O0):
define void @warp_image(float* %0, float* %1, float* %2, float* %3, i32 %4, i32 %5) {
%7 = alloca float*
%8 = alloca float*
%9 = alloca float*
%10 = alloca float*
%11 = alloca i32
%12 = alloca i32
%13 = alloca float
%14 = alloca float
%15 = alloca float
%16 = alloca float
%17 = alloca i32
%18 = alloca i32
%19 = alloca i32
%20 = alloca i32
%21 = alloca i32
%22 = alloca i32
%23 = alloca i32
store float* %0, float** %7
store float* %1, float** %8
store float* %2, float** %9
store float* %3, float** %10
store i32 %4, i32* %11
store i32 %5, i32* %12
store i32 0, i32* %18
br label %24
24: ; preds = %170, %6
%25 = load i32, i32* %18
%26 = load i32, i32* %12
%27 = icmp slt i32 %25, %26
br i1 %27, label %28, label %173
28: ; preds = %24
store i32 0, i32* %19
br label %29
29: ; preds = %166, %28
%30 = load i32, i32* %19
%31 = load i32, i32* %11
%32 = icmp slt i32 %30, %31
br i1 %32, label %33, label %169
33: ; preds = %29
%34 = load i32, i32* %18
%35 = load i32, i32* %11
%36 = mul nsw i32 %34, %35
%37 = load i32, i32* %19
%38 = add nsw i32 %36, %37
store i32 %38, i32* %17
%39 = load i32, i32* %18
%40 = sitofp i32 %39 to float
%41 = load float*, float** %10
%42 = load i32, i32* %17
%43 = sext i32 %42 to i64
%44 = getelementptr inbounds float, float* %41, i64 %43
%45 = load float, float* %44
%46 = fsub float %40, %45
store float %46, float* %13
%47 = load i32, i32* %19
%48 = sitofp i32 %47 to float
%49 = load float*, float** %9
%50 = load i32, i32* %17
%51 = sext i32 %50 to i64
%52 = getelementptr inbounds float, float* %49, i64 %51
%53 = load float, float* %52
%54 = fsub float %48, %53
store float %54, float* %14
%55 = load float, float* %14
%56 = fcmp olt float %55, 0.000000e+00
br i1 %56, label %72, label %57
57: ; preds = %33
%58 = load float, float* %14
%59 = load i32, i32* %11
%60 = sub nsw i32 %59, 1
%61 = sitofp i32 %60 to float
%62 = fcmp ogt float %58, %61
br i1 %62, label %72, label %63
63: ; preds = %57
%64 = load float, float* %13
%65 = fcmp olt float %64, 0.000000e+00
br i1 %65, label %72, label %66
66: ; preds = %63
%67 = load float, float* %13
%68 = load i32, i32* %12
%69 = sub nsw i32 %68, 1
%70 = sitofp i32 %69 to float
%71 = fcmp ogt float %67, %70
br i1 %71, label %72, label %77
72: ; preds = %66, %63, %57, %33
%73 = load float*, float** %8
%74 = load i32, i32* %17
%75 = sext i32 %74 to i64
%76 = getelementptr inbounds float, float* %73, i64 %75
store float 0.000000e+00, float* %76
br label %165
77: ; preds = %66
%78 = load float, float* %13
%79 = fpext float %78 to double
%80 = call double @llvm.floor.f64(double %79)
%81 = fptosi double %80 to i32
%82 = load i32, i32* %12
%83 = srem i32 %81, %82
store i32 %83, i32* %20
%84 = load float, float* %13
%85 = fpext float %84 to double
%86 = call double @llvm.ceil.f64(double %85)
%87 = fptosi double %86 to i32
%88 = load i32, i32* %12
%89 = srem i32 %87, %88
store i32 %89, i32* %22
%90 = load float, float* %13
%91 = load i32, i32* %20
%92 = sitofp i32 %91 to float
%93 = fsub float %90, %92
store float %93, float* %15
%94 = load float, float* %14
%95 = fpext float %94 to double
%96 = call double @llvm.floor.f64(double %95)
%97 = fptosi double %96 to i32
%98 = load i32, i32* %11
%99 = srem i32 %97, %98
store i32 %99, i32* %21
%100 = load float, float* %14
%101 = fpext float %100 to double
%102 = call double @llvm.ceil.f64(double %101)
%103 = fptosi double %102 to i32
%104 = load i32, i32* %11
%105 = srem i32 %103, %104
store i32 %105, i32* %23
%106 = load float, float* %14
%107 = load i32, i32* %21
%108 = sitofp i32 %107 to float
%109 = fsub float %106, %108
store float %109, float* %16
%110 = load float*, float** %7
%111 = load i32, i32* %20
%112 = load i32, i32* %11
%113 = mul nsw i32 %111, %112
%114 = load i32, i32* %21
%115 = add nsw i32 %113, %114
%116 = sext i32 %115 to i64
%117 = getelementptr inbounds float, float* %110, i64 %116
%118 = load float, float* %117
%119 = fpext float %118 to double
%120 = load float*, float** %7
%121 = load i32, i32* %22
%122 = load i32, i32* %11
%123 = mul nsw i32 %121, %122
%124 = load i32, i32* %21
%125 = add nsw i32 %123, %124
%126 = sext i32 %125 to i64
%127 = getelementptr inbounds float, float* %120, i64 %126
%128 = load float, float* %127
%129 = fpext float %128 to double
%130 = load float, float* %15
%131 = fpext float %130 to double
%132 = call float (double, double, double, ...) bitcast (float (...)* @linear_interp to float (double, double, double, ...)*)(double %119, double %129, double %131)
%133 = fpext float %132 to double
%134 = load float*, float** %7
%135 = load i32, i32* %20
%136 = load i32, i32* %11
%137 = mul nsw i32 %135, %136
%138 = load i32, i32* %23
%139 = add nsw i32 %137, %138
%140 = sext i32 %139 to i64
%141 = getelementptr inbounds float, float* %134, i64 %140
%142 = load float, float* %141
%143 = fpext float %142 to double
%144 = load float*, float** %7
%145 = load i32, i32* %22
%146 = load i32, i32* %11
%147 = mul nsw i32 %145, %146
%148 = load i32, i32* %23
%149 = add nsw i32 %147, %148
%150 = sext i32 %149 to i64
%151 = getelementptr inbounds float, float* %144, i64 %150
%152 = load float, float* %151
%153 = fpext float %152 to double
%154 = load float, float* %15
%155 = fpext float %154 to double
%156 = call float (double, double, double, ...) bitcast (float (...)* @linear_interp to float (double, double, double, ...)*)(double %143, double %153, double %155)
%157 = fpext float %156 to double
%158 = load float, float* %16
%159 = fpext float %158 to double
%160 = call float (double, double, double, ...) bitcast (float (...)* @linear_interp to float (double, double, double, ...)*)(double %133, double %157, double %159)
%161 = load float*, float** %8
%162 = load i32, i32* %17
%163 = sext i32 %162 to i64
%164 = getelementptr inbounds float, float* %161, i64 %163
store float %160, float* %164
br label %165
165: ; preds = %77, %72
br label %166
166: ; preds = %165
%167 = load i32, i32* %19
%168 = add nsw i32 %167, 1
store i32 %168, i32* %19
br label %29
169: ; preds = %29
br label %170
170: ; preds = %169
%171 = load i32, i32* %18
%172 = add nsw i32 %171, 1
store i32 %172, i32* %18
br label %24
173: ; preds = %24
ret void
}

LLVM IR (-O2):
define void @warp_image(float* nocapture %0, float* nocapture writeonly %1, float* nocapture %2, float* nocapture %3, i32 %4, i32 %5) {
%7 = icmp sgt i32 %5, 0
br i1 %7, label %8, label %97
8: ; preds = %6
%9 = icmp sgt i32 %4, 0
%10 = add nsw i32 %4, -1
%11 = sitofp i32 %10 to float
%12 = add nsw i32 %5, -1
%13 = sitofp i32 %12 to float
%14 = sext i32 %4 to i64
%15 = zext i32 %5 to i64
%16 = zext i32 %4 to i64
br label %17
17: ; preds = %8, %94
%18 = phi i64 [ 0, %8 ], [ %95, %94 ]
br i1 %9, label %19, label %94
19: ; preds = %17
%20 = mul nsw i64 %18, %14
%21 = trunc i64 %18 to i32
%22 = sitofp i32 %21 to float
br label %23
23: ; preds = %19, %89
%24 = phi i64 [ 0, %19 ], [ %92, %89 ]
%25 = add nsw i64 %24, %20
%26 = getelementptr inbounds float, float* %3, i64 %25
%27 = load float, float* %26
%28 = fsub float %22, %27
%29 = trunc i64 %24 to i32
%30 = sitofp i32 %29 to float
%31 = getelementptr inbounds float, float* %2, i64 %25
%32 = load float, float* %31
%33 = fsub float %30, %32
%34 = fcmp olt float %33, 0.000000e+00
br i1 %34, label %89, label %35
35: ; preds = %23
%36 = fcmp ogt float %33, %11
%37 = fcmp olt float %28, 0.000000e+00
%38 = select i1 %36, i1 true, i1 %37
%39 = fcmp ogt float %28, %13
%40 = select i1 %38, i1 true, i1 %39
br i1 %40, label %89, label %41
41: ; preds = %35
%42 = fpext float %28 to double
%43 = tail call double @llvm.floor.f64(double %42)
%44 = fptosi double %43 to i32
%45 = srem i32 %44, %5
%46 = tail call double @llvm.ceil.f64(double %42)
%47 = fptosi double %46 to i32
%48 = srem i32 %47, %5
%49 = sitofp i32 %45 to float
%50 = fsub float %28, %49
%51 = fpext float %33 to double
%52 = tail call double @llvm.floor.f64(double %51)
%53 = fptosi double %52 to i32
%54 = srem i32 %53, %4
%55 = tail call double @llvm.ceil.f64(double %51)
%56 = fptosi double %55 to i32
%57 = srem i32 %56, %4
%58 = sitofp i32 %54 to float
%59 = fsub float %33, %58
%60 = mul nsw i32 %45, %4
%61 = add nsw i32 %54, %60
%62 = sext i32 %61 to i64
%63 = getelementptr inbounds float, float* %0, i64 %62
%64 = load float, float* %63
%65 = fpext float %64 to double
%66 = mul nsw i32 %48, %4
%67 = add nsw i32 %66, %54
%68 = sext i32 %67 to i64
%69 = getelementptr inbounds float, float* %0, i64 %68
%70 = load float, float* %69
%71 = fpext float %70 to double
%72 = fpext float %50 to double
%73 = tail call float (double, double, double, ...) bitcast (float (...)* @linear_interp to float (double, double, double, ...)*)(double %65, double %71, double %72)
%74 = fpext float %73 to double
%75 = add nsw i32 %57, %60
%76 = sext i32 %75 to i64
%77 = getelementptr inbounds float, float* %0, i64 %76
%78 = load float, float* %77
%79 = fpext float %78 to double
%80 = add nsw i32 %57, %66
%81 = sext i32 %80 to i64
%82 = getelementptr inbounds float, float* %0, i64 %81
%83 = load float, float* %82
%84 = fpext float %83 to double
%85 = tail call float (double, double, double, ...) bitcast (float (...)* @linear_interp to float (double, double, double, ...)*)(double %79, double %84, double %72)
%86 = fpext float %85 to double
%87 = fpext float %59 to double
%88 = tail call float (double, double, double, ...) bitcast (float (...)* @linear_interp to float (double, double, double, ...)*)(double %74, double %86, double %87)
br label %89
89: ; preds = %23, %35, %41
%90 = phi float [ %88, %41 ], [ 0.000000e+00, %35 ], [ 0.000000e+00, %23 ]
%91 = getelementptr inbounds float, float* %1, i64 %25
store float %90, float* %91
%92 = add nuw nsw i64 %24, 1
%93 = icmp eq i64 %92, %16
br i1 %93, label %94, label %23
94: ; preds = %89, %17
%95 = add nuw nsw i64 %18, 1
%96 = icmp eq i64 %95, %15
br i1 %96, label %97, label %17
97: ; preds = %94, %6
ret void
}