// File: /workspace/github_cpp_makefile_star2_O3/Luthaf[P]zeopp/voro++/src/v_compute.cc
// Function: bool voro::voro_compute<voro::container_poly>::compute_cell<voro::voronoicell_neighbor>(voro::voronoicell_neighbor&, int, int, int, int, int)
bool voro_compute<c_class>::compute_cell(v_cell &c,int ijk,int s,int ci,int cj,int ck) {
	static const int count_list[8]={7,11,15,19,26,35,45,59},*count_e=count_list+8;
	double x,y,z,x1,y1,z1,qx=0,qy=0,qz=0;
	double xlo,ylo,zlo,xhi,yhi,zhi,x2,y2,z2,rs;
	int i,j,k,di,dj,dk,ei,ej,ek,f,g,l,disp;
	double fx,fy,fz,gxs,gys,gzs,*radp;
	unsigned int q,*e,*mijk;

	if(!con.initialize_voronoicell(c,ijk,s,ci,cj,ck,i,j,k,x,y,z,disp)) return false;
	con.r_init(ijk,s);

	// Initialize the Voronoi cell to fill the entire container
	double crs,mrs;

	int next_count=3,*count_p=(const_cast<int*> (count_list));

	// Test all particles in the particle's local region first
	for(l=0;l<s;l++) {
		x1=p[ijk][ps*l]-x;
		y1=p[ijk][ps*l+1]-y;
		z1=p[ijk][ps*l+2]-z;
		rs=con.r_scale(x1*x1+y1*y1+z1*z1,ijk,l);
		if(!c.nplane(x1,y1,z1,rs,id[ijk][l])) return false;
	}
	l++;
	while(l<co[ijk]) {
		x1=p[ijk][ps*l]-x;
		y1=p[ijk][ps*l+1]-y;
		z1=p[ijk][ps*l+2]-z;
		rs=con.r_scale(x1*x1+y1*y1+z1*z1,ijk,l);
		if(!c.nplane(x1,y1,z1,rs,id[ijk][l])) return false;
		l++;
	}

	// Now compute the maximum distance squared from the cell center to a
	// vertex. This is used to cut off the calculation since we only need
	// to test out to twice this range.
	mrs=c.max_radius_squared();

	// 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) {
		gxs=fx;
		m1=127+(3<<21);m2=1+(1<<21);di=wl_fgrid-1-di;if(di<0) di=0;
	} else {m1=m2=0;gxs=boxx-fx;}
	if(dj>=wl_hgrid) {
		gys=fy;
		m1|=(127<<7)+(3<<24);m2|=(1<<7)+(1<<24);dj=wl_fgrid-1-dj;if(dj<0) dj=0;
	} else gys=boxy-fy;
	if(dk>=wl_hgrid) {
		gzs=fz;
		m1|=(127<<14)+(3<<27);m2|=(1<<14)+(1<<27);dk=wl_fgrid-1-dk;if(dk<0) dk=0;
	} else gzs=boxz-fz;
	gxs*=gxs;gys*=gys;gzs*=gzs;

	// 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 {

		// At the intervals specified by count_list, we recompute the
		// maximum radius squared
		if(g==next_count) {
			mrs=c.max_radius_squared();
			if(count_p!=count_e) next_count=*(count_p++);
		}

		// If mrs is less than the minimum distance to any untested
		// block, then we are done
		if(con.r_ctest(radp[g],mrs)) return true;
		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_max_radius(di,dj,dk,fx,fy,fz,gxs,gys,gzs,crs,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.
		if(co[ijk]>0) {
			l=0;x2=x-qx;y2=y-qy;z2=z-qz;
			if(!con.r_ctest(crs,mrs)) {
				do {
					x1=p[ijk][ps*l]-x2;
					y1=p[ijk][ps*l+1]-y2;
					z1=p[ijk][ps*l+2]-z2;
					rs=con.r_scale(x1*x1+y1*y1+z1*z1,ijk,l);
					if(!c.nplane(x1,y1,z1,rs,id[ijk][l])) return false;
					l++;
				} while (l<co[ijk]);
			} else {
				do {
					x1=p[ijk][ps*l]-x2;
					y1=p[ijk][ps*l+1]-y2;
					z1=p[ijk][ps*l+2]-z2;
					rs=x1*x1+y1*y1+z1*z1;
					if(con.r_scale_check(rs,mrs,ijk,l)&&!c.nplane(x1,y1,z1,rs,id[ijk][l])) return false;
					l++;
				} while (l<co[ijk]);
			}
		}
	} while(g<f);

	// If we reach here, we were unable to compute the entire cell using
	// the first part of the worklist. This section of the algorithm
	// continues the worklist, but it now starts preparing the mask that we
	// need if we end up going block by block. We do the same as before,
	// but we put a mark down on the mask for every block that's tested.
	// The worklist also contains information about which neighbors of each
	// block are not also on the worklist, and we start storing those
	// points in a list in case we have to go block by block. Update the
	// mask counter, and if it wraps around then reset the whole mask; that
	// will only happen once every 2^32 tries.
	mv++;
	if(mv==0) {reset_mask();mv=1;}

	// Set the queue pointers
	int *qu_s=qu,*qu_e=qu;

	while(g<wl_seq_length-1) {

		// At the intervals specified by count_list, we recompute the
		// maximum radius squared
		if(g==next_count) {
			mrs=c.max_radius_squared();
			if(count_p!=count_e) next_count=*(count_p++);
		}

		// If mrs is less than the minimum distance to any untested
		// block, then we are done
		if(con.r_ctest(radp[g],mrs)) return true;
		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;

		// 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_max_radius(di,dj,dk,fx,fy,fz,gxs,gys,gzs,crs,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.
		if(co[ijk]>0) {
			l=0;x2=x-qx;y2=y-qy;z2=z-qz;
			if(!con.r_ctest(crs,mrs)) {
				do {
					x1=p[ijk][ps*l]-x2;
					y1=p[ijk][ps*l+1]-y2;
					z1=p[ijk][ps*l+2]-z2;
					rs=con.r_scale(x1*x1+y1*y1+z1*z1,ijk,l);
					if(!c.nplane(x1,y1,z1,rs,id[ijk][l])) return false;
					l++;
				} while (l<co[ijk]);
			} else {
				do {
					x1=p[ijk][ps*l]-x2;
					y1=p[ijk][ps*l+1]-y2;
					z1=p[ijk][ps*l+2]-z2;
					rs=x1*x1+y1*y1+z1*z1;
					if(con.r_scale_check(rs,mrs,ijk,l)&&!c.nplane(x1,y1,z1,rs,id[ijk][l])) return false;
					l++;
				} while (l<co[ijk]);
			}
		}

		// If there might not be enough memory on the list for these
		// additions, then add more
		if(qu_e>qu_l-18) add_list_memory(qu_s,qu_e);

		// Test the parts of the worklist element which tell us what
		// neighbors of this block are not on the worklist. Store them
		// on the block list, and mark the mask.
		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_ctest(radp[g],mrs)) return true;

	// 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) {

		// If we reached the end of the list memory loop back to the
		// start
		if(qu_s==qu_l) qu_s=qu;

		// Read in a block off the list, and compute the upper and lower
		// coordinates in each of the three dimensions
		ei=*(qu_s++);ej=*(qu_s++);ek=*(qu_s++);
		xlo=(ei-i)*boxx-fx;xhi=xlo+boxx;
		ylo=(ej-j)*boxy-fy;yhi=ylo+boxy;
		zlo=(ek-k)*boxz-fz;zhi=zlo+boxz;

		// Carry out plane tests to see if any particle in this block
		// could possibly intersect the cell
		if(ei>i) {
			if(ej>j) {
				if(ek>k) {if(corner_test(c,xlo,ylo,zlo,xhi,yhi,zhi)) continue;}
				else if(ek<k) {if(corner_test(c,xlo,ylo,zhi,xhi,yhi,zlo)) continue;}
				else {if(edge_z_test(c,xlo,ylo,zlo,xhi,yhi,zhi)) continue;}
			} else if(ej<j) {
				if(ek>k) {if(corner_test(c,xlo,yhi,zlo,xhi,ylo,zhi)) continue;}
				else if(ek<k) {if(corner_test(c,xlo,yhi,zhi,xhi,ylo,zlo)) continue;}
				else {if(edge_z_test(c,xlo,yhi,zlo,xhi,ylo,zhi)) continue;}
			} else {
				if(ek>k) {if(edge_y_test(c,xlo,ylo,zlo,xhi,yhi,zhi)) continue;}
				else if(ek<k) {if(edge_y_test(c,xlo,ylo,zhi,xhi,yhi,zlo)) continue;}
				else {if(face_x_test(c,xlo,ylo,zlo,yhi,zhi)) continue;}
			}
		} else if(ei<i) {
			if(ej>j) {
				if(ek>k) {if(corner_test(c,xhi,ylo,zlo,xlo,yhi,zhi)) continue;}
				else if(ek<k) {if(corner_test(c,xhi,ylo,zhi,xlo,yhi,zlo)) continue;}
				else {if(edge_z_test(c,xhi,ylo,zlo,xlo,yhi,zhi)) continue;}
			} else if(ej<j) {
				if(ek>k) {if(corner_test(c,xhi,yhi,zlo,xlo,ylo,zhi)) continue;}
				else if(ek<k) {if(corner_test(c,xhi,yhi,zhi,xlo,ylo,zlo)) continue;}
				else {if(edge_z_test(c,xhi,yhi,zlo,xlo,ylo,zhi)) continue;}
			} else {
				if(ek>k) {if(edge_y_test(c,xhi,ylo,zlo,xlo,yhi,zhi)) continue;}
				else if(ek<k) {if(edge_y_test(c,xhi,ylo,zhi,xlo,yhi,zlo)) continue;}
				else {if(face_x_test(c,xhi,ylo,zlo,yhi,zhi)) continue;}
			}
		} else {
			if(ej>j) {
				if(ek>k) {if(edge_x_test(c,xlo,ylo,zlo,xhi,yhi,zhi)) continue;}
				else if(ek<k) {if(edge_x_test(c,xlo,ylo,zhi,xhi,yhi,zlo)) continue;}
				else {if(face_y_test(c,xlo,ylo,zlo,xhi,zhi)) continue;}
			} else if(ej<j) {
				if(ek>k) {if(edge_x_test(c,xlo,yhi,zlo,xhi,ylo,zhi)) continue;}
				else if(ek<k) {if(edge_x_test(c,xlo,yhi,zhi,xhi,ylo,zlo)) continue;}
				else {if(face_y_test(c,xlo,yhi,zlo,xhi,zhi)) continue;}
			} else {
				if(ek>k) {if(face_z_test(c,xlo,ylo,zlo,xhi,yhi)) continue;}
				else if(ek<k) {if(face_z_test(c,xlo,ylo,zhi,xhi,yhi)) continue;}
				else voro_fatal_error("Compute cell routine revisiting central block, which should never\nhappen.",VOROPP_INTERNAL_ERROR);
			}
		}

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

		// Loop over all the elements in the block to test for cuts. It
		// would be possible to exclude some of these cases by testing
		// against mrs, but this will probably not save time.
		if(co[ijk]>0) {
			l=0;x2=x-qx;y2=y-qy;z2=z-qz;
			do {
				x1=p[ijk][ps*l]-x2;
				y1=p[ijk][ps*l+1]-y2;
				z1=p[ijk][ps*l+2]-z2;
				rs=con.r_scale(x1*x1+y1*y1+z1*z1,ijk,l);
				if(!c.nplane(x1,y1,z1,rs,id[ijk][l])) return false;
				l++;
			} while (l<co[ijk]);
		}

		// If there's not much memory on the block list then add more
		if((qu_s<=qu_e?(qu_l-qu_e)+(qu_s-qu):qu_s-qu_e)<18) add_list_memory(qu_s,qu_e);

		// Test the neighbors of the current block, and add them to the
		// block list if they haven't already been tested
		add_to_mask(ei,ej,ek,qu_e);
	}

	return true;
}

// Assembly:
pushq %rbp
pushq %r15
pushq %r14
pushq %r13
pushq %r12
pushq %rbx
subq $0x168, %rsp            # imm = 0x168
movq %r9, 0xb8(%rsp)
movl %ecx, %ebx
movq %rsi, %r15
movq %rdi, %r12
xorl %eax, %eax
movq %rax, 0xa0(%rsp)
movq %rax, 0x98(%rsp)
movq %rax, 0x90(%rsp)
movq (%rdi), %r13
movq 0xb8(%r13), %rax
movslq %edx, %r14
movq (%rax,%r14,8), %rax
movslq 0xd0(%r13), %rdx
movslq %ecx, %rcx
imulq %rdx, %rcx
movsd (%rax,%rcx,8), %xmm2
cmpb $0x1, 0xa8(%r13)
movq %r8, 0xd0(%rsp)
movsd %xmm2, 0x60(%rsp)
jne 0x1fa01
movsd 0x78(%r13), %xmm1
subsd 0x70(%r13), %xmm1
mulsd 0x17645(%rip), %xmm1    # 0x37030
movapd 0x1764d(%rip), %xmm0    # 0x37040
xorpd %xmm1, %xmm0
movl (%r13), %edx
movl %edx, 0x14(%rsp)
jmp 0x1fa1a
movsd 0x70(%r13), %xmm0
movsd 0x78(%r13), %xmm1
subsd %xmm2, %xmm0
subsd %xmm2, %xmm1
movl %r8d, 0x14(%rsp)
movsd 0x8(%rax,%rcx,8), %xmm4
cmpb $0x1, 0xa9(%r13)
jne 0x1fa56
movsd 0x88(%r13), %xmm3
subsd 0x80(%r13), %xmm3
mulsd 0x175ec(%rip), %xmm3    # 0x37030
movapd 0x175f4(%rip), %xmm2    # 0x37040
xorpd %xmm3, %xmm2
movl 0x4(%r13), %edx
jmp 0x1fa78
movsd 0x80(%r13), %xmm2
movsd 0x88(%r13), %xmm3
subsd %xmm4, %xmm2
subsd %xmm4, %xmm3
movq 0xb8(%rsp), %rdx
movl %edx, 0x10(%rsp)
movsd %xmm4, 0x30(%rsp)
movl 0x1a0(%rsp), %edx
movsd 0x10(%rax,%rcx,8), %xmm6
cmpb $0x1, 0xaa(%r13)
movsd %xmm6, 0x58(%rsp)
jne 0x1facf
movsd 0x98(%r13), %xmm5
subsd 0x90(%r13), %xmm5
mulsd 0x17577(%rip), %xmm5    # 0x37030
movapd 0x1757f(%rip), %xmm4    # 0x37040
xorpd %xmm5, %xmm4
movl 0x8(%r13), %eax
movl %eax, 0x4(%rsp)
jmp 0x1faed
movsd 0x90(%r13), %xmm4
movsd 0x98(%r13), %xmm5
subsd %xmm6, %xmm4
subsd %xmm6, %xmm5
movl %edx, 0x4(%rsp)
movq %r15, %rdi
callq 0xb6ec
movq 0x50(%r13), %rbp
cmpq 0x58(%r13), %rbp
jae 0x1fb2c
movq (%rbp), %rdi
movq (%rdi), %rax
movq %r15, %rsi
movsd 0x60(%rsp), %xmm0
movsd 0x30(%rsp), %xmm1
movsd 0x58(%rsp), %xmm2
callq *0x20(%rax)
testb %al, %al
je 0x1fd1d
addq $0x8, %rbp
jmp 0x1faf9
movl 0x4(%r13), %eax
imull 0x4(%rsp), %eax
addl 0x10(%rsp), %eax
imull (%r13), %eax
addl 0x14(%rsp), %eax
movl %r14d, %ecx
subl %eax, %ecx
movl %ecx, 0x6c(%rsp)
movq (%r12), %rax
movq 0xd8(%rax), %rcx
movq (%rcx,%r14,8), %rcx
leal 0x3(,%rbx,4), %edx
movslq %edx, %rdx
movsd (%rcx,%rdx,8), %xmm0
mulsd %xmm0, %xmm0
movsd %xmm0, 0xe8(%rax)
movsd 0xe0(%rax), %xmm1
mulsd %xmm1, %xmm1
subsd %xmm1, %xmm0
movsd %xmm0, 0xf0(%rax)
testl %ebx, %ebx
movsd 0x30(%rsp), %xmm2
jle 0x1fc57
movl %ebx, %ebx
movl $0x18, %r13d
xorl %ebp, %ebp
movq 0x58(%r12), %rax
movq (%rax,%r14,8), %rax
movslq 0x4c(%r12), %rcx
movslq %ebp, %rdx
imulq %rcx, %rdx
movsd (%rax,%rdx,8), %xmm0
movsd 0x8(%rax,%rdx,8), %xmm1
subsd 0x60(%rsp), %xmm0
subsd %xmm2, %xmm1
movsd 0x10(%rax,%rdx,8), %xmm2
subsd 0x58(%rsp), %xmm2
movapd %xmm1, %xmm3
mulsd %xmm1, %xmm3
movapd %xmm0, %xmm4
mulsd %xmm0, %xmm4
addsd %xmm3, %xmm4
movapd %xmm2, %xmm3
mulsd %xmm2, %xmm3
addsd %xmm4, %xmm3
movq (%r12), %rax
addsd 0xe8(%rax), %xmm3
movq 0x50(%r12), %rcx
movq 0xd8(%rax), %rax
movq (%rax,%r14,8), %rax
movsd (%rax,%r13), %xmm4
mulsd %xmm4, %xmm4
subsd %xmm4, %xmm3
movq (%rcx,%r14,8), %rax
movl (%rax,%rbp,4), %edx
movq %r15, %rdi
movq %r15, %rsi
callq 0xe636
testb %al, %al
je 0x1fd1d
incq %rbp
addq $0x20, %r13
cmpq %rbp, %rbx
movsd 0x30(%rsp), %xmm2
jne 0x1fba5
jmp 0x1fc59
xorl %ebx, %ebx
movq %rbx, %r13
shlq $0x5, %r13
incq %rbx
movq 0x60(%r12), %rax
movslq (%rax,%r14,4), %rax
cmpq %rax, %rbx
jge 0x1fd33
movq 0x58(%r12), %rax
movq (%rax,%r14,8), %rax
movslq 0x4c(%r12), %rcx
movslq %ebx, %rdx
imulq %rcx, %rdx
movsd (%rax,%rdx,8), %xmm0
movsd 0x8(%rax,%rdx,8), %xmm1
subsd 0x60(%rsp), %xmm0
subsd %xmm2, %xmm1
movsd 0x10(%rax,%rdx,8), %xmm2
subsd 0x58(%rsp), %xmm2
movapd %xmm1, %xmm3
mulsd %xmm1, %xmm3
movapd %xmm0, %xmm4
mulsd %xmm0, %xmm4
addsd %xmm3, %xmm4
movapd %xmm2, %xmm3
mulsd %xmm2, %xmm3
addsd %xmm4, %xmm3
movq (%r12), %rax
addsd 0xe8(%rax), %xmm3
movq 0x50(%r12), %rcx
movq 0xd8(%rax), %rax
movq (%rax,%r14,8), %rax
movsd 0x38(%rax,%r13), %xmm4
mulsd %xmm4, %xmm4
subsd %xmm4, %xmm3
movq (%rcx,%r14,8), %rax
movl (%rax,%rbx,4), %edx
movq %r15, %rdi
movq %r15, %rsi
callq 0xe636
movsd 0x30(%rsp), %xmm2
addq $0x20, %r13
incq %rbx
testb %al, %al
jne 0x1fc63
xorl %ebp, %ebp
movl %ebp, %eax
addq $0x168, %rsp            # imm = 0x168
popq %rbx
popq %r12
popq %r13
popq %r14
popq %r15
popq %rbp
retq
movq %r15, %rdi
callq 0x9a14
movapd %xmm0, %xmm6
xorps %xmm0, %xmm0
cvtsi2sdl 0xd0(%rsp), %xmm0
xorps %xmm1, %xmm1
cvtsi2sdl 0xb8(%rsp), %xmm1
movq (%r12), %rax
movl 0x1a0(%rsp), %ecx
xorps %xmm2, %xmm2
cvtsi2sd %ecx, %xmm2
movsd 0x60(%rsp), %xmm3
subsd 0x70(%rax), %xmm3
mulsd 0x18(%rax), %xmm0
subsd %xmm0, %xmm3
movsd 0x30(%rsp), %xmm4
subsd 0x80(%rax), %xmm4
mulsd 0x20(%rax), %xmm1
subsd %xmm1, %xmm4
movsd 0x58(%rsp), %xmm5
subsd 0x90(%rax), %xmm5
mulsd 0x28(%rax), %xmm2
subsd %xmm2, %xmm5
movsd 0x20(%r12), %xmm0
mulsd %xmm3, %xmm0
movsd 0x19572(%rip), %xmm1    # 0x39330
mulsd %xmm1, %xmm0
cvttsd2si %xmm0, %eax
movsd 0x28(%r12), %xmm0
movapd %xmm4, 0xf0(%rsp)
mulsd %xmm4, %xmm0
mulsd %xmm1, %xmm0
cvttsd2si %xmm0, %ecx
movsd 0x30(%r12), %xmm0
movapd %xmm5, 0xe0(%rsp)
mulsd %xmm5, %xmm0
mulsd %xmm1, %xmm0
cmpl $0x4, %eax
movapd %xmm3, 0x130(%rsp)
jl 0x1fe31
movl $0x7, %edx
subl %eax, %edx
xorl %eax, %eax
testl %edx, %edx
cmovgl %edx, %eax
movl $0x200001, 0x28(%rsp)   # imm = 0x200001
movl $0x60007f, 0x18(%rsp)   # imm = 0x60007F
movsd %xmm3, 0xb0(%rsp)
jmp 0x1fe55
movsd 0x8(%r12), %xmm1
subsd %xmm3, %xmm1
movsd %xmm1, 0xb0(%rsp)
movl $0x0, 0x18(%rsp)
movl $0x0, 0x28(%rsp)
cvttsd2si %xmm0, %edx
cmpl $0x4, %ecx
jl 0x1fe87
orl $0x3003f80, 0x18(%rsp)  # imm = 0x3003F80
orl $0x1000080, 0x28(%rsp)  # imm = 0x1000080
movl $0x7, %esi
subl %ecx, %esi
xorl %ecx, %ecx
testl %esi, %esi
cmovgl %esi, %ecx
movapd 0xf0(%rsp), %xmm0
jmp 0x1fe97
movsd 0x10(%r12), %xmm0
subsd 0xf0(%rsp), %xmm0
movsd %xmm0, 0xa8(%rsp)
cmpl $0x4, %edx
jl 0x1fece
orl $0x181fc000, 0x18(%rsp) # imm = 0x181FC000
orl $0x8004000, 0x28(%rsp)  # imm = 0x8004000
movl $0x7, %esi
subl %edx, %esi
xorl %edx, %edx
testl %esi, %esi
cmovgl %esi, %edx
movapd 0xe0(%rsp), %xmm1
jmp 0x1fede
movsd 0x18(%r12), %xmm1
subsd 0xe0(%rsp), %xmm1
movsd 0xb0(%rsp), %xmm0
mulsd %xmm0, %xmm0
movsd %xmm0, 0xb0(%rsp)
shll $0xa, %edx
shll $0x8, %ecx
shll $0x6, %eax
addl %ecx, %eax
addl %edx, %eax
cltq
leaq (,%rax,8), %rcx
addq 0x80(%r12), %rcx
movq %rcx, 0x78(%rsp)
movsd 0xa8(%rsp), %xmm0
mulsd %xmm0, %xmm0
movsd %xmm0, 0xa8(%rsp)
movq 0x78(%r12), %rcx
leaq (%rcx,%rax,4), %rax
movq %rax, 0x150(%rsp)
movl (%rax), %eax
cmpl $0x2, %eax
movl $0x1, %ecx
cmovgel %eax, %ecx
movq %rcx, 0x70(%rsp)
mulsd %xmm1, %xmm1
movsd %xmm1, 0x148(%rsp)
leaq 0x1952b(%rip), %rax     # 0x39490
movq %rax, 0xd8(%rsp)
movl $0x3, 0x80(%rsp)
movl $0x1, %eax
xorl %ecx, %ecx
movq %rcx, %r13
movq %rax, %rbx
movl 0x80(%rsp), %eax
cmpq %rax, %rcx
jne 0x1ffc3
movq %r15, %rdi
callq 0x9a14
movapd %xmm0, %xmm6
movq 0xd8(%rsp), %rax
cmpq 0x251bc(%rip), %rax     # 0x45168
je 0x1ffc3
movl (%rax), %ecx
movl %ecx, 0x80(%rsp)
addq $0x4, %rax
movq %rax, 0xd8(%rsp)
movq (%r12), %rax
movq 0x78(%rsp), %rcx
movsd (%rcx,%r13,8), %xmm0
movsd 0xf0(%rax), %xmm1
addsd %xmm0, %xmm1
mulsd %xmm6, %xmm0
ucomisd 0x18ab6(%rip), %xmm0    # 0x38aa0
movsd %xmm6, 0x8(%rsp)
movq %rbx, 0x50(%rsp)
jb 0x1fffd
sqrtsd %xmm0, %xmm0
jmp 0x20014
movsd %xmm1, 0x48(%rsp)
callq 0x20f0
movsd 0x48(%rsp), %xmm1
movsd 0x8(%rsp), %xmm6
ucomisd %xmm0, %xmm1
ja 0x20acf
movq 0x150(%rsp), %rax
movl 0x4(%rax,%r13,4), %ecx
xorl 0x18(%rsp), %ecx
addl 0x28(%rsp), %ecx
movl %ecx, %esi
andl $0x7f, %esi
addl $-0x40, %esi
movl %esi, %ebx
addl 0x14(%rsp), %ebx
js 0x204dd
cmpl 0x38(%r12), %ebx
jge 0x204dd
movl %ecx, %edx
shrl $0x7, %edx
andl $0x7f, %edx
addl $-0x40, %edx
movl %edx, %ebp
addl 0x10(%rsp), %ebp
js 0x204dd
cmpl 0x3c(%r12), %ebp
jge 0x204dd
shrl $0xe, %ecx
andl $0x7f, %ecx
addl $-0x40, %ecx
movl %ecx, %r14d
addl 0x4(%rsp), %r14d
js 0x204dd
cmpl 0x40(%r12), %r14d
jge 0x204dd
movq %r12, %rdi
movapd 0x130(%rsp), %xmm0
movaps 0xf0(%rsp), %xmm1
movapd 0xe0(%rsp), %xmm2
movsd 0xb0(%rsp), %xmm3
movsd 0xa8(%rsp), %xmm4
movsd 0x148(%rsp), %xmm5
leaq 0x128(%rsp), %r8
movsd 0x8(%rsp), %xmm6
callq 0x1d98a
movsd 0x8(%rsp), %xmm6
movsd 0x30(%rsp), %xmm1
testb %al, %al
jne 0x204dd
movq (%r12), %rax
cmpb $0x1, 0xa8(%rax)
jne 0x2014b
movq 0xd0(%rsp), %rcx
leal (%rbx,%rcx), %edx
movl (%rax), %ecx
cmpl %ecx, %edx
jge 0x2012b
addl %ebx, %ecx
movsd 0x78(%rax), %xmm0
subsd 0x70(%rax), %xmm0
xorpd 0x16f19(%rip), %xmm0    # 0x37040
movl %ecx, %ebx
jmp 0x20142
leal (%rcx,%rcx), %esi
xorpd %xmm0, %xmm0
cmpl %esi, %edx
jl 0x20142
subl %ecx, %ebx
movsd 0x78(%rax), %xmm0
subsd 0x70(%rax), %xmm0
movsd %xmm0, 0xa0(%rsp)
cmpb $0x1, 0xa9(%rax)
jne 0x20184
movq 0xb8(%rsp), %rcx
leal (%rcx,%rbp), %edx
movl 0x4(%rax), %ecx
cmpl %ecx, %edx
jge 0x2018d
addl %ebp, %ecx
movsd 0x88(%rax), %xmm0
subsd 0x80(%rax), %xmm0
xorpd 0x16ec0(%rip), %xmm0    # 0x37040
movl %ecx, %ebp
jmp 0x201aa
movl 0x1a0(%rsp), %ecx
jmp 0x201ba
leal (%rcx,%rcx), %esi
xorpd %xmm0, %xmm0
cmpl %esi, %edx
jl 0x201aa
subl %ecx, %ebp
movsd 0x88(%rax), %xmm0
subsd 0x80(%rax), %xmm0
movl 0x1a0(%rsp), %ecx
movsd %xmm0, 0x98(%rsp)
cmpb $0x1, 0xaa(%rax)
jne 0x20215
leal (%r14,%rcx), %edx
movl 0x8(%rax), %ecx
cmpl %ecx, %edx
jge 0x201ee
addl %r14d, %ecx
movsd 0x98(%rax), %xmm0
subsd 0x90(%rax), %xmm0
xorpd 0x16e57(%rip), %xmm0    # 0x37040
movl %ecx, %r14d
jmp 0x2020c
leal (%rcx,%rcx), %esi
xorpd %xmm0, %xmm0
cmpl %esi, %edx
jl 0x2020c
subl %ecx, %r14d
movsd 0x98(%rax), %xmm0
subsd 0x90(%rax), %xmm0
movsd %xmm0, 0x90(%rsp)
addl 0x6c(%rsp), %ebx
imull 0x4(%rax), %r14d
addl %ebp, %r14d
imull (%rax), %r14d
addl %ebx, %r14d
movq 0x60(%r12), %rcx
movslq %r14d, %rbx
cmpl $0x0, (%rcx,%rbx,4)
jle 0x204dd
movsd 0x60(%rsp), %xmm0
subsd 0xa0(%rsp), %xmm0
movsd %xmm0, 0x48(%rsp)
subsd 0x98(%rsp), %xmm1
movsd %xmm1, 0xc0(%rsp)
movsd 0x58(%rsp), %xmm0
subsd 0x90(%rsp), %xmm0
movsd %xmm0, 0x20(%rsp)
movsd 0x128(%rsp), %xmm0
movsd 0xf0(%rax), %xmm1
addsd %xmm0, %xmm1
mulsd %xmm6, %xmm0
ucomisd 0x18809(%rip), %xmm0    # 0x38aa0
jb 0x2029f
sqrtsd %xmm0, %xmm0
jmp 0x202b6
movsd %xmm1, 0x40(%rsp)
callq 0x20f0
movsd 0x40(%rsp), %xmm1
movsd 0x8(%rsp), %xmm6
movl $0x18, %r14d
xorl %ebp, %ebp
ucomisd %xmm0, %xmm1
jbe 0x2041f
movq (%r12), %rax
movq 0x58(%r12), %rcx
movq (%rcx,%rbx,8), %rcx
movslq 0x4c(%r12), %rdx
movslq %ebp, %rsi
imulq %rdx, %rsi
movsd (%rcx,%rsi,8), %xmm4
movsd 0x8(%rcx,%rsi,8), %xmm1
subsd 0x48(%rsp), %xmm4
subsd 0xc0(%rsp), %xmm1
movsd 0x10(%rcx,%rsi,8), %xmm2
subsd 0x20(%rsp), %xmm2
movapd %xmm1, %xmm0
mulsd %xmm1, %xmm0
movapd %xmm4, %xmm5
mulsd %xmm4, %xmm5
addsd %xmm0, %xmm5
movapd %xmm2, %xmm3
mulsd %xmm2, %xmm3
addsd %xmm5, %xmm3
movsd 0xe8(%rax), %xmm5
movq 0xd8(%rax), %rax
movq (%rax,%rbx,8), %rax
movsd (%rax,%r14), %xmm7
mulsd %xmm3, %xmm6
ucomisd 0x18754(%rip), %xmm6    # 0x38aa0
jb 0x20357
xorps %xmm0, %xmm0
sqrtsd %xmm6, %xmm0
jmp 0x203c6
movapd %xmm6, %xmm0
movsd %xmm4, 0x40(%rsp)
movsd %xmm1, 0x120(%rsp)
movsd %xmm2, 0x118(%rsp)
movsd %xmm3, 0x110(%rsp)
movsd %xmm5, 0x108(%rsp)
movsd %xmm7, 0x100(%rsp)
callq 0x20f0
movsd 0x100(%rsp), %xmm7
movsd 0x108(%rsp), %xmm5
movsd 0x110(%rsp), %xmm3
movsd 0x118(%rsp), %xmm2
movsd 0x120(%rsp), %xmm1
movsd 0x40(%rsp), %xmm4
mulsd %xmm7, %xmm7
subsd %xmm7, %xmm5
addsd %xmm5, %xmm3
ucomisd %xmm3, %xmm0
jbe 0x203fb
movq 0x50(%r12), %rax
movq (%rax,%rbx,8), %rax
movl (%rax,%rbp,4), %edx
movq %r15, %rdi
movq %r15, %rsi
movapd %xmm4, %xmm0
callq 0xe636
testb %al, %al
je 0x1fd1d
incq %rbp
movq 0x60(%r12), %rax
movslq (%rax,%rbx,4), %rax
addq $0x20, %r14
cmpq %rax, %rbp
movsd 0x8(%rsp), %xmm6
jl 0x202c8
jmp 0x204dd
movq 0x58(%r12), %rax
movq (%rax,%rbx,8), %rax
movslq 0x4c(%r12), %rcx
movslq %ebp, %rdx
imulq %rcx, %rdx
movsd (%rax,%rdx,8), %xmm0
movsd 0x8(%rax,%rdx,8), %xmm1
subsd 0x48(%rsp), %xmm0
subsd 0xc0(%rsp), %xmm1
movsd 0x10(%rax,%rdx,8), %xmm2
subsd 0x20(%rsp), %xmm2
movapd %xmm1, %xmm3
mulsd %xmm1, %xmm3
movapd %xmm0, %xmm4
mulsd %xmm0, %xmm4
addsd %xmm3, %xmm4
movapd %xmm2, %xmm3
mulsd %xmm2, %xmm3
addsd %xmm4, %xmm3
movq (%r12), %rax
addsd 0xe8(%rax), %xmm3
movq 0x50(%r12), %rcx
movq 0xd8(%rax), %rax
movq (%rax,%rbx,8), %rax
movsd (%rax,%r14), %xmm4
mulsd %xmm4, %xmm4
subsd %xmm4, %xmm3
movq (%rcx,%rbx,8), %rax
movl (%rax,%rbp,4), %edx
movq %r15, %rdi
movq %r15, %rsi
callq 0xe636
testb %al, %al
je 0x1fd1d
incq %rbp
movq 0x60(%r12), %rax
movslq (%rax,%rbx,4), %rax
addq $0x20, %r14
cmpq %rax, %rbp
movsd 0x8(%rsp), %xmm6
jl 0x2041f
leaq 0x1(%r13), %rcx
movq 0x50(%rsp), %r14
leaq 0x1(%r14), %rax
cmpq 0x70(%rsp), %rcx
jne 0x1ff7f
incl 0x70(%r12)
jne 0x20530
cmpl $0x0, 0x48(%r12)
jle 0x20527
movq 0x88(%r12), %rax
movq %rax, %rcx
movl $0x0, (%rcx)
addq $0x4, %rcx
movslq 0x48(%r12), %rdx
leaq (%rax,%rdx,4), %rdx
cmpq %rdx, %rcx
jb 0x2050f
movl $0x1, 0x70(%r12)
movq 0x90(%r12), %rax
movq %rax, 0x88(%rsp)
movq %rax, 0x38(%rsp)
cmpq $0x3d, %r13
movl 0x10(%rsp), %ebx
ja 0x20ae1
movl 0x80(%rsp), %eax
cmpq %rax, %r14
jne 0x20591
movq %r15, %rdi
callq 0x9a14
movapd %xmm0, %xmm6
movq 0xd8(%rsp), %rax
cmpq 0x24bee(%rip), %rax     # 0x45168
je 0x20591
movl (%rax), %ecx
movl %ecx, 0x80(%rsp)
addq $0x4, %rax
movq %rax, 0xd8(%rsp)
movq (%r12), %rax
movq 0x78(%rsp), %rcx
movsd (%rcx,%r14,8), %xmm0
movsd 0xf0(%rax), %xmm1
addsd %xmm0, %xmm1
mulsd %xmm6, %xmm0
ucomisd 0x184e8(%rip), %xmm0    # 0x38aa0
jb 0x205c0
sqrtsd %xmm0, %xmm0
jmp 0x205dd
movsd %xmm6, 0x8(%rsp)
movsd %xmm1, 0x50(%rsp)
callq 0x20f0
movsd 0x50(%rsp), %xmm1
movsd 0x8(%rsp), %xmm6
ucomisd %xmm0, %xmm1
ja 0x20acf
movq 0x150(%rsp), %rax
movl 0x4(%rax,%r14,4), %ebp
xorl 0x18(%rsp), %ebp
addl 0x28(%rsp), %ebp
movl %ebp, %esi
andl $0x7f, %esi
addl $-0x40, %esi
movl %esi, %r13d
addl 0x14(%rsp), %r13d
js 0x20aa9
movl 0x38(%r12), %eax
cmpl %eax, %r13d
jge 0x20aa9
movl %ebp, %edx
shrl $0x7, %edx
andl $0x7f, %edx
addl $-0x40, %edx
movl %edx, %r9d
addl %ebx, %r9d
js 0x20aa9
movl 0x3c(%r12), %edi
cmpl %edi, %r9d
jge 0x20aa9
movl %ebp, %ecx
shrl $0xe, %ecx
andl $0x7f, %ecx
addl $-0x40, %ecx
movl %ecx, %r10d
addl 0x4(%rsp), %r10d
js 0x20aa9
cmpl 0x40(%r12), %r10d
jge 0x20aa9
movl %r13d, %r8d
shlq $0x2, %r8
addq 0x88(%r12), %r8
imull %r10d, %edi
addl %r9d, %edi
imull %eax, %edi
movq %rdi, 0x160(%rsp)
movl 0x70(%r12), %eax
movq %r8, 0x158(%rsp)
movl %eax, (%r8,%rdi,4)
movq %r12, %rdi
movapd 0x130(%rsp), %xmm0
movapd 0xf0(%rsp), %xmm1
movapd 0xe0(%rsp), %xmm2
movsd 0xb0(%rsp), %xmm3
movsd 0xa8(%rsp), %xmm4
movsd 0x148(%rsp), %xmm5
leaq 0x128(%rsp), %r8
movsd %xmm6, 0x8(%rsp)
movsd 0x8(%rsp), %xmm6
movl %r9d, 0xcc(%rsp)
movq %r10, 0x140(%rsp)
callq 0x1d98a
movq 0x140(%rsp), %r10
movl 0xcc(%rsp), %r9d
movsd 0x8(%rsp), %xmm6
testb %al, %al
jne 0x20aa9
movq %r14, 0x50(%rsp)
movq (%r12), %rdi
subq $0x8, %rsp
movq 0xd8(%rsp), %rsi
movq 0xc0(%rsp), %rdx
movl 0x1a8(%rsp), %ecx
movl %r13d, %r8d
leaq 0x74(%rsp), %rax
pushq %rax
leaq 0xa0(%rsp), %rax
pushq %rax
leaq 0xb0(%rsp), %rax
pushq %rax
leaq 0xc0(%rsp), %rax
pushq %rax
pushq %r10
callq 0x184c4
addq $0x30, %rsp
movq 0x60(%r12), %rcx
movslq %eax, %rdi
cmpl $0x0, (%rcx,%rdi,4)
jle 0x20a48
movsd 0x60(%rsp), %xmm0
subsd 0xa0(%rsp), %xmm0
movsd %xmm0, 0x70(%rsp)
movsd 0x30(%rsp), %xmm0
subsd 0x98(%rsp), %xmm0
movsd %xmm0, 0x48(%rsp)
movsd 0x58(%rsp), %xmm0
subsd 0x90(%rsp), %xmm0
movsd %xmm0, 0xc0(%rsp)
movq (%r12), %rax
movsd 0x128(%rsp), %xmm0
movsd 0xf0(%rax), %xmm1
addsd %xmm0, %xmm1
mulsd 0x8(%rsp), %xmm0
ucomisd 0x182b7(%rip), %xmm0    # 0x38aa0
movq %rdi, 0x20(%rsp)
jb 0x207f6
sqrtsd %xmm0, %xmm0
jmp 0x2080c
movsd %xmm1, 0x40(%rsp)
callq 0x20f0
movsd 0x40(%rsp), %xmm1
movq 0x20(%rsp), %rdi
ucomisd %xmm0, %xmm1
jbe 0x20978
movl $0x18, %r14d
xorl %ebx, %ebx
movq (%r12), %rax
movq 0x58(%r12), %rcx
movq (%rcx,%rdi,8), %rcx
movslq 0x4c(%r12), %rdx
movslq %ebx, %rsi
imulq %rdx, %rsi
movsd (%rcx,%rsi,8), %xmm4
movsd 0x8(%rcx,%rsi,8), %xmm1
subsd 0x70(%rsp), %xmm4
subsd 0x48(%rsp), %xmm1
movsd 0x10(%rcx,%rsi,8), %xmm2
subsd 0xc0(%rsp), %xmm2
movapd %xmm1, %xmm0
mulsd %xmm1, %xmm0
movapd %xmm4, %xmm5
mulsd %xmm4, %xmm5
addsd %xmm0, %xmm5
movapd %xmm2, %xmm3
mulsd %xmm2, %xmm3
addsd %xmm5, %xmm3
movsd 0xe8(%rax), %xmm5
movq 0xd8(%rax), %rax
movq (%rax,%rdi,8), %rax
movsd (%rax,%r14), %xmm6
movsd 0x8(%rsp), %xmm0
mulsd %xmm3, %xmm0
ucomisd 0x181f8(%rip), %xmm0    # 0x38aa0
jb 0x208b0
sqrtsd %xmm0, %xmm0
jmp 0x20920
movsd %xmm4, 0x40(%rsp)
movsd %xmm1, 0x120(%rsp)
movsd %xmm2, 0x118(%rsp)
movsd %xmm3, 0x110(%rsp)
movsd %xmm5, 0x108(%rsp)
movsd %xmm6, 0x100(%rsp)
callq 0x20f0
movsd 0x100(%rsp), %xmm6
movsd 0x108(%rsp), %xmm5
movsd 0x110(%rsp), %xmm3
movsd 0x118(%rsp), %xmm2
movsd 0x120(%rsp), %xmm1
movsd 0x40(%rsp), %xmm4
movq 0x20(%rsp), %rdi
mulsd %xmm6, %xmm6
subsd %xmm6, %xmm5
addsd %xmm5, %xmm3
ucomisd %xmm3, %xmm0
jbe 0x2095a
movq 0x50(%r12), %rax
movq (%rax,%rdi,8), %rax
movl (%rax,%rbx,4), %edx
movq %r15, %rdi
movq %r15, %rsi
movapd %xmm4, %xmm0
callq 0xe636
movq 0x20(%rsp), %rdi
testb %al, %al
je 0x1fd1d
incq %rbx
movq 0x60(%r12), %rax
movslq (%rax,%rdi,4), %rax
addq $0x20, %r14
cmpq %rax, %rbx
jl 0x2081e
jmp 0x20a48
movl $0x18, %esi
xorl %r14d, %r14d
movq 0x58(%r12), %rax
movq (%rax,%rdi,8), %rax
movslq 0x4c(%r12), %rcx
movslq %r14d, %rdx
imulq %rcx, %rdx
movsd (%rax,%rdx,8), %xmm0
movsd 0x8(%rax,%rdx,8), %xmm1
subsd 0x70(%rsp), %xmm0
subsd 0x48(%rsp), %xmm1
movsd 0x10(%rax,%rdx,8), %xmm2
subsd 0xc0(%rsp), %xmm2
movapd %xmm1, %xmm3
mulsd %xmm1, %xmm3
movapd %xmm0, %xmm4
mulsd %xmm0, %xmm4
addsd %xmm3, %xmm4
movapd %xmm2, %xmm3
mulsd %xmm2, %xmm3
addsd %xmm4, %xmm3
movq (%r12), %rax
addsd 0xe8(%rax), %xmm3
movq 0x50(%r12), %rcx
movq 0xd8(%rax), %rax
movq (%rax,%rdi,8), %rax
movq %rsi, 0x20(%rsp)
movsd (%rax,%rsi), %xmm4
mulsd %xmm4, %xmm4
subsd %xmm4, %xmm3
movq (%rcx,%rdi,8), %rax
movl (%rax,%r14,4), %edx
movq %rdi, %rbx
movq %r15, %rdi
movq %r15, %rsi
callq 0xe636
testb %al, %al
je 0x1fd1d
incq %r14
movq 0x60(%r12), %rax
movslq (%rax,%rbx,4), %rax
movq 0x20(%rsp), %rsi
addq $0x20, %rsi
cmpq %rax, %r14
movq %rbx, %rdi
jl 0x20980
movq 0x98(%r12), %rax
addq $-0x48, %rax
cmpq %rax, 0x38(%rsp)
ja 0x20ab8
movq 0x160(%rsp), %rax
movq 0x158(%rsp), %rcx
leaq (%rcx,%rax,4), %rdx
subq $0x8, %rsp
movq %r12, %rdi
movl %ebp, %esi
movl %r13d, %ecx
movl 0xd4(%rsp), %r8d
movq 0x148(%rsp), %r9
leaq 0x40(%rsp), %rax
pushq %rax
callq 0x1ec3a
addq $0x10, %rsp
movl 0x10(%rsp), %ebx
movsd 0x8(%rsp), %xmm6
movq 0x50(%rsp), %r14
incq %r14
cmpq $0x3f, %r14
jne 0x20553
jmp 0x20ad7
movq %r12, %rdi
leaq 0x88(%rsp), %rsi
leaq 0x38(%rsp), %rdx
callq 0x2e82
jmp 0x20a5b
movb $0x1, %bpl
jmp 0x1fd1f
movl $0x3f, %eax
movq %rax, 0x70(%rsp)
movq (%r12), %rax
movl 0x70(%rsp), %ecx
movq 0x78(%rsp), %rdx
movsd (%rdx,%rcx,8), %xmm0
movsd 0xf0(%rax), %xmm1
addsd %xmm0, %xmm1
mulsd %xmm0, %xmm6
xorpd %xmm0, %xmm0
ucomisd %xmm0, %xmm6
jb 0x20b16
xorps %xmm0, %xmm0
sqrtsd %xmm6, %xmm0
jmp 0x20b2b
movapd %xmm6, %xmm0
movsd %xmm1, 0x8(%rsp)
callq 0x20f0
movsd 0x8(%rsp), %xmm1
movb $0x1, %bpl
ucomisd %xmm0, %xmm1
ja 0x1fd1f
movq 0x88(%rsp), %rax
cmpq 0x38(%rsp), %rax
je 0x1fd1f
movapd 0x164ed(%rip), %xmm0    # 0x37040
movapd 0x130(%rsp), %xmm1
xorpd %xmm0, %xmm1
movapd %xmm1, 0x130(%rsp)
movapd 0xf0(%rsp), %xmm1
xorpd %xmm0, %xmm1
movapd %xmm1, 0xf0(%rsp)
movapd 0xe0(%rsp), %xmm1
xorpd %xmm0, %xmm1
movapd %xmm1, 0xe0(%rsp)
cmpq 0x98(%r12), %rax
jne 0x20ba7
movq 0x90(%r12), %rax
movl (%rax), %r14d
movl 0x4(%rax), %r13d
leaq 0xc(%rax), %rcx
movq %rcx, 0x88(%rsp)
movl 0x8(%rax), %esi
movl %r13d, %eax
subl %ebx, %eax
movl %esi, %ecx
subl 0x4(%rsp), %ecx
movl %r14d, %edx
movsd 0x10(%r12), %xmm4
xorps %xmm7, %xmm7
cvtsi2sd %eax, %xmm7
mulsd %xmm4, %xmm7
addsd 0xf0(%rsp), %xmm7
xorps %xmm6, %xmm6
cvtsi2sd %ecx, %xmm6
addsd %xmm7, %xmm4
movsd 0x18(%r12), %xmm5
mulsd %xmm5, %xmm6
addsd 0xe0(%rsp), %xmm6
addsd %xmm6, %xmm5
subl 0x14(%rsp), %edx
xorps %xmm8, %xmm8
cvtsi2sd %edx, %xmm8
movsd 0x8(%r12), %xmm3
mulsd %xmm3, %xmm8
addsd 0x130(%rsp), %xmm8
addsd %xmm8, %xmm3
movq %rsi, 0x18(%rsp)
jle 0x20c5f
cmpl %ebx, %r13d
jle 0x20c90
cmpl 0x4(%rsp), %esi
jle 0x20ce0
movq %r12, %rdi
movq %r15, %rsi
movapd %xmm8, %xmm0
movapd %xmm7, %xmm1
movapd %xmm6, %xmm2
jmp 0x20e6e
jge 0x20cb5
cmpl %ebx, %r13d
jle 0x20cfe
cmpl 0x4(%rsp), %esi
jle 0x20d7c
movq %r12, %rdi
movq %r15, %rsi
movapd %xmm3, %xmm0
movapd %xmm7, %xmm1
movapd %xmm6, %xmm2
movapd %xmm8, %xmm3
jmp 0x20e6e
jge 0x20d2e
cmpl 0x4(%rsp), %esi
jle 0x20d9e
movq %r12, %rdi
movq %r15, %rsi
movapd %xmm8, %xmm0
movapd %xmm4, %xmm1
movapd %xmm6, %xmm2
jmp 0x20d25
cmpl %ebx, %r13d
jle 0x20d50
cmpl 0x4(%rsp), %esi
jle 0x20dbc
movq %r12, %rdi
movq %r15, %rsi
movapd %xmm8, %xmm0
movapd %xmm7, %xmm1
movapd %xmm6, %xmm2
jmp 0x20e95
movq %r12, %rdi
movq %r15, %rsi
movapd %xmm8, %xmm0
movapd %xmm7, %xmm1
jge 0x20dda
movapd %xmm5, %xmm2
jmp 0x20e6a
jge 0x20de3
cmpl 0x4(%rsp), %esi
jle 0x20e49
movq %r12, %rdi
movq %r15, %rsi
movapd %xmm3, %xmm0
movapd %xmm4, %xmm1
movapd %xmm6, %xmm2
movapd %xmm8, %xmm3
movapd %xmm7, %xmm4
jmp 0x20e6e
cmpl 0x4(%rsp), %esi
jle 0x20e09
movq %r12, %rdi
movq %r15, %rsi
movapd %xmm8, %xmm0
movapd %xmm7, %xmm1
movapd %xmm6, %xmm2
jmp 0x20edc
jge 0x20e27
cmpl 0x4(%rsp), %esi
jle 0x20e78
movq %r12, %rdi
movq %r15, %rsi
movapd %xmm8, %xmm0
movapd %xmm4, %xmm1
movapd %xmm6, %xmm2
movapd %xmm7, %xmm4
jmp 0x20e95
movq %r12, %rdi
movq %r15, %rsi
movapd %xmm3, %xmm0
movapd %xmm7, %xmm1
jge 0x20eae
movapd %xmm5, %xmm2
movapd %xmm8, %xmm3
jmp 0x20e6a
movq %r12, %rdi
movq %r15, %rsi
movapd %xmm8, %xmm0
movapd %xmm4, %xmm1
jge 0x20eb9
movapd %xmm5, %xmm2
jmp 0x20e66
movq %r12, %rdi
movq %r15, %rsi
movapd %xmm8, %xmm0
movapd %xmm7, %xmm1
jge 0x20e9f
movapd %xmm5, %xmm2
jmp 0x20e91
movapd %xmm6, %xmm2
jmp 0x20f27
cmpl 0x4(%rsp), %esi
jle 0x20ebf
movq %r12, %rdi
movq %r15, %rsi
movapd %xmm3, %xmm0
movapd %xmm7, %xmm1
movapd %xmm6, %xmm2
movapd %xmm8, %xmm3
jmp 0x20edc
movq %r12, %rdi
movq %r15, %rsi
movapd %xmm8, %xmm0
jge 0x20ee3
movapd %xmm7, %xmm1
movapd %xmm5, %xmm2
jmp 0x20ed8
cmpl 0x4(%rsp), %esi
jle 0x20efa
movq %r12, %rdi
movq %r15, %rsi
movapd %xmm8, %xmm0
movapd %xmm7, %xmm1
movapd %xmm6, %xmm2
jmp 0x20f13
movq %r12, %rdi
movq %r15, %rsi
movapd %xmm3, %xmm0
movapd %xmm4, %xmm1
jge 0x20f1a
movapd %xmm5, %xmm2
movapd %xmm8, %xmm3
movapd %xmm7, %xmm4
movapd %xmm6, %xmm5
callq 0x21254
jmp 0x20f2c
movq %r12, %rdi
movq %r15, %rsi
movapd %xmm8, %xmm0
movapd %xmm4, %xmm1
jge 0x20e9f
movapd %xmm5, %xmm2
movapd %xmm7, %xmm4
movapd %xmm6, %xmm5
callq 0x218ac
jmp 0x20f2c
movapd %xmm6, %xmm2
movapd %xmm5, %xmm4
callq 0x21a50
jmp 0x20f2c
movapd %xmm6, %xmm2
movapd %xmm8, %xmm3
jmp 0x20f27
movapd %xmm6, %xmm2
jmp 0x20f23
movq %r12, %rdi
movq %r15, %rsi
movapd %xmm3, %xmm0
jge 0x20ee3
movapd %xmm7, %xmm1
movapd %xmm5, %xmm2
movapd %xmm8, %xmm3
movapd %xmm6, %xmm5
callq 0x21608
jmp 0x20f2c
movapd %xmm7, %xmm1
movapd %xmm6, %xmm2
movapd %xmm4, %xmm3
movapd %xmm5, %xmm4
callq 0x217ac
jmp 0x20f2c
jge 0x2123d
movq %r12, %rdi
movq %r15, %rsi
movapd %xmm8, %xmm0
movapd %xmm7, %xmm1
movapd %xmm5, %xmm2
callq 0x21b50
jmp 0x20f2c
movapd %xmm6, %xmm2
movapd %xmm8, %xmm3
movapd %xmm7, %xmm4
callq 0x21460
testb %al, %al
jne 0x21209
movq (%r12), %rax
cmpb $0x1, 0xa8(%rax)
movl %r14d, %ecx
movq 0x18(%rsp), %r10
jne 0x20f9a
movq 0xd0(%rsp), %rcx
leal (%r14,%rcx), %esi
movl (%rax), %edx
cmpl %edx, %esi
jge 0x20f74
addl %r14d, %edx
movsd 0x78(%rax), %xmm0
subsd 0x70(%rax), %xmm0
xorpd 0x160d0(%rip), %xmm0    # 0x37040
movl %edx, %ecx
jmp 0x20f91
leal (%rdx,%rdx), %edi
xorpd %xmm0, %xmm0
movl %r14d, %ecx
cmpl %edi, %esi
jl 0x20f91
movl %r14d, %ecx
subl %edx, %ecx
movsd 0x78(%rax), %xmm0
subsd 0x70(%rax), %xmm0
movsd %xmm0, 0xa0(%rsp)
cmpb $0x1, 0xa9(%rax)
movl %r13d, %edx
jne 0x21006
movq 0xb8(%rsp), %rdx
leal (%rdx,%r13), %edi
movl 0x4(%rax), %esi
cmpl %esi, %edi
jge 0x20fd8
addl %r13d, %esi
movsd 0x88(%rax), %xmm0
subsd 0x80(%rax), %xmm0
xorpd 0x1606c(%rip), %xmm0    # 0x37040
movl %esi, %edx
jmp 0x20ffd
leal (%rsi,%rsi), %r8d
xorpd %xmm0, %xmm0
movl %r13d, %edx
cmpl %r8d, %edi
jl 0x20ffd
movl %r13d, %edx
subl %esi, %edx
movsd 0x88(%rax), %xmm0
subsd 0x80(%rax), %xmm0
movsd %xmm0, 0x98(%rsp)
cmpb $0x1, 0xaa(%rax)
movl %r10d, %esi
jne 0x21072
movl 0x1a0(%rsp), %esi
leal (%r10,%rsi), %r8d
movl 0x8(%rax), %edi
cmpl %edi, %r8d
jge 0x21044
addl %r10d, %edi
movsd 0x98(%rax), %xmm0
subsd 0x90(%rax), %xmm0
xorpd 0x16000(%rip), %xmm0    # 0x37040
movl %edi, %esi
jmp 0x21069
leal (%rdi,%rdi), %r9d
xorpd %xmm0, %xmm0
movl %r10d, %esi
cmpl %r9d, %r8d
jl 0x21069
movl %r10d, %esi
subl %edi, %esi
movsd 0x98(%rax), %xmm0
subsd 0x90(%rax), %xmm0
movsd %xmm0, 0x90(%rsp)
movq %r13, 0x78(%rsp)
movq %r14, 0x28(%rsp)
addl 0x6c(%rsp), %ecx
imull 0x4(%rax), %esi
addl %edx, %esi
imull (%rax), %esi
addl %ecx, %esi
movq 0x60(%r12), %rax
movslq %esi, %rbx
cmpl $0x0, (%rax,%rbx,4)
jle 0x211a1
movsd 0x60(%rsp), %xmm0
subsd 0xa0(%rsp), %xmm0
movsd %xmm0, 0x8(%rsp)
movsd 0x30(%rsp), %xmm0
subsd 0x98(%rsp), %xmm0
movsd %xmm0, 0x50(%rsp)
movsd 0x58(%rsp), %xmm0
subsd 0x90(%rsp), %xmm0
movsd %xmm0, 0x80(%rsp)
movl $0x18, %r14d
xorl %r13d, %r13d
movq 0x58(%r12), %rax
movq (%rax,%rbx,8), %rax
movslq 0x4c(%r12), %rcx
movslq %r13d, %rdx
imulq %rcx, %rdx
movsd (%rax,%rdx,8), %xmm0
subsd 0x8(%rsp), %xmm0
movsd 0x8(%rax,%rdx,8), %xmm1
subsd 0x50(%rsp), %xmm1
movsd 0x10(%rax,%rdx,8), %xmm2
subsd 0x80(%rsp), %xmm2
movapd %xmm1, %xmm3
mulsd %xmm1, %xmm3
movapd %xmm0, %xmm4
mulsd %xmm0, %xmm4
addsd %xmm3, %xmm4
movapd %xmm2, %xmm3
mulsd %xmm2, %xmm3
addsd %xmm4, %xmm3
movq (%r12), %rax
addsd 0xe8(%rax), %xmm3
movq 0x50(%r12), %rcx
movq 0xd8(%rax), %rax
movq (%rax,%rbx,8), %rax
movsd (%rax,%r14), %xmm4
mulsd %xmm4, %xmm4
subsd %xmm4, %xmm3
movq (%rcx,%rbx,8), %rax
movl (%rax,%r13,4), %edx
movq %r15, %rdi
movq %r15, %rsi
callq 0xe636
testb %al, %al
je 0x1fd1d
incq %r13
movq 0x60(%r12), %rax
movslq (%rax,%rbx,4), %rax
addq $0x20, %r14
cmpq %rax, %r13
jl 0x210e8
movq 0x88(%rsp), %rax
movq 0x38(%rsp), %rcx
movq %rax, %rdx
subq %rcx, %rdx
jbe 0x211bf
sarq $0x2, %rdx
movq %rdx, %rax
jmp 0x211dd
movq 0x98(%r12), %rdx
subq %rcx, %rdx
sarq $0x2, %rdx
subq 0x90(%r12), %rax
sarq $0x2, %rax
addq %rdx, %rax
movl 0x10(%rsp), %ebx
movq 0x28(%rsp), %r14
movq 0x78(%rsp), %r13
movq 0x18(%rsp), %rcx
cmpq $0x11, %rax
jle 0x21221
movq %r12, %rdi
movl %r14d, %esi
movl %r13d, %edx
leaq 0x38(%rsp), %r8
callq 0x1f78e
movq 0x88(%rsp), %rax
cmpq 0x38(%rsp), %rax
jne 0x20b95
jmp 0x1fd1f
movq %r12, %rdi
leaq 0x88(%rsp), %rsi
leaq 0x38(%rsp), %rdx
callq 0x2e82
movq 0x18(%rsp), %rcx
jmp 0x211f6
leaq 0x1811c(%rip), %rdi     # 0x39360
movl $0x3, %esi
callq 0x133c5
jmp 0x20f34
nop