#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 ti_zlema_start(double const *);
int ti_zlema(int size, double const *const *inputs, double const *options,
             double *const *outputs) {
  const double *input = inputs[0];
  const int period = (int)options[0];
  const int lag = (period - 1) / 2;
  double *output = outputs[0];
  if (period < 1)
    return 1;
  if (size <= ti_zlema_start(options))
    return 0;
  const double per = 2 / ((double)period + 1);
  double val = input[lag - 1];
  *output++ = val;
  int i;
  for (i = lag; i < size; ++i) {
    double c = input[i];
    double l = input[i - lag];
    val = ((c + (c - l)) - val) * per + val;
    *output++ = val;
  }
  assert(output - outputs[0] == size - ti_zlema_start(options));
  return 0;
}


x86 assembly (clang, -O0):
<ti_zlema>:
push   %rbp
mov    %rsp,%rbp
sub    $0x80,%rsp
mov    %edi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %rcx,-0x20(%rbp)
mov    -0x10(%rbp),%rax
mov    (%rax),%rax
mov    %rax,-0x28(%rbp)
mov    -0x18(%rbp),%rax
cvttsd2si (%rax),%eax
mov    %eax,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
sub    $0x1,%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
mov    %eax,-0x30(%rbp)
mov    -0x20(%rbp),%rax
mov    (%rax),%rax
mov    %rax,-0x38(%rbp)
cmpl   $0x1,-0x2c(%rbp)
jge    1182 <ti_zlema+0x62>
movl   $0x1,-0x4(%rbp)
jmp    12d5 <ti_zlema+0x1b5>
mov    -0x8(%rbp),%eax
mov    %eax,-0x64(%rbp)
mov    -0x18(%rbp),%rdi
call   1040 <ti_zlema_start@plt>
mov    %eax,%ecx
mov    -0x64(%rbp),%eax
cmp    %ecx,%eax
jg     11aa <ti_zlema+0x8a>
movl   $0x0,-0x4(%rbp)
jmp    12d5 <ti_zlema+0x1b5>
cvtsi2sdl -0x2c(%rbp),%xmm1
movsd  0xe51(%rip),%xmm0
addsd  %xmm0,%xmm1
movsd  0xe3d(%rip),%xmm0
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x40(%rbp)
mov    -0x28(%rbp),%rax
mov    -0x30(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x48(%rbp)
movsd  -0x48(%rbp),%xmm0
mov    -0x38(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x38(%rbp)
movsd  %xmm0,(%rax)
mov    -0x30(%rbp),%eax
mov    %eax,-0x4c(%rbp)
mov    -0x4c(%rbp),%eax
cmp    -0x8(%rbp),%eax
jge    128b <ti_zlema+0x16b>
mov    -0x28(%rbp),%rax
movslq -0x4c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x58(%rbp)
mov    -0x28(%rbp),%rax
mov    -0x4c(%rbp),%ecx
sub    -0x30(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x60(%rbp)
movsd  -0x58(%rbp),%xmm0
movsd  -0x58(%rbp),%xmm1
subsd  -0x60(%rbp),%xmm1
addsd  %xmm1,%xmm0
subsd  -0x48(%rbp),%xmm0
movsd  -0x40(%rbp),%xmm2
movsd  -0x48(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x48(%rbp)
movsd  -0x48(%rbp),%xmm0
mov    -0x38(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x38(%rbp)
movsd  %xmm0,(%rax)
mov    -0x4c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x4c(%rbp)
jmp    1201 <ti_zlema+0xe1>
mov    -0x38(%rbp),%rax
mov    -0x20(%rbp),%rcx
mov    (%rcx),%rcx
sub    %rcx,%rax
sar    $0x3,%rax
mov    %rax,-0x70(%rbp)
mov    -0x8(%rbp),%eax
mov    %eax,-0x74(%rbp)
mov    -0x18(%rbp),%rdi
call   1040 <ti_zlema_start@plt>
mov    -0x74(%rbp),%ecx
mov    %eax,%edx
mov    -0x70(%rbp),%rax
sub    %edx,%ecx
movslq %ecx,%rcx
cmp    %rcx,%rax
sete   %al
and    $0x1,%al
movzbl %al,%edi
call   1030 <assert@plt>
movl   $0x0,-0x4(%rbp)
mov    -0x4(%rbp),%eax
add    $0x80,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<ti_zlema>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
push   %rax
cvttsd2si (%rdx),%r14d
test   %r14d,%r14d
jle    11be <ti_zlema+0x9e>
mov    %rcx,%r15
mov    %rdx,%r12
mov    %edi,%r13d
mov    (%rsi),%rbx
mov    (%rcx),%rbp
mov    %rdx,%rdi
call   1040 <ti_zlema_start@plt>
xor    %ecx,%ecx
cmp    %r13d,%eax
jge    128c <ti_zlema+0x16c>
lea    -0x1(%r14),%eax
shr    $0x1f,%eax
add    %r14d,%eax
add    $0xffffffff,%eax
sar    %eax
movslq %eax,%rsi
movsd  -0x8(%rbx,%rsi,8),%xmm0
movsd  %xmm0,0x0(%rbp)
lea    0x8(%rbp),%rcx
cmp    %r13d,%esi
jge    1264 <ti_zlema+0x144>
cvtsi2sd %r14d,%xmm2
addsd  0xe6e(%rip),%xmm2
movsd  0xe6e(%rip),%xmm1
divsd  %xmm2,%xmm1
movslq %r13d,%rax
mov    %r13d,%edx
sub    %esi,%edx
mov    %rsi,%rdi
not    %rdi
test   $0x1,%dl
jne    11c8 <ti_zlema+0xa8>
mov    %rsi,%rdx
add    %rax,%rdi
jne    11fa <ti_zlema+0xda>
jmp    1267 <ti_zlema+0x147>
mov    $0x1,%ecx
jmp    128c <ti_zlema+0x16c>
movsd  (%rbx,%rsi,8),%xmm2
movapd %xmm2,%xmm3
subsd  (%rbx),%xmm3
addsd  %xmm2,%xmm3
subsd  %xmm0,%xmm3
mulsd  %xmm1,%xmm3
addsd  %xmm3,%xmm0
movsd  %xmm0,0x8(%rbp)
lea    0x1(%rsi),%rdx
add    $0x10,%rbp
mov    %rbp,%rcx
add    %rax,%rdi
je     1267 <ti_zlema+0x147>
mov    $0x1,%edi
sub    %rsi,%rdi
lea    (%rbx,%rdi,8),%rsi
cs nopw 0x0(%rax,%rax,1)
movsd  (%rbx,%rdx,8),%xmm2
movapd %xmm2,%xmm3
subsd  -0x8(%rsi,%rdx,8),%xmm3
addsd  %xmm2,%xmm3
subsd  %xmm0,%xmm3
mulsd  %xmm1,%xmm3
addsd  %xmm0,%xmm3
movsd  %xmm3,(%rcx)
movsd  0x8(%rbx,%rdx,8),%xmm2
movapd %xmm2,%xmm0
subsd  (%rsi,%rdx,8),%xmm0
addsd  %xmm2,%xmm0
subsd  %xmm3,%xmm0
mulsd  %xmm1,%xmm0
addsd  %xmm3,%xmm0
movsd  %xmm0,0x8(%rcx)
add    $0x2,%rdx
add    $0x10,%rcx
cmp    %rdx,%rax
jne    1210 <ti_zlema+0xf0>
mov    %rcx,%rbp
sub    (%r15),%rbp
sar    $0x3,%rbp
mov    %r12,%rdi
call   1040 <ti_zlema_start@plt>
sub    %eax,%r13d
movslq %r13d,%rax
xor    %edi,%edi
cmp    %rax,%rbp
sete   %dil
call   1030 <assert@plt>
xor    %ecx,%ecx
mov    %ecx,%eax
add    $0x8,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

LLVM IR (-O0):
define i32 @ti_zlema(i32 %0, double** %1, double* %2, double** %3) {
%5 = alloca i32
%6 = alloca i32
%7 = alloca double**
%8 = alloca double*
%9 = alloca double**
%10 = alloca double*
%11 = alloca i32
%12 = alloca i32
%13 = alloca double*
%14 = alloca double
%15 = alloca double
%16 = alloca i32
%17 = alloca double
%18 = alloca double
store i32 %0, i32* %6
store double** %1, double*** %7
store double* %2, double** %8
store double** %3, double*** %9
%19 = load double**, double*** %7
%20 = getelementptr inbounds double*, double** %19, i64 0
%21 = load double*, double** %20
store double* %21, double** %10
%22 = load double*, double** %8
%23 = getelementptr inbounds double, double* %22, i64 0
%24 = load double, double* %23
%25 = fptosi double %24 to i32
store i32 %25, i32* %11
%26 = load i32, i32* %11
%27 = sub nsw i32 %26, 1
%28 = sdiv i32 %27, 2
store i32 %28, i32* %12
%29 = load double**, double*** %9
%30 = getelementptr inbounds double*, double** %29, i64 0
%31 = load double*, double** %30
store double* %31, double** %13
%32 = load i32, i32* %11
%33 = icmp slt i32 %32, 1
br i1 %33, label %34, label %35
34: ; preds = %4
store i32 1, i32* %5
br label %106
35: ; preds = %4
%36 = load i32, i32* %6
%37 = load double*, double** %8
%38 = call i32 @ti_zlema_start(double* %37)
%39 = icmp sle i32 %36, %38
br i1 %39, label %40, label %41
40: ; preds = %35
store i32 0, i32* %5
br label %106
41: ; preds = %35
%42 = load i32, i32* %11
%43 = sitofp i32 %42 to double
%44 = fadd double %43, 1.000000e+00
%45 = fdiv double 2.000000e+00, %44
store double %45, double* %14
%46 = load double*, double** %10
%47 = load i32, i32* %12
%48 = sub nsw i32 %47, 1
%49 = sext i32 %48 to i64
%50 = getelementptr inbounds double, double* %46, i64 %49
%51 = load double, double* %50
store double %51, double* %15
%52 = load double, double* %15
%53 = load double*, double** %13
%54 = getelementptr inbounds double, double* %53, i32 1
store double* %54, double** %13
store double %52, double* %53
%55 = load i32, i32* %12
store i32 %55, i32* %16
br label %56
56: ; preds = %86, %41
%57 = load i32, i32* %16
%58 = load i32, i32* %6
%59 = icmp slt i32 %57, %58
br i1 %59, label %60, label %89
60: ; preds = %56
%61 = load double*, double** %10
%62 = load i32, i32* %16
%63 = sext i32 %62 to i64
%64 = getelementptr inbounds double, double* %61, i64 %63
%65 = load double, double* %64
store double %65, double* %17
%66 = load double*, double** %10
%67 = load i32, i32* %16
%68 = load i32, i32* %12
%69 = sub nsw i32 %67, %68
%70 = sext i32 %69 to i64
%71 = getelementptr inbounds double, double* %66, i64 %70
%72 = load double, double* %71
store double %72, double* %18
%73 = load double, double* %17
%74 = load double, double* %17
%75 = load double, double* %18
%76 = fsub double %74, %75
%77 = fadd double %73, %76
%78 = load double, double* %15
%79 = fsub double %77, %78
%80 = load double, double* %14
%81 = load double, double* %15
%82 = call double @llvm.fmuladd.f64(double %79, double %80, double %81)
store double %82, double* %15
%83 = load double, double* %15
%84 = load double*, double** %13
%85 = getelementptr inbounds double, double* %84, i32 1
store double* %85, double** %13
store double %83, double* %84
br label %86
86: ; preds = %60
%87 = load i32, i32* %16
%88 = add nsw i32 %87, 1
store i32 %88, i32* %16
br label %56
89: ; preds = %56
%90 = load double*, double** %13
%91 = load double**, double*** %9
%92 = getelementptr inbounds double*, double** %91, i64 0
%93 = load double*, double** %92
%94 = ptrtoint double* %90 to i64
%95 = ptrtoint double* %93 to i64
%96 = sub i64 %94, %95
%97 = sdiv exact i64 %96, 8
%98 = load i32, i32* %6
%99 = load double*, double** %8
%100 = call i32 @ti_zlema_start(double* %99)
%101 = sub nsw i32 %98, %100
%102 = sext i32 %101 to i64
%103 = icmp eq i64 %97, %102
%104 = zext i1 %103 to i32
%105 = call i32 @assert(i32 %104)
store i32 0, i32* %5
br label %106
106: ; preds = %89, %40, %34
%107 = load i32, i32* %5
ret i32 %107
}

LLVM IR (-O2):
define i32 @ti_zlema(i32 %0, double** nocapture %1, double* %2, double** nocapture %3) {
%5 = load double*, double** %1
%6 = load double, double* %2
%7 = fptosi double %6 to i32
%8 = add nsw i32 %7, -1
%9 = sdiv i32 %8, 2
%10 = load double*, double** %3
%11 = icmp slt i32 %7, 1
br i1 %11, label %89, label %12
12: ; preds = %4
%13 = tail call i32 @ti_zlema_start(double* %2)
%14 = icmp slt i32 %13, %0
br i1 %14, label %15, label %89
15: ; preds = %12
%16 = sitofp i32 %7 to double
%17 = fadd double %16, 1.000000e+00
%18 = fdiv double 2.000000e+00, %17
%19 = add nsw i32 %9, -1
%20 = sext i32 %19 to i64
%21 = getelementptr inbounds double, double* %5, i64 %20
%22 = load double, double* %21
store double %22, double* %10
%23 = getelementptr inbounds double, double* %10, i64 1
%24 = icmp slt i32 %9, %0
br i1 %24, label %25, label %76
25: ; preds = %15
%26 = sext i32 %9 to i64
%27 = sext i32 %0 to i64
%28 = sub nsw i64 %27, %26
%29 = xor i64 %26, -1
%30 = and i64 %28, 1
%31 = icmp eq i64 %30, 0
br i1 %31, label %42, label %32
32: ; preds = %25
%33 = getelementptr inbounds double, double* %5, i64 %26
%34 = load double, double* %33
%35 = load double, double* %5
%36 = fsub double %34, %35
%37 = fadd double %34, %36
%38 = fsub double %37, %22
%39 = tail call double @llvm.fmuladd.f64(double %38, double %18, double %22)
store double %39, double* %23
%40 = add nsw i64 %26, 1
%41 = getelementptr inbounds double, double* %10, i64 2
br label %42
42: ; preds = %32, %25
%43 = phi double* [ undef, %25 ], [ %41, %32 ]
%44 = phi i64 [ %26, %25 ], [ %40, %32 ]
%45 = phi double* [ %23, %25 ], [ %41, %32 ]
%46 = phi double [ %22, %25 ], [ %39, %32 ]
%47 = sub nsw i64 0, %27
%48 = icmp eq i64 %29, %47
br i1 %48, label %76, label %49
49: ; preds = %42, %49
%50 = phi i64 [ %73, %49 ], [ %44, %42 ]
%51 = phi double* [ %74, %49 ], [ %45, %42 ]
%52 = phi double [ %72, %49 ], [ %46, %42 ]
%53 = getelementptr inbounds double, double* %5, i64 %50
%54 = load double, double* %53
%55 = sub nsw i64 %50, %26
%56 = getelementptr inbounds double, double* %5, i64 %55
%57 = load double, double* %56
%58 = fsub double %54, %57
%59 = fadd double %54, %58
%60 = fsub double %59, %52
%61 = tail call double @llvm.fmuladd.f64(double %60, double %18, double %52)
store double %61, double* %51
%62 = add nsw i64 %50, 1
%63 = getelementptr inbounds double, double* %51, i64 1
%64 = getelementptr inbounds double, double* %5, i64 %62
%65 = load double, double* %64
%66 = sub nsw i64 %62, %26
%67 = getelementptr inbounds double, double* %5, i64 %66
%68 = load double, double* %67
%69 = fsub double %65, %68
%70 = fadd double %65, %69
%71 = fsub double %70, %61
%72 = tail call double @llvm.fmuladd.f64(double %71, double %18, double %61)
store double %72, double* %63
%73 = add nsw i64 %50, 2
%74 = getelementptr inbounds double, double* %51, i64 2
%75 = icmp eq i64 %73, %27
br i1 %75, label %76, label %49
76: ; preds = %42, %49, %15
%77 = phi double* [ %23, %15 ], [ %43, %42 ], [ %74, %49 ]
%78 = load double*, double** %3
%79 = ptrtoint double* %77 to i64
%80 = ptrtoint double* %78 to i64
%81 = sub i64 %79, %80
%82 = ashr exact i64 %81, 3
%83 = tail call i32 @ti_zlema_start(double* %2)
%84 = sub nsw i32 %0, %83
%85 = sext i32 %84 to i64
%86 = icmp eq i64 %82, %85
%87 = zext i1 %86 to i32
%88 = tail call i32 @assert(i32 %87)
br label %89
89: ; preds = %12, %4, %76
%90 = phi i32 [ 0, %76 ], [ 1, %4 ], [ 0, %12 ]
ret i32 %90
}