void opt_param_update(int *PassiveSet, int *NotGood, int *NonOptSet,
                      int *InfeaSet, double *X, double *Y, int *NotOptCols_nb,
                      int N, int K) {
  int i, k, prevNotGood;
  for (i = 0; i < N; i++) {
    prevNotGood = NotGood[i];
    NotGood[i] = 0;
    for (k = 0; k < K; k++) {
      if (prevNotGood && Y[k * N + i] < 0 && !PassiveSet[k * N + i]) {
        NonOptSet[k * N + i] = 1;
        NotGood[i]++;
      } else {
        NonOptSet[k * N + i] = 0;
      }
      if (prevNotGood && X[k * N + i] < 0 && PassiveSet[k * N + i]) {
        InfeaSet[k * N + i] = 1;
        NotGood[i]++;
      } else {
        InfeaSet[k * N + i] = 0;
      }
    }
  }
  *NotOptCols_nb = 0;
  for (i = 0; i < N; i++) {
    if (NotGood[i])
      (*NotOptCols_nb)++;
  }
}


x86 assembly (clang, -O0):
<opt_param_update>:
push   %rbp
mov    %rsp,%rbp
mov    0x20(%rbp),%eax
mov    0x18(%rbp),%eax
mov    0x10(%rbp),%rax
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %rcx,-0x20(%rbp)
mov    %r8,-0x28(%rbp)
mov    %r9,-0x30(%rbp)
movl   $0x0,-0x34(%rbp)
mov    -0x34(%rbp),%eax
cmp    0x18(%rbp),%eax
jge    12a3 <opt_param_update+0x1a3>
mov    -0x10(%rbp),%rax
movslq -0x34(%rbp),%rcx
mov    (%rax,%rcx,4),%eax
mov    %eax,-0x3c(%rbp)
mov    -0x10(%rbp),%rax
movslq -0x34(%rbp),%rcx
movl   $0x0,(%rax,%rcx,4)
movl   $0x0,-0x38(%rbp)
mov    -0x38(%rbp),%eax
cmp    0x20(%rbp),%eax
jge    1290 <opt_param_update+0x190>
cmpl   $0x0,-0x3c(%rbp)
je     11db <opt_param_update+0xdb>
mov    -0x30(%rbp),%rax
mov    -0x38(%rbp),%ecx
imul   0x18(%rbp),%ecx
add    -0x34(%rbp),%ecx
movslq %ecx,%rcx
xorps  %xmm0,%xmm0
ucomisd (%rax,%rcx,8),%xmm0
jbe    11db <opt_param_update+0xdb>
mov    -0x8(%rbp),%rax
mov    -0x38(%rbp),%ecx
imul   0x18(%rbp),%ecx
add    -0x34(%rbp),%ecx
movslq %ecx,%rcx
cmpl   $0x0,(%rax,%rcx,4)
jne    11db <opt_param_update+0xdb>
mov    -0x18(%rbp),%rax
mov    -0x38(%rbp),%ecx
imul   0x18(%rbp),%ecx
add    -0x34(%rbp),%ecx
movslq %ecx,%rcx
movl   $0x1,(%rax,%rcx,4)
mov    -0x10(%rbp),%rax
movslq -0x34(%rbp),%rcx
mov    (%rax,%rcx,4),%edx
add    $0x1,%edx
mov    %edx,(%rax,%rcx,4)
jmp    11f3 <opt_param_update+0xf3>
mov    -0x18(%rbp),%rax
mov    -0x38(%rbp),%ecx
imul   0x18(%rbp),%ecx
add    -0x34(%rbp),%ecx
movslq %ecx,%rcx
movl   $0x0,(%rax,%rcx,4)
cmpl   $0x0,-0x3c(%rbp)
je     1265 <opt_param_update+0x165>
mov    -0x28(%rbp),%rax
mov    -0x38(%rbp),%ecx
imul   0x18(%rbp),%ecx
add    -0x34(%rbp),%ecx
movslq %ecx,%rcx
xorps  %xmm0,%xmm0
ucomisd (%rax,%rcx,8),%xmm0
jbe    1265 <opt_param_update+0x165>
mov    -0x8(%rbp),%rax
mov    -0x38(%rbp),%ecx
imul   0x18(%rbp),%ecx
add    -0x34(%rbp),%ecx
movslq %ecx,%rcx
cmpl   $0x0,(%rax,%rcx,4)
je     1265 <opt_param_update+0x165>
mov    -0x20(%rbp),%rax
mov    -0x38(%rbp),%ecx
imul   0x18(%rbp),%ecx
add    -0x34(%rbp),%ecx
movslq %ecx,%rcx
movl   $0x1,(%rax,%rcx,4)
mov    -0x10(%rbp),%rax
movslq -0x34(%rbp),%rcx
mov    (%rax,%rcx,4),%edx
add    $0x1,%edx
mov    %edx,(%rax,%rcx,4)
jmp    127d <opt_param_update+0x17d>
mov    -0x20(%rbp),%rax
mov    -0x38(%rbp),%ecx
imul   0x18(%rbp),%ecx
add    -0x34(%rbp),%ecx
movslq %ecx,%rcx
movl   $0x0,(%rax,%rcx,4)
jmp    1282 <opt_param_update+0x182>
mov    -0x38(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x38(%rbp)
jmp    115d <opt_param_update+0x5d>
jmp    1295 <opt_param_update+0x195>
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    112d <opt_param_update+0x2d>
mov    0x10(%rbp),%rax
movl   $0x0,(%rax)
movl   $0x0,-0x34(%rbp)
mov    -0x34(%rbp),%eax
cmp    0x18(%rbp),%eax
jge    12f0 <opt_param_update+0x1f0>
mov    -0x10(%rbp),%rax
movslq -0x34(%rbp),%rcx
cmpl   $0x0,(%rax,%rcx,4)
je     12dd <opt_param_update+0x1dd>
mov    0x10(%rbp),%rax
mov    (%rax),%ecx
add    $0x1,%ecx
mov    %ecx,(%rax)
jmp    12e2 <opt_param_update+0x1e2>
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    12b4 <opt_param_update+0x1b4>
pop    %rbp
ret

x86 assembly (clang, -O2):
<opt_param_update>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
mov    0x40(%rsp),%eax
mov    0x38(%rsp),%r10
test   %eax,%eax
jle    1215 <opt_param_update+0x115>
mov    0x48(%rsp),%r14d
mov    %eax,%r12d
mov    %r14d,%r15d
lea    0x0(,%r12,4),%rbp
xor    %r13d,%r13d
xorpd  %xmm0,%xmm0
jmp    1161 <opt_param_update+0x61>
nopw   0x0(%rax,%rax,1)
add    $0x1,%r13
add    $0x4,%rcx
add    $0x4,%rdx
add    $0x4,%rdi
add    $0x8,%r9
add    $0x8,%r8
cmp    %r12,%r13
je     11e9 <opt_param_update+0xe9>
mov    (%rsi,%r13,4),%eax
movl   $0x0,(%rsi,%r13,4)
test   %r14d,%r14d
jle    1140 <opt_param_update+0x40>
mov    %r15,%rbx
xor    %r11d,%r11d
jmp    1199 <opt_param_update+0x99>
nopw   0x0(%rax,%rax,1)
movl   $0x0,(%rdx,%r11,1)
movl   $0x0,(%rcx,%r11,1)
add    %rbp,%r11
add    $0xffffffffffffffff,%rbx
je     1140 <opt_param_update+0x40>
test   %eax,%eax
je     1180 <opt_param_update+0x80>
ucomisd (%r9,%r11,2),%xmm0
jbe    11ac <opt_param_update+0xac>
cmpl   $0x0,(%rdi,%r11,1)
je     11d2 <opt_param_update+0xd2>
movl   $0x0,(%rdx,%r11,1)
ucomisd (%r8,%r11,2),%xmm0
jbe    1188 <opt_param_update+0x88>
cmpl   $0x0,(%rdi,%r11,1)
je     1188 <opt_param_update+0x88>
movl   $0x1,(%rcx,%r11,1)
addl   $0x1,(%rsi,%r13,4)
jmp    1190 <opt_param_update+0x90>
movl   $0x1,(%rdx,%r11,1)
addl   $0x1,(%rsi,%r13,4)
ucomisd (%r8,%r11,2),%xmm0
jbe    1188 <opt_param_update+0x88>
jmp    11bc <opt_param_update+0xbc>
movl   $0x0,(%r10)
mov    0x40(%rsp),%eax
test   %eax,%eax
jle    121c <opt_param_update+0x11c>
cmp    $0x1,%eax
jne    1227 <opt_param_update+0x127>
xor    %ecx,%ecx
xor    %eax,%eax
test   $0x1,%r12b
je     121c <opt_param_update+0x11c>
cmpl   $0x0,(%rsi,%rax,4)
je     121c <opt_param_update+0x11c>
add    $0x1,%ecx
mov    %ecx,(%r10)
jmp    121c <opt_param_update+0x11c>
movl   $0x0,(%r10)
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret
mov    %r12d,%edx
and    $0xfffffffe,%edx
xor    %ecx,%ecx
xor    %eax,%eax
jmp    1249 <opt_param_update+0x149>
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
add    $0x2,%rax
cmp    %rax,%rdx
je     1201 <opt_param_update+0x101>
cmpl   $0x0,(%rsi,%rax,4)
je     1255 <opt_param_update+0x155>
add    $0x1,%ecx
mov    %ecx,(%r10)
cmpl   $0x0,0x4(%rsi,%rax,4)
je     1240 <opt_param_update+0x140>
add    $0x1,%ecx
mov    %ecx,(%r10)
jmp    1240 <opt_param_update+0x140>

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

LLVM IR (-O2):
define void @opt_param_update(i32* nocapture %0, i32* nocapture %1, i32* nocapture writeonly %2, i32* nocapture writeonly %3, double* nocapture %4, double* nocapture %5, i32* nocapture %6, i32 %7, i32 %8) {
%10 = icmp sgt i32 %7, 0
br i1 %10, label %12, label %11
11: ; preds = %9
store i32 0, i32* %6
br label %101
12: ; preds = %9
%13 = icmp sgt i32 %8, 0
%14 = zext i32 %7 to i64
%15 = zext i32 %7 to i64
%16 = zext i32 %8 to i64
br label %17
17: ; preds = %12, %60
%18 = phi i64 [ 0, %12 ], [ %61, %60 ]
%19 = getelementptr inbounds i32, i32* %1, i64 %18
%20 = load i32, i32* %19
store i32 0, i32* %19
br i1 %13, label %21, label %60
21: ; preds = %17
%22 = icmp eq i32 %20, 0
br label %23
23: ; preds = %21, %57
%24 = phi i64 [ 0, %21 ], [ %58, %57 ]
%25 = mul nsw i64 %24, %14
%26 = add nuw nsw i64 %25, %18
br i1 %22, label %41, label %27
27: ; preds = %23
%28 = getelementptr inbounds double, double* %5, i64 %26
%29 = load double, double* %28
%30 = fcmp olt double %29, 0.000000e+00
br i1 %30, label %31, label %39
31: ; preds = %27
%32 = getelementptr inbounds i32, i32* %0, i64 %26
%33 = load i32, i32* %32
%34 = icmp eq i32 %33, 0
br i1 %34, label %35, label %39
35: ; preds = %31
%36 = getelementptr inbounds i32, i32* %2, i64 %26
store i32 1, i32* %36
%37 = load i32, i32* %19
%38 = add nsw i32 %37, 1
store i32 %38, i32* %19
br label %43
39: ; preds = %27, %31
%40 = getelementptr inbounds i32, i32* %2, i64 %26
store i32 0, i32* %40
br label %43
41: ; preds = %23
%42 = getelementptr inbounds i32, i32* %2, i64 %26
store i32 0, i32* %42
br label %55
43: ; preds = %39, %35
%44 = getelementptr inbounds double, double* %4, i64 %26
%45 = load double, double* %44
%46 = fcmp olt double %45, 0.000000e+00
br i1 %46, label %47, label %55
47: ; preds = %43
%48 = getelementptr inbounds i32, i32* %0, i64 %26
%49 = load i32, i32* %48
%50 = icmp eq i32 %49, 0
br i1 %50, label %55, label %51
51: ; preds = %47
%52 = getelementptr inbounds i32, i32* %3, i64 %26
store i32 1, i32* %52
%53 = load i32, i32* %19
%54 = add nsw i32 %53, 1
store i32 %54, i32* %19
br label %57
55: ; preds = %41, %47, %43
%56 = getelementptr inbounds i32, i32* %3, i64 %26
store i32 0, i32* %56
br label %57
57: ; preds = %51, %55
%58 = add nuw nsw i64 %24, 1
%59 = icmp eq i64 %58, %16
br i1 %59, label %60, label %23
60: ; preds = %57, %17
%61 = add nuw nsw i64 %18, 1
%62 = icmp eq i64 %61, %15
br i1 %62, label %63, label %17
63: ; preds = %60
store i32 0, i32* %6
br i1 %10, label %64, label %101
64: ; preds = %63
%65 = and i64 %15, 1
%66 = icmp eq i32 %7, 1
br i1 %66, label %91, label %67
67: ; preds = %64
%68 = and i64 %15, 4294967294
br label %69
69: ; preds = %86, %67
%70 = phi i32 [ 0, %67 ], [ %87, %86 ]
%71 = phi i64 [ 0, %67 ], [ %88, %86 ]
%72 = phi i64 [ 0, %67 ], [ %89, %86 ]
%73 = getelementptr inbounds i32, i32* %1, i64 %71
%74 = load i32, i32* %73
%75 = icmp eq i32 %74, 0
br i1 %75, label %78, label %76
76: ; preds = %69
%77 = add nsw i32 %70, 1
store i32 %77, i32* %6
br label %78
78: ; preds = %69, %76
%79 = phi i32 [ %70, %69 ], [ %77, %76 ]
%80 = or i64 %71, 1
%81 = getelementptr inbounds i32, i32* %1, i64 %80
%82 = load i32, i32* %81
%83 = icmp eq i32 %82, 0
br i1 %83, label %86, label %84
84: ; preds = %78
%85 = add nsw i32 %79, 1
store i32 %85, i32* %6
br label %86
86: ; preds = %84, %78
%87 = phi i32 [ %79, %78 ], [ %85, %84 ]
%88 = add nuw nsw i64 %71, 2
%89 = add i64 %72, 2
%90 = icmp eq i64 %89, %68
br i1 %90, label %91, label %69
91: ; preds = %86, %64
%92 = phi i32 [ 0, %64 ], [ %87, %86 ]
%93 = phi i64 [ 0, %64 ], [ %88, %86 ]
%94 = icmp eq i64 %65, 0
br i1 %94, label %101, label %95
95: ; preds = %91
%96 = getelementptr inbounds i32, i32* %1, i64 %93
%97 = load i32, i32* %96
%98 = icmp eq i32 %97, 0
br i1 %98, label %101, label %99
99: ; preds = %95
%100 = add nsw i32 %92, 1
store i32 %100, i32* %6
br label %101
101: ; preds = %91, %99, %95, %11, %63
ret void
}