// File: /workspace/github_cpp_makefile_star2_O3/Luthaf[P]zeopp/voro++/src/v_compute.cc
// Function: voro::voro_compute<voro::container_poly>::find_voronoi_cell(double, double, double, int, int, int, int, voro::particle_record&, double&)
void voro_compute<c_class>::find_voronoi_cell(double x,double y,double z,int ci,int cj,int ck,int ijk,particle_record &w,double &mrs) {
	double qx=0,qy=0,qz=0,rs;
	int i,j,k,di,dj,dk,ei,ej,ek,f,g,disp;
	double fx,fy,fz,mxs,mys,mzs,*radp;
	unsigned int q,*e,*mijk;

	// Init setup for parameters to return
	w.ijk=-1;mrs=large_number;

	con.initialize_search(ci,cj,ck,ijk,i,j,k,disp);

	// Test all particles in the particle's local region first
	scan_all(ijk,x,y,z,0,0,0,w,mrs);

	// Now compute the fractional position of the particle within its
	// region and store it in (fx,fy,fz). We use this to compute an index
	// (di,dj,dk) of which subregion the particle is within.
	unsigned int m1,m2;
	con.frac_pos(x,y,z,ci,cj,ck,fx,fy,fz);
	di=int(fx*xsp*wl_fgrid);dj=int(fy*ysp*wl_fgrid);dk=int(fz*zsp*wl_fgrid);

	// The indices (di,dj,dk) tell us which worklist to use, to test the
	// blocks in the optimal order. But we only store worklists for the
	// eighth of the region where di, dj, and dk are all less than half the
	// full grid. The rest of the cases are handled by symmetry. In this
	// section, we detect for these cases, by reflecting high values of di,
	// dj, and dk. For these cases, a mask is constructed in m1 and m2
	// which is used to flip the worklist information when it is loaded.
	if(di>=wl_hgrid) {
		mxs=boxx-fx;
		m1=127+(3<<21);m2=1+(1<<21);di=wl_fgrid-1-di;if(di<0) di=0;
	} else {m1=m2=0;mxs=fx;}
	if(dj>=wl_hgrid) {
		mys=boxy-fy;
		m1|=(127<<7)+(3<<24);m2|=(1<<7)+(1<<24);dj=wl_fgrid-1-dj;if(dj<0) dj=0;
	} else mys=fy;
	if(dk>=wl_hgrid) {
		mzs=boxz-fz;
		m1|=(127<<14)+(3<<27);m2|=(1<<14)+(1<<27);dk=wl_fgrid-1-dk;if(dk<0) dk=0;
	} else mzs=fz;

	// Do a quick test to account for the case when the minimum radius is
	// small enought that no other blocks need to be considered
	rs=con.r_max_add(mrs);
	if(mxs*mxs>rs&&mys*mys>rs&&mzs*mzs>rs) return;

	// Now compute which worklist we are going to use, and set radp and e to
	// point at the right offsets
	ijk=di+wl_hgrid*(dj+wl_hgrid*dk);
	radp=mrad+ijk*wl_seq_length;
	e=(const_cast<unsigned int*> (wl))+ijk*wl_seq_length;

	// Read in how many items in the worklist can be tested without having to
	// worry about writing to the mask
	f=e[0];g=0;
	do {

		// If mrs is less than the minimum distance to any untested
		// block, then we are done
		if(con.r_max_add(mrs)<radp[g]) return;
		g++;

		// Load in a block off the worklist, permute it with the
		// symmetry mask, and decode its position. These are all
		// integer bit operations so they should run very fast.
		q=e[g];q^=m1;q+=m2;
		di=q&127;di-=64;
		dj=(q>>7)&127;dj-=64;
		dk=(q>>14)&127;dk-=64;

		// Check that the worklist position is in range
		ei=di+i;if(ei<0||ei>=hx) continue;
		ej=dj+j;if(ej<0||ej>=hy) continue;
		ek=dk+k;if(ek<0||ek>=hz) continue;

		// Call the compute_min_max_radius() function. This returns
		// true if the minimum distance to the block is bigger than the
		// current mrs, in which case we skip this block and move on.
		// Otherwise, it computes the maximum distance to the block and
		// returns it in crs.
		if(compute_min_radius(di,dj,dk,fx,fy,fz,mrs)) continue;

		// Now compute which region we are going to loop over, adding a
		// displacement for the periodic cases
		ijk=con.region_index(ci,cj,ck,ei,ej,ek,qx,qy,qz,disp);

		// If mrs is bigger than the maximum distance to the block,
		// then we have to test all particles in the block for
		// intersections. Otherwise, we do additional checks and skip
		// those particles which can't possibly intersect the block.
		scan_all(ijk,x-qx,y-qy,z-qz,di,dj,dk,w,mrs);
	} while(g<f);

	// Update mask value and initialize queue
	mv++;
	if(mv==0) {reset_mask();mv=1;}
	int *qu_s=qu,*qu_e=qu;

	while(g<wl_seq_length-1) {

		// If mrs is less than the minimum distance to any untested
		// block, then we are done
		if(con.r_max_add(mrs)<radp[g]) return;
		g++;

		// Load in a block off the worklist, permute it with the
		// symmetry mask, and decode its position. These are all
		// integer bit operations so they should run very fast.
		q=e[g];q^=m1;q+=m2;
		di=q&127;di-=64;
		dj=(q>>7)&127;dj-=64;
		dk=(q>>14)&127;dk-=64;

		// Compute the position in the mask of the current block. If
		// this lies outside the mask, then skip it. Otherwise, mark
		// it.
		ei=di+i;if(ei<0||ei>=hx) continue;
		ej=dj+j;if(ej<0||ej>=hy) continue;
		ek=dk+k;if(ek<0||ek>=hz) continue;
		mijk=mask+ei+hx*(ej+hy*ek);
		*mijk=mv;

		// Skip this block if it is further away than the current
		// minimum radius
		if(compute_min_radius(di,dj,dk,fx,fy,fz,mrs)) continue;

		// Now compute which region we are going to loop over, adding a
		// displacement for the periodic cases
		ijk=con.region_index(ci,cj,ck,ei,ej,ek,qx,qy,qz,disp);
		scan_all(ijk,x-qx,y-qy,z-qz,di,dj,dk,w,mrs);

		if(qu_e>qu_l-18) add_list_memory(qu_s,qu_e);
		scan_bits_mask_add(q,mijk,ei,ej,ek,qu_e);
	}

	// Do a check to see if we've reached the radius cutoff
	if(con.r_max_add(mrs)<radp[g]) return;

	// We were unable to completely compute the cell based on the blocks in
	// the worklist, so now we have to go block by block, reading in items
	// off the list
	while(qu_s!=qu_e) {

		// Read the next entry of the queue
		if(qu_s==qu_l) qu_s=qu;
		ei=*(qu_s++);ej=*(qu_s++);ek=*(qu_s++);
		di=ei-i;dj=ej-j;dk=ek-k;
		if(compute_min_radius(di,dj,dk,fx,fy,fz,mrs)) continue;

		ijk=con.region_index(ci,cj,ck,ei,ej,ek,qx,qy,qz,disp);
		scan_all(ijk,x-qx,y-qy,z-qz,di,dj,dk,w,mrs);

		// Test the neighbors of the current block, and add them to the
		// block list if they haven't already been tested
		if((qu_s<=qu_e?(qu_l-qu_e)+(qu_s-qu):qu_s-qu_e)<18) add_list_memory(qu_s,qu_e);
		add_to_mask(ei,ej,ek,qu_e);
	}
}

// Assembly:
pushq %rbp
pushq %r15
pushq %r14
pushq %r13
pushq %r12
pushq %rbx
subq $0x128, %rsp            # imm = 0x128
movq %rdi, %r10
movq 0x160(%rsp), %rdi
movl $0xffffffff, (%r9)      # imm = 0xFFFFFFFF
movabsq $0x7fefffffffffffff, %rax # imm = 0x7FEFFFFFFFFFFFFF
movq %rax, (%rdi)
movq (%r10), %rax
cmpb $0x0, 0xa8(%rax)
movl (%rax), %r11d
movq %rsi, 0xd0(%rsp)
movl %r11d, 0x30(%rsp)
cmovnel %r11d, %esi
movl %esi, 0x10(%rsp)
cmpb $0x0, 0xa9(%rax)
movl 0x4(%rax), %ebp
movq %rdx, 0xd8(%rsp)
cmovnel %ebp, %edx
movq 0x60(%r10), %rsi
cmpb $0x0, 0xaa(%rax)
movq %rcx, 0xe0(%rsp)
movl %ecx, 0x3c(%rsp)
je 0x21cd7
movl 0x8(%rax), %ecx
movl %ecx, 0x3c(%rsp)
movl %r8d, 0x38(%rsp)
movslq %r8d, %rdi
movl (%rsi,%rdi,4), %r14d
testl %r14d, %r14d
movq %r10, 0x20(%rsp)
jle 0x21dde
movq 0x58(%r10), %r8
movq (%r8,%rdi,8), %r8
movslq 0x4c(%r10), %r13
movq 0xd8(%rax), %r10
movq (%r10,%rdi,8), %r10
addq $0x18, %r10
addq $0x10, %r8
shlq $0x3, %r13
movsd 0x1761c(%rip), %xmm7    # 0x39338
xorl %r11d, %r11d
xorl %ebx, %ebx
movq 0x160(%rsp), %rcx
movslq %r14d, %r14
movq %r11, %r15
shlq $0x5, %r15
addq %r10, %r15
movq %r13, %r12
imulq %r11, %r12
addq %r8, %r12
movsd -0x10(%r12), %xmm4
movsd -0x8(%r12), %xmm3
subsd %xmm0, %xmm4
subsd %xmm1, %xmm3
movsd (%r12), %xmm5
subsd %xmm2, %xmm5
mulsd %xmm3, %xmm3
mulsd %xmm4, %xmm4
addsd %xmm3, %xmm4
mulsd %xmm5, %xmm5
addsd %xmm4, %xmm5
movsd (%r15), %xmm3
mulsd %xmm3, %xmm3
subsd %xmm3, %xmm5
ucomisd %xmm5, %xmm7
ja 0x21d98
incq %r11
addq $0x20, %r15
addq %r13, %r12
cmpq %r14, %r11
jl 0x21d40
jmp 0x21db8
movsd %xmm5, (%rcx)
movl %r11d, 0x4(%r9)
movslq (%rsi,%rdi,4), %r14
incq %r11
movb $0x1, %bl
movapd %xmm5, %xmm7
cmpq %r14, %r11
jl 0x21d29
jmp 0x21dc1
movapd %xmm7, %xmm5
testb $0x1, %bl
je 0x21dee
movl 0x38(%rsp), %esi
movl %esi, (%r9)
movq $0x0, 0x8(%r9)
movl $0x0, 0x10(%r9)
movapd %xmm5, %xmm7
jmp 0x21dee
movsd 0x17552(%rip), %xmm7    # 0x39338
movq 0x160(%rsp), %rcx
xorps %xmm5, %xmm5
cvtsi2sdl 0xd0(%rsp), %xmm5
xorps %xmm4, %xmm4
cvtsi2sdl 0xd8(%rsp), %xmm4
xorps %xmm3, %xmm3
cvtsi2sdl 0xe0(%rsp), %xmm3
movapd %xmm0, %xmm8
subsd 0x70(%rax), %xmm8
mulsd 0x18(%rax), %xmm5
movapd %xmm1, %xmm6
subsd 0x80(%rax), %xmm6
mulsd 0x20(%rax), %xmm4
subsd %xmm5, %xmm8
movapd %xmm2, %xmm5
subsd 0x90(%rax), %xmm5
subsd %xmm4, %xmm6
mulsd 0x28(%rax), %xmm3
subsd %xmm3, %xmm5
movq 0x20(%rsp), %r11
movsd 0x20(%r11), %xmm4
mulsd %xmm8, %xmm4
movsd 0x174c7(%rip), %xmm3    # 0x39330
mulsd %xmm3, %xmm4
cvttsd2si %xmm4, %esi
movsd 0x28(%r11), %xmm4
movsd %xmm6, 0xc0(%rsp)
mulsd %xmm6, %xmm4
mulsd %xmm3, %xmm4
cvttsd2si %xmm4, %edi
movsd 0x30(%r11), %xmm4
movsd %xmm5, 0xb8(%rsp)
mulsd %xmm5, %xmm4
mulsd %xmm3, %xmm4
xorl %ebx, %ebx
cmpl $0x4, %esi
jl 0x21ed4
movsd 0x8(%r11), %xmm6
subsd %xmm8, %xmm6
movl $0x7, %r8d
subl %esi, %r8d
xorl %esi, %esi
testl %r8d, %r8d
cmovgl %r8d, %esi
movl $0x200001, %r14d        # imm = 0x200001
movl $0x60007f, %ebx         # imm = 0x60007F
jmp 0x21edc
movapd %xmm8, %xmm6
xorl %r14d, %r14d
cvttsd2si %xmm4, %r8d
cmpl $0x4, %edi
jl 0x21f16
movsd 0x10(%r11), %xmm5
subsd 0xc0(%rsp), %xmm5
orl $0x3003f80, %ebx        # imm = 0x3003F80
orl $0x1000080, %r14d       # imm = 0x1000080
movl $0x7, %r10d
subl %edi, %r10d
xorl %edi, %edi
testl %r10d, %r10d
cmovgl %r10d, %edi
jmp 0x21f1f
movsd 0xc0(%rsp), %xmm5
movl 0x10(%rsp), %r10d
cmpl $0x4, %r8d
movq %r9, 0x68(%rsp)
jl 0x21f60
movsd 0x18(%r11), %xmm3
subsd 0xb8(%rsp), %xmm3
orl $0x181fc000, %ebx       # imm = 0x181FC000
orl $0x8004000, %r14d       # imm = 0x8004000
movl $0x7, %r9d
subl %r8d, %r9d
xorl %r8d, %r8d
testl %r9d, %r9d
cmovgl %r9d, %r8d
jmp 0x21f69
movsd 0xb8(%rsp), %xmm3
movl %ebx, %r9d
movsd 0xe0(%rax), %xmm4
mulsd %xmm4, %xmm4
addsd %xmm7, %xmm4
mulsd %xmm6, %xmm6
ucomisd %xmm4, %xmm6
jbe 0x21fac
mulsd %xmm5, %xmm5
ucomisd %xmm4, %xmm5
jbe 0x21fac
mulsd %xmm3, %xmm3
ucomisd %xmm4, %xmm3
jbe 0x21fac
addq $0x128, %rsp            # imm = 0x128
popq %rbx
popq %r12
popq %r13
popq %r14
popq %r15
popq %rbp
retq
imull 0x3c(%rsp), %ebp
addl %edx, %ebp
imull 0x30(%rsp), %ebp
addl %r10d, %ebp
subl %ebp, 0x38(%rsp)
movq 0x78(%r11), %rbx
movq 0x80(%r11), %rax
shll $0xa, %r8d
shll $0x8, %edi
shll $0x6, %esi
addl %edi, %esi
addl %r8d, %esi
movslq %esi, %rdi
leaq (%rax,%rdi,8), %rsi
movq %rbx, 0xb0(%rsp)
movq %rdi, 0xc8(%rsp)
leaq (%rbx,%rdi,4), %r8
movl (%r8), %eax
cmpl $0x2, %eax
movl $0x1, %edi
cmovgel %eax, %edi
xorpd %xmm3, %xmm3
movapd %xmm3, 0x90(%rsp)
xorl %eax, %eax
movapd %xmm3, 0x80(%rsp)
movapd %xmm3, 0x70(%rsp)
movsd %xmm2, 0x100(%rsp)
movsd %xmm1, 0xf8(%rsp)
movsd %xmm0, 0xf0(%rsp)
movl %edx, 0x14(%rsp)
movsd %xmm8, 0x60(%rsp)
movl %r9d, 0x2c(%rsp)
movl %r14d, 0x30(%rsp)
movq %rsi, 0x50(%rsp)
movq %rdi, 0x58(%rsp)
movq %r8, 0x18(%rsp)
movq %rax, %r13
movq (%r11), %rax
movsd (%rcx), %xmm3
movsd 0xe0(%rax), %xmm0
mulsd %xmm0, %xmm0
addsd %xmm3, %xmm0
movsd (%rsi,%r13,8), %xmm1
ucomisd %xmm0, %xmm1
ja 0x21f9a
movl 0x4(%r8,%r13,4), %ebx
xorl %r9d, %ebx
addl %r14d, %ebx
movl %ebx, %ebp
andl $0x7f, %ebp
addl $-0x40, %ebp
movl %ebp, %r15d
addl %r10d, %r15d
js 0x2240d
cmpl 0x38(%r11), %r15d
jge 0x2240d
movl %ebx, %r14d
shrl $0x7, %r14d
andl $0x7f, %r14d
addl $-0x40, %r14d
movl %r14d, %r12d
addl %edx, %r12d
js 0x22137
cmpl 0x3c(%r11), %r12d
jge 0x22137
shrl $0xe, %ebx
andl $0x7f, %ebx
addl $-0x40, %ebx
movl %ebx, %eax
addl 0x3c(%rsp), %eax
js 0x22137
cmpl 0x40(%r11), %eax
jge 0x22137
movq %rax, 0x40(%rsp)
movq %r11, %rdi
movl %ebp, %esi
movl %r14d, %edx
movl %ebx, %ecx
movapd %xmm8, %xmm0
movsd 0xc0(%rsp), %xmm1
movsd 0xb8(%rsp), %xmm2
callq 0x2304e
testb %al, %al
je 0x22141
movq 0x20(%rsp), %r11
movq 0x160(%rsp), %rcx
movl 0x10(%rsp), %r10d
movl 0x14(%rsp), %edx
movsd 0x60(%rsp), %xmm8
jmp 0x223f4
movl 0x30(%rsp), %r14d
jmp 0x2240d
movq 0x20(%rsp), %r11
movq (%r11), %rax
cmpb $0x1, 0xa8(%rax)
movsd 0x100(%rsp), %xmm2
movsd 0xf8(%rsp), %xmm1
movsd 0xf0(%rsp), %xmm0
movq 0x18(%rsp), %r8
jne 0x221ab
movq 0xd0(%rsp), %rcx
leal (%r15,%rcx), %edx
movl (%rax), %ecx
cmpl %ecx, %edx
movsd 0x60(%rsp), %xmm8
jge 0x221fb
addl %r15d, %ecx
movsd 0x78(%rax), %xmm3
subsd 0x70(%rax), %xmm3
xorpd 0x14ea0(%rip), %xmm3    # 0x37040
movapd %xmm3, 0x70(%rsp)
movl %ecx, %r15d
jmp 0x221b2
movsd 0x60(%rsp), %xmm8
cmpb $0x1, 0xa9(%rax)
jne 0x22251
movq 0xd8(%rsp), %rcx
leal (%r12,%rcx), %edx
movl 0x4(%rax), %ecx
cmpl %ecx, %edx
jge 0x22221
addl %r12d, %ecx
movsd 0x88(%rax), %xmm3
subsd 0x80(%rax), %xmm3
xorpd 0x14e53(%rip), %xmm3    # 0x37040
movapd %xmm3, 0x80(%rsp)
movl %ecx, %r12d
jmp 0x22251
leal (%rcx,%rcx), %esi
xorpd %xmm3, %xmm3
movapd %xmm3, 0x70(%rsp)
cmpl %esi, %edx
jl 0x221b2
subl %ecx, %r15d
movsd 0x78(%rax), %xmm3
subsd 0x70(%rax), %xmm3
movapd %xmm3, 0x70(%rsp)
jmp 0x221b2
leal (%rcx,%rcx), %esi
xorpd %xmm3, %xmm3
movapd %xmm3, 0x80(%rsp)
cmpl %esi, %edx
jl 0x22251
subl %ecx, %r12d
movsd 0x88(%rax), %xmm3
subsd 0x80(%rax), %xmm3
movapd %xmm3, 0x80(%rsp)
cmpb $0x1, 0xaa(%rax)
jne 0x222b0
movq 0xe0(%rsp), %rcx
movq 0x40(%rsp), %rsi
leal (%rsi,%rcx), %edx
movl 0x8(%rax), %ecx
cmpl %ecx, %edx
jge 0x22298
addl %esi, %ecx
movsd 0x98(%rax), %xmm3
subsd 0x90(%rax), %xmm3
xorpd 0x14db5(%rip), %xmm3    # 0x37040
movapd %xmm3, 0x90(%rsp)
movl %ecx, %edx
jmp 0x222b5
leal (%rcx,%rcx), %esi
xorpd %xmm3, %xmm3
movapd %xmm3, 0x90(%rsp)
cmpl %esi, %edx
jge 0x22443
movq 0x40(%rsp), %rdx
addl 0x38(%rsp), %r15d
imull 0x4(%rax), %edx
addl %r12d, %edx
imull (%rax), %edx
addl %r15d, %edx
movq 0x60(%r11), %rcx
movq %rdx, %r12
movslq %edx, %rdx
movl (%rcx,%rdx,4), %r10d
testl %r10d, %r10d
jle 0x2241c
subsd 0x70(%rsp), %xmm0
subsd 0x80(%rsp), %xmm1
subsd 0x90(%rsp), %xmm2
movq 0x160(%rsp), %rsi
movsd (%rsi), %xmm3
movq 0x58(%r11), %rsi
movq (%rsi,%rdx,8), %rsi
movslq 0x4c(%r11), %rdi
movq 0xd8(%rax), %rax
movq (%rax,%rdx,8), %rax
addq $0x18, %rax
addq $0x10, %rsi
shlq $0x3, %rdi
xorl %r8d, %r8d
xorl %r9d, %r9d
movslq %r10d, %r10
movq %r8, %r11
shlq $0x5, %r11
addq %rax, %r11
movq %rdi, %r15
imulq %r8, %r15
addq %rsi, %r15
movsd -0x10(%r15), %xmm5
subsd %xmm0, %xmm5
movsd -0x8(%r15), %xmm6
subsd %xmm1, %xmm6
movsd (%r15), %xmm4
subsd %xmm2, %xmm4
mulsd %xmm6, %xmm6
mulsd %xmm5, %xmm5
addsd %xmm6, %xmm5
mulsd %xmm4, %xmm4
addsd %xmm5, %xmm4
movsd (%r11), %xmm5
mulsd %xmm5, %xmm5
subsd %xmm5, %xmm4
ucomisd %xmm4, %xmm3
ja 0x22397
incq %r8
addq $0x20, %r11
addq %rdi, %r15
cmpq %r10, %r8
jl 0x22342
jmp 0x223c6
movq 0x160(%rsp), %r9
movsd %xmm4, (%r9)
movq 0x68(%rsp), %r9
movl %r8d, 0x4(%r9)
movslq (%rcx,%rdx,4), %r10
incq %r8
movb $0x1, %r9b
movapd %xmm4, %xmm3
cmpq %r10, %r8
jl 0x2232b
jmp 0x223cc
testb $0x1, %r9b
je 0x223de
movq 0x68(%rsp), %rax
movl %r12d, (%rax)
movl %ebp, 0x8(%rax)
movl %r14d, 0xc(%rax)
movl %ebx, 0x10(%rax)
movq 0x20(%rsp), %r11
movq 0x160(%rsp), %rcx
movl 0x10(%rsp), %r10d
movl 0x14(%rsp), %edx
movl 0x2c(%rsp), %r9d
movl 0x30(%rsp), %r14d
movq 0x50(%rsp), %rsi
movq 0x58(%rsp), %rdi
movq 0x18(%rsp), %r8
leaq 0x1(%r13), %rax
cmpq %rdi, %rax
jne 0x2205f
jmp 0x22468
movq 0x160(%rsp), %rcx
movl 0x10(%rsp), %r10d
movl 0x14(%rsp), %edx
movl 0x2c(%rsp), %r9d
movl 0x30(%rsp), %r14d
movq 0x50(%rsp), %rsi
movq 0x58(%rsp), %rdi
jmp 0x2240d
movq 0x40(%rsp), %rdx
subl %ecx, %edx
movsd 0x98(%rax), %xmm3
subsd 0x90(%rax), %xmm3
movapd %xmm3, 0x90(%rsp)
jmp 0x222b5
incl 0x70(%r11)
jne 0x224aa
cmpl $0x0, 0x48(%r11)
jle 0x22496
movq 0x88(%r11), %rax
movq %rax, %rcx
movl $0x0, (%rcx)
addq $0x4, %rcx
movslq 0x48(%r11), %rdx
leaq (%rax,%rdx,4), %rdx
cmpq %rdx, %rcx
jb 0x2247f
movl $0x1, 0x70(%r11)
movq 0x160(%rsp), %rcx
movl 0x14(%rsp), %edx
movq 0x90(%r11), %rax
movq %rax, 0xe8(%rsp)
movq %rax, 0x48(%rsp)
cmpq $0x3d, %r13
ja 0x22b28
movq 0xb0(%rsp), %rax
movq 0xc8(%rsp), %r8
leaq (%rax,%r8,4), %rax
addq $0x4, %rax
movq %rax, 0x120(%rsp)
movq %rdi, 0x108(%rsp)
movslq %edi, %rdi
leaq (%rsi,%rdi,8), %r8
movq %r8, 0x18(%rsp)
leaq (%rax,%rdi,4), %rax
movq %rax, 0x40(%rsp)
xorl %r13d, %r13d
movq (%r11), %rax
movsd (%rcx), %xmm3
movsd 0xe0(%rax), %xmm0
mulsd %xmm0, %xmm0
addsd %xmm3, %xmm0
movq 0x18(%rsp), %rax
movsd (%rax,%r13,8), %xmm1
ucomisd %xmm0, %xmm1
ja 0x21f9a
movq 0x40(%rsp), %rax
movl (%rax,%r13,4), %ebp
xorl %r9d, %ebp
addl %r14d, %ebp
movl %ebp, %r14d
andl $0x7f, %r14d
addl $-0x40, %r14d
movl %r14d, %r8d
addl %r10d, %r8d
js 0x22681
movl 0x38(%r11), %eax
cmpl %eax, %r8d
jge 0x22681
movl %ebp, %ebx
shrl $0x7, %ebx
andl $0x7f, %ebx
addl $-0x40, %ebx
movl %ebx, 0x58(%rsp)
addl %edx, %ebx
js 0x22681
movl 0x3c(%r11), %r15d
movq %r15, 0xb0(%rsp)
cmpl %r15d, %ebx
jge 0x22681
movl %ebp, %r15d
shrl $0xe, %r15d
andl $0x7f, %r15d
addl $-0x40, %r15d
movl %r15d, 0xc8(%rsp)
addl 0x3c(%rsp), %r15d
js 0x22681
cmpl 0x40(%r11), %r15d
jge 0x22681
movl %r8d, %r12d
shlq $0x2, %r12
addq 0x88(%r11), %r12
movq 0xb0(%rsp), %rcx
imull %r15d, %ecx
addl %ebx, %ecx
imull %eax, %ecx
movl 0x70(%r11), %eax
movq %rcx, 0xb0(%rsp)
movl %eax, (%r12,%rcx,4)
movq %rdi, 0x110(%rsp)
movq %r11, %rdi
movl %r14d, %esi
movl 0x58(%rsp), %edx
movl 0xc8(%rsp), %ecx
movapd %xmm8, %xmm0
movsd 0xc0(%rsp), %xmm1
movsd 0xb8(%rsp), %xmm2
movq %r8, 0xa8(%rsp)
callq 0x2304e
movq 0xa8(%rsp), %r8
movq 0x110(%rsp), %rdi
movq 0x50(%rsp), %rsi
movl 0x2c(%rsp), %r9d
movsd 0x60(%rsp), %xmm8
movl 0x14(%rsp), %edx
movl 0x10(%rsp), %r10d
movq 0x160(%rsp), %rcx
movq 0x20(%rsp), %r11
movsd 0xf0(%rsp), %xmm0
movsd 0xf8(%rsp), %xmm1
movsd 0x100(%rsp), %xmm2
testb %al, %al
je 0x226a2
leaq (%rdi,%r13), %rax
incq %rax
incq %r13
decq %rax
cmpq $0x3e, %rax
movl 0x30(%rsp), %r14d
jl 0x22508
jmp 0x22b1b
movq %r12, 0x18(%rsp)
movq (%r11), %rax
cmpb $0x1, 0xa8(%rax)
movl %r8d, %edx
movq %rbx, %r10
jne 0x22713
movq 0xd0(%rsp), %rcx
leal (%r8,%rcx), %esi
movl (%rax), %ecx
cmpl %ecx, %esi
jge 0x226ea
addl %r8d, %ecx
movsd 0x78(%rax), %xmm3
subsd 0x70(%rax), %xmm3
xorpd 0x14960(%rip), %xmm3    # 0x37040
movapd %xmm3, 0x70(%rsp)
movl %ecx, %edx
jmp 0x22713
leal (%rcx,%rcx), %edi
xorpd %xmm3, %xmm3
movapd %xmm3, 0x70(%rsp)
movl %r8d, %edx
cmpl %edi, %esi
jl 0x22713
movl %r8d, %edx
subl %ecx, %edx
movsd 0x78(%rax), %xmm3
subsd 0x70(%rax), %xmm3
movapd %xmm3, 0x70(%rsp)
cmpb $0x1, 0xa9(%rax)
movl %r10d, %esi
jne 0x22791
movq 0xd8(%rsp), %rcx
leal (%r10,%rcx), %edi
movl 0x4(%rax), %ecx
cmpl %ecx, %edi
jge 0x2275a
addl %r10d, %ecx
movsd 0x88(%rax), %xmm3
subsd 0x80(%rax), %xmm3
xorpd 0x148f3(%rip), %xmm3    # 0x37040
movapd %xmm3, 0x80(%rsp)
movl %ecx, %esi
jmp 0x22791
leal (%rcx,%rcx), %r8d
xorpd %xmm3, %xmm3
movapd %xmm3, 0x80(%rsp)
movl %r10d, %esi
cmpl %r8d, %edi
jl 0x22791
movl %r10d, %esi
subl %ecx, %esi
movsd 0x88(%rax), %xmm3
subsd 0x80(%rax), %xmm3
movapd %xmm3, 0x80(%rsp)
cmpb $0x1, 0xaa(%rax)
movl %r15d, %ecx
jne 0x22810
movq 0xe0(%rsp), %rcx
leal (%r15,%rcx), %r8d
movl 0x8(%rax), %edi
cmpl %edi, %r8d
jge 0x227d9
addl %r15d, %edi
movsd 0x98(%rax), %xmm3
subsd 0x90(%rax), %xmm3
xorpd 0x14874(%rip), %xmm3    # 0x37040
movapd %xmm3, 0x90(%rsp)
movl %edi, %ecx
jmp 0x22810
leal (%rdi,%rdi), %r9d
xorpd %xmm3, %xmm3
movapd %xmm3, 0x90(%rsp)
movl %r15d, %ecx
cmpl %r9d, %r8d
jl 0x22810
movl %r15d, %ecx
subl %edi, %ecx
movsd 0x98(%rax), %xmm3
subsd 0x90(%rax), %xmm3
movapd %xmm3, 0x90(%rsp)
movq %r15, 0x118(%rsp)
movq %r10, 0x40(%rsp)
addl 0x38(%rsp), %edx
imull 0x4(%rax), %ecx
addl %esi, %ecx
imull (%rax), %ecx
addl %edx, %ecx
movq 0x60(%r11), %rdx
movslq %ecx, %rsi
movl (%rdx,%rsi,4), %r11d
testl %r11d, %r11d
jle 0x22951
subsd 0x70(%rsp), %xmm0
subsd 0x80(%rsp), %xmm1
subsd 0x90(%rsp), %xmm2
movq 0x160(%rsp), %rdi
movsd (%rdi), %xmm3
movq 0x20(%rsp), %r8
movq 0x58(%r8), %rdi
movq (%rdi,%rsi,8), %rdi
movslq 0x4c(%r8), %r8
movq 0xd8(%rax), %rax
movq (%rax,%rsi,8), %rax
addq $0x18, %rax
addq $0x10, %rdi
shlq $0x3, %r8
xorl %r9d, %r9d
xorl %r10d, %r10d
movslq %r11d, %r11
movq %r9, %rbx
shlq $0x5, %rbx
addq %rax, %rbx
movq %r8, %r12
imulq %r9, %r12
addq %rdi, %r12
movsd -0x10(%r12), %xmm5
subsd %xmm0, %xmm5
movsd -0x8(%r12), %xmm6
subsd %xmm1, %xmm6
movsd (%r12), %xmm4
subsd %xmm2, %xmm4
mulsd %xmm6, %xmm6
mulsd %xmm5, %xmm5
addsd %xmm6, %xmm5
mulsd %xmm4, %xmm4
addsd %xmm5, %xmm4
movsd (%rbx), %xmm5
mulsd %xmm5, %xmm5
subsd %xmm5, %xmm4
ucomisd %xmm4, %xmm3
ja 0x22900
incq %r9
addq $0x20, %rbx
addq %r8, %r12
cmpq %r11, %r9
jl 0x228a9
jmp 0x2292f
movq 0x160(%rsp), %r10
movsd %xmm4, (%r10)
movq 0x68(%rsp), %r10
movl %r9d, 0x4(%r10)
movslq (%rdx,%rsi,4), %r11
incq %r9
movb $0x1, %r10b
movapd %xmm4, %xmm3
cmpq %r11, %r9
jl 0x22892
jmp 0x22935
testb $0x1, %r10b
je 0x22951
movq 0x68(%rsp), %rax
movl %ecx, (%rax)
movl %r14d, 0x8(%rax)
movl 0x58(%rsp), %ecx
movl %ecx, 0xc(%rax)
movl 0xc8(%rsp), %ecx
movl %ecx, 0x10(%rax)
movq 0x20(%rsp), %rsi
movq 0x98(%rsi), %rax
addq $-0x48, %rax
cmpq %rax, 0x48(%rsp)
ja 0x229ff
movq 0xa8(%rsp), %rcx
movq 0x18(%rsp), %r15
movq 0xb0(%rsp), %rax
leaq (%r15,%rax,4), %rdx
movq 0x110(%rsp), %rbx
movl %ebx, %eax
leaq (%rax,%r13), %r14
incq %r14
addq %r13, %rbx
leaq 0x48(%rsp), %rax
movq %rax, (%rsp)
movq %rsi, %rdi
movl %ebp, %esi
movq 0x40(%rsp), %r8
movq 0x118(%rsp), %r9
callq 0x1ec3a
movq %r14, %rdi
movq 0x20(%rsp), %r11
cmpq $0x3e, %rbx
movq 0x160(%rsp), %rcx
movl 0x10(%rsp), %r10d
movl 0x14(%rsp), %edx
movsd 0x60(%rsp), %xmm8
movl 0x2c(%rsp), %r9d
movl 0x30(%rsp), %r14d
movq 0x50(%rsp), %rsi
movq 0x120(%rsp), %rax
jl 0x224e8
jmp 0x2303b
movl 0x74(%rsi), %eax
addl %eax, %eax
movl %eax, 0x74(%rsi)
movslq %eax, %r14
leaq (,%r14,4), %rdi
testl %r14d, %r14d
movq $-0x1, %rax
cmovsq %rax, %rdi
callq 0x2030
movq %rax, %rbx
movq 0x225b1(%rip), %rax     # 0x44fe0
movq (%rax), %rdi
leaq 0x169f0(%rip), %rsi     # 0x39429
movl %r14d, %edx
xorl %eax, %eax
callq 0x21b0
movq 0x48(%rsp), %rax
movq 0xe8(%rsp), %rdi
cmpq %rax, %rdi
jbe 0x22a78
movq 0x20(%rsp), %rsi
movq 0x98(%rsi), %rcx
movq %rbx, %r14
cmpq %rcx, %rdi
jae 0x22a9b
movl (%rdi), %edx
addq $0x4, %rdi
movl %edx, (%r14)
addq $0x4, %r14
jmp 0x22a64
movq %rbx, %r14
movq 0x20(%rsp), %rsi
movq 0x18(%rsp), %r15
jae 0x22abe
movl (%rdi), %ecx
addq $0x4, %rdi
movl %ecx, (%r14)
addq $0x4, %r14
cmpq %rax, %rdi
jb 0x22a87
jmp 0x22abe
movq 0x90(%rsi), %rdi
cmpq %rax, %rdi
jae 0x22acf
movq 0x18(%rsp), %r15
movl (%rdi), %ecx
addq $0x4, %rdi
movl %ecx, (%r14)
addq $0x4, %r14
cmpq %rax, %rdi
jb 0x22aac
movq 0x90(%rsi), %rdi
movq 0xa8(%rsp), %rcx
jmp 0x22adc
movq 0xa8(%rsp), %rcx
movq 0x18(%rsp), %r15
testq %rdi, %rdi
je 0x22af3
callq 0x2180
movq 0xa8(%rsp), %rcx
movq 0x20(%rsp), %rsi
movq %rbx, 0x90(%rsi)
movslq 0x74(%rsi), %rax
leaq (%rbx,%rax,4), %rax
movq %rax, 0x98(%rsi)
movq %r14, 0x48(%rsp)
movq %rbx, 0xe8(%rsp)
jmp 0x22979
movq 0x108(%rsp), %rax
addl %r13d, %eax
movl %eax, %edi
movq (%r11), %rax
movsd 0xe0(%rax), %xmm0
mulsd %xmm0, %xmm0
addsd (%rcx), %xmm0
movl %edi, %eax
movsd (%rsi,%rax,8), %xmm1
ucomisd %xmm0, %xmm1
ja 0x21f9a
movq 0xe8(%rsp), %r15
cmpq 0x48(%rsp), %r15
je 0x21f9a
cmpq 0x98(%r11), %r15
jne 0x22b6f
movq 0x90(%r11), %r15
movl (%r15), %r12d
movl 0x4(%r15), %eax
movl 0x8(%r15), %r14d
addq $0xc, %r15
movl %r12d, %ebp
subl %r10d, %ebp
movq %rax, 0x30(%rsp)
movl %eax, %r13d
subl %edx, %r13d
movl %r14d, %ebx
subl 0x3c(%rsp), %ebx
movsd (%rcx), %xmm3
movq %r11, %rdi
movl %ebp, %esi
movl %r13d, %edx
movl %ebx, %ecx
movapd %xmm8, %xmm0
movsd 0xc0(%rsp), %xmm1
movsd 0xb8(%rsp), %xmm2
callq 0x2304e
movsd 0x60(%rsp), %xmm8
movl 0x14(%rsp), %edx
movl 0x10(%rsp), %r10d
movq 0x160(%rsp), %rcx
movq 0x20(%rsp), %r11
testb %al, %al
jne 0x22b54
movq (%r11), %rax
cmpb $0x1, 0xa8(%rax)
movl %r12d, %edx
movsd 0x100(%rsp), %xmm2
movsd 0xf8(%rsp), %xmm1
movsd 0xf0(%rsp), %xmm0
movq %r12, %r10
movq 0x30(%rsp), %r12
movl %ebp, 0x50(%rsp)
movl %r13d, 0x2c(%rsp)
movl %ebx, 0x40(%rsp)
jne 0x22c7e
movq 0xd0(%rsp), %rcx
leal (%r10,%rcx), %esi
movl (%rax), %ecx
cmpl %ecx, %esi
jge 0x22c55
addl %r10d, %ecx
movsd 0x78(%rax), %xmm3
subsd 0x70(%rax), %xmm3
xorpd 0x143f5(%rip), %xmm3    # 0x37040
movapd %xmm3, 0x70(%rsp)
movl %ecx, %edx
jmp 0x22c7e
leal (%rcx,%rcx), %edi
xorpd %xmm3, %xmm3
movapd %xmm3, 0x70(%rsp)
movl %r10d, %edx
cmpl %edi, %esi
jl 0x22c7e
movl %r10d, %edx
subl %ecx, %edx
movsd 0x78(%rax), %xmm3
subsd 0x70(%rax), %xmm3
movapd %xmm3, 0x70(%rsp)
cmpb $0x1, 0xa9(%rax)
movl %r12d, %esi
jne 0x22cfc
movq 0xd8(%rsp), %rcx
leal (%r12,%rcx), %edi
movl 0x4(%rax), %ecx
cmpl %ecx, %edi
jge 0x22cc5
addl %r12d, %ecx
movsd 0x88(%rax), %xmm3
subsd 0x80(%rax), %xmm3
xorpd 0x14388(%rip), %xmm3    # 0x37040
movapd %xmm3, 0x80(%rsp)
movl %ecx, %esi
jmp 0x22cfc
leal (%rcx,%rcx), %r8d
xorpd %xmm3, %xmm3
movapd %xmm3, 0x80(%rsp)
movl %r12d, %esi
cmpl %r8d, %edi
jl 0x22cfc
movl %r12d, %esi
subl %ecx, %esi
movsd 0x88(%rax), %xmm3
subsd 0x80(%rax), %xmm3
movapd %xmm3, 0x80(%rsp)
cmpb $0x1, 0xaa(%rax)
movl %r14d, %ecx
jne 0x22d7b
movq 0xe0(%rsp), %rcx
leal (%r14,%rcx), %r8d
movl 0x8(%rax), %edi
cmpl %edi, %r8d
jge 0x22d44
addl %r14d, %edi
movsd 0x98(%rax), %xmm3
subsd 0x90(%rax), %xmm3
xorpd 0x14309(%rip), %xmm3    # 0x37040
movapd %xmm3, 0x90(%rsp)
movl %edi, %ecx
jmp 0x22d7b
leal (%rdi,%rdi), %r9d
xorpd %xmm3, %xmm3
movapd %xmm3, 0x90(%rsp)
movl %r14d, %ecx
cmpl %r9d, %r8d
jl 0x22d7b
movl %r14d, %ecx
subl %edi, %ecx
movsd 0x98(%rax), %xmm3
subsd 0x90(%rax), %xmm3
movapd %xmm3, 0x90(%rsp)
movq %r12, %rbp
movq %r14, %r12
movq %r10, 0x18(%rsp)
addl 0x38(%rsp), %edx
imull 0x4(%rax), %ecx
addl %esi, %ecx
imull (%rax), %ecx
addl %edx, %ecx
movq 0x60(%r11), %rdx
movslq %ecx, %rsi
movl (%rdx,%rsi,4), %r8d
testl %r8d, %r8d
jle 0x22ec5
movq %rdx, %r14
subsd 0x70(%rsp), %xmm0
subsd 0x80(%rsp), %xmm1
subsd 0x90(%rsp), %xmm2
movq 0x160(%rsp), %rdx
movsd (%rdx), %xmm3
movq 0x58(%r11), %rdx
movq (%rdx,%rsi,8), %rbx
movslq 0x4c(%r11), %rdi
movq 0xd8(%rax), %rax
movq (%rax,%rsi,8), %rax
addq $0x18, %rax
addq $0x10, %rbx
shlq $0x3, %rdi
xorl %edx, %edx
xorl %r10d, %r10d
movslq %r8d, %r11
movq %rdx, %r9
shlq $0x5, %r9
addq %rax, %r9
movq %rdi, %r8
imulq %rdx, %r8
movq %rbx, %r13
addq %rbx, %r8
movsd -0x10(%r8), %xmm5
subsd %xmm0, %xmm5
movsd -0x8(%r8), %xmm6
subsd %xmm1, %xmm6
movsd (%r8), %xmm4
subsd %xmm2, %xmm4
mulsd %xmm6, %xmm6
mulsd %xmm5, %xmm5
addsd %xmm6, %xmm5
mulsd %xmm4, %xmm4
addsd %xmm5, %xmm4
movsd (%r9), %xmm5
mulsd %xmm5, %xmm5
subsd %xmm5, %xmm4
ucomisd %xmm4, %xmm3
ja 0x22e67
incq %rdx
addq $0x20, %r9
addq %rdi, %r8
cmpq %r11, %rdx
jl 0x22e12
jmp 0x22e9e
movq 0x160(%rsp), %r8
movsd %xmm4, (%r8)
movq 0x68(%rsp), %r8
movl %edx, 0x4(%r8)
movslq (%r14,%rsi,4), %r8
incq %rdx
movb $0x1, %r10b
movapd %xmm4, %xmm3
cmpq %r8, %rdx
movq 0x20(%rsp), %r11
movq %r13, %rbx
jl 0x22df8
jmp 0x22ea9
testb $0x1, %r10b
movq 0x20(%rsp), %r11
je 0x22ec5
movq 0x68(%rsp), %rax
movl %ecx, (%rax)
movl 0x50(%rsp), %ecx
movl %ecx, 0x8(%rax)
movl 0x2c(%rsp), %ecx
movl %ecx, 0xc(%rax)
movl 0x40(%rsp), %ecx
movl %ecx, 0x10(%rax)
movq 0x48(%rsp), %rcx
movq %r15, %rax
subq %rcx, %rax
jbe 0x22ed8
sarq $0x2, %rax
jmp 0x22ef7
movq 0x98(%r11), %rdx
subq %rcx, %rdx
sarq $0x2, %rdx
movq %r15, %rax
subq 0x90(%r11), %rax
sarq $0x2, %rax
addq %rdx, %rax
movq 0x18(%rsp), %rsi
movq %r12, %r14
cmpq $0x11, %rax
movq %rbp, %r12
jle 0x22f3d
movq %r11, %rdi
movl %r12d, %edx
movl %r14d, %ecx
leaq 0x48(%rsp), %r8
callq 0x1f78e
movq 0x20(%rsp), %r11
movq 0x160(%rsp), %rcx
movl 0x10(%rsp), %r10d
movl 0x14(%rsp), %edx
movsd 0x60(%rsp), %xmm8
jmp 0x22b54
movl 0x74(%r11), %eax
addl %eax, %eax
movl %eax, 0x74(%r11)
movslq %eax, %r13
leaq (,%r13,4), %rdi
testl %r13d, %r13d
movq $-0x1, %rax
cmovsq %rax, %rdi
movq %r11, %rbp
callq 0x2030
movq %rax, %rbx
movq 0x2206e(%rip), %rax     # 0x44fe0
movq (%rax), %rdi
leaq 0x164ad(%rip), %rsi     # 0x39429
movl %r13d, %edx
xorl %eax, %eax
callq 0x21b0
movq 0x48(%rsp), %rax
cmpq %rax, %r15
jbe 0x22fdc
movq 0x98(%rbp), %rcx
movq %rbx, %r13
cmpq %rcx, %r15
movq %rbp, %r11
jae 0x22fb6
movl (%r15), %edx
addq $0x4, %r15
movl %edx, (%r13)
addq $0x4, %r13
cmpq %rcx, %r15
jb 0x22fa2
movq 0x90(%r11), %rdi
cmpq %rax, %rdi
movq 0x18(%rsp), %rsi
jae 0x23004
movl (%rdi), %ecx
addq $0x4, %rdi
movl %ecx, (%r13)
addq $0x4, %r13
cmpq %rax, %rdi
jb 0x22fc7
jmp 0x22ffd
movq %rbx, %r13
movq %rbp, %r11
movq 0x18(%rsp), %rsi
jae 0x22ffd
movl (%r15), %ecx
addq $0x4, %r15
movl %ecx, (%r13)
addq $0x4, %r13
cmpq %rax, %r15
jb 0x22fe9
movq 0x90(%r11), %rdi
testq %rdi, %rdi
je 0x23018
callq 0x2180
movq 0x18(%rsp), %rsi
movq 0x20(%rsp), %r11
movq %rbx, 0x90(%r11)
movslq 0x74(%r11), %rax
leaq (%rbx,%rax,4), %rax
movq %rax, 0x98(%r11)
movq %r13, 0x48(%rsp)
movq %rbx, %r15
jmp 0x22f08
movq 0x108(%rsp), %rax
leal (%rax,%r13), %edi
incl %edi
jmp 0x22b28