void get_max_Rfield_for_each_cluster(const double *Z, const double *R,
                                     double *max_rfs, int n, int rf) {
  int *curNode;
  int ndid, lid, rid, k;
  unsigned char *lvisited, *rvisited;
  const double *Zrow, *Rrow;
  double max_rf;
  const int bff = (((((double)(n)) / (double)(sizeof(unsigned char) * 8)) ==
                            ((double)(((n)) / ((sizeof(unsigned char) * 8))))
                        ? (((n)) / ((sizeof(unsigned char) * 8)))
                        : (((n)) / ((sizeof(unsigned char) * 8)) + 1)));
  k = 0;
  curNode = (int *)malloc(n * sizeof(int));
  lvisited = (unsigned char *)malloc(bff);
  rvisited = (unsigned char *)malloc(bff);
  curNode[k] = (n * 2) - 2;
  bzero(lvisited, bff);
  bzero(rvisited, bff);
  while (k >= 0) {
    ndid = curNode[k];
    Zrow = Z + ((ndid - n) * 4);
    Rrow = R + ((ndid - n) * 4);
    lid = (int)Zrow[0];
    rid = (int)Zrow[1];
    if (lid >= n && !(((lvisited)[(ndid - n) / (sizeof(unsigned char) * 8)] >>
                       (((sizeof(unsigned char) * 8) - 1) -
                        ((ndid - n) % (sizeof(unsigned char) * 8)))) &
                      0x1)) {
      ((lvisited)[(ndid - n) / (sizeof(unsigned char) * 8)] |=
       ((0x1) << (((sizeof(unsigned char) * 8) - 1) -
                  ((ndid - n) % (sizeof(unsigned char) * 8)))));
      curNode[k + 1] = lid;
      k++;
      continue;
    }
    if (rid >= n && !(((rvisited)[(ndid - n) / (sizeof(unsigned char) * 8)] >>
                       (((sizeof(unsigned char) * 8) - 1) -
                        ((ndid - n) % (sizeof(unsigned char) * 8)))) &
                      0x1)) {
      ((rvisited)[(ndid - n) / (sizeof(unsigned char) * 8)] |=
       ((0x1) << (((sizeof(unsigned char) * 8) - 1) -
                  ((ndid - n) % (sizeof(unsigned char) * 8)))));
      curNode[k + 1] = rid;
      k++;
      continue;
    }
    max_rf = Rrow[rf];
    if (lid >= n) {
      max_rf = ((max_rf > max_rfs[lid - n]) ? (max_rf) : (max_rfs[lid - n]));
    }
    if (rid >= n) {
      max_rf = ((max_rf > max_rfs[rid - n]) ? (max_rf) : (max_rfs[rid - n]));
    }
    max_rfs[ndid - n] = max_rf;
    k--;
  }
  free(curNode);
  free(lvisited);
  free(rvisited);
}


x86 assembly (clang, -O0):
<get_max_Rfield_for_each_cluster>:
push   %rbp
mov    %rsp,%rbp
sub    $0x80,%rsp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %ecx,-0x1c(%rbp)
mov    %r8d,-0x20(%rbp)
movslq -0x1c(%rbp),%rax
mov    %eax,%ecx
cvtsi2sd %ecx,%xmm0
movsd  0xea0(%rip),%xmm1
divsd  %xmm1,%xmm0
shr    $0x3,%rax
cvtsi2sd %rax,%xmm1
ucomisd %xmm1,%xmm0
jne    118e <get_max_Rfield_for_each_cluster+0x5e>
jp     118e <get_max_Rfield_for_each_cluster+0x5e>
movslq -0x1c(%rbp),%rax
shr    $0x3,%rax
mov    %rax,-0x70(%rbp)
jmp    119e <get_max_Rfield_for_each_cluster+0x6e>
movslq -0x1c(%rbp),%rax
shr    $0x3,%rax
add    $0x1,%rax
mov    %rax,-0x70(%rbp)
mov    -0x70(%rbp),%rax
mov    %eax,-0x64(%rbp)
movl   $0x0,-0x38(%rbp)
movslq -0x1c(%rbp),%rdi
shl    $0x2,%rdi
call   1050 <malloc@plt>
mov    %rax,-0x28(%rbp)
movslq -0x64(%rbp),%rdi
call   1050 <malloc@plt>
mov    %rax,-0x40(%rbp)
movslq -0x64(%rbp),%rdi
call   1050 <malloc@plt>
mov    %rax,-0x48(%rbp)
mov    -0x1c(%rbp),%edx
shl    $0x1,%edx
sub    $0x2,%edx
mov    -0x28(%rbp),%rax
movslq -0x38(%rbp),%rcx
mov    %edx,(%rax,%rcx,4)
mov    -0x40(%rbp),%rdi
movslq -0x64(%rbp),%rdx
xor    %esi,%esi
call   1040 <memset@plt>
mov    -0x48(%rbp),%rdi
movslq -0x64(%rbp),%rdx
xor    %esi,%esi
call   1040 <memset@plt>
cmpl   $0x0,-0x38(%rbp)
jl     1491 <get_max_Rfield_for_each_cluster+0x361>
mov    -0x28(%rbp),%rax
movslq -0x38(%rbp),%rcx
mov    (%rax,%rcx,4),%eax
mov    %eax,-0x2c(%rbp)
mov    -0x8(%rbp),%rax
mov    -0x2c(%rbp),%ecx
sub    -0x1c(%rbp),%ecx
shl    $0x2,%ecx
movslq %ecx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    %rax,-0x50(%rbp)
mov    -0x10(%rbp),%rax
mov    -0x2c(%rbp),%ecx
sub    -0x1c(%rbp),%ecx
shl    $0x2,%ecx
movslq %ecx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    %rax,-0x58(%rbp)
mov    -0x50(%rbp),%rax
cvttsd2si (%rax),%eax
mov    %eax,-0x30(%rbp)
mov    -0x50(%rbp),%rax
cvttsd2si 0x8(%rax),%eax
mov    %eax,-0x34(%rbp)
mov    -0x30(%rbp),%eax
cmp    -0x1c(%rbp),%eax
jl     130b <get_max_Rfield_for_each_cluster+0x1db>
mov    -0x40(%rbp),%rax
mov    -0x2c(%rbp),%ecx
sub    -0x1c(%rbp),%ecx
movslq %ecx,%rcx
shr    $0x3,%rcx
movzbl (%rax,%rcx,1),%eax
mov    -0x2c(%rbp),%ecx
sub    -0x1c(%rbp),%ecx
movslq %ecx,%rdx
and    $0x7,%rdx
mov    $0x7,%ecx
sub    %rdx,%rcx
sar    %cl,%eax
and    $0x1,%eax
cmp    $0x0,%eax
jne    130b <get_max_Rfield_for_each_cluster+0x1db>
mov    -0x2c(%rbp),%eax
sub    -0x1c(%rbp),%eax
movslq %eax,%rcx
and    $0x7,%rcx
mov    $0x7,%eax
sub    %rcx,%rax
mov    %eax,%ecx
mov    $0x1,%esi
shl    %cl,%esi
mov    -0x40(%rbp),%rax
mov    -0x2c(%rbp),%ecx
sub    -0x1c(%rbp),%ecx
movslq %ecx,%rcx
shr    $0x3,%rcx
movzbl (%rax,%rcx,1),%edx
or     %esi,%edx
mov    %dl,(%rax,%rcx,1)
mov    -0x30(%rbp),%edx
mov    -0x28(%rbp),%rax
mov    -0x38(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
mov    %edx,(%rax,%rcx,4)
mov    -0x38(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x38(%rbp)
jmp    1209 <get_max_Rfield_for_each_cluster+0xd9>
mov    -0x34(%rbp),%eax
cmp    -0x1c(%rbp),%eax
jl     13a8 <get_max_Rfield_for_each_cluster+0x278>
mov    -0x48(%rbp),%rax
mov    -0x2c(%rbp),%ecx
sub    -0x1c(%rbp),%ecx
movslq %ecx,%rcx
shr    $0x3,%rcx
movzbl (%rax,%rcx,1),%eax
mov    -0x2c(%rbp),%ecx
sub    -0x1c(%rbp),%ecx
movslq %ecx,%rdx
and    $0x7,%rdx
mov    $0x7,%ecx
sub    %rdx,%rcx
sar    %cl,%eax
and    $0x1,%eax
cmp    $0x0,%eax
jne    13a8 <get_max_Rfield_for_each_cluster+0x278>
mov    -0x2c(%rbp),%eax
sub    -0x1c(%rbp),%eax
movslq %eax,%rcx
and    $0x7,%rcx
mov    $0x7,%eax
sub    %rcx,%rax
mov    %eax,%ecx
mov    $0x1,%esi
shl    %cl,%esi
mov    -0x48(%rbp),%rax
mov    -0x2c(%rbp),%ecx
sub    -0x1c(%rbp),%ecx
movslq %ecx,%rcx
shr    $0x3,%rcx
movzbl (%rax,%rcx,1),%edx
or     %esi,%edx
mov    %dl,(%rax,%rcx,1)
mov    -0x34(%rbp),%edx
mov    -0x28(%rbp),%rax
mov    -0x38(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
mov    %edx,(%rax,%rcx,4)
mov    -0x38(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x38(%rbp)
jmp    1209 <get_max_Rfield_for_each_cluster+0xd9>
mov    -0x58(%rbp),%rax
movslq -0x20(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x60(%rbp)
mov    -0x30(%rbp),%eax
cmp    -0x1c(%rbp),%eax
jl     1413 <get_max_Rfield_for_each_cluster+0x2e3>
movsd  -0x60(%rbp),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x30(%rbp),%ecx
sub    -0x1c(%rbp),%ecx
movslq %ecx,%rcx
ucomisd (%rax,%rcx,8),%xmm0
jbe    13f2 <get_max_Rfield_for_each_cluster+0x2c2>
movsd  -0x60(%rbp),%xmm0
movsd  %xmm0,-0x78(%rbp)
jmp    1409 <get_max_Rfield_for_each_cluster+0x2d9>
mov    -0x18(%rbp),%rax
mov    -0x30(%rbp),%ecx
sub    -0x1c(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x78(%rbp)
movsd  -0x78(%rbp),%xmm0
movsd  %xmm0,-0x60(%rbp)
mov    -0x34(%rbp),%eax
cmp    -0x1c(%rbp),%eax
jl     146c <get_max_Rfield_for_each_cluster+0x33c>
movsd  -0x60(%rbp),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x34(%rbp),%ecx
sub    -0x1c(%rbp),%ecx
movslq %ecx,%rcx
ucomisd (%rax,%rcx,8),%xmm0
jbe    144b <get_max_Rfield_for_each_cluster+0x31b>
movsd  -0x60(%rbp),%xmm0
movsd  %xmm0,-0x80(%rbp)
jmp    1462 <get_max_Rfield_for_each_cluster+0x332>
mov    -0x18(%rbp),%rax
mov    -0x34(%rbp),%ecx
sub    -0x1c(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x80(%rbp)
movsd  -0x80(%rbp),%xmm0
movsd  %xmm0,-0x60(%rbp)
movsd  -0x60(%rbp),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x2c(%rbp),%ecx
sub    -0x1c(%rbp),%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x38(%rbp),%eax
add    $0xffffffff,%eax
mov    %eax,-0x38(%rbp)
jmp    1209 <get_max_Rfield_for_each_cluster+0xd9>
mov    -0x28(%rbp),%rdi
call   1030 <free@plt>
mov    -0x40(%rbp),%rdi
call   1030 <free@plt>
mov    -0x48(%rbp),%rdi
call   1030 <free@plt>
add    $0x80,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<get_max_Rfield_for_each_cluster>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x28,%rsp
mov    %r8d,0x14(%rsp)
mov    %ecx,%r13d
mov    %rdx,%r14
mov    %rsi,0x20(%rsp)
mov    %rdi,%r15
cvtsi2sd %ecx,%xmm0
mulsd  0xea3(%rip),%xmm0
movslq %ecx,%rbx
mov    %rbx,%rbp
shr    $0x3,%rbp
cvtsi2sd %rbp,%xmm1
cmpneqsd %xmm0,%xmm1
movq   %xmm1,%rax
sub    %eax,%ebp
lea    0x0(,%rbx,4),%rdi
call   1050 <malloc@plt>
mov    %rax,%r12
movslq %ebp,%rbp
mov    $0x1,%edi
mov    %rbp,%rsi
call   1040 <calloc@plt>
mov    %rax,0x18(%rsp)
mov    $0x1,%edi
mov    %rbp,%rsi
call   1040 <calloc@plt>
mov    0x18(%rsp),%r11
mov    %rax,0x8(%rsp)
lea    (%rbx,%rbx,1),%ecx
add    $0xfffffffe,%ecx
mov    %ecx,(%r12)
movslq 0x14(%rsp),%rax
mov    0x20(%rsp),%rdx
lea    (%rdx,%rax,8),%r8
xor    %r9d,%r9d
jmp    11ec <get_max_Rfield_for_each_cluster+0xbc>
add    %r11,%rdi
or     %bpl,%dl
mov    %dl,(%rdi)
mov    %r9d,%ecx
add    $0x1,%r9d
mov    %eax,0x4(%r12,%rcx,4)
mov    %r9d,%eax
mov    (%r12,%rax,4),%ecx
sub    %r13d,%ecx
lea    0x0(,%rcx,4),%eax
movslq %eax,%rbx
cvttsd2si (%r15,%rbx,8),%eax
movslq %ecx,%r10
cmp    %r13d,%eax
jl     122a <get_max_Rfield_for_each_cluster+0xfa>
mov    %r10,%rdi
mov    %r10d,%ecx
not    %cl
and    $0x7,%cl
mov    $0x1,%edx
shl    %cl,%edx
shr    $0x3,%rdi
movzbl (%r11,%rdi,1),%ebp
movzbl %cl,%ecx
bt     %ecx,%ebp
jae    11d1 <get_max_Rfield_for_each_cluster+0xa1>
cvttsd2si 0x8(%r15,%rbx,8),%esi
cmp    %r13d,%esi
jl     125e <get_max_Rfield_for_each_cluster+0x12e>
mov    %r10,%rdi
mov    %r10d,%ecx
not    %cl
and    $0x7,%cl
mov    $0x1,%edx
shl    %cl,%edx
shr    $0x3,%rdi
mov    0x8(%rsp),%rbp
movzbl 0x0(%rbp,%rdi,1),%ebp
movzbl %cl,%ecx
bt     %ecx,%ebp
jae    1293 <get_max_Rfield_for_each_cluster+0x163>
movsd  (%r8,%rbx,8),%xmm0
sub    %r13d,%eax
jl     1271 <get_max_Rfield_for_each_cluster+0x141>
cltq
maxsd  (%r14,%rax,8),%xmm0
sub    %r13d,%esi
jl     127f <get_max_Rfield_for_each_cluster+0x14f>
movslq %esi,%rax
maxsd  (%r14,%rax,8),%xmm0
movsd  %xmm0,(%r14,%r10,8)
test   %r9d,%r9d
jle    129f <get_max_Rfield_for_each_cluster+0x16f>
add    $0xffffffff,%r9d
jmp    11e5 <get_max_Rfield_for_each_cluster+0xb5>
add    0x8(%rsp),%rdi
mov    %esi,%eax
jmp    11d4 <get_max_Rfield_for_each_cluster+0xa4>
mov    %r12,%rdi
mov    %r11,%rbx
call   1030 <free@plt>
mov    %rbx,%rdi
call   1030 <free@plt>
mov    0x8(%rsp),%rdi
add    $0x28,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
jmp    1030 <free@plt>

LLVM IR (-O0):
define void @get_max_Rfield_for_each_cluster(double* %0, double* %1, double* %2, i32 %3, i32 %4) {
%6 = alloca double*
%7 = alloca double*
%8 = alloca double*
%9 = alloca i32
%10 = alloca i32
%11 = alloca i32*
%12 = alloca i32
%13 = alloca i32
%14 = alloca i32
%15 = alloca i32
%16 = alloca i8*
%17 = alloca i8*
%18 = alloca double*
%19 = alloca double*
%20 = alloca double
%21 = alloca i32
store double* %0, double** %6
store double* %1, double** %7
store double* %2, double** %8
store i32 %3, i32* %9
store i32 %4, i32* %10
%22 = load i32, i32* %9
%23 = sitofp i32 %22 to double
%24 = fdiv double %23, 8.000000e+00
%25 = load i32, i32* %9
%26 = sext i32 %25 to i64
%27 = udiv i64 %26, 8
%28 = uitofp i64 %27 to double
%29 = fcmp oeq double %24, %28
br i1 %29, label %30, label %34
30: ; preds = %5
%31 = load i32, i32* %9
%32 = sext i32 %31 to i64
%33 = udiv i64 %32, 8
br label %39
34: ; preds = %5
%35 = load i32, i32* %9
%36 = sext i32 %35 to i64
%37 = udiv i64 %36, 8
%38 = add i64 %37, 1
br label %39
39: ; preds = %34, %30
%40 = phi i64 [ %33, %30 ], [ %38, %34 ]
%41 = trunc i64 %40 to i32
store i32 %41, i32* %21
store i32 0, i32* %15
%42 = load i32, i32* %9
%43 = sext i32 %42 to i64
%44 = mul i64 %43, 4
%45 = call i8* @malloc(i64 %44)
%46 = bitcast i8* %45 to i32*
store i32* %46, i32** %11
%47 = load i32, i32* %21
%48 = sext i32 %47 to i64
%49 = call i8* @malloc(i64 %48)
store i8* %49, i8** %16
%50 = load i32, i32* %21
%51 = sext i32 %50 to i64
%52 = call i8* @malloc(i64 %51)
store i8* %52, i8** %17
%53 = load i32, i32* %9
%54 = mul nsw i32 %53, 2
%55 = sub nsw i32 %54, 2
%56 = load i32*, i32** %11
%57 = load i32, i32* %15
%58 = sext i32 %57 to i64
%59 = getelementptr inbounds i32, i32* %56, i64 %58
store i32 %55, i32* %59
%60 = load i8*, i8** %16
%61 = load i32, i32* %21
%62 = sext i32 %61 to i64
%63 = load i8*, i8** %17
%64 = load i32, i32* %21
%65 = sext i32 %64 to i64
br label %66
66: ; preds = %257, %172, %120, %39
%67 = load i32, i32* %15
%68 = icmp sge i32 %67, 0
br i1 %68, label %69, label %267
69: ; preds = %66
%70 = load i32*, i32** %11
%71 = load i32, i32* %15
%72 = sext i32 %71 to i64
%73 = getelementptr inbounds i32, i32* %70, i64 %72
%74 = load i32, i32* %73
store i32 %74, i32* %12
%75 = load double*, double** %6
%76 = load i32, i32* %12
%77 = load i32, i32* %9
%78 = sub nsw i32 %76, %77
%79 = mul nsw i32 %78, 4
%80 = sext i32 %79 to i64
%81 = getelementptr inbounds double, double* %75, i64 %80
store double* %81, double** %18
%82 = load double*, double** %7
%83 = load i32, i32* %12
%84 = load i32, i32* %9
%85 = sub nsw i32 %83, %84
%86 = mul nsw i32 %85, 4
%87 = sext i32 %86 to i64
%88 = getelementptr inbounds double, double* %82, i64 %87
store double* %88, double** %19
%89 = load double*, double** %18
%90 = getelementptr inbounds double, double* %89, i64 0
%91 = load double, double* %90
%92 = fptosi double %91 to i32
store i32 %92, i32* %13
%93 = load double*, double** %18
%94 = getelementptr inbounds double, double* %93, i64 1
%95 = load double, double* %94
%96 = fptosi double %95 to i32
store i32 %96, i32* %14
%97 = load i32, i32* %13
%98 = load i32, i32* %9
%99 = icmp sge i32 %97, %98
br i1 %99, label %100, label %148
100: ; preds = %69
%101 = load i8*, i8** %16
%102 = load i32, i32* %12
%103 = load i32, i32* %9
%104 = sub nsw i32 %102, %103
%105 = sext i32 %104 to i64
%106 = udiv i64 %105, 8
%107 = getelementptr inbounds i8, i8* %101, i64 %106
%108 = load i8, i8* %107
%109 = zext i8 %108 to i32
%110 = load i32, i32* %12
%111 = load i32, i32* %9
%112 = sub nsw i32 %110, %111
%113 = sext i32 %112 to i64
%114 = urem i64 %113, 8
%115 = sub i64 7, %114
%116 = trunc i64 %115 to i32
%117 = ashr i32 %109, %116
%118 = and i32 %117, 1
%119 = icmp ne i32 %118, 0
br i1 %119, label %148, label %120
120: ; preds = %100
%121 = load i32, i32* %12
%122 = load i32, i32* %9
%123 = sub nsw i32 %121, %122
%124 = sext i32 %123 to i64
%125 = urem i64 %124, 8
%126 = sub i64 7, %125
%127 = trunc i64 %126 to i32
%128 = shl i32 1, %127
%129 = load i8*, i8** %16
%130 = load i32, i32* %12
%131 = load i32, i32* %9
%132 = sub nsw i32 %130, %131
%133 = sext i32 %132 to i64
%134 = udiv i64 %133, 8
%135 = getelementptr inbounds i8, i8* %129, i64 %134
%136 = load i8, i8* %135
%137 = zext i8 %136 to i32
%138 = or i32 %137, %128
%139 = trunc i32 %138 to i8
store i8 %139, i8* %135
%140 = load i32, i32* %13
%141 = load i32*, i32** %11
%142 = load i32, i32* %15
%143 = add nsw i32 %142, 1
%144 = sext i32 %143 to i64
%145 = getelementptr inbounds i32, i32* %141, i64 %144
store i32 %140, i32* %145
%146 = load i32, i32* %15
%147 = add nsw i32 %146, 1
store i32 %147, i32* %15
br label %66
148: ; preds = %100, %69
%149 = load i32, i32* %14
%150 = load i32, i32* %9
%151 = icmp sge i32 %149, %150
br i1 %151, label %152, label %200
152: ; preds = %148
%153 = load i8*, i8** %17
%154 = load i32, i32* %12
%155 = load i32, i32* %9
%156 = sub nsw i32 %154, %155
%157 = sext i32 %156 to i64
%158 = udiv i64 %157, 8
%159 = getelementptr inbounds i8, i8* %153, i64 %158
%160 = load i8, i8* %159
%161 = zext i8 %160 to i32
%162 = load i32, i32* %12
%163 = load i32, i32* %9
%164 = sub nsw i32 %162, %163
%165 = sext i32 %164 to i64
%166 = urem i64 %165, 8
%167 = sub i64 7, %166
%168 = trunc i64 %167 to i32
%169 = ashr i32 %161, %168
%170 = and i32 %169, 1
%171 = icmp ne i32 %170, 0
br i1 %171, label %200, label %172
172: ; preds = %152
%173 = load i32, i32* %12
%174 = load i32, i32* %9
%175 = sub nsw i32 %173, %174
%176 = sext i32 %175 to i64
%177 = urem i64 %176, 8
%178 = sub i64 7, %177
%179 = trunc i64 %178 to i32
%180 = shl i32 1, %179
%181 = load i8*, i8** %17
%182 = load i32, i32* %12
%183 = load i32, i32* %9
%184 = sub nsw i32 %182, %183
%185 = sext i32 %184 to i64
%186 = udiv i64 %185, 8
%187 = getelementptr inbounds i8, i8* %181, i64 %186
%188 = load i8, i8* %187
%189 = zext i8 %188 to i32
%190 = or i32 %189, %180
%191 = trunc i32 %190 to i8
store i8 %191, i8* %187
%192 = load i32, i32* %14
%193 = load i32*, i32** %11
%194 = load i32, i32* %15
%195 = add nsw i32 %194, 1
%196 = sext i32 %195 to i64
%197 = getelementptr inbounds i32, i32* %193, i64 %196
store i32 %192, i32* %197
%198 = load i32, i32* %15
%199 = add nsw i32 %198, 1
store i32 %199, i32* %15
br label %66
200: ; preds = %152, %148
%201 = load double*, double** %19
%202 = load i32, i32* %10
%203 = sext i32 %202 to i64
%204 = getelementptr inbounds double, double* %201, i64 %203
%205 = load double, double* %204
store double %205, double* %20
%206 = load i32, i32* %13
%207 = load i32, i32* %9
%208 = icmp sge i32 %206, %207
br i1 %208, label %209, label %231
209: ; preds = %200
%210 = load double, double* %20
%211 = load double*, double** %8
%212 = load i32, i32* %13
%213 = load i32, i32* %9
%214 = sub nsw i32 %212, %213
%215 = sext i32 %214 to i64
%216 = getelementptr inbounds double, double* %211, i64 %215
%217 = load double, double* %216
%218 = fcmp ogt double %210, %217
br i1 %218, label %219, label %221
219: ; preds = %209
%220 = load double, double* %20
br label %229
221: ; preds = %209
%222 = load double*, double** %8
%223 = load i32, i32* %13
%224 = load i32, i32* %9
%225 = sub nsw i32 %223, %224
%226 = sext i32 %225 to i64
%227 = getelementptr inbounds double, double* %222, i64 %226
%228 = load double, double* %227
br label %229
229: ; preds = %221, %219
%230 = phi double [ %220, %219 ], [ %228, %221 ]
store double %230, double* %20
br label %231
231: ; preds = %229, %200
%232 = load i32, i32* %14
%233 = load i32, i32* %9
%234 = icmp sge i32 %232, %233
br i1 %234, label %235, label %257
235: ; preds = %231
%236 = load double, double* %20
%237 = load double*, double** %8
%238 = load i32, i32* %14
%239 = load i32, i32* %9
%240 = sub nsw i32 %238, %239
%241 = sext i32 %240 to i64
%242 = getelementptr inbounds double, double* %237, i64 %241
%243 = load double, double* %242
%244 = fcmp ogt double %236, %243
br i1 %244, label %245, label %247
245: ; preds = %235
%246 = load double, double* %20
br label %255
247: ; preds = %235
%248 = load double*, double** %8
%249 = load i32, i32* %14
%250 = load i32, i32* %9
%251 = sub nsw i32 %249, %250
%252 = sext i32 %251 to i64
%253 = getelementptr inbounds double, double* %248, i64 %252
%254 = load double, double* %253
br label %255
255: ; preds = %247, %245
%256 = phi double [ %246, %245 ], [ %254, %247 ]
store double %256, double* %20
br label %257
257: ; preds = %255, %231
%258 = load double, double* %20
%259 = load double*, double** %8
%260 = load i32, i32* %12
%261 = load i32, i32* %9
%262 = sub nsw i32 %260, %261
%263 = sext i32 %262 to i64
%264 = getelementptr inbounds double, double* %259, i64 %263
store double %258, double* %264
%265 = load i32, i32* %15
%266 = add nsw i32 %265, -1
store i32 %266, i32* %15
br label %66
267: ; preds = %66
%268 = load i32*, i32** %11
%269 = bitcast i32* %268 to i8*
%270 = load i8*, i8** %16
%271 = load i8*, i8** %17
ret void
}

LLVM IR (-O2):
define void @get_max_Rfield_for_each_cluster(double* nocapture %0, double* nocapture %1, double* nocapture %2, i32 %3, i32 %4) {
%6 = sitofp i32 %3 to double
%7 = fmul double %6, 1.250000e-01
%8 = sext i32 %3 to i64
%9 = lshr i64 %8, 3
%10 = uitofp i64 %9 to double
%11 = fcmp une double %7, %10
%12 = zext i1 %11 to i64
%13 = add nuw nsw i64 %9, %12
%14 = shl nsw i64 %8, 2
%15 = tail call i8* @malloc(i64 %14)
%16 = bitcast i8* %15 to i32*
%17 = shl i64 %13, 32
%18 = ashr exact i64 %17, 32
%19 = call i8* @calloc(i64 1, i64 %18)
%20 = call i8* @calloc(i64 1, i64 %18)
%21 = shl nsw i32 %3, 1
%22 = add nsw i32 %21, -2
store i32 %22, i32* %16
%23 = sext i32 %4 to i64
%24 = getelementptr inbounds double, double* %1, i64 %23
br label %25
25: ; preds = %65, %5
%26 = phi i32 [ %22, %5 ], [ %69, %65 ]
%27 = phi i32 [ 0, %5 ], [ %66, %65 ]
%28 = sub nsw i32 %26, %3
%29 = shl nsw i32 %28, 2
%30 = sext i32 %29 to i64
%31 = getelementptr inbounds double, double* %0, i64 %30
%32 = load double, double* %31
%33 = fptosi double %32 to i32
%34 = getelementptr inbounds double, double* %31, i64 1
%35 = load double, double* %34
%36 = fptosi double %35 to i32
%37 = icmp slt i32 %33, %3
br i1 %37, label %70, label %38
38: ; preds = %25
%39 = sext i32 %28 to i64
%40 = lshr i64 %39, 3
%41 = getelementptr inbounds i8, i8* %19, i64 %40
%42 = load i8, i8* %41
%43 = zext i8 %42 to i32
%44 = and i32 %28, 7
%45 = xor i32 %44, 7
%46 = shl nuw i32 1, %45
%47 = and i32 %46, %43
%48 = icmp eq i32 %47, 0
br i1 %48, label %55, label %70
49: ; preds = %93, %95
%50 = phi double [ %101, %95 ], [ %94, %93 ]
%51 = sext i32 %28 to i64
%52 = getelementptr inbounds double, double* %2, i64 %51
store double %50, double* %52
%53 = add nsw i32 %27, -1
%54 = icmp sgt i32 %27, 0
br i1 %54, label %65, label %102
55: ; preds = %38, %72
%56 = phi i32 [ %80, %72 ], [ %46, %38 ]
%57 = phi i8 [ %76, %72 ], [ %42, %38 ]
%58 = phi i8* [ %75, %72 ], [ %41, %38 ]
%59 = phi i32 [ %36, %72 ], [ %33, %38 ]
%60 = trunc i32 %56 to i8
%61 = or i8 %57, %60
store i8 %61, i8* %58
%62 = add nuw nsw i32 %27, 1
%63 = zext i32 %62 to i64
%64 = getelementptr inbounds i32, i32* %16, i64 %63
store i32 %59, i32* %64
br label %65
65: ; preds = %55, %49
%66 = phi i32 [ %53, %49 ], [ %62, %55 ]
%67 = zext i32 %66 to i64
%68 = getelementptr inbounds i32, i32* %16, i64 %67
%69 = load i32, i32* %68
br label %25
70: ; preds = %38, %25
%71 = icmp slt i32 %36, %3
br i1 %71, label %83, label %72
72: ; preds = %70
%73 = sext i32 %28 to i64
%74 = lshr i64 %73, 3
%75 = getelementptr inbounds i8, i8* %20, i64 %74
%76 = load i8, i8* %75
%77 = zext i8 %76 to i32
%78 = and i32 %28, 7
%79 = xor i32 %78, 7
%80 = shl nuw i32 1, %79
%81 = and i32 %80, %77
%82 = icmp eq i32 %81, 0
br i1 %82, label %55, label %83
83: ; preds = %72, %70
%84 = getelementptr inbounds double, double* %24, i64 %30
%85 = load double, double* %84
br i1 %37, label %93, label %86
86: ; preds = %83
%87 = sub nsw i32 %33, %3
%88 = sext i32 %87 to i64
%89 = getelementptr inbounds double, double* %2, i64 %88
%90 = load double, double* %89
%91 = fcmp ogt double %85, %90
%92 = select i1 %91, double %85, double %90
br label %93
93: ; preds = %86, %83
%94 = phi double [ %92, %86 ], [ %85, %83 ]
br i1 %71, label %49, label %95
95: ; preds = %93
%96 = sub nsw i32 %36, %3
%97 = sext i32 %96 to i64
%98 = getelementptr inbounds double, double* %2, i64 %97
%99 = load double, double* %98
%100 = fcmp ogt double %94, %99
%101 = select i1 %100, double %94, double %99
br label %49
102: ; preds = %49
tail call void @free(i8* %15)
tail call void @free(i8* %19)
tail call void @free(i8* %20)
ret void
}