#define NULL ((void *)0)
typedef unsigned long size_t; // Customize by platform.
typedef long intptr_t;
typedef unsigned long uintptr_t;
typedef long scalar_t__; // Either arithmetic or pointer type.
/* By default, we understand bool (as a convenience). */
typedef int bool;
#define false 0
#define true 1

/* Forward declarations */

/* Type definitions */

/* Variables and functions */
int /*<<< orphan*/ clr_HasBMS(unsigned char *);
int /*<<< orphan*/ exit(int);
int /*<<< orphan*/ fprintf(int /*<<< orphan*/, char *);
int /*<<< orphan*/ free(unsigned char *);
scalar_t__ malloc(int);
int /*<<< orphan*/ set_HasBMS(unsigned char *);
int /*<<< orphan*/ set_int3(unsigned char *, int);
int /*<<< orphan*/ stderr;
unsigned char *mk_BMS(unsigned char *pds, float *bindata, int *n,
                      float undef_low, float undef_hi) {
  int bms_size, i, nn, start;
  unsigned int imask, c;
  unsigned char *bms, *cbits;
  if (undef_low > undef_hi) {
    clr_HasBMS(pds);
    return NULL;
  }
  nn = *n;
  bms_size = 6 + (nn + 7) / 8;
  if (bms_size & 1)
    bms_size++;
  if ((bms = (unsigned char *)malloc(bms_size * sizeof(char))) == NULL) {
    fprintf(stderr, "*** not enough memory: mk_bms ***\n");
    exit(8);
  }
  set_int3(bms, bms_size);
  bms[3] = bms_size * 8 - 6 * 8 - nn;
  bms[4] = bms[5] = bms[bms_size - 1] = 0;
  cbits = bms + 6;
  c = start = 0;
  imask = 128;
  for (i = 0; i < nn; i++) {
    if (bindata[i] < undef_low || bindata[i] > undef_hi) {
      c += imask;
      bindata[start++] = bindata[i];
    }
    if ((imask >>= 1) == 0) {
      *cbits++ = c;
      c = 0;
      imask = 128;
    }
  }
  if (imask != 128)
    *cbits = c;
  if (nn == start) {
    free(bms);
    clr_HasBMS(pds);
    return NULL;
  }
  set_HasBMS(pds);
  *n = start;
  return bms;
}


x86 assembly (clang, -O0):
<mk_BMS>:
push   %rbp
mov    %rsp,%rbp
sub    $0x50,%rsp
mov    %rdi,-0x10(%rbp)
mov    %rsi,-0x18(%rbp)
mov    %rdx,-0x20(%rbp)
movss  %xmm0,-0x24(%rbp)
movss  %xmm1,-0x28(%rbp)
movss  -0x24(%rbp),%xmm0
ucomiss -0x28(%rbp),%xmm0
jbe    11b3 <mk_BMS+0x43>
mov    -0x10(%rbp),%rdi
call   1050 <clr_HasBMS@plt>
movq   $0x0,-0x8(%rbp)
jmp    139f <mk_BMS+0x22f>
mov    -0x20(%rbp),%rax
mov    (%rax),%eax
mov    %eax,-0x34(%rbp)
mov    -0x34(%rbp),%eax
add    $0x7,%eax
mov    $0x8,%ecx
cltd
idiv   %ecx
add    $0x6,%eax
mov    %eax,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
and    $0x1,%eax
cmp    $0x0,%eax
je     11e8 <mk_BMS+0x78>
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
movslq -0x2c(%rbp),%rax
shl    $0x0,%rax
mov    %eax,%edi
call   1070 <malloc@plt>
mov    %rax,-0x48(%rbp)
cmp    $0x0,%rax
jne    1224 <mk_BMS+0xb4>
mov    0x2dd4(%rip),%rax
mov    (%rax),%edi
lea    0xdeb(%rip),%rsi
call   1060 <fprintf@plt>
mov    $0x8,%edi
call   1090 <exit@plt>
mov    -0x48(%rbp),%rdi
mov    -0x2c(%rbp),%esi
call   1080 <set_int3@plt>
mov    -0x2c(%rbp),%eax
shl    $0x3,%eax
sub    $0x30,%eax
sub    -0x34(%rbp),%eax
mov    %al,%cl
mov    -0x48(%rbp),%rax
mov    %cl,0x3(%rax)
mov    -0x48(%rbp),%rax
mov    -0x2c(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movb   $0x0,(%rax,%rcx,1)
mov    -0x48(%rbp),%rax
movb   $0x0,0x5(%rax)
mov    -0x48(%rbp),%rax
movb   $0x0,0x4(%rax)
mov    -0x48(%rbp),%rax
add    $0x6,%rax
mov    %rax,-0x50(%rbp)
movl   $0x0,-0x38(%rbp)
movl   $0x0,-0x40(%rbp)
movl   $0x80,-0x3c(%rbp)
movl   $0x0,-0x30(%rbp)
mov    -0x30(%rbp),%eax
cmp    -0x34(%rbp),%eax
jge    1342 <mk_BMS+0x1d2>
mov    -0x18(%rbp),%rax
movslq -0x30(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm1
movss  -0x24(%rbp),%xmm0
ucomiss %xmm1,%xmm0
ja     12cc <mk_BMS+0x15c>
mov    -0x18(%rbp),%rax
movslq -0x30(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
ucomiss -0x28(%rbp),%xmm0
jbe    12f9 <mk_BMS+0x189>
mov    -0x3c(%rbp),%eax
add    -0x40(%rbp),%eax
mov    %eax,-0x40(%rbp)
mov    -0x18(%rbp),%rax
movslq -0x30(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x38(%rbp),%ecx
mov    %ecx,%edx
add    $0x1,%edx
mov    %edx,-0x38(%rbp)
movslq %ecx,%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x3c(%rbp),%eax
shr    $0x1,%eax
mov    %eax,-0x3c(%rbp)
cmp    $0x0,%eax
jne    132f <mk_BMS+0x1bf>
mov    -0x40(%rbp),%eax
mov    %al,%cl
mov    -0x50(%rbp),%rax
mov    %rax,%rdx
add    $0x1,%rdx
mov    %rdx,-0x50(%rbp)
mov    %cl,(%rax)
movl   $0x0,-0x40(%rbp)
movl   $0x80,-0x3c(%rbp)
jmp    1334 <mk_BMS+0x1c4>
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
jmp    128e <mk_BMS+0x11e>
cmpl   $0x80,-0x3c(%rbp)
je     135a <mk_BMS+0x1ea>
mov    -0x40(%rbp),%eax
mov    %al,%cl
mov    -0x50(%rbp),%rax
mov    %cl,(%rax)
mov    -0x34(%rbp),%eax
cmp    -0x38(%rbp),%eax
jne    1385 <mk_BMS+0x215>
mov    -0x48(%rbp),%rdi
call   1030 <free@plt>
mov    -0x10(%rbp),%rdi
call   1050 <clr_HasBMS@plt>
movq   $0x0,-0x8(%rbp)
jmp    139f <mk_BMS+0x22f>
mov    -0x10(%rbp),%rdi
call   1040 <set_HasBMS@plt>
mov    -0x38(%rbp),%ecx
mov    -0x20(%rbp),%rax
mov    %ecx,(%rax)
mov    -0x48(%rbp),%rax
mov    %rax,-0x8(%rbp)
mov    -0x8(%rbp),%rax
add    $0x50,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<mk_BMS>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x18,%rsp
mov    %rdi,%r15
ucomiss %xmm1,%xmm0
ja     1298 <mk_BMS+0x128>
mov    %rdx,%rbp
mov    %rsi,%rbx
movss  %xmm0,0x8(%rsp)
movss  %xmm1,0xc(%rsp)
mov    (%rdx),%r14d
lea    0x7(%r14),%eax
lea    0xe(%r14),%ecx
test   %eax,%eax
cmovns %eax,%ecx
sar    $0x3,%ecx
mov    %ecx,%eax
and    $0x1,%eax
lea    (%rcx,%rax,1),%r13d
add    $0x6,%r13d
mov    %r13d,%edi
call   1070 <malloc@plt>
test   %rax,%rax
je     12cd <mk_BMS+0x15d>
mov    %rax,%r12
mov    %rbp,0x10(%rsp)
mov    %rax,%rdi
mov    %r13d,%esi
call   1080 <set_int3@plt>
lea    0x0(,%r13,8),%eax
sub    %r14d,%eax
add    $0xd0,%al
mov    %al,0x3(%r12)
movslq %r13d,%rax
movb   $0x0,-0x1(%rax,%r12,1)
movw   $0x0,0x4(%r12)
test   %r14d,%r14d
jle    12a5 <mk_BMS+0x135>
mov    %r12,%rax
add    $0x6,%rax
mov    $0x80,%esi
xor    %ecx,%ecx
xor    %edx,%edx
xor    %ebp,%ebp
movss  0xc(%rsp),%xmm1
movss  0x8(%rsp),%xmm2
jmp    123b <mk_BMS+0xcb>
nopl   0x0(%rax)
shr    %esi
add    $0x1,%rcx
cmp    %rcx,%r14
je     1281 <mk_BMS+0x111>
movss  (%rbx,%rcx,4),%xmm0
ucomiss %xmm0,%xmm2
ja     1260 <mk_BMS+0xf0>
ucomiss %xmm1,%xmm0
ja     1260 <mk_BMS+0xf0>
cmp    $0x1,%esi
ja     1230 <mk_BMS+0xc0>
jmp    1272 <mk_BMS+0x102>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
add    %esi,%edx
movslq %ebp,%rdi
add    $0x1,%ebp
movss  %xmm0,(%rbx,%rdi,4)
cmp    $0x1,%esi
ja     1230 <mk_BMS+0xc0>
mov    %dl,(%rax)
add    $0x1,%rax
xor    %edx,%edx
mov    $0x80,%esi
jmp    1232 <mk_BMS+0xc2>
cmp    $0x80,%esi
je     128b <mk_BMS+0x11b>
mov    %dl,(%rax)
cmp    %ebp,%r14d
jne    12ac <mk_BMS+0x13c>
mov    %r12,%rdi
call   1030 <free@plt>
mov    %r15,%rdi
call   1050 <clr_HasBMS@plt>
xor    %r12d,%r12d
jmp    12bb <mk_BMS+0x14b>
xor    %ebp,%ebp
cmp    %ebp,%r14d
je     1290 <mk_BMS+0x120>
mov    %r15,%rdi
call   1040 <set_HasBMS@plt>
mov    0x10(%rsp),%rax
mov    %ebp,(%rax)
mov    %r12,%rax
add    $0x18,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret
mov    0x2d0c(%rip),%rax
mov    (%rax),%edi
lea    0xd23(%rip),%rsi
call   1060 <fprintf@plt>
mov    $0x8,%edi
call   1090 <exit@plt>

LLVM IR (-O0):
@stderr = global i32 0
@.str = [35 x i8] c"*** not enough memory: mk_bms ***\0A\00"
define i8* @mk_BMS(i8* %0, float* %1, i32* %2, float %3, float %4) {
%6 = alloca i8*
%7 = alloca i8*
%8 = alloca float*
%9 = alloca i32*
%10 = alloca float
%11 = alloca float
%12 = alloca i32
%13 = alloca i32
%14 = alloca i32
%15 = alloca i32
%16 = alloca i32
%17 = alloca i32
%18 = alloca i8*
%19 = alloca i8*
store i8* %0, i8** %7
store float* %1, float** %8
store i32* %2, i32** %9
store float %3, float* %10
store float %4, float* %11
%20 = load float, float* %10
%21 = load float, float* %11
%22 = fcmp ogt float %20, %21
br i1 %22, label %23, label %26
23: ; preds = %5
%24 = load i8*, i8** %7
%25 = call i32 @clr_HasBMS(i8* %24)
store i8* null, i8** %6
br label %143
26: ; preds = %5
%27 = load i32*, i32** %9
%28 = load i32, i32* %27
store i32 %28, i32* %14
%29 = load i32, i32* %14
%30 = add nsw i32 %29, 7
%31 = sdiv i32 %30, 8
%32 = add nsw i32 6, %31
store i32 %32, i32* %12
%33 = load i32, i32* %12
%34 = and i32 %33, 1
%35 = icmp ne i32 %34, 0
br i1 %35, label %36, label %39
36: ; preds = %26
%37 = load i32, i32* %12
%38 = add nsw i32 %37, 1
store i32 %38, i32* %12
br label %39
39: ; preds = %36, %26
%40 = load i32, i32* %12
%41 = sext i32 %40 to i64
%42 = mul i64 %41, 1
%43 = trunc i64 %42 to i32
%44 = call i64 @malloc(i32 %43)
%45 = inttoptr i64 %44 to i8*
store i8* %45, i8** %18
%46 = icmp eq i8* %45, null
br i1 %46, label %47, label %51
47: ; preds = %39
%48 = load i32, i32* @stderr
%49 = call i32 @fprintf(i32 %48, i8* getelementptr inbounds ([35 x i8], [35 x i8]* @.str, i64 0, i64 0))
%50 = call i32 @exit(i32 8)
unreachable
51: ; preds = %39
%52 = load i8*, i8** %18
%53 = load i32, i32* %12
%54 = call i32 @set_int3(i8* %52, i32 %53)
%55 = load i32, i32* %12
%56 = mul nsw i32 %55, 8
%57 = sub nsw i32 %56, 48
%58 = load i32, i32* %14
%59 = sub nsw i32 %57, %58
%60 = trunc i32 %59 to i8
%61 = load i8*, i8** %18
%62 = getelementptr inbounds i8, i8* %61, i64 3
store i8 %60, i8* %62
%63 = load i8*, i8** %18
%64 = load i32, i32* %12
%65 = sub nsw i32 %64, 1
%66 = sext i32 %65 to i64
%67 = getelementptr inbounds i8, i8* %63, i64 %66
store i8 0, i8* %67
%68 = load i8*, i8** %18
%69 = getelementptr inbounds i8, i8* %68, i64 5
store i8 0, i8* %69
%70 = load i8*, i8** %18
%71 = getelementptr inbounds i8, i8* %70, i64 4
store i8 0, i8* %71
%72 = load i8*, i8** %18
%73 = getelementptr inbounds i8, i8* %72, i64 6
store i8* %73, i8** %19
store i32 0, i32* %15
store i32 0, i32* %17
store i32 128, i32* %16
store i32 0, i32* %13
br label %74
74: ; preds = %118, %51
%75 = load i32, i32* %13
%76 = load i32, i32* %14
%77 = icmp slt i32 %75, %76
br i1 %77, label %78, label %121
78: ; preds = %74
%79 = load float*, float** %8
%80 = load i32, i32* %13
%81 = sext i32 %80 to i64
%82 = getelementptr inbounds float, float* %79, i64 %81
%83 = load float, float* %82
%84 = load float, float* %10
%85 = fcmp olt float %83, %84
br i1 %85, label %94, label %86
86: ; preds = %78
%87 = load float*, float** %8
%88 = load i32, i32* %13
%89 = sext i32 %88 to i64
%90 = getelementptr inbounds float, float* %87, i64 %89
%91 = load float, float* %90
%92 = load float, float* %11
%93 = fcmp ogt float %91, %92
br i1 %93, label %94, label %108
94: ; preds = %86, %78
%95 = load i32, i32* %16
%96 = load i32, i32* %17
%97 = add i32 %96, %95
store i32 %97, i32* %17
%98 = load float*, float** %8
%99 = load i32, i32* %13
%100 = sext i32 %99 to i64
%101 = getelementptr inbounds float, float* %98, i64 %100
%102 = load float, float* %101
%103 = load float*, float** %8
%104 = load i32, i32* %15
%105 = add nsw i32 %104, 1
store i32 %105, i32* %15
%106 = sext i32 %104 to i64
%107 = getelementptr inbounds float, float* %103, i64 %106
store float %102, float* %107
br label %108
108: ; preds = %94, %86
%109 = load i32, i32* %16
%110 = lshr i32 %109, 1
store i32 %110, i32* %16
%111 = icmp eq i32 %110, 0
br i1 %111, label %112, label %117
112: ; preds = %108
%113 = load i32, i32* %17
%114 = trunc i32 %113 to i8
%115 = load i8*, i8** %19
%116 = getelementptr inbounds i8, i8* %115, i32 1
store i8* %116, i8** %19
store i8 %114, i8* %115
store i32 0, i32* %17
store i32 128, i32* %16
br label %117
117: ; preds = %112, %108
br label %118
118: ; preds = %117
%119 = load i32, i32* %13
%120 = add nsw i32 %119, 1
store i32 %120, i32* %13
br label %74
121: ; preds = %74
%122 = load i32, i32* %16
%123 = icmp ne i32 %122, 128
br i1 %123, label %124, label %128
124: ; preds = %121
%125 = load i32, i32* %17
%126 = trunc i32 %125 to i8
%127 = load i8*, i8** %19
store i8 %126, i8* %127
br label %128
128: ; preds = %124, %121
%129 = load i32, i32* %14
%130 = load i32, i32* %15
%131 = icmp eq i32 %129, %130
br i1 %131, label %132, label %137
132: ; preds = %128
%133 = load i8*, i8** %18
%134 = call i32 @free(i8* %133)
%135 = load i8*, i8** %7
%136 = call i32 @clr_HasBMS(i8* %135)
store i8* null, i8** %6
br label %143
137: ; preds = %128
%138 = load i8*, i8** %7
%139 = call i32 @set_HasBMS(i8* %138)
%140 = load i32, i32* %15
%141 = load i32*, i32** %9
store i32 %140, i32* %141
%142 = load i8*, i8** %18
store i8* %142, i8** %6
br label %143
143: ; preds = %137, %132, %23
%144 = load i8*, i8** %6
ret i8* %144
}

LLVM IR (-O2):
@stderr = global i32 0
@.str = [35 x i8] c"*** not enough memory: mk_bms ***\0A\00"
define i8* @mk_BMS(i8* %0, float* nocapture %1, i32* nocapture %2, float %3, float %4) {
%6 = fcmp ogt float %3, %4
br i1 %6, label %7, label %9
7: ; preds = %5
%8 = tail call i32 @clr_HasBMS(i8* %0)
br label %81
9: ; preds = %5
%10 = load i32, i32* %2
%11 = add nsw i32 %10, 7
%12 = sdiv i32 %11, 8
%13 = and i32 %12, 1
%14 = or i32 %13, 6
%15 = add nsw i32 %14, %12
%16 = tail call i64 @malloc(i32 %15)
%17 = inttoptr i64 %16 to i8*
%18 = icmp eq i64 %16, 0
br i1 %18, label %19, label %23
19: ; preds = %9
%20 = load i32, i32* @stderr
%21 = tail call i32 @fprintf(i32 %20, i8* getelementptr inbounds ([35 x i8], [35 x i8]* @.str, i64 0, i64 0))
%22 = tail call i32 @exit(i32 8)
unreachable
23: ; preds = %9
%24 = tail call i32 @set_int3(i8* %17, i32 %15)
%25 = shl nsw i32 %15, 3
%26 = sub i32 208, %10
%27 = add i32 %26, %25
%28 = trunc i32 %27 to i8
%29 = getelementptr inbounds i8, i8* %17, i64 3
store i8 %28, i8* %29
%30 = add nsw i32 %15, -1
%31 = sext i32 %30 to i64
%32 = getelementptr inbounds i8, i8* %17, i64 %31
store i8 0, i8* %32
%33 = getelementptr inbounds i8, i8* %17, i64 5
store i8 0, i8* %33
%34 = getelementptr inbounds i8, i8* %17, i64 4
store i8 0, i8* %34
%35 = icmp sgt i32 %10, 0
br i1 %35, label %36, label %73
36: ; preds = %23
%37 = getelementptr inbounds i8, i8* %17, i64 6
%38 = zext i32 %10 to i64
br label %39
39: ; preds = %36, %63
%40 = phi i64 [ 0, %36 ], [ %67, %63 ]
%41 = phi i8* [ %37, %36 ], [ %66, %63 ]
%42 = phi i32 [ 0, %36 ], [ %65, %63 ]
%43 = phi i32 [ 128, %36 ], [ %64, %63 ]
%44 = phi i32 [ 0, %36 ], [ %56, %63 ]
%45 = getelementptr inbounds float, float* %1, i64 %40
%46 = load float, float* %45
%47 = fcmp olt float %46, %3
%48 = fcmp ogt float %46, %4
%49 = or i1 %47, %48
br i1 %49, label %50, label %55
50: ; preds = %39
%51 = add i32 %42, %43
%52 = add nsw i32 %44, 1
%53 = sext i32 %44 to i64
%54 = getelementptr inbounds float, float* %1, i64 %53
store float %46, float* %54
br label %55
55: ; preds = %39, %50
%56 = phi i32 [ %52, %50 ], [ %44, %39 ]
%57 = phi i32 [ %51, %50 ], [ %42, %39 ]
%58 = lshr i32 %43, 1
%59 = icmp ult i32 %43, 2
br i1 %59, label %60, label %63
60: ; preds = %55
%61 = trunc i32 %57 to i8
%62 = getelementptr inbounds i8, i8* %41, i64 1
store i8 %61, i8* %41
br label %63
63: ; preds = %55, %60
%64 = phi i32 [ 128, %60 ], [ %58, %55 ]
%65 = phi i32 [ 0, %60 ], [ %57, %55 ]
%66 = phi i8* [ %62, %60 ], [ %41, %55 ]
%67 = add nuw nsw i64 %40, 1
%68 = icmp eq i64 %67, %38
br i1 %68, label %69, label %39
69: ; preds = %63
%70 = icmp eq i32 %64, 128
br i1 %70, label %73, label %71
71: ; preds = %69
%72 = trunc i32 %65 to i8
store i8 %72, i8* %66
br label %73
73: ; preds = %23, %71, %69
%74 = phi i32 [ %56, %71 ], [ %56, %69 ], [ 0, %23 ]
%75 = icmp eq i32 %10, %74
br i1 %75, label %76, label %79
76: ; preds = %73
%77 = tail call i32 @free(i8* %17)
%78 = tail call i32 @clr_HasBMS(i8* %0)
br label %81
79: ; preds = %73
%80 = tail call i32 @set_HasBMS(i8* %0)
store i32 %74, i32* %2
br label %81
81: ; preds = %79, %76, %7
%82 = phi i8* [ null, %7 ], [ null, %76 ], [ %17, %79 ]
ret i8* %82
}