// File: /workspace/github_cpp_makefile_star2_O3/jamesliu516[P]amrdg2d/hllc.cpp
// Function: HLLC_flux(conservedvariables&, conservedvariables&, pointxy&, double*)
void HLLC_flux(SHBL &wl, SHBL &wr, PXYZ &normal, double flxface[4])
//sface ��������
{
	double s_l, s_r, s_m;
	double pStar;
	double aFace, ql, qr, al, ar, qFace; // speed

	double omgL, omgR;
//
	JBBL ul, ur, uFace;
	SHBL wLstar, wRstar;
	
	shblToJbbl(&wl, &ul);

	al = sqrt(GAMMA*ul.p/ul.q);

	ql = normal.x * ul.u + normal.y * ul.v;

	shblToJbbl(&wr, &ur);

	ar = sqrt(GAMMA*ur.p/ur.q);
	qr = normal.x * ur.u + normal.y * ur.v;
//
//
	RoeRiemannSolver(&ul, &ur, &uFace, &aFace);
	qFace = normal.x * uFace.u + normal.y * uFace.v;


	s_l = min2(ql-al, qFace-aFace);
	s_r = max2(qr+ar, qFace+aFace);


	s_m = (ur.q * qr * (s_r - qr) - ul.q * ql * (s_l - ql) + ul.p - ur.p)
			/ (ur.q * (s_r - qr) - ul.q * (s_l - ql));

	omgL = 1.0 / (s_l-s_m);
	omgR = 1.0 / (s_r-s_m);

	pStar = ul.q * (ql - s_l) * (ql - s_m) + ul.p;

	wLstar.q = omgL * ul.q * (s_l - ql);
	wLstar.qu = omgL * ((s_l - ql) * wl.qu + (pStar - ul.p) * normal.x);
	wLstar.qv = omgL * ((s_l - ql) * wl.qv + (pStar - ul.p) * normal.y);
	wLstar.te = omgL * ((s_l - ql) * wl.te - ul.p * ql + pStar * s_m);

	wRstar.q = omgR * ur.q * (s_r - qr);
	wRstar.qu = omgR * ((s_r - qr) * wr.qu + (pStar - ur.p) * normal.x);
	wRstar.qv = omgR * ((s_r - qr) * wr.qv + (pStar - ur.p) * normal.y);
	wRstar.te = omgR * ((s_r - qr) * wr.te - ur.p * qr + pStar * s_m);

	if(s_l > 0.0)
	{
		flxface[0] = ul.q * ql;
		flxface[1] = wl.qu * ql + ul.p * normal.x;
		flxface[2] = wl.qv * ql + ul.p * normal.y;
		flxface[3] = (wl.te + ul.p) * ql;

	}
	else if(s_r < 0.0)
	{
		flxface[0] = ur.q * qr;
		flxface[1] = wr.qu * qr + ur.p * normal.x;
		flxface[2] = wr.qv * qr + ur.p * normal.y;
		flxface[3] = (wr.te + ur.p) * qr;
	}
	else if(s_l <= 0.0 && s_m > 0.0)
	{
		flxface[0] = wLstar.q * s_m;
		flxface[1] = wLstar.qu * s_m + pStar * normal.x;
		flxface[2] = wLstar.qv * s_m + pStar * normal.y;
		flxface[3] = (wLstar.te + pStar) * s_m;
	}
	else if(s_m <= 0.0 && s_r >= 0.0)
	{
		flxface[0] = wRstar.q * s_m;
		flxface[1] = wRstar.qu * s_m + pStar * normal.x;
		flxface[2] = wRstar.qv * s_m + pStar * normal.y;
		flxface[3] = (wRstar.te + pStar) * s_m;
	}
}

// Assembly:
pushq %rbp
pushq %r15
pushq %r14
pushq %r13
pushq %r12
pushq %rbx
subq $0xc8, %rsp
movq %rdx, %r14
movq %rdi, %r15
xorps %xmm0, %xmm0
movaps %xmm0, 0xb0(%rsp)
movaps %xmm0, 0xa0(%rsp)
movsd (%rdi), %xmm1
movsd %xmm1, 0x78(%rsp)
movupd 0x8(%rdi), %xmm5
movapd %xmm1, %xmm0
unpcklpd %xmm1, %xmm0            # xmm0 = xmm0[0],xmm1[0]
divpd %xmm0, %xmm5
movupd %xmm5, 0x80(%rsp)
movapd %xmm5, %xmm3
unpckhpd %xmm5, %xmm3            # xmm3 = xmm3[1],xmm5[1]
movapd %xmm3, %xmm0
mulsd %xmm3, %xmm0
movapd %xmm5, %xmm2
mulsd %xmm5, %xmm2
addsd %xmm0, %xmm2
movsd 0xdf0b(%rip), %xmm0     # 0x1d140
mulsd %xmm1, %xmm0
mulsd %xmm2, %xmm0
addsd 0x18(%rdi), %xmm0
mulsd 0xdf5e(%rip), %xmm0     # 0x1d1a8
movq %rsi, %r12
movsd %xmm0, 0x90(%rsp)
mulsd 0xdf52(%rip), %xmm0     # 0x1d1b0
divsd %xmm1, %xmm0
xorpd %xmm6, %xmm6
ucomisd %xmm6, %xmm0
movq %rcx, 0x28(%rsp)
jb 0xf277
sqrtsd %xmm0, %xmm0
jmp 0xf296
movapd %xmm5, 0x30(%rsp)
movapd %xmm3, (%rsp)
callq 0x31b0
movapd (%rsp), %xmm3
xorpd %xmm6, %xmm6
movapd 0x30(%rsp), %xmm5
movsd %xmm0, 0x20(%rsp)
movsd 0x8(%r14), %xmm1
mulsd %xmm3, %xmm1
movsd (%r14), %xmm4
movsd (%r12), %xmm2
movsd %xmm2, 0x58(%rsp)
movupd 0x8(%r12), %xmm7
movapd %xmm2, %xmm0
unpcklpd %xmm2, %xmm0            # xmm0 = xmm0[0],xmm2[0]
divpd %xmm0, %xmm7
movapd %xmm7, %xmm8
unpckhpd %xmm7, %xmm8            # xmm8 = xmm8[1],xmm7[1]
movapd %xmm8, %xmm0
mulsd %xmm8, %xmm0
movapd %xmm7, %xmm3
mulsd %xmm7, %xmm3
addsd %xmm0, %xmm3
movsd 0xde4e(%rip), %xmm0     # 0x1d140
mulsd %xmm2, %xmm0
mulsd %xmm3, %xmm0
movupd %xmm7, 0x60(%rsp)
addsd 0x18(%r12), %xmm0
mulsd 0xde99(%rip), %xmm0     # 0x1d1a8
mulsd %xmm5, %xmm4
movsd %xmm0, 0x70(%rsp)
mulsd 0xde8f(%rip), %xmm0     # 0x1d1b0
addsd %xmm1, %xmm4
movsd %xmm4, 0x30(%rsp)
divsd %xmm2, %xmm0
ucomisd %xmm6, %xmm0
jb 0xf33b
sqrtsd %xmm0, %xmm0
jmp 0xf35a
movapd %xmm7, 0x10(%rsp)
movapd %xmm8, 0x40(%rsp)
callq 0x31b0
movapd 0x40(%rsp), %xmm8
movapd 0x10(%rsp), %xmm7
movsd %xmm0, (%rsp)
movsd 0x8(%r14), %xmm0
mulsd %xmm8, %xmm0
movsd %xmm0, 0x40(%rsp)
movsd (%r14), %xmm0
mulsd %xmm7, %xmm0
movsd %xmm0, 0x10(%rsp)
leaq 0x78(%rsp), %r13
leaq 0x58(%rsp), %rsi
leaq 0xa0(%rsp), %rbp
leaq 0x98(%rsp), %rbx
movq %r13, %rdi
movq %rbp, %rdx
movq %rbx, %rcx
callq 0xfa80
movsd (%r14), %xmm5
movsd 0x8(%r14), %xmm1
movsd 0x10(%rbp), %xmm2
mulsd %xmm1, %xmm2
movsd 0x8(%rbp), %xmm0
mulsd %xmm5, %xmm0
addsd %xmm2, %xmm0
movsd 0x30(%rsp), %xmm13
movapd %xmm13, %xmm10
subsd 0x20(%rsp), %xmm10
movsd (%rbx), %xmm3
movapd %xmm0, %xmm2
subsd %xmm3, %xmm2
minsd %xmm2, %xmm10
movsd (%r13), %xmm6
movsd 0x18(%r13), %xmm4
movapd %xmm13, %xmm2
mulsd %xmm6, %xmm2
xorpd %xmm8, %xmm8
ucomisd %xmm8, %xmm10
jbe 0xf461
movq 0x28(%rsp), %rax
movsd %xmm2, (%rax)
movsd (%r14), %xmm0
mulsd %xmm4, %xmm0
movsd 0x8(%r15), %xmm1
mulsd %xmm13, %xmm1
addsd %xmm0, %xmm1
movsd %xmm1, 0x8(%rax)
movsd 0x8(%r14), %xmm0
mulsd %xmm4, %xmm0
movsd 0x10(%r15), %xmm1
mulsd %xmm13, %xmm1
addsd %xmm0, %xmm1
movsd %xmm1, 0x10(%rax)
addsd 0x18(%r15), %xmm4
mulsd %xmm13, %xmm4
jmp 0xf71d
movsd 0x10(%rsp), %xmm11
addsd 0x40(%rsp), %xmm11
movsd (%rsp), %xmm9
addsd %xmm11, %xmm9
addsd %xmm0, %xmm3
maxsd %xmm3, %xmm9
movsd 0x58(%rsp), %xmm12
movsd 0x70(%rsp), %xmm7
movapd %xmm11, %xmm0
mulsd %xmm12, %xmm0
ucomisd %xmm9, %xmm8
movq 0x28(%rsp), %rax
movapd %xmm11, %xmm3
jbe 0xf4fe
movsd %xmm0, (%rax)
movsd (%r14), %xmm0
mulsd %xmm7, %xmm0
movsd 0x8(%r12), %xmm1
mulsd %xmm3, %xmm1
addsd %xmm0, %xmm1
movsd %xmm1, 0x8(%rax)
movsd 0x8(%r14), %xmm0
mulsd %xmm7, %xmm0
movsd 0x10(%r12), %xmm1
mulsd %xmm3, %xmm1
addsd %xmm0, %xmm1
movsd %xmm1, 0x10(%rax)
addsd 0x18(%r12), %xmm7
mulsd %xmm3, %xmm7
movapd %xmm7, %xmm4
jmp 0xf71d
movsd %xmm9, (%rsp)
movsd %xmm3, 0x10(%rsp)
subsd %xmm3, %xmm9
movapd %xmm10, %xmm11
subsd %xmm13, %xmm11
mulsd %xmm11, %xmm2
mulsd %xmm9, %xmm0
subsd %xmm2, %xmm0
addsd %xmm4, %xmm0
subsd %xmm7, %xmm0
movapd %xmm11, %xmm2
mulsd %xmm6, %xmm2
movapd %xmm12, %xmm3
mulsd %xmm9, %xmm3
subsd %xmm2, %xmm3
divsd %xmm3, %xmm0
movapd %xmm13, %xmm3
subsd %xmm10, %xmm3
mulsd %xmm6, %xmm3
movapd %xmm13, %xmm2
subsd %xmm0, %xmm2
mulsd %xmm3, %xmm2
addsd %xmm4, %xmm2
movapd %xmm2, %xmm3
mulsd %xmm0, %xmm3
movsd 0x8(%r12), %xmm14
mulsd %xmm9, %xmm14
movsd 0x10(%r12), %xmm13
mulsd %xmm9, %xmm13
movsd 0x18(%r12), %xmm15
mulsd %xmm9, %xmm15
ucomisd %xmm10, %xmm8
jb 0xf65d
ucomisd %xmm8, %xmm0
jbe 0xf65d
subsd %xmm0, %xmm10
movsd 0xdbe7(%rip), %xmm9     # 0x1d1a0
divsd %xmm10, %xmm9
mulsd %xmm9, %xmm6
mulsd %xmm11, %xmm6
movapd %xmm2, %xmm8
subsd %xmm4, %xmm8
mulsd %xmm8, %xmm5
movsd 0x8(%r15), %xmm7
mulsd %xmm11, %xmm7
addsd %xmm5, %xmm7
mulsd %xmm9, %xmm7
mulsd %xmm8, %xmm1
movsd 0x10(%r15), %xmm5
mulsd %xmm11, %xmm5
addsd %xmm1, %xmm5
mulsd %xmm9, %xmm5
movsd 0x30(%rsp), %xmm1
mulsd %xmm4, %xmm1
mulsd 0x18(%r15), %xmm11
subsd %xmm1, %xmm11
addsd %xmm11, %xmm3
mulsd %xmm9, %xmm3
mulsd %xmm0, %xmm6
movsd %xmm6, (%rax)
movsd (%r14), %xmm1
mulsd %xmm2, %xmm1
mulsd %xmm0, %xmm7
addsd %xmm1, %xmm7
movsd %xmm7, 0x8(%rax)
movsd 0x8(%r14), %xmm1
mulsd %xmm2, %xmm1
mulsd %xmm0, %xmm5
addsd %xmm1, %xmm5
movsd %xmm5, 0x10(%rax)
jmp 0xf711
movsd (%rsp), %xmm6
ucomisd %xmm8, %xmm6
jb 0xf722
xorpd %xmm4, %xmm4
ucomisd %xmm0, %xmm4
jb 0xf722
subsd %xmm0, %xmm6
movsd 0xdb18(%rip), %xmm8     # 0x1d1a0
divsd %xmm6, %xmm8
mulsd %xmm8, %xmm12
mulsd %xmm12, %xmm9
movapd %xmm2, %xmm4
subsd %xmm7, %xmm4
mulsd %xmm4, %xmm5
addsd %xmm5, %xmm14
mulsd %xmm8, %xmm14
mulsd %xmm4, %xmm1
addsd %xmm1, %xmm13
mulsd %xmm8, %xmm13
movsd 0x10(%rsp), %xmm1
mulsd %xmm7, %xmm1
subsd %xmm1, %xmm15
addsd %xmm15, %xmm3
mulsd %xmm8, %xmm3
mulsd %xmm0, %xmm9
movsd %xmm9, (%rax)
movsd (%r14), %xmm1
mulsd %xmm2, %xmm1
mulsd %xmm0, %xmm14
addsd %xmm1, %xmm14
movsd %xmm14, 0x8(%rax)
movsd 0x8(%r14), %xmm1
mulsd %xmm2, %xmm1
mulsd %xmm0, %xmm13
addsd %xmm1, %xmm13
movsd %xmm13, 0x10(%rax)
addsd %xmm2, %xmm3
mulsd %xmm0, %xmm3
movapd %xmm3, %xmm4
movsd %xmm4, 0x18(%rax)
addq $0xc8, %rsp
popq %rbx
popq %r12
popq %r13
popq %r14
popq %r15
popq %rbp
retq