#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 pow(double, double);
double getCustomSplitStatisticMultivariateRegression(
    unsigned int n, char *membership, double *time, double *event,
    unsigned int eventTypeSize, unsigned int eventTimeSize, double *eventTime,
    double *response, double mean, double variance, unsigned int maxLevel,
    double **feature, unsigned int featureCount) {
  double sumLeft, sumRght;
  double sumLeftSqr, sumRghtSqr;
  double delta;
  unsigned int leftSize, rghtSize;
  unsigned int i;
  sumLeft = sumRght = 0.0;
  leftSize = rghtSize = 0;
  delta = 0.0;
  for (i = 1; i <= n; i++) {
    if (membership[i] == 0x01) {
      sumLeft += response[i] - mean;
      leftSize++;
    } else {
      sumRght += response[i] - mean;
      rghtSize++;
    }
  }
  sumLeftSqr = pow(sumLeft, 2.0) / ((double)leftSize * variance);
  sumRghtSqr = pow(sumRght, 2.0) / ((double)rghtSize * variance);
  delta = sumLeftSqr + sumRghtSqr;
  return delta;
}


x86 assembly (clang, -O0):
<getCustomSplitStatisticMultivariateRegression>:
push   %rbp
mov    %rsp,%rbp
sub    $0x80,%rsp
mov    0x30(%rbp),%eax
mov    0x28(%rbp),%rax
mov    0x20(%rbp),%eax
mov    0x18(%rbp),%rax
mov    0x10(%rbp),%rax
mov    %edi,-0x4(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %rcx,-0x20(%rbp)
mov    %r8d,-0x24(%rbp)
mov    %r9d,-0x28(%rbp)
movsd  %xmm0,-0x30(%rbp)
movsd  %xmm1,-0x38(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x48(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x40(%rbp)
movl   $0x0,-0x68(%rbp)
movl   $0x0,-0x64(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x60(%rbp)
movl   $0x1,-0x6c(%rbp)
mov    -0x6c(%rbp),%eax
cmp    -0x4(%rbp),%eax
ja     11fb <getCustomSplitStatisticMultivariateRegression+0xeb>
mov    -0x10(%rbp),%rax
mov    -0x6c(%rbp),%ecx
movsbl (%rax,%rcx,1),%eax
cmp    $0x1,%eax
jne    11c4 <getCustomSplitStatisticMultivariateRegression+0xb4>
mov    0x18(%rbp),%rax
mov    -0x6c(%rbp),%ecx
movsd  (%rax,%rcx,8),%xmm0
subsd  -0x30(%rbp),%xmm0
addsd  -0x40(%rbp),%xmm0
movsd  %xmm0,-0x40(%rbp)
mov    -0x64(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x64(%rbp)
jmp    11e8 <getCustomSplitStatisticMultivariateRegression+0xd8>
mov    0x18(%rbp),%rax
mov    -0x6c(%rbp),%ecx
movsd  (%rax,%rcx,8),%xmm0
subsd  -0x30(%rbp),%xmm0
addsd  -0x48(%rbp),%xmm0
movsd  %xmm0,-0x48(%rbp)
mov    -0x68(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x68(%rbp)
jmp    11ed <getCustomSplitStatisticMultivariateRegression+0xdd>
mov    -0x6c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x6c(%rbp)
jmp    117b <getCustomSplitStatisticMultivariateRegression+0x6b>
movsd  -0x40(%rbp),%xmm0
movsd  0xdf8(%rip),%xmm1
movsd  %xmm1,-0x78(%rbp)
call   1030 <pow@plt>
movsd  -0x78(%rbp),%xmm1
mov    -0x64(%rbp),%eax
cvtsi2sd %rax,%xmm2
movsd  -0x38(%rbp),%xmm3
mulsd  %xmm3,%xmm2
divsd  %xmm2,%xmm0
movsd  %xmm0,-0x50(%rbp)
movsd  -0x48(%rbp),%xmm0
call   1030 <pow@plt>
mov    -0x68(%rbp),%eax
cvtsi2sd %rax,%xmm1
mulsd  -0x38(%rbp),%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x58(%rbp)
movsd  -0x50(%rbp),%xmm0
addsd  -0x58(%rbp),%xmm0
movsd  %xmm0,-0x60(%rbp)
movsd  -0x60(%rbp),%xmm0
add    $0x80,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<getCustomSplitStatisticMultivariateRegression>:
test   %edi,%edi
je     1150 <getCustomSplitStatisticMultivariateRegression+0x50>
movapd %xmm0,%xmm2
mov    0x10(%rsp),%r9
mov    %edi,%r8d
cmp    $0x1,%edi
jne    1161 <getCustomSplitStatisticMultivariateRegression+0x61>
xorpd  %xmm3,%xmm3
xor    %eax,%eax
mov    $0x1,%edi
xor    %ecx,%ecx
xorpd  %xmm0,%xmm0
test   $0x1,%r8b
je     11e9 <getCustomSplitStatisticMultivariateRegression+0xe9>
cmpb   $0x1,(%rsi,%rdi,1)
movsd  (%r9,%rdi,8),%xmm4
subsd  %xmm2,%xmm4
jne    11e2 <getCustomSplitStatisticMultivariateRegression+0xe2>
addsd  %xmm4,%xmm3
add    $0x1,%ecx
jmp    11e9 <getCustomSplitStatisticMultivariateRegression+0xe9>
xor    %ecx,%ecx
xorpd  %xmm0,%xmm0
xorpd  %xmm3,%xmm3
xor    %eax,%eax
jmp    11e9 <getCustomSplitStatisticMultivariateRegression+0xe9>
mov    %r8d,%r10d
and    $0xfffffffe,%r10d
neg    %r10
xorpd  %xmm3,%xmm3
xor    %eax,%eax
mov    $0x1,%edi
xor    %ecx,%ecx
xorpd  %xmm0,%xmm0
jmp    1199 <getCustomSplitStatisticMultivariateRegression+0x99>
xchg   %ax,%ax
addsd  %xmm4,%xmm3
add    $0x1,%ecx
lea    (%r10,%rdi,1),%rdx
add    $0x2,%rdx
add    $0x2,%rdi
cmp    $0x1,%rdx
je     1126 <getCustomSplitStatisticMultivariateRegression+0x26>
cmpb   $0x1,(%rsi,%rdi,1)
movsd  (%r9,%rdi,8),%xmm4
subsd  %xmm2,%xmm4
jne    11c0 <getCustomSplitStatisticMultivariateRegression+0xc0>
addsd  %xmm4,%xmm3
add    $0x1,%ecx
jmp    11c7 <getCustomSplitStatisticMultivariateRegression+0xc7>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
addsd  %xmm4,%xmm0
add    $0x1,%eax
cmpb   $0x1,0x1(%rsi,%rdi,1)
movsd  0x8(%r9,%rdi,8),%xmm4
subsd  %xmm2,%xmm4
je     1180 <getCustomSplitStatisticMultivariateRegression+0x80>
addsd  %xmm4,%xmm0
add    $0x1,%eax
jmp    1187 <getCustomSplitStatisticMultivariateRegression+0x87>
addsd  %xmm4,%xmm0
add    $0x1,%eax
mulsd  %xmm3,%xmm3
mov    %ecx,%ecx
xorps  %xmm2,%xmm2
cvtsi2sd %rcx,%xmm2
mulsd  %xmm1,%xmm2
divsd  %xmm2,%xmm3
mov    %eax,%eax
xorps  %xmm2,%xmm2
cvtsi2sd %rax,%xmm2
mulsd  %xmm0,%xmm0
mulsd  %xmm1,%xmm2
divsd  %xmm2,%xmm0
addsd  %xmm3,%xmm0
ret

LLVM IR (-O0):
define double @getCustomSplitStatisticMultivariateRegression(i32 %0, i8* %1, double* %2, double* %3, i32 %4, i32 %5, double* %6, double* %7, double %8, double %9, i32 %10, double** %11, i32 %12) {
%14 = alloca i32
%15 = alloca i8*
%16 = alloca double*
%17 = alloca double*
%18 = alloca i32
%19 = alloca i32
%20 = alloca double*
%21 = alloca double*
%22 = alloca double
%23 = alloca double
%24 = alloca i32
%25 = alloca double**
%26 = alloca i32
%27 = alloca double
%28 = alloca double
%29 = alloca double
%30 = alloca double
%31 = alloca double
%32 = alloca i32
%33 = alloca i32
%34 = alloca i32
store i32 %0, i32* %14
store i8* %1, i8** %15
store double* %2, double** %16
store double* %3, double** %17
store i32 %4, i32* %18
store i32 %5, i32* %19
store double* %6, double** %20
store double* %7, double** %21
store double %8, double* %22
store double %9, double* %23
store i32 %10, i32* %24
store double** %11, double*** %25
store i32 %12, i32* %26
store double 0.000000e+00, double* %28
store double 0.000000e+00, double* %27
store i32 0, i32* %33
store i32 0, i32* %32
store double 0.000000e+00, double* %31
store i32 1, i32* %34
br label %35
35: ; preds = %72, %13
%36 = load i32, i32* %34
%37 = load i32, i32* %14
%38 = icmp ule i32 %36, %37
br i1 %38, label %39, label %75
39: ; preds = %35
%40 = load i8*, i8** %15
%41 = load i32, i32* %34
%42 = zext i32 %41 to i64
%43 = getelementptr inbounds i8, i8* %40, i64 %42
%44 = load i8, i8* %43
%45 = sext i8 %44 to i32
%46 = icmp eq i32 %45, 1
br i1 %46, label %47, label %59
47: ; preds = %39
%48 = load double*, double** %21
%49 = load i32, i32* %34
%50 = zext i32 %49 to i64
%51 = getelementptr inbounds double, double* %48, i64 %50
%52 = load double, double* %51
%53 = load double, double* %22
%54 = fsub double %52, %53
%55 = load double, double* %27
%56 = fadd double %55, %54
store double %56, double* %27
%57 = load i32, i32* %32
%58 = add i32 %57, 1
store i32 %58, i32* %32
br label %71
59: ; preds = %39
%60 = load double*, double** %21
%61 = load i32, i32* %34
%62 = zext i32 %61 to i64
%63 = getelementptr inbounds double, double* %60, i64 %62
%64 = load double, double* %63
%65 = load double, double* %22
%66 = fsub double %64, %65
%67 = load double, double* %28
%68 = fadd double %67, %66
store double %68, double* %28
%69 = load i32, i32* %33
%70 = add i32 %69, 1
store i32 %70, i32* %33
br label %71
71: ; preds = %59, %47
br label %72
72: ; preds = %71
%73 = load i32, i32* %34
%74 = add i32 %73, 1
store i32 %74, i32* %34
br label %35
75: ; preds = %35
%76 = load double, double* %27
%77 = call double @pow(double %76, double 2.000000e+00)
%78 = load i32, i32* %32
%79 = uitofp i32 %78 to double
%80 = load double, double* %23
%81 = fmul double %79, %80
%82 = fdiv double %77, %81
store double %82, double* %29
%83 = load double, double* %28
%84 = call double @pow(double %83, double 2.000000e+00)
%85 = load i32, i32* %33
%86 = uitofp i32 %85 to double
%87 = load double, double* %23
%88 = fmul double %86, %87
%89 = fdiv double %84, %88
store double %89, double* %30
%90 = load double, double* %29
%91 = load double, double* %30
%92 = fadd double %90, %91
store double %92, double* %31
%93 = load double, double* %31
ret double %93
}

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