int stretch(the_image, out_image, x_stretch, y_stretch, bilinear, out_rows,
            out_cols, in_rows, in_cols)
float x_stretch, y_stretch;
int bilinear;
short **the_image, **out_image;
long out_cols, out_rows;
long in_cols, in_rows;
{
  double tmpx, tmpy;
  float fi, fj;
  int i, j, new_i, new_j;
  printf("\n");
  for (i = 0; i < out_rows; i++) {
    if ((i % 10) == 0)
      printf("%d ", i);
    for (j = 0; j < out_cols; j++) {
      fi = i;
      fj = j;
      tmpx = fj / x_stretch;
      tmpy = fi / y_stretch;
      new_i = tmpy;
      new_j = tmpx;
      if (bilinear == 0) {
        if (new_j < 0 || new_j >= in_cols || new_i < 0 || new_i >= in_rows)
          out_image[i][j] = 150;
        else
          out_image[i][j] = the_image[new_i][new_j];
      } else {
        out_image[i][j] =
            bilinear_interpolate(the_image, tmpx, tmpy, in_rows, in_cols);
      }
    }
  }
  return (1);
}


x86 assembly (clang, -O0):
<stretch>:
push   %rbp
mov    %rsp,%rbp
sub    $0x70,%rsp
movsd  %xmm1,-0x68(%rbp)
movaps %xmm0,%xmm1
movsd  -0x68(%rbp),%xmm0
mov    0x10(%rbp),%rax
cvtsd2ss %xmm1,%xmm1
cvtsd2ss %xmm0,%xmm0
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
movss  %xmm1,-0x14(%rbp)
movss  %xmm0,-0x18(%rbp)
mov    %edx,-0x1c(%rbp)
mov    %rcx,-0x28(%rbp)
mov    %r8,-0x30(%rbp)
mov    %r9,-0x38(%rbp)
lea    0xe97(%rip),%rdi
mov    $0x0,%al
call   1040 <printf@plt>
movl   $0x0,-0x54(%rbp)
movslq -0x54(%rbp),%rax
cmp    -0x28(%rbp),%rax
jge    12e5 <stretch+0x1c5>
mov    -0x54(%rbp),%eax
mov    $0xa,%ecx
cltd
idiv   %ecx
cmp    $0x0,%edx
jne    11aa <stretch+0x8a>
mov    -0x54(%rbp),%esi
lea    0xe5f(%rip),%rdi
mov    $0x0,%al
call   1040 <printf@plt>
movl   $0x0,-0x58(%rbp)
movslq -0x58(%rbp),%rax
cmp    -0x30(%rbp),%rax
jge    12d2 <stretch+0x1b2>
cvtsi2ssl -0x54(%rbp),%xmm0
movss  %xmm0,-0x4c(%rbp)
cvtsi2ssl -0x58(%rbp),%xmm0
movss  %xmm0,-0x50(%rbp)
movss  -0x50(%rbp),%xmm0
divss  -0x14(%rbp),%xmm0
cvtss2sd %xmm0,%xmm0
movsd  %xmm0,-0x40(%rbp)
movss  -0x4c(%rbp),%xmm0
divss  -0x18(%rbp),%xmm0
cvtss2sd %xmm0,%xmm0
movsd  %xmm0,-0x48(%rbp)
cvttsd2si -0x48(%rbp),%eax
mov    %eax,-0x5c(%rbp)
cvttsd2si -0x40(%rbp),%eax
mov    %eax,-0x60(%rbp)
cmpl   $0x0,-0x1c(%rbp)
jne    128b <stretch+0x16b>
cmpl   $0x0,-0x60(%rbp)
jl     1243 <stretch+0x123>
movslq -0x60(%rbp),%rax
cmp    0x10(%rbp),%rax
jge    1243 <stretch+0x123>
cmpl   $0x0,-0x5c(%rbp)
jl     1243 <stretch+0x123>
movslq -0x5c(%rbp),%rax
cmp    -0x38(%rbp),%rax
jl     125e <stretch+0x13e>
mov    -0x10(%rbp),%rax
movslq -0x54(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x58(%rbp),%rcx
movw   $0x96,(%rax,%rcx,2)
jmp    1286 <stretch+0x166>
mov    -0x8(%rbp),%rax
movslq -0x5c(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x60(%rbp),%rcx
mov    (%rax,%rcx,2),%dx
mov    -0x10(%rbp),%rax
movslq -0x54(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x58(%rbp),%rcx
mov    %dx,(%rax,%rcx,2)
jmp    12bf <stretch+0x19f>
mov    -0x8(%rbp),%rdi
movsd  -0x40(%rbp),%xmm0
movsd  -0x48(%rbp),%xmm1
mov    -0x38(%rbp),%rsi
mov    0x10(%rbp),%rdx
mov    $0x2,%al
call   1030 <bilinear_interpolate@plt>
mov    %ax,%dx
mov    -0x10(%rbp),%rax
movslq -0x54(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x58(%rbp),%rcx
mov    %dx,(%rax,%rcx,2)
jmp    12c4 <stretch+0x1a4>
mov    -0x58(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x58(%rbp)
jmp    11b1 <stretch+0x91>
jmp    12d7 <stretch+0x1b7>
mov    -0x54(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x54(%rbp)
jmp    1177 <stretch+0x57>
mov    $0x1,%eax
add    $0x70,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<stretch>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x38,%rsp
mov    %r9,0x28(%rsp)
mov    %r8,%rbx
mov    %rcx,%r14
mov    %edx,%r12d
movsd  %xmm1,0x8(%rsp)
movsd  %xmm0,0x10(%rsp)
mov    %rsi,%rbp
mov    %rdi,0x20(%rsp)
mov    $0xa,%edi
call   1030 <putchar@plt>
mov    %r14,0x30(%rsp)
test   %r14,%r14
jle    128a <stretch+0x15a>
movsd  0x10(%rsp),%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,0x10(%rsp)
movsd  0x8(%rsp),%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,0x1c(%rsp)
xor    %r15d,%r15d
jmp    11af <stretch+0x7f>
nopl   (%rax)
add    $0x1,%r15
cmp    0x30(%rsp),%r15
je     128a <stretch+0x15a>
imul   $0xcccccccd,%r15d,%eax
ror    %eax
cmp    $0x19999999,%eax
jbe    11d0 <stretch+0xa0>
test   %rbx,%rbx
jg     11e6 <stretch+0xb6>
jmp    11a0 <stretch+0x70>
cs nopw 0x0(%rax,%rax,1)
lea    0xe29(%rip),%rdi
mov    %r15d,%esi
xor    %eax,%eax
call   1050 <printf@plt>
test   %rbx,%rbx
jle    11a0 <stretch+0x70>
xorps  %xmm0,%xmm0
cvtsi2ss %r15d,%xmm0
divss  0x1c(%rsp),%xmm0
cvttss2si %xmm0,%r13d
cvtss2sd %xmm0,%xmm0
movsd  %xmm0,0x8(%rsp)
xor    %r14d,%r14d
jmp    1247 <stretch+0x117>
nopl   0x0(%rax,%rax,1)
cvtss2sd %xmm0,%xmm0
mov    0x20(%rsp),%rdi
movsd  0x8(%rsp),%xmm1
mov    0x28(%rsp),%rsi
mov    0x70(%rsp),%rdx
mov    $0x2,%al
call   1040 <bilinear_interpolate@plt>
mov    0x0(%rbp,%r15,8),%rcx
mov    %ax,(%rcx,%r14,2)
add    $0x1,%r14
cmp    %r14,%rbx
je     11a0 <stretch+0x70>
xorps  %xmm0,%xmm0
cvtsi2ss %r14d,%xmm0
divss  0x10(%rsp),%xmm0
test   %r12d,%r12d
jne    1210 <stretch+0xe0>
cvttss2si %xmm0,%ecx
mov    $0x96,%ax
test   %ecx,%ecx
js     1230 <stretch+0x100>
mov    %ecx,%ecx
cmp    0x70(%rsp),%rcx
jge    1230 <stretch+0x100>
test   %r13d,%r13d
js     1230 <stretch+0x100>
cmp    0x28(%rsp),%r13
jge    1230 <stretch+0x100>
mov    0x20(%rsp),%rax
mov    (%rax,%r13,8),%rax
movzwl (%rax,%rcx,2),%eax
jmp    1230 <stretch+0x100>
mov    $0x1,%eax
add    $0x38,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

LLVM IR (-O2):
@.str.1 = [4 x i8] c"%d \00"
define i32 @stretch(i16** %0, i16** nocapture %1, double %2, double %3, i32 %4, i64 %5, i64 %6, i64 %7, i64 %8) {
%10 = fptrunc double %2 to float
%11 = fptrunc double %3 to float
%12 = tail call i32 @putchar(i32 10)
%13 = icmp sgt i64 %5, 0
br i1 %13, label %14, label %65
14: ; preds = %9
%15 = icmp sgt i64 %6, 0
%16 = icmp eq i32 %4, 0
br label %17
17: ; preds = %14, %62
%18 = phi i64 [ 0, %14 ], [ %63, %62 ]
%19 = trunc i64 %18 to i32
%20 = urem i32 %19, 10
%21 = icmp eq i32 %20, 0
br i1 %21, label %22, label %24
22: ; preds = %17
%23 = tail call i32 (i8*, ...) @printf(i8* dereferenceable(1) getelementptr inbounds ([4 x i8], [4 x i8]* @.str.1, i64 0, i64 0), i32 %19)
br label %24
24: ; preds = %22, %17
br i1 %15, label %25, label %62
25: ; preds = %24
%26 = sitofp i32 %19 to float
%27 = fdiv float %26, %11
%28 = fptosi float %27 to i32
%29 = fpext float %27 to double
%30 = getelementptr inbounds i16*, i16** %1, i64 %18
%31 = icmp slt i32 %28, 0
%32 = zext i32 %28 to i64
%33 = icmp sge i64 %32, %7
%34 = getelementptr inbounds i16*, i16** %0, i64 %32
br label %35
35: ; preds = %25, %56
%36 = phi i64 [ 0, %25 ], [ %60, %56 ]
%37 = trunc i64 %36 to i32
%38 = sitofp i32 %37 to float
%39 = fdiv float %38, %10
%40 = fptosi float %39 to i32
br i1 %16, label %41, label %52
41: ; preds = %35
%42 = icmp slt i32 %40, 0
br i1 %42, label %56, label %43
43: ; preds = %41
%44 = zext i32 %40 to i64
%45 = icmp sge i64 %44, %8
%46 = select i1 %45, i1 true, i1 %31
%47 = select i1 %46, i1 true, i1 %33
br i1 %47, label %56, label %48
48: ; preds = %43
%49 = load i16*, i16** %34
%50 = getelementptr inbounds i16, i16* %49, i64 %44
%51 = load i16, i16* %50
br label %56
52: ; preds = %35
%53 = fpext float %39 to double
%54 = tail call i32 (i16**, double, double, i64, i64, ...) bitcast (i32 (...)* @bilinear_interpolate to i32 (i16**, double, double, i64, i64, ...)*)(i16** %0, double %53, double %29, i64 %7, i64 %8)
%55 = trunc i32 %54 to i16
br label %56
56: ; preds = %41, %43, %52, %48
%57 = phi i16 [ %55, %52 ], [ %51, %48 ], [ 150, %43 ], [ 150, %41 ]
%58 = load i16*, i16** %30
%59 = getelementptr inbounds i16, i16* %58, i64 %36
store i16 %57, i16* %59
%60 = add nuw nsw i64 %36, 1
%61 = icmp eq i64 %60, %6
br i1 %61, label %62, label %35
62: ; preds = %56, %24
%63 = add nuw nsw i64 %18, 1
%64 = icmp eq i64 %63, %5
br i1 %64, label %65, label %17
65: ; preds = %62, %9
ret i32 1
}