#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);
scalar_t__ find_parabolic_min_x(double *, double *, double *);
int /*<<< orphan*/ fprintf(int /*<<< orphan*/, char *, double, double, double,
                           ...);
int /*<<< orphan*/ stderr;
double optimum_center(double *x, double *y) {
  double result;
  assert(x[3] > x[2]);
  assert(x[2] > x[1]);
  assert(y[3] >= y[2]);
  assert(y[1] >= y[2]);
  if ((x[3] - x[1]) > 1E-7 && (x[2] - x[1]) > 1E-7 && (x[3] - x[2]) > 1E-7 &&
      find_parabolic_min_x(x, y, &result)) {
    assert(result > x[1] ||
           !fprintf(stderr, "\n\nx1=%.12le x2=%.12le x3=%.12le result=%.12le\n",
                    x[1], x[2], x[3], result) ||
           !fprintf(stderr, "y1=%.12le y2=%.12le y3=%.12le \n", y[1], y[2],
                    y[3]) ||
           !fprintf(stderr, "dy1=%.12le dy2=%.12le dy3=%.12le \n", y[1] - y[1],
                    y[2] - y[1], y[3] - y[1]) ||
           !fprintf(stderr, "x3-x1=%.12le x2-x1=%.12le result-x1=%.12le \n\n",
                    x[3] - x[2], x[2] - x[1], result - x[1]));
    assert(result < x[3]);
    return result;
  }
  return (x[3] + x[1]) / 2;
}


x86 assembly (clang, -O0):
<optimum_center>:
push   %rbp
mov    %rsp,%rbp
sub    $0x30,%rsp
mov    %rdi,-0x10(%rbp)
mov    %rsi,-0x18(%rbp)
mov    -0x10(%rbp),%rax
movsd  0x18(%rax),%xmm0
mov    -0x10(%rbp),%rax
ucomisd 0x10(%rax),%xmm0
seta   %al
and    $0x1,%al
movzbl %al,%edi
call   1040 <assert@plt>
mov    -0x10(%rbp),%rax
movsd  0x10(%rax),%xmm0
mov    -0x10(%rbp),%rax
ucomisd 0x8(%rax),%xmm0
seta   %al
and    $0x1,%al
movzbl %al,%edi
call   1040 <assert@plt>
mov    -0x18(%rbp),%rax
movsd  0x18(%rax),%xmm0
mov    -0x18(%rbp),%rax
ucomisd 0x10(%rax),%xmm0
setae  %al
and    $0x1,%al
movzbl %al,%edi
call   1040 <assert@plt>
mov    -0x18(%rbp),%rax
movsd  0x8(%rax),%xmm0
mov    -0x18(%rbp),%rax
ucomisd 0x10(%rax),%xmm0
setae  %al
and    $0x1,%al
movzbl %al,%edi
call   1040 <assert@plt>
mov    -0x10(%rbp),%rax
movsd  0x18(%rax),%xmm0
mov    -0x10(%rbp),%rax
subsd  0x8(%rax),%xmm0
movsd  0xe2a(%rip),%xmm1
ucomisd %xmm1,%xmm0
jbe    13d1 <optimum_center+0x2a1>
mov    -0x10(%rbp),%rax
movsd  0x10(%rax),%xmm0
mov    -0x10(%rbp),%rax
subsd  0x8(%rax),%xmm0
movsd  0xe06(%rip),%xmm1
ucomisd %xmm1,%xmm0
jbe    13d1 <optimum_center+0x2a1>
mov    -0x10(%rbp),%rax
movsd  0x18(%rax),%xmm0
mov    -0x10(%rbp),%rax
subsd  0x10(%rax),%xmm0
movsd  0xde2(%rip),%xmm1
ucomisd %xmm1,%xmm0
jbe    13d1 <optimum_center+0x2a1>
mov    -0x10(%rbp),%rdi
mov    -0x18(%rbp),%rsi
lea    -0x20(%rbp),%rdx
call   1030 <find_parabolic_min_x@plt>
cmp    $0x0,%rax
je     13d1 <optimum_center+0x2a1>
movsd  -0x20(%rbp),%xmm0
mov    -0x10(%rbp),%rcx
mov    $0x1,%al
ucomisd 0x8(%rcx),%xmm0
mov    %al,-0x21(%rbp)
ja     1396 <optimum_center+0x266>
mov    0x2d7d(%rip),%rax
mov    (%rax),%edi
mov    -0x10(%rbp),%rax
movsd  0x8(%rax),%xmm0
mov    -0x10(%rbp),%rax
movsd  0x10(%rax),%xmm1
mov    -0x10(%rbp),%rax
movsd  0x18(%rax),%xmm2
movsd  -0x20(%rbp),%xmm3
lea    0xd84(%rip),%rsi
mov    $0x4,%al
call   1050 <fprintf@plt>
mov    %eax,%ecx
mov    $0x1,%al
cmp    $0x0,%ecx
mov    %al,-0x21(%rbp)
je     1396 <optimum_center+0x266>
mov    0x2d36(%rip),%rax
mov    (%rax),%edi
mov    -0x18(%rbp),%rax
movsd  0x8(%rax),%xmm0
mov    -0x18(%rbp),%rax
movsd  0x10(%rax),%xmm1
mov    -0x18(%rbp),%rax
movsd  0x18(%rax),%xmm2
lea    0xd71(%rip),%rsi
mov    $0x3,%al
call   1050 <fprintf@plt>
mov    %eax,%ecx
mov    $0x1,%al
cmp    $0x0,%ecx
mov    %al,-0x21(%rbp)
je     1396 <optimum_center+0x266>
mov    0x2cf4(%rip),%rax
mov    (%rax),%edi
mov    -0x18(%rbp),%rax
movsd  0x8(%rax),%xmm0
mov    -0x18(%rbp),%rax
subsd  0x8(%rax),%xmm0
mov    -0x18(%rbp),%rax
movsd  0x10(%rax),%xmm1
mov    -0x18(%rbp),%rax
subsd  0x8(%rax),%xmm1
mov    -0x18(%rbp),%rax
movsd  0x18(%rax),%xmm2
mov    -0x18(%rbp),%rax
subsd  0x8(%rax),%xmm2
lea    0xd34(%rip),%rsi
mov    $0x3,%al
call   1050 <fprintf@plt>
mov    %eax,%ecx
mov    $0x1,%al
cmp    $0x0,%ecx
mov    %al,-0x21(%rbp)
je     1396 <optimum_center+0x266>
mov    0x2c97(%rip),%rax
mov    (%rax),%edi
mov    -0x10(%rbp),%rax
movsd  0x18(%rax),%xmm0
mov    -0x10(%rbp),%rax
subsd  0x10(%rax),%xmm0
mov    -0x10(%rbp),%rax
movsd  0x10(%rax),%xmm1
mov    -0x10(%rbp),%rax
subsd  0x8(%rax),%xmm1
movsd  -0x20(%rbp),%xmm2
mov    -0x10(%rbp),%rax
subsd  0x8(%rax),%xmm2
lea    0xcfe(%rip),%rsi
mov    $0x3,%al
call   1050 <fprintf@plt>
cmp    $0x0,%eax
setne  %al
xor    $0xff,%al
mov    %al,-0x21(%rbp)
mov    -0x21(%rbp),%al
and    $0x1,%al
movzbl %al,%edi
call   1040 <assert@plt>
movsd  -0x20(%rbp),%xmm1
mov    -0x10(%rbp),%rax
movsd  0x18(%rax),%xmm0
ucomisd %xmm1,%xmm0
seta   %al
and    $0x1,%al
movzbl %al,%edi
call   1040 <assert@plt>
movsd  -0x20(%rbp),%xmm0
movsd  %xmm0,-0x8(%rbp)
jmp    13f4 <optimum_center+0x2c4>
mov    -0x10(%rbp),%rax
movsd  0x18(%rax),%xmm0
mov    -0x10(%rbp),%rax
addsd  0x8(%rax),%xmm0
movsd  0xc1d(%rip),%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x8(%rbp)
movsd  -0x8(%rbp),%xmm0
add    $0x30,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<optimum_center>:
push   %rbp
push   %r15
push   %r14
push   %rbx
push   %rax
mov    %rdi,%rbx
movsd  0x18(%rdi),%xmm0
xor    %edi,%edi
ucomisd 0x10(%rbx),%xmm0
mov    %rsi,%r14
seta   %dil
call   1040 <assert@plt>
movsd  0x10(%rbx),%xmm0
xor    %edi,%edi
ucomisd 0x8(%rbx),%xmm0
seta   %dil
call   1040 <assert@plt>
movsd  0x18(%r14),%xmm0
xor    %edi,%edi
ucomisd 0x10(%r14),%xmm0
setae  %dil
call   1040 <assert@plt>
movsd  0x8(%r14),%xmm0
xor    %edi,%edi
ucomisd 0x10(%r14),%xmm0
setae  %dil
call   1040 <assert@plt>
movsd  0x8(%rbx),%xmm0
movsd  0x18(%rbx),%xmm1
movapd %xmm1,%xmm2
subsd  %xmm0,%xmm2
ucomisd 0xe51(%rip),%xmm2
jbe    1300 <optimum_center+0x1d0>
movsd  0x10(%rbx),%xmm2
movapd %xmm2,%xmm3
subsd  %xmm0,%xmm3
ucomisd 0xe36(%rip),%xmm3
jbe    1300 <optimum_center+0x1d0>
movapd %xmm1,%xmm3
subsd  %xmm2,%xmm3
ucomisd 0xe20(%rip),%xmm3
jbe    1300 <optimum_center+0x1d0>
mov    %rsp,%rdx
mov    %rbx,%rdi
mov    %r14,%rsi
call   1030 <find_parabolic_min_x@plt>
test   %rax,%rax
je     12f6 <optimum_center+0x1c6>
movsd  (%rsp),%xmm3
movsd  0x8(%rbx),%xmm0
mov    $0x1,%ebp
ucomisd %xmm0,%xmm3
ja     12d3 <optimum_center+0x1a3>
mov    0x2dc3(%rip),%r15
mov    (%r15),%edi
movsd  0x10(%rbx),%xmm1
movsd  0x18(%rbx),%xmm2
lea    0xddf(%rip),%rsi
mov    $0x4,%al
call   1050 <fprintf@plt>
test   %eax,%eax
je     12d3 <optimum_center+0x1a3>
mov    (%r15),%edi
movsd  0x8(%r14),%xmm0
movsd  0x10(%r14),%xmm1
movsd  0x18(%r14),%xmm2
lea    0xde3(%rip),%rsi
mov    $0x3,%al
call   1050 <fprintf@plt>
test   %eax,%eax
je     12d3 <optimum_center+0x1a3>
mov    (%r15),%edi
movsd  0x8(%r14),%xmm0
movsd  0x10(%r14),%xmm1
subsd  %xmm0,%xmm1
movsd  0x18(%r14),%xmm2
subsd  %xmm0,%xmm2
subsd  %xmm0,%xmm0
lea    0xdd0(%rip),%rsi
mov    $0x3,%al
call   1050 <fprintf@plt>
test   %eax,%eax
je     12d3 <optimum_center+0x1a3>
mov    (%r15),%edi
movsd  0x18(%rbx),%xmm0
movsd  0x8(%rbx),%xmm3
movsd  0x10(%rbx),%xmm1
subsd  %xmm1,%xmm0
subsd  %xmm3,%xmm1
movsd  (%rsp),%xmm2
subsd  %xmm3,%xmm2
lea    0xdbe(%rip),%rsi
mov    $0x3,%al
call   1050 <fprintf@plt>
xor    %ebp,%ebp
test   %eax,%eax
sete   %bpl
mov    %ebp,%edi
call   1040 <assert@plt>
movsd  0x18(%rbx),%xmm0
xor    %edi,%edi
ucomisd (%rsp),%xmm0
seta   %dil
call   1040 <assert@plt>
movsd  (%rsp),%xmm0
jmp    130c <optimum_center+0x1dc>
movsd  0x8(%rbx),%xmm0
movsd  0x18(%rbx),%xmm1
addsd  %xmm1,%xmm0
mulsd  0xcfc(%rip),%xmm0
add    $0x8,%rsp
pop    %rbx
pop    %r14
pop    %r15
pop    %rbp
ret

LLVM IR (-O0):
@stderr = global i32 0
@.str = [47 x i8] c"\0A\0Ax1=%.12le x2=%.12le x3=%.12le result=%.12le\0A\00"
@.str.1 = [32 x i8] c"y1=%.12le y2=%.12le y3=%.12le \0A\00"
@.str.2 = [35 x i8] c"dy1=%.12le dy2=%.12le dy3=%.12le \0A\00"
@.str.3 = [46 x i8] c"x3-x1=%.12le x2-x1=%.12le result-x1=%.12le \0A\0A\00"
define double @optimum_center(double* %0, double* %1) {
%3 = alloca double
%4 = alloca double*
%5 = alloca double*
%6 = alloca double
store double* %0, double** %4
store double* %1, double** %5
%7 = load double*, double** %4
%8 = getelementptr inbounds double, double* %7, i64 3
%9 = load double, double* %8
%10 = load double*, double** %4
%11 = getelementptr inbounds double, double* %10, i64 2
%12 = load double, double* %11
%13 = fcmp ogt double %9, %12
%14 = zext i1 %13 to i32
%15 = call i32 @assert(i32 %14)
%16 = load double*, double** %4
%17 = getelementptr inbounds double, double* %16, i64 2
%18 = load double, double* %17
%19 = load double*, double** %4
%20 = getelementptr inbounds double, double* %19, i64 1
%21 = load double, double* %20
%22 = fcmp ogt double %18, %21
%23 = zext i1 %22 to i32
%24 = call i32 @assert(i32 %23)
%25 = load double*, double** %5
%26 = getelementptr inbounds double, double* %25, i64 3
%27 = load double, double* %26
%28 = load double*, double** %5
%29 = getelementptr inbounds double, double* %28, i64 2
%30 = load double, double* %29
%31 = fcmp oge double %27, %30
%32 = zext i1 %31 to i32
%33 = call i32 @assert(i32 %32)
%34 = load double*, double** %5
%35 = getelementptr inbounds double, double* %34, i64 1
%36 = load double, double* %35
%37 = load double*, double** %5
%38 = getelementptr inbounds double, double* %37, i64 2
%39 = load double, double* %38
%40 = fcmp oge double %36, %39
%41 = zext i1 %40 to i32
%42 = call i32 @assert(i32 %41)
%43 = load double*, double** %4
%44 = getelementptr inbounds double, double* %43, i64 3
%45 = load double, double* %44
%46 = load double*, double** %4
%47 = getelementptr inbounds double, double* %46, i64 1
%48 = load double, double* %47
%49 = fsub double %45, %48
%50 = fcmp ogt double %49, 0x3E7AD7F29ABCAF48
br i1 %50, label %51, label %168
51: ; preds = %2
%52 = load double*, double** %4
%53 = getelementptr inbounds double, double* %52, i64 2
%54 = load double, double* %53
%55 = load double*, double** %4
%56 = getelementptr inbounds double, double* %55, i64 1
%57 = load double, double* %56
%58 = fsub double %54, %57
%59 = fcmp ogt double %58, 0x3E7AD7F29ABCAF48
br i1 %59, label %60, label %168
60: ; preds = %51
%61 = load double*, double** %4
%62 = getelementptr inbounds double, double* %61, i64 3
%63 = load double, double* %62
%64 = load double*, double** %4
%65 = getelementptr inbounds double, double* %64, i64 2
%66 = load double, double* %65
%67 = fsub double %63, %66
%68 = fcmp ogt double %67, 0x3E7AD7F29ABCAF48
br i1 %68, label %69, label %168
69: ; preds = %60
%70 = load double*, double** %4
%71 = load double*, double** %5
%72 = call i64 @find_parabolic_min_x(double* %70, double* %71, double* %6)
%73 = icmp ne i64 %72, 0
br i1 %73, label %74, label %168
74: ; preds = %69
%75 = load double, double* %6
%76 = load double*, double** %4
%77 = getelementptr inbounds double, double* %76, i64 1
%78 = load double, double* %77
%79 = fcmp ogt double %75, %78
br i1 %79, label %156, label %80
80: ; preds = %74
%81 = load i32, i32* @stderr
%82 = load double*, double** %4
%83 = getelementptr inbounds double, double* %82, i64 1
%84 = load double, double* %83
%85 = load double*, double** %4
%86 = getelementptr inbounds double, double* %85, i64 2
%87 = load double, double* %86
%88 = load double*, double** %4
%89 = getelementptr inbounds double, double* %88, i64 3
%90 = load double, double* %89
%91 = load double, double* %6
%92 = call i32 (i32, i8*, double, double, double, ...) @fprintf(i32 %81, i8* getelementptr inbounds ([47 x i8], [47 x i8]* @.str, i64 0, i64 0), double %84, double %87, double %90, double %91)
%93 = icmp ne i32 %92, 0
br i1 %93, label %94, label %156
94: ; preds = %80
%95 = load i32, i32* @stderr
%96 = load double*, double** %5
%97 = getelementptr inbounds double, double* %96, i64 1
%98 = load double, double* %97
%99 = load double*, double** %5
%100 = getelementptr inbounds double, double* %99, i64 2
%101 = load double, double* %100
%102 = load double*, double** %5
%103 = getelementptr inbounds double, double* %102, i64 3
%104 = load double, double* %103
%105 = call i32 (i32, i8*, double, double, double, ...) @fprintf(i32 %95, i8* getelementptr inbounds ([32 x i8], [32 x i8]* @.str.1, i64 0, i64 0), double %98, double %101, double %104)
%106 = icmp ne i32 %105, 0
br i1 %106, label %107, label %156
107: ; preds = %94
%108 = load i32, i32* @stderr
%109 = load double*, double** %5
%110 = getelementptr inbounds double, double* %109, i64 1
%111 = load double, double* %110
%112 = load double*, double** %5
%113 = getelementptr inbounds double, double* %112, i64 1
%114 = load double, double* %113
%115 = fsub double %111, %114
%116 = load double*, double** %5
%117 = getelementptr inbounds double, double* %116, i64 2
%118 = load double, double* %117
%119 = load double*, double** %5
%120 = getelementptr inbounds double, double* %119, i64 1
%121 = load double, double* %120
%122 = fsub double %118, %121
%123 = load double*, double** %5
%124 = getelementptr inbounds double, double* %123, i64 3
%125 = load double, double* %124
%126 = load double*, double** %5
%127 = getelementptr inbounds double, double* %126, i64 1
%128 = load double, double* %127
%129 = fsub double %125, %128
%130 = call i32 (i32, i8*, double, double, double, ...) @fprintf(i32 %108, i8* getelementptr inbounds ([35 x i8], [35 x i8]* @.str.2, i64 0, i64 0), double %115, double %122, double %129)
%131 = icmp ne i32 %130, 0
br i1 %131, label %132, label %156
132: ; preds = %107
%133 = load i32, i32* @stderr
%134 = load double*, double** %4
%135 = getelementptr inbounds double, double* %134, i64 3
%136 = load double, double* %135
%137 = load double*, double** %4
%138 = getelementptr inbounds double, double* %137, i64 2
%139 = load double, double* %138
%140 = fsub double %136, %139
%141 = load double*, double** %4
%142 = getelementptr inbounds double, double* %141, i64 2
%143 = load double, double* %142
%144 = load double*, double** %4
%145 = getelementptr inbounds double, double* %144, i64 1
%146 = load double, double* %145
%147 = fsub double %143, %146
%148 = load double, double* %6
%149 = load double*, double** %4
%150 = getelementptr inbounds double, double* %149, i64 1
%151 = load double, double* %150
%152 = fsub double %148, %151
%153 = call i32 (i32, i8*, double, double, double, ...) @fprintf(i32 %133, i8* getelementptr inbounds ([46 x i8], [46 x i8]* @.str.3, i64 0, i64 0), double %140, double %147, double %152)
%154 = icmp ne i32 %153, 0
%155 = xor i1 %154, true
br label %156
156: ; preds = %132, %107, %94, %80, %74
%157 = phi i1 [ true, %107 ], [ true, %94 ], [ true, %80 ], [ true, %74 ], [ %155, %132 ]
%158 = zext i1 %157 to i32
%159 = call i32 @assert(i32 %158)
%160 = load double, double* %6
%161 = load double*, double** %4
%162 = getelementptr inbounds double, double* %161, i64 3
%163 = load double, double* %162
%164 = fcmp olt double %160, %163
%165 = zext i1 %164 to i32
%166 = call i32 @assert(i32 %165)
%167 = load double, double* %6
store double %167, double* %3
br label %177
168: ; preds = %69, %60, %51, %2
%169 = load double*, double** %4
%170 = getelementptr inbounds double, double* %169, i64 3
%171 = load double, double* %170
%172 = load double*, double** %4
%173 = getelementptr inbounds double, double* %172, i64 1
%174 = load double, double* %173
%175 = fadd double %171, %174
%176 = fdiv double %175, 2.000000e+00
store double %176, double* %3
br label %177
177: ; preds = %168, %156
%178 = load double, double* %3
ret double %178
}

LLVM IR (-O2):
@stderr = global i32 0
@.str = [47 x i8] c"\0A\0Ax1=%.12le x2=%.12le x3=%.12le result=%.12le\0A\00"
@.str.1 = [32 x i8] c"y1=%.12le y2=%.12le y3=%.12le \0A\00"
@.str.2 = [35 x i8] c"dy1=%.12le dy2=%.12le dy3=%.12le \0A\00"
@.str.3 = [46 x i8] c"x3-x1=%.12le x2-x1=%.12le result-x1=%.12le \0A\0A\00"
define double @optimum_center(double* %0, double* %1) {
%3 = alloca double
%4 = bitcast double* %3 to i8*
%5 = getelementptr inbounds double, double* %0, i64 3
%6 = load double, double* %5
%7 = getelementptr inbounds double, double* %0, i64 2
%8 = load double, double* %7
%9 = fcmp ogt double %6, %8
%10 = zext i1 %9 to i32
%11 = tail call i32 @assert(i32 %10)
%12 = load double, double* %7
%13 = getelementptr inbounds double, double* %0, i64 1
%14 = load double, double* %13
%15 = fcmp ogt double %12, %14
%16 = zext i1 %15 to i32
%17 = tail call i32 @assert(i32 %16)
%18 = getelementptr inbounds double, double* %1, i64 3
%19 = load double, double* %18
%20 = getelementptr inbounds double, double* %1, i64 2
%21 = load double, double* %20
%22 = fcmp oge double %19, %21
%23 = zext i1 %22 to i32
%24 = tail call i32 @assert(i32 %23)
%25 = getelementptr inbounds double, double* %1, i64 1
%26 = load double, double* %25
%27 = load double, double* %20
%28 = fcmp oge double %26, %27
%29 = zext i1 %28 to i32
%30 = tail call i32 @assert(i32 %29)
%31 = load double, double* %5
%32 = load double, double* %13
%33 = fsub double %31, %32
%34 = fcmp ogt double %33, 0x3E7AD7F29ABCAF48
br i1 %34, label %35, label %96
35: ; preds = %2
%36 = load double, double* %7
%37 = fsub double %36, %32
%38 = fcmp ogt double %37, 0x3E7AD7F29ABCAF48
%39 = fsub double %31, %36
%40 = fcmp ogt double %39, 0x3E7AD7F29ABCAF48
%41 = select i1 %38, i1 %40, i1 false
br i1 %41, label %42, label %96
42: ; preds = %35
%43 = call i64 @find_parabolic_min_x(double* %0, double* %1, double* %3)
%44 = icmp eq i64 %43, 0
br i1 %44, label %45, label %48
45: ; preds = %42
%46 = load double, double* %5
%47 = load double, double* %13
br label %96
48: ; preds = %42
%49 = load double, double* %3
%50 = load double, double* %13
%51 = fcmp ogt double %49, %50
br i1 %51, label %87, label %52
52: ; preds = %48
%53 = load i32, i32* @stderr
%54 = load double, double* %7
%55 = load double, double* %5
%56 = call i32 (i32, i8*, double, double, double, ...) @fprintf(i32 %53, i8* getelementptr inbounds ([47 x i8], [47 x i8]* @.str, i64 0, i64 0), double %50, double %54, double %55, double %49)
%57 = icmp eq i32 %56, 0
br i1 %57, label %87, label %58
58: ; preds = %52
%59 = load i32, i32* @stderr
%60 = load double, double* %25
%61 = load double, double* %20
%62 = load double, double* %18
%63 = call i32 (i32, i8*, double, double, double, ...) @fprintf(i32 %59, i8* getelementptr inbounds ([32 x i8], [32 x i8]* @.str.1, i64 0, i64 0), double %60, double %61, double %62)
%64 = icmp eq i32 %63, 0
br i1 %64, label %87, label %65
65: ; preds = %58
%66 = load i32, i32* @stderr
%67 = load double, double* %25
%68 = fsub double %67, %67
%69 = load double, double* %20
%70 = fsub double %69, %67
%71 = load double, double* %18
%72 = fsub double %71, %67
%73 = call i32 (i32, i8*, double, double, double, ...) @fprintf(i32 %66, i8* getelementptr inbounds ([35 x i8], [35 x i8]* @.str.2, i64 0, i64 0), double %68, double %70, double %72)
%74 = icmp eq i32 %73, 0
br i1 %74, label %87, label %75
75: ; preds = %65
%76 = load i32, i32* @stderr
%77 = load double, double* %5
%78 = load double, double* %7
%79 = fsub double %77, %78
%80 = load double, double* %13
%81 = fsub double %78, %80
%82 = load double, double* %3
%83 = fsub double %82, %80
%84 = call i32 (i32, i8*, double, double, double, ...) @fprintf(i32 %76, i8* getelementptr inbounds ([46 x i8], [46 x i8]* @.str.3, i64 0, i64 0), double %79, double %81, double %83)
%85 = icmp eq i32 %84, 0
%86 = zext i1 %85 to i32
br label %87
87: ; preds = %75, %65, %58, %52, %48
%88 = phi i32 [ 1, %65 ], [ 1, %58 ], [ 1, %52 ], [ 1, %48 ], [ %86, %75 ]
%89 = call i32 @assert(i32 %88)
%90 = load double, double* %3
%91 = load double, double* %5
%92 = fcmp olt double %90, %91
%93 = zext i1 %92 to i32
%94 = call i32 @assert(i32 %93)
%95 = load double, double* %3
br label %101
96: ; preds = %45, %35, %2
%97 = phi double [ %47, %45 ], [ %32, %35 ], [ %32, %2 ]
%98 = phi double [ %46, %45 ], [ %31, %35 ], [ %31, %2 ]
%99 = fadd double %98, %97
%100 = fmul double %99, 5.000000e-01
br label %101
101: ; preds = %96, %87
%102 = phi double [ %95, %87 ], [ %100, %96 ]
ret double %102
}