// File: /workspace/github_c_makefile_star2_O0/samtools[P]htslib/kfunc.c
// Function: kf_erfc
double kf_erfc(double x)
{
	const double p0 = 220.2068679123761;
	const double p1 = 221.2135961699311;
	const double p2 = 112.0792914978709;
	const double p3 = 33.912866078383;
	const double p4 = 6.37396220353165;
	const double p5 = .7003830644436881;
	const double p6 = .03526249659989109;
	const double q0 = 440.4137358247522;
	const double q1 = 793.8265125199484;
	const double q2 = 637.3336333788311;
	const double q3 = 296.5642487796737;
	const double q4 = 86.78073220294608;
	const double q5 = 16.06417757920695;
	const double q6 = 1.755667163182642;
	const double q7 = .08838834764831844;
	double expntl, z, p;
	z = fabs(x) * M_SQRT2;
	if (z > 37.) return x > 0.? 0. : 2.;
	expntl = exp(z * z * - .5);
	if (z < 10. / M_SQRT2) // for small z
	    p = expntl * ((((((p6 * z + p5) * z + p4) * z + p3) * z + p2) * z + p1) * z + p0)
			/ (((((((q7 * z + q6) * z + q5) * z + q4) * z + q3) * z + q2) * z + q1) * z + q0);
	else p = expntl / 2.506628274631001 / (z + 1. / (z + 2. / (z + 3. / (z + 4. / (z + .65)))));
	return x > 0.? 2. * p : 2. * (1. - p);
}

// Assembly:
pushq %rbp
movq %rsp, %rbp
subq $0xb0, %rsp
movsd %xmm0, -0x10(%rbp)
movsd 0x1980(%rip), %xmm0     # 0x31f8
movsd %xmm0, -0x18(%rbp)
movsd 0x196b(%rip), %xmm0     # 0x31f0
movsd %xmm0, -0x20(%rbp)
movsd 0x1956(%rip), %xmm0     # 0x31e8
movsd %xmm0, -0x28(%rbp)
movsd 0x1941(%rip), %xmm0     # 0x31e0
movsd %xmm0, -0x30(%rbp)
movsd 0x192c(%rip), %xmm0     # 0x31d8
movsd %xmm0, -0x38(%rbp)
movsd 0x1917(%rip), %xmm0     # 0x31d0
movsd %xmm0, -0x40(%rbp)
movsd 0x1902(%rip), %xmm0     # 0x31c8
movsd %xmm0, -0x48(%rbp)
movsd 0x18ed(%rip), %xmm0     # 0x31c0
movsd %xmm0, -0x50(%rbp)
movsd 0x18d8(%rip), %xmm0     # 0x31b8
movsd %xmm0, -0x58(%rbp)
movsd 0x18c3(%rip), %xmm0     # 0x31b0
movsd %xmm0, -0x60(%rbp)
movsd 0x18ae(%rip), %xmm0     # 0x31a8
movsd %xmm0, -0x68(%rbp)
movsd 0x1899(%rip), %xmm0     # 0x31a0
movsd %xmm0, -0x70(%rbp)
movsd 0x1884(%rip), %xmm0     # 0x3198
movsd %xmm0, -0x78(%rbp)
movsd 0x186f(%rip), %xmm0     # 0x3190
movsd %xmm0, -0x80(%rbp)
movsd 0x185a(%rip), %xmm0     # 0x3188
movsd %xmm0, -0x88(%rbp)
movsd -0x10(%rbp), %xmm0
movaps 0x173e(%rip), %xmm1     # 0x3080
pand %xmm1, %xmm0
movsd 0x1832(%rip), %xmm1     # 0x3180
mulsd %xmm1, %xmm0
movsd %xmm0, -0x98(%rbp)
movsd -0x98(%rbp), %xmm0
movsd 0x180e(%rip), %xmm1     # 0x3178
ucomisd %xmm1, %xmm0
jbe 0x19a1
movsd -0x10(%rbp), %xmm1
xorps %xmm3, %xmm3
movsd 0x17a0(%rip), %xmm2     # 0x3120
movaps %xmm3, %xmm0
cmpltsd %xmm1, %xmm0
movaps %xmm0, %xmm1
andpd %xmm3, %xmm1
andnpd %xmm2, %xmm0
orpd %xmm1, %xmm0
movsd %xmm0, -0x8(%rbp)
jmp 0x1c2e
movsd -0x98(%rbp), %xmm0
mulsd -0x98(%rbp), %xmm0
movsd 0x184f(%rip), %xmm1     # 0x3208
mulsd %xmm1, %xmm0
callq 0x1060
movsd %xmm0, -0x90(%rbp)
movsd 0x182e(%rip), %xmm0     # 0x3200
ucomisd -0x98(%rbp), %xmm0
jbe 0x1b45
movsd -0x90(%rbp), %xmm0
movsd -0x98(%rbp), %xmm1
movsd 0x17d0(%rip), %xmm2     # 0x31c8
mulsd %xmm2, %xmm1
movsd 0x17cc(%rip), %xmm2     # 0x31d0
addsd %xmm2, %xmm1
movsd -0x98(%rbp), %xmm2
mulsd %xmm2, %xmm1
movsd 0x17bc(%rip), %xmm2     # 0x31d8
addsd %xmm2, %xmm1
movsd -0x98(%rbp), %xmm2
mulsd %xmm2, %xmm1
movsd 0x17ac(%rip), %xmm2     # 0x31e0
addsd %xmm2, %xmm1
movsd -0x98(%rbp), %xmm2
mulsd %xmm2, %xmm1
movsd 0x179c(%rip), %xmm2     # 0x31e8
addsd %xmm2, %xmm1
movsd -0x98(%rbp), %xmm2
mulsd %xmm2, %xmm1
movsd 0x178c(%rip), %xmm2     # 0x31f0
addsd %xmm2, %xmm1
movsd -0x98(%rbp), %xmm2
mulsd %xmm2, %xmm1
movsd 0x177c(%rip), %xmm2     # 0x31f8
addsd %xmm2, %xmm1
mulsd %xmm1, %xmm0
movsd -0x98(%rbp), %xmm1
movsd 0x16f4(%rip), %xmm2     # 0x3188
mulsd %xmm2, %xmm1
movsd 0x16f0(%rip), %xmm2     # 0x3190
addsd %xmm2, %xmm1
movsd -0x98(%rbp), %xmm2
mulsd %xmm2, %xmm1
movsd 0x16e0(%rip), %xmm2     # 0x3198
addsd %xmm2, %xmm1
movsd -0x98(%rbp), %xmm2
mulsd %xmm2, %xmm1
movsd 0x16d0(%rip), %xmm2     # 0x31a0
addsd %xmm2, %xmm1
movsd -0x98(%rbp), %xmm2
mulsd %xmm2, %xmm1
movsd 0x16c0(%rip), %xmm2     # 0x31a8
addsd %xmm2, %xmm1
movsd -0x98(%rbp), %xmm2
mulsd %xmm2, %xmm1
movsd 0x16b0(%rip), %xmm2     # 0x31b0
addsd %xmm2, %xmm1
movsd -0x98(%rbp), %xmm2
mulsd %xmm2, %xmm1
movsd 0x16a0(%rip), %xmm2     # 0x31b8
addsd %xmm2, %xmm1
movsd -0x98(%rbp), %xmm2
mulsd %xmm2, %xmm1
movsd 0x1690(%rip), %xmm2     # 0x31c0
addsd %xmm2, %xmm1
divsd %xmm1, %xmm0
movsd %xmm0, -0xa0(%rbp)
jmp 0x1bd5
movsd -0x90(%rbp), %xmm0
movsd 0x16c3(%rip), %xmm1     # 0x3218
divsd %xmm1, %xmm0
movsd -0x98(%rbp), %xmm1
movsd -0x98(%rbp), %xmm3
movsd -0x98(%rbp), %xmm4
movsd -0x98(%rbp), %xmm5
movsd 0x168f(%rip), %xmm6     # 0x3210
addsd -0x98(%rbp), %xmm6
movsd 0x15af(%rip), %xmm2     # 0x3140
divsd %xmm6, %xmm2
addsd %xmm2, %xmm5
movsd 0x158f(%rip), %xmm2     # 0x3130
divsd %xmm5, %xmm2
addsd %xmm2, %xmm4
movsd 0x156f(%rip), %xmm2     # 0x3120
divsd %xmm4, %xmm2
addsd %xmm2, %xmm3
movsd 0x144f(%rip), %xmm2     # 0x3010
divsd %xmm3, %xmm2
addsd %xmm2, %xmm1
divsd %xmm1, %xmm0
movsd %xmm0, -0xa0(%rbp)
movsd -0x10(%rbp), %xmm0
xorps %xmm1, %xmm1
ucomisd %xmm1, %xmm0
jbe 0x1bfd
movsd 0x1535(%rip), %xmm0     # 0x3120
mulsd -0xa0(%rbp), %xmm0
movsd %xmm0, -0xa8(%rbp)
jmp 0x1c21
movsd 0x140b(%rip), %xmm1     # 0x3010
subsd -0xa0(%rbp), %xmm1
movsd 0x150b(%rip), %xmm0     # 0x3120
mulsd %xmm1, %xmm0
movsd %xmm0, -0xa8(%rbp)
movsd -0xa8(%rbp), %xmm0
movsd %xmm0, -0x8(%rbp)
movsd -0x8(%rbp), %xmm0
addq $0xb0, %rsp
popq %rbp
retq
nopl (%rax)