void decorr_star_sysrem(double *star_coeff, double **mag, double **sig,
                        double *mag_ave, int numlc, int numJD, int **clip,
                        double sigma_clip2, double *image_coeff,
                        double saturation) {
  int j, k, n;
  long double val1, val2;
  for (j = 0; j < numlc; j++) {
    n = 0;
    val1 = 0.0;
    val2 = 0.0;
    for (k = 0; k < numJD; k++) {
      if (clip[j][k]) {
        val1 += (long double)(mag[j][k] - mag_ave[j]) * image_coeff[k] /
                (sig[j][k] * sig[j][k]);
        val2 += (long double)(image_coeff[k] * image_coeff[k]) /
                (sig[j][k] * sig[j][k]);
        n++;
      }
    }
    if (n > 0 && val2 > 0.0)
      star_coeff[j] = (double)val1 / val2;
    else
      star_coeff[j] = 0.0;
  }
}


x86 assembly (clang, -O0):
<decorr_star_sysrem>:
push   %rbp
mov    %rsp,%rbp
sub    $0x30,%rsp
mov    0x18(%rbp),%rax
mov    0x10(%rbp),%rax
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %rcx,-0x20(%rbp)
mov    %r8d,-0x24(%rbp)
mov    %r9d,-0x28(%rbp)
movsd  %xmm0,-0x30(%rbp)
movsd  %xmm1,-0x38(%rbp)
movl   $0x0,-0x3c(%rbp)
mov    -0x3c(%rbp),%eax
cmp    -0x24(%rbp),%eax
jge    12c8 <decorr_star_sysrem+0x1c8>
movl   $0x0,-0x44(%rbp)
fldz
fld    %st(0)
fstpt  -0x60(%rbp)
fstpt  -0x70(%rbp)
movl   $0x0,-0x40(%rbp)
mov    -0x40(%rbp),%eax
cmp    -0x28(%rbp),%eax
jge    124d <decorr_star_sysrem+0x14d>
mov    0x10(%rbp),%rax
movslq -0x3c(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x40(%rbp),%rcx
cmpl   $0x0,(%rax,%rcx,4)
je     123a <decorr_star_sysrem+0x13a>
mov    -0x10(%rbp),%rax
movslq -0x3c(%rbp),%rdx
mov    (%rax,%rdx,8),%rax
movslq -0x40(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x20(%rbp),%rax
movsd  (%rax,%rdx,8),%xmm1
subsd  %xmm1,%xmm0
movsd  %xmm0,-0x98(%rbp)
ff
fldl   -0x98(%rbp)
mov    0x18(%rbp),%rax
fldl   (%rax,%rcx,8)
fmulp  %st,%st(1)
mov    -0x18(%rbp),%rax
mov    (%rax,%rdx,8),%rax
movsd  (%rax,%rcx,8),%xmm0
mulsd  %xmm0,%xmm0
movsd  %xmm0,-0x90(%rbp)
ff
fldl   -0x90(%rbp)
fdivrp %st,%st(1)
fldt   -0x60(%rbp)
faddp  %st,%st(1)
fstpt  -0x60(%rbp)
mov    0x18(%rbp),%rax
movslq -0x40(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mulsd  %xmm0,%xmm0
movsd  %xmm0,-0xa0(%rbp)
ff
fldl   -0xa0(%rbp)
mov    -0x18(%rbp),%rax
movslq -0x3c(%rbp),%rdx
mov    (%rax,%rdx,8),%rax
movsd  (%rax,%rcx,8),%xmm0
mulsd  %xmm0,%xmm0
movsd  %xmm0,-0xa8(%rbp)
ff
fldl   -0xa8(%rbp)
fdivrp %st,%st(1)
fldt   -0x70(%rbp)
faddp  %st,%st(1)
fstpt  -0x70(%rbp)
mov    -0x44(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x44(%rbp)
jmp    123f <decorr_star_sysrem+0x13f>
mov    -0x40(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x40(%rbp)
jmp    115d <decorr_star_sysrem+0x5d>
cmpl   $0x0,-0x44(%rbp)
jle    12a5 <decorr_star_sysrem+0x1a5>
fldt   -0x70(%rbp)
fldz
fxch   %st(1)
fucomip %st(1),%st
fstp   %st(0)
jbe    12a5 <decorr_star_sysrem+0x1a5>
jmp    126d <decorr_star_sysrem+0x16d>
fldt   -0x60(%rbp)
fstpl  -0x78(%rbp)
movsd  -0x78(%rbp),%xmm0
movsd  %xmm0,-0x80(%rbp)
fldl   -0x80(%rbp)
fldt   -0x70(%rbp)
fdivrp %st,%st(1)
fstpl  -0x88(%rbp)
movsd  -0x88(%rbp),%xmm0
ff
mov    -0x8(%rbp),%rax
movslq -0x3c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
jmp    12b5 <decorr_star_sysrem+0x1b5>
mov    -0x8(%rbp),%rax
movslq -0x3c(%rbp),%rcx
xorps  %xmm0,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
jmp    12ba <decorr_star_sysrem+0x1ba>
mov    -0x3c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x3c(%rbp)
jmp    1139 <decorr_star_sysrem+0x39>
add    $0x30,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<decorr_star_sysrem>:
push   %rbp
push   %r15
push   %r14
push   %r12
push   %rbx
test   %r8d,%r8d
fldz
jle    1218 <decorr_star_sysrem+0x118>
fstp   %st(0)
mov    0x38(%rsp),%r11
mov    0x30(%rsp),%r10
mov    %r8d,%r8d
mov    %r9d,%r14d
xor    %ebx,%ebx
fldz
jmp    114c <decorr_star_sysrem+0x4c>
fstp   %st(1)
fstp   %st(0)
nop
fldz
fldz
fxch   %st(1)
fstp   %st(0)
fstp   %st(0)
movsd  %xmm0,(%rdi,%rbx,8)
add    $0x1,%rbx
cmp    %r8,%rbx
je     1218 <decorr_star_sysrem+0x118>
xorps  %xmm0,%xmm0
test   %r9d,%r9d
jle    1130 <decorr_star_sysrem+0x30>
mov    (%r10,%rbx,8),%r15
xor    %eax,%eax
fld    %st(0)
fld    %st(1)
xor    %r12d,%r12d
jmp    1179 <decorr_star_sysrem+0x79>
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
add    $0x1,%rax
cmp    %rax,%r14
je     11e0 <decorr_star_sysrem+0xe0>
cmpl   $0x0,(%r15,%rax,4)
je     1170 <decorr_star_sysrem+0x70>
mov    (%rsi,%rbx,8),%rbp
movsd  0x0(%rbp,%rax,8),%xmm1
subsd  (%rcx,%rbx,8),%xmm1
movsd  %xmm1,-0x18(%rsp)
fldl   -0x18(%rsp)
movsd  (%r11,%rax,8),%xmm1
movsd  %xmm1,-0x8(%rsp)
fmull  -0x8(%rsp)
mov    (%rdx,%rbx,8),%rbp
movsd  0x0(%rbp,%rax,8),%xmm2
mulsd  %xmm2,%xmm2
movsd  %xmm2,-0x20(%rsp)
fldl   -0x20(%rsp)
fdivr  %st,%st(1)
fxch   %st(2)
faddp  %st,%st(1)
mulsd  %xmm1,%xmm1
movsd  %xmm1,-0x10(%rsp)
fxch   %st(1)
fdivrl -0x10(%rsp)
faddp  %st,%st(2)
add    $0x1,%r12d
jmp    1170 <decorr_star_sysrem+0x70>
nop
test   %r12d,%r12d
jle    112b <decorr_star_sysrem+0x2b>
fxch   %st(1)
fucomi %st(2),%st
jbe    1136 <decorr_star_sysrem+0x36>
fxch   %st(1)
fstpl  -0x28(%rsp)
movsd  -0x28(%rsp),%xmm0
movsd  %xmm0,-0x30(%rsp)
fdivrl -0x30(%rsp)
fstpl  -0x38(%rsp)
movsd  -0x38(%rsp),%xmm0
jmp    1130 <decorr_star_sysrem+0x30>
fstp   %st(0)
pop    %rbx
pop    %r12
pop    %r14
pop    %r15
pop    %rbp
ret

LLVM IR (-O0):
define void @decorr_star_sysrem(double* %0, double** %1, double** %2, double* %3, i32 %4, i32 %5, i32** %6, double %7, double* %8, double %9) {
%11 = alloca double*
%12 = alloca double**
%13 = alloca double**
%14 = alloca double*
%15 = alloca i32
%16 = alloca i32
%17 = alloca i32**
%18 = alloca double
%19 = alloca double*
%20 = alloca double
%21 = alloca i32
%22 = alloca i32
%23 = alloca i32
%24 = alloca x86_fp80
%25 = alloca x86_fp80
store double* %0, double** %11
store double** %1, double*** %12
store double** %2, double*** %13
store double* %3, double** %14
store i32 %4, i32* %15
store i32 %5, i32* %16
store i32** %6, i32*** %17
store double %7, double* %18
store double* %8, double** %19
store double %9, double* %20
store i32 0, i32* %21
br label %26
26: ; preds = %157, %10
%27 = load i32, i32* %21
%28 = load i32, i32* %15
%29 = icmp slt i32 %27, %28
br i1 %29, label %30, label %160
30: ; preds = %26
store i32 0, i32* %23
store x86_fp80 0xK00000000000000000000, x86_fp80* %24
store x86_fp80 0xK00000000000000000000, x86_fp80* %25
store i32 0, i32* %22
br label %31
31: ; preds = %131, %30
%32 = load i32, i32* %22
%33 = load i32, i32* %16
%34 = icmp slt i32 %32, %33
br i1 %34, label %35, label %134
35: ; preds = %31
%36 = load i32**, i32*** %17
%37 = load i32, i32* %21
%38 = sext i32 %37 to i64
%39 = getelementptr inbounds i32*, i32** %36, i64 %38
%40 = load i32*, i32** %39
%41 = load i32, i32* %22
%42 = sext i32 %41 to i64
%43 = getelementptr inbounds i32, i32* %40, i64 %42
%44 = load i32, i32* %43
%45 = icmp ne i32 %44, 0
br i1 %45, label %46, label %130
46: ; preds = %35
%47 = load double**, double*** %12
%48 = load i32, i32* %21
%49 = sext i32 %48 to i64
%50 = getelementptr inbounds double*, double** %47, i64 %49
%51 = load double*, double** %50
%52 = load i32, i32* %22
%53 = sext i32 %52 to i64
%54 = getelementptr inbounds double, double* %51, i64 %53
%55 = load double, double* %54
%56 = load double*, double** %14
%57 = load i32, i32* %21
%58 = sext i32 %57 to i64
%59 = getelementptr inbounds double, double* %56, i64 %58
%60 = load double, double* %59
%61 = fsub double %55, %60
%62 = fpext double %61 to x86_fp80
%63 = load double*, double** %19
%64 = load i32, i32* %22
%65 = sext i32 %64 to i64
%66 = getelementptr inbounds double, double* %63, i64 %65
%67 = load double, double* %66
%68 = fpext double %67 to x86_fp80
%69 = fmul x86_fp80 %62, %68
%70 = load double**, double*** %13
%71 = load i32, i32* %21
%72 = sext i32 %71 to i64
%73 = getelementptr inbounds double*, double** %70, i64 %72
%74 = load double*, double** %73
%75 = load i32, i32* %22
%76 = sext i32 %75 to i64
%77 = getelementptr inbounds double, double* %74, i64 %76
%78 = load double, double* %77
%79 = load double**, double*** %13
%80 = load i32, i32* %21
%81 = sext i32 %80 to i64
%82 = getelementptr inbounds double*, double** %79, i64 %81
%83 = load double*, double** %82
%84 = load i32, i32* %22
%85 = sext i32 %84 to i64
%86 = getelementptr inbounds double, double* %83, i64 %85
%87 = load double, double* %86
%88 = fmul double %78, %87
%89 = fpext double %88 to x86_fp80
%90 = fdiv x86_fp80 %69, %89
%91 = load x86_fp80, x86_fp80* %24
%92 = fadd x86_fp80 %91, %90
store x86_fp80 %92, x86_fp80* %24
%93 = load double*, double** %19
%94 = load i32, i32* %22
%95 = sext i32 %94 to i64
%96 = getelementptr inbounds double, double* %93, i64 %95
%97 = load double, double* %96
%98 = load double*, double** %19
%99 = load i32, i32* %22
%100 = sext i32 %99 to i64
%101 = getelementptr inbounds double, double* %98, i64 %100
%102 = load double, double* %101
%103 = fmul double %97, %102
%104 = fpext double %103 to x86_fp80
%105 = load double**, double*** %13
%106 = load i32, i32* %21
%107 = sext i32 %106 to i64
%108 = getelementptr inbounds double*, double** %105, i64 %107
%109 = load double*, double** %108
%110 = load i32, i32* %22
%111 = sext i32 %110 to i64
%112 = getelementptr inbounds double, double* %109, i64 %111
%113 = load double, double* %112
%114 = load double**, double*** %13
%115 = load i32, i32* %21
%116 = sext i32 %115 to i64
%117 = getelementptr inbounds double*, double** %114, i64 %116
%118 = load double*, double** %117
%119 = load i32, i32* %22
%120 = sext i32 %119 to i64
%121 = getelementptr inbounds double, double* %118, i64 %120
%122 = load double, double* %121
%123 = fmul double %113, %122
%124 = fpext double %123 to x86_fp80
%125 = fdiv x86_fp80 %104, %124
%126 = load x86_fp80, x86_fp80* %25
%127 = fadd x86_fp80 %126, %125
store x86_fp80 %127, x86_fp80* %25
%128 = load i32, i32* %23
%129 = add nsw i32 %128, 1
store i32 %129, i32* %23
br label %130
130: ; preds = %46, %35
br label %131
131: ; preds = %130
%132 = load i32, i32* %22
%133 = add nsw i32 %132, 1
store i32 %133, i32* %22
br label %31
134: ; preds = %31
%135 = load i32, i32* %23
%136 = icmp sgt i32 %135, 0
br i1 %136, label %137, label %151
137: ; preds = %134
%138 = load x86_fp80, x86_fp80* %25
%139 = fcmp ogt x86_fp80 %138, 0xK00000000000000000000
br i1 %139, label %140, label %151
140: ; preds = %137
%141 = load x86_fp80, x86_fp80* %24
%142 = fptrunc x86_fp80 %141 to double
%143 = fpext double %142 to x86_fp80
%144 = load x86_fp80, x86_fp80* %25
%145 = fdiv x86_fp80 %143, %144
%146 = fptrunc x86_fp80 %145 to double
%147 = load double*, double** %11
%148 = load i32, i32* %21
%149 = sext i32 %148 to i64
%150 = getelementptr inbounds double, double* %147, i64 %149
store double %146, double* %150
br label %156
151: ; preds = %137, %134
%152 = load double*, double** %11
%153 = load i32, i32* %21
%154 = sext i32 %153 to i64
%155 = getelementptr inbounds double, double* %152, i64 %154
store double 0.000000e+00, double* %155
br label %156
156: ; preds = %151, %140
br label %157
157: ; preds = %156
%158 = load i32, i32* %21
%159 = add nsw i32 %158, 1
store i32 %159, i32* %21
br label %26
160: ; preds = %26
ret void
}

LLVM IR (-O2):
define void @decorr_star_sysrem(double* nocapture writeonly %0, double** nocapture %1, double** nocapture %2, double* nocapture %3, i32 %4, i32 %5, i32** nocapture %6, double %7, double* nocapture %8, double %9) {
%11 = icmp sgt i32 %4, 0
br i1 %11, label %12, label %75
12: ; preds = %10
%13 = icmp sgt i32 %5, 0
%14 = zext i32 %4 to i64
%15 = zext i32 %5 to i64
br label %16
16: ; preds = %12, %70
%17 = phi i64 [ 0, %12 ], [ %73, %70 ]
br i1 %13, label %18, label %70
18: ; preds = %16
%19 = getelementptr inbounds i32*, i32** %6, i64 %17
%20 = load i32*, i32** %19
%21 = getelementptr inbounds double*, double** %1, i64 %17
%22 = getelementptr inbounds double, double* %3, i64 %17
%23 = getelementptr inbounds double*, double** %2, i64 %17
br label %24
24: ; preds = %18, %55
%25 = phi i64 [ 0, %18 ], [ %59, %55 ]
%26 = phi x86_fp80 [ 0xK00000000000000000000, %18 ], [ %58, %55 ]
%27 = phi x86_fp80 [ 0xK00000000000000000000, %18 ], [ %57, %55 ]
%28 = phi i32 [ 0, %18 ], [ %56, %55 ]
%29 = getelementptr inbounds i32, i32* %20, i64 %25
%30 = load i32, i32* %29
%31 = icmp eq i32 %30, 0
br i1 %31, label %55, label %32
32: ; preds = %24
%33 = load double*, double** %21
%34 = getelementptr inbounds double, double* %33, i64 %25
%35 = load double, double* %34
%36 = load double, double* %22
%37 = fsub double %35, %36
%38 = fpext double %37 to x86_fp80
%39 = getelementptr inbounds double, double* %8, i64 %25
%40 = load double, double* %39
%41 = fpext double %40 to x86_fp80
%42 = fmul x86_fp80 %38, %41
%43 = load double*, double** %23
%44 = getelementptr inbounds double, double* %43, i64 %25
%45 = load double, double* %44
%46 = fmul double %45, %45
%47 = fpext double %46 to x86_fp80
%48 = fdiv x86_fp80 %42, %47
%49 = fadd x86_fp80 %27, %48
%50 = fmul double %40, %40
%51 = fpext double %50 to x86_fp80
%52 = fdiv x86_fp80 %51, %47
%53 = fadd x86_fp80 %26, %52
%54 = add nsw i32 %28, 1
br label %55
55: ; preds = %24, %32
%56 = phi i32 [ %54, %32 ], [ %28, %24 ]
%57 = phi x86_fp80 [ %49, %32 ], [ %27, %24 ]
%58 = phi x86_fp80 [ %53, %32 ], [ %26, %24 ]
%59 = add nuw nsw i64 %25, 1
%60 = icmp eq i64 %59, %15
br i1 %60, label %61, label %24
61: ; preds = %55
%62 = icmp sgt i32 %56, 0
%63 = fcmp ogt x86_fp80 %58, 0xK00000000000000000000
%64 = select i1 %62, i1 %63, i1 false
br i1 %64, label %65, label %70
65: ; preds = %61
%66 = fptrunc x86_fp80 %57 to double
%67 = fpext double %66 to x86_fp80
%68 = fdiv x86_fp80 %67, %58
%69 = fptrunc x86_fp80 %68 to double
br label %70
70: ; preds = %16, %61, %65
%71 = phi double [ %69, %65 ], [ 0.000000e+00, %61 ], [ 0.000000e+00, %16 ]
%72 = getelementptr inbounds double, double* %0, i64 %17
store double %71, double* %72
%73 = add nuw nsw i64 %17, 1
%74 = icmp eq i64 %73, %14
br i1 %74, label %75, label %16
75: ; preds = %70, %10
ret void
}