#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 */
double fabs(double);
void ep1R(double *x, double *t, double *u, double *v, int n) {
  int j;
  for (j = 0; j < n; j++) {
    if (fabs(u[j]) > fabs(v[j])) {
      t[j] = (fabs(u[j]) + v[j]) / 2;
      if (u[j] > 0)
        x[j] = t[j];
      else
        x[j] = -t[j];
    } else if (fabs(u[j]) <= v[j]) {
      t[j] = v[j];
      x[j] = u[j];
    } else {
      t[j] = x[j] = 0;
    }
  }
}


x86 assembly (clang, -O0):
<ep1R>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %rcx,-0x20(%rbp)
mov    %r8d,-0x24(%rbp)
movl   $0x0,-0x28(%rbp)
mov    -0x28(%rbp),%eax
cmp    -0x24(%rbp),%eax
jge    12b2 <ep1R+0x1b2>
mov    -0x18(%rbp),%rax
movslq -0x28(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movaps 0xec1(%rip),%xmm1
pand   %xmm1,%xmm0
mov    -0x20(%rbp),%rax
movslq -0x28(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
movaps 0xea9(%rip),%xmm2
pand   %xmm2,%xmm1
ucomisd %xmm1,%xmm0
jbe    1212 <ep1R+0x112>
mov    -0x18(%rbp),%rax
movslq -0x28(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movaps 0xe87(%rip),%xmm1
pand   %xmm1,%xmm0
mov    -0x20(%rbp),%rax
movslq -0x28(%rbp),%rcx
addsd  (%rax,%rcx,8),%xmm0
movsd  0xe7e(%rip),%xmm1
divsd  %xmm1,%xmm0
mov    -0x10(%rbp),%rax
movslq -0x28(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x18(%rbp),%rax
movslq -0x28(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jbe    11dc <ep1R+0xdc>
mov    -0x10(%rbp),%rax
movslq -0x28(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x8(%rbp),%rax
movslq -0x28(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
jmp    120d <ep1R+0x10d>
mov    -0x10(%rbp),%rax
movslq -0x28(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mov    -0x8(%rbp),%rax
movslq -0x28(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
jmp    129f <ep1R+0x19f>
mov    -0x18(%rbp),%rax
movslq -0x28(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
movaps 0xdda(%rip),%xmm0
pand   %xmm0,%xmm1
mov    -0x20(%rbp),%rax
movslq -0x28(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
ucomisd %xmm1,%xmm0
jb     127a <ep1R+0x17a>
mov    -0x20(%rbp),%rax
movslq -0x28(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x28(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x18(%rbp),%rax
movslq -0x28(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x8(%rbp),%rax
movslq -0x28(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
jmp    129a <ep1R+0x19a>
mov    -0x8(%rbp),%rax
movslq -0x28(%rbp),%rcx
xorps  %xmm0,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x10(%rbp),%rax
movslq -0x28(%rbp),%rcx
xorps  %xmm0,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
jmp    129f <ep1R+0x19f>
jmp    12a4 <ep1R+0x1a4>
mov    -0x28(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x28(%rbp)
jmp    111f <ep1R+0x1f>
pop    %rbp
ret

x86 assembly (clang, -O2):
<ep1R>:
test   %r8d,%r8d
jle    11b9 <ep1R+0xb9>
mov    %r8d,%r8d
xor    %eax,%eax
movapd 0xeea(%rip),%xmm0
movsd  0xf02(%rip),%xmm1
xorpd  %xmm2,%xmm2
movapd 0xee6(%rip),%xmm3
jmp    114e <ep1R+0x4e>
nopl   0x0(%rax)
movq   $0x0,(%rdi,%rax,8)
xorpd  %xmm4,%xmm4
mov    %rsi,%r9
movsd  %xmm4,(%r9,%rax,8)
add    $0x1,%rax
cmp    %rax,%r8
je     11b9 <ep1R+0xb9>
movsd  (%rdx,%rax,8),%xmm5
andpd  %xmm0,%xmm5
movsd  (%rcx,%rax,8),%xmm4
movapd %xmm4,%xmm6
andpd  %xmm0,%xmm6
ucomisd %xmm6,%xmm5
jbe    1190 <ep1R+0x90>
addsd  %xmm5,%xmm4
mulsd  %xmm1,%xmm4
movsd  %xmm4,(%rsi,%rax,8)
movsd  (%rdx,%rax,8),%xmm5
ucomisd %xmm2,%xmm5
jbe    11b0 <ep1R+0xb0>
mov    %rdi,%r9
jmp    113f <ep1R+0x3f>
nopw   0x0(%rax,%rax,1)
ucomisd %xmm5,%xmm4
jb     1130 <ep1R+0x30>
movsd  %xmm4,(%rsi,%rax,8)
movsd  (%rdx,%rax,8),%xmm4
mov    %rdi,%r9
jmp    113f <ep1R+0x3f>
cs nopw 0x0(%rax,%rax,1)
nop
xorpd  %xmm3,%xmm4
mov    %rdi,%r9
jmp    113f <ep1R+0x3f>
ret

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

LLVM IR (-O2):
define void @ep1R(double* nocapture writeonly %0, double* nocapture writeonly %1, double* nocapture %2, double* nocapture %3, i32 %4) {
%6 = icmp sgt i32 %4, 0
br i1 %6, label %7, label %39
7: ; preds = %5
%8 = zext i32 %4 to i64
br label %9
9: ; preds = %7, %33
%10 = phi i64 [ 0, %7 ], [ %37, %33 ]
%11 = getelementptr inbounds double, double* %2, i64 %10
%12 = load double, double* %11
%13 = tail call double @llvm.fabs.f64(double %12)
%14 = getelementptr inbounds double, double* %3, i64 %10
%15 = load double, double* %14
%16 = tail call double @llvm.fabs.f64(double %15)
%17 = fcmp ogt double %13, %16
br i1 %17, label %18, label %26
18: ; preds = %9
%19 = fadd double %13, %15
%20 = fmul double %19, 5.000000e-01
%21 = getelementptr inbounds double, double* %1, i64 %10
store double %20, double* %21
%22 = load double, double* %11
%23 = fcmp ogt double %22, 0.000000e+00
br i1 %23, label %33, label %24
24: ; preds = %18
%25 = fneg double %20
br label %33
26: ; preds = %9
%27 = fcmp ugt double %13, %15
br i1 %27, label %31, label %28
28: ; preds = %26
%29 = getelementptr inbounds double, double* %1, i64 %10
store double %15, double* %29
%30 = load double, double* %11
br label %33
31: ; preds = %26
%32 = getelementptr inbounds double, double* %0, i64 %10
store double 0.000000e+00, double* %32
br label %33
33: ; preds = %18, %24, %31, %28
%34 = phi double* [ %0, %24 ], [ %1, %31 ], [ %0, %28 ], [ %0, %18 ]
%35 = phi double [ %25, %24 ], [ 0.000000e+00, %31 ], [ %30, %28 ], [ %20, %18 ]
%36 = getelementptr inbounds double, double* %34, i64 %10
store double %35, double* %36
%37 = add nuw nsw i64 %10, 1
%38 = icmp eq i64 %37, %8
br i1 %38, label %39, label %9
39: ; preds = %33, %5
ret void
}