// File: /workspace/github_cpp_makefile_star2_O2/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 0xfffc
movq 0x18(%rbx), %rdi
movq %r14, -0x50(%rbp)
movq %r14, %rsi
callq 0xc9d2
movq 0x8(%rbx), %rcx
leaq -0x258(%rbp), %rdx
movq %rcx, -0x10(%rdx)
andl $0x0, (%rdx)
andq $0x0, 0x8(%rdx)
movq %rdx, 0x10(%rdx)
movq %rdx, 0x18(%rdx)
andq $0x0, 0x20(%rdx)
leaq -0x200(%rbp), %rdx
movq %rdx, -0x30(%rdx)
pushq $0x1
popq %rsi
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
andl $0x0, (%rdx)
andq $0x0, 0x8(%rdx)
movq %rdx, 0x10(%rdx)
movq %rdx, 0x18(%rdx)
andq $0x0, 0x20(%rdx)
leaq -0x1b8(%rbp), %rdx
movq %rcx, -0x10(%rdx)
andl $0x0, (%rdx)
andq $0x0, 0x8(%rdx)
movq %rdx, 0x10(%rdx)
movq %rdx, 0x18(%rdx)
andq $0x0, 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
andl $0x0, (%rdx)
andq $0x0, 0x8(%rdx)
movq %rdx, 0x10(%rdx)
movq %rdx, 0x18(%rdx)
andq $0x0, 0x20(%rdx)
movq %rcx, -0x70(%rbp)
leaq -0x268(%rbp), %rdx
movq %rdx, -0x68(%rbp)
leaq 0xf0ea(%rip), %rdx      # 0x1ec28
movq %rdx, -0x78(%rbp)
andq $0x0, -0x60(%rbp)
andl $0x0, -0x58(%rbp)
leaq -0x128(%rbp), %r15
movq %rcx, 0x8(%r15)
leaq -0x1c8(%rbp), %rcx
movq %rcx, 0x10(%r15)
movq %rdx, (%r15)
andq $0x0, 0x18(%r15)
movq %rax, %r12
andl $0x0, 0x20(%r15)
leaq -0x40(%rbp), %rcx
andq $0x0, 0x10(%rcx)
leaq -0xf8(%rbp), %rax
andl $0x0, (%rax)
andq $0x0, 0x8(%rax)
movq %rax, 0x10(%rax)
movq %rax, 0x18(%rax)
andq $0x0, 0x20(%rax)
movapd %xmm0, (%rcx)
movq 0x60(%r12), %r14
cmpq 0x68(%r12), %r14
je 0xfc23
movq (%r14), %rsi
movq %r15, %rdi
callq 0xf740
leaq -0x40(%rbp), %rdi
movq %r14, %rsi
callq 0x10046
movq 0x8(%rbx), %r13
movq (%r14), %rsi
movq %r13, %rdi
callq 0x9b48
movq %r13, %rdi
movq %rax, %rsi
callq 0x9bec
cmpq $0x2, %rax
jb 0xfc1a
movq 0x8(%rbx), %rdi
movq (%r14), %rsi
callq 0x9b48
movq %rax, -0xd0(%rbp)
leaq -0x100(%rbp), %rdi
leaq -0xd0(%rbp), %rsi
callq 0x103b6
movq 0x8(%rbx), %rdi
movq (%r14), %rsi
callq 0x9b48
movq %r15, %rdi
movq %rax, %rsi
callq 0xf7fc
addq $0x10, %r14
jmp 0xfb9f
movq -0x40(%rbp), %rdi
movq -0x38(%rbp), %rsi
leaq -0xc42(%rip), %rdx      # 0xeff0
callq 0x1048c
leaq -0xd0(%rbp), %rcx
orq $-0x1, 0x8(%rcx)
movabsq $-0x4010000000000000, %rax # imm = 0xBFF0000000000000
movq %rax, 0x10(%rcx)
leaq 0xebd8(%rip), %rax      # 0x1e830
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), %r12
leaq -0x128(%rbp), %r13
movq -0x38(%rbp), %rax
cmpq %rax, %r12
je 0xfd7d
movsd 0x8(%r12), %xmm0
movsd %xmm0, -0x48(%rbp)
cmpq %rax, %r12
je 0xfd59
movsd 0x8(%r12), %xmm0
ucomisd -0x48(%rbp), %xmm0
jne 0xfd59
jp 0xfd59
movq 0x8(%rbx), %rdi
movq (%r12), %rsi
callq 0x9b48
movq %rax, %r14
movq 0x8(%rbx), %rdi
movq %rax, %rsi
callq 0x9bec
movq %rax, %r15
movq (%r12), %rsi
leaq -0x78(%rbp), %rdi
callq 0xf740
movq (%r12), %rsi
movq %r13, %rdi
callq 0xf79e
movq -0x118(%rbp), %rdi
movq %r14, %rsi
callq 0xf638
cmpq $0x2, %r15
setb %cl
xorb $0x1, %al
orb %cl, %al
jne 0xfd4c
movq -0x68(%rbp), %rdi
movq %r14, %rsi
callq 0xea14
addq %rax, %rax
cmpq %r15, %rax
jbe 0xfd4c
leaq -0x78(%rbp), %rdi
movq %r14, %rsi
callq 0xf7fc
movq %r13, %rdi
movq %r14, %rsi
callq 0xf872
addq $0x10, %r12
movq -0x38(%rbp), %rax
jmp 0xfcac
movq %rbx, %rdi
leaq -0x78(%rbp), %rsi
movq %r13, %rdx
movq -0x50(%rbp), %rcx
movsd -0x48(%rbp), %xmm0
leaq -0xd0(%rbp), %r8
callq 0xfe68
jmp 0xfc93
leaq 0xb8(%rbx), %rsi
leaq -0x88(%rbp), %rdi
callq 0x1006a
movq 0xa8(%rbx), %rsi
leaq -0xd0(%rbp), %rdi
callq 0x749a
testb %al, %al
je 0xfdbb
movq 0xa8(%rbx), %rdi
movq (%rdi), %rax
leaq -0xd0(%rbp), %rsi
callq *0x8(%rax)
movq 0xe0(%rbx), %rdi
callq 0x1008a
leaq -0x88(%rbp), %rdi
callq 0x100c4
leaq -0x100(%rbp), %rdi
callq 0x7602
leaq -0x40(%rbp), %rdi
callq 0xc57e
leaq -0x1c8(%rbp), %rdi
callq 0x6d3e
leaq -0x268(%rbp), %rdi
callq 0x6d3e
addq $0x248, %rsp            # imm = 0x248
popq %rbx
popq %r12
popq %r13
popq %r14
popq %r15
popq %rbp
retq
jmp 0xfe2f
jmp 0xfe2f
movq %rax, %rbx
leaq -0x88(%rbp), %rdi
callq 0x100c4
jmp 0xfe32
jmp 0xfe2f
jmp 0xfe2f
jmp 0xfe2f
jmp 0xfe2f
movq %rax, %rbx
leaq -0x100(%rbp), %rdi
callq 0x7602
leaq -0x40(%rbp), %rdi
callq 0xc57e
leaq -0x1c8(%rbp), %rdi
callq 0x6d3e
leaq -0x268(%rbp), %rdi
callq 0x6d3e
movq %rbx, %rdi
callq 0x4660
nop