// File: /workspace/github_cpp_makefile_star2_O2/CK7[P]abawaca/src/ClusterSeparatorSplitScafs.cpp
// Function: ClusterSeparatorSplitScafs::separate_dimension(unsigned long)
void ClusterSeparatorSplitScafs::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);
	ClusterSeparatorSplitScafs::IsScafIn		in_operator(fraction_dps_to_be_considered_in);
	ClusterSeparatorSplitScafs::ClusterStats	curr_cluster_stats1(curr_cluster1, scaf_db, in_operator);
	ClusterSeparatorSplitScafs::ClusterStats	curr_cluster_stats2(curr_cluster2, scaf_db, in_operator);
	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 += vit->dp;
		dps.push_back(*vit);
		if(scaf_db.ndps(scaf_db.dp2scaf(vit->dp)) > 1)
			useful_scafs.insert(scaf_db.dp2scaf(vit->dp));
	}

	if(curr_cluster_stats2.num_scafs_in_this_cluster() != useful_scafs.size()) {
			cerr << "Error: # useful scafs=" << useful_scafs.size() << " but in cluster 2: " << curr_cluster_stats2.num_scafs_in_this_cluster() << endl;
		semaphore.notify();
		return;
	}
	sort(dps.begin(), dps.end(), comp_by_value);
	ClusterSeparatorSplitScafs::ClusteringResult	best_dimension_separation;

	for(auto vit=dps.begin(); vit!=dps.end();) {
		double value = vit->value;
		while(vit!=dps.end() && vit->value==value) {
			curr_cluster_stats1 += vit->dp;
			curr_cluster_stats2 -= vit->dp;
			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, useful_scafs.size(), 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 $0x238, %rsp            # imm = 0x238
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 -0x248(%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 -0x1f0(%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 -0x1e0(%rbp), %rdx
andl $0x0, (%rdx)
andq $0x0, 0x8(%rdx)
movq %rdx, 0x10(%rdx)
movq %rdx, 0x18(%rdx)
andq $0x0, 0x20(%rdx)
leaq -0x1a8(%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)
movq %rax, %r15
leaq -0x150(%rbp), %rax
movq %rax, -0x30(%rax)
movq %rsi, -0x28(%rax)
movupd %xmm0, -0x20(%rax)
movl %edi, -0x10(%rax)
movupd %xmm0, -0x8(%rax)
leaq -0x140(%rbp), %rax
andl $0x0, (%rax)
andq $0x0, 0x8(%rax)
movq %rax, 0x10(%rax)
movq %rax, 0x18(%rax)
andq $0x0, 0x20(%rax)
movsd 0xe8(%rbx), %xmm1
leaq -0xc0(%rbp), %rax
movq %rcx, -0x88(%rbp)
leaq -0x258(%rbp), %rdx
movq %rdx, -0x80(%rbp)
andq $0x0, -0x78(%rbp)
movsd %xmm1, (%rax)
andl $0x0, -0x70(%rbp)
movq %rax, -0x68(%rbp)
leaq -0x118(%rbp), %r12
movq %rcx, (%r12)
leaq -0x1b8(%rbp), %rcx
movq %rcx, 0x8(%r12)
andq $0x0, 0x10(%r12)
andl $0x0, 0x18(%r12)
movq %rax, 0x20(%r12)
leaq -0x40(%rbp), %rcx
andq $0x0, 0x10(%rcx)
leaq -0xb0(%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(%r15), %r14
cmpq 0x68(%r15), %r14
je 0x11129
movq (%r14), %rsi
movq %r12, %rdi
callq 0x10b44
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 0x11123
movq 0x8(%rbx), %rdi
movq (%r14), %rsi
callq 0x9b48
movq %rax, -0xf0(%rbp)
leaq -0xb8(%rbp), %rdi
leaq -0xf0(%rbp), %rsi
callq 0x103b6
addq $0x10, %r14
jmp 0x110c0
movslq -0x100(%rbp), %rax
cmpq %rax, -0x90(%rbp)
jne 0x11218
movq -0x40(%rbp), %rdi
movq -0x38(%rbp), %rsi
leaq -0x215c(%rip), %rdx     # 0xeff0
callq 0x1048c
leaq -0xf0(%rbp), %r14
orq $-0x1, 0x8(%r14)
movabsq $-0x4010000000000000, %rax # imm = 0xBFF0000000000000
movq %rax, 0x10(%r14)
leaq 0xdb6e(%rip), %rax      # 0x1ece0
movq %rax, (%r14)
movabsq $0x412e848000000000, %rax # imm = 0x412E848000000000
movq %rax, 0x18(%r14)
movq %rax, 0x20(%r14)
movabsq $0x3fb999999999999a, %rax # imm = 0x3FB999999999999A
movq %rax, 0x28(%r14)
movq -0x40(%rbp), %r13
leaq -0x88(%rbp), %r15
leaq -0x118(%rbp), %r12
movq -0x38(%rbp), %rax
cmpq %rax, %r13
je 0x1126d
movsd 0x8(%r13), %xmm0
movsd %xmm0, -0x48(%rbp)
cmpq %rax, %r13
je 0x111f5
movsd 0x8(%r13), %xmm0
ucomisd -0x48(%rbp), %xmm0
jne 0x111f5
jp 0x111f5
movq (%r13), %rsi
movq %r15, %rdi
callq 0x10b44
movq (%r13), %rsi
movq %r12, %rdi
callq 0x10cb0
addq $0x10, %r13
movq -0x38(%rbp), %rax
jmp 0x111bf
movq -0x90(%rbp), %r8
movq %rbx, %rdi
movq %r15, %rsi
movq %r12, %rdx
movq -0x50(%rbp), %rcx
movsd -0x48(%rbp), %xmm0
movq %r14, %r9
callq 0x1134c
jmp 0x111a7
movq 0xddc9(%rip), %rdi      # 0x1efe8
leaq 0x5142(%rip), %rsi      # 0x16368
callq 0x4340
movq -0x90(%rbp), %rsi
movq %rax, %rdi
callq 0x4240
leaq 0x513e(%rip), %rsi      # 0x1637f
movq %rax, %rdi
callq 0x4340
movl -0x100(%rbp), %esi
movq %rax, %rdi
callq 0x4630
movq %rax, %rdi
callq 0x41a0
movq 0xe0(%rbx), %rdi
callq 0x1008a
jmp 0x112bd
leaq 0xb8(%rbx), %rsi
leaq -0x60(%rbp), %rdi
callq 0x1006a
movq 0xa8(%rbx), %rsi
leaq -0xf0(%rbp), %rdi
callq 0x1146e
testb %al, %al
je 0x112a8
movq 0xa8(%rbx), %rdi
movq (%rdi), %rax
leaq -0xf0(%rbp), %rsi
callq *0x8(%rax)
movq 0xe0(%rbx), %rdi
callq 0x1008a
leaq -0x60(%rbp), %rdi
callq 0x100c4
leaq -0xb8(%rbp), %rdi
callq 0x7602
leaq -0x40(%rbp), %rdi
callq 0xc57e
leaq -0x1b8(%rbp), %rdi
callq 0x6d3e
leaq -0x258(%rbp), %rdi
callq 0x6d3e
addq $0x238, %rsp            # imm = 0x238
popq %rbx
popq %r12
popq %r13
popq %r14
popq %r15
popq %rbp
retq
jmp 0x11314
movq %rax, %rbx
leaq -0x60(%rbp), %rdi
callq 0x100c4
jmp 0x11317
jmp 0x11314
jmp 0x11314
jmp 0x11314
jmp 0x11314
movq %rax, %rbx
leaq -0xb8(%rbp), %rdi
callq 0x7602
leaq -0x40(%rbp), %rdi
callq 0xc57e
leaq -0x1b8(%rbp), %rdi
callq 0x6d3e
leaq -0x258(%rbp), %rdi
callq 0x6d3e
movq %rbx, %rdi
callq 0x4660