#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*/ isinf(double);
int /*<<< orphan*/ isnan(double);
double chisqr_sysrem(double p[], int numlc, int numJD, double **mag,
                     double **sig, double *mag_ave, double sigma_clip,
                     double saturation) {
  int j, k, n;
  long double S, dum;
  S = 0.0;
  n = 0;
  for (j = 0; j < numlc; j++)
    for (k = 0; k < numJD; k++)
      if (!isnan(mag[j][k]) && !isinf(mag[j][k]) && mag[j][k] > saturation &&
          (((mag[j][k] - mag_ave[j]) > 0 ? (mag[j][k] - mag_ave[j])
                                         : (-(mag[j][k] - mag_ave[j]))) /
                   sig[j][k] <
               sigma_clip ||
           sigma_clip < 0)) {
        dum = (long double)(mag[j][k] - mag_ave[j] -
                            p[j + 1] * p[numlc + k + 1]) /
              sig[j][k];
        S += dum * dum;
      }
  return ((double)(S / (double)(numJD * numlc)));
}


x86 assembly (clang, -O0):
<chisqr_sysrem>:
push   %rbp
mov    %rsp,%rbp
sub    $0x90,%rsp
mov    %rdi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
mov    %edx,-0x10(%rbp)
mov    %rcx,-0x18(%rbp)
mov    %r8,-0x20(%rbp)
mov    %r9,-0x28(%rbp)
movsd  %xmm0,-0x30(%rbp)
movsd  %xmm1,-0x38(%rbp)
fldz
fstpt  -0x60(%rbp)
movl   $0x0,-0x44(%rbp)
movl   $0x0,-0x3c(%rbp)
mov    -0x3c(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    1353 <chisqr_sysrem+0x233>
movl   $0x0,-0x40(%rbp)
mov    -0x40(%rbp),%eax
cmp    -0x10(%rbp),%eax
jge    1340 <chisqr_sysrem+0x220>
mov    -0x18(%rbp),%rax
movslq -0x3c(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x40(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
call   1030 <isnan@plt>
cmp    $0x0,%eax
jne    132d <chisqr_sysrem+0x20d>
mov    -0x18(%rbp),%rax
movslq -0x3c(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x40(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
call   1040 <isinf@plt>
cmp    $0x0,%eax
jne    132d <chisqr_sysrem+0x20d>
mov    -0x18(%rbp),%rax
movslq -0x3c(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x40(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
ucomisd -0x38(%rbp),%xmm0
jbe    132d <chisqr_sysrem+0x20d>
mov    -0x18(%rbp),%rax
movslq -0x3c(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x40(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x28(%rbp),%rax
movslq -0x3c(%rbp),%rcx
subsd  (%rax,%rcx,8),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jbe    1241 <chisqr_sysrem+0x121>
mov    -0x18(%rbp),%rax
movslq -0x3c(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x40(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x28(%rbp),%rax
movslq -0x3c(%rbp),%rcx
subsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x90(%rbp)
ff
jmp    1282 <chisqr_sysrem+0x162>
mov    -0x18(%rbp),%rax
movslq -0x3c(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x40(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x28(%rbp),%rax
movslq -0x3c(%rbp),%rcx
subsd  (%rax,%rcx,8),%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
movsd  %xmm0,-0x90(%rbp)
ff
movsd  -0x90(%rbp),%xmm1
ff
mov    -0x20(%rbp),%rax
movslq -0x3c(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x40(%rbp),%rcx
divsd  (%rax,%rcx,8),%xmm1
movsd  -0x30(%rbp),%xmm0
ucomisd %xmm1,%xmm0
ja     12bc <chisqr_sysrem+0x19c>
xorps  %xmm0,%xmm0
ucomisd -0x30(%rbp),%xmm0
jbe    132d <chisqr_sysrem+0x20d>
mov    -0x18(%rbp),%rax
movslq -0x3c(%rbp),%rdx
mov    (%rax,%rdx,8),%rax
movslq -0x40(%rbp),%rcx
mov    %ecx,%edi
movsd  (%rax,%rcx,8),%xmm0
mov    -0x28(%rbp),%rax
movsd  (%rax,%rdx,8),%xmm1
subsd  %xmm1,%xmm0
mov    -0x8(%rbp),%rax
movsd  0x8(%rax,%rdx,8),%xmm1
mov    -0xc(%rbp),%esi
add    %edi,%esi
movslq %esi,%rsi
movsd  0x8(%rax,%rsi,8),%xmm2
mulsd  %xmm2,%xmm1
subsd  %xmm1,%xmm0
movsd  %xmm0,-0x88(%rbp)
ff
fldl   -0x88(%rbp)
mov    -0x20(%rbp),%rax
mov    (%rax,%rdx,8),%rax
fldl   (%rax,%rcx,8)
fdivrp %st,%st(1)
fstpt  -0x70(%rbp)
fldt   -0x70(%rbp)
fldt   -0x60(%rbp)
fxch   %st(1)
fmul   %st,%st(0)
faddp  %st,%st(1)
fstpt  -0x60(%rbp)
jmp    1332 <chisqr_sysrem+0x212>
mov    -0x40(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x40(%rbp)
jmp    1171 <chisqr_sysrem+0x51>
jmp    1345 <chisqr_sysrem+0x225>
mov    -0x3c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x3c(%rbp)
jmp    115e <chisqr_sysrem+0x3e>
fldt   -0x60(%rbp)
mov    -0x10(%rbp),%eax
mov    -0xc(%rbp),%ecx
imul   %ecx,%eax
cvtsi2sd %eax,%xmm0
movsd  %xmm0,-0x78(%rbp)
fldl   -0x78(%rbp)
fdivrp %st,%st(1)
fstpl  -0x80(%rbp)
movsd  -0x80(%rbp),%xmm0
add    $0x90,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<chisqr_sysrem>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x68,%rsp
movsd  %xmm1,0x48(%rsp)
mov    %r9,0x40(%rsp)
mov    %r8,0x38(%rsp)
mov    %edx,0x4(%rsp)
fldz
fstpt  0x8(%rsp)
mov    %rsi,0x18(%rsp)
test   %esi,%esi
jle    12b2 <chisqr_sysrem+0x192>
mov    %rcx,%rbx
mov    %rdi,%rbp
mov    0x18(%rsp),%rax
mov    %eax,%ecx
mov    %rcx,0x28(%rsp)
mov    0x4(%rsp),%r14d
add    $0x1,%eax
mov    %rax,0x20(%rsp)
fldz
fstpt  0x8(%rsp)
xor    %r15d,%r15d
movsd  %xmm0,0x30(%rsp)
jmp    119e <chisqr_sysrem+0x7e>
cs nopw 0x0(%rax,%rax,1)
nop
mov    %r12,%r15
cmp    0x28(%rsp),%r12
je     12b2 <chisqr_sysrem+0x192>
cmpl   $0x0,0x4(%rsp)
jle    12a0 <chisqr_sysrem+0x180>
lea    0x1(%r15),%r12
xor    %r13d,%r13d
jmp    11f9 <chisqr_sysrem+0xd9>
movsd  0x8(%rbp,%r15,8),%xmm2
mov    0x20(%rsp),%rax
add    %r13d,%eax
cltq
mulsd  0x0(%rbp,%rax,8),%xmm2
subsd  %xmm2,%xmm5
movsd  %xmm5,0x58(%rsp)
fldl   0x58(%rsp)
movsd  %xmm1,0x60(%rsp)
fdivl  0x60(%rsp)
fmul   %st,%st(0)
fldt   0x8(%rsp)
faddp  %st,%st(1)
fstpt  0x8(%rsp)
nopl   (%rax)
add    $0x1,%r13
cmp    %r13,%r14
je     1190 <chisqr_sysrem+0x70>
mov    (%rbx,%r15,8),%rax
movsd  (%rax,%r13,8),%xmm0
call   1030 <isnan@plt>
test   %eax,%eax
jne    11f0 <chisqr_sysrem+0xd0>
mov    (%rbx,%r15,8),%rax
movsd  (%rax,%r13,8),%xmm0
call   1040 <isinf@plt>
test   %eax,%eax
jne    11f0 <chisqr_sysrem+0xd0>
mov    (%rbx,%r15,8),%rax
movsd  (%rax,%r13,8),%xmm5
ucomisd 0x48(%rsp),%xmm5
jbe    11f0 <chisqr_sysrem+0xd0>
movsd  0x30(%rsp),%xmm0
xorpd  %xmm1,%xmm1
ucomisd %xmm0,%xmm1
mov    0x40(%rsp),%rax
subsd  (%rax,%r15,8),%xmm5
mov    0x38(%rsp),%rax
mov    (%rax,%r15,8),%rax
movsd  (%rax,%r13,8),%xmm1
ja     11b2 <chisqr_sysrem+0x92>
movapd %xmm5,%xmm2
xorpd  0xd95(%rip),%xmm2
xorpd  %xmm3,%xmm3
cmpltsd %xmm5,%xmm3
movapd %xmm3,%xmm4
andpd  %xmm5,%xmm4
andnpd %xmm2,%xmm3
orpd   %xmm4,%xmm3
divsd  %xmm1,%xmm3
ucomisd %xmm3,%xmm0
jbe    11f0 <chisqr_sysrem+0xd0>
jmp    11b2 <chisqr_sysrem+0x92>
nopw   0x0(%rax,%rax,1)
add    $0x1,%r15
mov    %r15,%r12
cmp    0x28(%rsp),%r12
jne    119e <chisqr_sysrem+0x7e>
mov    0x4(%rsp),%eax
imul   0x18(%rsp),%eax
mov    %eax,0x14(%rsp)
fldt   0x8(%rsp)
fidivl 0x14(%rsp)
fstpl  0x50(%rsp)
movsd  0x50(%rsp),%xmm0
add    $0x68,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

LLVM IR (-O0):
define double @chisqr_sysrem(double* %0, i32 %1, i32 %2, double** %3, double** %4, double* %5, double %6, double %7) {
%9 = alloca double*
%10 = alloca i32
%11 = alloca i32
%12 = alloca double**
%13 = alloca double**
%14 = alloca double*
%15 = alloca double
%16 = alloca double
%17 = alloca i32
%18 = alloca i32
%19 = alloca i32
%20 = alloca x86_fp80
%21 = alloca x86_fp80
store double* %0, double** %9
store i32 %1, i32* %10
store i32 %2, i32* %11
store double** %3, double*** %12
store double** %4, double*** %13
store double* %5, double** %14
store double %6, double* %15
store double %7, double* %16
store x86_fp80 0xK00000000000000000000, x86_fp80* %20
store i32 0, i32* %19
store i32 0, i32* %17
br label %22
22: ; preds = %187, %8
%23 = load i32, i32* %17
%24 = load i32, i32* %10
%25 = icmp slt i32 %23, %24
br i1 %25, label %26, label %190
26: ; preds = %22
store i32 0, i32* %18
br label %27
27: ; preds = %183, %26
%28 = load i32, i32* %18
%29 = load i32, i32* %11
%30 = icmp slt i32 %28, %29
br i1 %30, label %31, label %186
31: ; preds = %27
%32 = load double**, double*** %12
%33 = load i32, i32* %17
%34 = sext i32 %33 to i64
%35 = getelementptr inbounds double*, double** %32, i64 %34
%36 = load double*, double** %35
%37 = load i32, i32* %18
%38 = sext i32 %37 to i64
%39 = getelementptr inbounds double, double* %36, i64 %38
%40 = load double, double* %39
%41 = call i32 @isnan(double %40)
%42 = icmp ne i32 %41, 0
br i1 %42, label %182, label %43
43: ; preds = %31
%44 = load double**, double*** %12
%45 = load i32, i32* %17
%46 = sext i32 %45 to i64
%47 = getelementptr inbounds double*, double** %44, i64 %46
%48 = load double*, double** %47
%49 = load i32, i32* %18
%50 = sext i32 %49 to i64
%51 = getelementptr inbounds double, double* %48, i64 %50
%52 = load double, double* %51
%53 = call i32 @isinf(double %52)
%54 = icmp ne i32 %53, 0
br i1 %54, label %182, label %55
55: ; preds = %43
%56 = load double**, double*** %12
%57 = load i32, i32* %17
%58 = sext i32 %57 to i64
%59 = getelementptr inbounds double*, double** %56, i64 %58
%60 = load double*, double** %59
%61 = load i32, i32* %18
%62 = sext i32 %61 to i64
%63 = getelementptr inbounds double, double* %60, i64 %62
%64 = load double, double* %63
%65 = load double, double* %16
%66 = fcmp ogt double %64, %65
br i1 %66, label %67, label %182
67: ; preds = %55
%68 = load double**, double*** %12
%69 = load i32, i32* %17
%70 = sext i32 %69 to i64
%71 = getelementptr inbounds double*, double** %68, i64 %70
%72 = load double*, double** %71
%73 = load i32, i32* %18
%74 = sext i32 %73 to i64
%75 = getelementptr inbounds double, double* %72, i64 %74
%76 = load double, double* %75
%77 = load double*, double** %14
%78 = load i32, i32* %17
%79 = sext i32 %78 to i64
%80 = getelementptr inbounds double, double* %77, i64 %79
%81 = load double, double* %80
%82 = fsub double %76, %81
%83 = fcmp ogt double %82, 0.000000e+00
br i1 %83, label %84, label %100
84: ; preds = %67
%85 = load double**, double*** %12
%86 = load i32, i32* %17
%87 = sext i32 %86 to i64
%88 = getelementptr inbounds double*, double** %85, i64 %87
%89 = load double*, double** %88
%90 = load i32, i32* %18
%91 = sext i32 %90 to i64
%92 = getelementptr inbounds double, double* %89, i64 %91
%93 = load double, double* %92
%94 = load double*, double** %14
%95 = load i32, i32* %17
%96 = sext i32 %95 to i64
%97 = getelementptr inbounds double, double* %94, i64 %96
%98 = load double, double* %97
%99 = fsub double %93, %98
br label %117
100: ; preds = %67
%101 = load double**, double*** %12
%102 = load i32, i32* %17
%103 = sext i32 %102 to i64
%104 = getelementptr inbounds double*, double** %101, i64 %103
%105 = load double*, double** %104
%106 = load i32, i32* %18
%107 = sext i32 %106 to i64
%108 = getelementptr inbounds double, double* %105, i64 %107
%109 = load double, double* %108
%110 = load double*, double** %14
%111 = load i32, i32* %17
%112 = sext i32 %111 to i64
%113 = getelementptr inbounds double, double* %110, i64 %112
%114 = load double, double* %113
%115 = fsub double %109, %114
%116 = fneg double %115
br label %117
117: ; preds = %100, %84
%118 = phi double [ %99, %84 ], [ %116, %100 ]
%119 = load double**, double*** %13
%120 = load i32, i32* %17
%121 = sext i32 %120 to i64
%122 = getelementptr inbounds double*, double** %119, i64 %121
%123 = load double*, double** %122
%124 = load i32, i32* %18
%125 = sext i32 %124 to i64
%126 = getelementptr inbounds double, double* %123, i64 %125
%127 = load double, double* %126
%128 = fdiv double %118, %127
%129 = load double, double* %15
%130 = fcmp olt double %128, %129
br i1 %130, label %134, label %131
131: ; preds = %117
%132 = load double, double* %15
%133 = fcmp olt double %132, 0.000000e+00
br i1 %133, label %134, label %182
134: ; preds = %131, %117
%135 = load double**, double*** %12
%136 = load i32, i32* %17
%137 = sext i32 %136 to i64
%138 = getelementptr inbounds double*, double** %135, i64 %137
%139 = load double*, double** %138
%140 = load i32, i32* %18
%141 = sext i32 %140 to i64
%142 = getelementptr inbounds double, double* %139, i64 %141
%143 = load double, double* %142
%144 = load double*, double** %14
%145 = load i32, i32* %17
%146 = sext i32 %145 to i64
%147 = getelementptr inbounds double, double* %144, i64 %146
%148 = load double, double* %147
%149 = fsub double %143, %148
%150 = load double*, double** %9
%151 = load i32, i32* %17
%152 = add nsw i32 %151, 1
%153 = sext i32 %152 to i64
%154 = getelementptr inbounds double, double* %150, i64 %153
%155 = load double, double* %154
%156 = load double*, double** %9
%157 = load i32, i32* %10
%158 = load i32, i32* %18
%159 = add nsw i32 %157, %158
%160 = add nsw i32 %159, 1
%161 = sext i32 %160 to i64
%162 = getelementptr inbounds double, double* %156, i64 %161
%163 = load double, double* %162
%164 = fneg double %155
%165 = call double @llvm.fmuladd.f64(double %164, double %163, double %149)
%166 = fpext double %165 to x86_fp80
%167 = load double**, double*** %13
%168 = load i32, i32* %17
%169 = sext i32 %168 to i64
%170 = getelementptr inbounds double*, double** %167, i64 %169
%171 = load double*, double** %170
%172 = load i32, i32* %18
%173 = sext i32 %172 to i64
%174 = getelementptr inbounds double, double* %171, i64 %173
%175 = load double, double* %174
%176 = fpext double %175 to x86_fp80
%177 = fdiv x86_fp80 %166, %176
store x86_fp80 %177, x86_fp80* %21
%178 = load x86_fp80, x86_fp80* %21
%179 = load x86_fp80, x86_fp80* %21
%180 = load x86_fp80, x86_fp80* %20
%181 = call x86_fp80 @llvm.fmuladd.f80(x86_fp80 %178, x86_fp80 %179, x86_fp80 %180)
store x86_fp80 %181, x86_fp80* %20
br label %182
182: ; preds = %134, %131, %55, %43, %31
br label %183
183: ; preds = %182
%184 = load i32, i32* %18
%185 = add nsw i32 %184, 1
store i32 %185, i32* %18
br label %27
186: ; preds = %27
br label %187
187: ; preds = %186
%188 = load i32, i32* %17
%189 = add nsw i32 %188, 1
store i32 %189, i32* %17
br label %22
190: ; preds = %22
%191 = load x86_fp80, x86_fp80* %20
%192 = load i32, i32* %11
%193 = load i32, i32* %10
%194 = mul nsw i32 %192, %193
%195 = sitofp i32 %194 to double
%196 = fpext double %195 to x86_fp80
%197 = fdiv x86_fp80 %191, %196
%198 = fptrunc x86_fp80 %197 to double
ret double %198
}

LLVM IR (-O2):
define double @chisqr_sysrem(double* nocapture %0, i32 %1, i32 %2, double** nocapture %3, double** nocapture %4, double* nocapture %5, double %6, double %7) {
%9 = icmp sgt i32 %1, 0
br i1 %9, label %10, label %79
10: ; preds = %8
%11 = icmp sgt i32 %2, 0
%12 = fcmp olt double %6, 0.000000e+00
%13 = zext i32 %1 to i64
%14 = zext i32 %2 to i64
br label %15
15: ; preds = %10, %75
%16 = phi i64 [ 0, %10 ], [ %76, %75 ]
%17 = phi x86_fp80 [ 0xK00000000000000000000, %10 ], [ %77, %75 ]
br i1 %11, label %20, label %18
18: ; preds = %15
%19 = add nuw nsw i64 %16, 1
br label %75
20: ; preds = %15
%21 = getelementptr inbounds double*, double** %3, i64 %16
%22 = getelementptr inbounds double, double* %5, i64 %16
%23 = getelementptr inbounds double*, double** %4, i64 %16
%24 = add nuw nsw i64 %16, 1
%25 = getelementptr inbounds double, double* %0, i64 %24
br label %26
26: ; preds = %20, %71
%27 = phi i64 [ 0, %20 ], [ %73, %71 ]
%28 = phi x86_fp80 [ %17, %20 ], [ %72, %71 ]
%29 = load double*, double** %21
%30 = getelementptr inbounds double, double* %29, i64 %27
%31 = load double, double* %30
%32 = tail call i32 @isnan(double %31)
%33 = icmp eq i32 %32, 0
br i1 %33, label %34, label %71
34: ; preds = %26
%35 = load double*, double** %21
%36 = getelementptr inbounds double, double* %35, i64 %27
%37 = load double, double* %36
%38 = tail call i32 @isinf(double %37)
%39 = icmp eq i32 %38, 0
br i1 %39, label %40, label %71
40: ; preds = %34
%41 = load double*, double** %21
%42 = getelementptr inbounds double, double* %41, i64 %27
%43 = load double, double* %42
%44 = fcmp ogt double %43, %7
br i1 %44, label %45, label %71
45: ; preds = %40
%46 = load double, double* %22
%47 = fsub double %43, %46
%48 = fcmp ogt double %47, 0.000000e+00
%49 = fneg double %47
%50 = select i1 %48, double %47, double %49
%51 = load double*, double** %23
%52 = getelementptr inbounds double, double* %51, i64 %27
%53 = load double, double* %52
%54 = fdiv double %50, %53
%55 = fcmp olt double %54, %6
%56 = or i1 %12, %55
br i1 %56, label %57, label %71
57: ; preds = %45
%58 = load double, double* %25
%59 = trunc i64 %27 to i32
%60 = add i32 %59, 1
%61 = add i32 %60, %1
%62 = sext i32 %61 to i64
%63 = getelementptr inbounds double, double* %0, i64 %62
%64 = load double, double* %63
%65 = fneg double %58
%66 = tail call double @llvm.fmuladd.f64(double %65, double %64, double %47)
%67 = fpext double %66 to x86_fp80
%68 = fpext double %53 to x86_fp80
%69 = fdiv x86_fp80 %67, %68
%70 = tail call x86_fp80 @llvm.fmuladd.f80(x86_fp80 %69, x86_fp80 %69, x86_fp80 %28)
br label %71
71: ; preds = %26, %34, %40, %57, %45
%72 = phi x86_fp80 [ %28, %26 ], [ %28, %34 ], [ %70, %57 ], [ %28, %45 ], [ %28, %40 ]
%73 = add nuw nsw i64 %27, 1
%74 = icmp eq i64 %73, %14
br i1 %74, label %75, label %26
75: ; preds = %71, %18
%76 = phi i64 [ %19, %18 ], [ %24, %71 ]
%77 = phi x86_fp80 [ %17, %18 ], [ %72, %71 ]
%78 = icmp eq i64 %76, %13
br i1 %78, label %79, label %15
79: ; preds = %75, %8
%80 = phi x86_fp80 [ 0xK00000000000000000000, %8 ], [ %77, %75 ]
%81 = mul nsw i32 %2, %1
%82 = sitofp i32 %81 to x86_fp80
%83 = fdiv x86_fp80 %80, %82
%84 = fptrunc x86_fp80 %83 to double
ret double %84
}