int _variable_cols_double(int *cols, double *X, int M, int N,
                          int min_constant) {
  int i, j;
  int ro = (M - 1) * N;
  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++) {
      if (X[j] != X[ro + j]) {
        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>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x10(%rbp)
mov    %rsi,-0x18(%rbp)
mov    %edx,-0x1c(%rbp)
mov    %ecx,-0x20(%rbp)
mov    %r8d,-0x24(%rbp)
mov    -0x1c(%rbp),%eax
sub    $0x1,%eax
imul   -0x20(%rbp),%eax
mov    %eax,-0x30(%rbp)
mov    -0x20(%rbp),%eax
mov    %eax,-0x34(%rbp)
movl   $0x0,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0x20(%rbp),%eax
jge    1159 <_variable_cols_double+0x59>
mov    -0x10(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movl   $0x0,(%rax,%rcx,4)
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    1130 <_variable_cols_double+0x30>
movl   $0x0,-0x28(%rbp)
mov    -0x28(%rbp),%eax
cmp    -0x1c(%rbp),%eax
jge    1241 <_variable_cols_double+0x141>
mov    -0x28(%rbp),%eax
imul   -0x20(%rbp),%eax
mov    %eax,-0x30(%rbp)
movl   $0x0,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0x20(%rbp),%eax
jge    122e <_variable_cols_double+0x12e>
mov    -0x18(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x30(%rbp),%ecx
add    -0x2c(%rbp),%ecx
movslq %ecx,%rcx
ucomisd (%rax,%rcx,8),%xmm0
jne    11b9 <_variable_cols_double+0xb9>
jp     11b9 <_variable_cols_double+0xb9>
jmp    121b <_variable_cols_double+0x11b>
mov    -0x10(%rbp),%rax
movslq -0x2c(%rbp),%rcx
cmpl   $0x0,(%rax,%rcx,4)
jne    1216 <_variable_cols_double+0x116>
mov    -0x10(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movl   $0x1,(%rax,%rcx,4)
mov    -0x34(%rbp),%eax
add    $0xffffffff,%eax
mov    %eax,-0x34(%rbp)
mov    -0x34(%rbp),%eax
cmp    -0x24(%rbp),%eax
jge    11fb <_variable_cols_double+0xfb>
movl   $0x0,-0x4(%rbp)
jmp    1248 <_variable_cols_double+0x148>
cmpl   $0x0,-0x34(%rbp)
jne    1211 <_variable_cols_double+0x111>
movl   $0x1,-0x4(%rbp)
jmp    1248 <_variable_cols_double+0x148>
jmp    1216 <_variable_cols_double+0x116>
jmp    121b <_variable_cols_double+0x11b>
jmp    1220 <_variable_cols_double+0x120>
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    117d <_variable_cols_double+0x7d>
jmp    1233 <_variable_cols_double+0x133>
mov    -0x28(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x28(%rbp)
jmp    1160 <_variable_cols_double+0x60>
movl   $0x1,-0x4(%rbp)
mov    -0x4(%rbp),%eax
pop    %rbp
ret

x86 assembly (clang, -O2):
<_variable_cols_double>:
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    1142 <_variable_cols_double+0x32>
lea    0x0(,%rbx,4),%rdx
mov    %r13,%rdi
xor    %esi,%esi
call   1030 <memset@plt>
mov    $0x1,%eax
test   %r12d,%r12d
jle    11b7 <_variable_cols_double+0xa7>
movslq %r14d,%r9
mov    %r12d,%r8d
shl    $0x3,%r9
xor    %esi,%esi
mov    %rbp,%rdi
mov    %r14d,%edx
jmp    116c <_variable_cols_double+0x5c>
add    $0x1,%rsi
add    %r9,%rdi
cmp    %r8,%rsi
je     11b7 <_variable_cols_double+0xa7>
test   %r14d,%r14d
jle    1160 <_variable_cols_double+0x50>
xor    %ecx,%ecx
jmp    1189 <_variable_cols_double+0x79>
cs nopw 0x0(%rax,%rax,1)
nop
add    $0x1,%rcx
cmp    %rcx,%rbx
je     1160 <_variable_cols_double+0x50>
movsd  0x0(%rbp,%rcx,8),%xmm0
ucomisd (%rdi,%rcx,8),%xmm0
jne    1198 <_variable_cols_double+0x88>
jnp    1180 <_variable_cols_double+0x70>
cmpl   $0x0,0x0(%r13,%rcx,4)
jne    1180 <_variable_cols_double+0x70>
movl   $0x1,0x0(%r13,%rcx,4)
cmp    %r15d,%edx
jle    11b5 <_variable_cols_double+0xa5>
add    $0xffffffff,%edx
jne    1180 <_variable_cols_double+0x70>
jmp    11b7 <_variable_cols_double+0xa7>
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(i32* %0, double* %1, i32 %2, i32 %3, i32 %4) {
%6 = alloca i32
%7 = alloca i32*
%8 = alloca double*
%9 = alloca i32
%10 = alloca i32
%11 = alloca i32
%12 = alloca i32
%13 = alloca i32
%14 = alloca i32
%15 = alloca i32
store i32* %0, i32** %7
store double* %1, double** %8
store i32 %2, i32* %9
store i32 %3, i32* %10
store i32 %4, i32* %11
%16 = load i32, i32* %9
%17 = sub nsw i32 %16, 1
%18 = load i32, i32* %10
%19 = mul nsw i32 %17, %18
store i32 %19, i32* %14
%20 = load i32, i32* %10
store i32 %20, i32* %15
store i32 0, i32* %13
br label %21
21: ; preds = %30, %5
%22 = load i32, i32* %13
%23 = load i32, i32* %10
%24 = icmp slt i32 %22, %23
br i1 %24, label %25, label %33
25: ; preds = %21
%26 = load i32*, i32** %7
%27 = load i32, i32* %13
%28 = sext i32 %27 to i64
%29 = getelementptr inbounds i32, i32* %26, i64 %28
store i32 0, i32* %29
br label %30
30: ; preds = %25
%31 = load i32, i32* %13
%32 = add nsw i32 %31, 1
store i32 %32, i32* %13
br label %21
33: ; preds = %21
store i32 0, i32* %12
br label %34
34: ; preds = %89, %33
%35 = load i32, i32* %12
%36 = load i32, i32* %9
%37 = icmp slt i32 %35, %36
br i1 %37, label %38, label %92
38: ; preds = %34
%39 = load i32, i32* %12
%40 = load i32, i32* %10
%41 = mul nsw i32 %39, %40
store i32 %41, i32* %14
store i32 0, i32* %13
br label %42
42: ; preds = %85, %38
%43 = load i32, i32* %13
%44 = load i32, i32* %10
%45 = icmp slt i32 %43, %44
br i1 %45, label %46, label %88
46: ; preds = %42
%47 = load double*, double** %8
%48 = load i32, i32* %13
%49 = sext i32 %48 to i64
%50 = getelementptr inbounds double, double* %47, i64 %49
%51 = load double, double* %50
%52 = load double*, double** %8
%53 = load i32, i32* %14
%54 = load i32, i32* %13
%55 = add nsw i32 %53, %54
%56 = sext i32 %55 to i64
%57 = getelementptr inbounds double, double* %52, i64 %56
%58 = load double, double* %57
%59 = fcmp une double %51, %58
br i1 %59, label %60, label %84
60: ; preds = %46
%61 = load i32*, i32** %7
%62 = load i32, i32* %13
%63 = sext i32 %62 to i64
%64 = getelementptr inbounds i32, i32* %61, i64 %63
%65 = load i32, i32* %64
%66 = icmp eq i32 %65, 0
br i1 %66, label %67, label %83
67: ; preds = %60
%68 = load i32*, i32** %7
%69 = load i32, i32* %13
%70 = sext i32 %69 to i64
%71 = getelementptr inbounds i32, i32* %68, i64 %70
store i32 1, i32* %71
%72 = load i32, i32* %15
%73 = add nsw i32 %72, -1
store i32 %73, i32* %15
%74 = load i32, i32* %15
%75 = load i32, i32* %11
%76 = icmp slt i32 %74, %75
br i1 %76, label %77, label %78
77: ; preds = %67
store i32 0, i32* %6
br label %93
78: ; preds = %67
%79 = load i32, i32* %15
%80 = icmp eq i32 %79, 0
br i1 %80, label %81, label %82
81: ; preds = %78
store i32 1, i32* %6
br label %93
82: ; preds = %78
br label %83
83: ; preds = %82, %60
br label %84
84: ; preds = %83, %46
br label %85
85: ; preds = %84
%86 = load i32, i32* %13
%87 = add nsw i32 %86, 1
store i32 %87, i32* %13
br label %42
88: ; preds = %42
br label %89
89: ; preds = %88
%90 = load i32, i32* %12
%91 = add nsw i32 %90, 1
store i32 %91, i32* %12
br label %34
92: ; preds = %34
store i32 1, i32* %6
br label %93
93: ; preds = %92, %81, %77
%94 = load i32, i32* %6
ret i32 %94
}

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