#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 */
typedef scalar_t__ uint;

/* Variables and functions */
double fabs(double);
double pow(double, double);
void DC(double amplitude, double *output, uint *index, uint Ni, double expon,
        double srate, double dt) {
  uint i;
  if (amplitude == 0.)
    (*index) += Ni;
  else {
    if ((expon == -1) || (expon == 0))
      amplitude = fabs(amplitude);
    else if (expon != 1)
      amplitude = pow(amplitude, expon);
    for (i = 0; i < Ni; i++)
      output[(*index)++] = amplitude;
  }
  return;
}


x86 assembly (clang, -O0):
<DC>:
push   %rbp
mov    %rsp,%rbp
sub    $0x40,%rsp
movsd  %xmm0,-0x8(%rbp)
mov    %rdi,-0x10(%rbp)
mov    %rsi,-0x18(%rbp)
mov    %rdx,-0x20(%rbp)
movsd  %xmm1,-0x28(%rbp)
movsd  %xmm2,-0x30(%rbp)
movsd  %xmm3,-0x38(%rbp)
movsd  -0x8(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    1163 <DC+0x53>
jp     1163 <DC+0x53>
mov    -0x20(%rbp),%rcx
mov    -0x18(%rbp),%rax
add    (%rax),%rcx
mov    %rcx,(%rax)
jmp    123d <DC+0x12d>
movsd  -0x28(%rbp),%xmm0
movsd  0xe90(%rip),%xmm1
ucomisd %xmm1,%xmm0
jne    1185 <DC+0x75>
jp     1185 <DC+0x75>
jmp    119d <DC+0x8d>
movsd  -0x28(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    11b7 <DC+0xa7>
jp     11b7 <DC+0xa7>
movsd  -0x8(%rbp),%xmm0
movaps 0xe67(%rip),%xmm1
pand   %xmm1,%xmm0
movsd  %xmm0,-0x8(%rbp)
jmp    11f2 <DC+0xe2>
movsd  -0x28(%rbp),%xmm0
movsd  0xe44(%rip),%xmm1
ucomisd %xmm1,%xmm0
jne    11d9 <DC+0xc9>
jp     11d9 <DC+0xc9>
jmp    11ed <DC+0xdd>
movsd  -0x8(%rbp),%xmm0
movsd  -0x28(%rbp),%xmm1
call   1030 <pow@plt>
movsd  %xmm0,-0x8(%rbp)
jmp    11f2 <DC+0xe2>
movq   $0x0,-0x40(%rbp)
mov    -0x40(%rbp),%rax
cmp    -0x20(%rbp),%rax
jge    1238 <DC+0x128>
movsd  -0x8(%rbp),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x18(%rbp),%rdx
mov    (%rdx),%rcx
mov    %rcx,%rsi
add    $0x1,%rsi
mov    %rsi,(%rdx)
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x40(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x40(%rbp)
jmp    11fa <DC+0xea>
jmp    123d <DC+0x12d>
add    $0x40,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<DC>:
push   %r15
push   %r14
push   %rbx
mov    %rdx,%rbx
mov    %rsi,%r14
xorpd  %xmm2,%xmm2
ucomisd %xmm2,%xmm0
jne    1132 <DC+0x22>
jp     1132 <DC+0x22>
add    (%r14),%rbx
mov    %rbx,%rax
jmp    1240 <DC+0x130>
mov    %rdi,%r15
ucomisd 0xec3(%rip),%xmm1
jne    1141 <DC+0x31>
jnp    1160 <DC+0x50>
xorpd  %xmm2,%xmm2
ucomisd %xmm2,%xmm1
jne    114d <DC+0x3d>
jnp    1160 <DC+0x50>
ucomisd 0xeb3(%rip),%xmm1
jne    1159 <DC+0x49>
jnp    1168 <DC+0x58>
call   1030 <pow@plt>
jmp    1168 <DC+0x58>
andpd  0xea8(%rip),%xmm0
test   %rbx,%rbx
jle    1243 <DC+0x133>
mov    (%r14),%rax
cmp    $0x4,%rbx
jae    1182 <DC+0x72>
xor    %r8d,%r8d
jmp    122c <DC+0x11c>
mov    %rbx,%r8
and    $0xfffffffffffffffc,%r8
movapd %xmm0,%xmm1
unpcklpd %xmm0,%xmm1
lea    -0x4(%r8),%rcx
mov    %rcx,%rsi
shr    $0x2,%rsi
add    $0x1,%rsi
mov    %esi,%edx
and    $0x3,%edx
cmp    $0xc,%rcx
jae    11af <DC+0x9f>
xor    %edi,%edi
jmp    11f9 <DC+0xe9>
and    $0xfffffffffffffffc,%rsi
lea    (%r15,%rax,8),%rcx
add    $0x70,%rcx
xor    %edi,%edi
nopl   (%rax)
movupd %xmm1,-0x70(%rcx,%rdi,8)
movupd %xmm1,-0x60(%rcx,%rdi,8)
movupd %xmm1,-0x50(%rcx,%rdi,8)
movupd %xmm1,-0x40(%rcx,%rdi,8)
movupd %xmm1,-0x30(%rcx,%rdi,8)
movupd %xmm1,-0x20(%rcx,%rdi,8)
movupd %xmm1,-0x10(%rcx,%rdi,8)
movupd %xmm1,(%rcx,%rdi,8)
add    $0x10,%rdi
add    $0xfffffffffffffffc,%rsi
jne    11c0 <DC+0xb0>
test   %rdx,%rdx
je     1224 <DC+0x114>
add    %rax,%rdi
lea    (%r15,%rdi,8),%rcx
add    $0x10,%rcx
shl    $0x5,%rdx
xor    %esi,%esi
nop
movupd %xmm1,-0x10(%rcx,%rsi,1)
movupd %xmm1,(%rcx,%rsi,1)
add    $0x20,%rsi
cmp    %rsi,%rdx
jne    1210 <DC+0x100>
add    %r8,%rax
cmp    %rbx,%r8
je     1240 <DC+0x130>
sub    %r8,%rbx
nop
movsd  %xmm0,(%r15,%rax,8)
add    $0x1,%rax
add    $0xffffffffffffffff,%rbx
jne    1230 <DC+0x120>
mov    %rax,(%r14)
pop    %rbx
pop    %r14
pop    %r15
ret

LLVM IR (-O0):
define void @DC(double %0, double* %1, i64* %2, i64 %3, double %4, double %5, double %6) {
%8 = alloca double
%9 = alloca double*
%10 = alloca i64*
%11 = alloca i64
%12 = alloca double
%13 = alloca double
%14 = alloca double
%15 = alloca i64
store double %0, double* %8
store double* %1, double** %9
store i64* %2, i64** %10
store i64 %3, i64* %11
store double %4, double* %12
store double %5, double* %13
store double %6, double* %14
%16 = load double, double* %8
%17 = fcmp oeq double %16, 0.000000e+00
br i1 %17, label %18, label %23
18: ; preds = %7
%19 = load i64, i64* %11
%20 = load i64*, i64** %10
%21 = load i64, i64* %20
%22 = add nsw i64 %21, %19
store i64 %22, i64* %20
br label %56
23: ; preds = %7
%24 = load double, double* %12
%25 = fcmp oeq double %24, -1.000000e+00
br i1 %25, label %29, label %26
26: ; preds = %23
%27 = load double, double* %12
%28 = fcmp oeq double %27, 0.000000e+00
br i1 %28, label %29, label %32
29: ; preds = %26, %23
%30 = load double, double* %8
%31 = call double @llvm.fabs.f64(double %30)
store double %31, double* %8
br label %40
32: ; preds = %26
%33 = load double, double* %12
%34 = fcmp une double %33, 1.000000e+00
br i1 %34, label %35, label %39
35: ; preds = %32
%36 = load double, double* %8
%37 = load double, double* %12
%38 = call double @pow(double %36, double %37)
store double %38, double* %8
br label %39
39: ; preds = %35, %32
br label %40
40: ; preds = %39, %29
store i64 0, i64* %15
br label %41
41: ; preds = %52, %40
%42 = load i64, i64* %15
%43 = load i64, i64* %11
%44 = icmp slt i64 %42, %43
br i1 %44, label %45, label %55
45: ; preds = %41
%46 = load double, double* %8
%47 = load double*, double** %9
%48 = load i64*, i64** %10
%49 = load i64, i64* %48
%50 = add nsw i64 %49, 1
store i64 %50, i64* %48
%51 = getelementptr inbounds double, double* %47, i64 %49
store double %46, double* %51
br label %52
52: ; preds = %45
%53 = load i64, i64* %15
%54 = add nsw i64 %53, 1
store i64 %54, i64* %15
br label %41
55: ; preds = %41
br label %56
56: ; preds = %55, %18
ret void
}

LLVM IR (-O2):
define void @DC(double %0, double* nocapture writeonly %1, i64* nocapture %2, i64 %3, double %4, double %5, double %6) {
%8 = fcmp oeq double %0, 0.000000e+00
br i1 %8, label %9, label %12
9: ; preds = %7
%10 = load i64, i64* %2
%11 = add nsw i64 %10, %3
br label %97
12: ; preds = %7
%13 = fcmp oeq double %4, -1.000000e+00
%14 = fcmp oeq double %4, 0.000000e+00
%15 = or i1 %13, %14
br i1 %15, label %16, label %18
16: ; preds = %12
%17 = tail call double @llvm.fabs.f64(double %0)
br label %22
18: ; preds = %12
%19 = fcmp une double %4, 1.000000e+00
br i1 %19, label %20, label %22
20: ; preds = %18
%21 = tail call double @pow(double %0, double %4)
br label %22
22: ; preds = %18, %20, %16
%23 = phi double [ %17, %16 ], [ %21, %20 ], [ %0, %18 ]
%24 = icmp sgt i64 %3, 0
br i1 %24, label %25, label %99
25: ; preds = %22
%26 = load i64, i64* %2
%27 = icmp ult i64 %3, 4
br i1 %27, label %87, label %28
28: ; preds = %25
%29 = and i64 %3, -4
%30 = add i64 %26, %29
%31 = insertelement <2 x double> poison, double %23, i64 0
%32 = shufflevector <2 x double> %31, <2 x double> poison, <2 x i32> zeroinitializer
%33 = insertelement <2 x double> poison, double %23, i64 0
%34 = shufflevector <2 x double> %33, <2 x double> poison, <2 x i32> zeroinitializer
%35 = add i64 %29, -4
%36 = lshr exact i64 %35, 2
%37 = add nuw nsw i64 %36, 1
%38 = and i64 %37, 3
%39 = icmp ult i64 %35, 12
br i1 %39, label %71, label %40
40: ; preds = %28
%41 = and i64 %37, 9223372036854775804
br label %42
42: ; preds = %42, %40
%43 = phi i64 [ 0, %40 ], [ %68, %42 ]
%44 = phi i64 [ 0, %40 ], [ %69, %42 ]
%45 = add i64 %26, %43
%46 = getelementptr inbounds double, double* %1, i64 %45
%47 = bitcast double* %46 to <2 x double>*
store <2 x double> %32, <2 x double>* %47
%48 = getelementptr inbounds double, double* %46, i64 2
%49 = bitcast double* %48 to <2 x double>*
store <2 x double> %34, <2 x double>* %49
%50 = or i64 %43, 4
%51 = add i64 %26, %50
%52 = getelementptr inbounds double, double* %1, i64 %51
%53 = bitcast double* %52 to <2 x double>*
store <2 x double> %32, <2 x double>* %53
%54 = getelementptr inbounds double, double* %52, i64 2
%55 = bitcast double* %54 to <2 x double>*
store <2 x double> %34, <2 x double>* %55
%56 = or i64 %43, 8
%57 = add i64 %26, %56
%58 = getelementptr inbounds double, double* %1, i64 %57
%59 = bitcast double* %58 to <2 x double>*
store <2 x double> %32, <2 x double>* %59
%60 = getelementptr inbounds double, double* %58, i64 2
%61 = bitcast double* %60 to <2 x double>*
store <2 x double> %34, <2 x double>* %61
%62 = or i64 %43, 12
%63 = add i64 %26, %62
%64 = getelementptr inbounds double, double* %1, i64 %63
%65 = bitcast double* %64 to <2 x double>*
store <2 x double> %32, <2 x double>* %65
%66 = getelementptr inbounds double, double* %64, i64 2
%67 = bitcast double* %66 to <2 x double>*
store <2 x double> %34, <2 x double>* %67
%68 = add nuw i64 %43, 16
%69 = add nuw i64 %44, 4
%70 = icmp eq i64 %69, %41
br i1 %70, label %71, label %42
71: ; preds = %42, %28
%72 = phi i64 [ 0, %28 ], [ %68, %42 ]
%73 = icmp eq i64 %38, 0
br i1 %73, label %85, label %74
74: ; preds = %71, %74
%75 = phi i64 [ %82, %74 ], [ %72, %71 ]
%76 = phi i64 [ %83, %74 ], [ 0, %71 ]
%77 = add i64 %26, %75
%78 = getelementptr inbounds double, double* %1, i64 %77
%79 = bitcast double* %78 to <2 x double>*
store <2 x double> %32, <2 x double>* %79
%80 = getelementptr inbounds double, double* %78, i64 2
%81 = bitcast double* %80 to <2 x double>*
store <2 x double> %34, <2 x double>* %81
%82 = add nuw i64 %75, 4
%83 = add i64 %76, 1
%84 = icmp eq i64 %83, %38
br i1 %84, label %85, label %74
85: ; preds = %74, %71
%86 = icmp eq i64 %29, %3
br i1 %86, label %97, label %87
87: ; preds = %25, %85
%88 = phi i64 [ %26, %25 ], [ %30, %85 ]
%89 = phi i64 [ 0, %25 ], [ %29, %85 ]
br label %90
90: ; preds = %87, %90
%91 = phi i64 [ %93, %90 ], [ %88, %87 ]
%92 = phi i64 [ %95, %90 ], [ %89, %87 ]
%93 = add nsw i64 %91, 1
%94 = getelementptr inbounds double, double* %1, i64 %91
store double %23, double* %94
%95 = add nuw nsw i64 %92, 1
%96 = icmp eq i64 %95, %3
br i1 %96, label %97, label %90
97: ; preds = %90, %85, %9
%98 = phi i64 [ %11, %9 ], [ %30, %85 ], [ %93, %90 ]
store i64 %98, i64* %2
br label %99
99: ; preds = %97, %22
ret void
}