#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*/ isnan(double);
double doalarm(int N, double *resid, double *err) {
  double a, ave;
  long double sum1, sum2;
  int sign, i, n;
  sum1 = 0.;
  sum2 = 0.;
  for (i = 0; i < N; i++) {
    if (!isnan(resid[i]) && !isnan(err[i]) && err[i] > 0.) {
      sum1 += (long double)(resid[i] / (err[i] * err[i]));
      sum2 += (long double)(1. / (err[i] * err[i]));
    }
  }
  a = 0.;
  if (N > 0) {
    ave = sum1 / sum2;
    resid[0] -= ave;
    sum1 = resid[0] * resid[0] / (err[0] * err[0]);
    sum2 = resid[0] / err[0];
    n = 1;
    sign = ((resid[0]) >= 0 ? (1) : (-1));
    for (i = 1; i < N; i++) {
      if (!isnan(resid[i])) {
        resid[i] -= ave;
        if (((resid[i]) >= 0 ? (1) : (-1)) != sign) {
          a += (sum2 * sum2);
          sum2 = resid[i] / err[i];
          n = 1;
          sign = ((resid[i]) >= 0 ? (1) : (-1));
        } else {
          sum2 += resid[i] / err[i];
          n++;
        }
        sum1 += resid[i] * resid[i] / (err[i] * err[i]);
      }
    }
    a += (sum2 * sum2);
    a = (a / sum1) - 2.2732395;
  }
  return (a);
}


x86 assembly (clang, -O0):
<doalarm>:
push   %rbp
mov    %rsp,%rbp
sub    $0xa0,%rsp
mov    %edi,-0x4(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
fldz
fld    %st(0)
fstpt  -0x40(%rbp)
fstpt  -0x50(%rbp)
movl   $0x0,-0x58(%rbp)
mov    -0x58(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    120d <doalarm+0xfd>
mov    -0x10(%rbp),%rax
movslq -0x58(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
call   1030 <isnan@plt>
cmp    $0x0,%eax
jne    11fa <doalarm+0xea>
mov    -0x18(%rbp),%rax
movslq -0x58(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
call   1030 <isnan@plt>
cmp    $0x0,%eax
jne    11fa <doalarm+0xea>
mov    -0x18(%rbp),%rax
movslq -0x58(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jbe    11fa <doalarm+0xea>
mov    -0x10(%rbp),%rax
movslq -0x58(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x18(%rbp),%rax
movsd  (%rax,%rcx,8),%xmm1
mulsd  %xmm1,%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x98(%rbp)
ff
fldl   -0x98(%rbp)
fldt   -0x40(%rbp)
faddp  %st,%st(1)
fstpt  -0x40(%rbp)
mov    -0x18(%rbp),%rax
movslq -0x58(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
mulsd  %xmm1,%xmm1
movsd  0xe28(%rip),%xmm0
divsd  %xmm1,%xmm0
movsd  %xmm0,-0xa0(%rbp)
ff
fldl   -0xa0(%rbp)
fldt   -0x50(%rbp)
faddp  %st,%st(1)
fstpt  -0x50(%rbp)
jmp    11ff <doalarm+0xef>
mov    -0x58(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x58(%rbp)
jmp    1137 <doalarm+0x27>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x20(%rbp)
cmpl   $0x0,-0x4(%rbp)
jle    1432 <doalarm+0x322>
fldt   -0x40(%rbp)
fldt   -0x50(%rbp)
fdivrp %st,%st(1)
fstpl  -0x28(%rbp)
movsd  -0x28(%rbp),%xmm1
mov    -0x10(%rbp),%rax
movsd  (%rax),%xmm0
subsd  %xmm1,%xmm0
movsd  %xmm0,(%rax)
mov    -0x10(%rbp),%rax
movsd  (%rax),%xmm0
mulsd  %xmm0,%xmm0
mov    -0x18(%rbp),%rax
movsd  (%rax),%xmm1
mulsd  %xmm1,%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x68(%rbp)
fldl   -0x68(%rbp)
fstpt  -0x40(%rbp)
mov    -0x10(%rbp),%rax
movsd  (%rax),%xmm0
mov    -0x18(%rbp),%rax
movsd  (%rax),%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x70(%rbp)
fldl   -0x70(%rbp)
fstpt  -0x50(%rbp)
movl   $0x1,-0x5c(%rbp)
mov    -0x10(%rbp),%rax
movsd  (%rax),%xmm0
xorps  %xmm1,%xmm1
mov    $0xffffffff,%eax
mov    $0x1,%ecx
ucomisd %xmm1,%xmm0
cmovae %ecx,%eax
mov    %eax,-0x54(%rbp)
movl   $0x1,-0x58(%rbp)
mov    -0x58(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    1406 <doalarm+0x2f6>
mov    -0x10(%rbp),%rax
movslq -0x58(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
call   1030 <isnan@plt>
cmp    $0x0,%eax
jne    13f3 <doalarm+0x2e3>
movsd  -0x28(%rbp),%xmm1
mov    -0x10(%rbp),%rax
movslq -0x58(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
subsd  %xmm1,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x10(%rbp),%rax
movslq -0x58(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
xorps  %xmm1,%xmm1
mov    $0xffffffff,%eax
mov    $0x1,%ecx
ucomisd %xmm1,%xmm0
cmovae %ecx,%eax
cmp    -0x54(%rbp),%eax
je     1388 <doalarm+0x278>
fldt   -0x50(%rbp)
fldl   -0x20(%rbp)
fxch   %st(1)
fmul   %st,%st(0)
faddp  %st,%st(1)
fstpl  -0x20(%rbp)
mov    -0x10(%rbp),%rax
movslq -0x58(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x18(%rbp),%rax
movsd  (%rax,%rcx,8),%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x88(%rbp)
ff
fldl   -0x88(%rbp)
fstpt  -0x50(%rbp)
movl   $0x1,-0x5c(%rbp)
mov    -0x10(%rbp),%rax
movslq -0x58(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
xorps  %xmm1,%xmm1
mov    $0xffffffff,%eax
mov    $0x1,%ecx
ucomisd %xmm1,%xmm0
cmovae %ecx,%eax
mov    %eax,-0x54(%rbp)
jmp    13bb <doalarm+0x2ab>
mov    -0x10(%rbp),%rax
movslq -0x58(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x18(%rbp),%rax
movsd  (%rax,%rcx,8),%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x80(%rbp)
fldl   -0x80(%rbp)
fldt   -0x50(%rbp)
faddp  %st,%st(1)
fstpt  -0x50(%rbp)
mov    -0x5c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x5c(%rbp)
mov    -0x10(%rbp),%rax
movslq -0x58(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mulsd  %xmm0,%xmm0
mov    -0x18(%rbp),%rax
movsd  (%rax,%rcx,8),%xmm1
mulsd  %xmm1,%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x90(%rbp)
ff
fldl   -0x90(%rbp)
fldt   -0x40(%rbp)
faddp  %st,%st(1)
fstpt  -0x40(%rbp)
jmp    13f8 <doalarm+0x2e8>
mov    -0x58(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x58(%rbp)
jmp    12b2 <doalarm+0x1a2>
fldt   -0x50(%rbp)
fldl   -0x20(%rbp)
fxch   %st(1)
fmul   %st,%st(0)
faddp  %st,%st(1)
fstpl  -0x20(%rbp)
fldl   -0x20(%rbp)
fldt   -0x40(%rbp)
fdivrp %st,%st(1)
fldl   0xbdd(%rip)
faddp  %st,%st(1)
fstpl  -0x78(%rbp)
movsd  -0x78(%rbp),%xmm0
movsd  %xmm0,-0x20(%rbp)
movsd  -0x20(%rbp),%xmm0
add    $0xa0,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<doalarm>:
push   %rbp
push   %r15
push   %r14
push   %r12
push   %rbx
sub    $0xa0,%rsp
test   %edi,%edi
jle    1298 <doalarm+0x188>
mov    %rdx,%r14
mov    %rsi,%rbx
mov    %edi,%r15d
mov    %edi,%r12d
fldz
xor    %ebp,%ebp
fld    %st(0)
fstpt  0x14(%rsp)
fstpt  0x8(%rsp)
jmp    1159 <doalarm+0x49>
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
add    $0x1,%rbp
cmp    %rbp,%r12
je     11d4 <doalarm+0xc4>
movsd  (%rbx,%rbp,8),%xmm0
call   1030 <isnan@plt>
test   %eax,%eax
jne    1150 <doalarm+0x40>
movsd  (%r14,%rbp,8),%xmm0
call   1030 <isnan@plt>
test   %eax,%eax
jne    1150 <doalarm+0x40>
movsd  (%r14,%rbp,8),%xmm0
ucomisd 0xe8c(%rip),%xmm0
jbe    1150 <doalarm+0x40>
movsd  (%rbx,%rbp,8),%xmm1
mulsd  %xmm0,%xmm0
divsd  %xmm0,%xmm1
movsd  %xmm1,0x90(%rsp)
fldt   0x8(%rsp)
faddl  0x90(%rsp)
fstpt  0x8(%rsp)
movsd  0xe4d(%rip),%xmm1
divsd  %xmm0,%xmm1
movsd  %xmm1,0x98(%rsp)
fldt   0x14(%rsp)
faddl  0x98(%rsp)
fstpt  0x14(%rsp)
jmp    1150 <doalarm+0x40>
test   %r15d,%r15d
jle    1298 <doalarm+0x188>
fldt   0x14(%rsp)
fldt   0x8(%rsp)
fdivp  %st,%st(1)
fstpl  0x88(%rsp)
movsd  0x88(%rsp),%xmm1
movsd  (%rbx),%xmm0
movsd  %xmm1,0x28(%rsp)
subsd  %xmm1,%xmm0
movsd  %xmm0,(%rbx)
movapd %xmm0,%xmm1
mulsd  %xmm0,%xmm1
movsd  (%r14),%xmm2
movapd %xmm0,%xmm3
divsd  %xmm2,%xmm3
mulsd  %xmm2,%xmm2
divsd  %xmm2,%xmm1
movsd  %xmm1,0x78(%rsp)
fldl   0x78(%rsp)
fstpt  0x14(%rsp)
movsd  %xmm3,0x80(%rsp)
fldl   0x80(%rsp)
fstpt  0x8(%rsp)
xorpd  %xmm1,%xmm1
movsd  %xmm1,0x20(%rsp)
cmp    $0x1,%r15d
jne    12ac <doalarm+0x19c>
movsd  0x20(%rsp),%xmm0
movsd  %xmm0,0x48(%rsp)
fldt   0x8(%rsp)
fmul   %st,%st(0)
faddl  0x48(%rsp)
fstpl  0x40(%rsp)
movsd  0x40(%rsp),%xmm0
movsd  %xmm0,0x38(%rsp)
fldt   0x14(%rsp)
fdivrl 0x38(%rsp)
faddl  0xd7c(%rip)
fstpl  0x30(%rsp)
movsd  0x30(%rsp),%xmm0
jmp    129c <doalarm+0x18c>
xorpd  %xmm0,%xmm0
add    $0xa0,%rsp
pop    %rbx
pop    %r12
pop    %r14
pop    %r15
pop    %rbp
ret
xorpd  %xmm1,%xmm1
xor    %eax,%eax
ucomisd %xmm1,%xmm0
setae  %al
lea    (%rax,%rax,1),%r15d
add    $0xffffffff,%r15d
mov    $0x1,%ebp
xorpd  %xmm0,%xmm0
movsd  %xmm0,0x20(%rsp)
jmp    131d <doalarm+0x20d>
movsd  (%r14,%rbp,8),%xmm1
movapd %xmm0,%xmm2
divsd  %xmm1,%xmm2
movsd  %xmm2,0x58(%rsp)
fldt   0x8(%rsp)
faddl  0x58(%rsp)
fstpt  0x8(%rsp)
mulsd  %xmm0,%xmm0
mulsd  %xmm1,%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,0x50(%rsp)
fldt   0x14(%rsp)
faddl  0x50(%rsp)
fstpt  0x14(%rsp)
add    $0x1,%rbp
cmp    %rbp,%r12
je     1258 <doalarm+0x148>
movsd  (%rbx,%rbp,8),%xmm0
call   1030 <isnan@plt>
test   %eax,%eax
jne    1310 <doalarm+0x200>
movsd  (%rbx,%rbp,8),%xmm0
subsd  0x28(%rsp),%xmm0
xor    %eax,%eax
ucomisd 0xcd0(%rip),%xmm0
movsd  %xmm0,(%rbx,%rbp,8)
setae  %al
add    %eax,%eax
add    $0xffffffff,%eax
cmp    %r15d,%eax
je     12d2 <doalarm+0x1c2>
movsd  0x20(%rsp),%xmm1
movsd  %xmm1,0x70(%rsp)
fldt   0x8(%rsp)
fmul   %st,%st(0)
faddl  0x70(%rsp)
fstpl  0x68(%rsp)
movsd  0x68(%rsp),%xmm1
movsd  %xmm1,0x20(%rsp)
movsd  (%r14,%rbp,8),%xmm1
movapd %xmm0,%xmm2
divsd  %xmm1,%xmm2
movsd  %xmm2,0x60(%rsp)
fldl   0x60(%rsp)
fstpt  0x8(%rsp)
mov    %eax,%r15d
jmp    12f2 <doalarm+0x1e2>

LLVM IR (-O0):
define double @doalarm(i32 %0, double* %1, double* %2) {
%4 = alloca i32
%5 = alloca double*
%6 = alloca double*
%7 = alloca double
%8 = alloca double
%9 = alloca x86_fp80
%10 = alloca x86_fp80
%11 = alloca i32
%12 = alloca i32
%13 = alloca i32
store i32 %0, i32* %4
store double* %1, double** %5
store double* %2, double** %6
store x86_fp80 0xK00000000000000000000, x86_fp80* %9
store x86_fp80 0xK00000000000000000000, x86_fp80* %10
store i32 0, i32* %12
br label %14
14: ; preds = %78, %3
%15 = load i32, i32* %12
%16 = load i32, i32* %4
%17 = icmp slt i32 %15, %16
br i1 %17, label %18, label %81
18: ; preds = %14
%19 = load double*, double** %5
%20 = load i32, i32* %12
%21 = sext i32 %20 to i64
%22 = getelementptr inbounds double, double* %19, i64 %21
%23 = load double, double* %22
%24 = call i32 @isnan(double %23)
%25 = icmp ne i32 %24, 0
br i1 %25, label %77, label %26
26: ; preds = %18
%27 = load double*, double** %6
%28 = load i32, i32* %12
%29 = sext i32 %28 to i64
%30 = getelementptr inbounds double, double* %27, i64 %29
%31 = load double, double* %30
%32 = call i32 @isnan(double %31)
%33 = icmp ne i32 %32, 0
br i1 %33, label %77, label %34
34: ; preds = %26
%35 = load double*, double** %6
%36 = load i32, i32* %12
%37 = sext i32 %36 to i64
%38 = getelementptr inbounds double, double* %35, i64 %37
%39 = load double, double* %38
%40 = fcmp ogt double %39, 0.000000e+00
br i1 %40, label %41, label %77
41: ; preds = %34
%42 = load double*, double** %5
%43 = load i32, i32* %12
%44 = sext i32 %43 to i64
%45 = getelementptr inbounds double, double* %42, i64 %44
%46 = load double, double* %45
%47 = load double*, double** %6
%48 = load i32, i32* %12
%49 = sext i32 %48 to i64
%50 = getelementptr inbounds double, double* %47, i64 %49
%51 = load double, double* %50
%52 = load double*, double** %6
%53 = load i32, i32* %12
%54 = sext i32 %53 to i64
%55 = getelementptr inbounds double, double* %52, i64 %54
%56 = load double, double* %55
%57 = fmul double %51, %56
%58 = fdiv double %46, %57
%59 = fpext double %58 to x86_fp80
%60 = load x86_fp80, x86_fp80* %9
%61 = fadd x86_fp80 %60, %59
store x86_fp80 %61, x86_fp80* %9
%62 = load double*, double** %6
%63 = load i32, i32* %12
%64 = sext i32 %63 to i64
%65 = getelementptr inbounds double, double* %62, i64 %64
%66 = load double, double* %65
%67 = load double*, double** %6
%68 = load i32, i32* %12
%69 = sext i32 %68 to i64
%70 = getelementptr inbounds double, double* %67, i64 %69
%71 = load double, double* %70
%72 = fmul double %66, %71
%73 = fdiv double 1.000000e+00, %72
%74 = fpext double %73 to x86_fp80
%75 = load x86_fp80, x86_fp80* %10
%76 = fadd x86_fp80 %75, %74
store x86_fp80 %76, x86_fp80* %10
br label %77
77: ; preds = %41, %34, %26, %18
br label %78
78: ; preds = %77
%79 = load i32, i32* %12
%80 = add nsw i32 %79, 1
store i32 %80, i32* %12
br label %14
81: ; preds = %14
store double 0.000000e+00, double* %7
%82 = load i32, i32* %4
%83 = icmp sgt i32 %82, 0
br i1 %83, label %84, label %242
84: ; preds = %81
%85 = load x86_fp80, x86_fp80* %9
%86 = load x86_fp80, x86_fp80* %10
%87 = fdiv x86_fp80 %85, %86
%88 = fptrunc x86_fp80 %87 to double
store double %88, double* %8
%89 = load double, double* %8
%90 = load double*, double** %5
%91 = getelementptr inbounds double, double* %90, i64 0
%92 = load double, double* %91
%93 = fsub double %92, %89
store double %93, double* %91
%94 = load double*, double** %5
%95 = getelementptr inbounds double, double* %94, i64 0
%96 = load double, double* %95
%97 = load double*, double** %5
%98 = getelementptr inbounds double, double* %97, i64 0
%99 = load double, double* %98
%100 = fmul double %96, %99
%101 = load double*, double** %6
%102 = getelementptr inbounds double, double* %101, i64 0
%103 = load double, double* %102
%104 = load double*, double** %6
%105 = getelementptr inbounds double, double* %104, i64 0
%106 = load double, double* %105
%107 = fmul double %103, %106
%108 = fdiv double %100, %107
%109 = fpext double %108 to x86_fp80
store x86_fp80 %109, x86_fp80* %9
%110 = load double*, double** %5
%111 = getelementptr inbounds double, double* %110, i64 0
%112 = load double, double* %111
%113 = load double*, double** %6
%114 = getelementptr inbounds double, double* %113, i64 0
%115 = load double, double* %114
%116 = fdiv double %112, %115
%117 = fpext double %116 to x86_fp80
store x86_fp80 %117, x86_fp80* %10
store i32 1, i32* %13
%118 = load double*, double** %5
%119 = getelementptr inbounds double, double* %118, i64 0
%120 = load double, double* %119
%121 = fcmp oge double %120, 0.000000e+00
%122 = zext i1 %121 to i64
%123 = select i1 %121, i32 1, i32 -1
store i32 %123, i32* %11
store i32 1, i32* %12
br label %124
124: ; preds = %226, %84
%125 = load i32, i32* %12
%126 = load i32, i32* %4
%127 = icmp slt i32 %125, %126
br i1 %127, label %128, label %229
128: ; preds = %124
%129 = load double*, double** %5
%130 = load i32, i32* %12
%131 = sext i32 %130 to i64
%132 = getelementptr inbounds double, double* %129, i64 %131
%133 = load double, double* %132
%134 = call i32 @isnan(double %133)
%135 = icmp ne i32 %134, 0
br i1 %135, label %225, label %136
136: ; preds = %128
%137 = load double, double* %8
%138 = load double*, double** %5
%139 = load i32, i32* %12
%140 = sext i32 %139 to i64
%141 = getelementptr inbounds double, double* %138, i64 %140
%142 = load double, double* %141
%143 = fsub double %142, %137
store double %143, double* %141
%144 = load double*, double** %5
%145 = load i32, i32* %12
%146 = sext i32 %145 to i64
%147 = getelementptr inbounds double, double* %144, i64 %146
%148 = load double, double* %147
%149 = fcmp oge double %148, 0.000000e+00
%150 = zext i1 %149 to i64
%151 = select i1 %149, i32 1, i32 -1
%152 = load i32, i32* %11
%153 = icmp ne i32 %151, %152
br i1 %153, label %154, label %181
154: ; preds = %136
%155 = load x86_fp80, x86_fp80* %10
%156 = load x86_fp80, x86_fp80* %10
%157 = load double, double* %7
%158 = fpext double %157 to x86_fp80
%159 = call x86_fp80 @llvm.fmuladd.f80(x86_fp80 %155, x86_fp80 %156, x86_fp80 %158)
%160 = fptrunc x86_fp80 %159 to double
store double %160, double* %7
%161 = load double*, double** %5
%162 = load i32, i32* %12
%163 = sext i32 %162 to i64
%164 = getelementptr inbounds double, double* %161, i64 %163
%165 = load double, double* %164
%166 = load double*, double** %6
%167 = load i32, i32* %12
%168 = sext i32 %167 to i64
%169 = getelementptr inbounds double, double* %166, i64 %168
%170 = load double, double* %169
%171 = fdiv double %165, %170
%172 = fpext double %171 to x86_fp80
store x86_fp80 %172, x86_fp80* %10
store i32 1, i32* %13
%173 = load double*, double** %5
%174 = load i32, i32* %12
%175 = sext i32 %174 to i64
%176 = getelementptr inbounds double, double* %173, i64 %175
%177 = load double, double* %176
%178 = fcmp oge double %177, 0.000000e+00
%179 = zext i1 %178 to i64
%180 = select i1 %178, i32 1, i32 -1
store i32 %180, i32* %11
br label %198
181: ; preds = %136
%182 = load double*, double** %5
%183 = load i32, i32* %12
%184 = sext i32 %183 to i64
%185 = getelementptr inbounds double, double* %182, i64 %184
%186 = load double, double* %185
%187 = load double*, double** %6
%188 = load i32, i32* %12
%189 = sext i32 %188 to i64
%190 = getelementptr inbounds double, double* %187, i64 %189
%191 = load double, double* %190
%192 = fdiv double %186, %191
%193 = fpext double %192 to x86_fp80
%194 = load x86_fp80, x86_fp80* %10
%195 = fadd x86_fp80 %194, %193
store x86_fp80 %195, x86_fp80* %10
%196 = load i32, i32* %13
%197 = add nsw i32 %196, 1
store i32 %197, i32* %13
br label %198
198: ; preds = %181, %154
%199 = load double*, double** %5
%200 = load i32, i32* %12
%201 = sext i32 %200 to i64
%202 = getelementptr inbounds double, double* %199, i64 %201
%203 = load double, double* %202
%204 = load double*, double** %5
%205 = load i32, i32* %12
%206 = sext i32 %205 to i64
%207 = getelementptr inbounds double, double* %204, i64 %206
%208 = load double, double* %207
%209 = fmul double %203, %208
%210 = load double*, double** %6
%211 = load i32, i32* %12
%212 = sext i32 %211 to i64
%213 = getelementptr inbounds double, double* %210, i64 %212
%214 = load double, double* %213
%215 = load double*, double** %6
%216 = load i32, i32* %12
%217 = sext i32 %216 to i64
%218 = getelementptr inbounds double, double* %215, i64 %217
%219 = load double, double* %218
%220 = fmul double %214, %219
%221 = fdiv double %209, %220
%222 = fpext double %221 to x86_fp80
%223 = load x86_fp80, x86_fp80* %9
%224 = fadd x86_fp80 %223, %222
store x86_fp80 %224, x86_fp80* %9
br label %225
225: ; preds = %198, %128
br label %226
226: ; preds = %225
%227 = load i32, i32* %12
%228 = add nsw i32 %227, 1
store i32 %228, i32* %12
br label %124
229: ; preds = %124
%230 = load x86_fp80, x86_fp80* %10
%231 = load x86_fp80, x86_fp80* %10
%232 = load double, double* %7
%233 = fpext double %232 to x86_fp80
%234 = call x86_fp80 @llvm.fmuladd.f80(x86_fp80 %230, x86_fp80 %231, x86_fp80 %233)
%235 = fptrunc x86_fp80 %234 to double
store double %235, double* %7
%236 = load double, double* %7
%237 = fpext double %236 to x86_fp80
%238 = load x86_fp80, x86_fp80* %9
%239 = fdiv x86_fp80 %237, %238
%240 = fsub x86_fp80 %239, 0xK4000917CC1871E6CD000
%241 = fptrunc x86_fp80 %240 to double
store double %241, double* %7
br label %242
242: ; preds = %229, %81
%243 = load double, double* %7
ret double %243
}

LLVM IR (-O2):
define double @doalarm(i32 %0, double* nocapture %1, double* nocapture %2) {
%4 = icmp sgt i32 %0, 0
br i1 %4, label %5, label %113
5: ; preds = %3
%6 = zext i32 %0 to i64
br label %7
7: ; preds = %5, %32
%8 = phi i64 [ 0, %5 ], [ %35, %32 ]
%9 = phi x86_fp80 [ 0xK00000000000000000000, %5 ], [ %34, %32 ]
%10 = phi x86_fp80 [ 0xK00000000000000000000, %5 ], [ %33, %32 ]
%11 = getelementptr inbounds double, double* %1, i64 %8
%12 = load double, double* %11
%13 = tail call i32 @isnan(double %12)
%14 = icmp eq i32 %13, 0
br i1 %14, label %15, label %32
15: ; preds = %7
%16 = getelementptr inbounds double, double* %2, i64 %8
%17 = load double, double* %16
%18 = tail call i32 @isnan(double %17)
%19 = icmp eq i32 %18, 0
br i1 %19, label %20, label %32
20: ; preds = %15
%21 = load double, double* %16
%22 = fcmp ogt double %21, 0.000000e+00
br i1 %22, label %23, label %32
23: ; preds = %20
%24 = load double, double* %11
%25 = fmul double %21, %21
%26 = fdiv double %24, %25
%27 = fpext double %26 to x86_fp80
%28 = fadd x86_fp80 %10, %27
%29 = fdiv double 1.000000e+00, %25
%30 = fpext double %29 to x86_fp80
%31 = fadd x86_fp80 %9, %30
br label %32
32: ; preds = %7, %15, %20, %23
%33 = phi x86_fp80 [ %10, %7 ], [ %10, %15 ], [ %28, %23 ], [ %10, %20 ]
%34 = phi x86_fp80 [ %9, %7 ], [ %9, %15 ], [ %31, %23 ], [ %9, %20 ]
%35 = add nuw nsw i64 %8, 1
%36 = icmp eq i64 %35, %6
br i1 %36, label %37, label %7
37: ; preds = %32
br i1 %4, label %38, label %113
38: ; preds = %37
%39 = fdiv x86_fp80 %33, %34
%40 = fptrunc x86_fp80 %39 to double
%41 = load double, double* %1
%42 = fsub double %41, %40
store double %42, double* %1
%43 = fmul double %42, %42
%44 = load double, double* %2
%45 = fmul double %44, %44
%46 = fdiv double %43, %45
%47 = fpext double %46 to x86_fp80
%48 = fdiv double %42, %44
%49 = fpext double %48 to x86_fp80
%50 = icmp eq i32 %0, 1
br i1 %50, label %102, label %51
51: ; preds = %38
%52 = fcmp oge double %42, 0.000000e+00
%53 = select i1 %52, i32 1, i32 -1
%54 = zext i32 %0 to i64
br label %55
55: ; preds = %51, %95
%56 = phi i64 [ 1, %51 ], [ %100, %95 ]
%57 = phi i32 [ %53, %51 ], [ %99, %95 ]
%58 = phi x86_fp80 [ %49, %51 ], [ %98, %95 ]
%59 = phi double [ 0.000000e+00, %51 ], [ %97, %95 ]
%60 = phi x86_fp80 [ %47, %51 ], [ %96, %95 ]
%61 = getelementptr inbounds double, double* %1, i64 %56
%62 = load double, double* %61
%63 = tail call i32 @isnan(double %62)
%64 = icmp eq i32 %63, 0
br i1 %64, label %65, label %95
65: ; preds = %55
%66 = load double, double* %61
%67 = fsub double %66, %40
store double %67, double* %61
%68 = fcmp oge double %67, 0.000000e+00
%69 = select i1 %68, i32 1, i32 -1
%70 = icmp eq i32 %69, %57
br i1 %70, label %79, label %71
71: ; preds = %65
%72 = fpext double %59 to x86_fp80
%73 = tail call x86_fp80 @llvm.fmuladd.f80(x86_fp80 %58, x86_fp80 %58, x86_fp80 %72)
%74 = fptrunc x86_fp80 %73 to double
%75 = getelementptr inbounds double, double* %2, i64 %56
%76 = load double, double* %75
%77 = fdiv double %67, %76
%78 = fpext double %77 to x86_fp80
br label %85
79: ; preds = %65
%80 = getelementptr inbounds double, double* %2, i64 %56
%81 = load double, double* %80
%82 = fdiv double %67, %81
%83 = fpext double %82 to x86_fp80
%84 = fadd x86_fp80 %58, %83
br label %85
85: ; preds = %79, %71
%86 = phi double [ %76, %71 ], [ %81, %79 ]
%87 = phi double [ %74, %71 ], [ %59, %79 ]
%88 = phi x86_fp80 [ %78, %71 ], [ %84, %79 ]
%89 = phi i32 [ %69, %71 ], [ %57, %79 ]
%90 = fmul double %67, %67
%91 = fmul double %86, %86
%92 = fdiv double %90, %91
%93 = fpext double %92 to x86_fp80
%94 = fadd x86_fp80 %60, %93
br label %95
95: ; preds = %55, %85
%96 = phi x86_fp80 [ %60, %55 ], [ %94, %85 ]
%97 = phi double [ %59, %55 ], [ %87, %85 ]
%98 = phi x86_fp80 [ %58, %55 ], [ %88, %85 ]
%99 = phi i32 [ %57, %55 ], [ %89, %85 ]
%100 = add nuw nsw i64 %56, 1
%101 = icmp eq i64 %100, %54
br i1 %101, label %102, label %55
102: ; preds = %95, %38
%103 = phi x86_fp80 [ %47, %38 ], [ %96, %95 ]
%104 = phi double [ 0.000000e+00, %38 ], [ %97, %95 ]
%105 = phi x86_fp80 [ %49, %38 ], [ %98, %95 ]
%106 = fpext double %104 to x86_fp80
%107 = tail call x86_fp80 @llvm.fmuladd.f80(x86_fp80 %105, x86_fp80 %105, x86_fp80 %106)
%108 = fptrunc x86_fp80 %107 to double
%109 = fpext double %108 to x86_fp80
%110 = fdiv x86_fp80 %109, %103
%111 = fadd x86_fp80 %110, 0xKC000917CC1871E6CD000
%112 = fptrunc x86_fp80 %111 to double
br label %113
113: ; preds = %3, %102, %37
%114 = phi double [ %112, %102 ], [ 0.000000e+00, %37 ], [ 0.000000e+00, %3 ]
ret double %114
}