#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 */
float bessi0(float);
scalar_t__ pow(float, int);
float sqrt(scalar_t__);
void wdwksr(float *window, int isize, int n1, int n2, float theta) {
  if (isize != 0) {
    float bestht = 1.0 / bessi0(theta);
    float point = 1.0 / (float)isize;
    int i, i1 = (((n1) > (1)) ? (n1) : (1));
    int i2 = (((n2) < (isize)) ? (n2) : (isize));
    for (i = i1; i <= i2; i++)
      window[i - 1] =
          bestht * bessi0(theta * (sqrt(1.0 - pow((float)i * point, 2))));
  }
}


x86 assembly (clang, -O0):
<wdwksr>:
push   %rbp
mov    %rsp,%rbp
sub    $0x40,%rsp
mov    %rdi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
mov    %edx,-0x10(%rbp)
mov    %ecx,-0x14(%rbp)
movss  %xmm0,-0x18(%rbp)
cmpl   $0x0,-0xc(%rbp)
je     1277 <wdwksr+0x147>
movss  -0x18(%rbp),%xmm0
call   1030 <bessi0@plt>
cvtss2sd %xmm0,%xmm1
movsd  0xe96(%rip),%xmm0
divsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,-0x1c(%rbp)
cvtsi2ssl -0xc(%rbp),%xmm0
cvtss2sd %xmm0,%xmm1
movsd  0xe78(%rip),%xmm0
divsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,-0x20(%rbp)
cmpl   $0x1,-0x10(%rbp)
jle    11aa <wdwksr+0x7a>
mov    -0x10(%rbp),%eax
mov    %eax,-0x30(%rbp)
jmp    11b7 <wdwksr+0x87>
mov    $0x1,%eax
mov    %eax,-0x30(%rbp)
jmp    11b7 <wdwksr+0x87>
mov    -0x30(%rbp),%eax
mov    %eax,-0x28(%rbp)
mov    -0x14(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    11d4 <wdwksr+0xa4>
mov    -0x14(%rbp),%eax
mov    %eax,-0x34(%rbp)
jmp    11da <wdwksr+0xaa>
mov    -0xc(%rbp),%eax
mov    %eax,-0x34(%rbp)
mov    -0x34(%rbp),%eax
mov    %eax,-0x2c(%rbp)
mov    -0x28(%rbp),%eax
mov    %eax,-0x24(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0x2c(%rbp),%eax
jg     1272 <wdwksr+0x142>
movss  -0x1c(%rbp),%xmm0
movss  %xmm0,-0x38(%rbp)
movss  -0x18(%rbp),%xmm0
movss  %xmm0,-0x3c(%rbp)
cvtsi2ssl -0x24(%rbp),%xmm0
mulss  -0x20(%rbp),%xmm0
mov    $0x2,%edi
call   1040 <pow@plt>
cvtsi2sd %rax,%xmm1
movsd  0xdd9(%rip),%xmm0
subsd  %xmm1,%xmm0
cvttsd2si %xmm0,%rdi
call   1050 <sqrt@plt>
movaps %xmm0,%xmm1
movss  -0x3c(%rbp),%xmm0
mulss  %xmm1,%xmm0
call   1030 <bessi0@plt>
movaps %xmm0,%xmm1
movss  -0x38(%rbp),%xmm0
mulss  %xmm1,%xmm0
mov    -0x8(%rbp),%rax
mov    -0x24(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    11e6 <wdwksr+0xb6>
jmp    1277 <wdwksr+0x147>
add    $0x40,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<wdwksr>:
push   %rbp
push   %r15
push   %r14
push   %rbx
sub    $0x18,%rsp
movss  %xmm0,0xc(%rsp)
test   %esi,%esi
je     1206 <wdwksr+0xd6>
mov    %ecx,%ebx
mov    %edx,%r15d
mov    %esi,%ebp
mov    %rdi,%r14
movss  0xc(%rsp),%xmm0
call   1030 <bessi0@plt>
test   %r15d,%r15d
mov    $0x1,%eax
cmovg  %r15d,%eax
cmp    %ebp,%ebx
cmovge %ebp,%ebx
cmp    %ebx,%eax
jg     1206 <wdwksr+0xd6>
movss  0xe82(%rip),%xmm1
movaps %xmm1,%xmm2
divss  %xmm0,%xmm2
movss  %xmm2,0x10(%rsp)
xorps  %xmm0,%xmm0
cvtsi2ss %ebp,%xmm0
divss  %xmm0,%xmm1
movss  %xmm1,0x14(%rsp)
mov    %eax,%ebp
add    $0x1,%ebx
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
xorps  %xmm0,%xmm0
cvtsi2ss %ebp,%xmm0
mulss  0x14(%rsp),%xmm0
mov    $0x2,%edi
call   1040 <pow@plt>
xorps  %xmm0,%xmm0
cvtsi2sd %rax,%xmm0
movsd  0xe31(%rip),%xmm1
subsd  %xmm0,%xmm1
cvttsd2si %xmm1,%rdi
call   1050 <sqrt@plt>
mulss  0xc(%rsp),%xmm0
call   1030 <bessi0@plt>
mulss  0x10(%rsp),%xmm0
movss  %xmm0,-0x4(%r14,%rbp,4)
add    $0x1,%rbp
cmp    %rbp,%rbx
jne    11b0 <wdwksr+0x80>
add    $0x18,%rsp
pop    %rbx
pop    %r14
pop    %r15
pop    %rbp
ret

LLVM IR (-O0):
define void @wdwksr(float* %0, i32 %1, i32 %2, i32 %3, float %4) {
%6 = alloca float*
%7 = alloca i32
%8 = alloca i32
%9 = alloca i32
%10 = alloca float
%11 = alloca float
%12 = alloca float
%13 = alloca i32
%14 = alloca i32
%15 = alloca i32
store float* %0, float** %6
store i32 %1, i32* %7
store i32 %2, i32* %8
store i32 %3, i32* %9
store float %4, float* %10
%16 = load i32, i32* %7
%17 = icmp ne i32 %16, 0
br i1 %17, label %18, label %74
18: ; preds = %5
%19 = load float, float* %10
%20 = call float @bessi0(float %19)
%21 = fpext float %20 to double
%22 = fdiv double 1.000000e+00, %21
%23 = fptrunc double %22 to float
store float %23, float* %11
%24 = load i32, i32* %7
%25 = sitofp i32 %24 to float
%26 = fpext float %25 to double
%27 = fdiv double 1.000000e+00, %26
%28 = fptrunc double %27 to float
store float %28, float* %12
%29 = load i32, i32* %8
%30 = icmp sgt i32 %29, 1
br i1 %30, label %31, label %33
31: ; preds = %18
%32 = load i32, i32* %8
br label %34
33: ; preds = %18
br label %34
34: ; preds = %33, %31
%35 = phi i32 [ %32, %31 ], [ 1, %33 ]
store i32 %35, i32* %14
%36 = load i32, i32* %9
%37 = load i32, i32* %7
%38 = icmp slt i32 %36, %37
br i1 %38, label %39, label %41
39: ; preds = %34
%40 = load i32, i32* %9
br label %43
41: ; preds = %34
%42 = load i32, i32* %7
br label %43
43: ; preds = %41, %39
%44 = phi i32 [ %40, %39 ], [ %42, %41 ]
store i32 %44, i32* %15
%45 = load i32, i32* %14
store i32 %45, i32* %13
br label %46
46: ; preds = %70, %43
%47 = load i32, i32* %13
%48 = load i32, i32* %15
%49 = icmp sle i32 %47, %48
br i1 %49, label %50, label %73
50: ; preds = %46
%51 = load float, float* %11
%52 = load float, float* %10
%53 = load i32, i32* %13
%54 = sitofp i32 %53 to float
%55 = load float, float* %12
%56 = fmul float %54, %55
%57 = call i64 @pow(float %56, i32 2)
%58 = sitofp i64 %57 to double
%59 = fsub double 1.000000e+00, %58
%60 = fptosi double %59 to i64
%61 = call float @sqrt(i64 %60)
%62 = fmul float %52, %61
%63 = call float @bessi0(float %62)
%64 = fmul float %51, %63
%65 = load float*, float** %6
%66 = load i32, i32* %13
%67 = sub nsw i32 %66, 1
%68 = sext i32 %67 to i64
%69 = getelementptr inbounds float, float* %65, i64 %68
store float %64, float* %69
br label %70
70: ; preds = %50
%71 = load i32, i32* %13
%72 = add nsw i32 %71, 1
store i32 %72, i32* %13
br label %46
73: ; preds = %46
br label %74
74: ; preds = %73, %5
ret void
}

LLVM IR (-O2):
define void @wdwksr(float* nocapture writeonly %0, i32 %1, i32 %2, i32 %3, float %4) {
%6 = icmp eq i32 %1, 0
br i1 %6, label %38, label %7
7: ; preds = %5
%8 = tail call float @bessi0(float %4)
%9 = fdiv float 1.000000e+00, %8
%10 = sitofp i32 %1 to float
%11 = fdiv float 1.000000e+00, %10
%12 = icmp sgt i32 %2, 1
%13 = select i1 %12, i32 %2, i32 1
%14 = icmp slt i32 %3, %1
%15 = select i1 %14, i32 %3, i32 %1
%16 = icmp sgt i32 %13, %15
br i1 %16, label %38, label %17
17: ; preds = %7
%18 = zext i32 %13 to i64
%19 = add i32 %15, 1
%20 = zext i32 %19 to i64
br label %21
21: ; preds = %17, %21
%22 = phi i64 [ %18, %17 ], [ %36, %21 ]
%23 = trunc i64 %22 to i32
%24 = sitofp i32 %23 to float
%25 = fmul float %11, %24
%26 = tail call i64 @pow(float %25, i32 2)
%27 = sitofp i64 %26 to double
%28 = fsub double 1.000000e+00, %27
%29 = fptosi double %28 to i64
%30 = tail call float @sqrt(i64 %29)
%31 = fmul float %30, %4
%32 = tail call float @bessi0(float %31)
%33 = fmul float %9, %32
%34 = add nsw i64 %22, -1
%35 = getelementptr inbounds float, float* %0, i64 %34
store float %33, float* %35
%36 = add nuw nsw i64 %22, 1
%37 = icmp eq i64 %36, %20
br i1 %37, label %38, label %21
38: ; preds = %21, %7, %5
ret void
}