#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 log10(int);
double get_block_likelihood(int branch_genome_size, int number_of_branch_snps,
                            int block_genome_size_without_gaps,
                            int number_of_block_snps) {
  double part1, part2, part3, part4;
  if (block_genome_size_without_gaps == 0) {
    return 0.0;
  }
  if (number_of_block_snps == 0) {
    return 0.0;
  }
  part1 =
      log10(number_of_block_snps * 1.0 / block_genome_size_without_gaps * 1.0) *
      number_of_block_snps;
  if ((block_genome_size_without_gaps - number_of_block_snps) == 0) {
    part2 = 0.0;
  } else {
    part2 =
        log10(((block_genome_size_without_gaps - number_of_block_snps) * 1.0) /
              block_genome_size_without_gaps * 1.0) *
        (block_genome_size_without_gaps - number_of_block_snps);
  }
  if ((number_of_branch_snps - number_of_block_snps) == 0) {
    part3 = 0.0;
  } else {
    part3 =
        log10(((number_of_branch_snps - number_of_block_snps) * 1.0) /
              ((branch_genome_size - block_genome_size_without_gaps) * 1.0)) *
        (number_of_branch_snps - number_of_block_snps);
  }
  if (((branch_genome_size - block_genome_size_without_gaps) -
       (number_of_branch_snps - number_of_block_snps)) == 0) {
    part4 = 0.0;
  } else {
    part4 =
        log10(((branch_genome_size - block_genome_size_without_gaps) -
               (number_of_branch_snps - number_of_block_snps) * 1.0) /
              ((branch_genome_size - block_genome_size_without_gaps) * 1.0)) *
        ((branch_genome_size - block_genome_size_without_gaps) -
         (number_of_branch_snps - number_of_block_snps));
  }
  return (part1 + part2 + part3 + part4) * -1;
}


x86 assembly (clang, -O0):
<get_block_likelihood>:
push   %rbp
mov    %rsp,%rbp
sub    $0x40,%rsp
mov    %edi,-0xc(%rbp)
mov    %esi,-0x10(%rbp)
mov    %edx,-0x14(%rbp)
mov    %ecx,-0x18(%rbp)
cmpl   $0x0,-0x14(%rbp)
jne    113b <get_block_likelihood+0x2b>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x8(%rbp)
jmp    130c <get_block_likelihood+0x1fc>
cmpl   $0x0,-0x18(%rbp)
jne    1152 <get_block_likelihood+0x42>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x8(%rbp)
jmp    130c <get_block_likelihood+0x1fc>
cvtsi2sdl -0x18(%rbp),%xmm0
movsd  0xea1(%rip),%xmm1
mulsd  %xmm1,%xmm0
cvtsi2sdl -0x14(%rbp),%xmm1
divsd  %xmm1,%xmm0
movsd  0xe8c(%rip),%xmm1
mulsd  %xmm1,%xmm0
cvttsd2si %xmm0,%edi
call   1030 <log10@plt>
imul   -0x18(%rbp),%eax
cvtsi2sd %eax,%xmm0
movsd  %xmm0,-0x20(%rbp)
mov    -0x14(%rbp),%eax
sub    -0x18(%rbp),%eax
cmp    $0x0,%eax
jne    11aa <get_block_likelihood+0x9a>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x28(%rbp)
jmp    11f0 <get_block_likelihood+0xe0>
mov    -0x14(%rbp),%eax
sub    -0x18(%rbp),%eax
cvtsi2sd %eax,%xmm0
movsd  0xe44(%rip),%xmm1
mulsd  %xmm1,%xmm0
cvtsi2sdl -0x14(%rbp),%xmm1
divsd  %xmm1,%xmm0
movsd  0xe2f(%rip),%xmm1
mulsd  %xmm1,%xmm0
cvttsd2si %xmm0,%edi
call   1030 <log10@plt>
mov    -0x14(%rbp),%ecx
sub    -0x18(%rbp),%ecx
imul   %ecx,%eax
cvtsi2sd %eax,%xmm0
movsd  %xmm0,-0x28(%rbp)
mov    -0x10(%rbp),%eax
sub    -0x18(%rbp),%eax
cmp    $0x0,%eax
jne    120c <get_block_likelihood+0xfc>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x30(%rbp)
jmp    1257 <get_block_likelihood+0x147>
mov    -0x10(%rbp),%eax
sub    -0x18(%rbp),%eax
cvtsi2sd %eax,%xmm0
movsd  0xde2(%rip),%xmm1
mulsd  %xmm1,%xmm0
mov    -0xc(%rbp),%eax
sub    -0x14(%rbp),%eax
cvtsi2sd %eax,%xmm1
movsd  0xdcc(%rip),%xmm2
mulsd  %xmm2,%xmm1
divsd  %xmm1,%xmm0
cvttsd2si %xmm0,%edi
call   1030 <log10@plt>
mov    -0x10(%rbp),%ecx
sub    -0x18(%rbp),%ecx
imul   %ecx,%eax
cvtsi2sd %eax,%xmm0
movsd  %xmm0,-0x30(%rbp)
mov    -0xc(%rbp),%eax
sub    -0x14(%rbp),%eax
mov    -0x10(%rbp),%ecx
sub    -0x18(%rbp),%ecx
sub    %ecx,%eax
cmp    $0x0,%eax
jne    127b <get_block_likelihood+0x16b>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x38(%rbp)
jmp    12e7 <get_block_likelihood+0x1d7>
mov    -0xc(%rbp),%eax
sub    -0x14(%rbp),%eax
cvtsi2sd %eax,%xmm1
mov    -0x10(%rbp),%eax
sub    -0x18(%rbp),%eax
cvtsi2sd %eax,%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
addsd  %xmm1,%xmm0
mov    -0xc(%rbp),%eax
sub    -0x14(%rbp),%eax
cvtsi2sd %eax,%xmm1
movsd  0xd44(%rip),%xmm2
mulsd  %xmm2,%xmm1
divsd  %xmm1,%xmm0
cvttsd2si %xmm0,%edi
call   1030 <log10@plt>
mov    -0xc(%rbp),%ecx
sub    -0x14(%rbp),%ecx
mov    -0x10(%rbp),%edx
sub    -0x18(%rbp),%edx
sub    %edx,%ecx
imul   %ecx,%eax
cvtsi2sd %eax,%xmm0
movsd  %xmm0,-0x38(%rbp)
movsd  -0x20(%rbp),%xmm0
addsd  -0x28(%rbp),%xmm0
addsd  -0x30(%rbp),%xmm0
addsd  -0x38(%rbp),%xmm0
movsd  0xd05(%rip),%xmm1
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0x8(%rbp)
movsd  -0x8(%rbp),%xmm0
add    $0x40,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<get_block_likelihood>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x28,%rsp
xorpd  %xmm0,%xmm0
test   %edx,%edx
je     123f <get_block_likelihood+0x12f>
mov    %ecx,%ebp
test   %ecx,%ecx
je     123f <get_block_likelihood+0x12f>
mov    %edx,%r13d
mov    %esi,%r12d
mov    %edi,%r15d
xorps  %xmm0,%xmm0
cvtsi2sd %ebp,%xmm0
cvtsi2sd %edx,%xmm1
movsd  %xmm1,0x10(%rsp)
divsd  %xmm1,%xmm0
cvttsd2si %xmm0,%edi
call   1030 <log10@plt>
mov    %eax,%r14d
mov    %r13d,%ebx
sub    %ebp,%ebx
je     1187 <get_block_likelihood+0x77>
xorps  %xmm0,%xmm0
cvtsi2sd %ebx,%xmm0
divsd  0x10(%rsp),%xmm0
cvttsd2si %xmm0,%edi
call   1030 <log10@plt>
imul   %ebx,%eax
xorps  %xmm0,%xmm0
cvtsi2sd %eax,%xmm0
jmp    118b <get_block_likelihood+0x7b>
xorpd  %xmm0,%xmm0
imul   %ebp,%r14d
sub    %ebp,%r12d
xorps  %xmm1,%xmm1
cvtsi2sd %r12d,%xmm1
movapd %xmm0,0x10(%rsp)
movsd  %xmm1,0x8(%rsp)
je     11ca <get_block_likelihood+0xba>
sub    %r13d,%r15d
xorps  %xmm0,%xmm0
cvtsi2sd %r15d,%xmm0
divsd  %xmm0,%xmm1
cvttsd2si %xmm1,%edi
call   1030 <log10@plt>
imul   %r12d,%eax
cvtsi2sd %eax,%xmm2
jmp    11d1 <get_block_likelihood+0xc1>
sub    %r13d,%r15d
xorpd  %xmm2,%xmm2
cvtsi2sd %r14d,%xmm3
mov    %r15d,%ebx
sub    %r12d,%ebx
je     1221 <get_block_likelihood+0x111>
xorps  %xmm0,%xmm0
cvtsi2sd %r15d,%xmm0
movapd %xmm0,%xmm1
subsd  0x8(%rsp),%xmm1
divsd  %xmm0,%xmm1
cvttsd2si %xmm1,%edi
movsd  %xmm2,0x8(%rsp)
movsd  %xmm3,0x20(%rsp)
call   1030 <log10@plt>
movsd  0x20(%rsp),%xmm3
movsd  0x8(%rsp),%xmm2
imul   %ebx,%eax
xorps  %xmm1,%xmm1
cvtsi2sd %eax,%xmm1
jmp    1225 <get_block_likelihood+0x115>
xorpd  %xmm1,%xmm1
movapd 0x10(%rsp),%xmm0
addsd  %xmm3,%xmm0
addsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
xorpd  0xdc1(%rip),%xmm0
add    $0x28,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

LLVM IR (-O0):
define double @get_block_likelihood(i32 %0, i32 %1, i32 %2, i32 %3) {
%5 = alloca double
%6 = alloca i32
%7 = alloca i32
%8 = alloca i32
%9 = alloca i32
%10 = alloca double
%11 = alloca double
%12 = alloca double
%13 = alloca double
store i32 %0, i32* %6
store i32 %1, i32* %7
store i32 %2, i32* %8
store i32 %3, i32* %9
%14 = load i32, i32* %8
%15 = icmp eq i32 %14, 0
br i1 %15, label %16, label %17
16: ; preds = %4
store double 0.000000e+00, double* %5
br label %128
17: ; preds = %4
%18 = load i32, i32* %9
%19 = icmp eq i32 %18, 0
br i1 %19, label %20, label %21
20: ; preds = %17
store double 0.000000e+00, double* %5
br label %128
21: ; preds = %17
%22 = load i32, i32* %9
%23 = sitofp i32 %22 to double
%24 = fmul double %23, 1.000000e+00
%25 = load i32, i32* %8
%26 = sitofp i32 %25 to double
%27 = fdiv double %24, %26
%28 = fmul double %27, 1.000000e+00
%29 = fptosi double %28 to i32
%30 = call i32 @log10(i32 %29)
%31 = load i32, i32* %9
%32 = mul nsw i32 %30, %31
%33 = sitofp i32 %32 to double
store double %33, double* %10
%34 = load i32, i32* %8
%35 = load i32, i32* %9
%36 = sub nsw i32 %34, %35
%37 = icmp eq i32 %36, 0
br i1 %37, label %38, label %39
38: ; preds = %21
store double 0.000000e+00, double* %11
br label %56
39: ; preds = %21
%40 = load i32, i32* %8
%41 = load i32, i32* %9
%42 = sub nsw i32 %40, %41
%43 = sitofp i32 %42 to double
%44 = fmul double %43, 1.000000e+00
%45 = load i32, i32* %8
%46 = sitofp i32 %45 to double
%47 = fdiv double %44, %46
%48 = fmul double %47, 1.000000e+00
%49 = fptosi double %48 to i32
%50 = call i32 @log10(i32 %49)
%51 = load i32, i32* %8
%52 = load i32, i32* %9
%53 = sub nsw i32 %51, %52
%54 = mul nsw i32 %50, %53
%55 = sitofp i32 %54 to double
store double %55, double* %11
br label %56
56: ; preds = %39, %38
%57 = load i32, i32* %7
%58 = load i32, i32* %9
%59 = sub nsw i32 %57, %58
%60 = icmp eq i32 %59, 0
br i1 %60, label %61, label %62
61: ; preds = %56
store double 0.000000e+00, double* %12
br label %81
62: ; preds = %56
%63 = load i32, i32* %7
%64 = load i32, i32* %9
%65 = sub nsw i32 %63, %64
%66 = sitofp i32 %65 to double
%67 = fmul double %66, 1.000000e+00
%68 = load i32, i32* %6
%69 = load i32, i32* %8
%70 = sub nsw i32 %68, %69
%71 = sitofp i32 %70 to double
%72 = fmul double %71, 1.000000e+00
%73 = fdiv double %67, %72
%74 = fptosi double %73 to i32
%75 = call i32 @log10(i32 %74)
%76 = load i32, i32* %7
%77 = load i32, i32* %9
%78 = sub nsw i32 %76, %77
%79 = mul nsw i32 %75, %78
%80 = sitofp i32 %79 to double
store double %80, double* %12
br label %81
81: ; preds = %62, %61
%82 = load i32, i32* %6
%83 = load i32, i32* %8
%84 = sub nsw i32 %82, %83
%85 = load i32, i32* %7
%86 = load i32, i32* %9
%87 = sub nsw i32 %85, %86
%88 = sub nsw i32 %84, %87
%89 = icmp eq i32 %88, 0
br i1 %89, label %90, label %91
90: ; preds = %81
store double 0.000000e+00, double* %13
br label %119
91: ; preds = %81
%92 = load i32, i32* %6
%93 = load i32, i32* %8
%94 = sub nsw i32 %92, %93
%95 = sitofp i32 %94 to double
%96 = load i32, i32* %7
%97 = load i32, i32* %9
%98 = sub nsw i32 %96, %97
%99 = sitofp i32 %98 to double
%100 = fneg double %99
%101 = call double @llvm.fmuladd.f64(double %100, double 1.000000e+00, double %95)
%102 = load i32, i32* %6
%103 = load i32, i32* %8
%104 = sub nsw i32 %102, %103
%105 = sitofp i32 %104 to double
%106 = fmul double %105, 1.000000e+00
%107 = fdiv double %101, %106
%108 = fptosi double %107 to i32
%109 = call i32 @log10(i32 %108)
%110 = load i32, i32* %6
%111 = load i32, i32* %8
%112 = sub nsw i32 %110, %111
%113 = load i32, i32* %7
%114 = load i32, i32* %9
%115 = sub nsw i32 %113, %114
%116 = sub nsw i32 %112, %115
%117 = mul nsw i32 %109, %116
%118 = sitofp i32 %117 to double
store double %118, double* %13
br label %119
119: ; preds = %91, %90
%120 = load double, double* %10
%121 = load double, double* %11
%122 = fadd double %120, %121
%123 = load double, double* %12
%124 = fadd double %122, %123
%125 = load double, double* %13
%126 = fadd double %124, %125
%127 = fmul double %126, -1.000000e+00
store double %127, double* %5
br label %128
128: ; preds = %119, %20, %16
%129 = load double, double* %5
ret double %129
}

LLVM IR (-O2):
define double @get_block_likelihood(i32 %0, i32 %1, i32 %2, i32 %3) {
%5 = icmp eq i32 %2, 0
%6 = icmp eq i32 %3, 0
%7 = or i1 %5, %6
br i1 %7, label %60, label %8
8: ; preds = %4
%9 = sitofp i32 %3 to double
%10 = sitofp i32 %2 to double
%11 = fdiv double %9, %10
%12 = fptosi double %11 to i32
%13 = tail call i32 @log10(i32 %12)
%14 = mul nsw i32 %13, %3
%15 = sitofp i32 %14 to double
%16 = sub nsw i32 %2, %3
%17 = icmp eq i32 %16, 0
br i1 %17, label %25, label %18
18: ; preds = %8
%19 = sitofp i32 %16 to double
%20 = fdiv double %19, %10
%21 = fptosi double %20 to i32
%22 = tail call i32 @log10(i32 %21)
%23 = mul nsw i32 %22, %16
%24 = sitofp i32 %23 to double
br label %25
25: ; preds = %8, %18
%26 = phi double [ %24, %18 ], [ 0.000000e+00, %8 ]
%27 = sub nsw i32 %1, %3
%28 = icmp eq i32 %27, 0
br i1 %28, label %29, label %31
29: ; preds = %25
%30 = sub nsw i32 %0, %2
br label %40
31: ; preds = %25
%32 = sitofp i32 %27 to double
%33 = sub nsw i32 %0, %2
%34 = sitofp i32 %33 to double
%35 = fdiv double %32, %34
%36 = fptosi double %35 to i32
%37 = tail call i32 @log10(i32 %36)
%38 = mul nsw i32 %37, %27
%39 = sitofp i32 %38 to double
br label %40
40: ; preds = %29, %31
%41 = phi i32 [ %30, %29 ], [ %33, %31 ]
%42 = phi double [ 0.000000e+00, %29 ], [ %39, %31 ]
%43 = sub nsw i32 %41, %27
%44 = icmp eq i32 %43, 0
br i1 %44, label %54, label %45
45: ; preds = %40
%46 = sitofp i32 %41 to double
%47 = sitofp i32 %27 to double
%48 = fsub double %46, %47
%49 = fdiv double %48, %46
%50 = fptosi double %49 to i32
%51 = tail call i32 @log10(i32 %50)
%52 = mul nsw i32 %51, %43
%53 = sitofp i32 %52 to double
br label %54
54: ; preds = %40, %45
%55 = phi double [ %53, %45 ], [ 0.000000e+00, %40 ]
%56 = fadd double %26, %15
%57 = fadd double %56, %42
%58 = fadd double %57, %55
%59 = fneg double %58
br label %60
60: ; preds = %4, %54
%61 = phi double [ %59, %54 ], [ 0.000000e+00, %4 ]
ret double %61
}