// File: Luthaf[P]zeopp/voro++/src/v_compute.cc
// Function: voro::voro_compute<voro::container>::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
movq %rsp, %rbp
pushq %r14
pushq %rbx
subq $0x130, %rsp            # imm = 0x130
movq 0x10(%rbp), %rax
movq %rdi, -0x18(%rbp)
movsd %xmm0, -0x20(%rbp)
movsd %xmm1, -0x28(%rbp)
movsd %xmm2, -0x30(%rbp)
movl %esi, -0x34(%rbp)
movl %edx, -0x38(%rbp)
movl %ecx, -0x3c(%rbp)
movl %r8d, -0x40(%rbp)
movq %r9, -0x48(%rbp)
movq -0x18(%rbp), %rax
movq %rax, -0x108(%rbp)
xorps %xmm0, %xmm0
movsd %xmm0, -0x50(%rbp)
xorps %xmm0, %xmm0
movsd %xmm0, -0x58(%rbp)
xorps %xmm0, %xmm0
movsd %xmm0, -0x60(%rbp)
movq -0x48(%rbp), %rcx
movl $0xffffffff, (%rcx)     # imm = 0xFFFFFFFF
movq 0x10(%rbp), %rcx
movsd 0x2ec21(%rip), %xmm0    # 0x58000
movsd %xmm0, (%rcx)
movq (%rax), %rdi
movl -0x34(%rbp), %esi
movl -0x38(%rbp), %edx
movl -0x3c(%rbp), %ecx
movl -0x40(%rbp), %r8d
leaq -0x6c(%rbp), %r9
leaq -0x70(%rbp), %r11
leaq -0x74(%rbp), %r10
leaq -0x98(%rbp), %rax
movq %r11, (%rsp)
movq %r10, 0x8(%rsp)
movq %rax, 0x10(%rsp)
callq 0x29f40
movq -0x108(%rbp), %rdi
movl -0x40(%rbp), %esi
movsd -0x20(%rbp), %xmm0
movsd -0x28(%rbp), %xmm1
movsd -0x30(%rbp), %xmm2
movq -0x48(%rbp), %r9
movq 0x10(%rbp), %rax
xorl %r8d, %r8d
movl %r8d, %edx
movl %r8d, %ecx
movq %rax, (%rsp)
callq 0x2a010
movq -0x108(%rbp), %rax
movq (%rax), %rdi
movsd -0x20(%rbp), %xmm0
movsd -0x28(%rbp), %xmm1
movsd -0x30(%rbp), %xmm2
cvtsi2sdl -0x34(%rbp), %xmm3
cvtsi2sdl -0x38(%rbp), %xmm4
cvtsi2sdl -0x3c(%rbp), %xmm5
leaq -0xa0(%rbp), %rsi
leaq -0xa8(%rbp), %rdx
leaq -0xb0(%rbp), %rcx
callq 0x21f10
movq -0x108(%rbp), %rax
movsd -0xa0(%rbp), %xmm0
mulsd 0x20(%rax), %xmm0
movsd 0x2eb4e(%rip), %xmm1    # 0x57ff8
mulsd %xmm1, %xmm0
cvttsd2si %xmm0, %ecx
movl %ecx, -0x78(%rbp)
movsd -0xa8(%rbp), %xmm0
mulsd 0x28(%rax), %xmm0
movsd 0x2eb2e(%rip), %xmm1    # 0x57ff8
mulsd %xmm1, %xmm0
cvttsd2si %xmm0, %ecx
movl %ecx, -0x7c(%rbp)
movsd -0xb0(%rbp), %xmm0
mulsd 0x30(%rax), %xmm0
movsd 0x2eb0e(%rip), %xmm1    # 0x57ff8
mulsd %xmm1, %xmm0
cvttsd2si %xmm0, %eax
movl %eax, -0x80(%rbp)
cmpl $0x4, -0x78(%rbp)
jl 0x29545
movq -0x108(%rbp), %rax
movsd 0x8(%rax), %xmm0
subsd -0xa0(%rbp), %xmm0
movsd %xmm0, -0xb8(%rbp)
movl $0x60007f, -0xec(%rbp)  # imm = 0x60007F
movl $0x200001, -0xf0(%rbp)  # imm = 0x200001
movl $0x7, %eax
subl -0x78(%rbp), %eax
movl %eax, -0x78(%rbp)
cmpl $0x0, -0x78(%rbp)
jge 0x29543
movl $0x0, -0x78(%rbp)
jmp 0x29569
movl $0x0, -0xf0(%rbp)
movl $0x0, -0xec(%rbp)
movsd -0xa0(%rbp), %xmm0
movsd %xmm0, -0xb8(%rbp)
cmpl $0x4, -0x7c(%rbp)
jl 0x295c7
movq -0x108(%rbp), %rax
movsd 0x10(%rax), %xmm0
subsd -0xa8(%rbp), %xmm0
movsd %xmm0, -0xc0(%rbp)
movl -0xec(%rbp), %eax
orl $0x3003f80, %eax        # imm = 0x3003F80
movl %eax, -0xec(%rbp)
movl -0xf0(%rbp), %eax
orl $0x1000080, %eax        # imm = 0x1000080
movl %eax, -0xf0(%rbp)
movl $0x7, %eax
subl -0x7c(%rbp), %eax
movl %eax, -0x7c(%rbp)
cmpl $0x0, -0x7c(%rbp)
jge 0x295c5
movl $0x0, -0x7c(%rbp)
jmp 0x295d7
movsd -0xa8(%rbp), %xmm0
movsd %xmm0, -0xc0(%rbp)
cmpl $0x4, -0x80(%rbp)
jl 0x29635
movq -0x108(%rbp), %rax
movsd 0x18(%rax), %xmm0
subsd -0xb0(%rbp), %xmm0
movsd %xmm0, -0xc8(%rbp)
movl -0xec(%rbp), %eax
orl $0x181fc000, %eax       # imm = 0x181FC000
movl %eax, -0xec(%rbp)
movl -0xf0(%rbp), %eax
orl $0x8004000, %eax        # imm = 0x8004000
movl %eax, -0xf0(%rbp)
movl $0x7, %eax
subl -0x80(%rbp), %eax
movl %eax, -0x80(%rbp)
cmpl $0x0, -0x80(%rbp)
jge 0x29633
movl $0x0, -0x80(%rbp)
jmp 0x29645
movsd -0xb0(%rbp), %xmm0
movsd %xmm0, -0xc8(%rbp)
movq -0x108(%rbp), %rax
movq (%rax), %rdi
movq 0x10(%rbp), %rax
movsd (%rax), %xmm0
callq 0x2a1a0
movsd %xmm0, -0x68(%rbp)
movsd -0xb8(%rbp), %xmm0
mulsd -0xb8(%rbp), %xmm0
ucomisd -0x68(%rbp), %xmm0
jbe 0x296ab
movsd -0xc0(%rbp), %xmm0
mulsd -0xc0(%rbp), %xmm0
ucomisd -0x68(%rbp), %xmm0
jbe 0x296ab
movsd -0xc8(%rbp), %xmm0
mulsd -0xc8(%rbp), %xmm0
ucomisd -0x68(%rbp), %xmm0
jbe 0x296ab
jmp 0x29f2a
movq -0x108(%rbp), %rax
movl -0x78(%rbp), %ecx
movl -0x7c(%rbp), %edx
movl -0x80(%rbp), %esi
shll $0x2, %esi
addl %esi, %edx
shll $0x2, %edx
addl %edx, %ecx
movl %ecx, -0x40(%rbp)
movq 0x80(%rax), %rcx
movl -0x40(%rbp), %edx
shll $0x6, %edx
movslq %edx, %rdx
shlq $0x3, %rdx
addq %rdx, %rcx
movq %rcx, -0xd0(%rbp)
movq 0x78(%rax), %rax
movl -0x40(%rbp), %ecx
shll $0x6, %ecx
movslq %ecx, %rcx
shlq $0x2, %rcx
addq %rcx, %rax
movq %rax, -0xe0(%rbp)
movq -0xe0(%rbp), %rax
movl (%rax), %eax
movl %eax, -0x90(%rbp)
movl $0x0, -0x94(%rbp)
movq -0x108(%rbp), %rax
movq (%rax), %rdi
movq 0x10(%rbp), %rax
movsd (%rax), %xmm0
callq 0x2a1a0
movaps %xmm0, %xmm1
movq -0xd0(%rbp), %rax
movslq -0x94(%rbp), %rcx
movsd (%rax,%rcx,8), %xmm0
ucomisd %xmm1, %xmm0
jbe 0x29752
jmp 0x29f2a
movl -0x94(%rbp), %eax
addl $0x1, %eax
movl %eax, -0x94(%rbp)
movq -0xe0(%rbp), %rax
movslq -0x94(%rbp), %rcx
movl (%rax,%rcx,4), %eax
movl %eax, -0xd4(%rbp)
movl -0xec(%rbp), %eax
xorl -0xd4(%rbp), %eax
movl %eax, -0xd4(%rbp)
movl -0xf0(%rbp), %eax
addl -0xd4(%rbp), %eax
movl %eax, -0xd4(%rbp)
movl -0xd4(%rbp), %eax
andl $0x7f, %eax
movl %eax, -0x78(%rbp)
movl -0x78(%rbp), %eax
subl $0x40, %eax
movl %eax, -0x78(%rbp)
movl -0xd4(%rbp), %eax
shrl $0x7, %eax
andl $0x7f, %eax
movl %eax, -0x7c(%rbp)
movl -0x7c(%rbp), %eax
subl $0x40, %eax
movl %eax, -0x7c(%rbp)
movl -0xd4(%rbp), %eax
shrl $0xe, %eax
andl $0x7f, %eax
movl %eax, -0x80(%rbp)
movl -0x80(%rbp), %eax
subl $0x40, %eax
movl %eax, -0x80(%rbp)
movl -0x78(%rbp), %eax
addl -0x6c(%rbp), %eax
movl %eax, -0x84(%rbp)
cmpl $0x0, -0x84(%rbp)
jl 0x29808
movq -0x108(%rbp), %rcx
movl -0x84(%rbp), %eax
cmpl 0x38(%rcx), %eax
jl 0x2980d
jmp 0x29943
movl -0x7c(%rbp), %eax
addl -0x70(%rbp), %eax
movl %eax, -0x88(%rbp)
cmpl $0x0, -0x88(%rbp)
jl 0x29834
movq -0x108(%rbp), %rcx
movl -0x88(%rbp), %eax
cmpl 0x3c(%rcx), %eax
jl 0x29839
jmp 0x29943
movl -0x80(%rbp), %eax
addl -0x74(%rbp), %eax
movl %eax, -0x8c(%rbp)
cmpl $0x0, -0x8c(%rbp)
jl 0x29860
movq -0x108(%rbp), %rcx
movl -0x8c(%rbp), %eax
cmpl 0x40(%rcx), %eax
jl 0x29865
jmp 0x29943
movq -0x108(%rbp), %rdi
movl -0x78(%rbp), %esi
movl -0x7c(%rbp), %edx
movl -0x80(%rbp), %ecx
movsd -0xa0(%rbp), %xmm0
movsd -0xa8(%rbp), %xmm1
movsd -0xb0(%rbp), %xmm2
movq 0x10(%rbp), %rax
movsd (%rax), %xmm3
callq 0x2a1c0
testb $0x1, %al
jne 0x298a0
jmp 0x298a5
jmp 0x29943
movq -0x108(%rbp), %rax
movq (%rax), %rdi
movl -0x34(%rbp), %esi
movl -0x38(%rbp), %edx
movl -0x3c(%rbp), %ecx
movl -0x84(%rbp), %r8d
movl -0x88(%rbp), %r9d
movl -0x8c(%rbp), %r14d
leaq -0x50(%rbp), %rbx
leaq -0x58(%rbp), %r11
leaq -0x60(%rbp), %r10
leaq -0x98(%rbp), %rax
movl %r14d, (%rsp)
movq %rbx, 0x8(%rsp)
movq %r11, 0x10(%rsp)
movq %r10, 0x18(%rsp)
movq %rax, 0x20(%rsp)
callq 0x23640
movq -0x108(%rbp), %rdi
movl %eax, -0x40(%rbp)
movl -0x40(%rbp), %esi
movsd -0x20(%rbp), %xmm0
subsd -0x50(%rbp), %xmm0
movsd -0x28(%rbp), %xmm1
subsd -0x58(%rbp), %xmm1
movsd -0x30(%rbp), %xmm2
subsd -0x60(%rbp), %xmm2
movl -0x78(%rbp), %edx
movl -0x7c(%rbp), %ecx
movl -0x80(%rbp), %r8d
movq -0x48(%rbp), %r9
movq 0x10(%rbp), %rax
movq %rax, (%rsp)
callq 0x2a010
movl -0x94(%rbp), %eax
cmpl -0x90(%rbp), %eax
jl 0x2971a
movq -0x108(%rbp), %rax
movl 0x70(%rax), %ecx
addl $0x1, %ecx
movl %ecx, 0x70(%rax)
cmpl $0x0, 0x70(%rax)
jne 0x29985
movq -0x108(%rbp), %rdi
callq 0x1f710
movq -0x108(%rbp), %rax
movl $0x1, 0x70(%rax)
movq -0x108(%rbp), %rax
movq 0x90(%rax), %rcx
movq %rcx, -0xf8(%rbp)
movq 0x90(%rax), %rax
movq %rax, -0x100(%rbp)
cmpl $0x3f, -0x94(%rbp)
jge 0x29c9d
movq -0x108(%rbp), %rax
movq (%rax), %rdi
movq 0x10(%rbp), %rax
movsd (%rax), %xmm0
callq 0x2a1a0
movaps %xmm0, %xmm1
movq -0xd0(%rbp), %rax
movslq -0x94(%rbp), %rcx
movsd (%rax,%rcx,8), %xmm0
ucomisd %xmm1, %xmm0
jbe 0x299ed
jmp 0x29f2a
movl -0x94(%rbp), %eax
addl $0x1, %eax
movl %eax, -0x94(%rbp)
movq -0xe0(%rbp), %rax
movslq -0x94(%rbp), %rcx
movl (%rax,%rcx,4), %eax
movl %eax, -0xd4(%rbp)
movl -0xec(%rbp), %eax
xorl -0xd4(%rbp), %eax
movl %eax, -0xd4(%rbp)
movl -0xf0(%rbp), %eax
addl -0xd4(%rbp), %eax
movl %eax, -0xd4(%rbp)
movl -0xd4(%rbp), %eax
andl $0x7f, %eax
movl %eax, -0x78(%rbp)
movl -0x78(%rbp), %eax
subl $0x40, %eax
movl %eax, -0x78(%rbp)
movl -0xd4(%rbp), %eax
shrl $0x7, %eax
andl $0x7f, %eax
movl %eax, -0x7c(%rbp)
movl -0x7c(%rbp), %eax
subl $0x40, %eax
movl %eax, -0x7c(%rbp)
movl -0xd4(%rbp), %eax
shrl $0xe, %eax
andl $0x7f, %eax
movl %eax, -0x80(%rbp)
movl -0x80(%rbp), %eax
subl $0x40, %eax
movl %eax, -0x80(%rbp)
movl -0x78(%rbp), %eax
addl -0x6c(%rbp), %eax
movl %eax, -0x84(%rbp)
cmpl $0x0, -0x84(%rbp)
jl 0x29aa3
movq -0x108(%rbp), %rcx
movl -0x84(%rbp), %eax
cmpl 0x38(%rcx), %eax
jl 0x29aa8
jmp 0x299a8
movl -0x7c(%rbp), %eax
addl -0x70(%rbp), %eax
movl %eax, -0x88(%rbp)
cmpl $0x0, -0x88(%rbp)
jl 0x29acf
movq -0x108(%rbp), %rcx
movl -0x88(%rbp), %eax
cmpl 0x3c(%rcx), %eax
jl 0x29ad4
jmp 0x299a8
movl -0x80(%rbp), %eax
addl -0x74(%rbp), %eax
movl %eax, -0x8c(%rbp)
cmpl $0x0, -0x8c(%rbp)
jl 0x29afb
movq -0x108(%rbp), %rcx
movl -0x8c(%rbp), %eax
cmpl 0x40(%rcx), %eax
jl 0x29b00
jmp 0x299a8
movq -0x108(%rbp), %rdi
movq 0x88(%rdi), %rax
movslq -0x84(%rbp), %rcx
shlq $0x2, %rcx
addq %rcx, %rax
movl 0x38(%rdi), %ecx
movl -0x88(%rbp), %edx
movl 0x3c(%rdi), %esi
imull -0x8c(%rbp), %esi
addl %esi, %edx
imull %edx, %ecx
movslq %ecx, %rcx
shlq $0x2, %rcx
addq %rcx, %rax
movq %rax, -0xe8(%rbp)
movl 0x70(%rdi), %ecx
movq -0xe8(%rbp), %rax
movl %ecx, (%rax)
movl -0x78(%rbp), %esi
movl -0x7c(%rbp), %edx
movl -0x80(%rbp), %ecx
movsd -0xa0(%rbp), %xmm0
movsd -0xa8(%rbp), %xmm1
movsd -0xb0(%rbp), %xmm2
movq 0x10(%rbp), %rax
movsd (%rax), %xmm3
callq 0x2a1c0
testb $0x1, %al
jne 0x29b85
jmp 0x29b8a
jmp 0x299a8
movq -0x108(%rbp), %rax
movq (%rax), %rdi
movl -0x34(%rbp), %esi
movl -0x38(%rbp), %edx
movl -0x3c(%rbp), %ecx
movl -0x84(%rbp), %r8d
movl -0x88(%rbp), %r9d
movl -0x8c(%rbp), %r14d
leaq -0x50(%rbp), %rbx
leaq -0x58(%rbp), %r11
leaq -0x60(%rbp), %r10
leaq -0x98(%rbp), %rax
movl %r14d, (%rsp)
movq %rbx, 0x8(%rsp)
movq %r11, 0x10(%rsp)
movq %r10, 0x18(%rsp)
movq %rax, 0x20(%rsp)
callq 0x23640
movq -0x108(%rbp), %rdi
movl %eax, -0x40(%rbp)
movl -0x40(%rbp), %esi
movsd -0x20(%rbp), %xmm0
subsd -0x50(%rbp), %xmm0
movsd -0x28(%rbp), %xmm1
subsd -0x58(%rbp), %xmm1
movsd -0x30(%rbp), %xmm2
subsd -0x60(%rbp), %xmm2
movl -0x78(%rbp), %edx
movl -0x7c(%rbp), %ecx
movl -0x80(%rbp), %r8d
movq -0x48(%rbp), %r9
movq 0x10(%rbp), %rax
movq %rax, (%rsp)
callq 0x2a010
movq -0x108(%rbp), %rcx
movq -0x100(%rbp), %rax
movq 0x98(%rcx), %rcx
addq $-0x48, %rcx
cmpq %rcx, %rax
jbe 0x29c60
movq -0x108(%rbp), %rdi
leaq -0xf8(%rbp), %rsi
leaq -0x100(%rbp), %rdx
callq 0x238b0
movq -0x108(%rbp), %rdi
movl -0xd4(%rbp), %esi
movq -0xe8(%rbp), %rdx
movl -0x84(%rbp), %ecx
movl -0x88(%rbp), %r8d
movl -0x8c(%rbp), %r9d
leaq -0x100(%rbp), %rax
movq %rax, (%rsp)
callq 0x23a60
jmp 0x299a8
movq -0x108(%rbp), %rax
movq (%rax), %rdi
movq 0x10(%rbp), %rax
movsd (%rax), %xmm0
callq 0x2a1a0
movaps %xmm0, %xmm1
movq -0xd0(%rbp), %rax
movslq -0x94(%rbp), %rcx
movsd (%rax,%rcx,8), %xmm0
ucomisd %xmm1, %xmm0
jbe 0x29cd5
jmp 0x29f2a
jmp 0x29cd7
movq -0xf8(%rbp), %rax
cmpq -0x100(%rbp), %rax
je 0x29f2a
movq -0x108(%rbp), %rcx
movq -0xf8(%rbp), %rax
cmpq 0x98(%rcx), %rax
jne 0x29d17
movq -0x108(%rbp), %rax
movq 0x90(%rax), %rax
movq %rax, -0xf8(%rbp)
movq -0x108(%rbp), %rdi
movq -0xf8(%rbp), %rax
movq %rax, %rcx
addq $0x4, %rcx
movq %rcx, -0xf8(%rbp)
movl (%rax), %eax
movl %eax, -0x84(%rbp)
movq -0xf8(%rbp), %rax
movq %rax, %rcx
addq $0x4, %rcx
movq %rcx, -0xf8(%rbp)
movl (%rax), %eax
movl %eax, -0x88(%rbp)
movq -0xf8(%rbp), %rax
movq %rax, %rcx
addq $0x4, %rcx
movq %rcx, -0xf8(%rbp)
movl (%rax), %eax
movl %eax, -0x8c(%rbp)
movl -0x84(%rbp), %eax
subl -0x6c(%rbp), %eax
movl %eax, -0x78(%rbp)
movl -0x88(%rbp), %eax
subl -0x70(%rbp), %eax
movl %eax, -0x7c(%rbp)
movl -0x8c(%rbp), %eax
subl -0x74(%rbp), %eax
movl %eax, -0x80(%rbp)
movl -0x78(%rbp), %esi
movl -0x7c(%rbp), %edx
movl -0x80(%rbp), %ecx
movsd -0xa0(%rbp), %xmm0
movsd -0xa8(%rbp), %xmm1
movsd -0xb0(%rbp), %xmm2
movq 0x10(%rbp), %rax
movsd (%rax), %xmm3
callq 0x2a1c0
testb $0x1, %al
jne 0x29dcd
jmp 0x29dd2
jmp 0x29cd7
movq -0x108(%rbp), %rax
movq (%rax), %rdi
movl -0x34(%rbp), %esi
movl -0x38(%rbp), %edx
movl -0x3c(%rbp), %ecx
movl -0x84(%rbp), %r8d
movl -0x88(%rbp), %r9d
movl -0x8c(%rbp), %r14d
leaq -0x50(%rbp), %rbx
leaq -0x58(%rbp), %r11
leaq -0x60(%rbp), %r10
leaq -0x98(%rbp), %rax
movl %r14d, (%rsp)
movq %rbx, 0x8(%rsp)
movq %r11, 0x10(%rsp)
movq %r10, 0x18(%rsp)
movq %rax, 0x20(%rsp)
callq 0x23640
movq -0x108(%rbp), %rdi
movl %eax, -0x40(%rbp)
movl -0x40(%rbp), %esi
movsd -0x20(%rbp), %xmm0
subsd -0x50(%rbp), %xmm0
movsd -0x28(%rbp), %xmm1
subsd -0x58(%rbp), %xmm1
movsd -0x30(%rbp), %xmm2
subsd -0x60(%rbp), %xmm2
movl -0x78(%rbp), %edx
movl -0x7c(%rbp), %ecx
movl -0x80(%rbp), %r8d
movq -0x48(%rbp), %r9
movq 0x10(%rbp), %rax
movq %rax, (%rsp)
callq 0x2a010
movq -0xf8(%rbp), %rax
cmpq -0x100(%rbp), %rax
ja 0x29ebd
movq -0x108(%rbp), %rdx
movq 0x98(%rdx), %rax
movq -0x100(%rbp), %rcx
subq %rcx, %rax
sarq $0x2, %rax
movq -0xf8(%rbp), %rcx
movq 0x90(%rdx), %rdx
subq %rdx, %rcx
sarq $0x2, %rcx
addq %rcx, %rax
movq %rax, -0x110(%rbp)
jmp 0x29ed9
movq -0xf8(%rbp), %rax
movq -0x100(%rbp), %rcx
subq %rcx, %rax
sarq $0x2, %rax
movq %rax, -0x110(%rbp)
movq -0x110(%rbp), %rax
cmpq $0x12, %rax
jge 0x29f00
movq -0x108(%rbp), %rdi
leaq -0xf8(%rbp), %rsi
leaq -0x100(%rbp), %rdx
callq 0x238b0
movq -0x108(%rbp), %rdi
movl -0x84(%rbp), %esi
movl -0x88(%rbp), %edx
movl -0x8c(%rbp), %ecx
leaq -0x100(%rbp), %r8
callq 0x25470
jmp 0x29cd7
addq $0x130, %rsp            # imm = 0x130
popq %rbx
popq %r14
popq %rbp
retq
nopw %cs:(%rax,%rax)