#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 *calloc(int, int);
int /*<<< orphan*/ equal_clusters(int const, int *, int *);
int /*<<< orphan*/ free(int *);
int *malloc(int const);
void emalg(int nclusters, int nrows, int ncolumns, double **data, int **mask,
           double weight[], int transpose,
           void getclustercenter(int, int, int, double **, int **, int[],
                                 double **, int **, int),
           double metric(int, double **, double **, int **, int **,
                         const double[], int, int, int),
           int clusterid[], double **cdata, int **cmask) {
  const int nitems = (transpose == 0) ? nrows : ncolumns;
  const int ndata = (transpose == 0) ? ncolumns : nrows;
  int *cn = calloc(nclusters, sizeof(int));
  int *savedids = malloc(nitems * sizeof(int));
  int changed;
  int iteration = 0;
  int period = 10;
  int i, j;
  for (i = 0; i < nitems; i++)
    cn[clusterid[i]]++;
  do {
    if (iteration % period == 0) {
      for (i = 0; i < nitems; i++)
        savedids[i] = clusterid[i];
      period *= 2;
    }
    iteration++;
    getclustercenter(nclusters, nrows, ncolumns, data, mask, clusterid, cdata,
                     cmask, transpose);
    changed = 0;
    for (i = 0; i < nitems; i++) {
      double distance;
      int jnow = clusterid[i];
      if (cn[jnow] == 1)
        continue;
      distance =
          metric(ndata, data, cdata, mask, cmask, weight, i, jnow, transpose);
      for (j = 0; j < nclusters; j++) {
        double tdistance;
        if (j == jnow)
          continue;
        tdistance =
            metric(ndata, data, cdata, mask, cmask, weight, i, j, transpose);
        if (tdistance < distance) {
          distance = tdistance;
          cn[clusterid[i]]--;
          clusterid[i] = j;
          cn[j]++;
          changed = 1;
        }
      }
    }
  } while (changed && !equal_clusters(nitems, savedids, clusterid));
  free(savedids);
  free(cn);
  return;
}


x86 assembly (clang, -O0):
<emalg>:
push   %rbp
mov    %rsp,%rbp
push   %rbx
sub    $0x98,%rsp
mov    0x38(%rbp),%rax
mov    0x30(%rbp),%rax
mov    0x28(%rbp),%rax
mov    0x20(%rbp),%rax
mov    0x18(%rbp),%rax
mov    0x10(%rbp),%eax
mov    %edi,-0xc(%rbp)
mov    %esi,-0x10(%rbp)
mov    %edx,-0x14(%rbp)
mov    %rcx,-0x20(%rbp)
mov    %r8,-0x28(%rbp)
mov    %r9,-0x30(%rbp)
cmpl   $0x0,0x10(%rbp)
jne    118d <emalg+0x4d>
mov    -0x10(%rbp),%eax
mov    %eax,-0x7c(%rbp)
jmp    1193 <emalg+0x53>
mov    -0x14(%rbp),%eax
mov    %eax,-0x7c(%rbp)
mov    -0x7c(%rbp),%eax
mov    %eax,-0x34(%rbp)
cmpl   $0x0,0x10(%rbp)
jne    11ae <emalg+0x6e>
mov    -0x14(%rbp),%eax
mov    %eax,-0x80(%rbp)
jmp    11b4 <emalg+0x74>
mov    -0x10(%rbp),%eax
mov    %eax,-0x80(%rbp)
mov    -0x80(%rbp),%eax
mov    %eax,-0x38(%rbp)
mov    -0xc(%rbp),%edi
mov    $0x4,%esi
call   1040 <calloc@plt>
mov    %rax,-0x40(%rbp)
movslq -0x34(%rbp),%rax
shl    $0x2,%rax
mov    %eax,%edi
call   1050 <malloc@plt>
mov    %rax,-0x48(%rbp)
movl   $0x0,-0x50(%rbp)
movl   $0xa,-0x54(%rbp)
movl   $0x0,-0x58(%rbp)
mov    -0x58(%rbp),%eax
cmp    -0x34(%rbp),%eax
jge    1226 <emalg+0xe6>
mov    -0x40(%rbp),%rax
mov    0x28(%rbp),%rcx
movslq -0x58(%rbp),%rdx
movslq (%rcx,%rdx,4),%rcx
mov    (%rax,%rcx,4),%edx
add    $0x1,%edx
mov    %edx,(%rax,%rcx,4)
mov    -0x58(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x58(%rbp)
jmp    11f3 <emalg+0xb3>
jmp    122b <emalg+0xeb>
mov    -0x50(%rbp),%eax
cltd
idivl  -0x54(%rbp)
cmp    $0x0,%edx
jne    127b <emalg+0x13b>
movl   $0x0,-0x58(%rbp)
mov    -0x58(%rbp),%eax
cmp    -0x34(%rbp),%eax
jge    1272 <emalg+0x132>
mov    0x28(%rbp),%rax
movslq -0x58(%rbp),%rcx
mov    (%rax,%rcx,4),%edx
mov    -0x48(%rbp),%rax
movslq -0x58(%rbp),%rcx
mov    %edx,(%rax,%rcx,4)
mov    -0x58(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x58(%rbp)
jmp    1242 <emalg+0x102>
mov    -0x54(%rbp),%eax
shl    $0x1,%eax
mov    %eax,-0x54(%rbp)
mov    -0x50(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x50(%rbp)
mov    0x18(%rbp),%rax
mov    -0xc(%rbp),%edi
mov    -0x10(%rbp),%esi
mov    -0x14(%rbp),%edx
mov    -0x20(%rbp),%rcx
mov    -0x28(%rbp),%r8
mov    0x28(%rbp),%r9
mov    0x30(%rbp),%rbx
mov    0x38(%rbp),%r11
mov    0x10(%rbp),%r10d
mov    %rbx,(%rsp)
mov    %r11,0x8(%rsp)
mov    %r10d,0x10(%rsp)
call   *%rax
movl   $0x0,-0x4c(%rbp)
movl   $0x0,-0x58(%rbp)
mov    -0x58(%rbp),%eax
cmp    -0x34(%rbp),%eax
jge    1413 <emalg+0x2d3>
mov    0x28(%rbp),%rax
movslq -0x58(%rbp),%rcx
mov    (%rax,%rcx,4),%eax
mov    %eax,-0x6c(%rbp)
mov    -0x40(%rbp),%rax
movslq -0x6c(%rbp),%rcx
cmpl   $0x1,(%rax,%rcx,4)
jne    12f8 <emalg+0x1b8>
jmp    1405 <emalg+0x2c5>
mov    0x20(%rbp),%rax
mov    -0x38(%rbp),%edi
mov    -0x20(%rbp),%rsi
mov    0x30(%rbp),%rdx
mov    -0x28(%rbp),%rcx
mov    0x38(%rbp),%r8
mov    -0x30(%rbp),%r9
mov    -0x58(%rbp),%ebx
mov    -0x6c(%rbp),%r11d
mov    0x10(%rbp),%r10d
mov    %ebx,(%rsp)
mov    %r11d,0x8(%rsp)
mov    %r10d,0x10(%rsp)
call   *%rax
movsd  %xmm0,-0x68(%rbp)
movl   $0x0,-0x5c(%rbp)
mov    -0x5c(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    1400 <emalg+0x2c0>
mov    -0x5c(%rbp),%eax
cmp    -0x6c(%rbp),%eax
jne    1356 <emalg+0x216>
jmp    13f2 <emalg+0x2b2>
mov    0x20(%rbp),%rax
mov    -0x38(%rbp),%edi
mov    -0x20(%rbp),%rsi
mov    0x30(%rbp),%rdx
mov    -0x28(%rbp),%rcx
mov    0x38(%rbp),%r8
mov    -0x30(%rbp),%r9
mov    -0x58(%rbp),%ebx
mov    -0x5c(%rbp),%r11d
mov    0x10(%rbp),%r10d
mov    %ebx,(%rsp)
mov    %r11d,0x8(%rsp)
mov    %r10d,0x10(%rsp)
call   *%rax
movsd  %xmm0,-0x78(%rbp)
movsd  -0x78(%rbp),%xmm1
movsd  -0x68(%rbp),%xmm0
ucomisd %xmm1,%xmm0
jbe    13ed <emalg+0x2ad>
movsd  -0x78(%rbp),%xmm0
movsd  %xmm0,-0x68(%rbp)
mov    -0x40(%rbp),%rax
mov    0x28(%rbp),%rcx
movslq -0x58(%rbp),%rdx
movslq (%rcx,%rdx,4),%rcx
mov    (%rax,%rcx,4),%edx
add    $0xffffffff,%edx
mov    %edx,(%rax,%rcx,4)
mov    -0x5c(%rbp),%edx
mov    0x28(%rbp),%rax
movslq -0x58(%rbp),%rcx
mov    %edx,(%rax,%rcx,4)
mov    -0x40(%rbp),%rax
movslq -0x5c(%rbp),%rcx
mov    (%rax,%rcx,4),%edx
add    $0x1,%edx
mov    %edx,(%rax,%rcx,4)
movl   $0x1,-0x4c(%rbp)
jmp    13f2 <emalg+0x2b2>
mov    -0x5c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x5c(%rbp)
jmp    1339 <emalg+0x1f9>
jmp    1405 <emalg+0x2c5>
mov    -0x58(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x58(%rbp)
jmp    12c7 <emalg+0x187>
jmp    1418 <emalg+0x2d8>
xor    %eax,%eax
cmpl   $0x0,-0x4c(%rbp)
mov    %al,-0x81(%rbp)
je     1448 <emalg+0x308>
mov    -0x34(%rbp),%edi
mov    -0x48(%rbp),%rsi
mov    0x28(%rbp),%rdx
call   1060 <equal_clusters@plt>
cmp    $0x0,%eax
setne  %al
xor    $0xff,%al
mov    %al,-0x81(%rbp)
mov    -0x81(%rbp),%al
test   $0x1,%al
jne    122b <emalg+0xeb>
mov    -0x48(%rbp),%rdi
call   1030 <free@plt>
mov    -0x40(%rbp),%rdi
call   1030 <free@plt>
add    $0x98,%rsp
pop    %rbx
pop    %rbp
ret

x86 assembly (clang, -O2):
<emalg>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x88,%rsp
mov    %r9,0x58(%rsp)
mov    %r8,0x40(%rsp)
mov    %rcx,0x38(%rsp)
mov    0xe8(%rsp),%rax
mov    %rax,0x30(%rsp)
mov    0xe0(%rsp),%rax
mov    %rax,0x28(%rsp)
mov    0xd8(%rsp),%rbp
mov    0xd0(%rsp),%rax
mov    %rax,0x50(%rsp)
mov    0xc8(%rsp),%rax
mov    %rax,0x70(%rsp)
mov    0xc0(%rsp),%eax
mov    %rax,0x20(%rsp)
test   %eax,%eax
mov    %edx,%ebx
cmove  %esi,%ebx
mov    %esi,0x48(%rsp)
mov    %edx,0x4c(%rsp)
cmove  %edx,%esi
mov    %esi,0x14(%rsp)
mov    %edi,0xc(%rsp)
mov    $0x4,%esi
call   1040 <calloc@plt>
mov    %rax,%r15
lea    0x0(,%rbx,4),%edi
call   1060 <malloc@plt>
mov    %rax,0x18(%rsp)
test   %ebx,%ebx
jle    1272 <emalg+0x122>
lea    -0x1(%rbx),%rcx
mov    %ebx,%eax
and    $0x3,%eax
cmp    $0x3,%rcx
jae    1207 <emalg+0xb7>
xor    %ecx,%ecx
jmp    1241 <emalg+0xf1>
mov    %ebx,%edx
and    $0xfffffffc,%edx
xor    %ecx,%ecx
xchg   %ax,%ax
movslq 0x0(%rbp,%rcx,4),%rsi
addl   $0x1,(%r15,%rsi,4)
movslq 0x4(%rbp,%rcx,4),%rsi
addl   $0x1,(%r15,%rsi,4)
movslq 0x8(%rbp,%rcx,4),%rsi
addl   $0x1,(%r15,%rsi,4)
movslq 0xc(%rbp,%rcx,4),%rsi
addl   $0x1,(%r15,%rsi,4)
add    $0x4,%rcx
cmp    %rcx,%rdx
jne    1210 <emalg+0xc0>
test   %rax,%rax
je     1272 <emalg+0x122>
lea    0x0(,%rcx,4),%rcx
add    %rbp,%rcx
xor    %edx,%edx
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
movslq (%rcx,%rdx,4),%rsi
addl   $0x1,(%r15,%rsi,4)
add    $0x1,%rdx
cmp    %rdx,%rax
jne    1260 <emalg+0x110>
lea    0x0(,%rbx,4),%rax
mov    %rax,0x68(%rsp)
mov    0xc(%rsp),%r12d
xor    %r14d,%r14d
movl   $0xa,0x8(%rsp)
mov    %rbx,0x80(%rsp)
nopw   0x0(%rax,%rax,1)
mov    %r14d,%eax
cltd
idivl  0x8(%rsp)
test   %edx,%edx
jne    12cc <emalg+0x17c>
test   %ebx,%ebx
jle    12c2 <emalg+0x172>
mov    0x18(%rsp),%rdi
mov    %rbp,%rsi
mov    0x68(%rsp),%rdx
call   1050 <memcpy@plt>
mov    0x8(%rsp),%eax
add    %eax,%eax
mov    %eax,0x8(%rsp)
sub    $0x8,%rsp
mov    0x14(%rsp),%edi
mov    0x50(%rsp),%esi
mov    0x54(%rsp),%edx
mov    0x40(%rsp),%rcx
mov    0x48(%rsp),%r8
mov    %rbp,%r9
push   0x28(%rsp)
push   0x40(%rsp)
push   0x40(%rsp)
call   *0x90(%rsp)
add    $0x20,%rsp
test   %ebx,%ebx
jle    1423 <emalg+0x2d3>
add    $0x1,%r14d
mov    %r14,0x78(%rsp)
xor    %r13d,%r13d
movl   $0x0,0x10(%rsp)
jmp    1335 <emalg+0x1e5>
xchg   %ax,%ax
add    $0x1,%r13
mov    0x80(%rsp),%rbx
cmp    %rbx,%r13
je     1400 <emalg+0x2b0>
movslq 0x0(%rbp,%r13,4),%rax
cmpl   $0x1,(%r15,%rax,4)
je     1320 <emalg+0x1d0>
mov    %eax,%r14d
sub    $0x8,%rsp
mov    0x1c(%rsp),%edi
mov    0x40(%rsp),%rsi
mov    0x30(%rsp),%rdx
mov    0x48(%rsp),%rcx
mov    0x38(%rsp),%r8
mov    0x60(%rsp),%r9
push   0x28(%rsp)
push   %r14
push   %r13
call   *0x70(%rsp)
add    $0x20,%rsp
movsd  %xmm0,0x60(%rsp)
cmpl   $0x0,0xc(%rsp)
jle    1320 <emalg+0x1d0>
xor    %ebx,%ebx
jmp    1399 <emalg+0x249>
cs nopw 0x0(%rax,%rax,1)
add    $0x1,%rbx
cmp    %rbx,%r12
je     1320 <emalg+0x1d0>
cmp    %rbx,%r14
je     1390 <emalg+0x240>
sub    $0x8,%rsp
mov    0x1c(%rsp),%edi
mov    0x40(%rsp),%rsi
mov    0x30(%rsp),%rdx
mov    0x48(%rsp),%rcx
mov    0x38(%rsp),%r8
mov    0x60(%rsp),%r9
push   0x28(%rsp)
push   %rbx
push   %r13
call   *0x70(%rsp)
add    $0x20,%rsp
movsd  0x60(%rsp),%xmm1
ucomisd %xmm0,%xmm1
jbe    1390 <emalg+0x240>
movslq 0x0(%rbp,%r13,4),%rax
addl   $0xffffffff,(%r15,%rax,4)
mov    %ebx,0x0(%rbp,%r13,4)
addl   $0x1,(%r15,%rbx,4)
movl   $0x1,0x10(%rsp)
movsd  %xmm0,0x60(%rsp)
jmp    1390 <emalg+0x240>
xchg   %ax,%ax
cmpl   $0x0,0x10(%rsp)
mov    0x78(%rsp),%r14
je     1423 <emalg+0x2d3>
mov    %ebx,%edi
mov    0x18(%rsp),%rsi
mov    %rbp,%rdx
call   1070 <equal_clusters@plt>
test   %eax,%eax
je     12a0 <emalg+0x150>
mov    0x18(%rsp),%rdi
call   1030 <free@plt>
mov    %r15,%rdi
add    $0x88,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
jmp    1030 <free@plt>

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

LLVM IR (-O2):
define void @emalg(i32 %0, i32 %1, i32 %2, double** %3, i32** %4, double* %5, i32 %6, void (i32, i32, i32, double**, i32**, i32*, double**, i32**, i32)* nocapture %7, double (i32, double**, double**, i32**, i32**, double*, i32, i32, i32)* nocapture %8, i32* %9, double** %10, i32** %11) {
%13 = bitcast i32* %9 to i8*
%14 = icmp eq i32 %6, 0
%15 = select i1 %14, i32 %1, i32 %2
%16 = select i1 %14, i32 %2, i32 %1
%17 = tail call i32* @calloc(i32 %0, i32 4)
%18 = shl i32 %15, 2
%19 = tail call i32* @malloc(i32 %18)
%20 = bitcast i32* %19 to i8*
%21 = icmp sgt i32 %15, 0
%22 = zext i32 %15 to i64
br i1 %21, label %23, label %44
23: ; preds = %12
%24 = add nsw i64 %22, -1
%25 = and i64 %22, 3
%26 = icmp ult i64 %24, 3
br i1 %26, label %29, label %27
27: ; preds = %23
%28 = and i64 %22, 4294967292
br label %48
29: ; preds = %48, %23
%30 = phi i64 [ 0, %23 ], [ %78, %48 ]
%31 = icmp eq i64 %25, 0
br i1 %31, label %44, label %32
32: ; preds = %29, %32
%33 = phi i64 [ %41, %32 ], [ %30, %29 ]
%34 = phi i64 [ %42, %32 ], [ 0, %29 ]
%35 = getelementptr inbounds i32, i32* %9, i64 %33
%36 = load i32, i32* %35
%37 = sext i32 %36 to i64
%38 = getelementptr inbounds i32, i32* %17, i64 %37
%39 = load i32, i32* %38
%40 = add nsw i32 %39, 1
store i32 %40, i32* %38
%41 = add nuw nsw i64 %33, 1
%42 = add i64 %34, 1
%43 = icmp eq i64 %42, %25
br i1 %43, label %44, label %32
44: ; preds = %29, %32, %12
%45 = icmp sgt i32 %0, 0
%46 = shl nuw nsw i64 %22, 2
%47 = zext i32 %0 to i64
br label %81
48: ; preds = %48, %27
%49 = phi i64 [ 0, %27 ], [ %78, %48 ]
%50 = phi i64 [ 0, %27 ], [ %79, %48 ]
%51 = getelementptr inbounds i32, i32* %9, i64 %49
%52 = load i32, i32* %51
%53 = sext i32 %52 to i64
%54 = getelementptr inbounds i32, i32* %17, i64 %53
%55 = load i32, i32* %54
%56 = add nsw i32 %55, 1
store i32 %56, i32* %54
%57 = or i64 %49, 1
%58 = getelementptr inbounds i32, i32* %9, i64 %57
%59 = load i32, i32* %58
%60 = sext i32 %59 to i64
%61 = getelementptr inbounds i32, i32* %17, i64 %60
%62 = load i32, i32* %61
%63 = add nsw i32 %62, 1
store i32 %63, i32* %61
%64 = or i64 %49, 2
%65 = getelementptr inbounds i32, i32* %9, i64 %64
%66 = load i32, i32* %65
%67 = sext i32 %66 to i64
%68 = getelementptr inbounds i32, i32* %17, i64 %67
%69 = load i32, i32* %68
%70 = add nsw i32 %69, 1
store i32 %70, i32* %68
%71 = or i64 %49, 3
%72 = getelementptr inbounds i32, i32* %9, i64 %71
%73 = load i32, i32* %72
%74 = sext i32 %73 to i64
%75 = getelementptr inbounds i32, i32* %17, i64 %74
%76 = load i32, i32* %75
%77 = add nsw i32 %76, 1
store i32 %77, i32* %75
%78 = add nuw nsw i64 %49, 4
%79 = add i64 %50, 4
%80 = icmp eq i64 %79, %28
br i1 %80, label %29, label %48
81: ; preds = %44, %136
%82 = phi i32 [ %91, %136 ], [ 10, %44 ]
%83 = phi i32 [ %92, %136 ], [ 0, %44 ]
%84 = srem i32 %83, %82
%85 = icmp eq i32 %84, 0
br i1 %85, label %86, label %90
86: ; preds = %81
br i1 %21, label %87, label %88
87: ; preds = %86
br label %88
88: ; preds = %87, %86
%89 = shl nsw i32 %82, 1
br label %90
90: ; preds = %88, %81
%91 = phi i32 [ %89, %88 ], [ %82, %81 ]
%92 = add nuw nsw i32 %83, 1
tail call void %7(i32 %0, i32 %1, i32 %2, double** %3, i32** %4, i32* %9, double** %10, i32** %11, i32 %6)
br i1 %21, label %93, label %139
93: ; preds = %90, %130
%94 = phi i64 [ %132, %130 ], [ 0, %90 ]
%95 = phi i32 [ %131, %130 ], [ 0, %90 ]
%96 = getelementptr inbounds i32, i32* %9, i64 %94
%97 = load i32, i32* %96
%98 = sext i32 %97 to i64
%99 = getelementptr inbounds i32, i32* %17, i64 %98
%100 = load i32, i32* %99
%101 = icmp eq i32 %100, 1
br i1 %101, label %130, label %102
102: ; preds = %93
%103 = trunc i64 %94 to i32
%104 = tail call double %8(i32 %16, double** %3, double** %10, i32** %4, i32** %11, double* %5, i32 %103, i32 %97, i32 %6)
br i1 %45, label %105, label %130
105: ; preds = %102
%106 = zext i32 %97 to i64
br label %107
107: ; preds = %105, %125
%108 = phi i64 [ 0, %105 ], [ %128, %125 ]
%109 = phi i32 [ %95, %105 ], [ %127, %125 ]
%110 = phi double [ %104, %105 ], [ %126, %125 ]
%111 = icmp eq i64 %108, %106
br i1 %111, label %125, label %112
112: ; preds = %107
%113 = trunc i64 %108 to i32
%114 = tail call double %8(i32 %16, double** %3, double** %10, i32** %4, i32** %11, double* %5, i32 %103, i32 %113, i32 %6)
%115 = fcmp olt double %114, %110
br i1 %115, label %116, label %125
116: ; preds = %112
%117 = load i32, i32* %96
%118 = sext i32 %117 to i64
%119 = getelementptr inbounds i32, i32* %17, i64 %118
%120 = load i32, i32* %119
%121 = add nsw i32 %120, -1
store i32 %121, i32* %119
store i32 %113, i32* %96
%122 = getelementptr inbounds i32, i32* %17, i64 %108
%123 = load i32, i32* %122
%124 = add nsw i32 %123, 1
store i32 %124, i32* %122
br label %125
125: ; preds = %112, %116, %107
%126 = phi double [ %110, %107 ], [ %114, %116 ], [ %110, %112 ]
%127 = phi i32 [ %109, %107 ], [ 1, %116 ], [ %109, %112 ]
%128 = add nuw nsw i64 %108, 1
%129 = icmp eq i64 %128, %47
br i1 %129, label %130, label %107
130: ; preds = %125, %102, %93
%131 = phi i32 [ %95, %93 ], [ %95, %102 ], [ %127, %125 ]
%132 = add nuw nsw i64 %94, 1
%133 = icmp eq i64 %132, %22
br i1 %133, label %134, label %93
134: ; preds = %130
%135 = icmp eq i32 %131, 0
br i1 %135, label %139, label %136
136: ; preds = %134
%137 = tail call i32 @equal_clusters(i32 %15, i32* %19, i32* %9)
%138 = icmp eq i32 %137, 0
br i1 %138, label %81, label %139
139: ; preds = %90, %134, %136
%140 = tail call i32 @free(i32* %19)
%141 = tail call i32 @free(i32* %17)
ret void
}