// File: /workspace/github_cpp_makefile_star2_O1/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:
pushq %r15
pushq %r14
pushq %r12
pushq %rbx
subq $0x138, %rsp            # imm = 0x138
movq %rcx, %rbx
movq %rsi, %r14
xorps %xmm0, %xmm0
movaps %xmm0, 0xa0(%rsp)
movaps %xmm0, 0x90(%rsp)
movaps %xmm0, 0x80(%rsp)
movaps %xmm0, 0x70(%rsp)
movaps %xmm0, 0x120(%rsp)
movaps %xmm0, 0x110(%rsp)
movupd 0x10(%rdi), %xmm4
movapd %xmm4, 0xa0(%rsp)
movupd (%rdx), %xmm1
movupd (%rdi), %xmm5
movapd %xmm5, 0x90(%rsp)
movapd %xmm5, %xmm8
unpckhpd %xmm5, %xmm8            # xmm8 = xmm8[1],xmm5[1]
movapd %xmm4, %xmm6
unpckhpd %xmm4, %xmm6            # xmm6 = xmm6[1],xmm4[1]
movsd 0xf16f(%rip), %xmm2     # 0x1e148
mulsd %xmm5, %xmm2
movapd %xmm4, %xmm0
mulsd %xmm4, %xmm0
mulsd %xmm8, %xmm8
addsd %xmm0, %xmm8
movsd 0xf1b9(%rip), %xmm0     # 0x1e1b0
mulsd %xmm6, %xmm0
divsd %xmm5, %xmm0
xorpd %xmm7, %xmm7
ucomisd %xmm7, %xmm0
movapd %xmm1, 0xc0(%rsp)
jb 0xf018
sqrtsd %xmm0, %xmm0
jmp 0xf068
movsd %xmm2, 0x8(%rsp)
movapd %xmm5, 0x10(%rsp)
movapd %xmm6, 0x30(%rsp)
movapd %xmm4, 0x20(%rsp)
movapd %xmm8, 0x50(%rsp)
callq 0x31b0
movapd 0x50(%rsp), %xmm8
movapd 0x20(%rsp), %xmm4
xorpd %xmm7, %xmm7
movapd 0x30(%rsp), %xmm6
movapd 0x10(%rsp), %xmm5
movsd 0x8(%rsp), %xmm2
movapd 0xc0(%rsp), %xmm1
movsd %xmm0, 0x48(%rsp)
movapd %xmm5, %xmm3
shufpd $0x1, %xmm5, %xmm3      # xmm3 = xmm3[1],xmm5[0]
divsd 0xf129(%rip), %xmm6     # 0x1e1a8
unpcklpd %xmm4, %xmm5            # xmm5 = xmm5[0],xmm4[0]
mulsd %xmm8, %xmm2
movapd %xmm1, %xmm4
unpckhpd %xmm1, %xmm4            # xmm4 = xmm4[1],xmm1[1]
movsd 0xa0(%rsp), %xmm0
movapd %xmm4, 0xf0(%rsp)
mulsd %xmm4, %xmm0
movsd 0x98(%rsp), %xmm4
mulsd %xmm1, %xmm4
addsd %xmm0, %xmm4
movapd %xmm4, 0x50(%rsp)
movupd 0x10(%r14), %xmm1
movapd %xmm1, 0x80(%rsp)
movupd (%r14), %xmm4
movapd %xmm4, 0x70(%rsp)
movapd %xmm1, 0xd0(%rsp)
unpckhpd %xmm1, %xmm1            # xmm1 = xmm1[1,1]
movsd 0xf0c4(%rip), %xmm0     # 0x1e1b0
movapd %xmm1, 0xe0(%rsp)
mulsd %xmm1, %xmm0
movapd %xmm4, 0x60(%rsp)
divsd %xmm4, %xmm0
ucomisd %xmm7, %xmm0
jb 0xf10f
sqrtsd %xmm0, %xmm0
jmp 0xf144
movsd %xmm2, 0x8(%rsp)
movapd %xmm3, 0x20(%rsp)
movapd %xmm5, 0x10(%rsp)
movapd %xmm6, 0x30(%rsp)
callq 0x31b0
movapd 0x30(%rsp), %xmm6
movapd 0x10(%rsp), %xmm5
movapd 0x20(%rsp), %xmm3
movsd 0x8(%rsp), %xmm2
movapd %xmm0, 0xb0(%rsp)
mulpd %xmm5, %xmm3
movapd %xmm3, 0x20(%rsp)
addsd %xmm6, %xmm2
movsd %xmm2, 0x8(%rsp)
leaq 0x70(%rsp), %rsi
movups 0x8(%rsi), %xmm0
movaps %xmm0, 0x10(%rsp)
movaps 0x10(%rsi), %xmm0
movaps %xmm0, 0x30(%rsp)
leaq 0x90(%rsp), %r14
leaq 0x110(%rsp), %r15
leaq 0x108(%rsp), %r12
movq %r14, %rdi
movq %r15, %rdx
movq %r12, %rcx
callq 0xfe08
movapd 0x10(%rsp), %xmm10
movapd 0x50(%rsp), %xmm11
movapd %xmm11, %xmm5
subsd 0x48(%rsp), %xmm5
movsd (%r12), %xmm0
movapd 0x30(%rsp), %xmm1
movhpd 0x10(%r15), %xmm1       # xmm1 = xmm1[0],mem[0]
movapd 0xf0(%rsp), %xmm2
mulpd %xmm1, %xmm2
movapd 0xc0(%rsp), %xmm12
movapd %xmm12, %xmm1
unpcklpd %xmm12, %xmm1           # xmm1 = xmm1[0],xmm12[0]
movhpd 0x8(%r15), %xmm10       # xmm10 = xmm10[0],mem[0]
mulpd %xmm1, %xmm10
addpd %xmm2, %xmm10
movapd %xmm10, %xmm1
unpckhpd %xmm10, %xmm1           # xmm1 = xmm1[1],xmm10[1]
subsd %xmm0, %xmm1
minsd %xmm1, %xmm5
movsd (%r14), %xmm1
movsd 0x18(%r14), %xmm4
movapd %xmm11, %xmm2
mulsd %xmm1, %xmm2
xorpd %xmm3, %xmm3
ucomisd %xmm3, %xmm5
jbe 0xf269
movsd %xmm2, (%rbx)
movsd 0x8(%rsp), %xmm3
addsd %xmm4, %xmm3
unpcklpd %xmm4, %xmm4            # xmm4 = xmm4[0,0]
mulpd %xmm4, %xmm12
movapd %xmm11, %xmm0
unpcklpd %xmm11, %xmm0           # xmm0 = xmm0[0],xmm11[0]
movapd 0x20(%rsp), %xmm1
mulpd %xmm0, %xmm1
addpd %xmm12, %xmm1
movupd %xmm1, 0x8(%rbx)
mulsd %xmm11, %xmm3
jmp 0xf537
movapd 0x20(%rsp), %xmm15
movapd 0x60(%rsp), %xmm9
movapd %xmm9, %xmm6
unpckhpd %xmm9, %xmm6            # xmm6 = xmm6[1],xmm9[1]
movapd 0xd0(%rsp), %xmm8
movapd %xmm8, %xmm7
mulsd %xmm8, %xmm7
mulsd %xmm6, %xmm6
addsd %xmm7, %xmm6
movsd 0xeea3(%rip), %xmm7     # 0x1e148
mulsd %xmm9, %xmm7
mulsd %xmm6, %xmm7
movapd %xmm9, %xmm6
shufpd $0x1, %xmm9, %xmm6      # xmm6 = xmm6[1],xmm9[0]
unpcklpd %xmm8, %xmm9            # xmm9 = xmm9[0],xmm8[0]
mulpd %xmm6, %xmm9
movapd 0xe0(%rsp), %xmm14
divsd 0xeed2(%rip), %xmm14    # 0x1e1a8
addsd %xmm7, %xmm14
movapd 0xb0(%rsp), %xmm6
unpcklpd %xmm0, %xmm6            # xmm6 = xmm6[0],xmm0[0]
addpd %xmm10, %xmm6
movapd %xmm6, %xmm0
unpckhpd %xmm6, %xmm0            # xmm0 = xmm0[1],xmm6[1]
maxsd %xmm0, %xmm6
movsd 0x70(%rsp), %xmm7
movsd 0x88(%rsp), %xmm8
movapd %xmm10, %xmm0
mulsd %xmm7, %xmm0
ucomisd %xmm6, %xmm3
jbe 0xf354
movsd %xmm0, (%rbx)
addsd %xmm8, %xmm14
unpcklpd %xmm8, %xmm8            # xmm8 = xmm8[0,0]
mulpd %xmm8, %xmm12
movapd %xmm10, %xmm0
unpcklpd %xmm10, %xmm0           # xmm0 = xmm0[0],xmm10[0]
mulpd %xmm0, %xmm9
addpd %xmm12, %xmm9
movupd %xmm9, 0x8(%rbx)
mulsd %xmm10, %xmm14
movapd %xmm14, %xmm3
jmp 0xf537
movapd %xmm9, 0x60(%rsp)
movapd %xmm6, %xmm13
movapd %xmm6, %xmm9
movapd %xmm10, 0x10(%rsp)
subsd %xmm10, %xmm9
movapd %xmm5, %xmm10
subsd %xmm11, %xmm10
mulsd %xmm10, %xmm2
mulsd %xmm9, %xmm0
subsd %xmm2, %xmm0
addsd %xmm4, %xmm0
subsd %xmm8, %xmm0
movapd %xmm10, %xmm2
mulsd %xmm1, %xmm2
movapd %xmm7, %xmm3
mulsd %xmm9, %xmm3
subsd %xmm2, %xmm3
divsd %xmm3, %xmm0
movsd 0xedec(%rip), %xmm6     # 0x1e1a0
movapd %xmm11, %xmm3
subsd %xmm5, %xmm3
mulsd %xmm1, %xmm3
movapd %xmm11, %xmm2
subsd %xmm0, %xmm2
mulsd %xmm3, %xmm2
addsd %xmm4, %xmm2
movapd %xmm2, %xmm3
mulsd %xmm0, %xmm3
xorpd %xmm11, %xmm11
ucomisd %xmm5, %xmm11
jb 0xf489
ucomisd %xmm11, %xmm0
jbe 0xf489
subsd %xmm0, %xmm5
divsd %xmm5, %xmm6
mulsd %xmm6, %xmm1
mulsd %xmm10, %xmm1
movapd %xmm2, %xmm5
subsd %xmm4, %xmm5
unpcklpd %xmm5, %xmm5            # xmm5 = xmm5[0,0]
mulpd %xmm12, %xmm5
movapd %xmm10, %xmm7
unpcklpd %xmm10, %xmm7           # xmm7 = xmm7[0],xmm10[0]
mulpd %xmm15, %xmm7
addpd %xmm5, %xmm7
movapd %xmm6, %xmm5
unpcklpd %xmm6, %xmm5            # xmm5 = xmm5[0],xmm6[0]
mulpd %xmm7, %xmm5
movapd 0x50(%rsp), %xmm7
mulsd %xmm4, %xmm7
movsd 0x8(%rsp), %xmm13
mulsd %xmm10, %xmm13
subsd %xmm7, %xmm13
addsd %xmm13, %xmm3
mulsd %xmm6, %xmm3
mulsd %xmm0, %xmm1
movsd %xmm1, (%rbx)
movapd %xmm2, %xmm1
unpcklpd %xmm2, %xmm1            # xmm1 = xmm1[0],xmm2[0]
mulpd %xmm12, %xmm1
movapd %xmm0, %xmm4
unpcklpd %xmm0, %xmm4            # xmm4 = xmm4[0],xmm0[0]
mulpd %xmm5, %xmm4
addpd %xmm1, %xmm4
movupd %xmm4, 0x8(%rbx)
jmp 0xf52f
ucomisd %xmm11, %xmm13
jb 0xf53c
xorpd %xmm1, %xmm1
ucomisd %xmm0, %xmm1
jb 0xf53c
subsd %xmm0, %xmm13
divsd %xmm13, %xmm6
mulsd %xmm6, %xmm7
mulsd %xmm9, %xmm7
movapd %xmm2, %xmm1
subsd %xmm8, %xmm1
unpcklpd %xmm1, %xmm1            # xmm1 = xmm1[0,0]
mulpd %xmm12, %xmm1
movapd %xmm9, %xmm4
unpcklpd %xmm9, %xmm4            # xmm4 = xmm4[0],xmm9[0]
mulpd 0x60(%rsp), %xmm4
addpd %xmm1, %xmm4
movapd %xmm6, %xmm1
unpcklpd %xmm6, %xmm1            # xmm1 = xmm1[0],xmm6[0]
mulpd %xmm4, %xmm1
movapd 0x10(%rsp), %xmm4
mulsd %xmm8, %xmm4
mulsd %xmm9, %xmm14
subsd %xmm4, %xmm14
addsd %xmm14, %xmm3
mulsd %xmm6, %xmm3
mulsd %xmm0, %xmm7
movsd %xmm7, (%rbx)
movapd %xmm2, %xmm4
unpcklpd %xmm2, %xmm4            # xmm4 = xmm4[0],xmm2[0]
mulpd %xmm12, %xmm4
movapd %xmm0, %xmm5
unpcklpd %xmm0, %xmm5            # xmm5 = xmm5[0],xmm0[0]
mulpd %xmm1, %xmm5
addpd %xmm4, %xmm5
movupd %xmm5, 0x8(%rbx)
addsd %xmm2, %xmm3
mulsd %xmm0, %xmm3
movsd %xmm3, 0x18(%rbx)
addq $0x138, %rsp            # imm = 0x138
popq %rbx
popq %r12
popq %r14
popq %r15
retq