#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 */
double fmax(double, double);
double fmin(double, double);
int find_node_split_dim(double **data, int *node_indices, int n_features,
                        int n_points) {
  double min_val, max_val, val, spread, max_spread;
  int i, j, j_max;
  j_max = 0;
  max_spread = 0;
  for (j = 0; j < n_features; j++) {
    max_val = data[node_indices[0]][j];
    min_val = max_val;
    for (i = 1; i < n_points; i++) {
      val = data[node_indices[i]][j];
      max_val = fmax(max_val, val);
      min_val = fmin(min_val, val);
    }
    spread = max_val - min_val;
    if (spread > max_spread) {
      max_spread = spread;
      j_max = j;
    }
  }
  return (j_max);
}


x86 assembly (clang, -O0):
<find_node_split_dim>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %edx,-0x14(%rbp)
mov    %ecx,-0x18(%rbp)
movl   $0x0,-0x4c(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x40(%rbp)
movl   $0x0,-0x48(%rbp)
mov    -0x48(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    123a <find_node_split_dim+0x13a>
mov    -0x8(%rbp),%rax
mov    -0x10(%rbp),%rcx
movslq (%rcx),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x48(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x28(%rbp)
movsd  -0x28(%rbp),%xmm0
movsd  %xmm0,-0x20(%rbp)
movl   $0x1,-0x44(%rbp)
mov    -0x44(%rbp),%eax
cmp    -0x18(%rbp),%eax
jge    11f8 <find_node_split_dim+0xf8>
mov    -0x8(%rbp),%rax
mov    -0x10(%rbp),%rcx
movslq -0x44(%rbp),%rdx
movslq (%rcx,%rdx,4),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x48(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x30(%rbp)
movsd  -0x28(%rbp),%xmm3
movsd  -0x30(%rbp),%xmm2
movaps %xmm2,%xmm4
movaps %xmm3,%xmm0
cmpunordsd %xmm0,%xmm0
movaps %xmm0,%xmm1
pand   %xmm4,%xmm1
maxsd  %xmm3,%xmm2
pandn  %xmm2,%xmm0
por    %xmm1,%xmm0
movsd  %xmm0,-0x28(%rbp)
movsd  -0x20(%rbp),%xmm3
movsd  -0x30(%rbp),%xmm2
movaps %xmm2,%xmm4
movaps %xmm3,%xmm0
cmpunordsd %xmm0,%xmm0
movaps %xmm0,%xmm1
pand   %xmm4,%xmm1
minsd  %xmm3,%xmm2
pandn  %xmm2,%xmm0
por    %xmm1,%xmm0
movsd  %xmm0,-0x20(%rbp)
mov    -0x44(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x44(%rbp)
jmp    1162 <find_node_split_dim+0x62>
movsd  -0x28(%rbp),%xmm0
subsd  -0x20(%rbp),%xmm0
movsd  %xmm0,-0x38(%rbp)
movsd  -0x38(%rbp),%xmm0
ucomisd -0x40(%rbp),%xmm0
jbe    1227 <find_node_split_dim+0x127>
movsd  -0x38(%rbp),%xmm0
movsd  %xmm0,-0x40(%rbp)
mov    -0x48(%rbp),%eax
mov    %eax,-0x4c(%rbp)
jmp    122c <find_node_split_dim+0x12c>
mov    -0x48(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x48(%rbp)
jmp    1128 <find_node_split_dim+0x28>
mov    -0x4c(%rbp),%eax
pop    %rbp
ret

x86 assembly (clang, -O2):
<find_node_split_dim>:
push   %r14
push   %rbx
test   %edx,%edx
jle    12be <find_node_split_dim+0x1be>
movslq (%rsi),%rax
mov    (%rdi,%rax,8),%r8
mov    %edx,%r10d
mov    %ecx,%r9d
add    $0xffffffffffffffff,%r9
mov    %r9,%r11
and    $0xfffffffffffffffe,%r11
neg    %r11
xorpd  %xmm0,%xmm0
xor    %r14d,%r14d
xor    %eax,%eax
jmp    1165 <find_node_split_dim+0x65>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
movapd %xmm2,%xmm1
subsd  %xmm1,%xmm2
ucomisd %xmm0,%xmm2
cmova  %r14d,%eax
maxsd  %xmm0,%xmm2
add    $0x1,%r14
movapd %xmm2,%xmm0
cmp    %r10,%r14
je     12c0 <find_node_split_dim+0x1c0>
movsd  (%r8,%r14,8),%xmm2
cmp    $0x2,%ecx
jl     1140 <find_node_split_dim+0x40>
mov    $0x1,%ebx
jne    1190 <find_node_split_dim+0x90>
movapd %xmm2,%xmm1
movapd %xmm2,%xmm3
jmp    125b <find_node_split_dim+0x15b>
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
mov    $0x1,%ebx
movapd %xmm2,%xmm4
movapd %xmm2,%xmm3
movapd %xmm2,%xmm1
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
movslq (%rsi,%rbx,4),%rdx
mov    (%rdi,%rdx,8),%rdx
movsd  (%rdx,%r14,8),%xmm5
cmpunordsd %xmm3,%xmm3
movapd %xmm3,%xmm6
andpd  %xmm5,%xmm6
movapd %xmm5,%xmm7
maxsd  %xmm4,%xmm7
andnpd %xmm7,%xmm3
orpd   %xmm6,%xmm3
cmpunordsd %xmm1,%xmm1
movapd %xmm1,%xmm4
andpd  %xmm5,%xmm4
minsd  %xmm2,%xmm5
andnpd %xmm5,%xmm1
orpd   %xmm4,%xmm1
movslq 0x4(%rsi,%rbx,4),%rdx
mov    (%rdi,%rdx,8),%rdx
movsd  (%rdx,%r14,8),%xmm2
movapd %xmm2,%xmm4
maxsd  %xmm3,%xmm4
cmpunordsd %xmm3,%xmm3
movapd %xmm3,%xmm5
andnpd %xmm4,%xmm5
andpd  %xmm2,%xmm3
orpd   %xmm5,%xmm3
movapd %xmm2,%xmm4
minsd  %xmm1,%xmm4
cmpunordsd %xmm1,%xmm1
movapd %xmm1,%xmm5
andnpd %xmm4,%xmm5
andpd  %xmm2,%xmm1
orpd   %xmm5,%xmm1
lea    (%r11,%rbx,1),%rdx
add    $0x2,%rdx
add    $0x2,%rbx
movapd %xmm1,%xmm2
movapd %xmm3,%xmm4
cmp    $0x1,%rdx
jne    11b0 <find_node_split_dim+0xb0>
movapd %xmm3,%xmm2
test   $0x1,%r9b
je     1144 <find_node_split_dim+0x44>
movslq (%rsi,%rbx,4),%rdx
mov    (%rdi,%rdx,8),%rdx
movsd  (%rdx,%r14,8),%xmm4
movapd %xmm3,%xmm2
cmpunordsd %xmm3,%xmm2
movapd %xmm2,%xmm5
andpd  %xmm4,%xmm5
movapd %xmm4,%xmm6
maxsd  %xmm3,%xmm6
andnpd %xmm6,%xmm2
orpd   %xmm5,%xmm2
movapd %xmm1,%xmm3
cmpunordsd %xmm1,%xmm3
movapd %xmm3,%xmm5
andpd  %xmm4,%xmm5
minsd  %xmm1,%xmm4
andnpd %xmm4,%xmm3
orpd   %xmm5,%xmm3
movapd %xmm3,%xmm1
jmp    1144 <find_node_split_dim+0x44>
xor    %eax,%eax
pop    %rbx
pop    %r14
ret

LLVM IR (-O0):
define i32 @find_node_split_dim(double** %0, i32* %1, i32 %2, i32 %3) {
%5 = alloca double**
%6 = alloca i32*
%7 = alloca i32
%8 = alloca i32
%9 = alloca double
%10 = alloca double
%11 = alloca double
%12 = alloca double
%13 = alloca double
%14 = alloca i32
%15 = alloca i32
%16 = alloca i32
store double** %0, double*** %5
store i32* %1, i32** %6
store i32 %2, i32* %7
store i32 %3, i32* %8
store i32 0, i32* %16
store double 0.000000e+00, double* %13
store i32 0, i32* %15
br label %17
17: ; preds = %72, %4
%18 = load i32, i32* %15
%19 = load i32, i32* %7
%20 = icmp slt i32 %18, %19
br i1 %20, label %21, label %75
21: ; preds = %17
%22 = load double**, double*** %5
%23 = load i32*, i32** %6
%24 = getelementptr inbounds i32, i32* %23, i64 0
%25 = load i32, i32* %24
%26 = sext i32 %25 to i64
%27 = getelementptr inbounds double*, double** %22, i64 %26
%28 = load double*, double** %27
%29 = load i32, i32* %15
%30 = sext i32 %29 to i64
%31 = getelementptr inbounds double, double* %28, i64 %30
%32 = load double, double* %31
store double %32, double* %10
%33 = load double, double* %10
store double %33, double* %9
store i32 1, i32* %14
br label %34
34: ; preds = %58, %21
%35 = load i32, i32* %14
%36 = load i32, i32* %8
%37 = icmp slt i32 %35, %36
br i1 %37, label %38, label %61
38: ; preds = %34
%39 = load double**, double*** %5
%40 = load i32*, i32** %6
%41 = load i32, i32* %14
%42 = sext i32 %41 to i64
%43 = getelementptr inbounds i32, i32* %40, i64 %42
%44 = load i32, i32* %43
%45 = sext i32 %44 to i64
%46 = getelementptr inbounds double*, double** %39, i64 %45
%47 = load double*, double** %46
%48 = load i32, i32* %15
%49 = sext i32 %48 to i64
%50 = getelementptr inbounds double, double* %47, i64 %49
%51 = load double, double* %50
store double %51, double* %11
%52 = load double, double* %10
%53 = load double, double* %11
%54 = call double @llvm.maxnum.f64(double %52, double %53)
store double %54, double* %10
%55 = load double, double* %9
%56 = load double, double* %11
%57 = call double @llvm.minnum.f64(double %55, double %56)
store double %57, double* %9
br label %58
58: ; preds = %38
%59 = load i32, i32* %14
%60 = add nsw i32 %59, 1
store i32 %60, i32* %14
br label %34
61: ; preds = %34
%62 = load double, double* %10
%63 = load double, double* %9
%64 = fsub double %62, %63
store double %64, double* %12
%65 = load double, double* %12
%66 = load double, double* %13
%67 = fcmp ogt double %65, %66
br i1 %67, label %68, label %71
68: ; preds = %61
%69 = load double, double* %12
store double %69, double* %13
%70 = load i32, i32* %15
store i32 %70, i32* %16
br label %71
71: ; preds = %68, %61
br label %72
72: ; preds = %71
%73 = load i32, i32* %15
%74 = add nsw i32 %73, 1
store i32 %74, i32* %15
br label %17
75: ; preds = %17
%76 = load i32, i32* %16
ret i32 %76
}

LLVM IR (-O2):
define i32 @find_node_split_dim(double** nocapture %0, i32* nocapture %1, i32 %2, i32 %3) {
%5 = icmp sgt i32 %2, 0
br i1 %5, label %6, label %79
6: ; preds = %4
%7 = load i32, i32* %1
%8 = sext i32 %7 to i64
%9 = getelementptr inbounds double*, double** %0, i64 %8
%10 = load double*, double** %9
%11 = icmp sgt i32 %3, 1
%12 = zext i32 %2 to i64
%13 = zext i32 %3 to i64
%14 = add nsw i64 %13, -1
%15 = and i64 %14, 1
%16 = icmp eq i32 %3, 2
%17 = and i64 %14, -2
%18 = icmp eq i64 %15, 0
br label %19
19: ; preds = %6, %69
%20 = phi i64 [ 0, %6 ], [ %77, %69 ]
%21 = phi i32 [ 0, %6 ], [ %76, %69 ]
%22 = phi double [ 0.000000e+00, %6 ], [ %74, %69 ]
%23 = getelementptr inbounds double, double* %10, i64 %20
%24 = load double, double* %23
br i1 %11, label %25, label %69
25: ; preds = %19
br i1 %16, label %53, label %26
26: ; preds = %25, %26
%27 = phi i64 [ %50, %26 ], [ 1, %25 ]
%28 = phi double [ %49, %26 ], [ %24, %25 ]
%29 = phi double [ %48, %26 ], [ %24, %25 ]
%30 = phi i64 [ %51, %26 ], [ 0, %25 ]
%31 = getelementptr inbounds i32, i32* %1, i64 %27
%32 = load i32, i32* %31
%33 = sext i32 %32 to i64
%34 = getelementptr inbounds double*, double** %0, i64 %33
%35 = load double*, double** %34
%36 = getelementptr inbounds double, double* %35, i64 %20
%37 = load double, double* %36
%38 = tail call double @llvm.maxnum.f64(double %29, double %37)
%39 = tail call double @llvm.minnum.f64(double %28, double %37)
%40 = add nuw nsw i64 %27, 1
%41 = getelementptr inbounds i32, i32* %1, i64 %40
%42 = load i32, i32* %41
%43 = sext i32 %42 to i64
%44 = getelementptr inbounds double*, double** %0, i64 %43
%45 = load double*, double** %44
%46 = getelementptr inbounds double, double* %45, i64 %20
%47 = load double, double* %46
%48 = tail call double @llvm.maxnum.f64(double %38, double %47)
%49 = tail call double @llvm.minnum.f64(double %39, double %47)
%50 = add nuw nsw i64 %27, 2
%51 = add i64 %30, 2
%52 = icmp eq i64 %51, %17
br i1 %52, label %53, label %26
53: ; preds = %26, %25
%54 = phi double [ undef, %25 ], [ %48, %26 ]
%55 = phi double [ undef, %25 ], [ %49, %26 ]
%56 = phi i64 [ 1, %25 ], [ %50, %26 ]
%57 = phi double [ %24, %25 ], [ %49, %26 ]
%58 = phi double [ %24, %25 ], [ %48, %26 ]
br i1 %18, label %69, label %59
59: ; preds = %53
%60 = getelementptr inbounds i32, i32* %1, i64 %56
%61 = load i32, i32* %60
%62 = sext i32 %61 to i64
%63 = getelementptr inbounds double*, double** %0, i64 %62
%64 = load double*, double** %63
%65 = getelementptr inbounds double, double* %64, i64 %20
%66 = load double, double* %65
%67 = tail call double @llvm.maxnum.f64(double %58, double %66)
%68 = tail call double @llvm.minnum.f64(double %57, double %66)
br label %69
69: ; preds = %59, %53, %19
%70 = phi double [ %24, %19 ], [ %54, %53 ], [ %67, %59 ]
%71 = phi double [ %24, %19 ], [ %55, %53 ], [ %68, %59 ]
%72 = fsub double %70, %71
%73 = fcmp ogt double %72, %22
%74 = select i1 %73, double %72, double %22
%75 = trunc i64 %20 to i32
%76 = select i1 %73, i32 %75, i32 %21
%77 = add nuw nsw i64 %20, 1
%78 = icmp eq i64 %77, %12
br i1 %78, label %79, label %19
79: ; preds = %69, %4
%80 = phi i32 [ 0, %4 ], [ %76, %69 ]
ret i32 %80
}