double itersolve_chisqr_sysrem(double *star_coeff, double *image_coeff,
                               int numlc, int numJD, double **mag, double **sig,
                               double *mag_ave, double sigma_clip,
                               double upperdiff, double saturation,
                               int **clip) {
  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 (clip[j][k]) {
        dum = (long double)(mag[j][k] - mag_ave[j] -
                            star_coeff[j] * image_coeff[k]) /
              sig[j][k];
        S += dum * dum;
        n++;
      }
  return ((double)(S / (double)n));
}


x86 assembly (clang, -O0):
<itersolve_chisqr_sysrem>:
push   %rbp
mov    %rsp,%rbp
sub    $0x10,%rsp
mov    0x18(%rbp),%rax
mov    0x10(%rbp),%rax
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %edx,-0x14(%rbp)
mov    %ecx,-0x18(%rbp)
mov    %r8,-0x20(%rbp)
mov    %r9,-0x28(%rbp)
movsd  %xmm0,-0x30(%rbp)
movsd  %xmm1,-0x38(%rbp)
movsd  %xmm2,-0x40(%rbp)
fldz
fstpt  -0x60(%rbp)
movl   $0x0,-0x4c(%rbp)
movl   $0x0,-0x44(%rbp)
mov    -0x44(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    1219 <itersolve_chisqr_sysrem+0x119>
movl   $0x0,-0x48(%rbp)
mov    -0x48(%rbp),%eax
cmp    -0x18(%rbp),%eax
jge    1206 <itersolve_chisqr_sysrem+0x106>
mov    0x18(%rbp),%rax
movslq -0x44(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x48(%rbp),%rcx
cmpl   $0x0,(%rax,%rcx,4)
je     11f3 <itersolve_chisqr_sysrem+0xf3>
mov    -0x20(%rbp),%rax
movslq -0x44(%rbp),%rdx
mov    (%rax,%rdx,8),%rax
movslq -0x48(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    0x10(%rbp),%rax
movsd  (%rax,%rdx,8),%xmm1
subsd  %xmm1,%xmm0
mov    -0x8(%rbp),%rax
movsd  (%rax,%rdx,8),%xmm1
mov    -0x10(%rbp),%rax
movsd  (%rax,%rcx,8),%xmm2
mulsd  %xmm2,%xmm1
subsd  %xmm1,%xmm0
movsd  %xmm0,-0x88(%rbp)
ff
fldl   -0x88(%rbp)
mov    -0x28(%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)
mov    -0x4c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x4c(%rbp)
jmp    11f8 <itersolve_chisqr_sysrem+0xf8>
mov    -0x48(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x48(%rbp)
jmp    115b <itersolve_chisqr_sysrem+0x5b>
jmp    120b <itersolve_chisqr_sysrem+0x10b>
mov    -0x44(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x44(%rbp)
jmp    1148 <itersolve_chisqr_sysrem+0x48>
fldt   -0x60(%rbp)
cvtsi2sdl -0x4c(%rbp),%xmm0
movsd  %xmm0,-0x78(%rbp)
fldl   -0x78(%rbp)
fdivrp %st,%st(1)
fstpl  -0x80(%rbp)
movsd  -0x80(%rbp),%xmm0
add    $0x10,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<itersolve_chisqr_sysrem>:
push   %rbp
push   %r15
push   %r14
push   %r12
push   %rbx
fldz
test   %edx,%edx
jle    1199 <itersolve_chisqr_sysrem+0x99>
fstp   %st(0)
mov    0x38(%rsp),%r10
mov    0x30(%rsp),%r14
mov    %edx,%r11d
mov    %ecx,%r15d
xor    %eax,%eax
fldz
xor    %r12d,%r12d
jmp    1139 <itersolve_chisqr_sysrem+0x39>
nopl   (%rax)
add    $0x1,%rax
cmp    %r11,%rax
je     119c <itersolve_chisqr_sysrem+0x9c>
test   %ecx,%ecx
jle    1130 <itersolve_chisqr_sysrem+0x30>
mov    (%r10,%rax,8),%rdx
xor    %ebx,%ebx
jmp    1159 <itersolve_chisqr_sysrem+0x59>
cs nopw 0x0(%rax,%rax,1)
nop
add    $0x1,%rbx
cmp    %rbx,%r15
je     1130 <itersolve_chisqr_sysrem+0x30>
cmpl   $0x0,(%rdx,%rbx,4)
je     1150 <itersolve_chisqr_sysrem+0x50>
mov    (%r8,%rax,8),%rbp
movsd  0x0(%rbp,%rbx,8),%xmm0
subsd  (%r14,%rax,8),%xmm0
movsd  (%rdi,%rax,8),%xmm1
mulsd  (%rsi,%rbx,8),%xmm1
subsd  %xmm1,%xmm0
movsd  %xmm0,-0x8(%rsp)
fldl   -0x8(%rsp)
mov    (%r9,%rax,8),%rbp
fdivl  0x0(%rbp,%rbx,8)
fmul   %st,%st(0)
faddp  %st,%st(1)
add    $0x1,%r12d
jmp    1150 <itersolve_chisqr_sysrem+0x50>
xor    %r12d,%r12d
mov    %r12d,-0x14(%rsp)
fidivl -0x14(%rsp)
fstpl  -0x10(%rsp)
movsd  -0x10(%rsp),%xmm0
pop    %rbx
pop    %r12
pop    %r14
pop    %r15
pop    %rbp
ret

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

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