// File: melop[P]admixem/./newran02/newran.cpp
// Function: ChiSq::ChiSq(int, double)
ChiSq::ChiSq(int df, Real noncen)
{
  if (df<=0 || noncen<0.0) Throw(Logic_error("Newran: illegal parameters"));
  else if (df==1) { version=0; c1=new ChiSq1(noncen); }
  else if (noncen==0.0)
  {
     if (df==2) { version=1; c1=new Exponential(); }
     else { version=2; c1=new Gamma2(0.5*df); }
  }
  else if (df==2) { version=3; c1=new ChiSq1(noncen/2.0); }
  else if (df==3) { version=4; c1=new Exponential(); c2=new ChiSq1(noncen); }
  else { version=5; c1=new Gamma2(0.5*(df-1)); c2=new ChiSq1(noncen); }
  if (!c1 || (version>3 && !c2)) ErrorNoSpace();
  mean=df+noncen; var=2*df+4.0*noncen;
}

// Assembly:
pushq %rbp
movq %rsp, %rbp
subq $0x80, %rsp
movq %rdi, -0x8(%rbp)
movl %esi, -0xc(%rbp)
movsd %xmm0, -0x18(%rbp)
movq -0x8(%rbp), %rdi
movq %rdi, -0x30(%rbp)
callq 0x198a0
movq -0x30(%rbp), %rax
leaq 0x94a51(%rip), %rcx     # 0xa8c80
addq $0x10, %rcx
movq %rcx, (%rax)
cmpl $0x0, -0xc(%rbp)
jle 0x14246
xorps %xmm0, %xmm0
ucomisd -0x18(%rbp), %xmm0
jbe 0x142ac
movl $0x1, %edi
callq 0x8230
movq %rax, %rdi
movq %rdi, %rax
movq %rax, -0x38(%rbp)
leaq 0x627d2(%rip), %rsi     # 0x76a33
callq 0x12d90
jmp 0x14268
movq -0x38(%rbp), %rdi
leaq 0x94755(%rip), %rsi     # 0xa89c8
xorl %eax, %eax
movl %eax, %edx
callq 0x8820
jmp 0x14658
movq -0x38(%rbp), %rdi
movq %rax, %rcx
movl %edx, %eax
movq %rcx, -0x20(%rbp)
movl %eax, -0x24(%rbp)
callq 0x8370
jmp 0x14646
movq %rax, %rcx
movl %edx, %eax
movq %rcx, -0x20(%rbp)
movl %eax, -0x24(%rbp)
jmp 0x14646
cmpl $0x1, -0xc(%rbp)
jne 0x1430d
movq -0x30(%rbp), %rax
movl $0x0, 0x18(%rax)
movl $0x38, %edi
callq 0x84c0
movq %rax, -0x40(%rbp)
jmp 0x142cd
movq -0x40(%rbp), %rdi
movsd -0x18(%rbp), %xmm0
callq 0x14160
jmp 0x142dd
movq -0x30(%rbp), %rax
movq -0x40(%rbp), %rcx
movq %rcx, 0x8(%rax)
jmp 0x145dc
movq -0x40(%rbp), %rdi
movq %rax, %rcx
movl %edx, %eax
movq %rcx, -0x20(%rbp)
movl %eax, -0x24(%rbp)
movl $0x38, %esi
callq 0x84d0
jmp 0x14646
movsd -0x18(%rbp), %xmm0
xorps %xmm1, %xmm1
ucomisd %xmm1, %xmm0
jne 0x143fc
jp 0x143fc
cmpl $0x2, -0xc(%rbp)
jne 0x14393
movq -0x30(%rbp), %rax
movl $0x1, 0x18(%rax)
movl $0x28, %edi
callq 0x84c0
movq %rax, -0x48(%rbp)
jmp 0x14346
movq -0x48(%rbp), %rdi
xorps %xmm0, %xmm0
movaps %xmm0, 0x10(%rdi)
movaps %xmm0, (%rdi)
movq $0x0, 0x20(%rdi)
callq 0x19900
jmp 0x14363
movq -0x30(%rbp), %rax
movq -0x48(%rbp), %rcx
movq %rcx, 0x8(%rax)
jmp 0x143f7
movq -0x48(%rbp), %rdi
movq %rax, %rcx
movl %edx, %eax
movq %rcx, -0x20(%rbp)
movl %eax, -0x24(%rbp)
movl $0x28, %esi
callq 0x84d0
jmp 0x14646
movq -0x30(%rbp), %rax
movl $0x2, 0x18(%rax)
movl $0x40, %edi
callq 0x84c0
movq %rax, -0x50(%rbp)
jmp 0x143ae
movq -0x50(%rbp), %rdi
cvtsi2sdl -0xc(%rbp), %xmm0
movsd 0x622d1(%rip), %xmm1    # 0x76690
mulsd %xmm1, %xmm0
callq 0x15ba0
jmp 0x143ca
movq -0x30(%rbp), %rax
movq -0x50(%rbp), %rcx
movq %rcx, 0x8(%rax)
jmp 0x143f7
movq -0x50(%rbp), %rdi
movq %rax, %rcx
movl %edx, %eax
movq %rcx, -0x20(%rbp)
movl %eax, -0x24(%rbp)
movl $0x40, %esi
callq 0x84d0
jmp 0x14646
jmp 0x145da
cmpl $0x2, -0xc(%rbp)
jne 0x14469
movq -0x30(%rbp), %rax
movl $0x3, 0x18(%rax)
movl $0x38, %edi
callq 0x84c0
movq %rax, -0x58(%rbp)
jmp 0x1441d
movq -0x58(%rbp), %rdi
movsd -0x18(%rbp), %xmm0
movsd 0x62262(%rip), %xmm1    # 0x76690
mulsd %xmm1, %xmm0
callq 0x14160
jmp 0x14439
movq -0x30(%rbp), %rax
movq -0x58(%rbp), %rcx
movq %rcx, 0x8(%rax)
jmp 0x145d8
movq -0x58(%rbp), %rdi
movq %rax, %rcx
movl %edx, %eax
movq %rcx, -0x20(%rbp)
movl %eax, -0x24(%rbp)
movl $0x38, %esi
callq 0x84d0
jmp 0x14646
cmpl $0x3, -0xc(%rbp)
jne 0x14526
movq -0x30(%rbp), %rax
movl $0x4, 0x18(%rax)
movl $0x28, %edi
callq 0x84c0
movq %rax, -0x60(%rbp)
jmp 0x1448e
movq -0x60(%rbp), %rdi
xorps %xmm0, %xmm0
movaps %xmm0, 0x10(%rdi)
movaps %xmm0, (%rdi)
movq $0x0, 0x20(%rdi)
callq 0x19900
jmp 0x144ab
movq -0x30(%rbp), %rax
movq -0x60(%rbp), %rcx
movq %rcx, 0x8(%rax)
movl $0x38, %edi
callq 0x84c0
movq %rax, -0x68(%rbp)
jmp 0x144c7
movq -0x68(%rbp), %rdi
movsd -0x18(%rbp), %xmm0
callq 0x14160
jmp 0x144d7
movq -0x30(%rbp), %rax
movq -0x68(%rbp), %rcx
movq %rcx, 0x10(%rax)
jmp 0x145d6
movq -0x60(%rbp), %rdi
movq %rax, %rcx
movl %edx, %eax
movq %rcx, -0x20(%rbp)
movl %eax, -0x24(%rbp)
movl $0x28, %esi
callq 0x84d0
jmp 0x14646
movq -0x68(%rbp), %rdi
movq %rax, %rcx
movl %edx, %eax
movq %rcx, -0x20(%rbp)
movl %eax, -0x24(%rbp)
movl $0x38, %esi
callq 0x84d0
jmp 0x14646
movq -0x30(%rbp), %rax
movl $0x5, 0x18(%rax)
movl $0x40, %edi
callq 0x84c0
movq %rax, -0x70(%rbp)
jmp 0x14541
movq -0x70(%rbp), %rdi
movl -0xc(%rbp), %eax
decl %eax
cvtsi2sd %eax, %xmm0
movsd 0x6213a(%rip), %xmm1    # 0x76690
mulsd %xmm1, %xmm0
callq 0x15ba0
jmp 0x14561
movq -0x30(%rbp), %rax
movq -0x70(%rbp), %rcx
movq %rcx, 0x8(%rax)
movl $0x38, %edi
callq 0x84c0
movq %rax, -0x78(%rbp)
jmp 0x1457d
movq -0x78(%rbp), %rdi
movsd -0x18(%rbp), %xmm0
callq 0x14160
jmp 0x1458d
movq -0x30(%rbp), %rax
movq -0x78(%rbp), %rcx
movq %rcx, 0x10(%rax)
jmp 0x145d6
movq -0x70(%rbp), %rdi
movq %rax, %rcx
movl %edx, %eax
movq %rcx, -0x20(%rbp)
movl %eax, -0x24(%rbp)
movl $0x40, %esi
callq 0x84d0
jmp 0x14646
movq -0x78(%rbp), %rdi
movq %rax, %rcx
movl %edx, %eax
movq %rcx, -0x20(%rbp)
movl %eax, -0x24(%rbp)
movl $0x38, %esi
callq 0x84d0
jmp 0x14646
jmp 0x145d8
jmp 0x145da
jmp 0x145dc
jmp 0x145de
movq -0x30(%rbp), %rax
cmpq $0x0, 0x8(%rax)
je 0x145fe
movq -0x30(%rbp), %rax
cmpl $0x3, 0x18(%rax)
jle 0x14607
movq -0x30(%rbp), %rax
cmpq $0x0, 0x10(%rax)
jne 0x14607
callq 0x136a0
jmp 0x14605
jmp 0x14607
movq -0x30(%rbp), %rax
cvtsi2sdl -0xc(%rbp), %xmm0
addsd -0x18(%rbp), %xmm0
movsd %xmm0, 0x20(%rax)
movl -0xc(%rbp), %ecx
shll %ecx
cvtsi2sd %ecx, %xmm1
movsd -0x18(%rbp), %xmm0
movsd 0x62088(%rip), %xmm2    # 0x766b8
mulsd %xmm2, %xmm0
addsd %xmm1, %xmm0
movsd %xmm0, 0x28(%rax)
addq $0x80, %rsp
popq %rbp
retq
movq -0x30(%rbp), %rdi
callq 0x19e90
movq -0x20(%rbp), %rdi
callq 0x8860
nopl (%rax,%rax)