LLM4Binary/decompile-bench-makefile huggingface dataset:
// File: /workspace/github_cpp_makefile_star2_O3/melop[P]admixem/./newran02/newran.cpp
// Function: PosGen::Build(bool)
void PosGen::Build(bool sym)                 // set up arrays
{
   #ifdef MONITOR
      cout << "building PosGen arrays\n";
   #endif
   int i;
   NotReady=false;
   sx=new Real[60]; sfx=new Real[60];
   if (!sx || !sfx) ErrorNoSpace();
   Real sxi=0.0; Real inc = sym ? 0.01 : 0.02;
   for (i=0; i<60; i++)
   {
      sx[i]=sxi; Real f1=Density(sxi); sfx[i]=f1;
      if (f1<=0.0) goto L20;
      sxi+=inc/f1;
   }
   Throw(Runtime_error("Newran: area too large"));
L20:
   if (i<50) Throw(Runtime_error("Newran: area too small"));
   xi = sym ? 2*i : i;
   return;
}

// Assembly:
pushq %rbp
pushq %r14
pushq %rbx
subq $0x10, %rsp
movl %esi, %ebp
movq %rdi, %rbx
movb $0x0, 0x20(%rdi)
movl $0x1e0, %edi            # imm = 0x1E0
callq 0x7030
movq %rax, 0x10(%rbx)
movl $0x1e0, %edi            # imm = 0x1E0
callq 0x7030
movq %rax, 0x18(%rbx)
testl %ebp, %ebp
jne 0xbe5d
movsd 0x2857d(%rip), %xmm0    # 0x343d8
jmp 0xbe65
movsd 0x2856b(%rip), %xmm0    # 0x343d0
movsd %xmm0, (%rsp)
xorl %r14d, %r14d
xorpd %xmm0, %xmm0
movq 0x10(%rbx), %rax
movsd %xmm0, (%rax,%r14,8)
movq (%rbx), %rax
movq %rbx, %rdi
movsd %xmm0, 0x8(%rsp)
callq *0x10(%rax)
movq 0x18(%rbx), %rax
movsd %xmm0, (%rax,%r14,8)
xorpd %xmm1, %xmm1
ucomisd %xmm0, %xmm1
jae 0xbee7
movsd (%rsp), %xmm1
divsd %xmm0, %xmm1
movsd 0x8(%rsp), %xmm0
addsd %xmm1, %xmm0
incq %r14
cmpq $0x3c, %r14
jne 0xbe71
movl $0x1, %edi
callq 0x7180
movq %rax, %rbx
leaq 0x288e9(%rip), %rsi     # 0x347b7
movq %rax, %rdi
callq 0xb924
leaq 0x3bdab(%rip), %rsi     # 0x47c88
movq %rbx, %rdi
xorl %edx, %edx
callq 0x7660
cmpq $0x31, %r14
jbe 0xbf08
movl %ebp, %ecx
shll %cl, %r14d
xorps %xmm0, %xmm0
cvtsi2sd %r14d, %xmm0
movsd %xmm0, 0x8(%rbx)
addq $0x10, %rsp
popq %rbx
popq %r14
popq %rbp
retq
movl $0x1, %edi
callq 0x7180
movq %rax, %rbx
leaq 0x288b2(%rip), %rsi     # 0x347ce
movq %rax, %rdi
callq 0xb924
jmp 0xbed6
jmp 0xbf28
movq %rax, %r14
movq %rbx, %rdi
callq 0x72a0
movq %r14, %rdi
callq 0x76a0
// File: /workspace/github_cpp_makefile_star2_O3/melop[P]admixem/./newran02/newran.cpp
// Function: PosGen::Next()
Real PosGen::Next()
{
   Real ak,y; int ir;
   if (NotReady) Build(false);
   do
   {
      Real r1=Random::Next();
      ir = (int)(r1*xi); Real sxi=sx[ir];
      ak=sxi+(sx[ir+1]-sxi)*Random::Next();
      y=sfx[ir]*Random::Next();
   }
   while ( y>=sfx[ir+1] && y>=Density(ak) );
   return ak;
}

// Assembly:
pushq %r14
pushq %rbx
subq $0x18, %rsp
movq %rdi, %rbx
cmpb $0x1, 0x20(%rdi)
jne 0xbf9a
movq %rbx, %rdi
xorl %esi, %esi
callq 0xbe22
callq 0xbba0
mulsd 0x8(%rbx), %xmm0
cvttsd2si %xmm0, %eax
movq 0x10(%rbx), %rcx
movslq %eax, %r14
movsd (%rcx,%r14,8), %xmm0
movsd %xmm0, 0x10(%rsp)
movsd 0x8(%rcx,%r14,8), %xmm1
subsd %xmm0, %xmm1
movsd %xmm1, 0x8(%rsp)
callq 0xbba0
mulsd 0x8(%rsp), %xmm0
addsd 0x10(%rsp), %xmm0
movsd %xmm0, 0x10(%rsp)
movq 0x18(%rbx), %rax
movsd (%rax,%r14,8), %xmm0
movsd %xmm0, 0x8(%rsp)
callq 0xbba0
mulsd 0x8(%rsp), %xmm0
movq 0x18(%rbx), %rax
ucomisd 0x8(%rax,%r14,8), %xmm0
jb 0xc030
movq (%rbx), %rax
movq %rbx, %rdi
movsd %xmm0, 0x8(%rsp)
movsd 0x10(%rsp), %xmm0
callq *0x10(%rax)
movsd 0x8(%rsp), %xmm1
ucomisd %xmm0, %xmm1
jae 0xbf9a
movsd 0x10(%rsp), %xmm0
addq $0x18, %rsp
popq %rbx
popq %r14
retq
// File: /workspace/github_cpp_makefile_star2_O3/melop[P]admixem/./newran02/newran.cpp
// Function: SymGen::Next()
Real SymGen::Next()
{
   Real s,ak,y; int ir;
   if (NotReady) Build(true);
   do
   {
      s=1.0;
      Real r1=Random::Next();
      if (r1>0.5) { s=-1.0; r1=1.0-r1; }
      ir = (int)(r1*xi); Real sxi=sx[ir];
      ak=sxi+(sx[ir+1]-sxi)*Random::Next();
      y=sfx[ir]*Random::Next();
   }
   while ( y>=sfx[ir+1] && y>=Density(ak) );
   return s*ak;
}

// Assembly:
pushq %r14
pushq %rbx
subq $0x38, %rsp
movq %rdi, %rbx
cmpb $0x1, 0x20(%rdi)
jne 0xc05b
movq %rbx, %rdi
movl $0x1, %esi
callq 0xbe22
callq 0xbba0
movsd 0x27fb0(%rip), %xmm1    # 0x34018
subsd %xmm0, %xmm1
movsd 0x2836c(%rip), %xmm2    # 0x343e0
cmpltsd %xmm0, %xmm2
andpd %xmm2, %xmm1
movapd %xmm2, 0x20(%rsp)
andnpd %xmm0, %xmm2
orpd %xmm1, %xmm2
mulsd 0x8(%rbx), %xmm2
cvttsd2si %xmm2, %eax
movq 0x10(%rbx), %rcx
movslq %eax, %r14
movsd (%rcx,%r14,8), %xmm0
movsd %xmm0, 0x10(%rsp)
movsd 0x8(%rcx,%r14,8), %xmm1
subsd %xmm0, %xmm1
movsd %xmm1, 0x8(%rsp)
callq 0xbba0
mulsd 0x8(%rsp), %xmm0
addsd 0x10(%rsp), %xmm0
movapd %xmm0, 0x10(%rsp)
movq 0x18(%rbx), %rax
movsd (%rax,%r14,8), %xmm0
movsd %xmm0, 0x8(%rsp)
callq 0xbba0
mulsd 0x8(%rsp), %xmm0
movq 0x18(%rbx), %rax
ucomisd 0x8(%rax,%r14,8), %xmm0
jb 0xc11c
movq (%rbx), %rax
movq %rbx, %rdi
movsd %xmm0, 0x8(%rsp)
movapd 0x10(%rsp), %xmm0
callq *0x10(%rax)
movsd 0x8(%rsp), %xmm1
ucomisd %xmm0, %xmm1
jae 0xc05b
movaps 0x27f6d(%rip), %xmm1    # 0x34090
movaps 0x10(%rsp), %xmm2
xorps %xmm2, %xmm1
movaps 0x20(%rsp), %xmm0
andps %xmm0, %xmm1
andnps %xmm2, %xmm0
orps %xmm1, %xmm0
addq $0x38, %rsp
popq %rbx
popq %r14
retq
nop
// File: /workspace/github_cpp_makefile_star2_O3/melop[P]admixem/./newran02/newran.cpp
// Function: AsymGen::Build()
void AsymGen::Build()                        // set up arrays
{
   #ifdef MONITOR
      cout << "building AsymGen arrays\n";
   #endif
   int i;
   NotReady=false;
   sx=new Real[121]; sfx=new Real[121];
   if (!sx || !sfx)  ErrorNoSpace();
   Real sxi=mode;
   for (i=0; i<120; i++)
   {
      sx[i]=sxi; Real f1=Density(sxi); sfx[i]=f1;
      if (f1<=0.0) goto L20;
      sxi+=0.01/f1;
   }
   Throw(Runtime_error("Newran: area too large (a)"));
L20:
   ic=i-1; sx[120]=sxi; sfx[120]=0.0;
   sxi=mode;
   for (; i<120; i++)
   {
      sx[i]=sxi; Real f1=Density(sxi); sfx[i]=f1;
      if (f1<=0.0) goto L30;
      sxi-=0.01/f1;
   }
   Throw(Runtime_error("Newran: area too large (b)"));
L30:
   if (i<100)  Throw(Runtime_error("Newran: area too small"));
   xi=i;
   return;
}

// Assembly:
pushq %r14
pushq %rbx
pushq %rax
movq %rdi, %rbx
movb $0x0, 0x24(%rdi)
movl $0x3c8, %edi            # imm = 0x3C8
callq 0x7030
movq %rax, 0x10(%rbx)
movl $0x3c8, %edi            # imm = 0x3C8
callq 0x7030
movq %rax, 0x18(%rbx)
movsd 0x28(%rbx), %xmm0
xorl %r14d, %r14d
movq 0x10(%rbx), %rax
movsd %xmm0, (%rax,%r14,8)
movq (%rbx), %rax
movq %rbx, %rdi
movsd %xmm0, (%rsp)
callq *0x10(%rax)
movq 0x18(%rbx), %rax
movsd %xmm0, (%rax,%r14,8)
xorpd %xmm1, %xmm1
ucomisd %xmm0, %xmm1
jae 0xc1f0
movsd 0x28217(%rip), %xmm1    # 0x343d0
divsd %xmm0, %xmm1
movsd (%rsp), %xmm0
addsd %xmm1, %xmm0
incq %r14
cmpq $0x78, %r14
jne 0xc185
movl $0x1, %edi
callq 0x7180
movq %rax, %rbx
leaq 0x28602(%rip), %rsi     # 0x347e5
movq %rax, %rdi
callq 0xb924
jmp 0xc284
leal -0x1(%r14), %ecx
movl %ecx, 0x20(%rbx)
movq 0x10(%rbx), %rcx
movsd (%rsp), %xmm0
movsd %xmm0, 0x3c0(%rcx)
movq $0x0, 0x3c0(%rax)
cmpq $0x77, %r14
ja 0xc268
movsd 0x28(%rbx), %xmm0
movq 0x10(%rbx), %rax
movsd %xmm0, (%rax,%r14,8)
movq (%rbx), %rax
movq %rbx, %rdi
movsd %xmm0, (%rsp)
callq *0x10(%rax)
movq 0x18(%rbx), %rax
movsd %xmm0, (%rax,%r14,8)
xorpd %xmm1, %xmm1
ucomisd %xmm0, %xmm1
jae 0xc295
movsd 0x28196(%rip), %xmm1    # 0x343e8
divsd %xmm0, %xmm1
movsd (%rsp), %xmm0
addsd %xmm1, %xmm0
incq %r14
cmpq $0x78, %r14
jne 0xc21e
movl $0x1, %edi
callq 0x7180
movq %rax, %rbx
leaq 0x28584(%rip), %rsi     # 0x34800
movq %rax, %rdi
callq 0xb924
leaq 0x3b9fd(%rip), %rsi     # 0x47c88
movq %rbx, %rdi
xorl %edx, %edx
callq 0x7660
cmpq $0x63, %r14
jbe 0xc2b0
xorps %xmm0, %xmm0
cvtsi2sd %r14d, %xmm0
movsd %xmm0, 0x8(%rbx)
addq $0x8, %rsp
popq %rbx
popq %r14
retq
movl $0x1, %edi
callq 0x7180
movq %rax, %rbx
leaq 0x2850a(%rip), %rsi     # 0x347ce
movq %rax, %rdi
callq 0xb924
jmp 0xc284
jmp 0xc2d2
jmp 0xc2d2
movq %rax, %r14
movq %rbx, %rdi
callq 0x72a0
movq %r14, %rdi
callq 0x76a0
nop
// File: /workspace/github_cpp_makefile_star2_O3/melop[P]admixem/./newran02/newran.cpp
// Function: AsymGen::Next()
Real AsymGen::Next()
{
   Real ak,y; int ir1;
   if (NotReady) Build();
   do
   {
      Real r1=Random::Next();
      int ir=(int)(r1*xi); Real sxi=sx[ir];
      ir1 = (ir==ic) ? 120 : ir+1;
      ak=sxi+(sx[ir1]-sxi)*Random::Next();
      y=sfx[ir]*Random::Next();
   }
   while ( y>=sfx[ir1] && y>=Density(ak) );
   return ak;
}

// Assembly:
pushq %rbp
pushq %r15
pushq %r14
pushq %rbx
subq $0x18, %rsp
movq %rdi, %rbx
cmpb $0x1, 0x24(%rdi)
jne 0xc301
movq %rbx, %rdi
callq 0xc156
movl $0x78, %ebp
callq 0xbba0
mulsd 0x8(%rbx), %xmm0
cvttsd2si %xmm0, %eax
movq 0x10(%rbx), %rcx
movslq %eax, %r14
leal 0x1(%r14), %eax
cmpl %r14d, 0x20(%rbx)
cmovel %ebp, %eax
movsd (%rcx,%r14,8), %xmm0
movsd %xmm0, 0x10(%rsp)
movslq %eax, %r15
movsd (%rcx,%r15,8), %xmm1
subsd %xmm0, %xmm1
movsd %xmm1, 0x8(%rsp)
callq 0xbba0
mulsd 0x8(%rsp), %xmm0
addsd 0x10(%rsp), %xmm0
movsd %xmm0, 0x10(%rsp)
movq 0x18(%rbx), %rax
movsd (%rax,%r14,8), %xmm0
movsd %xmm0, 0x8(%rsp)
callq 0xbba0
mulsd 0x8(%rsp), %xmm0
movq 0x18(%rbx), %rax
ucomisd (%rax,%r15,8), %xmm0
jb 0xc3a8
movq (%rbx), %rax
movq %rbx, %rdi
movsd %xmm0, 0x8(%rsp)
movsd 0x10(%rsp), %xmm0
callq *0x10(%rax)
movsd 0x8(%rsp), %xmm1
ucomisd %xmm0, %xmm1
jae 0xc306
movsd 0x10(%rsp), %xmm0
addq $0x18, %rsp
popq %rbx
popq %r14
popq %r15
popq %rbp
retq
nop
// File: /workspace/github_cpp_makefile_star2_O3/melop[P]admixem/./newran02/newran.cpp
// Function: ChiSq::Next()
Real ChiSq::Next()
{
   switch(version)
   {
   case 0: return c1->Next();
   case 1: case 2: return c1->Next()*2.0;
   case 3: return c1->Next() + c1->Next();
   case 4: case 5: Real s1 = c1->Next()*2.0; Real s2 = c2->Next();
	   return s1 + s2; // this is to make it work with Microsoft VC5
   }
   return 0;
}

// Assembly:
pushq %rbx
subq $0x10, %rsp
movl 0x18(%rdi), %eax
cmpq $0x5, %rax
ja 0xc9a9
movq %rdi, %rbx
leaq 0x27bc8(%rip), %rcx     # 0x34520
movslq (%rcx,%rax,4), %rax
addq %rcx, %rax
jmpq *%rax
movq 0x8(%rbx), %rdi
movq (%rdi), %rax
callq *(%rax)
addsd %xmm0, %xmm0
jmp 0xc9ad
movq 0x8(%rbx), %rdi
movq (%rdi), %rax
callq *(%rax)
addsd %xmm0, %xmm0
movsd %xmm0, 0x8(%rsp)
movq 0x10(%rbx), %rdi
jmp 0xc99c
movq 0x8(%rbx), %rdi
movq (%rdi), %rax
callq *(%rax)
movsd %xmm0, 0x8(%rsp)
movq 0x8(%rbx), %rdi
movq (%rdi), %rax
callq *(%rax)
addsd 0x8(%rsp), %xmm0
jmp 0xc9ad
xorpd %xmm0, %xmm0
addq $0x10, %rsp
popq %rbx
retq
movq 0x8(%rbx), %rdi
movq (%rdi), %rax
addq $0x10, %rsp
popq %rbx
jmpq *(%rax)
nop
// File: /workspace/github_cpp_makefile_star2_O3/melop[P]admixem/./newran02/newran.cpp
// Function: Poisson1::Density(double) const
Real Poisson1::Density(Real x) const             // Poisson density function
{
   if (x < 0.0) return 0.0;
   double ix = floor(x);                         // use integer part
   double l = ln_mu * ix - mu - ln_gamma(1.0 + ix);
   return  (l < -40.0) ? 0.0 : exp(l);
}

// Assembly:
xorpd %xmm1, %xmm1
ucomisd %xmm0, %xmm1
ja 0xcb60
pushq %rbx
subq $0x10, %rsp
movq %rdi, %rbx
callq 0x7740
movsd 0x38(%rbx), %xmm1
mulsd %xmm0, %xmm1
subsd 0x30(%rbx), %xmm1
movsd %xmm1, 0x8(%rsp)
addsd 0x274db(%rip), %xmm0    # 0x34018
callq 0xcb6e
movsd 0x8(%rsp), %xmm1
subsd %xmm0, %xmm1
movsd 0x278e4(%rip), %xmm0    # 0x34438
ucomisd %xmm1, %xmm0
leaq 0x10(%rsp), %rsp
popq %rbx
jbe 0xcb65
xorpd %xmm0, %xmm0
retq
movapd %xmm1, %xmm0
jmp 0x7748
// File: /workspace/github_cpp_makefile_star2_O3/melop[P]admixem/./newran02/newran.cpp
// Function: ln_gamma(double)
Real ln_gamma(Real xx)
{
   // log gamma function adapted from numerical recipes in C

   if (xx<1.0)                           // Use reflection formula
   {
      double piz = 3.14159265359 * (1.0-xx);
      return log(piz/sin(piz))-ln_gamma(2.0-xx);
   }
   else
   {
      static double cof[6]={76.18009173,-86.50532033,24.01409822,
         -1.231739516,0.120858003e-2,-0.536382e-5};

      double x=xx-1.0; double tmp=x+5.5;
      tmp -= (x+0.5)*log(tmp); double ser=1.0;
      for (int j=0;j<=5;j++) { x += 1.0; ser += cof[j]/x; }
      return -tmp+log(2.50662827465*ser);
   }
}

// Assembly:
subq $0x28, %rsp
movapd %xmm0, %xmm1
movsd 0x2749a(%rip), %xmm0    # 0x34018
ucomisd %xmm1, %xmm0
jbe 0xcbe4
movsd 0x2748c(%rip), %xmm0    # 0x34018
subsd %xmm1, %xmm0
mulsd 0x278b8(%rip), %xmm0    # 0x34450
movsd %xmm0, (%rsp)
movsd %xmm1, 0x18(%rsp)
callq 0x7718
movsd (%rsp), %xmm1
divsd %xmm0, %xmm1
movapd %xmm1, %xmm0
callq 0x7730
movsd %xmm0, (%rsp)
movsd 0x27851(%rip), %xmm0    # 0x34418
subsd 0x18(%rsp), %xmm0
callq 0xcb6e
movsd (%rsp), %xmm1
subsd %xmm0, %xmm1
movapd %xmm1, %xmm0
jmp 0xcc82
addsd 0x2781c(%rip), %xmm1    # 0x34408
movsd %xmm1, 0x18(%rsp)
movsd 0x27846(%rip), %xmm0    # 0x34440
addsd %xmm1, %xmm0
movsd 0x277da(%rip), %xmm2    # 0x343e0
addsd %xmm1, %xmm2
movapd %xmm2, (%rsp)
movsd %xmm0, 0x20(%rsp)
callq 0x7730
movsd 0x273f6(%rip), %xmm5    # 0x34018
movsd 0x18(%rsp), %xmm3
movapd (%rsp), %xmm4
xorpd 0x2745b(%rip), %xmm4    # 0x34090
xorl %eax, %eax
leaq 0x27902(%rip), %rcx     # 0x34540
movapd %xmm5, %xmm1
addsd %xmm5, %xmm3
movsd (%rax,%rcx), %xmm2
divsd %xmm3, %xmm2
addsd %xmm2, %xmm1
addq $0x8, %rax
cmpq $0x30, %rax
jne 0xcc42
mulsd %xmm0, %xmm4
addsd 0x20(%rsp), %xmm4
movapd %xmm4, (%rsp)
mulsd 0x277d4(%rip), %xmm1    # 0x34448
movapd %xmm1, %xmm0
callq 0x7730
subsd (%rsp), %xmm0
addq $0x28, %rsp
retq
nop
// File: /workspace/github_cpp_makefile_star2_O3/melop[P]admixem/./newran02/newran.cpp
// Function: Binomial1::Density(double) const
Real Binomial1::Density(Real x) const            // Binomial density function
{
   double ix = floor(x);                         // use integer part
   if (ix < 0.0 || ix > n) return 0.0;
   double l = ln_n_fac - ln_gamma(ix+1) - ln_gamma(n-ix+1)
      + ix * ln_p + (n-ix) * ln_q;
   return  (l < -40.0) ? 0.0 : exp(l);
}

// Assembly:
pushq %rbx
subq $0x10, %rsp
movq %rdi, %rbx
callq 0x7740
movapd %xmm0, %xmm1
xorpd %xmm0, %xmm0
ucomisd %xmm1, %xmm0
ja 0xcdae
xorps %xmm0, %xmm0
cvtsi2sdl 0x58(%rbx), %xmm0
ucomisd %xmm0, %xmm1
ja 0xcdae
movsd 0x50(%rbx), %xmm0
movsd %xmm0, (%rsp)
movapd %xmm1, %xmm0
addsd 0x272dd(%rip), %xmm0    # 0x34018
movsd %xmm1, 0x8(%rsp)
callq 0xcb6e
movsd (%rsp), %xmm1
subsd %xmm0, %xmm1
movsd %xmm1, (%rsp)
xorps %xmm0, %xmm0
cvtsi2sdl 0x58(%rbx), %xmm0
subsd 0x8(%rsp), %xmm0
addsd 0x272ae(%rip), %xmm0    # 0x34018
callq 0xcb6e
movsd (%rsp), %xmm2
subsd %xmm0, %xmm2
movsd 0x40(%rbx), %xmm1
movsd 0x8(%rsp), %xmm3
mulsd %xmm3, %xmm1
xorps %xmm0, %xmm0
cvtsi2sdl 0x58(%rbx), %xmm0
addsd %xmm2, %xmm1
subsd %xmm3, %xmm0
mulsd 0x48(%rbx), %xmm0
addsd %xmm1, %xmm0
movsd 0x27690(%rip), %xmm1    # 0x34438
ucomisd %xmm0, %xmm1
jbe 0xcdb8
xorpd %xmm0, %xmm0
addq $0x10, %rsp
popq %rbx
retq
addq $0x10, %rsp
popq %rbx
jmp 0x7748
// File: /workspace/github_cpp_makefile_star2_O3/melop[P]admixem/./newran02/newran.cpp
// Function: Poisson2::Poisson2(double)
Poisson2::Poisson2(Real mux)
{
   Real probs[40];
   probs[0]=exp(-mux);
   for (int i=1; i<40; i++) probs[i]=probs[i-1]*mux/i;
   dg=new DiscreteGen(40,probs);
   if (!dg) ErrorNoSpace();
}

// Assembly:
pushq %r14
pushq %rbx
subq $0x158, %rsp            # imm = 0x158
movapd %xmm0, %xmm1
movapd %xmm0, (%rsp)
movq %rdi, %rbx
leaq 0x3a021(%rip), %rax     # 0x46e00
movq %rax, (%rdi)
movapd 0x272a6(%rip), %xmm0    # 0x34090
xorpd %xmm1, %xmm0
callq 0x7748
movapd (%rsp), %xmm2
movsd %xmm0, 0x10(%rsp)
movl $0x1, %eax
mulsd %xmm2, %xmm0
xorps %xmm1, %xmm1
cvtsi2sd %eax, %xmm1
divsd %xmm1, %xmm0
movsd %xmm0, 0x10(%rsp,%rax,8)
incq %rax
cmpq $0x28, %rax
jne 0xce03
movl $0x38, %edi
callq 0x7380
movq %rax, %r14
leaq 0x10(%rsp), %rdx
movq %rax, %rdi
movl $0x28, %esi
callq 0xd5d6
movq %r14, 0x8(%rbx)
addq $0x158, %rsp            # imm = 0x158
popq %rbx
popq %r14
retq
movq %rax, %rbx
movl $0x38, %esi
movq %r14, %rdi
callq 0x7390
movq %rbx, %rdi
callq 0x76a0
nop
// File: /workspace/github_cpp_makefile_star2_O3/melop[P]admixem/./newran02/newran.cpp
// Function: Binomial2::Binomial2(int, double)
Binomial2::Binomial2(int nx, Real px)
{
   Real qx = 1.0 - px;
   Real probs[40];
   int k = (int)(nx * px);
   probs[k] = exp(ln_gamma(nx+1) - ln_gamma(k+1) - ln_gamma(nx-k+1)
      + k * log(px) + (nx-k) * log(qx));
   int i;
   int m = (nx >= 40) ? 39 : nx;
   for (i=k+1; i<=m; i++) probs[i]=probs[i-1] * px * (nx-i+1) / qx / i;
   for (i=k-1; i>=0; i--) probs[i]=probs[i+1] * qx * (i+1) / px / (nx-i);
   dg = new DiscreteGen(m + 1, probs);
   if (!dg) ErrorNoSpace();
}

// Assembly:
pushq %rbp
pushq %r15
pushq %r14
pushq %r13
pushq %r12
pushq %rbx
subq $0x188, %rsp            # imm = 0x188
movsd %xmm0, 0x28(%rsp)
movl %esi, %ebx
leaq 0x39f92(%rip), %rax     # 0x46e60
movq %rdi, 0x30(%rsp)
movq %rax, (%rdi)
cvtsi2sd %esi, %xmm2
movsd 0x27136(%rip), %xmm1    # 0x34018
mulsd %xmm0, %xmm2
movapd %xmm2, 0x10(%rsp)
cvttsd2si %xmm2, %r12d
subsd %xmm0, %xmm1
movsd %xmm1, 0x20(%rsp)
leal 0x1(%rbx), %eax
xorps %xmm0, %xmm0
cvtsi2sd %eax, %xmm0
callq 0xcb6e
movsd %xmm0, 0x8(%rsp)
leal 0x1(%r12), %r13d
xorps %xmm0, %xmm0
cvtsi2sd %r13d, %xmm0
callq 0xcb6e
movsd %xmm0, 0x38(%rsp)
movl %ebx, %r14d
subl %r12d, %r14d
leal 0x1(%r14), %r15d
xorps %xmm0, %xmm0
cvtsi2sd %r15d, %xmm0
callq 0xcb6e
movsd 0x8(%rsp), %xmm1
subsd 0x38(%rsp), %xmm1
subsd %xmm0, %xmm1
movsd %xmm1, 0x8(%rsp)
cvttpd2dq 0x10(%rsp), %xmm0
cvtdq2pd %xmm0, %xmm0
movaps %xmm0, 0x10(%rsp)
movsd 0x28(%rsp), %xmm0
callq 0x7730
mulsd 0x10(%rsp), %xmm0
addsd 0x8(%rsp), %xmm0
movsd %xmm0, 0x8(%rsp)
xorps %xmm0, %xmm0
cvtsi2sd %r14d, %xmm0
movsd %xmm0, 0x10(%rsp)
movsd 0x20(%rsp), %xmm0
callq 0x7730
mulsd 0x10(%rsp), %xmm0
addsd 0x8(%rsp), %xmm0
callq 0x7748
movsd 0x20(%rsp), %xmm3
movsd 0x28(%rsp), %xmm2
movslq %r12d, %rax
movsd %xmm0, 0x40(%rsp,%rax,8)
cmpl $0x27, %ebx
movl $0x27, %ebp
cmovll %ebx, %ebp
cmpl %eax, %ebp
jle 0xd01d
leaq (%rsp,%rax,8), %rcx
addq $0x48, %rcx
movl %ebp, %edx
negl %edx
mulsd %xmm2, %xmm0
xorps %xmm1, %xmm1
cvtsi2sd %r14d, %xmm1
mulsd %xmm0, %xmm1
divsd %xmm3, %xmm1
xorps %xmm0, %xmm0
cvtsi2sd %r13d, %xmm0
divsd %xmm0, %xmm1
movsd %xmm1, (%rcx)
decl %r14d
addq $0x8, %rcx
leal (%rdx,%r13), %esi
incl %esi
movl %r13d, %edi
incl %edi
movl %edi, %r13d
movapd %xmm1, %xmm0
cmpl $0x1, %esi
jne 0xcfdb
testl %r12d, %r12d
jle 0xd061
movsd 0x40(%rsp,%rax,8), %xmm0
xorps %xmm1, %xmm1
cvtsi2sd %eax, %xmm1
mulsd %xmm3, %xmm0
mulsd %xmm0, %xmm1
xorps %xmm0, %xmm0
cvtsi2sd %r15d, %xmm0
divsd %xmm2, %xmm1
divsd %xmm0, %xmm1
movsd %xmm1, 0x38(%rsp,%rax,8)
leaq -0x1(%rax), %rcx
incl %r15d
movapd %xmm1, %xmm0
cmpq $0x1, %rax
movq %rcx, %rax
ja 0xd028
movl $0x38, %edi
callq 0x7380
movq %rax, %r14
incl %ebp
leaq 0x40(%rsp), %rdx
movq %rax, %rdi
movl %ebp, %esi
callq 0xd5d6
movq 0x30(%rsp), %rax
movq %r14, 0x8(%rax)
addq $0x188, %rsp            # imm = 0x188
popq %rbx
popq %r12
popq %r13
popq %r14
popq %r15
popq %rbp
retq
movq %rax, %rbx
movl $0x38, %esi
movq %r14, %rdi
callq 0x7390
movq %rbx, %rdi
callq 0x76a0
// File: /workspace/github_cpp_makefile_star2_O3/melop[P]admixem/./newran02/newran.cpp
// Function: NegativeBinomial::Density(double) const
Real NegativeBinomial::Density(Real x) const
{
   if (x < 0.0) return 0.0;
   Real ix = floor(x);
   Real l = c + ln_gamma(ix + N) - ln_gamma(ix + 1) + ix * ln_q;
   return  (l < -40.0) ? 0.0 : exp(l);
}

// Assembly:
xorpd %xmm1, %xmm1
ucomisd %xmm0, %xmm1
ja 0xd392
pushq %rbx
subq $0x10, %rsp
movq %rdi, %rbx
callq 0x7740
movapd %xmm0, %xmm1
movsd %xmm0, 0x8(%rsp)
movsd 0x58(%rbx), %xmm0
movsd %xmm0, (%rsp)
movsd 0x30(%rbx), %xmm0
addsd %xmm1, %xmm0
callq 0xcb6e
addsd (%rsp), %xmm0
movsd %xmm0, (%rsp)
movsd 0x26cbd(%rip), %xmm0    # 0x34018
addsd 0x8(%rsp), %xmm0
callq 0xcb6e
movsd (%rsp), %xmm2
subsd %xmm0, %xmm2
movsd 0x8(%rsp), %xmm1
mulsd 0x50(%rbx), %xmm1
addsd %xmm2, %xmm1
movsd 0x270b2(%rip), %xmm0    # 0x34438
ucomisd %xmm1, %xmm0
leaq 0x10(%rsp), %rsp
popq %rbx
jbe 0xd397
xorpd %xmm0, %xmm0
retq
movapd %xmm1, %xmm0
jmp 0x7748
// File: /workspace/github_cpp_makefile_star2_O3/melop[P]admixem/./newran02/newran.cpp
// Function: Gamma::Gamma(double)
Gamma::Gamma(Real alpha)                         // general gamma generator
{
   if (alpha<1.0) method = new Gamma1(alpha);
   else if (alpha==1.0) method = new Exponential();
   else method = new Gamma2(alpha);
   if (!method)  ErrorNoSpace();
}

// Assembly:
pushq %r14
pushq %rbx
pushq %rax
movq %rdi, %rbx
leaq 0x39bc4(%rip), %rax     # 0x470a0
movq %rax, (%rdi)
movsd 0x26b31(%rip), %xmm1    # 0x34018
ucomisd %xmm0, %xmm1
jbe 0xd541
movl $0x40, %edi
movsd %xmm0, (%rsp)
callq 0x7380
movq %rax, %r14
movb $0x1, 0x20(%rax)
leaq 0x39ad6(%rip), %rax     # 0x46fe0
movq %rax, (%r14)
movsd (%rsp), %xmm1
movapd %xmm1, %xmm0
movsd 0x26afa(%rip), %xmm2    # 0x34018
addsd %xmm2, %xmm0
divsd %xmm1, %xmm2
movsd %xmm2, 0x30(%r14)
callq 0xcb6e
movsd %xmm0, 0x28(%r14)
movsd (%rsp), %xmm0
jmp 0xd5c4
ucomisd 0x26acf(%rip), %xmm0    # 0x34018
jne 0xd57f
jp 0xd57f
movl $0x28, %edi
callq 0x7380
movq %rax, %r14
xorpd %xmm0, %xmm0
movapd %xmm0, (%rax)
movapd %xmm0, 0x10(%rax)
movq $0x0, 0x20(%rax)
movb $0x1, 0x20(%rax)
leaq 0x3a53e(%rip), %rax     # 0x47ab8
movq %rax, (%r14)
jmp 0xd5ca
movl $0x40, %edi
movsd %xmm0, (%rsp)
callq 0x7380
movq %rax, %r14
movsd 0x26e6f(%rip), %xmm0    # 0x34408
movsd (%rsp), %xmm1
addsd %xmm1, %xmm0
movsd %xmm0, 0x28(%rax)
movb $0x1, 0x24(%rax)
leaq 0x39a8e(%rip), %rax     # 0x47040
movq %rax, (%r14)
movsd %xmm1, 0x30(%r14)
movapd %xmm1, %xmm0
callq 0xcb6e
movsd %xmm0, 0x38(%r14)
movq %r14, 0x8(%rbx)
addq $0x8, %rsp
popq %rbx
popq %r14
retq
// File: /workspace/github_cpp_makefile_star2_O3/melop[P]admixem/./newran02/newran.cpp
// Function: DiscreteGen::DiscreteGen(int, double*)
DiscreteGen::DiscreteGen(int n1, Real* prob)     // discrete generator
						 // values on 0,...,n1-1
{
   #ifdef MONITOR
      cout << "constructing DiscreteGen\n";
   #endif
   Gen(n1, prob); val=0;
   mean=0.0; var=0.0;
   { for (int i=0; i<n; i++) mean = mean + i*prob[i]; }
   { for (int i=0; i<n; i++) var = var + square(i-mean) * prob[i]; }
}

// Assembly:
pushq %r14
pushq %rbx
pushq %rax
movq %rdx, %rbx
movq %rdi, %r14
leaq 0x39b19(%rip), %rax     # 0x47100
movq %rax, (%rdi)
callq 0xd666
xorpd %xmm0, %xmm0
movupd %xmm0, 0x20(%r14)
movq $0x0, 0x30(%r14)
movslq 0x18(%r14), %rax
testq %rax, %rax
jle 0xd65e
xorpd %xmm0, %xmm0
xorl %ecx, %ecx
xorps %xmm1, %xmm1
cvtsi2sd %ecx, %xmm1
mulsd (%rbx,%rcx,8), %xmm1
addsd %xmm1, %xmm0
movsd %xmm0, 0x28(%r14)
incq %rcx
cmpq %rcx, %rax
jne 0xd610
testl %eax, %eax
jle 0xd65e
xorpd %xmm1, %xmm1
xorl %ecx, %ecx
xorps %xmm2, %xmm2
cvtsi2sd %ecx, %xmm2
subsd %xmm0, %xmm2
mulsd %xmm2, %xmm2
mulsd (%rbx,%rcx,8), %xmm2
addsd %xmm2, %xmm1
movsd %xmm1, 0x30(%r14)
incq %rcx
cmpq %rcx, %rax
jne 0xd638
addq $0x8, %rsp
popq %rbx
popq %r14
retq
// File: /workspace/github_cpp_makefile_star2_O3/melop[P]admixem/./newran02/newran.cpp
// Function: DiscreteGen::DiscreteGen(int, double*, double*)
DiscreteGen::DiscreteGen(int n1, Real* prob, Real* val1)
                                                 // discrete generator
                                                 // values on *val
{
   #ifdef MONITOR
      cout << "constructing DiscreteGen\n";
   #endif
   Gen(n1, prob); val = new Real[n1];
   if (!val)  ErrorNoSpace();
   for (int i=0; i<n1; i++) val[i]=val1[i];
   mean=0.0; var=0.0;
   { for (int i=0; i<n; i++) mean = mean + val[i]*prob[i]; }
   { for (int i=0; i<n; i++) var = var + square(val[i]-mean)*prob[i]; }
}

// Assembly:
pushq %rbp
pushq %r15
pushq %r14
pushq %r13
pushq %r12
pushq %rbx
pushq %rax
movq %rcx, %r12
movq %rdx, %rbx
movl %esi, %ebp
movq %rdi, %r14
leaq 0x39883(%rip), %rax     # 0x47100
movq %rax, (%rdi)
callq 0xd666
movl %ebp, %r13d
shlq $0x3, %r13
testl %ebp, %ebp
movq $-0x1, %rdi
cmovnsq %r13, %rdi
callq 0x7030
movq %rax, %r15
movq %rax, 0x20(%r14)
testl %ebp, %ebp
jle 0xd8b7
movq %r15, %rdi
movq %r12, %rsi
movq %r13, %rdx
callq 0x72e0
xorpd %xmm0, %xmm0
movupd %xmm0, 0x28(%r14)
movslq 0x18(%r14), %rax
testq %rax, %rax
jle 0xd91e
leaq 0x28(%r14), %rcx
xorpd %xmm0, %xmm0
xorl %edx, %edx
movsd (%r15,%rdx,8), %xmm1
mulsd (%rbx,%rdx,8), %xmm1
addsd %xmm1, %xmm0
movsd %xmm0, (%rcx)
incq %rdx
cmpq %rdx, %rax
jne 0xd8d4
testl %eax, %eax
jle 0xd91e
xorpd %xmm1, %xmm1
xorl %ecx, %ecx
movsd (%r15,%rcx,8), %xmm2
subsd %xmm0, %xmm2
mulsd %xmm2, %xmm2
mulsd (%rbx,%rcx,8), %xmm2
addsd %xmm2, %xmm1
movsd %xmm1, 0x30(%r14)
incq %rcx
cmpq %rcx, %rax
jne 0xd8f9
addq $0x8, %rsp
popq %rbx
popq %r12
popq %r13
popq %r14
popq %r15
popq %rbp
retq
nop
// File: /workspace/github_cpp_makefile_star2_O3/melop[P]admixem/./newran02/newran.cpp
// Function: MultipliedRandom::Variance() const
ExtReal MultipliedRandom::Variance() const
{
   ExtReal e1 = square(rv1->Mean()); ExtReal e2 = square(rv2->Mean());
   ExtReal v1 = rv1->Variance(); ExtReal v2 = rv2->Variance();
   ExtReal r=v1*v2+v1*e2+e1*v2;
   return r;
}

// Assembly:
pushq %r15
pushq %r14
pushq %r12
pushq %rbx
subq $0x88, %rsp
movq %rdi, %r14
movq 0x8(%rdi), %rdi
movq (%rdi), %rax
callq *0x28(%rax)
leaq 0x18(%rsp), %r15
movsd %xmm0, (%r15)
movl %eax, 0x8(%r15)
movq %r15, %rdi
movq %r15, %rsi
callq 0xaf88
leaq 0x78(%rsp), %rbx
movsd %xmm0, (%rbx)
movl %eax, 0x8(%rbx)
movq 0x10(%r14), %rdi
movq (%rdi), %rax
callq *0x28(%rax)
leaq 0x8(%rsp), %r12
movsd %xmm0, (%r12)
movl %eax, 0x8(%r12)
movq %r12, %rdi
movq %r12, %rsi
callq 0xaf88
movsd %xmm0, (%r15)
movl %eax, 0x8(%r15)
movq 0x8(%r14), %rdi
movq (%rdi), %rax
callq *0x30(%rax)
movsd %xmm0, (%r12)
movl %eax, 0x8(%r12)
movq 0x10(%r14), %rdi
movq (%rdi), %rax
callq *0x30(%rax)
leaq 0x68(%rsp), %r14
movsd %xmm0, (%r14)
movl %eax, 0x8(%r14)
leaq 0x8(%rsp), %r15
movq %r15, %rdi
movq %r14, %rsi
callq 0xaf88
leaq 0x48(%rsp), %r12
movsd %xmm0, (%r12)
movl %eax, 0x8(%r12)
leaq 0x18(%rsp), %rsi
movq %r15, %rdi
callq 0xaf88
leaq 0x38(%rsp), %rsi
movsd %xmm0, (%rsi)
movl %eax, 0x8(%rsi)
movq %r12, %rdi
callq 0xae9c
leaq 0x58(%rsp), %r15
movsd %xmm0, (%r15)
movl %eax, 0x8(%r15)
movq %rbx, %rdi
movq %r14, %rsi
callq 0xaf88
leaq 0x28(%rsp), %rsi
movsd %xmm0, (%rsi)
movl %eax, 0x8(%rsi)
movq %r15, %rdi
callq 0xae9c
addq $0x88, %rsp
popq %rbx
popq %r12
popq %r14
popq %r15
retq
// File: /workspace/github_cpp_makefile_star2_O3/melop[P]admixem/./newran02/newran.cpp
// Function: MixedRandom::MixedRandom(int, double*, Random**)
MixedRandom::MixedRandom(int nx, Real* probs, Random** rvx)
{
   n = nx;
   rv = new Random*[n]; if (!rv) ErrorNoSpace();
   for (int i=0; i<n; i++) rv[i]=rvx[i];
   Build(probs);
}

// Assembly:
pushq %rbp
pushq %r15
pushq %r14
pushq %r12
pushq %rbx
movq %rcx, %r14
movq %rdx, %rbx
movl %esi, %ebp
movq %rdi, %r15
leaq 0x38de4(%rip), %rax     # 0x47160
movq %rax, (%rdi)
movl $0x4, %eax
movl %eax, 0x28(%rdi)
movl %eax, 0x38(%rdi)
movl %esi, 0x8(%rdi)
movl %esi, %r12d
leaq (,%r12,8), %rax
testl %esi, %esi
movq $-0x1, %rdi
cmovnsq %rax, %rdi
callq 0x7030
movq %rax, 0x18(%r15)
testl %ebp, %ebp
jle 0xe3c8
xorl %eax, %eax
movq (%r14,%rax,8), %rcx
movq 0x18(%r15), %rdx
movq %rcx, (%rdx,%rax,8)
incq %rax
cmpq %rax, %r12
jne 0xe3b4
movq %r15, %rdi
movq %rbx, %rsi
popq %rbx
popq %r12
popq %r14
popq %r15
popq %rbp
jmp 0xe3dc
nop
// File: /workspace/github_cpp_makefile_star2_O3/melop[P]admixem/./newran02/newran.cpp
// Function: MixedRandom::Build(double*)
void MixedRandom::Build(Real* probs)
{
   int i;
   dg=new DiscreteGen(n,probs);
   if (!dg) ErrorNoSpace();
   mean=0.0; var=0.0;
   for (i=0; i<n; i++) mean = mean + (rv[i])->Mean()*probs[i];
   for (i=0; i<n; i++)
   {
      ExtReal sigsq=(rv[i])->Variance();
      ExtReal mudif=(rv[i])->Mean()-mean;
      var = var + (sigsq+square(mudif))*probs[i];
   }

}

// Assembly:
pushq %rbp
pushq %r15
pushq %r14
pushq %r13
pushq %r12
pushq %rbx
subq $0x68, %rsp
movq %rsi, %rbx
movq %rdi, %r14
movl $0x38, %edi
callq 0x7380
movq %rax, %r15
movl 0x8(%r14), %esi
movq %rax, %rdi
movq %rbx, %rdx
callq 0xd5d6
movq %r15, 0x10(%r14)
xorl %eax, %eax
movq %rax, 0x20(%r14)
movl %eax, 0x28(%r14)
movq %rax, 0x30(%r14)
movl %eax, 0x38(%r14)
cmpl %eax, 0x8(%r14)
jle 0xe57f
leaq 0x20(%r14), %r15
leaq 0x30(%r14), %rax
movq %rax, 0x30(%rsp)
xorl %r13d, %r13d
movq %rsp, %rbp
leaq 0x20(%rsp), %r12
movq 0x18(%r14), %rax
movq (%rax,%r13,8), %rdi
movq (%rdi), %rax
callq *0x28(%rax)
movsd %xmm0, 0x10(%rsp)
movl %eax, 0x18(%rsp)
movsd (%rbx,%r13,8), %xmm0
movl $0x0, 0x8(%rsp)
movsd %xmm0, (%rsp)
leaq 0x10(%rsp), %rdi
movq %rbp, %rsi
callq 0xaf88
movsd %xmm0, 0x20(%rsp)
movl %eax, 0x28(%rsp)
movq %r15, %rdi
movq %r12, %rsi
callq 0xae9c
movsd %xmm0, 0x20(%r14)
movl %eax, 0x28(%r14)
incq %r13
movslq 0x8(%r14), %rax
cmpq %rax, %r13
jl 0xe444
testl %eax, %eax
jle 0xe57f
xorl %r12d, %r12d
movq %rsp, %r13
leaq 0x10(%rsp), %rbp
movq 0x18(%r14), %rax
movq (%rax,%r12,8), %rdi
movq (%rdi), %rax
callq *0x30(%rax)
movsd %xmm0, 0x20(%rsp)
movl %eax, 0x28(%rsp)
movq 0x18(%r14), %rax
movq (%rax,%r12,8), %rdi
movq (%rdi), %rax
callq *0x28(%rax)
movsd %xmm0, (%rsp)
movl %eax, 0x8(%rsp)
movq %r13, %rdi
movq %r15, %rsi
callq 0xaf12
movsd %xmm0, 0x10(%rsp)
movl %eax, 0x18(%rsp)
movq %rbp, %rdi
movq %rbp, %rsi
callq 0xaf88
movsd %xmm0, 0x48(%rsp)
movl %eax, 0x50(%rsp)
leaq 0x20(%rsp), %rdi
leaq 0x48(%rsp), %rsi
callq 0xae9c
movsd %xmm0, 0x58(%rsp)
movl %eax, 0x60(%rsp)
movsd (%rbx,%r12,8), %xmm0
movl $0x0, 0x40(%rsp)
movsd %xmm0, 0x38(%rsp)
leaq 0x58(%rsp), %rdi
leaq 0x38(%rsp), %rsi
callq 0xaf88
movsd %xmm0, (%rsp)
movl %eax, 0x8(%rsp)
movq 0x30(%rsp), %rdi
movq %r13, %rsi
callq 0xae9c
movsd %xmm0, 0x30(%r14)
movl %eax, 0x38(%r14)
incq %r12
movslq 0x8(%r14), %rax
cmpq %rax, %r12
jl 0xe4ba
addq $0x68, %rsp
popq %rbx
popq %r12
popq %r13
popq %r14
popq %r15
popq %rbp
retq
movq %rax, %rbx
movl $0x38, %esi
movq %r15, %rdi
callq 0x7390
movq %rbx, %rdi
callq 0x76a0
// File: /workspace/github_cpp_makefile_star2_O3/melop[P]admixem/./newran02/newran.cpp
// Function: MixedRandom::MixedRandom(AddedSelectedRandom&)
MixedRandom::MixedRandom(AddedSelectedRandom& sr)
{
   n = sr.nelems();                       // number of terms;
   Real* probs = new Real[n]; rv = new Random*[n];
   if (!probs || !rv) ErrorNoSpace();
   int i=0; sr.load(&i,probs,rv);
   Build(probs); delete [] probs;
}

// Assembly:
pushq %r15
pushq %r14
pushq %r12
pushq %rbx
pushq %rax
movq %rsi, %r14
movq %rdi, %rbx
leaq 0x38ba5(%rip), %rax     # 0x47160
movq %rax, (%rdi)
movl $0x4, %eax
movl %eax, 0x28(%rdi)
movl %eax, 0x38(%rdi)
movq (%rsi), %rax
movq %rsi, %rdi
callq *0x40(%rax)
movl %eax, 0x8(%rbx)
movl %eax, %ecx
shlq $0x3, %rcx
testl %eax, %eax
movq $-0x1, %r15
cmovnsq %rcx, %r15
movq %r15, %rdi
callq 0x7030
movq %rax, %r12
movq %r15, %rdi
callq 0x7030
movq %rax, 0x18(%rbx)
leaq 0x4(%rsp), %rsi
movl $0x0, (%rsi)
movq (%r14), %r8
movq %r14, %rdi
movq %r12, %rdx
movq %rax, %rcx
callq *0x48(%r8)
movq %rbx, %rdi
movq %r12, %rsi
callq 0xe3dc
movq %r12, %rdi
callq 0x7490
addq $0x8, %rsp
popq %rbx
popq %r12
popq %r14
popq %r15
retq
nop
// File: /workspace/github_cpp_makefile_star2_O3/melop[P]admixem/./newran02/newran.cpp
// Function: RandomPermutation::Next(int, int, int*, int)
void RandomPermutation::Next(int N, int M, int p[], int start)
{
   // N = size of urn; M = number of draws
   if (N < M) Throw(Logic_error("Newran: N < M in RandomPermutation"));
   int i;
   int* q = new int [N];
   if (!q) ErrorNoSpace();
   for (i = 0; i < N; i++) q[i] = start + i;
   for (i = 0; i < M; i++)
   {
      int k = i + (int)(U.Next() * (N - i));       // uniform on i ... N-1
      p[i] = q[k]; q[k] = q[i];                    // swap i and k terms
                                                   // but put i-th term into p
   }
   delete [] q;
}

// Assembly:
pushq %rbp
pushq %r15
pushq %r14
pushq %r13
pushq %r12
pushq %rbx
pushq %rax
cmpl %edx, %esi
jl 0xea0f
movl %r8d, %r12d
movq %rcx, %rbx
movl %edx, %r15d
movl %esi, %ebp
movl %esi, %r13d
leaq (,%r13,4), %rax
testl %esi, %esi
movq $-0x1, %rdi
cmovnsq %rax, %rdi
callq 0x7030
movq %rax, %r14
testl %ebp, %ebp
jle 0xe9bb
leal 0x1(%r13), %eax
decq %r13
movq %r13, %xmm0
pshufd $0x44, %xmm0, %xmm0     # xmm0 = xmm0[0,1,0,1]
shrl %eax
shlq $0x3, %rax
movdqa 0x25b8d(%rip), %xmm1    # 0x344d0
xorl %ecx, %ecx
movdqa 0x25b93(%rip), %xmm2    # 0x344e0
pxor %xmm2, %xmm0
pcmpeqd %xmm3, %xmm3
movdqa 0x25b93(%rip), %xmm4    # 0x344f0
movdqa %xmm1, %xmm5
pxor %xmm2, %xmm5
movdqa %xmm5, %xmm6
pcmpgtd %xmm0, %xmm6
pcmpeqd %xmm0, %xmm5
pshufd $0xf5, %xmm5, %xmm7     # xmm7 = xmm5[1,1,3,3]
pand %xmm6, %xmm7
pshufd $0xf5, %xmm6, %xmm5     # xmm5 = xmm6[1,1,3,3]
por %xmm7, %xmm5
movd %xmm5, %edx
notl %edx
testb $0x1, %dl
je 0xe992
movl %r12d, (%r14,%rcx)
pxor %xmm3, %xmm5
pextrw $0x4, %xmm5, %edx
testb $0x1, %dl
je 0xe9aa
leal 0x1(%r12), %edx
movl %edx, 0x4(%r14,%rcx)
paddq %xmm4, %xmm1
addq $0x8, %rcx
addl $0x2, %r12d
cmpq %rcx, %rax
jne 0xe95d
testl %r15d, %r15d
jle 0xe9f9
movl %r15d, %r15d
xorl %r12d, %r12d
callq 0xbba0
xorps %xmm1, %xmm1
cvtsi2sd %ebp, %xmm1
mulsd %xmm0, %xmm1
cvttsd2si %xmm1, %eax
cltq
addq %r12, %rax
movl (%r14,%rax,4), %ecx
movl %ecx, (%rbx,%r12,4)
movl (%r14,%r12,4), %ecx
movl %ecx, (%r14,%rax,4)
incq %r12
decl %ebp
cmpq %r12, %r15
jne 0xe9c6
movq %r14, %rdi
addq $0x8, %rsp
popq %rbx
popq %r12
popq %r13
popq %r14
popq %r15
popq %rbp
jmp 0x7490
movl $0x1, %edi
callq 0x7180
movq %rax, %rbx
leaq 0x25f61(%rip), %rsi     # 0x34984
movq %rax, %rdi
callq 0xb8ea
leaq 0x37f86(%rip), %rsi     # 0x469b8
movq %rbx, %rdi
xorl %edx, %edx
callq 0x7660
movq %rax, %r14
movq %rbx, %rdi
callq 0x72a0
movq %r14, %rdi
callq 0x76a0
// File: /workspace/github_cpp_makefile_star2_O3/melop[P]admixem/./newran02/newran.cpp
// Function: ShellSortAscending(int*, int)
void ShellSortAscending(int* a, int N)
{
   int h;
   for (h = 1; h <= N / 9; h = 3*h + 1) {}
   for (; h > 0; h /= 3) for (int i = h; i < N; ++i)
   {
      int v = a[i]; int j = i;
      while (j>=h && a[j-h] > v) { a[j] = a[j-h]; j-= h; }
      a[j] = v;
   }
}

// Assembly:
movl $0x1, %eax
cmpl $0x9, %esi
jl 0xea89
movslq %esi, %rax
imulq $0x38e38e39, %rax, %rcx # imm = 0x38E38E39
movq %rcx, %rax
shrq $0x3f, %rax
sarq $0x21, %rcx
addl %eax, %ecx
movl $0x1, %eax
movl %eax, %edx
leal (%rdx,%rdx,2), %r8d
leal (%rdx,%rdx,2), %eax
incl %eax
cmpl %ecx, %r8d
jl 0xea75
testl %edx, %edx
js 0xeaf2
pushq %rbp
pushq %r14
pushq %rbx
movl $0xaaaaaaab, %ecx       # imm = 0xAAAAAAAB
movl %eax, %edx
cmpl %esi, %eax
jge 0xeadf
xorl %r8d, %r8d
movq %rdx, %r9
movl (%rdi,%r9,4), %r10d
movl %r9d, %ebx
cmpq %rdx, %r9
jb 0xeacd
movl %r8d, %r11d
leal (%rax,%r11), %ebx
movslq %r11d, %r14
movl (%rdi,%r14,4), %ebp
cmpl %r10d, %ebp
jle 0xeacd
movl %ebx, %ebx
movl %ebp, (%rdi,%rbx,4)
subl %eax, %r11d
jge 0xeaad
addl %eax, %r11d
movl %r11d, %ebx
movslq %ebx, %r11
movl %r10d, (%rdi,%r11,4)
incq %r9
incl %r8d
cmpl %esi, %r9d
jl 0xea9e
imulq %rcx, %rdx
shrq $0x21, %rdx
cmpl $0x3, %eax
movl %edx, %eax
jae 0xea92
popq %rbx
popq %r14
popq %rbp
retq
nop
// File: /workspace/github_cpp_makefile_star2_O3/melop[P]admixem/./newran02/newran.cpp
// Function: VariPoisson::iNext_very_small(double)
int VariPoisson::iNext_very_small(Real mu)
{
   // mu is expected value of Poisson random number
   // Efficient when mu is small
   // generate a Poisson variable with mean mu
   if (mu == 0) return 0;
   Real t = exp(-mu); int k = 0; Real u = U.Next();
   for (Real s = t; s < u; s += t) { ++k; t *= mu / k; }
   return k;
}

// Assembly:
pushq %rbx
subq $0x20, %rsp
movapd %xmm0, %xmm1
xorl %ebx, %ebx
xorpd %xmm0, %xmm0
ucomisd %xmm0, %xmm1
jne 0xebfd
jnp 0xec56
movapd 0x2548b(%rip), %xmm0    # 0x34090
xorpd %xmm1, %xmm0
movapd %xmm1, 0x10(%rsp)
callq 0x7748
movsd %xmm0, 0x8(%rsp)
callq 0xbba0
movsd 0x8(%rsp), %xmm5
movapd 0x10(%rsp), %xmm4
ucomisd %xmm5, %xmm0
jbe 0xec56
xorl %ebx, %ebx
movapd %xmm5, %xmm1
incl %ebx
xorps %xmm2, %xmm2
cvtsi2sd %ebx, %xmm2
movapd %xmm4, %xmm3
divsd %xmm2, %xmm3
mulsd %xmm3, %xmm1
addsd %xmm1, %xmm5
ucomisd %xmm5, %xmm0
ja 0xec37
movl %ebx, %eax
addq $0x20, %rsp
popq %rbx
retq
// File: /workspace/github_cpp_makefile_star2_O3/melop[P]admixem/./newran02/newran.cpp
// Function: VariPoisson::iNext_small(double)
int VariPoisson::iNext_small(Real mu)
{
   // mu is expected value of Poisson random number
   // Efficient when mu is reasonably small
   // generate a Poisson variable with mean mu
   // start iteration at mode
   if (mu < 10) return iNext_very_small(mu);
   int k_start = (int)mu;  Real u = U.Next();
   Real t1 = exp(k_start * log(mu) - mu - ln_gamma(k_start+1));
   if (t1 > u) return k_start;
   int k1 = k_start; int k2 = k_start; Real t2 = t1; Real s = t1;
   for(;;)
   {
      ++k1; t1 *= mu / k1; s += t1; if (s > u) return k1;
      if (k2 > 0) { t2 *= k2 / mu; --k2; s += t2; if (s > u) return k2; }
   }
}

// Assembly:
movsd 0x25812(%rip), %xmm1    # 0x34478
ucomisd %xmm0, %xmm1
ja 0xebe6
pushq %rbp
pushq %rbx
subq $0x38, %rsp
cvttsd2si %xmm0, %ebx
movapd %xmm0, 0x10(%rsp)
callq 0xbba0
movsd %xmm0, 0x30(%rsp)
movapd 0x10(%rsp), %xmm0
cvttpd2dq %xmm0, %xmm1
cvtdq2pd %xmm1, %xmm1
movaps %xmm1, (%rsp)
callq 0x7730
mulsd (%rsp), %xmm0
subsd 0x10(%rsp), %xmm0
movsd %xmm0, (%rsp)
leal 0x1(%rbx), %ebp
xorps %xmm0, %xmm0
cvtsi2sd %ebp, %xmm0
movsd %xmm0, 0x28(%rsp)
callq 0xcb6e
movsd (%rsp), %xmm1
subsd %xmm0, %xmm1
movapd %xmm1, %xmm0
callq 0x7748
movsd 0x30(%rsp), %xmm6
ucomisd %xmm6, %xmm0
jbe 0xece9
movl %ebx, %eax
jmp 0xed51
movapd 0x10(%rsp), %xmm5
movapd %xmm5, %xmm1
divsd 0x28(%rsp), %xmm1
mulsd %xmm0, %xmm1
movapd %xmm0, %xmm2
addsd %xmm1, %xmm2
ucomisd %xmm6, %xmm2
jbe 0xed0f
movl %ebp, %eax
jmp 0xed51
testl %ebx, %ebx
jle 0xed30
xorps %xmm3, %xmm3
cvtsi2sd %ebx, %xmm3
divsd %xmm5, %xmm3
mulsd %xmm3, %xmm0
decl %ebx
addsd %xmm0, %xmm2
ucomisd %xmm6, %xmm2
movl %ebx, %eax
ja 0xed51
incl %ebp
xorps %xmm3, %xmm3
cvtsi2sd %ebp, %xmm3
movapd %xmm5, %xmm4
divsd %xmm3, %xmm4
mulsd %xmm4, %xmm1
addsd %xmm1, %xmm2
ucomisd %xmm6, %xmm2
movl %ebp, %eax
jbe 0xed0f
addq $0x38, %rsp
popq %rbx
popq %rbp
retq
// File: /workspace/github_cpp_makefile_star2_O3/melop[P]admixem/./newran02/newran.cpp
// Function: VariPoisson::iNext_large(double)
int VariPoisson::iNext_large(Real mu)
{
   // mu is expected value of Poisson random number
   // reasonably accurate when mu is large, but should try to improve
   // generate a Poisson variable with mean mu
   // method is to start with normal variable X and use X with prob 1/3
   // and X**2 with prob 2/3. In each case X has mean and variance to
   // agree with Poisson rv after adjusting with Sheppard's correction

   const Real sc = 1.0 / 12.0;        // Sheppard correction
   int k;
   if (U.Next() > 1.0 / 3.0)
   {
      Real musc = 0.5 * (mu - sc); Real meansq = sqrt(mu * mu - musc);
      Real sigma = sqrt(musc / (mu + meansq));
      k = (int)floor(square(N.Next() * sigma + sqrt(meansq)) + 0.5);
   }
   else k = (int)floor(N.Next() * sqrt(mu - sc) + mu + 0.5);
   return k < 0 ? 0 : k;
}

// Assembly:
pushq %rbx
subq $0x20, %rsp
movsd %xmm0, 0x8(%rsp)
movq %rdi, %rbx
callq 0xbba0
ucomisd 0x2570d(%rip), %xmm0    # 0x34480
jbe 0xedab
movsd 0x2570b(%rip), %xmm2    # 0x34488
movsd 0x8(%rsp), %xmm1
addsd %xmm1, %xmm2
mulsd 0x25651(%rip), %xmm2    # 0x343e0
movapd %xmm1, %xmm0
mulsd %xmm1, %xmm0
subsd %xmm2, %xmm0
xorpd %xmm3, %xmm3
ucomisd %xmm3, %xmm0
jb 0xede1
sqrtsd %xmm0, %xmm0
jmp 0xedfc
addq $0x8, %rbx
movq %rbx, %rdi
callq 0xc03e
movapd %xmm0, %xmm2
movsd 0x256c5(%rip), %xmm0    # 0x34488
movsd 0x8(%rsp), %xmm3
addsd %xmm3, %xmm0
xorpd %xmm1, %xmm1
ucomisd %xmm1, %xmm0
jb 0xee19
movapd %xmm3, %xmm1
sqrtsd %xmm0, %xmm0
jmp 0xee30
movsd %xmm2, 0x10(%rsp)
callq 0x7728
xorpd %xmm3, %xmm3
movsd 0x10(%rsp), %xmm2
movsd 0x8(%rsp), %xmm1
movsd %xmm0, 0x8(%rsp)
addsd %xmm0, %xmm1
divsd %xmm1, %xmm2
ucomisd %xmm3, %xmm2
jb 0xee3a
xorps %xmm0, %xmm0
sqrtsd %xmm2, %xmm0
jmp 0xee43
movsd %xmm2, 0x10(%rsp)
callq 0x7728
movsd 0x10(%rsp), %xmm2
movsd 0x8(%rsp), %xmm1
mulsd %xmm0, %xmm2
addsd %xmm1, %xmm2
jmp 0xee92
movapd %xmm2, %xmm0
callq 0x7728
movsd %xmm0, 0x18(%rsp)
addq $0x8, %rbx
movq %rbx, %rdi
callq 0xc03e
movapd %xmm0, %xmm2
xorpd %xmm0, %xmm0
movsd 0x8(%rsp), %xmm1
ucomisd %xmm0, %xmm1
movapd %xmm1, %xmm0
jb 0xee73
sqrtsd %xmm0, %xmm0
jmp 0xee84
movsd %xmm2, 0x10(%rsp)
callq 0x7728
movsd 0x10(%rsp), %xmm2
mulsd 0x18(%rsp), %xmm2
addsd %xmm0, %xmm2
mulsd %xmm2, %xmm2
addsd 0x25546(%rip), %xmm2    # 0x343e0
movapd %xmm2, %xmm0
callq 0x7740
cvttsd2si %xmm0, %ecx
xorl %eax, %eax
testl %ecx, %ecx
cmovgl %ecx, %eax
addq $0x20, %rsp
popq %rbx
retq
// File: /workspace/github_cpp_makefile_star2_O3/melop[P]admixem/./newran02/newran.cpp
// Function: VariPoisson::iNext(double)
int VariPoisson::iNext(Real mu)
{
   if (mu <= 0)
   {
      if (mu == 0) return 0;
      else Throw(Logic_error("Newran: illegal parameters"));
   }
   if (mu < 10) return iNext_very_small(mu);
   else if (mu < 100) return iNext_small(mu);
   else if (mu < 200)
   {
      // do in two statements so order of evaluation is preserved
      int i = (int)P100.Next();
      return i + iNext_small(mu - 100);
   }
   else if (mu < 300)
   {
      // do in two statements so order of evaluation is preserved
      int i = (int)P200.Next();
      return i + iNext_small(mu - 200);
   }
   else return iNext_large(mu);
}

// Assembly:
pushq %r14
pushq %rbx
pushq %rax
xorpd %xmm1, %xmm1
ucomisd %xmm0, %xmm1
jae 0xeedc
movsd 0x255ae(%rip), %xmm1    # 0x34478
ucomisd %xmm0, %xmm1
jbe 0xeef6
addq $0x8, %rsp
popq %rbx
popq %r14
jmp 0xebe6
xorl %eax, %eax
ucomisd %xmm1, %xmm0
jne 0xef84
jp 0xef84
addq $0x8, %rsp
popq %rbx
popq %r14
retq
movsd 0x25112(%rip), %xmm1    # 0x34010
ucomisd %xmm0, %xmm1
jbe 0xef10
addq $0x8, %rsp
popq %rbx
popq %r14
jmp 0xec5e
movsd 0x25578(%rip), %xmm1    # 0x34490
ucomisd %xmm0, %xmm1
jbe 0xef3f
movq 0x38(%rdi), %rdi
movq (%rdi), %rax
movsd %xmm0, (%rsp)
callq *(%rax)
cvttsd2si %xmm0, %ebx
movsd (%rsp), %xmm0
addsd 0x2556b(%rip), %xmm0    # 0x344a8
jmp 0xef6c
movsd 0x25551(%rip), %xmm1    # 0x34498
ucomisd %xmm0, %xmm1
jbe 0xef78
movq 0x48(%rdi), %rdi
movq (%rdi), %rax
movsd %xmm0, (%rsp)
callq *(%rax)
cvttsd2si %xmm0, %ebx
movsd (%rsp), %xmm0
addsd 0x25534(%rip), %xmm0    # 0x344a0
callq 0xec5e
addl %ebx, %eax
jmp 0xeeee
addq $0x8, %rsp
popq %rbx
popq %r14
jmp 0xed58
movl $0x1, %edi
callq 0x7180
movq %rax, %rbx
leaq 0x25883(%rip), %rsi     # 0x3481b
movq %rax, %rdi
callq 0xb8ea
leaq 0x37a11(%rip), %rsi     # 0x469b8
movq %rbx, %rdi
xorl %edx, %edx
callq 0x7660
movq %rax, %r14
movq %rbx, %rdi
callq 0x72a0
movq %r14, %rdi
callq 0x76a0
// File: /workspace/github_cpp_makefile_star2_O3/melop[P]admixem/./newran02/newran.cpp
// Function: VariBinomial::iNext_very_small(int, double)
int VariBinomial::iNext_very_small(int n, Real p)
{
   // Efficient when n * p is small
   // generate a Binomial variable with parameters n and p

   Real q = 1 - p; if (q == 0) return n;
   Real r = p / q; Real t = pow(q, n); Real s = t;
   Real u = U.Next();
   for (int k = 0; k <= n; ++k)
   {
      if (s >= u) return k;
      t *= r * (Real)(n - k) / (Real)(k + 1);
      s += t;
   }
   return 0;    // can happen only if we have round-off error
}

// Assembly:
pushq %rbx
subq $0x20, %rsp
movl %esi, %ebx
movsd 0x25031(%rip), %xmm2    # 0x34018
subsd %xmm0, %xmm2
xorpd %xmm1, %xmm1
ucomisd %xmm1, %xmm2
jne 0xeffe
jp 0xeffe
movl %ebx, %eax
jmp 0xf081
movsd %xmm0, 0x10(%rsp)
xorps %xmm1, %xmm1
cvtsi2sd %ebx, %xmm1
movsd %xmm2, 0x8(%rsp)
movapd %xmm2, %xmm0
callq 0x7690
movsd %xmm0, 0x18(%rsp)
callq 0xbba0
xorl %eax, %eax
testl %ebx, %ebx
js 0xf081
movsd 0x18(%rsp), %xmm5
ucomisd %xmm0, %xmm5
jae 0xf081
movsd 0x10(%rsp), %xmm4
divsd 0x8(%rsp), %xmm4
xorl %ecx, %ecx
movl %ebx, %eax
movapd %xmm5, %xmm1
movl %ecx, %edx
leal 0x1(%rdx), %ecx
cmpl %ebx, %edx
jge 0xf07a
xorps %xmm2, %xmm2
cvtsi2sd %eax, %xmm2
mulsd %xmm4, %xmm2
xorps %xmm3, %xmm3
cvtsi2sd %ecx, %xmm3
divsd %xmm3, %xmm2
mulsd %xmm2, %xmm1
addsd %xmm1, %xmm5
decl %eax
ucomisd %xmm0, %xmm5
jb 0xf04b
xorl %eax, %eax
cmpl %ebx, %edx
cmovll %ecx, %eax
addq $0x20, %rsp
popq %rbx
retq
nop
// File: /workspace/github_cpp_makefile_star2_O3/melop[P]admixem/./newran02/newran.cpp
// Function: VariBinomial::iNext(int, double)
int VariBinomial::iNext(int n, Real p)
{
   if (p > 0.5) return n - iNext(n, 1.0 - p);

   if (n <= 0 || p <= 0)
   {
      if (n < 0 || p < 0) Throw(Logic_error("Newran: illegal parameters"));
      else return 0;
   }

   Real mean = n * p;
   if (mean <= 10) return iNext_very_small(n, p);
   else if (mean <= 200) return iNext_small(n, p);
   else return iNext_large(n,p);

}

// Assembly:
pushq %rbp
pushq %r14
pushq %rbx
subq $0x10, %rsp
movl %esi, %ebp
movq %rdi, %rbx
ucomisd 0x25167(%rip), %xmm0    # 0x343e0
jbe 0xf29c
movsd 0x24d95(%rip), %xmm1    # 0x34018
subsd %xmm0, %xmm1
movq %rbx, %rdi
movl %ebp, %esi
movapd %xmm1, %xmm0
callq 0xf264
subl %eax, %ebp
jmp 0xf39c
testl %ebp, %ebp
jle 0xf307
xorpd %xmm1, %xmm1
ucomisd %xmm0, %xmm1
jae 0xf307
cvtsi2sd %ebp, %xmm2
mulsd %xmm0, %xmm2
movsd 0x251be(%rip), %xmm1    # 0x34478
ucomisd %xmm2, %xmm1
jae 0xf348
movsd 0x251c4(%rip), %xmm1    # 0x34490
ucomisd %xmm2, %xmm1
jae 0xf357
movsd 0x24d3a(%rip), %xmm1    # 0x34018
subsd %xmm0, %xmm1
movsd %xmm2, 0x8(%rsp)
mulsd %xmm2, %xmm1
addsd 0x25194(%rip), %xmm1    # 0x34488
xorpd %xmm0, %xmm0
ucomisd %xmm0, %xmm1
jb 0xf366
xorps %xmm0, %xmm0
sqrtsd %xmm1, %xmm0
jmp 0xf36f
testl %ebp, %ebp
js 0xf31b
xorl %ebp, %ebp
xorpd %xmm1, %xmm1
ucomisd %xmm0, %xmm1
jbe 0xf39c
movl $0x1, %edi
callq 0x7180
movq %rax, %rbx
leaq 0x254ec(%rip), %rsi     # 0x3481b
movq %rax, %rdi
callq 0xb8ea
leaq 0x3767a(%rip), %rsi     # 0x469b8
movq %rbx, %rdi
xorl %edx, %edx
callq 0x7660
movl %ebp, %esi
addq $0x10, %rsp
popq %rbx
popq %r14
popq %rbp
jmp 0xefd8
movl %ebp, %esi
addq $0x10, %rsp
popq %rbx
popq %r14
popq %rbp
jmp 0xf088
movapd %xmm1, %xmm0
callq 0x7728
movsd %xmm0, (%rsp)
addq $0x8, %rbx
movq %rbx, %rdi
callq 0xc03e
mulsd (%rsp), %xmm0
addsd 0x8(%rsp), %xmm0
addsd 0x2504d(%rip), %xmm0    # 0x343e0
callq 0x7740
cvttsd2si %xmm0, %ebp
movl %ebp, %eax
addq $0x10, %rsp
popq %rbx
popq %r14
popq %rbp
retq
movq %rax, %r14
movq %rbx, %rdi
callq 0x72a0
movq %r14, %rdi
callq 0x76a0
// File: /workspace/github_cpp_makefile_star2_O3/melop[P]admixem/./newran02/newran.cpp
// Function: VariBinomial::iNext_large(int, double)
int VariBinomial::iNext_large(int n, Real p)
{
   const Real sc = 1.0 / 12.0;        // Sheppard correction
   Real mean = n * p;
   Real sd = sqrt(mean * (1.0 - p) - sc);
   return (int)floor(N.Next() * sd + mean + 0.5);
}

// Assembly:
pushq %rbx
subq $0x10, %rsp
movapd %xmm0, %xmm1
cvtsi2sd %esi, %xmm2
mulsd %xmm0, %xmm2
movsd 0x24c45(%rip), %xmm0    # 0x34018
subsd %xmm1, %xmm0
movsd %xmm2, 0x8(%rsp)
mulsd %xmm2, %xmm0
addsd 0x2509f(%rip), %xmm0    # 0x34488
movq %rdi, %rbx
xorpd %xmm1, %xmm1
ucomisd %xmm1, %xmm0
jb 0xf3fc
sqrtsd %xmm0, %xmm0
jmp 0xf401
callq 0x7728
movsd %xmm0, (%rsp)
addq $0x8, %rbx
movq %rbx, %rdi
callq 0xc03e
mulsd (%rsp), %xmm0
addsd 0x8(%rsp), %xmm0
addsd 0x24fbb(%rip), %xmm0    # 0x343e0
callq 0x7740
cvttsd2si %xmm0, %eax
addq $0x10, %rsp
popq %rbx
retq
// File: /workspace/github_cpp_makefile_star2_O3/melop[P]admixem/config.cpp
// Function: SimulationConfigurations::GetConfig(std::__cxx11::basic_string<char, std::char_traits<char>, std::allocator<char>>)
const string SimulationConfigurations::GetConfig(string sKey) {

	//printf("Size of map: %d\n", this->_mpConfigs.size());

	

	if (  this->_mpConfigs.find(sKey) == this ->_mpConfigs.end()) { // the requested key is not found in the config

		//printf("The key %s was not found in the configuration file.\n", sKey.c_str());
		return "";

	}
	

	return  this->_mpConfigs[sKey];


}

// Assembly:
pushq %r15
pushq %r14
pushq %r12
pushq %rbx
pushq %rax
movq %rdx, %r14
movq %rsi, %r12
movq %rdi, %rbx
leaq 0x58(%rsi), %r15
movq %r15, %rdi
movq %rdx, %rsi
callq 0x21cca
addq $0x60, %r12
cmpq %r12, %rax
je 0x1028d
movq %r15, %rdi
movq %r14, %rsi
callq 0x19f5e
leaq 0x10(%rbx), %rcx
movq %rcx, (%rbx)
movq (%rax), %rsi
movq 0x8(%rax), %rdx
addq %rsi, %rdx
movq %rbx, %rdi
callq 0xa9b2
jmp 0x102a6
leaq 0x10(%rbx), %rax
movq %rax, (%rbx)
leaq 0x24d9e(%rip), %rdx     # 0x35039
movq %rbx, %rdi
movq %rdx, %rsi
callq 0xa91c
movq %rbx, %rax
addq $0x8, %rsp
popq %rbx
popq %r12
popq %r14
popq %r15
retq
nop