// File: 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
subq $0x80, %rsp
movsd %xmm0, -0x10(%rbp)
movl %edi, -0x14(%rbp)
xorps %xmm0, %xmm0
ucomisd -0x10(%rbp), %xmm0
jae 0x5283
cmpl $0x1, -0x14(%rbp)
jge 0x5295
movsd 0xcd7d(%rip), %xmm0     # 0x12008
movsd %xmm0, -0x8(%rbp)
jmp 0x553b
movsd 0xcd93(%rip), %xmm0     # 0x12030
mulsd -0x10(%rbp), %xmm0
movsd %xmm0, -0x20(%rbp)
movl -0x14(%rbp), %eax
andl $0x1, %eax
cmpl $0x0, %eax
sete %al
andb $0x1, %al
movb %al, -0x41(%rbp)
xorps %xmm0, %xmm0
movsd %xmm0, -0x50(%rbp)
cmpl $0x1, -0x14(%rbp)
jle 0x52ec
movsd -0x20(%rbp), %xmm0
movq %xmm0, %rax
movabsq $-0x8000000000000000, %rcx # imm = 0x8000000000000000
xorq %rcx, %rax
movq %rax, %xmm0
callq 0x5550
movsd %xmm0, -0x50(%rbp)
testb $0x1, -0x41(%rbp)
je 0x52fe
movsd -0x50(%rbp), %xmm0
movsd %xmm0, -0x58(%rbp)
jmp 0x5338
movsd -0x10(%rbp), %xmm0
callq 0x2130
movq %xmm0, %rax
movabsq $-0x8000000000000000, %rcx # imm = 0x8000000000000000
xorq %rcx, %rax
movq %rax, %xmm0
callq 0x55a0
movaps %xmm0, %xmm1
movsd 0xcd21(%rip), %xmm0     # 0x12050
mulsd %xmm1, %xmm0
movsd %xmm0, -0x58(%rbp)
movsd -0x58(%rbp), %xmm0
movsd %xmm0, -0x28(%rbp)
cmpl $0x2, -0x14(%rbp)
jle 0x5531
cvtsi2sdl -0x14(%rbp), %xmm1
movsd 0xccaf(%rip), %xmm0     # 0x12008
subsd %xmm0, %xmm1
movsd 0xcccb(%rip), %xmm0     # 0x12030
mulsd %xmm1, %xmm0
movsd %xmm0, -0x10(%rbp)
movb -0x41(%rbp), %al
movsd 0xcc8f(%rip), %xmm0     # 0x12008
movsd 0xccaf(%rip), %xmm1     # 0x12030
movsd %xmm1, -0x68(%rbp)
testb $0x1, %al
movsd %xmm0, -0x60(%rbp)
jne 0x5399
movsd -0x68(%rbp), %xmm0
movsd %xmm0, -0x60(%rbp)
movsd -0x60(%rbp), %xmm0
movsd %xmm0, -0x40(%rbp)
movsd -0x20(%rbp), %xmm0
cvtsi2sdl 0x14e10(%rip), %xmm1 # 0x1a1c0
ucomisd %xmm1, %xmm0
jbe 0x5479
testb $0x1, -0x41(%rbp)
je 0x53ca
xorps %xmm0, %xmm0
movsd %xmm0, -0x70(%rbp)
jmp 0x53d7
movsd 0x14e4e(%rip), %xmm0    # 0x1a220
movsd %xmm0, -0x70(%rbp)
movsd -0x70(%rbp), %xmm0
movsd %xmm0, -0x30(%rbp)
movsd -0x20(%rbp), %xmm0
callq 0x2140
movsd %xmm0, -0x38(%rbp)
movsd -0x40(%rbp), %xmm1
movsd -0x10(%rbp), %xmm0
ucomisd %xmm1, %xmm0
jb 0x546a
movsd -0x40(%rbp), %xmm0
callq 0x2140
addsd -0x30(%rbp), %xmm0
movsd %xmm0, -0x30(%rbp)
movsd -0x38(%rbp), %xmm0
movsd -0x40(%rbp), %xmm2
movsd -0x20(%rbp), %xmm1
movq %xmm1, %rax
movabsq $-0x8000000000000000, %rcx # imm = 0x8000000000000000
xorq %rcx, %rax
movq %rax, %xmm1
mulsd %xmm2, %xmm0
addsd %xmm1, %xmm0
subsd -0x30(%rbp), %xmm0
callq 0x5550
addsd -0x28(%rbp), %xmm0
movsd %xmm0, -0x28(%rbp)
movsd 0xcbaa(%rip), %xmm0     # 0x12008
addsd -0x40(%rbp), %xmm0
movsd %xmm0, -0x40(%rbp)
jmp 0x53f0
movsd -0x28(%rbp), %xmm0
movsd %xmm0, -0x8(%rbp)
jmp 0x553b
testb $0x1, -0x41(%rbp)
je 0x548e
movsd 0xcb81(%rip), %xmm0     # 0x12008
movsd %xmm0, -0x78(%rbp)
jmp 0x54b6
movsd 0x14d92(%rip), %xmm0    # 0x1a228
movsd %xmm0, -0x80(%rbp)
movsd -0x20(%rbp), %xmm0
callq 0x2130
movaps %xmm0, %xmm1
movsd -0x80(%rbp), %xmm0
divsd %xmm1, %xmm0
movsd %xmm0, -0x78(%rbp)
movsd -0x78(%rbp), %xmm0
movsd %xmm0, -0x30(%rbp)
xorps %xmm0, %xmm0
movsd %xmm0, -0x38(%rbp)
movsd -0x40(%rbp), %xmm1
movsd -0x10(%rbp), %xmm0
ucomisd %xmm1, %xmm0
jb 0x5513
movsd -0x30(%rbp), %xmm0
movsd -0x20(%rbp), %xmm1
divsd -0x40(%rbp), %xmm1
mulsd %xmm1, %xmm0
movsd %xmm0, -0x30(%rbp)
movsd -0x38(%rbp), %xmm0
addsd -0x30(%rbp), %xmm0
movsd %xmm0, -0x38(%rbp)
movsd 0xcb01(%rip), %xmm0     # 0x12008
addsd -0x40(%rbp), %xmm0
movsd %xmm0, -0x40(%rbp)
jmp 0x54c8
movsd -0x38(%rbp), %xmm0
movsd -0x50(%rbp), %xmm2
movsd -0x28(%rbp), %xmm1
mulsd %xmm2, %xmm0
addsd %xmm1, %xmm0
movsd %xmm0, -0x8(%rbp)
jmp 0x553b
movsd -0x28(%rbp), %xmm0
movsd %xmm0, -0x8(%rbp)
movsd -0x8(%rbp), %xmm0
addq $0x80, %rsp
popq %rbp
retq
nopl (%rax)