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


x86 assembly (clang, -O0):
<dedisperse_lagged>:
push   %rbp
mov    %rsp,%rbp
sub    $0x50,%rsp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %edx,-0x14(%rbp)
mov    %ecx,-0x18(%rbp)
cmpl   $0x2,-0x14(%rbp)
setge  %al
and    $0x1,%al
movzbl %al,%edi
call   1040 <assert@plt>
cmpl   $0x1,-0x18(%rbp)
setge  %al
and    $0x1,%al
movzbl %al,%edi
call   1040 <assert@plt>
mov    -0x8(%rbp),%rax
mov    0x8(%rax),%rax
mov    -0x8(%rbp),%rcx
mov    (%rcx),%rcx
sub    %rcx,%rax
sar    $0x2,%rax
mov    -0x14(%rbp),%ecx
add    -0x18(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
cmp    %rcx,%rax
setge  %al
and    $0x1,%al
movzbl %al,%edi
call   1040 <assert@plt>
mov    -0x10(%rbp),%rax
mov    0x8(%rax),%rax
mov    -0x10(%rbp),%rcx
mov    (%rcx),%rcx
sub    %rcx,%rax
sar    $0x2,%rax
mov    -0x14(%rbp),%ecx
add    -0x18(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
cmp    %rcx,%rax
setge  %al
and    $0x1,%al
movzbl %al,%edi
call   1040 <assert@plt>
mov    -0x14(%rbp),%edi
call   1030 <get_npass@plt>
mov    %eax,-0x1c(%rbp)
mov    -0x14(%rbp),%eax
mov    -0x1c(%rbp),%ecx
mov    $0x1,%edx
shl    %cl,%edx
mov    %edx,%ecx
cmp    %ecx,%eax
sete   %al
and    $0x1,%al
movzbl %al,%edi
call   1040 <assert@plt>
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    12c5 <dedisperse_lagged+0x185>
movl   $0x0,-0x38(%rbp)
mov    -0x38(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    1291 <dedisperse_lagged+0x151>
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),%eax
mov    %eax,-0x48(%rbp)
mov    -0x18(%rbp),%ecx
mov    -0x38(%rbp),%eax
cltd
idivl  -0x20(%rbp)
add    %eax,%ecx
mov    -0x38(%rbp),%eax
cltd
idivl  -0x20(%rbp)
mov    -0x48(%rbp),%edx
mov    %eax,%r8d
call   1060 <dedisperse_kernel_lagged@plt>
mov    -0x20(%rbp),%eax
add    -0x38(%rbp),%eax
mov    %eax,-0x38(%rbp)
jmp    1235 <dedisperse_lagged+0xf5>
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    -0x20(%rbp),%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
mov    %eax,-0x20(%rbp)
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    1222 <dedisperse_lagged+0xe2>
movl   $0x0,-0x44(%rbp)
mov    -0x44(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    1311 <dedisperse_lagged+0x1d1>
mov    -0x30(%rbp),%rax
movslq -0x44(%rbp),%rcx
mov    (%rax,%rcx,8),%rdi
mov    -0x28(%rbp),%rax
movslq -0x44(%rbp),%rcx
mov    (%rax,%rcx,8),%rsi
mov    -0x18(%rbp),%eax
add    -0x44(%rbp),%eax
cltq
shl    $0x2,%rax
mov    %eax,%edx
call   1050 <memcpy@plt>
mov    -0x44(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x44(%rbp)
jmp    12cc <dedisperse_lagged+0x18c>
add    $0x50,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<dedisperse_lagged>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x38,%rsp
mov    %ecx,%r14d
mov    %edx,%ebp
mov    %rsi,%r15
mov    %rdi,%r12
xor    %edi,%edi
cmp    $0x2,%edx
setge  %dil
call   1040 <assert@plt>
xor    %edi,%edi
test   %r14d,%r14d
setg   %dil
call   1040 <assert@plt>
mov    0x8(%r12),%rax
sub    (%r12),%rax
sar    $0x2,%rax
mov    %r14,0x8(%rsp)
lea    (%r14,%rbp,1),%ecx
add    $0xffffffff,%ecx
movslq %ecx,%r14
xor    %edi,%edi
cmp    %r14,%rax
setge  %dil
call   1040 <assert@plt>
mov    0x8(%r15),%rax
sub    (%r15),%rax
sar    $0x2,%rax
xor    %edi,%edi
cmp    %r14,%rax
setge  %dil
call   1040 <assert@plt>
mov    %ebp,%edi
call   1030 <get_npass@plt>
mov    %eax,%ebx
mov    $0x1,%eax
mov    %ebx,%ecx
shl    %cl,%eax
xor    %edi,%edi
cmp    %ebp,%eax
sete   %dil
call   1040 <assert@plt>
mov    %ebx,0x4(%rsp)
test   %ebx,%ebx
jle    1292 <dedisperse_lagged+0x152>
movslq %ebp,%r14
xor    %edx,%edx
mov    %ebp,%eax
mov    %ebp,(%rsp)
mov    %rbp,0x20(%rsp)
mov    0x8(%rsp),%r13
jmp    122f <dedisperse_lagged+0xef>
nopl   0x0(%rax)
mov    (%rsp),%edx
mov    %edx,%eax
shr    $0x1f,%eax
add    %edx,%eax
sar    %eax
mov    0x28(%rsp),%rdx
add    $0x1,%edx
mov    %eax,(%rsp)
mov    0x10(%rsp),%rbx
mov    %rbx,%r15
cmp    0x4(%rsp),%edx
mov    0x18(%rsp),%r12
mov    0x20(%rsp),%rbp
je     1295 <dedisperse_lagged+0x155>
mov    %rdx,0x28(%rsp)
mov    %r12,0x10(%rsp)
mov    %r15,0x18(%rsp)
test   %ebp,%ebp
jle    1200 <dedisperse_lagged+0xc0>
movslq (%rsp),%rbp
lea    0x0(,%rbp,8),%rax
mov    %rax,0x30(%rsp)
mov    0x18(%rsp),%rbx
mov    0x10(%rsp),%r15
xor    %r12d,%r12d
mov    %r12d,%eax
cltd
idiv   %ebp
lea    (%rax,%r13,1),%ecx
mov    %r15,%rdi
mov    %rbx,%rsi
mov    %ebp,%edx
mov    %eax,%r8d
call   1060 <dedisperse_kernel_lagged@plt>
add    %rbp,%r12
mov    0x30(%rsp),%rax
add    %rax,%r15
add    %rax,%rbx
cmp    %r14,%r12
jl     1260 <dedisperse_lagged+0x120>
jmp    1200 <dedisperse_lagged+0xc0>
mov    %r15,%rbx
test   %ebp,%ebp
mov    0x8(%rsp),%r14
jle    12ca <dedisperse_lagged+0x18a>
mov    %ebp,%r15d
shl    $0x2,%r14d
shl    $0x2,%r15
xor    %ebp,%ebp
nopl   0x0(%rax,%rax,1)
mov    (%rbx,%rbp,2),%rdi
mov    (%r12,%rbp,2),%rsi
lea    (%r14,%rbp,1),%edx
call   1050 <memcpy@plt>
add    $0x4,%rbp
cmp    %rbp,%r15
jne    12b0 <dedisperse_lagged+0x170>
add    $0x38,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

LLVM IR (-O0):
define void @dedisperse_lagged(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**
%16 = alloca i32
store float** %0, float*** %5
store float** %1, float*** %6
store i32 %2, i32* %7
store i32 %3, i32* %8
%17 = load i32, i32* %7
%18 = icmp sge i32 %17, 2
%19 = zext i1 %18 to i32
%20 = call i32 @assert(i32 %19)
%21 = load i32, i32* %8
%22 = icmp sge i32 %21, 1
%23 = zext i1 %22 to i32
%24 = call i32 @assert(i32 %23)
%25 = load float**, float*** %5
%26 = getelementptr inbounds float*, float** %25, i64 1
%27 = load float*, float** %26
%28 = load float**, float*** %5
%29 = getelementptr inbounds float*, float** %28, i64 0
%30 = load float*, float** %29
%31 = ptrtoint float* %27 to i64
%32 = ptrtoint float* %30 to i64
%33 = sub i64 %31, %32
%34 = sdiv exact i64 %33, 4
%35 = load i32, i32* %7
%36 = load i32, i32* %8
%37 = add nsw i32 %35, %36
%38 = sub nsw i32 %37, 1
%39 = sext i32 %38 to i64
%40 = icmp sge i64 %34, %39
%41 = zext i1 %40 to i32
%42 = call i32 @assert(i32 %41)
%43 = load float**, float*** %6
%44 = getelementptr inbounds float*, float** %43, i64 1
%45 = load float*, float** %44
%46 = load float**, float*** %6
%47 = getelementptr inbounds float*, float** %46, i64 0
%48 = load float*, float** %47
%49 = ptrtoint float* %45 to i64
%50 = ptrtoint float* %48 to i64
%51 = sub i64 %49, %50
%52 = sdiv exact i64 %51, 4
%53 = load i32, i32* %7
%54 = load i32, i32* %8
%55 = add nsw i32 %53, %54
%56 = sub nsw i32 %55, 1
%57 = sext i32 %56 to i64
%58 = icmp sge i64 %52, %57
%59 = zext i1 %58 to i32
%60 = call i32 @assert(i32 %59)
%61 = load i32, i32* %7
%62 = call i32 @get_npass(i32 %61)
store i32 %62, i32* %9
%63 = load i32, i32* %7
%64 = load i32, i32* %9
%65 = shl i32 1, %64
%66 = icmp eq i32 %63, %65
%67 = zext i1 %66 to i32
%68 = call i32 @assert(i32 %67)
%69 = load i32, i32* %7
store i32 %69, i32* %10
%70 = load float**, float*** %5
store float** %70, float*** %11
%71 = load float**, float*** %6
store float** %71, float*** %12
store i32 0, i32* %13
br label %72
72: ; preds = %110, %4
%73 = load i32, i32* %13
%74 = load i32, i32* %9
%75 = icmp slt i32 %73, %74
br i1 %75, label %76, label %113
76: ; preds = %72
store i32 0, i32* %14
br label %77
77: ; preds = %100, %76
%78 = load i32, i32* %14
%79 = load i32, i32* %7
%80 = icmp slt i32 %78, %79
br i1 %80, label %81, label %104
81: ; preds = %77
%82 = load float**, float*** %11
%83 = load i32, i32* %14
%84 = sext i32 %83 to i64
%85 = getelementptr inbounds float*, float** %82, i64 %84
%86 = load float**, float*** %12
%87 = load i32, i32* %14
%88 = sext i32 %87 to i64
%89 = getelementptr inbounds float*, float** %86, i64 %88
%90 = load i32, i32* %10
%91 = load i32, i32* %8
%92 = load i32, i32* %14
%93 = load i32, i32* %10
%94 = sdiv i32 %92, %93
%95 = add nsw i32 %91, %94
%96 = load i32, i32* %14
%97 = load i32, i32* %10
%98 = sdiv i32 %96, %97
%99 = call i32 @dedisperse_kernel_lagged(float** %85, float** %89, i32 %90, i32 %95, i32 %98)
br label %100
100: ; preds = %81
%101 = load i32, i32* %10
%102 = load i32, i32* %14
%103 = add nsw i32 %102, %101
store i32 %103, i32* %14
br label %77
104: ; preds = %77
%105 = load float**, float*** %11
store float** %105, float*** %15
%106 = load float**, float*** %12
store float** %106, float*** %11
%107 = load float**, float*** %15
store float** %107, float*** %12
%108 = load i32, i32* %10
%109 = sdiv i32 %108, 2
store i32 %109, i32* %10
br label %110
110: ; preds = %104
%111 = load i32, i32* %13
%112 = add nsw i32 %111, 1
store i32 %112, i32* %13
br label %72
113: ; preds = %72
store i32 0, i32* %16
br label %114
114: ; preds = %136, %113
%115 = load i32, i32* %16
%116 = load i32, i32* %7
%117 = icmp slt i32 %115, %116
br i1 %117, label %118, label %139
118: ; preds = %114
%119 = load float**, float*** %12
%120 = load i32, i32* %16
%121 = sext i32 %120 to i64
%122 = getelementptr inbounds float*, float** %119, i64 %121
%123 = load float*, float** %122
%124 = load float**, float*** %11
%125 = load i32, i32* %16
%126 = sext i32 %125 to i64
%127 = getelementptr inbounds float*, float** %124, i64 %126
%128 = load float*, float** %127
%129 = load i32, i32* %8
%130 = load i32, i32* %16
%131 = add nsw i32 %129, %130
%132 = sext i32 %131 to i64
%133 = mul i64 %132, 4
%134 = trunc i64 %133 to i32
%135 = call i32 @memcpy(float* %123, float* %128, i32 %134)
br label %136
136: ; preds = %118
%137 = load i32, i32* %16
%138 = add nsw i32 %137, 1
store i32 %138, i32* %16
br label %114
139: ; preds = %114
ret void
}

LLVM IR (-O2):
define void @dedisperse_lagged(float** %0, float** %1, i32 %2, i32 %3) {
%5 = icmp sgt i32 %2, 1
%6 = zext i1 %5 to i32
%7 = tail call i32 @assert(i32 %6)
%8 = icmp sgt i32 %3, 0
%9 = zext i1 %8 to i32
%10 = tail call i32 @assert(i32 %9)
%11 = getelementptr inbounds float*, float** %0, i64 1
%12 = load float*, float** %11
%13 = load float*, float** %0
%14 = ptrtoint float* %12 to i64
%15 = ptrtoint float* %13 to i64
%16 = sub i64 %14, %15
%17 = ashr exact i64 %16, 2
%18 = add i32 %2, -1
%19 = add i32 %18, %3
%20 = sext i32 %19 to i64
%21 = icmp sge i64 %17, %20
%22 = zext i1 %21 to i32
%23 = tail call i32 @assert(i32 %22)
%24 = getelementptr inbounds float*, float** %1, i64 1
%25 = load float*, float** %24
%26 = load float*, float** %1
%27 = ptrtoint float* %25 to i64
%28 = ptrtoint float* %26 to i64
%29 = sub i64 %27, %28
%30 = ashr exact i64 %29, 2
%31 = icmp sge i64 %30, %20
%32 = zext i1 %31 to i32
%33 = tail call i32 @assert(i32 %32)
%34 = tail call i32 @get_npass(i32 %2)
%35 = shl nuw i32 1, %34
%36 = icmp eq i32 %35, %2
%37 = zext i1 %36 to i32
%38 = tail call i32 @assert(i32 %37)
%39 = icmp sgt i32 %34, 0
br i1 %39, label %40, label %50
40: ; preds = %4
%41 = icmp sgt i32 %2, 0
%42 = sext i32 %2 to i64
br label %43
43: ; preds = %40, %56
%44 = phi i32 [ %2, %40 ], [ %57, %56 ]
%45 = phi i32 [ 0, %40 ], [ %58, %56 ]
%46 = phi float** [ %1, %40 ], [ %47, %56 ]
%47 = phi float** [ %0, %40 ], [ %46, %56 ]
br i1 %41, label %48, label %56
48: ; preds = %43
%49 = sext i32 %44 to i64
br label %60
50: ; preds = %56, %4
%51 = phi float** [ %0, %4 ], [ %46, %56 ]
%52 = phi float** [ %1, %4 ], [ %47, %56 ]
%53 = icmp sgt i32 %2, 0
br i1 %53, label %54, label %70
54: ; preds = %50
%55 = zext i32 %2 to i64
br label %71
56: ; preds = %60, %43
%57 = sdiv i32 %44, 2
%58 = add nuw nsw i32 %45, 1
%59 = icmp eq i32 %58, %34
br i1 %59, label %50, label %43
60: ; preds = %48, %60
%61 = phi i64 [ 0, %48 ], [ %68, %60 ]
%62 = getelementptr inbounds float*, float** %47, i64 %61
%63 = getelementptr inbounds float*, float** %46, i64 %61
%64 = trunc i64 %61 to i32
%65 = sdiv i32 %64, %44
%66 = add nsw i32 %65, %3
%67 = tail call i32 @dedisperse_kernel_lagged(float** %62, float** %63, i32 %44, i32 %66, i32 %65)
%68 = add i64 %61, %49
%69 = icmp slt i64 %68, %42
br i1 %69, label %60, label %56
70: ; preds = %71, %50
ret void
71: ; preds = %54, %71
%72 = phi i64 [ 0, %54 ], [ %81, %71 ]
%73 = getelementptr inbounds float*, float** %52, i64 %72
%74 = load float*, float** %73
%75 = getelementptr inbounds float*, float** %51, i64 %72
%76 = load float*, float** %75
%77 = trunc i64 %72 to i32
%78 = add i32 %77, %3
%79 = shl i32 %78, 2
%80 = tail call i32 @memcpy(float* %74, float* %76, i32 %79)
%81 = add nuw nsw i64 %72, 1
%82 = icmp eq i64 %81, %55
br i1 %82, label %70, label %71
}