// File: /workspace/github_cpp_makefile_star2_O3/ashlall[P]StreamStats/ChiSquareExperiments/./../SourceCode/ChiSquare/Chi_Square_Distribution.cpp
// Function: pochisq(double, int)
double pochisq(double x, int df)
   {
        double a, s;
        double e, c, z;
 
        if (x <= 0.0 || df < 1) {
            return 1.0;
        }
        a = 0.5 * x;
        bool even = (df & 1) == 0;
        double y = 0;
        if (df > 1) {
            y = ex(-a);
        }
        s = (even ? y : (2.0 * poz(-sqrt(x))));
        if (df > 2) {
            x = 0.5 * (df - 1.0);
            z = (even ? 1.0 : 0.5);
            if (a > MAX_X) {
                e = (even ? 0.0 : LOG_SQRT_PI);
                c = log(a);
                while (z <= x) {
                    e = log(z) + e;
                    s += ex(c * z - a - e);
                    z += 1.0;
                }
                return s;
            } else {
                e = (even ? 1.0 : (I_SQRT_PI / sqrt(a)));
                c = 0.0;
                while (z <= x) {
                    e = e * (a / z);
                    c = c + e;
                    z += 1.0;
                }
                return c * y + s;
            }
        } else {
            return s;
        }
    }

// Assembly:
pushq %rbp
movq %rsp, %rbp
pushq %rbx
subq $0x38, %rsp
xorpd %xmm5, %xmm5
ucomisd %xmm0, %xmm5
setae %al
testl %edi, %edi
setle %cl
orb %al, %cl
je 0x40a9
movsd 0x4f64(%rip), %xmm4     # 0x9008
jmp 0x428f
movl %edi, %ebx
movapd %xmm0, %xmm2
movsd 0x4f89(%rip), %xmm3     # 0x9040
mulsd %xmm0, %xmm3
cmpl $0x2, %edi
movapd %xmm3, -0x40(%rbp)
jb 0x40fb
movapd 0x5103(%rip), %xmm0     # 0x91d0
xorpd %xmm3, %xmm0
xorl %eax, %eax
subl 0x90ff(%rip), %eax      # 0xd1d8
cvtsi2sd %eax, %xmm1
ucomisd %xmm0, %xmm1
ja 0x40fb
movsd %xmm2, -0x18(%rbp)
callq 0x2380
movsd -0x18(%rbp), %xmm2
movapd -0x40(%rbp), %xmm3
movapd %xmm0, %xmm5
movapd %xmm5, %xmm4
testb $0x1, %bl
je 0x4144
movsd %xmm5, -0x10(%rbp)
xorpd %xmm0, %xmm0
ucomisd %xmm0, %xmm2
jb 0x411c
xorps %xmm0, %xmm0
sqrtsd %xmm2, %xmm0
jmp 0x4125
movapd %xmm2, %xmm0
callq 0x2190
xorpd 0x50a3(%rip), %xmm0     # 0x91d0
callq 0x4355
movapd %xmm0, %xmm4
addsd %xmm0, %xmm4
movapd -0x40(%rbp), %xmm3
movsd -0x10(%rbp), %xmm5
cmpl $0x3, %ebx
jb 0x428f
xorl %eax, %eax
testb $0x1, %bl
sete %al
cvtsi2sd %ebx, %xmm6
addsd 0x4eb7(%rip), %xmm6     # 0x9018
mulsd 0x4ed7(%rip), %xmm6     # 0x9040
leaq 0x5070(%rip), %rcx      # 0x91e0
movsd (%rcx,%rax,8), %xmm7
xorps %xmm0, %xmm0
cvtsi2sdl 0x9058(%rip), %xmm0 # 0xd1d8
ucomisd %xmm0, %xmm3
movsd %xmm6, -0x30(%rbp)
jbe 0x41ca
movsd %xmm4, -0x20(%rbp)
movsd 0x9068(%rip), %xmm0     # 0xd200
movsd %xmm0, -0x10(%rbp)
movapd %xmm3, %xmm0
movsd %xmm7, -0x18(%rbp)
callq 0x21a0
movsd -0x18(%rbp), %xmm3
movsd %xmm0, -0x28(%rbp)
movsd -0x30(%rbp), %xmm0
ucomisd %xmm3, %xmm0
jae 0x41e1
movsd -0x20(%rbp), %xmm4
jmp 0x428f
testb $0x1, %bl
jne 0x429a
movsd 0x4e2c(%rip), %xmm8     # 0x9008
jmp 0x42f8
xorpd %xmm5, %xmm5
testb $0x1, %bl
movsd -0x20(%rbp), %xmm4
je 0x41f4
movsd -0x10(%rbp), %xmm5
movapd -0x40(%rbp), %xmm0
xorpd 0x4fcf(%rip), %xmm0     # 0x91d0
movapd %xmm0, -0x40(%rbp)
movsd %xmm3, -0x18(%rbp)
movsd %xmm4, -0x20(%rbp)
movsd %xmm5, -0x10(%rbp)
movapd %xmm3, %xmm0
callq 0x21a0
movsd -0x18(%rbp), %xmm3
movapd %xmm0, %xmm5
addsd -0x10(%rbp), %xmm5
movsd -0x28(%rbp), %xmm0
mulsd %xmm3, %xmm0
addsd -0x40(%rbp), %xmm0
xorl %eax, %eax
subl 0x8f96(%rip), %eax      # 0xd1d8
xorps %xmm2, %xmm2
cvtsi2sd %eax, %xmm2
subsd %xmm5, %xmm0
xorpd %xmm1, %xmm1
ucomisd %xmm0, %xmm2
ja 0x426f
movsd %xmm5, -0x10(%rbp)
callq 0x2380
movsd -0x10(%rbp), %xmm5
movsd -0x18(%rbp), %xmm3
movapd %xmm0, %xmm1
movsd -0x20(%rbp), %xmm4
addsd %xmm1, %xmm4
addsd 0x4d88(%rip), %xmm3     # 0x9008
movsd -0x30(%rbp), %xmm0
ucomisd %xmm3, %xmm0
jae 0x4206
movapd %xmm4, %xmm0
addq $0x38, %rsp
popq %rbx
popq %rbp
retq
movsd 0x8f65(%rip), %xmm8     # 0xd208
xorpd %xmm0, %xmm0
ucomisd %xmm0, %xmm3
jb 0x42b6
xorps %xmm0, %xmm0
sqrtsd %xmm3, %xmm0
jmp 0x42f3
movapd %xmm3, %xmm0
movsd %xmm4, -0x20(%rbp)
movsd %xmm5, -0x10(%rbp)
movsd %xmm7, -0x18(%rbp)
movsd %xmm8, -0x28(%rbp)
callq 0x2190
movsd -0x28(%rbp), %xmm8
movsd -0x18(%rbp), %xmm7
movsd -0x30(%rbp), %xmm6
movsd -0x10(%rbp), %xmm5
movsd -0x20(%rbp), %xmm4
movapd -0x40(%rbp), %xmm3
divsd %xmm0, %xmm8
xorpd %xmm0, %xmm0
ucomisd %xmm7, %xmm6
jb 0x4326
movsd 0x4cfe(%rip), %xmm1     # 0x9008
movapd %xmm3, %xmm2
divsd %xmm7, %xmm2
mulsd %xmm2, %xmm8
addsd %xmm8, %xmm0
addsd %xmm1, %xmm7
ucomisd %xmm7, %xmm6
jae 0x430a
mulsd %xmm5, %xmm0
addsd %xmm4, %xmm0
jmp 0x4293