// File: Kingsford-Group[P]ghost2/src/gzReader.hpp
// Function: ipsenMikhailovVector(double*, int, double)
vector<double> ipsenMikhailovVector(double* spectrum, int s, double gamma)
{
  double c = 0;
  vector<double> dens;
  double ub = gamma>.1?2.01:(int)((2 + .1*gamma)*100)/100.0;
  double step = gamma<.1?.01:(int)(10*gamma)/100.0;
  double sigsq = 2.0*gamma*gamma;
  double norm = 1.0/sqrt(PI * sigsq);
  for(double i = 0; i < ub; i += step){
    double d = gaussianDensity(spectrum, s, i, sigsq, norm);
    c += d;
    dens.push_back(d);
  }
  for(int i=0; i < dens.size(); i++) dens[i] /= c;
  return dens;
}

// Assembly:
pushq %rbp
movq %rsp, %rbp
subq $0xb0, %rsp
movq %rdi, -0x80(%rbp)
movq %rdi, %rax
movq %rax, -0x78(%rbp)
movq %rdi, -0x8(%rbp)
movq %rsi, -0x10(%rbp)
movl %edx, -0x14(%rbp)
movsd %xmm0, -0x20(%rbp)
xorps %xmm0, %xmm0
movsd %xmm0, -0x28(%rbp)
movb $0x0, -0x29(%rbp)
callq 0x251f0
movsd -0x20(%rbp), %xmm0
movsd 0xa494c(%rip), %xmm1    # 0xbc090
ucomisd %xmm1, %xmm0
jbe 0x1775c
movsd 0xa494e(%rip), %xmm0    # 0xbc0a0
movsd %xmm0, -0x88(%rbp)
jmp 0x177a1
movsd -0x20(%rbp), %xmm0
movsd 0xa4927(%rip), %xmm1    # 0xbc090
mulsd %xmm1, %xmm0
movsd 0xa4923(%rip), %xmm1    # 0xbc098
addsd %xmm1, %xmm0
movsd 0xa48d7(%rip), %xmm1    # 0xbc058
mulsd %xmm1, %xmm0
cvttsd2si %xmm0, %eax
cvtsi2sd %eax, %xmm0
movsd 0xa48c3(%rip), %xmm1    # 0xbc058
divsd %xmm1, %xmm0
movsd %xmm0, -0x88(%rbp)
movsd -0x88(%rbp), %xmm0
movsd %xmm0, -0x38(%rbp)
movsd 0xa48da(%rip), %xmm0    # 0xbc090
ucomisd -0x20(%rbp), %xmm0
jbe 0x177cf
movsd 0xa48a3(%rip), %xmm0    # 0xbc068
movsd %xmm0, -0x90(%rbp)
jmp 0x177f8
movsd 0xa48d1(%rip), %xmm0    # 0xbc0a8
mulsd -0x20(%rbp), %xmm0
cvttsd2si %xmm0, %eax
cvtsi2sd %eax, %xmm0
movsd 0xa486c(%rip), %xmm1    # 0xbc058
divsd %xmm1, %xmm0
movsd %xmm0, -0x90(%rbp)
movsd -0x90(%rbp), %xmm0
movsd %xmm0, -0x40(%rbp)
movsd 0xa488b(%rip), %xmm0    # 0xbc098
mulsd -0x20(%rbp), %xmm0
mulsd -0x20(%rbp), %xmm0
movsd %xmm0, -0x48(%rbp)
movsd 0xa488c(%rip), %xmm0    # 0xbc0b0
mulsd -0x48(%rbp), %xmm0
callq 0x10560
movaps %xmm0, %xmm1
movsd 0xa4837(%rip), %xmm0    # 0xbc070
divsd %xmm1, %xmm0
movsd %xmm0, -0x50(%rbp)
xorps %xmm0, %xmm0
movsd %xmm0, -0x58(%rbp)
movsd -0x58(%rbp), %xmm1
movsd -0x38(%rbp), %xmm0
ucomisd %xmm1, %xmm0
jbe 0x178df
movq -0x10(%rbp), %rdi
movl -0x14(%rbp), %esi
movsd -0x58(%rbp), %xmm0
movsd -0x48(%rbp), %xmm1
movsd -0x50(%rbp), %xmm2
callq 0x17530
movsd %xmm0, -0x98(%rbp)
jmp 0x17883
movq -0x80(%rbp), %rdi
movsd -0x98(%rbp), %xmm0
movsd %xmm0, -0x60(%rbp)
movsd -0x60(%rbp), %xmm1
movsd -0x28(%rbp), %xmm0
addsd %xmm1, %xmm0
movsd %xmm0, -0x28(%rbp)
leaq -0x60(%rbp), %rsi
callq 0x253d0
jmp 0x178b2
jmp 0x178b4
movsd -0x40(%rbp), %xmm0
addsd -0x58(%rbp), %xmm0
movsd %xmm0, -0x58(%rbp)
jmp 0x1784a
movq -0x80(%rbp), %rdi
movq %rax, %rcx
movl %edx, %eax
movq %rcx, -0x68(%rbp)
movl %eax, -0x6c(%rbp)
callq 0x25480
jmp 0x17962
movl $0x0, -0x70(%rbp)
movq -0x80(%rbp), %rdi
movslq -0x70(%rbp), %rax
movq %rax, -0xa0(%rbp)
callq 0x254e0
movq %rax, %rcx
movq -0xa0(%rbp), %rax
cmpq %rcx, %rax
jae 0x17942
movq -0x80(%rbp), %rdi
movsd -0x28(%rbp), %xmm0
movsd %xmm0, -0xa8(%rbp)
movslq -0x70(%rbp), %rsi
callq 0x25500
movsd -0xa8(%rbp), %xmm1
movsd (%rax), %xmm0
divsd %xmm1, %xmm0
movsd %xmm0, (%rax)
movl -0x70(%rbp), %eax
addl $0x1, %eax
movl %eax, -0x70(%rbp)
jmp 0x178e6
movb $0x1, -0x29(%rbp)
testb $0x1, -0x29(%rbp)
jne 0x17955
movq -0x80(%rbp), %rdi
callq 0x25480
movq -0x78(%rbp), %rax
addq $0xb0, %rsp
popq %rbp
retq
movq -0x68(%rbp), %rdi
callq 0x10cc0
nopl (%rax,%rax)