// File: jamesliu516[P]amrdg2d/hllc.cpp
// Function: HLLC_flux(privariables*, privariables*, pointxy*, double*)
void HLLC_flux(JBBL *ujbL, JBBL *ujbR, PXYZ *nml, 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, wl, wr;

	PXYZ normal;

	//Flux flxL, flxR, flxLs, flxRs;
	
	normal.x = nml->x;
	normal.y = nml->y;


	ul.q = ujbL->q;
	ul.u = ujbL->u;
	ul.v = ujbL->v;
	ul.p = ujbL->p;

	jbblToShbl(&ul, &wl);

	al = sqrt(GAMMA*ul.p/ul.q);

	ql = normal.x * ul.u + normal.y * ul.v;

	ur.q = ujbR->q;
	ur.u = ujbR->u;
	ur.v = ujbR->v;
	ur.p = ujbR->p;


	jbblToShbl(&ur, &wr);

	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;


	//	flxface->fq = flxL.fq;
	//	flxface->fqu = flxL.fqu;
	//	flxface->fqv = flxL.fqv;
	//	flxface->fte = flxL.fte;

	}
	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;


	//	flxface->fq = flxR.fq;
	//	flxface->fqu = flxR.fqu;
	//	flxface->fqv = flxR.fqv;
	//	flxface->fte = flxR.fte;
	}
	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;

	//	flxface->fq = flxLs.fq;
	//	flxface->fqu = flxLs.fqu;
	//	flxface->fqv = flxLs.fqv;
	//	flxface->fte = flxLs.fte;

	}
	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;

		//flxface->fq = flxRs.fq;
		//flxface->fqu = flxRs.fqu;
		//flxface->fqv = flxRs.fqv;
		//flxface->fte = flxRs.fte;	

	}
}

// Assembly:
subq $0x188, %rsp            # imm = 0x188
movq %rdi, 0x180(%rsp)
movq %rsi, 0x178(%rsp)
movq %rdx, 0x170(%rsp)
movq %rcx, 0x168(%rsp)
leaq 0xe8(%rsp), %rdi
callq 0x6530
leaq 0xc8(%rsp), %rdi
callq 0x6530
leaq 0xa8(%rsp), %rdi
callq 0x6530
leaq 0x88(%rsp), %rdi
callq 0x21aa0
leaq 0x68(%rsp), %rdi
callq 0x21aa0
leaq 0x48(%rsp), %rdi
callq 0x21aa0
leaq 0x28(%rsp), %rdi
callq 0x21aa0
leaq 0x18(%rsp), %rdi
callq 0x6310
movq 0x170(%rsp), %rax
movsd (%rax), %xmm0
movsd %xmm0, 0x18(%rsp)
movq 0x170(%rsp), %rax
movsd 0x8(%rax), %xmm0
movsd %xmm0, 0x20(%rsp)
movq 0x180(%rsp), %rax
movsd (%rax), %xmm0
movsd %xmm0, 0xe8(%rsp)
movq 0x180(%rsp), %rax
movsd 0x8(%rax), %xmm0
movsd %xmm0, 0xf0(%rsp)
movq 0x180(%rsp), %rax
movsd 0x10(%rax), %xmm0
movsd %xmm0, 0xf8(%rsp)
movq 0x180(%rsp), %rax
movsd 0x18(%rax), %xmm0
movsd %xmm0, 0x100(%rsp)
leaq 0xe8(%rsp), %rdi
leaq 0x48(%rsp), %rsi
callq 0x23570
movsd 0x1d50f(%rip), %xmm0    # 0x3f1d8
mulsd 0x100(%rsp), %xmm0
divsd 0xe8(%rsp), %xmm0
callq 0x31c0
movsd %xmm0, 0x128(%rsp)
movsd 0x18(%rsp), %xmm0
movsd 0xf0(%rsp), %xmm2
movsd 0x20(%rsp), %xmm1
mulsd 0xf8(%rsp), %xmm1
mulsd %xmm2, %xmm0
addsd %xmm1, %xmm0
movsd %xmm0, 0x138(%rsp)
movq 0x178(%rsp), %rax
movsd (%rax), %xmm0
movsd %xmm0, 0xc8(%rsp)
movq 0x178(%rsp), %rax
movsd 0x8(%rax), %xmm0
movsd %xmm0, 0xd0(%rsp)
movq 0x178(%rsp), %rax
movsd 0x10(%rax), %xmm0
movsd %xmm0, 0xd8(%rsp)
movq 0x178(%rsp), %rax
movsd 0x18(%rax), %xmm0
movsd %xmm0, 0xe0(%rsp)
leaq 0xc8(%rsp), %rdi
leaq 0x28(%rsp), %rsi
callq 0x23570
movsd 0x1d44f(%rip), %xmm0    # 0x3f1d8
mulsd 0xe0(%rsp), %xmm0
divsd 0xc8(%rsp), %xmm0
callq 0x31c0
movsd %xmm0, 0x120(%rsp)
movsd 0x18(%rsp), %xmm0
movsd 0xd0(%rsp), %xmm2
movsd 0x20(%rsp), %xmm1
mulsd 0xd8(%rsp), %xmm1
mulsd %xmm2, %xmm0
addsd %xmm1, %xmm0
movsd %xmm0, 0x130(%rsp)
leaq 0xe8(%rsp), %rdi
leaq 0xc8(%rsp), %rsi
leaq 0xa8(%rsp), %rdx
leaq 0x140(%rsp), %rcx
callq 0x23630
movsd 0x18(%rsp), %xmm0
movsd 0xb0(%rsp), %xmm2
movsd 0x20(%rsp), %xmm1
mulsd 0xb8(%rsp), %xmm1
mulsd %xmm2, %xmm0
addsd %xmm1, %xmm0
movsd %xmm0, 0x118(%rsp)
movsd 0x138(%rsp), %xmm1
subsd 0x128(%rsp), %xmm1
movsd 0x118(%rsp), %xmm0
subsd 0x140(%rsp), %xmm0
ucomisd %xmm1, %xmm0
jbe 0x21e70
movsd 0x138(%rsp), %xmm0
subsd 0x128(%rsp), %xmm0
movsd %xmm0, 0x10(%rsp)
jmp 0x21e88
movsd 0x118(%rsp), %xmm0
subsd 0x140(%rsp), %xmm0
movsd %xmm0, 0x10(%rsp)
movsd 0x10(%rsp), %xmm0
movsd %xmm0, 0x160(%rsp)
movsd 0x130(%rsp), %xmm0
addsd 0x120(%rsp), %xmm0
movsd 0x118(%rsp), %xmm1
addsd 0x140(%rsp), %xmm1
ucomisd %xmm1, %xmm0
jbe 0x21edb
movsd 0x130(%rsp), %xmm0
addsd 0x120(%rsp), %xmm0
movsd %xmm0, 0x8(%rsp)
jmp 0x21ef3
movsd 0x118(%rsp), %xmm0
addsd 0x140(%rsp), %xmm0
movsd %xmm0, 0x8(%rsp)
movsd 0x8(%rsp), %xmm0
movsd %xmm0, 0x158(%rsp)
movsd 0xc8(%rsp), %xmm0
mulsd 0x130(%rsp), %xmm0
movsd 0x158(%rsp), %xmm2
subsd 0x130(%rsp), %xmm2
movsd 0xe8(%rsp), %xmm1
mulsd 0x138(%rsp), %xmm1
movsd 0x160(%rsp), %xmm3
subsd 0x138(%rsp), %xmm3
mulsd %xmm3, %xmm1
movq %xmm1, %rax
movabsq $-0x8000000000000000, %rcx # imm = 0x8000000000000000
xorq %rcx, %rax
movq %rax, %xmm1
mulsd %xmm2, %xmm0
addsd %xmm1, %xmm0
addsd 0x100(%rsp), %xmm0
subsd 0xe0(%rsp), %xmm0
movsd 0xc8(%rsp), %xmm1
movsd 0x158(%rsp), %xmm3
subsd 0x130(%rsp), %xmm3
movsd 0xe8(%rsp), %xmm2
movsd 0x160(%rsp), %xmm4
subsd 0x138(%rsp), %xmm4
mulsd %xmm4, %xmm2
movq %xmm2, %rax
movabsq $-0x8000000000000000, %rcx # imm = 0x8000000000000000
xorq %rcx, %rax
movq %rax, %xmm2
mulsd %xmm3, %xmm1
addsd %xmm2, %xmm1
divsd %xmm1, %xmm0
movsd %xmm0, 0x150(%rsp)
movsd 0x160(%rsp), %xmm1
subsd 0x150(%rsp), %xmm1
movsd 0x1d059(%rip), %xmm0    # 0x3f058
divsd %xmm1, %xmm0
movsd %xmm0, 0x110(%rsp)
movsd 0x158(%rsp), %xmm1
subsd 0x150(%rsp), %xmm1
movsd 0x1d032(%rip), %xmm0    # 0x3f058
divsd %xmm1, %xmm0
movsd %xmm0, 0x108(%rsp)
movsd 0xe8(%rsp), %xmm0
movsd 0x138(%rsp), %xmm1
subsd 0x160(%rsp), %xmm1
mulsd %xmm1, %xmm0
movsd 0x138(%rsp), %xmm2
subsd 0x150(%rsp), %xmm2
movsd 0x100(%rsp), %xmm1
mulsd %xmm2, %xmm0
addsd %xmm1, %xmm0
movsd %xmm0, 0x148(%rsp)
movsd 0x110(%rsp), %xmm0
mulsd 0xe8(%rsp), %xmm0
movsd 0x160(%rsp), %xmm1
subsd 0x138(%rsp), %xmm1
mulsd %xmm1, %xmm0
movsd %xmm0, 0x88(%rsp)
movsd 0x110(%rsp), %xmm0
movsd 0x160(%rsp), %xmm1
subsd 0x138(%rsp), %xmm1
movsd 0x50(%rsp), %xmm3
movsd 0x148(%rsp), %xmm2
subsd 0x100(%rsp), %xmm2
mulsd 0x18(%rsp), %xmm2
mulsd %xmm3, %xmm1
addsd %xmm2, %xmm1
mulsd %xmm1, %xmm0
movsd %xmm0, 0x90(%rsp)
movsd 0x110(%rsp), %xmm0
movsd 0x160(%rsp), %xmm1
subsd 0x138(%rsp), %xmm1
movsd 0x58(%rsp), %xmm3
movsd 0x148(%rsp), %xmm2
subsd 0x100(%rsp), %xmm2
mulsd 0x20(%rsp), %xmm2
mulsd %xmm3, %xmm1
addsd %xmm2, %xmm1
mulsd %xmm1, %xmm0
movsd %xmm0, 0x98(%rsp)
movsd 0x110(%rsp), %xmm0
movsd 0x160(%rsp), %xmm2
subsd 0x138(%rsp), %xmm2
movsd 0x60(%rsp), %xmm3
movsd 0x100(%rsp), %xmm1
mulsd 0x138(%rsp), %xmm1
movq %xmm1, %rax
movabsq $-0x8000000000000000, %rcx # imm = 0x8000000000000000
xorq %rcx, %rax
movq %rax, %xmm1
mulsd %xmm3, %xmm2
addsd %xmm1, %xmm2
movsd 0x148(%rsp), %xmm1
movsd 0x150(%rsp), %xmm3
mulsd %xmm3, %xmm1
addsd %xmm2, %xmm1
mulsd %xmm1, %xmm0
movsd %xmm0, 0xa0(%rsp)
movsd 0x108(%rsp), %xmm0
mulsd 0xc8(%rsp), %xmm0
movsd 0x158(%rsp), %xmm1
subsd 0x130(%rsp), %xmm1
mulsd %xmm1, %xmm0
movsd %xmm0, 0x68(%rsp)
movsd 0x108(%rsp), %xmm0
movsd 0x158(%rsp), %xmm1
subsd 0x130(%rsp), %xmm1
movsd 0x30(%rsp), %xmm3
movsd 0x148(%rsp), %xmm2
subsd 0xe0(%rsp), %xmm2
mulsd 0x18(%rsp), %xmm2
mulsd %xmm3, %xmm1
addsd %xmm2, %xmm1
mulsd %xmm1, %xmm0
movsd %xmm0, 0x70(%rsp)
movsd 0x108(%rsp), %xmm0
movsd 0x158(%rsp), %xmm1
subsd 0x130(%rsp), %xmm1
movsd 0x38(%rsp), %xmm3
movsd 0x148(%rsp), %xmm2
subsd 0xe0(%rsp), %xmm2
mulsd 0x20(%rsp), %xmm2
mulsd %xmm3, %xmm1
addsd %xmm2, %xmm1
mulsd %xmm1, %xmm0
movsd %xmm0, 0x78(%rsp)
movsd 0x108(%rsp), %xmm0
movsd 0x158(%rsp), %xmm2
subsd 0x130(%rsp), %xmm2
movsd 0x40(%rsp), %xmm3
movsd 0xe0(%rsp), %xmm1
mulsd 0x130(%rsp), %xmm1
movq %xmm1, %rax
movabsq $-0x8000000000000000, %rcx # imm = 0x8000000000000000
xorq %rcx, %rax
movq %rax, %xmm1
mulsd %xmm3, %xmm2
addsd %xmm1, %xmm2
movsd 0x148(%rsp), %xmm1
movsd 0x150(%rsp), %xmm3
mulsd %xmm3, %xmm1
addsd %xmm2, %xmm1
mulsd %xmm1, %xmm0
movsd %xmm0, 0x80(%rsp)
movsd 0x160(%rsp), %xmm0
xorps %xmm1, %xmm1
ucomisd %xmm1, %xmm0
jbe 0x223c5
movsd 0xe8(%rsp), %xmm0
mulsd 0x138(%rsp), %xmm0
movq 0x168(%rsp), %rax
movsd %xmm0, (%rax)
movsd 0x50(%rsp), %xmm0
movsd 0x138(%rsp), %xmm2
movsd 0x100(%rsp), %xmm1
mulsd 0x18(%rsp), %xmm1
mulsd %xmm2, %xmm0
addsd %xmm1, %xmm0
movq 0x168(%rsp), %rax
movsd %xmm0, 0x8(%rax)
movsd 0x58(%rsp), %xmm0
movsd 0x138(%rsp), %xmm2
movsd 0x100(%rsp), %xmm1
mulsd 0x20(%rsp), %xmm1
mulsd %xmm2, %xmm0
addsd %xmm1, %xmm0
movq 0x168(%rsp), %rax
movsd %xmm0, 0x10(%rax)
movsd 0x60(%rsp), %xmm0
addsd 0x100(%rsp), %xmm0
mulsd 0x138(%rsp), %xmm0
movq 0x168(%rsp), %rax
movsd %xmm0, 0x18(%rax)
jmp 0x2263b
xorps %xmm0, %xmm0
ucomisd 0x158(%rsp), %xmm0
jbe 0x22485
movsd 0xc8(%rsp), %xmm0
mulsd 0x130(%rsp), %xmm0
movq 0x168(%rsp), %rax
movsd %xmm0, (%rax)
movsd 0x30(%rsp), %xmm0
movsd 0x130(%rsp), %xmm2
movsd 0xe0(%rsp), %xmm1
mulsd 0x18(%rsp), %xmm1
mulsd %xmm2, %xmm0
addsd %xmm1, %xmm0
movq 0x168(%rsp), %rax
movsd %xmm0, 0x8(%rax)
movsd 0x38(%rsp), %xmm0
movsd 0x130(%rsp), %xmm2
movsd 0xe0(%rsp), %xmm1
mulsd 0x20(%rsp), %xmm1
mulsd %xmm2, %xmm0
addsd %xmm1, %xmm0
movq 0x168(%rsp), %rax
movsd %xmm0, 0x10(%rax)
movsd 0x40(%rsp), %xmm0
addsd 0xe0(%rsp), %xmm0
mulsd 0x130(%rsp), %xmm0
movq 0x168(%rsp), %rax
movsd %xmm0, 0x18(%rax)
jmp 0x22639
xorps %xmm0, %xmm0
ucomisd 0x160(%rsp), %xmm0
jb 0x22564
movsd 0x150(%rsp), %xmm0
xorps %xmm1, %xmm1
ucomisd %xmm1, %xmm0
jbe 0x22564
movsd 0x88(%rsp), %xmm0
mulsd 0x150(%rsp), %xmm0
movq 0x168(%rsp), %rax
movsd %xmm0, (%rax)
movsd 0x90(%rsp), %xmm0
movsd 0x150(%rsp), %xmm2
movsd 0x148(%rsp), %xmm1
mulsd 0x18(%rsp), %xmm1
mulsd %xmm2, %xmm0
addsd %xmm1, %xmm0
movq 0x168(%rsp), %rax
movsd %xmm0, 0x8(%rax)
movsd 0x98(%rsp), %xmm0
movsd 0x150(%rsp), %xmm2
movsd 0x148(%rsp), %xmm1
mulsd 0x20(%rsp), %xmm1
mulsd %xmm2, %xmm0
addsd %xmm1, %xmm0
movq 0x168(%rsp), %rax
movsd %xmm0, 0x10(%rax)
movsd 0xa0(%rsp), %xmm0
addsd 0x148(%rsp), %xmm0
mulsd 0x150(%rsp), %xmm0
movq 0x168(%rsp), %rax
movsd %xmm0, 0x18(%rax)
jmp 0x22637
xorps %xmm0, %xmm0
ucomisd 0x150(%rsp), %xmm0
jb 0x22635
movsd 0x158(%rsp), %xmm0
xorps %xmm1, %xmm1
ucomisd %xmm1, %xmm0
jb 0x22635
movsd 0x68(%rsp), %xmm0
mulsd 0x150(%rsp), %xmm0
movq 0x168(%rsp), %rax
movsd %xmm0, (%rax)
movsd 0x70(%rsp), %xmm0
movsd 0x150(%rsp), %xmm2
movsd 0x148(%rsp), %xmm1
mulsd 0x18(%rsp), %xmm1
mulsd %xmm2, %xmm0
addsd %xmm1, %xmm0
movq 0x168(%rsp), %rax
movsd %xmm0, 0x8(%rax)
movsd 0x78(%rsp), %xmm0
movsd 0x150(%rsp), %xmm2
movsd 0x148(%rsp), %xmm1
mulsd 0x20(%rsp), %xmm1
mulsd %xmm2, %xmm0
addsd %xmm1, %xmm0
movq 0x168(%rsp), %rax
movsd %xmm0, 0x10(%rax)
movsd 0x80(%rsp), %xmm0
addsd 0x148(%rsp), %xmm0
mulsd 0x150(%rsp), %xmm0
movq 0x168(%rsp), %rax
movsd %xmm0, 0x18(%rax)
jmp 0x22637
jmp 0x22639
jmp 0x2263b
addq $0x188, %rsp            # imm = 0x188
retq
nopw %cs:(%rax,%rax)