#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*/ dedisperse_kernel(float **, float **, int, int);
int get_npass(int);
int /*<<< orphan*/ memcpy(float *, float *, int);
void dedisperse(float **inin, float **outout, int nchan, int ndat) {
  int npass = get_npass(nchan);
  int bs = nchan;
  float **in = inin;
  float **out = outout;
  for (int i = 0; i < npass; i++) {
    for (int j = 0; j < nchan; j += bs) {
      dedisperse_kernel(in + j, out + j, bs, ndat);
    }
    bs /= 2;
    float **tmp = in;
    in = out;
    out = tmp;
  }
  memcpy(out[0], in[0], nchan * ndat * sizeof(float));
}


x86 assembly (clang, -O0):
<dedisperse>:
push   %rbp
mov    %rsp,%rbp
sub    $0x40,%rsp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %edx,-0x14(%rbp)
mov    %ecx,-0x18(%rbp)
mov    -0x14(%rbp),%edi
call   1030 <get_npass@plt>
mov    %eax,-0x1c(%rbp)
mov    -0x14(%rbp),%eax
mov    %eax,-0x20(%rbp)
mov    -0x8(%rbp),%rax
mov    %rax,-0x28(%rbp)
mov    -0x10(%rbp),%rax
mov    %rax,-0x30(%rbp)
movl   $0x0,-0x34(%rbp)
mov    -0x34(%rbp),%eax
cmp    -0x1c(%rbp),%eax
jge    11f8 <dedisperse+0xc8>
movl   $0x0,-0x38(%rbp)
mov    -0x38(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    11c4 <dedisperse+0x94>
mov    -0x28(%rbp),%rdi
movslq -0x38(%rbp),%rax
shl    $0x3,%rax
add    %rax,%rdi
mov    -0x30(%rbp),%rsi
movslq -0x38(%rbp),%rax
shl    $0x3,%rax
add    %rax,%rsi
mov    -0x20(%rbp),%edx
mov    -0x18(%rbp),%ecx
call   1050 <dedisperse_kernel@plt>
mov    -0x20(%rbp),%eax
add    -0x38(%rbp),%eax
mov    %eax,-0x38(%rbp)
jmp    1181 <dedisperse+0x51>
mov    -0x20(%rbp),%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
mov    %eax,-0x20(%rbp)
mov    -0x28(%rbp),%rax
mov    %rax,-0x40(%rbp)
mov    -0x30(%rbp),%rax
mov    %rax,-0x28(%rbp)
mov    -0x40(%rbp),%rax
mov    %rax,-0x30(%rbp)
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    116e <dedisperse+0x3e>
mov    -0x30(%rbp),%rax
mov    (%rax),%rdi
mov    -0x28(%rbp),%rax
mov    (%rax),%rsi
mov    -0x14(%rbp),%eax
imul   -0x18(%rbp),%eax
cltq
shl    $0x2,%rax
mov    %eax,%edx
call   1040 <memcpy@plt>
add    $0x40,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<dedisperse>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x28,%rsp
mov    %ecx,(%rsp)
mov    %edx,%ebp
mov    %rsi,%rbx
mov    %rdi,0x8(%rsp)
mov    %edx,%edi
call   1030 <get_npass@plt>
mov    %eax,0x14(%rsp)
test   %eax,%eax
jle    11f6 <dedisperse+0xc6>
movslq %ebp,%r14
xor    %esi,%esi
mov    %ebp,0x4(%rsp)
mov    (%rsp),%ecx
mov    %ebp,0x10(%rsp)
jmp    119b <dedisperse+0x6b>
mov    0x4(%rsp),%esi
mov    %esi,%eax
shr    $0x1f,%eax
add    %esi,%eax
sar    %eax
mov    0x20(%rsp),%rsi
add    $0x1,%esi
mov    %eax,0x4(%rsp)
mov    0x18(%rsp),%rax
mov    %rax,%rbx
cmp    0x14(%rsp),%esi
mov    0x10(%rsp),%ebp
je     11fc <dedisperse+0xcc>
mov    %rsi,0x20(%rsp)
mov    0x8(%rsp),%rax
mov    %rax,0x18(%rsp)
mov    %rbx,0x8(%rsp)
test   %ebp,%ebp
jle    1170 <dedisperse+0x40>
movslq 0x4(%rsp),%rbx
lea    0x0(,%rbx,8),%r12
mov    0x8(%rsp),%rbp
mov    0x18(%rsp),%r15
xor    %r13d,%r13d
nopl   (%rax)
mov    %r15,%rdi
mov    %rbp,%rsi
mov    %ebx,%edx
mov    (%rsp),%ecx
call   1050 <dedisperse_kernel@plt>
mov    (%rsp),%ecx
add    %rbx,%r13
add    %r12,%r15
add    %r12,%rbp
cmp    %r14,%r13
jl     11d0 <dedisperse+0xa0>
jmp    1170 <dedisperse+0x40>
mov    %rbx,%rax
mov    (%rsp),%ecx
mov    (%rax),%rdi
mov    0x8(%rsp),%rax
mov    (%rax),%rsi
imul   %ecx,%ebp
shl    $0x2,%ebp
mov    %ebp,%edx
add    $0x28,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
jmp    1040 <memcpy@plt>

LLVM IR (-O0):
define void @dedisperse(float** %0, float** %1, i32 %2, i32 %3) {
%5 = alloca float**
%6 = alloca float**
%7 = alloca i32
%8 = alloca i32
%9 = alloca i32
%10 = alloca i32
%11 = alloca float**
%12 = alloca float**
%13 = alloca i32
%14 = alloca i32
%15 = alloca float**
store float** %0, float*** %5
store float** %1, float*** %6
store i32 %2, i32* %7
store i32 %3, i32* %8
%16 = load i32, i32* %7
%17 = call i32 @get_npass(i32 %16)
store i32 %17, i32* %9
%18 = load i32, i32* %7
store i32 %18, i32* %10
%19 = load float**, float*** %5
store float** %19, float*** %11
%20 = load float**, float*** %6
store float** %20, float*** %12
store i32 0, i32* %13
br label %21
21: ; preds = %52, %4
%22 = load i32, i32* %13
%23 = load i32, i32* %9
%24 = icmp slt i32 %22, %23
br i1 %24, label %25, label %55
25: ; preds = %21
store i32 0, i32* %14
br label %26
26: ; preds = %42, %25
%27 = load i32, i32* %14
%28 = load i32, i32* %7
%29 = icmp slt i32 %27, %28
br i1 %29, label %30, label %46
30: ; preds = %26
%31 = load float**, float*** %11
%32 = load i32, i32* %14
%33 = sext i32 %32 to i64
%34 = getelementptr inbounds float*, float** %31, i64 %33
%35 = load float**, float*** %12
%36 = load i32, i32* %14
%37 = sext i32 %36 to i64
%38 = getelementptr inbounds float*, float** %35, i64 %37
%39 = load i32, i32* %10
%40 = load i32, i32* %8
%41 = call i32 @dedisperse_kernel(float** %34, float** %38, i32 %39, i32 %40)
br label %42
42: ; preds = %30
%43 = load i32, i32* %10
%44 = load i32, i32* %14
%45 = add nsw i32 %44, %43
store i32 %45, i32* %14
br label %26
46: ; preds = %26
%47 = load i32, i32* %10
%48 = sdiv i32 %47, 2
store i32 %48, i32* %10
%49 = load float**, float*** %11
store float** %49, float*** %15
%50 = load float**, float*** %12
store float** %50, float*** %11
%51 = load float**, float*** %15
store float** %51, float*** %12
br label %52
52: ; preds = %46
%53 = load i32, i32* %13
%54 = add nsw i32 %53, 1
store i32 %54, i32* %13
br label %21
55: ; preds = %21
%56 = load float**, float*** %12
%57 = getelementptr inbounds float*, float** %56, i64 0
%58 = load float*, float** %57
%59 = load float**, float*** %11
%60 = getelementptr inbounds float*, float** %59, i64 0
%61 = load float*, float** %60
%62 = load i32, i32* %7
%63 = load i32, i32* %8
%64 = mul nsw i32 %62, %63
%65 = sext i32 %64 to i64
%66 = mul i64 %65, 4
%67 = trunc i64 %66 to i32
%68 = call i32 @memcpy(float* %58, float* %61, i32 %67)
ret void
}

LLVM IR (-O2):
define void @dedisperse(float** %0, float** %1, i32 %2, i32 %3) {
%5 = tail call i32 @get_npass(i32 %2)
%6 = icmp sgt i32 %5, 0
br i1 %6, label %7, label %17
7: ; preds = %4
%8 = icmp sgt i32 %2, 0
%9 = sext i32 %2 to i64
br label %10
10: ; preds = %7, %25
%11 = phi i32 [ %2, %7 ], [ %26, %25 ]
%12 = phi i32 [ 0, %7 ], [ %27, %25 ]
%13 = phi float** [ %1, %7 ], [ %14, %25 ]
%14 = phi float** [ %0, %7 ], [ %13, %25 ]
br i1 %8, label %15, label %25
15: ; preds = %10
%16 = sext i32 %11 to i64
br label %29
17: ; preds = %25, %4
%18 = phi float** [ %0, %4 ], [ %13, %25 ]
%19 = phi float** [ %1, %4 ], [ %14, %25 ]
%20 = load float*, float** %19
%21 = load float*, float** %18
%22 = shl i32 %2, 2
%23 = mul i32 %22, %3
%24 = tail call i32 @memcpy(float* %20, float* %21, i32 %23)
ret void
25: ; preds = %29, %10
%26 = sdiv i32 %11, 2
%27 = add nuw nsw i32 %12, 1
%28 = icmp eq i32 %27, %5
br i1 %28, label %17, label %10
29: ; preds = %15, %29
%30 = phi i64 [ 0, %15 ], [ %34, %29 ]
%31 = getelementptr inbounds float*, float** %14, i64 %30
%32 = getelementptr inbounds float*, float** %13, i64 %30
%33 = tail call i32 @dedisperse_kernel(float** %31, float** %32, i32 %11, i32 %3)
%34 = add i64 %30, %16
%35 = icmp slt i64 %34, %9
br i1 %35, label %29, label %25
}