int active_(int n, double *l, double *u, int *nbd, double *x, int *iwhere,
            int *prjctd, int *cnstnd, int *boxed) {
  int i, nbdd;
  nbdd = 0;
  *prjctd = 0;
  *cnstnd = 0;
  *boxed = 1;
  for (i = 0; i < n; ++i) {
    if (nbd[i] > 0) {
      if (nbd[i] <= 2 && x[i] <= l[i]) {
        if (x[i] < l[i]) {
          *prjctd = 1;
          x[i] = l[i];
        }
        ++nbdd;
      } else if (nbd[i] >= 2 && x[i] >= u[i]) {
        if (x[i] > u[i]) {
          *prjctd = 1;
          x[i] = u[i];
        }
        ++nbdd;
      }
    }
  }
  for (i = 0; i < n; ++i) {
    if (nbd[i] != 2) {
      *boxed = 0;
    }
    if (nbd[i] == 0) {
      iwhere[i] = -1;
    } else {
      *cnstnd = 1;
      if (nbd[i] == 2 && u[i] - l[i] <= 0.) {
        iwhere[i] = 3;
      } else {
        iwhere[i] = 0;
      }
    }
  }
  return 0;
}


x86 assembly (clang, -O0):
<active_>:
push   %rbp
mov    %rsp,%rbp
mov    0x20(%rbp),%rax
mov    0x18(%rbp),%rax
mov    0x10(%rbp),%rax
mov    %edi,-0x4(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %rcx,-0x20(%rbp)
mov    %r8,-0x28(%rbp)
mov    %r9,-0x30(%rbp)
movl   $0x0,-0x38(%rbp)
mov    0x10(%rbp),%rax
movl   $0x0,(%rax)
mov    0x18(%rbp),%rax
movl   $0x0,(%rax)
mov    0x20(%rbp),%rax
movl   $0x1,(%rax)
movl   $0x0,-0x34(%rbp)
mov    -0x34(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    1299 <active_+0x199>
mov    -0x20(%rbp),%rax
movslq -0x34(%rbp),%rcx
cmpl   $0x0,(%rax,%rcx,4)
jle    1286 <active_+0x186>
mov    -0x20(%rbp),%rax
movslq -0x34(%rbp),%rcx
cmpl   $0x2,(%rax,%rcx,4)
jg     11fd <active_+0xfd>
mov    -0x28(%rbp),%rax
movslq -0x34(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
mov    -0x10(%rbp),%rax
movslq -0x34(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
ucomisd %xmm1,%xmm0
jb     11fd <active_+0xfd>
mov    -0x28(%rbp),%rax
movslq -0x34(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
mov    -0x10(%rbp),%rax
movslq -0x34(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
ucomisd %xmm1,%xmm0
jbe    11ef <active_+0xef>
mov    0x10(%rbp),%rax
movl   $0x1,(%rax)
mov    -0x10(%rbp),%rax
movslq -0x34(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x28(%rbp),%rax
movslq -0x34(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x38(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x38(%rbp)
jmp    1281 <active_+0x181>
mov    -0x20(%rbp),%rax
movslq -0x34(%rbp),%rcx
cmpl   $0x2,(%rax,%rcx,4)
jl     127c <active_+0x17c>
mov    -0x28(%rbp),%rax
movslq -0x34(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x18(%rbp),%rax
movslq -0x34(%rbp),%rcx
ucomisd (%rax,%rcx,8),%xmm0
jb     127c <active_+0x17c>
mov    -0x28(%rbp),%rax
movslq -0x34(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x18(%rbp),%rax
movslq -0x34(%rbp),%rcx
ucomisd (%rax,%rcx,8),%xmm0
jbe    1273 <active_+0x173>
mov    0x10(%rbp),%rax
movl   $0x1,(%rax)
mov    -0x18(%rbp),%rax
movslq -0x34(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x28(%rbp),%rax
movslq -0x34(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x38(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x38(%rbp)
jmp    1281 <active_+0x181>
jmp    1286 <active_+0x186>
jmp    128b <active_+0x18b>
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    1153 <active_+0x53>
movl   $0x0,-0x34(%rbp)
mov    -0x34(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    136c <active_+0x26c>
mov    -0x20(%rbp),%rax
movslq -0x34(%rbp),%rcx
cmpl   $0x2,(%rax,%rcx,4)
je     12c8 <active_+0x1c8>
mov    0x20(%rbp),%rax
movl   $0x0,(%rax)
mov    -0x20(%rbp),%rax
movslq -0x34(%rbp),%rcx
cmpl   $0x0,(%rax,%rcx,4)
jne    12ee <active_+0x1ee>
mov    -0x30(%rbp),%rax
movslq -0x34(%rbp),%rcx
movl   $0xffffffff,(%rax,%rcx,4)
jmp    1359 <active_+0x259>
mov    0x18(%rbp),%rax
movl   $0x1,(%rax)
mov    -0x20(%rbp),%rax
movslq -0x34(%rbp),%rcx
cmpl   $0x2,(%rax,%rcx,4)
jne    1345 <active_+0x245>
mov    -0x18(%rbp),%rax
movslq -0x34(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
mov    -0x10(%rbp),%rax
movslq -0x34(%rbp),%rcx
subsd  (%rax,%rcx,8),%xmm1
xorps  %xmm0,%xmm0
ucomisd %xmm1,%xmm0
jb     1345 <active_+0x245>
mov    -0x30(%rbp),%rax
movslq -0x34(%rbp),%rcx
movl   $0x3,(%rax,%rcx,4)
jmp    1354 <active_+0x254>
mov    -0x30(%rbp),%rax
movslq -0x34(%rbp),%rcx
movl   $0x0,(%rax,%rcx,4)
jmp    1359 <active_+0x259>
jmp    135e <active_+0x25e>
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    12a0 <active_+0x1a0>
xor    %eax,%eax
pop    %rbp
ret

x86 assembly (clang, -O2):
<active_>:
push   %rbp
push   %r14
push   %rbx
mov    0x30(%rsp),%r10
mov    0x28(%rsp),%r11
mov    0x20(%rsp),%r14
movl   $0x0,(%r14)
movl   $0x0,(%r11)
movl   $0x1,(%r10)
test   %edi,%edi
jle    11ed <active_+0xed>
mov    %edi,%eax
xor    %ebx,%ebx
jmp    1149 <active_+0x49>
cs nopw 0x0(%rax,%rax,1)
add    $0x1,%rbx
cmp    %rbx,%rax
je     118a <active_+0x8a>
mov    (%rcx,%rbx,4),%ebp
test   %ebp,%ebp
jle    1140 <active_+0x40>
movsd  (%r8,%rbx,8),%xmm0
cmp    $0x2,%ebp
ja     1170 <active_+0x70>
movsd  (%rsi,%rbx,8),%xmm1
ucomisd %xmm0,%xmm1
jae    1179 <active_+0x79>
cmp    $0x2,%ebp
jb     1140 <active_+0x40>
nopl   0x0(%rax,%rax,1)
movsd  (%rdx,%rbx,8),%xmm1
ucomisd %xmm1,%xmm0
jbe    1140 <active_+0x40>
movl   $0x1,(%r14)
movsd  %xmm1,(%r8,%rbx,8)
jmp    1140 <active_+0x40>
test   %edi,%edi
jle    11ed <active_+0xed>
xor    %edi,%edi
xorpd  %xmm0,%xmm0
jmp    11af <active_+0xaf>
cs nopw 0x0(%rax,%rax,1)
xor    %ebp,%ebp
mov    %ebp,(%r9,%rdi,4)
add    $0x1,%rdi
cmp    %rdi,%rax
je     11ed <active_+0xed>
cmpl   $0x2,(%rcx,%rdi,4)
je     11c2 <active_+0xc2>
movl   $0x0,(%r10)
cmpl   $0x0,(%rcx,%rdi,4)
je     11e6 <active_+0xe6>
movl   $0x1,(%r11)
cmpl   $0x2,(%rcx,%rdi,4)
jne    11a0 <active_+0xa0>
movsd  (%rdx,%rdi,8),%xmm1
subsd  (%rsi,%rdi,8),%xmm1
mov    $0x3,%ebp
ucomisd %xmm1,%xmm0
jb     11a0 <active_+0xa0>
jmp    11a2 <active_+0xa2>
mov    $0xffffffff,%ebp
jmp    11a2 <active_+0xa2>
xor    %eax,%eax
pop    %rbx
pop    %r14
pop    %rbp
ret

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

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