int _variable_cols_double_approx(int *cols, double *X, int M, int N, double tol,
                                 int min_constant) {
  int i, j;
  int ro = (M - 1) * N;
  double diff;
  int nconstant = N;
  for (j = 0; j < N; j++)
    cols[j] = 0;
  for (i = 0; i < M; i++) {
    ro = i * N;
    for (j = 0; j < N; j++) {
      diff = X[j] - X[ro + j];
      if (diff > tol || -diff > tol) {
        if (cols[j] == 0) {
          cols[j] = 1;
          nconstant--;
          if (nconstant < min_constant)
            return 0;
          if (nconstant == 0)
            return 1;
        }
      }
    }
  }
  return 1;
}


x86 assembly (clang, -O0):
<_variable_cols_double_approx>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x10(%rbp)
mov    %rsi,-0x18(%rbp)
mov    %edx,-0x1c(%rbp)
mov    %ecx,-0x20(%rbp)
movsd  %xmm0,-0x28(%rbp)
mov    %r8d,-0x2c(%rbp)
mov    -0x1c(%rbp),%eax
sub    $0x1,%eax
imul   -0x20(%rbp),%eax
mov    %eax,-0x38(%rbp)
mov    -0x20(%rbp),%eax
mov    %eax,-0x44(%rbp)
movl   $0x0,-0x34(%rbp)
mov    -0x34(%rbp),%eax
cmp    -0x20(%rbp),%eax
jge    115e <_variable_cols_double_approx+0x5e>
mov    -0x10(%rbp),%rax
movslq -0x34(%rbp),%rcx
movl   $0x0,(%rax,%rcx,4)
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    1135 <_variable_cols_double_approx+0x35>
movl   $0x0,-0x30(%rbp)
mov    -0x30(%rbp),%eax
cmp    -0x1c(%rbp),%eax
jge    1271 <_variable_cols_double_approx+0x171>
mov    -0x30(%rbp),%eax
imul   -0x20(%rbp),%eax
mov    %eax,-0x38(%rbp)
movl   $0x0,-0x34(%rbp)
mov    -0x34(%rbp),%eax
cmp    -0x20(%rbp),%eax
jge    125e <_variable_cols_double_approx+0x15e>
mov    -0x18(%rbp),%rax
movslq -0x34(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x38(%rbp),%ecx
add    -0x34(%rbp),%ecx
movslq %ecx,%rcx
subsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x40(%rbp)
movsd  -0x40(%rbp),%xmm0
ucomisd -0x28(%rbp),%xmm0
ja     11e9 <_variable_cols_double_approx+0xe9>
movsd  -0x40(%rbp),%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
ucomisd -0x28(%rbp),%xmm0
jbe    124b <_variable_cols_double_approx+0x14b>
mov    -0x10(%rbp),%rax
movslq -0x34(%rbp),%rcx
cmpl   $0x0,(%rax,%rcx,4)
jne    1246 <_variable_cols_double_approx+0x146>
mov    -0x10(%rbp),%rax
movslq -0x34(%rbp),%rcx
movl   $0x1,(%rax,%rcx,4)
mov    -0x44(%rbp),%eax
add    $0xffffffff,%eax
mov    %eax,-0x44(%rbp)
mov    -0x44(%rbp),%eax
cmp    -0x2c(%rbp),%eax
jge    122b <_variable_cols_double_approx+0x12b>
movl   $0x0,-0x4(%rbp)
jmp    1278 <_variable_cols_double_approx+0x178>
cmpl   $0x0,-0x44(%rbp)
jne    1241 <_variable_cols_double_approx+0x141>
movl   $0x1,-0x4(%rbp)
jmp    1278 <_variable_cols_double_approx+0x178>
jmp    1246 <_variable_cols_double_approx+0x146>
jmp    124b <_variable_cols_double_approx+0x14b>
jmp    1250 <_variable_cols_double_approx+0x150>
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    1182 <_variable_cols_double_approx+0x82>
jmp    1263 <_variable_cols_double_approx+0x163>
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
jmp    1165 <_variable_cols_double_approx+0x65>
movl   $0x1,-0x4(%rbp)
mov    -0x4(%rbp),%eax
pop    %rbp
ret

x86 assembly (clang, -O2):
<_variable_cols_double_approx>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
push   %rax
mov    %r8d,%r15d
mov    %ecx,%r14d
mov    %edx,%r12d
mov    %rsi,%rbp
mov    %rdi,%r13
mov    %ecx,%ebx
test   %ecx,%ecx
jle    114c <_variable_cols_double_approx+0x3c>
lea    0x0(,%rbx,4),%rdx
mov    %r13,%rdi
xor    %esi,%esi
movsd  %xmm0,(%rsp)
call   1030 <memset@plt>
movsd  (%rsp),%xmm0
mov    $0x1,%eax
test   %r12d,%r12d
jle    11e5 <_variable_cols_double_approx+0xd5>
movslq %r14d,%r9
mov    %r12d,%r8d
shl    $0x3,%r9
xor    %esi,%esi
movapd 0xe92(%rip),%xmm2
mov    %rbp,%rdi
mov    %r14d,%edx
jmp    118c <_variable_cols_double_approx+0x7c>
cs nopw 0x0(%rax,%rax,1)
add    $0x1,%rsi
add    %r9,%rdi
cmp    %r8,%rsi
je     11e5 <_variable_cols_double_approx+0xd5>
test   %r14d,%r14d
jle    1180 <_variable_cols_double_approx+0x70>
xor    %ecx,%ecx
jmp    11b1 <_variable_cols_double_approx+0xa1>
cs nopw 0x0(%rax,%rax,1)
nop
cmpl   $0x0,0x0(%r13,%rcx,4)
je     11ce <_variable_cols_double_approx+0xbe>
add    $0x1,%rcx
cmp    %rcx,%rbx
je     1180 <_variable_cols_double_approx+0x70>
movsd  0x0(%rbp,%rcx,8),%xmm1
subsd  (%rdi,%rcx,8),%xmm1
ucomisd %xmm0,%xmm1
ja     11a0 <_variable_cols_double_approx+0x90>
xorpd  %xmm2,%xmm1
ucomisd %xmm0,%xmm1
jbe    11a8 <_variable_cols_double_approx+0x98>
jmp    11a0 <_variable_cols_double_approx+0x90>
movl   $0x1,0x0(%r13,%rcx,4)
cmp    %r15d,%edx
jle    11e3 <_variable_cols_double_approx+0xd3>
add    $0xffffffff,%edx
jne    11a8 <_variable_cols_double_approx+0x98>
jmp    11e5 <_variable_cols_double_approx+0xd5>
xor    %eax,%eax
add    $0x8,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

LLVM IR (-O2):
define i32 @_variable_cols_double_approx(i32* nocapture %0, double* nocapture %1, i32 %2, i32 %3, double %4, i32 %5) {
%7 = icmp sgt i32 %3, 0
br i1 %7, label %8, label %12
8: ; preds = %6
%9 = bitcast i32* %0 to i8*
%10 = zext i32 %3 to i64
%11 = shl nuw nsw i64 %10, 2
br label %12
12: ; preds = %8, %6
%13 = icmp sgt i32 %2, 0
br i1 %13, label %14, label %52
14: ; preds = %12
%15 = sext i32 %3 to i64
%16 = zext i32 %2 to i64
%17 = zext i32 %3 to i64
br label %18
18: ; preds = %14, %48
%19 = phi i64 [ 0, %14 ], [ %50, %48 ]
%20 = phi i32 [ %3, %14 ], [ %49, %48 ]
%21 = mul nsw i64 %19, %15
br i1 %7, label %22, label %48
22: ; preds = %18, %44
%23 = phi i64 [ %46, %44 ], [ 0, %18 ]
%24 = phi i32 [ %45, %44 ], [ %20, %18 ]
%25 = getelementptr inbounds double, double* %1, i64 %23
%26 = load double, double* %25
%27 = add nsw i64 %23, %21
%28 = getelementptr inbounds double, double* %1, i64 %27
%29 = load double, double* %28
%30 = fsub double %26, %29
%31 = fcmp ogt double %30, %4
%32 = fneg double %30
%33 = fcmp ogt double %32, %4
%34 = or i1 %31, %33
br i1 %34, label %35, label %44
35: ; preds = %22
%36 = getelementptr inbounds i32, i32* %0, i64 %23
%37 = load i32, i32* %36
%38 = icmp eq i32 %37, 0
br i1 %38, label %39, label %44
39: ; preds = %35
store i32 1, i32* %36
%40 = add nsw i32 %24, -1
%41 = icmp sgt i32 %24, %5
br i1 %41, label %42, label %52
42: ; preds = %39
%43 = icmp eq i32 %40, 0
br i1 %43, label %52, label %44
44: ; preds = %22, %42, %35
%45 = phi i32 [ %40, %42 ], [ %24, %35 ], [ %24, %22 ]
%46 = add nuw nsw i64 %23, 1
%47 = icmp eq i64 %46, %17
br i1 %47, label %48, label %22
48: ; preds = %44, %18
%49 = phi i32 [ %20, %18 ], [ %45, %44 ]
%50 = add nuw nsw i64 %19, 1
%51 = icmp eq i64 %50, %16
br i1 %51, label %52, label %18
52: ; preds = %48, %42, %39, %12
%53 = phi i32 [ 1, %12 ], [ 1, %42 ], [ 0, %39 ], [ 1, %48 ]
ret i32 %53
}