#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*/ bacdic(float *, int, double *, int, float *, int);
double *dvector(int, int);
int /*<<< orphan*/ free_dvector(double *, int);
int /*<<< orphan*/ invrse(float *, int);
int /*<<< orphan*/ memcof(float *, int, int, double *, double *);
void wafil3(float *x, int n, int kc, float *fcv, float *tmp) {
  double *cof, dum;
  int i, j, nd, m, off1, off2;
  for (i = kc + 1; i <= n - kc; i++) {
    tmp[i] = 0.0;
    for (j = 1; j <= 2 * kc + 1; j++)
      tmp[i] += fcv[j] * x[i + j - kc - 1];
  }
  nd = (((n - kc - kc) < (kc * 2)) ? (n - kc - kc) : (kc * 2));
  j = (((nd) < (kc)) ? (nd) : (kc));
  m = (((256) < (j)) ? (256) : (j));
  cof = dvector(1, m);
  off1 = n - nd - kc;
  off2 = n - kc;
  memcof(tmp + off1, nd, m, &dum, cof);
  bacdic(tmp + off1, nd, cof, m, tmp + off2, kc);
  invrse(tmp + 1, n);
  memcof(tmp + off1, nd, m, &dum, cof);
  bacdic(tmp + off1, nd, cof, m, tmp + off2, kc);
  invrse(tmp + 1, n);
  for (i = 1; i <= n; i++)
    x[i] -= tmp[i];
  free_dvector(cof, 1);
}


x86 assembly (clang, -O0):
<wafil3>:
push   %rbp
mov    %rsp,%rbp
sub    $0x60,%rsp
mov    %rdi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
mov    %edx,-0x10(%rbp)
mov    %rcx,-0x18(%rbp)
mov    %r8,-0x20(%rbp)
mov    -0x10(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
mov    -0x34(%rbp),%eax
mov    -0xc(%rbp),%ecx
sub    -0x10(%rbp),%ecx
cmp    %ecx,%eax
jg     120f <wafil3+0xbf>
mov    -0x20(%rbp),%rax
movslq -0x34(%rbp),%rcx
xorps  %xmm0,%xmm0
movss  %xmm0,(%rax,%rcx,4)
movl   $0x1,-0x38(%rbp)
mov    -0x38(%rbp),%eax
mov    -0x10(%rbp),%ecx
shl    $0x1,%ecx
add    $0x1,%ecx
cmp    %ecx,%eax
jg     11fc <wafil3+0xac>
mov    -0x18(%rbp),%rax
movslq -0x38(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x34(%rbp),%ecx
add    -0x38(%rbp),%ecx
sub    -0x10(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm2
mov    -0x20(%rbp),%rax
movslq -0x34(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,(%rax,%rcx,4)
mov    -0x38(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x38(%rbp)
jmp    119b <wafil3+0x4b>
jmp    1201 <wafil3+0xb1>
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    1173 <wafil3+0x23>
mov    -0xc(%rbp),%eax
sub    -0x10(%rbp),%eax
sub    -0x10(%rbp),%eax
mov    -0x10(%rbp),%ecx
shl    $0x1,%ecx
cmp    %ecx,%eax
jge    1237 <wafil3+0xe7>
mov    -0xc(%rbp),%eax
sub    -0x10(%rbp),%eax
sub    -0x10(%rbp),%eax
mov    %eax,-0x4c(%rbp)
jmp    1240 <wafil3+0xf0>
mov    -0x10(%rbp),%eax
shl    $0x1,%eax
mov    %eax,-0x4c(%rbp)
mov    -0x4c(%rbp),%eax
mov    %eax,-0x3c(%rbp)
mov    -0x3c(%rbp),%eax
cmp    -0x10(%rbp),%eax
jge    125d <wafil3+0x10d>
mov    -0x3c(%rbp),%eax
mov    %eax,-0x50(%rbp)
jmp    1263 <wafil3+0x113>
mov    -0x10(%rbp),%eax
mov    %eax,-0x50(%rbp)
mov    -0x50(%rbp),%eax
mov    %eax,-0x38(%rbp)
mov    $0x100,%eax
cmp    -0x38(%rbp),%eax
jge    1284 <wafil3+0x134>
mov    $0x100,%eax
mov    %eax,-0x54(%rbp)
jmp    128a <wafil3+0x13a>
mov    -0x38(%rbp),%eax
mov    %eax,-0x54(%rbp)
mov    -0x54(%rbp),%eax
mov    %eax,-0x40(%rbp)
mov    -0x40(%rbp),%esi
mov    $0x1,%edi
call   1040 <dvector@plt>
mov    %rax,-0x28(%rbp)
mov    -0xc(%rbp),%eax
sub    -0x3c(%rbp),%eax
sub    -0x10(%rbp),%eax
mov    %eax,-0x44(%rbp)
mov    -0xc(%rbp),%eax
sub    -0x10(%rbp),%eax
mov    %eax,-0x48(%rbp)
mov    -0x20(%rbp),%rdi
movslq -0x44(%rbp),%rax
shl    $0x2,%rax
add    %rax,%rdi
mov    -0x3c(%rbp),%esi
mov    -0x40(%rbp),%edx
mov    -0x28(%rbp),%r8
lea    -0x30(%rbp),%rcx
call   1050 <memcof@plt>
mov    -0x20(%rbp),%rdi
movslq -0x44(%rbp),%rax
shl    $0x2,%rax
add    %rax,%rdi
mov    -0x3c(%rbp),%esi
mov    -0x28(%rbp),%rdx
mov    -0x40(%rbp),%ecx
mov    -0x20(%rbp),%r8
movslq -0x48(%rbp),%rax
shl    $0x2,%rax
add    %rax,%r8
mov    -0x10(%rbp),%r9d
call   1030 <bacdic@plt>
mov    -0x20(%rbp),%rdi
add    $0x4,%rdi
mov    -0xc(%rbp),%esi
call   1060 <invrse@plt>
mov    -0x20(%rbp),%rdi
movslq -0x44(%rbp),%rax
shl    $0x2,%rax
add    %rax,%rdi
mov    -0x3c(%rbp),%esi
mov    -0x40(%rbp),%edx
mov    -0x28(%rbp),%r8
lea    -0x30(%rbp),%rcx
call   1050 <memcof@plt>
mov    -0x20(%rbp),%rdi
movslq -0x44(%rbp),%rax
shl    $0x2,%rax
add    %rax,%rdi
mov    -0x3c(%rbp),%esi
mov    -0x28(%rbp),%rdx
mov    -0x40(%rbp),%ecx
mov    -0x20(%rbp),%r8
movslq -0x48(%rbp),%rax
shl    $0x2,%rax
add    %rax,%r8
mov    -0x10(%rbp),%r9d
call   1030 <bacdic@plt>
mov    -0x20(%rbp),%rdi
add    $0x4,%rdi
mov    -0xc(%rbp),%esi
call   1060 <invrse@plt>
movl   $0x1,-0x34(%rbp)
mov    -0x34(%rbp),%eax
cmp    -0xc(%rbp),%eax
jg     13c0 <wafil3+0x270>
mov    -0x20(%rbp),%rax
movslq -0x34(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm1
mov    -0x8(%rbp),%rax
movslq -0x34(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
subss  %xmm1,%xmm0
movss  %xmm0,(%rax,%rcx,4)
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    1383 <wafil3+0x233>
mov    -0x28(%rbp),%rdi
mov    $0x1,%esi
call   1070 <free_dvector@plt>
add    $0x60,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<wafil3>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x38,%rsp
mov    %r8,%r14
mov    %rdi,%r15
mov    %esi,0xc(%rsp)
mov    %esi,%eax
sub    %edx,%eax
lea    (%rdx,%rdx,1),%r13d
mov    %eax,%r8d
sub    %edx,%r8d
jle    123f <wafil3+0xef>
mov    %r13d,%r9d
or     $0x1,%r9d
lea    0x2(%r13),%ebp
movslq %edx,%rsi
movslq %eax,%r12
mov    $0x2,%edi
sub    %rbp,%rdi
xor    %eax,%eax
jmp    11c4 <wafil3+0x74>
nopl   0x0(%rax)
movss  (%rcx,%rbx,4),%xmm1
add    %ebx,%r10d
movslq %r10d,%rbp
mulss  (%r15,%rbp,4),%xmm1
addss  %xmm0,%xmm1
movss  %xmm1,(%r14,%rsi,4)
add    $0x1,%rax
cmp    %r12,%rsi
je     1242 <wafil3+0xf2>
mov    %rsi,%r10
add    $0x1,%rsi
movl   $0x0,0x4(%r14,%r10,4)
test   %r9d,%r9d
jle    11bb <wafil3+0x6b>
sub    %edx,%r10d
mov    $0x1,%ebx
xorps  %xmm0,%xmm0
test   %r13d,%r13d
je     11a0 <wafil3+0x50>
nopl   0x0(%rax)
movss  (%rcx,%rbx,4),%xmm1
lea    (%rax,%rbx,1),%rbp
movslq %ebp,%rbp
mulss  (%r15,%rbp,4),%xmm1
addss  %xmm0,%xmm1
movss  %xmm1,(%r14,%rsi,4)
movss  0x4(%rcx,%rbx,4),%xmm0
add    $0x1,%ebp
movslq %ebp,%rbp
mulss  (%r15,%rbp,4),%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,(%r14,%rsi,4)
lea    (%rdi,%rbx,1),%rbp
add    $0x2,%rbp
add    $0x2,%rbx
cmp    $0x1,%rbp
jne    11f0 <wafil3+0xa0>
jmp    11a0 <wafil3+0x50>
movslq %eax,%r12
cmp    %r13d,%r8d
cmovl  %r8d,%r13d
cmp    %edx,%r13d
mov    %edx,%eax
cmovl  %r13d,%eax
cmp    $0x100,%eax
mov    $0x100,%ebx
cmovl  %eax,%ebx
mov    $0x1,%edi
mov    %ebx,%esi
mov    %rdx,%rbp
mov    %rdx,0x10(%rsp)
call   1040 <dvector@plt>
mov    0xc(%rsp),%ecx
sub    %r13d,%ecx
sub    %ebp,%ecx
movslq %ecx,%rcx
lea    (%r14,%rcx,4),%rdi
mov    %rdi,0x18(%rsp)
lea    0x30(%rsp),%rcx
mov    %r13d,%esi
mov    %ebx,%edx
mov    %rax,%rbp
mov    %rax,%r8
call   1050 <memcof@plt>
lea    (%r14,%r12,4),%r8
mov    %r8,0x28(%rsp)
mov    0x18(%rsp),%r12
mov    %r12,%rdi
mov    %r13d,%esi
mov    %rbp,%rdx
mov    %ebx,%ecx
mov    0x10(%rsp),%r9
call   1030 <bacdic@plt>
lea    0x4(%r14),%rdi
mov    %rdi,0x20(%rsp)
mov    0xc(%rsp),%esi
call   1060 <invrse@plt>
mov    %r12,%rdi
mov    %r13d,%esi
mov    %ebx,%edx
lea    0x30(%rsp),%rcx
mov    %rbp,%r8
call   1050 <memcof@plt>
mov    %r12,%rdi
mov    %r13d,%esi
mov    %rbp,%r13
mov    %rbp,%rdx
mov    %ebx,%ecx
mov    0x20(%rsp),%rbx
mov    0x28(%rsp),%r8
mov    0x10(%rsp),%r9
call   1030 <bacdic@plt>
mov    %rbx,%rdi
mov    0xc(%rsp),%ebp
mov    %ebp,%esi
call   1060 <invrse@plt>
test   %ebp,%ebp
jle    14ba <wafil3+0x36a>
mov    %ebp,%edx
lea    0x1(%rdx),%rax
mov    $0x1,%ecx
cmp    $0x8,%ebp
jb     1415 <wafil3+0x2c5>
lea    0x4(%r15),%rsi
lea    (%r14,%rax,4),%rdi
cmp    %rdi,%rsi
jae    134d <wafil3+0x1fd>
lea    (%r15,%rax,4),%rsi
cmp    %rsi,%rbx
jb     1415 <wafil3+0x2c5>
mov    %edx,%ecx
and    $0xfffffff8,%ecx
lea    -0x8(%rcx),%rsi
mov    %rsi,%rbp
shr    $0x3,%rbp
add    $0x1,%rbp
test   %rsi,%rsi
je     14d6 <wafil3+0x386>
mov    %rbp,%rdi
and    $0xfffffffffffffffe,%rdi
mov    $0x1,%esi
cs nopw 0x0(%rax,%rax,1)
movups (%r14,%rsi,4),%xmm0
movups 0x10(%r14,%rsi,4),%xmm1
movups (%r15,%rsi,4),%xmm2
subps  %xmm0,%xmm2
movups 0x10(%r15,%rsi,4),%xmm0
subps  %xmm1,%xmm0
movups 0x20(%r15,%rsi,4),%xmm1
movups 0x30(%r15,%rsi,4),%xmm3
movups %xmm2,(%r15,%rsi,4)
movups %xmm0,0x10(%r15,%rsi,4)
movups 0x20(%r14,%rsi,4),%xmm0
subps  %xmm0,%xmm1
movups 0x30(%r14,%rsi,4),%xmm0
subps  %xmm0,%xmm3
movups %xmm1,0x20(%r15,%rsi,4)
movups %xmm3,0x30(%r15,%rsi,4)
add    $0x10,%rsi
add    $0xfffffffffffffffe,%rdi
jne    1380 <wafil3+0x230>
test   $0x1,%bpl
je     1408 <wafil3+0x2b8>
movups (%r14,%rsi,4),%xmm0
movups 0x10(%r14,%rsi,4),%xmm1
movups (%r15,%rsi,4),%xmm2
subps  %xmm0,%xmm2
movups 0x10(%r15,%rsi,4),%xmm0
subps  %xmm1,%xmm0
movups %xmm2,(%r15,%rsi,4)
movups %xmm0,0x10(%r15,%rsi,4)
cmp    %rcx,%rdx
je     14ba <wafil3+0x36a>
or     $0x1,%rcx
mov    %eax,%esi
sub    %ecx,%esi
mov    %rcx,%rdx
not    %rdx
add    %rax,%rdx
and    $0x3,%rsi
je     144c <wafil3+0x2fc>
nopl   0x0(%rax,%rax,1)
movss  (%r15,%rcx,4),%xmm0
subss  (%r14,%rcx,4),%xmm0
movss  %xmm0,(%r15,%rcx,4)
add    $0x1,%rcx
add    $0xffffffffffffffff,%rsi
jne    1430 <wafil3+0x2e0>
cmp    $0x3,%rdx
jb     14ba <wafil3+0x36a>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
movss  (%r15,%rcx,4),%xmm0
subss  (%r14,%rcx,4),%xmm0
movss  0x4(%r15,%rcx,4),%xmm1
movss  %xmm0,(%r15,%rcx,4)
subss  0x4(%r14,%rcx,4),%xmm1
movss  %xmm1,0x4(%r15,%rcx,4)
movss  0x8(%r15,%rcx,4),%xmm0
subss  0x8(%r14,%rcx,4),%xmm0
movss  %xmm0,0x8(%r15,%rcx,4)
movss  0xc(%r15,%rcx,4),%xmm0
subss  0xc(%r14,%rcx,4),%xmm0
movss  %xmm0,0xc(%r15,%rcx,4)
add    $0x4,%rcx
cmp    %rcx,%rax
jne    1460 <wafil3+0x310>
mov    %r13,%rdi
mov    $0x1,%esi
call   1070 <free_dvector@plt>
add    $0x38,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret
mov    $0x1,%esi
test   $0x1,%bpl
jne    13e1 <wafil3+0x291>
jmp    1408 <wafil3+0x2b8>

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

LLVM IR (-O2):
define void @wafil3(float* nocapture %0, i32 %1, i32 %2, float* nocapture %3, float* %4) {
%6 = alloca double
%7 = bitcast double* %6 to i8*
%8 = sub i32 %1, %2
%9 = icmp sgt i32 %8, %2
%10 = shl i32 %2, 1
br i1 %9, label %13, label %11
11: ; preds = %5
%12 = sext i32 %8 to i64
br label %66
13: ; preds = %5
%14 = or i32 %10, 1
%15 = icmp slt i32 %14, 1
%16 = add i32 %10, 2
%17 = sext i32 %2 to i64
%18 = sext i32 %8 to i64
%19 = zext i32 %16 to i64
%20 = icmp eq i32 %10, 0
%21 = add nsw i64 %19, -2
br label %35
22: ; preds = %42, %39
%23 = phi float [ 0.000000e+00, %39 ], [ %62, %42 ]
%24 = phi i64 [ 1, %39 ], [ %63, %42 ]
%25 = getelementptr inbounds float, float* %3, i64 %24
%26 = load float, float* %25
%27 = trunc i64 %24 to i32
%28 = add i32 %41, %27
%29 = sext i32 %28 to i64
%30 = getelementptr inbounds float, float* %0, i64 %29
%31 = load float, float* %30
%32 = tail call float @llvm.fmuladd.f32(float %26, float %31, float %23)
store float %32, float* %38
br label %33
33: ; preds = %22, %35
%34 = icmp eq i64 %37, %18
br i1 %34, label %66, label %35
35: ; preds = %13, %33
%36 = phi i64 [ %17, %13 ], [ %37, %33 ]
%37 = add nsw i64 %36, 1
%38 = getelementptr inbounds float, float* %4, i64 %37
store float 0.000000e+00, float* %38
br i1 %15, label %33, label %39
39: ; preds = %35
%40 = trunc i64 %36 to i32
%41 = sub i32 %40, %2
br i1 %20, label %22, label %42
42: ; preds = %39, %42
%43 = phi float [ %62, %42 ], [ 0.000000e+00, %39 ]
%44 = phi i64 [ %63, %42 ], [ 1, %39 ]
%45 = phi i64 [ %64, %42 ], [ 0, %39 ]
%46 = getelementptr inbounds float, float* %3, i64 %44
%47 = load float, float* %46
%48 = trunc i64 %44 to i32
%49 = add i32 %41, %48
%50 = sext i32 %49 to i64
%51 = getelementptr inbounds float, float* %0, i64 %50
%52 = load float, float* %51
%53 = tail call float @llvm.fmuladd.f32(float %47, float %52, float %43)
store float %53, float* %38
%54 = add nuw nsw i64 %44, 1
%55 = getelementptr inbounds float, float* %3, i64 %54
%56 = load float, float* %55
%57 = trunc i64 %54 to i32
%58 = add i32 %41, %57
%59 = sext i32 %58 to i64
%60 = getelementptr inbounds float, float* %0, i64 %59
%61 = load float, float* %60
%62 = tail call float @llvm.fmuladd.f32(float %56, float %61, float %53)
store float %62, float* %38
%63 = add nuw nsw i64 %44, 2
%64 = add i64 %45, 2
%65 = icmp eq i64 %64, %21
br i1 %65, label %22, label %42
66: ; preds = %33, %11
%67 = phi i64 [ %12, %11 ], [ %18, %33 ]
%68 = sub nsw i32 %8, %2
%69 = icmp slt i32 %68, %10
%70 = select i1 %69, i32 %68, i32 %10
%71 = icmp slt i32 %70, %2
%72 = select i1 %71, i32 %70, i32 %2
%73 = icmp slt i32 %72, 256
%74 = select i1 %73, i32 %72, i32 256
%75 = tail call double* @dvector(i32 1, i32 %74)
%76 = add i32 %70, %2
%77 = sub i32 %1, %76
%78 = sext i32 %77 to i64
%79 = getelementptr inbounds float, float* %4, i64 %78
%80 = call i32 @memcof(float* %79, i32 %70, i32 %74, double* %6, double* %75)
%81 = getelementptr inbounds float, float* %4, i64 %67
%82 = call i32 @bacdic(float* %79, i32 %70, double* %75, i32 %74, float* %81, i32 %2)
%83 = getelementptr inbounds float, float* %4, i64 1
%84 = call i32 @invrse(float* %83, i32 %1)
%85 = call i32 @memcof(float* %79, i32 %70, i32 %74, double* %6, double* %75)
%86 = call i32 @bacdic(float* %79, i32 %70, double* %75, i32 %74, float* %81, i32 %2)
%87 = call i32 @invrse(float* %83, i32 %1)
%88 = icmp slt i32 %1, 1
br i1 %88, label %224, label %89
89: ; preds = %66
%90 = add nuw i32 %1, 1
%91 = zext i32 %90 to i64
%92 = add nsw i64 %91, -1
%93 = icmp ult i64 %92, 8
br i1 %93, label %176, label %94
94: ; preds = %89
%95 = getelementptr float, float* %0, i64 1
%96 = getelementptr float, float* %0, i64 %91
%97 = getelementptr float, float* %4, i64 1
%98 = getelementptr float, float* %4, i64 %91
%99 = icmp ult float* %95, %98
%100 = icmp ult float* %97, %96
%101 = and i1 %99, %100
br i1 %101, label %176, label %102
102: ; preds = %94
%103 = and i64 %92, -8
%104 = or i64 %103, 1
%105 = add nsw i64 %103, -8
%106 = lshr exact i64 %105, 3
%107 = add nuw nsw i64 %106, 1
%108 = and i64 %107, 1
%109 = icmp eq i64 %105, 0
br i1 %109, label %154, label %110
110: ; preds = %102
%111 = and i64 %107, 4611686018427387902
br label %112
112: ; preds = %112, %110
%113 = phi i64 [ 0, %110 ], [ %149, %112 ]
%114 = phi i64 [ 0, %110 ], [ %150, %112 ]
%115 = or i64 %113, 1
%116 = getelementptr inbounds float, float* %4, i64 %115
%117 = bitcast float* %116 to <4 x float>*
%118 = load <4 x float>, <4 x float>* %117
%119 = getelementptr inbounds float, float* %116, i64 4
%120 = bitcast float* %119 to <4 x float>*
%121 = load <4 x float>, <4 x float>* %120
%122 = getelementptr inbounds float, float* %0, i64 %115
%123 = bitcast float* %122 to <4 x float>*
%124 = load <4 x float>, <4 x float>* %123
%125 = getelementptr inbounds float, float* %122, i64 4
%126 = bitcast float* %125 to <4 x float>*
%127 = load <4 x float>, <4 x float>* %126
%128 = fsub <4 x float> %124, %118
%129 = fsub <4 x float> %127, %121
%130 = bitcast float* %122 to <4 x float>*
store <4 x float> %128, <4 x float>* %130
%131 = bitcast float* %125 to <4 x float>*
store <4 x float> %129, <4 x float>* %131
%132 = or i64 %113, 9
%133 = getelementptr inbounds float, float* %4, i64 %132
%134 = bitcast float* %133 to <4 x float>*
%135 = load <4 x float>, <4 x float>* %134
%136 = getelementptr inbounds float, float* %133, i64 4
%137 = bitcast float* %136 to <4 x float>*
%138 = load <4 x float>, <4 x float>* %137
%139 = getelementptr inbounds float, float* %0, i64 %132
%140 = bitcast float* %139 to <4 x float>*
%141 = load <4 x float>, <4 x float>* %140
%142 = getelementptr inbounds float, float* %139, i64 4
%143 = bitcast float* %142 to <4 x float>*
%144 = load <4 x float>, <4 x float>* %143
%145 = fsub <4 x float> %141, %135
%146 = fsub <4 x float> %144, %138
%147 = bitcast float* %139 to <4 x float>*
store <4 x float> %145, <4 x float>* %147
%148 = bitcast float* %142 to <4 x float>*
store <4 x float> %146, <4 x float>* %148
%149 = add nuw i64 %113, 16
%150 = add i64 %114, 2
%151 = icmp eq i64 %150, %111
br i1 %151, label %152, label %112
152: ; preds = %112
%153 = or i64 %149, 1
br label %154
154: ; preds = %152, %102
%155 = phi i64 [ 1, %102 ], [ %153, %152 ]
%156 = icmp eq i64 %108, 0
br i1 %156, label %174, label %157
157: ; preds = %154
%158 = getelementptr inbounds float, float* %4, i64 %155
%159 = bitcast float* %158 to <4 x float>*
%160 = load <4 x float>, <4 x float>* %159
%161 = getelementptr inbounds float, float* %158, i64 4
%162 = bitcast float* %161 to <4 x float>*
%163 = load <4 x float>, <4 x float>* %162
%164 = getelementptr inbounds float, float* %0, i64 %155
%165 = bitcast float* %164 to <4 x float>*
%166 = load <4 x float>, <4 x float>* %165
%167 = getelementptr inbounds float, float* %164, i64 4
%168 = bitcast float* %167 to <4 x float>*
%169 = load <4 x float>, <4 x float>* %168
%170 = fsub <4 x float> %166, %160
%171 = fsub <4 x float> %169, %163
%172 = bitcast float* %164 to <4 x float>*
store <4 x float> %170, <4 x float>* %172
%173 = bitcast float* %167 to <4 x float>*
store <4 x float> %171, <4 x float>* %173
br label %174
174: ; preds = %154, %157
%175 = icmp eq i64 %92, %103
br i1 %175, label %224, label %176
176: ; preds = %94, %89, %174
%177 = phi i64 [ 1, %94 ], [ 1, %89 ], [ %104, %174 ]
%178 = sub nsw i64 %91, %177
%179 = xor i64 %177, -1
%180 = add nsw i64 %179, %91
%181 = and i64 %178, 3
%182 = icmp eq i64 %181, 0
br i1 %182, label %194, label %183
183: ; preds = %176, %183
%184 = phi i64 [ %191, %183 ], [ %177, %176 ]
%185 = phi i64 [ %192, %183 ], [ 0, %176 ]
%186 = getelementptr inbounds float, float* %4, i64 %184
%187 = load float, float* %186
%188 = getelementptr inbounds float, float* %0, i64 %184
%189 = load float, float* %188
%190 = fsub float %189, %187
store float %190, float* %188
%191 = add nuw nsw i64 %184, 1
%192 = add i64 %185, 1
%193 = icmp eq i64 %192, %181
br i1 %193, label %194, label %183
194: ; preds = %183, %176
%195 = phi i64 [ %177, %176 ], [ %191, %183 ]
%196 = icmp ult i64 %180, 3
br i1 %196, label %224, label %197
197: ; preds = %194, %197
%198 = phi i64 [ %222, %197 ], [ %195, %194 ]
%199 = getelementptr inbounds float, float* %4, i64 %198
%200 = load float, float* %199
%201 = getelementptr inbounds float, float* %0, i64 %198
%202 = load float, float* %201
%203 = fsub float %202, %200
store float %203, float* %201
%204 = add nuw nsw i64 %198, 1
%205 = getelementptr inbounds float, float* %4, i64 %204
%206 = load float, float* %205
%207 = getelementptr inbounds float, float* %0, i64 %204
%208 = load float, float* %207
%209 = fsub float %208, %206
store float %209, float* %207
%210 = add nuw nsw i64 %198, 2
%211 = getelementptr inbounds float, float* %4, i64 %210
%212 = load float, float* %211
%213 = getelementptr inbounds float, float* %0, i64 %210
%214 = load float, float* %213
%215 = fsub float %214, %212
store float %215, float* %213
%216 = add nuw nsw i64 %198, 3
%217 = getelementptr inbounds float, float* %4, i64 %216
%218 = load float, float* %217
%219 = getelementptr inbounds float, float* %0, i64 %216
%220 = load float, float* %219
%221 = fsub float %220, %218
store float %221, float* %219
%222 = add nuw nsw i64 %198, 4
%223 = icmp eq i64 %222, %91
br i1 %223, label %224, label %197
224: ; preds = %194, %197, %174, %66
%225 = call i32 @free_dvector(double* %75, i32 1)
ret void
}