int get_changedRxns_by_one_met(int *IA, double *VA, int lenA, double *lb,
                               double *ub, int *tmprun, int *changedRxn) {
  int i, j, k;
  k = 0;
  for (i = 1; i < lenA + 1; i++) {
    if ((lb[IA[i] - 1] > -1e-6 && VA[i] < -1e-6) ||
        (ub[IA[i] - 1] < 1e-6 && VA[i] > 1e-6)) {
      if (tmprun[IA[i] - 1] == 0) {
        changedRxn[k] = IA[i] - 1;
        k = k + 1;
        tmprun[IA[i] - 1] = 1;
      }
    }
  }
  for (i = 0; i < k; i++) {
    tmprun[changedRxn[i]] = 0;
  }
  return k;
}


x86 assembly (clang, -O0):
<get_changedRxns_by_one_met>:
push   %rbp
mov    %rsp,%rbp
mov    0x10(%rbp),%rax
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %edx,-0x14(%rbp)
mov    %rcx,-0x20(%rbp)
mov    %r8,-0x28(%rbp)
mov    %r9,-0x30(%rbp)
movl   $0x0,-0x3c(%rbp)
movl   $0x1,-0x34(%rbp)
mov    -0x34(%rbp),%eax
mov    -0x14(%rbp),%ecx
add    $0x1,%ecx
cmp    %ecx,%eax
jge    1241 <get_changedRxns_by_one_met+0x141>
mov    -0x20(%rbp),%rax
mov    -0x8(%rbp),%rcx
movslq -0x34(%rbp),%rdx
mov    (%rcx,%rdx,4),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  0xea0(%rip),%xmm1
ucomisd %xmm1,%xmm0
jbe    1185 <get_changedRxns_by_one_met+0x85>
mov    -0x10(%rbp),%rax
movslq -0x34(%rbp),%rcx
movsd  0xe86(%rip),%xmm0
ucomisd (%rax,%rcx,8),%xmm0
ja     11cc <get_changedRxns_by_one_met+0xcc>
mov    -0x28(%rbp),%rax
mov    -0x8(%rbp),%rcx
movslq -0x34(%rbp),%rdx
mov    (%rcx,%rdx,4),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  0xe66(%rip),%xmm0
ucomisd (%rax,%rcx,8),%xmm0
jbe    122e <get_changedRxns_by_one_met+0x12e>
mov    -0x10(%rbp),%rax
movslq -0x34(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  0xe46(%rip),%xmm1
ucomisd %xmm1,%xmm0
jbe    122e <get_changedRxns_by_one_met+0x12e>
mov    -0x30(%rbp),%rax
mov    -0x8(%rbp),%rcx
movslq -0x34(%rbp),%rdx
mov    (%rcx,%rdx,4),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
cmpl   $0x0,(%rax,%rcx,4)
jne    1229 <get_changedRxns_by_one_met+0x129>
mov    -0x8(%rbp),%rax
movslq -0x34(%rbp),%rcx
mov    (%rax,%rcx,4),%edx
sub    $0x1,%edx
mov    0x10(%rbp),%rax
movslq -0x3c(%rbp),%rcx
mov    %edx,(%rax,%rcx,4)
mov    -0x3c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x3c(%rbp)
mov    -0x30(%rbp),%rax
mov    -0x8(%rbp),%rcx
movslq -0x34(%rbp),%rdx
mov    (%rcx,%rdx,4),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movl   $0x1,(%rax,%rcx,4)
jmp    122e <get_changedRxns_by_one_met+0x12e>
jmp    1233 <get_changedRxns_by_one_met+0x133>
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    112d <get_changedRxns_by_one_met+0x2d>
movl   $0x0,-0x34(%rbp)
mov    -0x34(%rbp),%eax
cmp    -0x3c(%rbp),%eax
jge    1279 <get_changedRxns_by_one_met+0x179>
mov    -0x30(%rbp),%rax
mov    0x10(%rbp),%rcx
movslq -0x34(%rbp),%rdx
movslq (%rcx,%rdx,4),%rcx
movl   $0x0,(%rax,%rcx,4)
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    1248 <get_changedRxns_by_one_met+0x148>
mov    -0x3c(%rbp),%eax
pop    %rbp
ret

x86 assembly (clang, -O2):
<get_changedRxns_by_one_met>:
push   %rbx
test   %edx,%edx
jle    11a9 <get_changedRxns_by_one_met+0xa9>
mov    0x10(%rsp),%r10
add    $0x1,%edx
xor    %eax,%eax
mov    $0x1,%r11d
movsd  0xedf(%rip),%xmm0
movsd  0xedf(%rip),%xmm1
jmp    1141 <get_changedRxns_by_one_met+0x41>
nopl   0x0(%rax,%rax,1)
cmpl   $0x0,-0x4(%r9,%rbx,4)
je     1170 <get_changedRxns_by_one_met+0x70>
add    $0x1,%r11
cmp    %r11,%rdx
je     118c <get_changedRxns_by_one_met+0x8c>
movslq (%rdi,%r11,4),%rbx
movsd  -0x8(%rcx,%rbx,8),%xmm2
ucomisd %xmm0,%xmm2
jbe    1159 <get_changedRxns_by_one_met+0x59>
ucomisd (%rsi,%r11,8),%xmm0
ja     1130 <get_changedRxns_by_one_met+0x30>
ucomisd -0x8(%r8,%rbx,8),%xmm1
jbe    1138 <get_changedRxns_by_one_met+0x38>
movsd  (%rsi,%r11,8),%xmm2
ucomisd %xmm1,%xmm2
ja     1130 <get_changedRxns_by_one_met+0x30>
jmp    1138 <get_changedRxns_by_one_met+0x38>
add    $0xffffffffffffffff,%rbx
cltq
mov    %ebx,(%r10,%rax,4)
add    $0x1,%eax
movslq (%rdi,%r11,4),%rbx
movl   $0x1,-0x4(%r9,%rbx,4)
jmp    1138 <get_changedRxns_by_one_met+0x38>
test   %eax,%eax
jle    1225 <get_changedRxns_by_one_met+0x125>
mov    %eax,%edx
lea    -0x1(%rdx),%rsi
mov    %edx,%ecx
and    $0x3,%ecx
cmp    $0x3,%rsi
jae    11ad <get_changedRxns_by_one_met+0xad>
xor    %esi,%esi
jmp    11fc <get_changedRxns_by_one_met+0xfc>
xor    %eax,%eax
jmp    1225 <get_changedRxns_by_one_met+0x125>
and    $0xfffffffc,%edx
xor    %esi,%esi
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
movslq (%r10,%rsi,4),%rdi
movl   $0x0,(%r9,%rdi,4)
movslq 0x4(%r10,%rsi,4),%rdi
movl   $0x0,(%r9,%rdi,4)
movslq 0x8(%r10,%rsi,4),%rdi
movl   $0x0,(%r9,%rdi,4)
movslq 0xc(%r10,%rsi,4),%rdi
movl   $0x0,(%r9,%rdi,4)
add    $0x4,%rsi
cmp    %rsi,%rdx
jne    11c0 <get_changedRxns_by_one_met+0xc0>
test   %rcx,%rcx
je     1225 <get_changedRxns_by_one_met+0x125>
lea    (%r10,%rsi,4),%rdx
xor    %esi,%esi
nopw   0x0(%rax,%rax,1)
movslq (%rdx,%rsi,4),%rdi
movl   $0x0,(%r9,%rdi,4)
add    $0x1,%rsi
cmp    %rsi,%rcx
jne    1210 <get_changedRxns_by_one_met+0x110>
pop    %rbx
ret

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

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