#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 /*<<< orphan*/ fill_poles(double *, int);
int /*<<< orphan*/ invspline1D(float *, int, int, double *, int);
int /*<<< orphan*/ isfinite(float);
bool prepare_spline(float *img, int w, int h, int pd, int order) {
  if (order < 3)
    return true;
  for (int i = 0; i < w * h * pd; i++)
    if (!isfinite(img[i]))
      img[i] = 0;
  double z[5];
  if (!fill_poles(z, order))
    return false;
  int npoles = order / 2;
  for (int k = 0; k < pd; k++) {
    for (int y = 0; y < h; y++)
      invspline1D(img + (y * w + 0) * pd + k, pd * 1, w, z, npoles);
    for (int x = 0; x < w; x++)
      invspline1D(img + (0 * w + x) * pd + k, pd * w, h, z, npoles);
  }
  return true;
}


x86 assembly (clang, -O0):
<prepare_spline>:
push   %rbp
mov    %rsp,%rbp
sub    $0x60,%rsp
mov    %rdi,-0x10(%rbp)
mov    %esi,-0x14(%rbp)
mov    %edx,-0x18(%rbp)
mov    %ecx,-0x1c(%rbp)
mov    %r8d,-0x20(%rbp)
cmpl   $0x3,-0x20(%rbp)
jge    115f <prepare_spline+0x2f>
movl   $0x1,-0x4(%rbp)
jmp    12ce <prepare_spline+0x19e>
movl   $0x0,-0x24(%rbp)
mov    -0x24(%rbp),%eax
mov    -0x14(%rbp),%ecx
imul   -0x18(%rbp),%ecx
imul   -0x1c(%rbp),%ecx
cmp    %ecx,%eax
jge    11ba <prepare_spline+0x8a>
mov    -0x10(%rbp),%rax
movslq -0x24(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
call   1040 <isfinite@plt>
cmp    $0x0,%eax
jne    11a7 <prepare_spline+0x77>
mov    -0x10(%rbp),%rax
movslq -0x24(%rbp),%rcx
xorps  %xmm0,%xmm0
movss  %xmm0,(%rax,%rcx,4)
jmp    11ac <prepare_spline+0x7c>
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    1166 <prepare_spline+0x36>
lea    -0x50(%rbp),%rdi
mov    -0x20(%rbp),%esi
call   1050 <fill_poles@plt>
cmp    $0x0,%eax
jne    11db <prepare_spline+0xab>
movl   $0x0,-0x4(%rbp)
jmp    12ce <prepare_spline+0x19e>
mov    -0x20(%rbp),%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
mov    %eax,-0x54(%rbp)
movl   $0x0,-0x58(%rbp)
mov    -0x58(%rbp),%eax
cmp    -0x1c(%rbp),%eax
jge    12c7 <prepare_spline+0x197>
movl   $0x0,-0x5c(%rbp)
mov    -0x5c(%rbp),%eax
cmp    -0x18(%rbp),%eax
jge    1259 <prepare_spline+0x129>
mov    -0x10(%rbp),%rdi
mov    -0x5c(%rbp),%eax
imul   -0x14(%rbp),%eax
add    $0x0,%eax
imul   -0x1c(%rbp),%eax
cltq
shl    $0x2,%rax
add    %rax,%rdi
movslq -0x58(%rbp),%rax
shl    $0x2,%rax
add    %rax,%rdi
mov    -0x1c(%rbp),%esi
shl    $0x0,%esi
mov    -0x14(%rbp),%edx
lea    -0x50(%rbp),%rcx
mov    -0x54(%rbp),%r8d
call   1030 <invspline1D@plt>
mov    -0x5c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x5c(%rbp)
jmp    1203 <prepare_spline+0xd3>
movl   $0x0,-0x60(%rbp)
mov    -0x60(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    12b4 <prepare_spline+0x184>
mov    -0x10(%rbp),%rdi
imul   $0x0,-0x14(%rbp),%eax
add    -0x60(%rbp),%eax
imul   -0x1c(%rbp),%eax
cltq
shl    $0x2,%rax
add    %rax,%rdi
movslq -0x58(%rbp),%rax
shl    $0x2,%rax
add    %rax,%rdi
mov    -0x1c(%rbp),%esi
imul   -0x14(%rbp),%esi
mov    -0x18(%rbp),%edx
lea    -0x50(%rbp),%rcx
mov    -0x54(%rbp),%r8d
call   1030 <invspline1D@plt>
mov    -0x60(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x60(%rbp)
jmp    1260 <prepare_spline+0x130>
jmp    12b9 <prepare_spline+0x189>
mov    -0x58(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x58(%rbp)
jmp    11f0 <prepare_spline+0xc0>
movl   $0x1,-0x4(%rbp)
mov    -0x4(%rbp),%eax
add    $0x60,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<prepare_spline>:
mov    $0x1,%eax
cmp    $0x3,%r8d
jl     1311 <prepare_spline+0x1e1>
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x88,%rsp
mov    %r8d,%r15d
mov    %ecx,%r14d
mov    %esi,%r13d
mov    %rdi,%r12
mov    %ecx,%eax
imul   %esi,%eax
mov    %eax,0x24(%rsp)
mov    %edx,0xc(%rsp)
imul   %edx,%eax
test   %eax,%eax
jle    11a2 <prepare_spline+0x72>
mov    %eax,%ebx
xor    %ebp,%ebp
jmp    1189 <prepare_spline+0x59>
cs nopw 0x0(%rax,%rax,1)
add    $0x1,%rbp
cmp    %rbp,%rbx
je     11a2 <prepare_spline+0x72>
movss  (%r12,%rbp,4),%xmm0
call   1040 <isfinite@plt>
test   %eax,%eax
jne    1180 <prepare_spline+0x50>
movl   $0x0,(%r12,%rbp,4)
jmp    1180 <prepare_spline+0x50>
lea    0x60(%rsp),%rdi
mov    %r15d,%esi
call   1050 <fill_poles@plt>
test   %eax,%eax
je     12f7 <prepare_spline+0x1c7>
mov    $0x1,%eax
test   %r14d,%r14d
jle    1300 <prepare_spline+0x1d0>
mov    %r13d,%eax
mov    %r15d,%r13d
shr    $0x1f,%r13d
add    %r15d,%r13d
sar    %r13d
mov    %r14d,0x20(%rsp)
mov    %r14d,%ecx
mov    0xc(%rsp),%r15d
mov    %r15d,%edx
mov    %rdx,0x30(%rsp)
mov    %eax,0x1c(%rsp)
mov    %eax,%eax
mov    %rax,0x48(%rsp)
mov    %rcx,0x50(%rsp)
lea    0x0(,%rcx,4),%rax
mov    %rax,0x40(%rsp)
xor    %eax,%eax
mov    %rax,0x10(%rsp)
lea    0x60(%rsp),%r8
mov    %r12,0x38(%rsp)
mov    %r12,0x28(%rsp)
mov    0x24(%rsp),%r12d
jmp    1252 <prepare_spline+0x122>
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
mov    0x10(%rsp),%rcx
add    $0x1,%rcx
addq   $0x4,0x28(%rsp)
mov    %rcx,%rax
mov    %rcx,0x10(%rsp)
cmp    0x50(%rsp),%rcx
je     12fb <prepare_spline+0x1cb>
test   %r15d,%r15d
mov    0x20(%rsp),%r14d
mov    0x1c(%rsp),%r15d
jle    12ab <prepare_spline+0x17b>
mov    0x38(%rsp),%rax
mov    0x10(%rsp),%rcx
lea    (%rax,%rcx,4),%rax
mov    %rax,0x58(%rsp)
xor    %ebp,%ebp
mov    0x30(%rsp),%rbx
nopl   0x0(%rax,%rax,1)
movslq %ebp,%rbp
mov    0x58(%rsp),%rax
lea    (%rax,%rbp,4),%rdi
mov    %r14d,%esi
mov    %r15d,%edx
mov    %r8,%rcx
mov    %r13d,%r8d
call   1030 <invspline1D@plt>
lea    0x60(%rsp),%r8
add    %r12d,%ebp
add    $0xffffffffffffffff,%rbx
jne    1280 <prepare_spline+0x150>
mov    0x48(%rsp),%rbx
mov    0x28(%rsp),%r14
test   %r15d,%r15d
mov    0xc(%rsp),%r15d
mov    0x40(%rsp),%rbp
jle    1230 <prepare_spline+0x100>
nopl   0x0(%rax,%rax,1)
mov    %r14,%rdi
mov    %r12d,%esi
mov    %r15d,%edx
mov    %r8,%rcx
mov    %r13d,%r8d
call   1030 <invspline1D@plt>
lea    0x60(%rsp),%r8
add    %rbp,%r14
add    $0xffffffffffffffff,%rbx
jne    12d0 <prepare_spline+0x1a0>
jmp    1230 <prepare_spline+0x100>
xor    %eax,%eax
jmp    1300 <prepare_spline+0x1d0>
mov    $0x1,%eax
add    $0x88,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

LLVM IR (-O0):
define i32 @prepare_spline(float* %0, i32 %1, i32 %2, i32 %3, i32 %4) {
%6 = alloca i32
%7 = alloca float*
%8 = alloca i32
%9 = alloca i32
%10 = alloca i32
%11 = alloca i32
%12 = alloca i32
%13 = alloca [5 x double]
%14 = alloca i32
%15 = alloca i32
%16 = alloca i32
%17 = alloca i32
store float* %0, float** %7
store i32 %1, i32* %8
store i32 %2, i32* %9
store i32 %3, i32* %10
store i32 %4, i32* %11
%18 = load i32, i32* %11
%19 = icmp slt i32 %18, 3
br i1 %19, label %20, label %21
20: ; preds = %5
store i32 1, i32* %6
br label %120
21: ; preds = %5
store i32 0, i32* %12
br label %22
22: ; preds = %44, %21
%23 = load i32, i32* %12
%24 = load i32, i32* %8
%25 = load i32, i32* %9
%26 = mul nsw i32 %24, %25
%27 = load i32, i32* %10
%28 = mul nsw i32 %26, %27
%29 = icmp slt i32 %23, %28
br i1 %29, label %30, label %47
30: ; preds = %22
%31 = load float*, float** %7
%32 = load i32, i32* %12
%33 = sext i32 %32 to i64
%34 = getelementptr inbounds float, float* %31, i64 %33
%35 = load float, float* %34
%36 = call i32 @isfinite(float %35)
%37 = icmp ne i32 %36, 0
br i1 %37, label %43, label %38
38: ; preds = %30
%39 = load float*, float** %7
%40 = load i32, i32* %12
%41 = sext i32 %40 to i64
%42 = getelementptr inbounds float, float* %39, i64 %41
store float 0.000000e+00, float* %42
br label %43
43: ; preds = %38, %30
br label %44
44: ; preds = %43
%45 = load i32, i32* %12
%46 = add nsw i32 %45, 1
store i32 %46, i32* %12
br label %22
47: ; preds = %22
%48 = getelementptr inbounds [5 x double], [5 x double]* %13, i64 0, i64 0
%49 = load i32, i32* %11
%50 = call i32 @fill_poles(double* %48, i32 %49)
%51 = icmp ne i32 %50, 0
br i1 %51, label %53, label %52
52: ; preds = %47
store i32 0, i32* %6
br label %120
53: ; preds = %47
%54 = load i32, i32* %11
%55 = sdiv i32 %54, 2
store i32 %55, i32* %14
store i32 0, i32* %15
br label %56
56: ; preds = %116, %53
%57 = load i32, i32* %15
%58 = load i32, i32* %10
%59 = icmp slt i32 %57, %58
br i1 %59, label %60, label %119
60: ; preds = %56
store i32 0, i32* %16
br label %61
61: ; preds = %84, %60
%62 = load i32, i32* %16
%63 = load i32, i32* %9
%64 = icmp slt i32 %62, %63
br i1 %64, label %65, label %87
65: ; preds = %61
%66 = load float*, float** %7
%67 = load i32, i32* %16
%68 = load i32, i32* %8
%69 = mul nsw i32 %67, %68
%70 = add nsw i32 %69, 0
%71 = load i32, i32* %10
%72 = mul nsw i32 %70, %71
%73 = sext i32 %72 to i64
%74 = getelementptr inbounds float, float* %66, i64 %73
%75 = load i32, i32* %15
%76 = sext i32 %75 to i64
%77 = getelementptr inbounds float, float* %74, i64 %76
%78 = load i32, i32* %10
%79 = mul nsw i32 %78, 1
%80 = load i32, i32* %8
%81 = getelementptr inbounds [5 x double], [5 x double]* %13, i64 0, i64 0
%82 = load i32, i32* %14
%83 = call i32 @invspline1D(float* %77, i32 %79, i32 %80, double* %81, i32 %82)
br label %84
84: ; preds = %65
%85 = load i32, i32* %16
%86 = add nsw i32 %85, 1
store i32 %86, i32* %16
br label %61
87: ; preds = %61
store i32 0, i32* %17
br label %88
88: ; preds = %112, %87
%89 = load i32, i32* %17
%90 = load i32, i32* %8
%91 = icmp slt i32 %89, %90
br i1 %91, label %92, label %115
92: ; preds = %88
%93 = load float*, float** %7
%94 = load i32, i32* %8
%95 = mul nsw i32 0, %94
%96 = load i32, i32* %17
%97 = add nsw i32 %95, %96
%98 = load i32, i32* %10
%99 = mul nsw i32 %97, %98
%100 = sext i32 %99 to i64
%101 = getelementptr inbounds float, float* %93, i64 %100
%102 = load i32, i32* %15
%103 = sext i32 %102 to i64
%104 = getelementptr inbounds float, float* %101, i64 %103
%105 = load i32, i32* %10
%106 = load i32, i32* %8
%107 = mul nsw i32 %105, %106
%108 = load i32, i32* %9
%109 = getelementptr inbounds [5 x double], [5 x double]* %13, i64 0, i64 0
%110 = load i32, i32* %14
%111 = call i32 @invspline1D(float* %104, i32 %107, i32 %108, double* %109, i32 %110)
br label %112
112: ; preds = %92
%113 = load i32, i32* %17
%114 = add nsw i32 %113, 1
store i32 %114, i32* %17
br label %88
115: ; preds = %88
br label %116
116: ; preds = %115
%117 = load i32, i32* %15
%118 = add nsw i32 %117, 1
store i32 %118, i32* %15
br label %56
119: ; preds = %56
store i32 1, i32* %6
br label %120
120: ; preds = %119, %52, %20
%121 = load i32, i32* %6
ret i32 %121
}

LLVM IR (-O2):
define i32 @prepare_spline(float* %0, i32 %1, i32 %2, i32 %3, i32 %4) {
%6 = alloca [5 x double]
%7 = icmp slt i32 %4, 3
br i1 %7, label %67, label %8
8: ; preds = %5
%9 = mul i32 %3, %1
%10 = mul i32 %9, %2
%11 = icmp sgt i32 %10, 0
br i1 %11, label %12, label %14
12: ; preds = %8
%13 = zext i32 %10 to i64
br label %19
14: ; preds = %26, %8
%15 = bitcast [5 x double]* %6 to i8*
%16 = getelementptr inbounds [5 x double], [5 x double]* %6, i64 0, i64 0
%17 = call i32 @fill_poles(double* %16, i32 %4)
%18 = icmp eq i32 %17, 0
br i1 %18, label %65, label %29
19: ; preds = %12, %26
%20 = phi i64 [ 0, %12 ], [ %27, %26 ]
%21 = getelementptr inbounds float, float* %0, i64 %20
%22 = load float, float* %21
%23 = tail call i32 @isfinite(float %22)
%24 = icmp eq i32 %23, 0
br i1 %24, label %25, label %26
25: ; preds = %19
store float 0.000000e+00, float* %21
br label %26
26: ; preds = %19, %25
%27 = add nuw nsw i64 %20, 1
%28 = icmp eq i64 %27, %13
br i1 %28, label %14, label %19
29: ; preds = %14
%30 = sdiv i32 %4, 2
%31 = icmp sgt i32 %3, 0
br i1 %31, label %32, label %65
32: ; preds = %29
%33 = icmp sgt i32 %2, 0
%34 = icmp sgt i32 %1, 0
%35 = zext i32 %3 to i64
%36 = zext i32 %3 to i64
%37 = zext i32 %2 to i64
%38 = zext i32 %1 to i64
br label %39
39: ; preds = %32, %55
%40 = phi i64 [ 0, %32 ], [ %56, %55 ]
br i1 %33, label %41, label %43
41: ; preds = %39
%42 = getelementptr inbounds float, float* %0, i64 %40
br label %46
43: ; preds = %46, %39
br i1 %34, label %44, label %55
44: ; preds = %43
%45 = getelementptr inbounds float, float* %0, i64 %40
br label %58
46: ; preds = %41, %46
%47 = phi i64 [ 0, %41 ], [ %53, %46 ]
%48 = trunc i64 %47 to i32
%49 = mul i32 %9, %48
%50 = sext i32 %49 to i64
%51 = getelementptr inbounds float, float* %42, i64 %50
%52 = call i32 @invspline1D(float* %51, i32 %3, i32 %1, double* %16, i32 %30)
%53 = add nuw nsw i64 %47, 1
%54 = icmp eq i64 %53, %37
br i1 %54, label %43, label %46
55: ; preds = %58, %43
%56 = add nuw nsw i64 %40, 1
%57 = icmp eq i64 %56, %36
br i1 %57, label %65, label %39
58: ; preds = %44, %58
%59 = phi i64 [ 0, %44 ], [ %63, %58 ]
%60 = mul nsw i64 %59, %35
%61 = getelementptr inbounds float, float* %45, i64 %60
%62 = call i32 @invspline1D(float* %61, i32 %9, i32 %2, double* %16, i32 %30)
%63 = add nuw nsw i64 %59, 1
%64 = icmp eq i64 %63, %38
br i1 %64, label %55, label %58
65: ; preds = %55, %29, %14
%66 = phi i32 [ 0, %14 ], [ 1, %29 ], [ 1, %55 ]
br label %67
67: ; preds = %5, %65
%68 = phi i32 [ %66, %65 ], [ 1, %5 ]
ret i32 %68
}