// File: /workspace/github_cpp_makefile_star2_O1/CK7[P]abawaca/src/ClusterSeparatorBySensitivitySpecificity.cpp
// Function: ClusterSeparatorBySensitivitySpecificity::separate_dimension(unsigned long)
void ClusterSeparatorBySensitivitySpecificity::separate_dimension(size_t dimension_index)
{
	semaphore.wait();
	const Dimension& dimension = clustering_db.get_dimension(dimension_index);
	const set<double>& values = dimension.get_values();
	int values_checked = 0;
	// Initialize - curr_cluster1 will have no points while curr_cluster2 will have all dps
	Cluster 						curr_cluster1(scaf_db);
	Cluster 						curr_cluster2(scaf_db);
	ClusterSeparatorBySensitivitySpecificity::ClusterStats	curr_cluster_stats1(curr_cluster1, scaf_db);
	ClusterSeparatorBySensitivitySpecificity::ClusterStats	curr_cluster_stats2(curr_cluster2, scaf_db);
	vector<Dimension::DataPoint>				dps;
	set<size_t> 						useful_scafs;

	for(vector<Dimension::DataPoint>::const_iterator vit = dimension.begin(); vit!=dimension.end(); vit++) {
		curr_cluster_stats2.add_dp(vit->dp);
		dps.push_back(*vit);
		if(scaf_db.ndps(scaf_db.dp2scaf(vit->dp)) > 1) {
			useful_scafs.insert(scaf_db.dp2scaf(vit->dp));
			curr_cluster_stats2.assign_scaf(scaf_db.dp2scaf(vit->dp));
		}
	}

	sort(dps.begin(), dps.end(), comp_by_value);
	ClusterSeparatorBySensitivitySpecificity::ClusteringResult	best_dimension_separation;

	for(auto vit=dps.begin(); vit!=dps.end();) {
		double value = vit->value;
		while(vit!=dps.end() && vit->value==value) {
			size_t scaf = scaf_db.dp2scaf(vit->dp);
			size_t scaf_total_ndps = scaf_db.ndps(scaf);

			curr_cluster_stats1.add_dp(vit->dp);
			curr_cluster_stats2.remove_dp(vit->dp);
			if(curr_cluster_stats2.get_cluster().assigned(scaf) && (scaf_total_ndps >= 2) && (2*curr_cluster_stats1.get_cluster().ndps(scaf) > scaf_total_ndps)) {
				curr_cluster_stats1.assign_scaf(scaf);
				curr_cluster_stats2.unassign_scaf(scaf);
			}
			vit++;
		}

		// Now do what we actually need to do in this specific class
		this->separate_specific(curr_cluster_stats1, curr_cluster_stats2, dimension_index, value, best_dimension_separation);
	}

	// Only a single access to the variables common to all threads
	unique_lock<mutex> lck(update_mtx);

	if(best_dimension_separation.is_better(*best_separation))
		best_separation->copy(best_dimension_separation);

	semaphore.notify();
}

// Assembly:
pushq %rbp
movq %rsp, %rbp
pushq %r15
pushq %r14
pushq %r13
pushq %r12
pushq %rbx
subq $0x248, %rsp            # imm = 0x248
movq %rsi, %r14
movq %rdi, %rbx
movq 0xe0(%rdi), %rdi
callq 0x13e50
movq 0x18(%rbx), %rdi
movq %r14, -0x50(%rbp)
movq %r14, %rsi
callq 0xff94
movq %rax, %r15
movq 0x8(%rbx), %rcx
leaq -0x258(%rbp), %rdx
movq %rcx, -0x10(%rdx)
xorl %eax, %eax
movl %eax, (%rdx)
movq %rax, 0x8(%rdx)
movq %rdx, 0x10(%rdx)
movq %rdx, 0x18(%rdx)
movq %rax, 0x20(%rdx)
leaq -0x200(%rbp), %rdx
movq %rdx, -0x30(%rdx)
movl $0x1, %esi
movq %rsi, -0x28(%rdx)
xorpd %xmm0, %xmm0
movupd %xmm0, -0x20(%rdx)
movl $0x3f800000, %edi       # imm = 0x3F800000
movl %edi, -0x10(%rdx)
movupd %xmm0, -0x8(%rdx)
leaq -0x1f0(%rbp), %rdx
movl %eax, (%rdx)
movq %rax, 0x8(%rdx)
movq %rdx, 0x10(%rdx)
movq %rdx, 0x18(%rdx)
movq %rax, 0x20(%rdx)
leaq -0x1b8(%rbp), %rdx
movq %rcx, -0x10(%rdx)
movl %eax, (%rdx)
movq %rax, 0x8(%rdx)
movq %rdx, 0x10(%rdx)
movq %rdx, 0x18(%rdx)
movq %rax, 0x20(%rdx)
leaq -0x160(%rbp), %rdx
movq %rdx, -0x30(%rdx)
movq %rsi, -0x28(%rdx)
movupd %xmm0, -0x20(%rdx)
movl %edi, -0x10(%rdx)
movupd %xmm0, -0x8(%rdx)
leaq -0x150(%rbp), %rdx
movl %eax, (%rdx)
movq %rax, 0x8(%rdx)
movq %rdx, 0x10(%rdx)
movq %rdx, 0x18(%rdx)
movq %rax, 0x20(%rdx)
movq %rcx, -0x80(%rbp)
leaq -0x268(%rbp), %rdx
movq %rdx, -0x78(%rbp)
leaq 0xf2e5(%rip), %rdx      # 0x22be8
movq %rdx, -0x88(%rbp)
movq %rax, -0x70(%rbp)
movl %eax, -0x68(%rbp)
movq %rcx, -0xa8(%rbp)
leaq -0x1c8(%rbp), %rcx
movq %rcx, -0xa0(%rbp)
movq %rdx, -0xb0(%rbp)
movq %rax, -0x98(%rbp)
movl %eax, -0x90(%rbp)
movq %rax, -0x30(%rbp)
movapd %xmm0, -0x40(%rbp)
leaq -0x120(%rbp), %rcx
movl %eax, (%rcx)
movq %rax, 0x8(%rcx)
movq %rcx, 0x10(%rcx)
movq %rcx, 0x18(%rcx)
movq %rax, 0x20(%rcx)
movq 0x60(%r15), %r12
cmpq 0x68(%r15), %r12
je 0x13a12
leaq -0xb0(%rbp), %r14
movq (%r12), %rsi
movq %r14, %rdi
callq 0x13340
movq -0x38(%rbp), %rsi
cmpq -0x30(%rbp), %rsi
je 0x13998
movupd (%r12), %xmm0
movupd %xmm0, (%rsi)
addq $0x10, -0x38(%rbp)
jmp 0x139a4
leaq -0x40(%rbp), %rdi
movq %r12, %rdx
callq 0x140be
movq 0x8(%rbx), %r13
movq (%r12), %rsi
movq %r13, %rdi
callq 0xcd08
movq %r13, %rdi
movq %rax, %rsi
callq 0xcdbe
cmpq $0x2, %rax
jb 0x13a04
movq 0x8(%rbx), %rdi
movq (%r12), %rsi
callq 0xcd08
movq %rax, -0xf8(%rbp)
leaq -0x128(%rbp), %rdi
leaq -0xf8(%rbp), %rsi
callq 0x141a0
movq 0x8(%rbx), %rdi
movq (%r12), %rsi
callq 0xcd08
movq %r14, %rdi
movq %rax, %rsi
callq 0x13474
addq $0x10, %r12
cmpq 0x68(%r15), %r12
jne 0x13971
movq -0x40(%rbp), %r14
movq -0x38(%rbp), %r15
cmpq %r15, %r14
je 0x13a5a
movq %r15, %rax
subq %r14, %rax
sarq $0x4, %rax
bsrq %rax, %rdx
xorl $0x3f, %edx
addl %edx, %edx
xorq $0x7e, %rdx
leaq -0xfa5(%rip), %rcx      # 0x12a98
movq %r14, %rdi
movq %r15, %rsi
callq 0x1425a
leaq -0xfb7(%rip), %rdx      # 0x12a98
movq %r14, %rdi
movq %r15, %rsi
callq 0x14392
leaq -0xf8(%rbp), %rcx
movq $-0x1, 0x8(%rcx)
movabsq $-0x4010000000000000, %rax # imm = 0xBFF0000000000000
movq %rax, 0x10(%rcx)
leaq 0xed3a(%rip), %rax      # 0x227b8
movq %rax, (%rcx)
xorpd %xmm0, %xmm0
movupd %xmm0, 0x18(%rcx)
movabsq $0x3fe999999999999a, %rax # imm = 0x3FE999999999999A
movq %rax, 0x28(%rcx)
movq %rax, 0x30(%rcx)
movq %rax, 0x38(%rcx)
movabsq $0x3ff999999999999a, %rax # imm = 0x3FF999999999999A
movq %rax, 0x40(%rcx)
movq -0x40(%rbp), %r15
leaq -0xb0(%rbp), %r12
cmpq -0x38(%rbp), %r15
je 0x13bd8
movsd 0x8(%r15), %xmm0
movsd %xmm0, -0x48(%rbp)
cmpq -0x38(%rbp), %r15
je 0x13bb1
movsd 0x8(%r15), %xmm0
ucomisd -0x48(%rbp), %xmm0
jne 0x13bb1
jp 0x13bb1
movq 0x8(%rbx), %rdi
movq (%r15), %rsi
callq 0xcd08
movq %rax, %r14
movq 0x8(%rbx), %rdi
movq %rax, %rsi
callq 0xcdbe
movq %rax, %r13
movq (%r15), %rsi
leaq -0x88(%rbp), %rdi
callq 0x13340
movq (%r15), %rsi
movq %r12, %rdi
callq 0x133da
movq -0xa0(%rbp), %rax
movq 0x80(%rax), %rdx
addq $0x78, %rax
movq %rax, %rcx
testq %rdx, %rdx
je 0x13b5c
movq %rax, %rcx
xorl %esi, %esi
cmpq %r14, 0x20(%rdx)
setb %sil
cmovaeq %rdx, %rcx
movq 0x10(%rdx,%rsi,8), %rdx
testq %rdx, %rdx
jne 0x13b44
movq %rax, %rdx
cmpq %rax, %rcx
je 0x13b6f
cmpq %r14, 0x20(%rcx)
cmovaq %rax, %rcx
movq %rcx, %rdx
cmpq $0x2, %r13
jb 0x13ba8
cmpq %rax, %rdx
je 0x13ba8
movq -0x78(%rbp), %rdi
movq %r14, %rsi
callq 0x1243a
addq %rax, %rax
cmpq %r13, %rax
jbe 0x13ba8
leaq -0x88(%rbp), %rdi
movq %r14, %rsi
callq 0x13474
movq %r12, %rdi
movq %r14, %rsi
callq 0x1355e
addq $0x10, %r15
jmp 0x13ace
movq %rbx, %rdi
leaq -0x88(%rbp), %rsi
movq %r12, %rdx
movq -0x50(%rbp), %rcx
movsd -0x48(%rbp), %xmm0
leaq -0xf8(%rbp), %r8
callq 0x13cfe
jmp 0x13ab9
leaq 0xb8(%rbx), %rax
leaq -0x60(%rbp), %rdi
movq %rax, (%rdi)
movb $0x0, 0x8(%rdi)
callq 0x14046
movb $0x1, -0x58(%rbp)
movq 0xa8(%rbx), %rsi
leaq -0xf8(%rbp), %rdi
callq 0x9af6
testb %al, %al
je 0x13c1e
movq 0xa8(%rbx), %rdi
movq (%rdi), %rax
leaq -0xf8(%rbp), %rsi
callq *0x8(%rax)
movq 0xe0(%rbx), %rdi
callq 0x13ea6
leaq -0x60(%rbp), %rdi
callq 0x13eee
leaq -0x128(%rbp), %rdi
callq 0x9c80
movq -0x40(%rbp), %rdi
testq %rdi, %rdi
je 0x13c4d
callq 0x5340
leaq -0x158(%rbp), %rdi
callq 0x9c80
leaq -0x190(%rbp), %rdi
callq 0x9cd8
leaq -0x1c0(%rbp), %rdi
callq 0x9d4e
leaq -0x1f8(%rbp), %rdi
callq 0x9c80
leaq -0x230(%rbp), %rdi
callq 0x9cd8
leaq -0x260(%rbp), %rdi
callq 0x9d4e
addq $0x248, %rsp            # imm = 0x248
popq %rbx
popq %r12
popq %r13
popq %r14
popq %r15
popq %rbp
retq
jmp 0x13cc1
jmp 0x13cc1
movq %rax, %rbx
leaq -0x60(%rbp), %rdi
callq 0x13eee
jmp 0x13cc4
jmp 0x13cc1
jmp 0x13cc1
jmp 0x13cc1
jmp 0x13cc1
movq %rax, %rbx
leaq -0x128(%rbp), %rdi
callq 0x9c80
movq -0x40(%rbp), %rdi
testq %rdi, %rdi
je 0x13cde
callq 0x5340
leaq -0x1c8(%rbp), %rdi
callq 0x13f08
leaq -0x268(%rbp), %rdi
callq 0x13f08
movq %rbx, %rdi
callq 0x56a0