// File: jingdao[P]PointCloudApp/kdtree.cpp
// Function: KdTree::searchNode(KdTree::KdTreeNode*, std::vector<int, std::allocator<int>>*, KdTree::Cube*, KdTree::Cube, int)
void KdTree::searchNode(KdTreeNode* node,std::vector<int> *indices,Cube* range,Cube cell,int depth){
	if (!node) return; //nothing here
	if (node->isLeaf) {
		float x = float_data[node->index * 4];
		float y = float_data[node->index * 4 + 1];
		float z = float_data[node->index * 4 + 2];
		if (x >= range->x1 && x <= range->x2 && 
			y >= range->y1 && y <= range->y2 && 
			z >= range->z1 && z <= range->z2)
			indices->push_back(node->index);
	} else {
		if (cell.x1 >= range->x1 && cell.x2 <= range->x2 &&
			cell.y1 >= range->y1 && cell.y2 <= range->y2 &&
			cell.z1 >= range->z1 && cell.z2 <= range->z2) {
			//region fully contained in range
			getIndexFromSubTree(node->left,indices);
			getIndexFromSubTree(node->right,indices);
		} else if ((cell.x2 >= range->x1 && cell.x1 <= range->x2) &&
			(cell.y2 >= range->y1 && cell.y1 <= range->y2) &&
			(cell.z2 >= range->z1 && cell.z1 <= range->z2)) {
			//region intersects range
			int dimension = depth % 3;
			Cube leftCell = cell;
			Cube rightCell = cell;
			//trim cells according to splitting plane
			switch (dimension) {
				case 0:
					leftCell.x2 = node->value;
					rightCell.x1 = node->value;
					break;
				case 1:
					leftCell.y2 = node->value;
					rightCell.y1 = node->value;
					break;
				case 2:
					leftCell.z2 = node->value;
					rightCell.z1 = node->value;
					break;
			}
			searchNode(node->left,indices,range,leftCell,depth+1);
			searchNode(node->right,indices,range,rightCell,depth+1);
		}
	}
}

// Assembly:
pushq %rbp
movq %rsp, %rbp
subq $0xd0, %rsp
leaq 0x10(%rbp), %rax
movq %rax, -0xb0(%rbp)
movq %rdi, -0x8(%rbp)
movq %rsi, -0x10(%rbp)
movq %rdx, -0x18(%rbp)
movq %rcx, -0x20(%rbp)
movl %r8d, -0x24(%rbp)
movq -0x8(%rbp), %rax
movq %rax, -0xa8(%rbp)
cmpq $0x0, -0x10(%rbp)
jne 0xe771
jmp 0xeb5a
movq -0x10(%rbp), %rax
testb $0x1, (%rax)
je 0xe857
movq -0xa8(%rbp), %rax
movq 0x30(%rax), %rcx
movq -0x10(%rbp), %rdx
movl 0x4(%rdx), %edx
shll $0x2, %edx
movslq %edx, %rdx
movss (%rcx,%rdx,4), %xmm0
movss %xmm0, -0x28(%rbp)
movq 0x30(%rax), %rcx
movq -0x10(%rbp), %rdx
movl 0x4(%rdx), %edx
shll $0x2, %edx
addl $0x1, %edx
movslq %edx, %rdx
movss (%rcx,%rdx,4), %xmm0
movss %xmm0, -0x2c(%rbp)
movq 0x30(%rax), %rax
movq -0x10(%rbp), %rcx
movl 0x4(%rcx), %ecx
shll $0x2, %ecx
addl $0x2, %ecx
movslq %ecx, %rcx
movss (%rax,%rcx,4), %xmm0
movss %xmm0, -0x30(%rbp)
movss -0x28(%rbp), %xmm0
movq -0x20(%rbp), %rax
ucomiss (%rax), %xmm0
jb 0xe852
movss -0x28(%rbp), %xmm1
movq -0x20(%rbp), %rax
movss 0xc(%rax), %xmm0
ucomiss %xmm1, %xmm0
jb 0xe852
movss -0x2c(%rbp), %xmm0
movq -0x20(%rbp), %rax
ucomiss 0x4(%rax), %xmm0
jb 0xe852
movss -0x2c(%rbp), %xmm1
movq -0x20(%rbp), %rax
movss 0x10(%rax), %xmm0
ucomiss %xmm1, %xmm0
jb 0xe852
movss -0x30(%rbp), %xmm0
movq -0x20(%rbp), %rax
ucomiss 0x8(%rax), %xmm0
jb 0xe852
movss -0x30(%rbp), %xmm1
movq -0x20(%rbp), %rax
movss 0x14(%rax), %xmm0
ucomiss %xmm1, %xmm0
jb 0xe852
movq -0x18(%rbp), %rdi
movq -0x10(%rbp), %rsi
addq $0x4, %rsi
callq 0xae30
jmp 0xeb5a
movq -0xb0(%rbp), %rax
movss (%rax), %xmm0
movq -0x20(%rbp), %rax
ucomiss (%rax), %xmm0
jb 0xe922
movq -0xb0(%rbp), %rax
movss 0xc(%rax), %xmm1
movq -0x20(%rbp), %rax
movss 0xc(%rax), %xmm0
ucomiss %xmm1, %xmm0
jb 0xe922
movq -0xb0(%rbp), %rax
movss 0x4(%rax), %xmm0
movq -0x20(%rbp), %rax
ucomiss 0x4(%rax), %xmm0
jb 0xe922
movq -0xb0(%rbp), %rax
movss 0x10(%rax), %xmm1
movq -0x20(%rbp), %rax
movss 0x10(%rax), %xmm0
ucomiss %xmm1, %xmm0
jb 0xe922
movq -0xb0(%rbp), %rax
movss 0x8(%rax), %xmm0
movq -0x20(%rbp), %rax
ucomiss 0x8(%rax), %xmm0
jb 0xe922
movq -0xb0(%rbp), %rax
movss 0x14(%rax), %xmm1
movq -0x20(%rbp), %rax
movss 0x14(%rax), %xmm0
ucomiss %xmm1, %xmm0
jb 0xe922
movq -0xa8(%rbp), %rdi
movq -0x10(%rbp), %rax
movq 0x8(%rax), %rsi
movq -0x18(%rbp), %rdx
callq 0xeb70
movq -0xa8(%rbp), %rdi
movq -0x10(%rbp), %rax
movq 0x10(%rax), %rsi
movq -0x18(%rbp), %rdx
callq 0xeb70
jmp 0xeb58
movq -0xb0(%rbp), %rax
movss 0xc(%rax), %xmm0
movq -0x20(%rbp), %rax
ucomiss (%rax), %xmm0
jb 0xeb56
movq -0xb0(%rbp), %rax
movss (%rax), %xmm1
movq -0x20(%rbp), %rax
movss 0xc(%rax), %xmm0
ucomiss %xmm1, %xmm0
jb 0xeb56
movq -0xb0(%rbp), %rax
movss 0x10(%rax), %xmm0
movq -0x20(%rbp), %rax
ucomiss 0x4(%rax), %xmm0
jb 0xeb56
movq -0xb0(%rbp), %rax
movss 0x4(%rax), %xmm1
movq -0x20(%rbp), %rax
movss 0x10(%rax), %xmm0
ucomiss %xmm1, %xmm0
jb 0xeb56
movq -0xb0(%rbp), %rax
movss 0x14(%rax), %xmm0
movq -0x20(%rbp), %rax
ucomiss 0x8(%rax), %xmm0
jb 0xeb56
movq -0xb0(%rbp), %rax
movss 0x8(%rax), %xmm1
movq -0x20(%rbp), %rax
movss 0x14(%rax), %xmm0
ucomiss %xmm1, %xmm0
jb 0xeb56
movq -0xb0(%rbp), %rax
movslq -0x24(%rbp), %rdx
movl %edx, %ecx
imulq $0x55555556, %rdx, %rsi # imm = 0x55555556
movq %rsi, %rdx
shrq $0x3f, %rdx
shrq $0x20, %rsi
addl %edx, %esi
movl %esi, %edx
leal (%rdx,%rdx,2), %edx
subl %edx, %ecx
movl %ecx, -0x34(%rbp)
movq 0x10(%rax), %rcx
movq %rcx, -0x40(%rbp)
movups (%rax), %xmm0
movaps %xmm0, -0x50(%rbp)
movq 0x10(%rax), %rcx
movq %rcx, -0x60(%rbp)
movups (%rax), %xmm0
movaps %xmm0, -0x70(%rbp)
movl -0x34(%rbp), %eax
movl %eax, -0xb4(%rbp)
testl %eax, %eax
je 0xea3a
jmp 0xea20
movl -0xb4(%rbp), %eax
subl $0x1, %eax
je 0xea58
jmp 0xea2d
movl -0xb4(%rbp), %eax
subl $0x2, %eax
je 0xea76
jmp 0xea92
movq -0x10(%rbp), %rax
movss 0x4(%rax), %xmm0
movss %xmm0, -0x44(%rbp)
movq -0x10(%rbp), %rax
movss 0x4(%rax), %xmm0
movss %xmm0, -0x70(%rbp)
jmp 0xea92
movq -0x10(%rbp), %rax
movss 0x4(%rax), %xmm0
movss %xmm0, -0x40(%rbp)
movq -0x10(%rbp), %rax
movss 0x4(%rax), %xmm0
movss %xmm0, -0x6c(%rbp)
jmp 0xea92
movq -0x10(%rbp), %rax
movss 0x4(%rax), %xmm0
movss %xmm0, -0x3c(%rbp)
movq -0x10(%rbp), %rax
movss 0x4(%rax), %xmm0
movss %xmm0, -0x68(%rbp)
movq -0xa8(%rbp), %rdi
movq -0x10(%rbp), %rax
movq 0x8(%rax), %rsi
movq -0x18(%rbp), %rdx
movq -0x20(%rbp), %rcx
movq -0x50(%rbp), %rax
movq %rax, -0x88(%rbp)
movq -0x48(%rbp), %rax
movq %rax, -0x80(%rbp)
movq -0x40(%rbp), %rax
movq %rax, -0x78(%rbp)
movl -0x24(%rbp), %r8d
addl $0x1, %r8d
leaq -0x88(%rbp), %rax
movq (%rax), %r9
movq %r9, (%rsp)
movq 0x8(%rax), %r9
movq %r9, 0x8(%rsp)
movq 0x10(%rax), %rax
movq %rax, 0x10(%rsp)
callq 0xe730
movq -0xa8(%rbp), %rdi
movq -0x10(%rbp), %rax
movq 0x10(%rax), %rsi
movq -0x18(%rbp), %rdx
movq -0x20(%rbp), %rcx
movq -0x70(%rbp), %rax
movq %rax, -0xa0(%rbp)
movq -0x68(%rbp), %rax
movq %rax, -0x98(%rbp)
movq -0x60(%rbp), %rax
movq %rax, -0x90(%rbp)
movl -0x24(%rbp), %r8d
addl $0x1, %r8d
leaq -0xa0(%rbp), %rax
movq (%rax), %r9
movq %r9, (%rsp)
movq 0x8(%rax), %r9
movq %r9, 0x8(%rsp)
movq 0x10(%rax), %rax
movq %rax, 0x10(%rsp)
callq 0xe730
jmp 0xeb58
jmp 0xeb5a
addq $0xd0, %rsp
popq %rbp
retq
nopw %cs:(%rax,%rax)