user: What would the C source look like for the assembly code provided here?

```asm
<Rates4Sites>:
  endbr64
  push   %r15
  movapd %xmm0,%xmm1
  push   %r14
  mov    %r8,%r14
  push   %r13
  mov    %esi,%r13d
  push   %r12
  mov    %edx,%r12d
  push   %rbp
  mov    %rdi,%rbp
  push   %rbx
  mov    %ecx,%ebx
  sub    $0x18,%rsp
  ucomisd D11(%rip),%xmm0
  jp     L94
  jne    L94
  test   %rdi,%rdi
  je     L101
  test   %edx,%edx
  jle    L101
  cmp    $0x1,%edx
  je     L187
  shr    %edx
  movapd D14(%rip),%xmm0
  mov    %rdi,%rax
  sub    $0x1,%edx
  shl    $0x4,%rdx
  lea    0x10(%rdi,%rdx,1),%rdx
  nopl   0x0(%rax,%rax,1)
L169:
  movups %xmm0,(%rax)
  add    $0x10,%rax
  cmp    %rdx,%rax
  jne    L169
  mov    %r12d,%eax
  and    $0xfffffffe,%eax
  and    $0x1,%r12d
  je     L101
L168:
  mov    D14(%rip),%rbx
  cltq
  mov    %rbx,0x0(%rbp,%rax,8)
L101:
  add    $0x18,%rsp
  xor    %eax,%eax
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
  nopw   0x0(%rax,%rax,1)
L94:
  cmp    $0x1,%r13d
  jg     L114
  test   %r12d,%r12d
  jle    L37
  lea    -0x1(%r12),%eax
  mov    %rbp,%r13
  lea    0x8(%rbp,%rax,8),%r14
  nopl   0x0(%rax)
L200:
  movapd %xmm1,%xmm0
  movsd  %xmm1,(%rsp)
  add    $0x8,%r13
  call   <rndgamma@plt>
  movsd  (%rsp),%xmm1
  divsd  %xmm1,%xmm0
  movsd  %xmm0,-0x8(%r13)
  cmp    %r14,%r13
  jne    L200
L37:
  test   %ebx,%ebx
  je     L101
  cmp    $0x1,%r12d
  jle    L31
  lea    -0x2(%r12),%edx
  movsd  0x0(%rbp),%xmm0
  lea    0x8(%rbp),%rax
  lea    0x10(%rbp,%rdx,8),%rdx
  nopl   0x0(%rax,%rax,1)
L196:
  addsd  (%rax),%xmm0
  add    $0x8,%rax
  movsd  %xmm0,-0x8(%rax)
  cmp    %rdx,%rax
  jne    L196
L31:
  movslq %r12d,%rax
  movsd  D14(%rip),%xmm0
  mov    %r12d,%edx
  mov    %rbp,%rsi
  divsd  -0x8(%rbp,%rax,8),%xmm0
  cvttsd2si %xmm0,%edi
  call   <abyx@plt>
  jmp    L101
  nop
L114:
  lea    0x0(%r13,%r13,1),%r15d
  movslq %r13d,%rax
  xor    %ecx,%ecx
  mov    %r13d,%edx
  movslq %r15d,%r15
  shl    $0x3,%rax
  movapd %xmm1,%xmm0
  mov    %r14,%rsi
  lea    0x0(,%r15,8),%r8
  lea    (%r14,%r8,1),%r15
  mov    %r8,0x8(%rsp)
  mov    %r15,%r10
  sub    %rax,%r10
  mov    %r10,%rdi
  mov    %r10,(%rsp)
  call   <DiscreteGamma@plt>
  mov    (%rsp),%r10
  mov    %r15,%rcx
  mov    %r13d,%esi
  mov    0x8(%rsp),%r8
  mov    %r12d,%edi
  mov    %r10,%rdx
  add    %r10,%r8
  call   <MultiNomial@plt>
  lea    -0x1(%r13),%eax
  mov    %r15,%rcx
  mov    %r14,%r8
  lea    0x8(%r14,%rax,8),%r9
  xor    %esi,%esi
  lea    0x8(%rbp),%r10
L59:
  mov    (%rcx),%edx
  lea    0x8(%r8),%rdi
  test   %edx,%edx
  jle    L190
  movslq %esi,%r11
  movslq %edx,%rdi
  lea    -0x1(%rdx),%r13d
  add    %r11,%rdi
  lea    0x0(%rbp,%r11,8),%rax
  lea    0x0(%rbp,%rdi,8),%rdi
  cmp    %rdi,%r8
  lea    0x8(%r8),%rdi
  setae  %r14b
  cmp    %rdi,%rax
  setae  %r15b
  or     %r15b,%r14b
  je     L237
  cmp    $0x2,%r13d
  jbe    L237
  mov    %edx,%r11d
  movsd  (%r8),%xmm0
  shr    %r11d
  sub    $0x1,%r11d
  unpcklpd %xmm0,%xmm0
  shl    $0x4,%r11
  lea    0x10(%rax,%r11,1),%r11
  xchg   %ax,%ax
L148:
  movups %xmm0,(%rax)
  add    $0x10,%rax
  cmp    %r11,%rax
  jne    L148
  mov    %edx,%r11d
  and    $0xfffffffe,%r11d
  lea    (%r11,%rsi,1),%eax
  cmp    %r11d,%edx
  je     L73
  movsd  (%r8),%xmm0
  cltq
  movsd  %xmm0,0x0(%rbp,%rax,8)
L73:
  add    %edx,%esi
L190:
  add    $0x4,%rcx
  mov    %rdi,%r8
  cmp    %r9,%rdi
  jne    L59
  jmp    L37
L237:
  add    %r11,%r13
  lea    (%r10,%r13,8),%r11
  nopw   0x0(%rax,%rax,1)
L160:
  movsd  (%r8),%xmm0
  add    $0x8,%rax
  movsd  %xmm0,-0x8(%rax)
  cmp    %r11,%rax
  jne    L160
  jmp    L73
L187:
  xor    %eax,%eax
  jmp    L168
```

tool_calls: name: unpack_data
arguments: data_tag: D11
type_infered: byte[8]
name: unpack_data
arguments: data_tag: D14
type_infered: byte[16]

tool: {"D11": [0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00]}

tool: {"D14": [0x00,0x00,0x00,0x00,0x00,0x00,0xf0,0x3f,0x00,0x00,0x00,0x00,0x00,0x00,0xf0,0x3f]}

assistant: ```c
int Rates4Sites(double rates[], double alpha, int ncatG, int ls, int cdf, double space[]) {
  int h, ir, j, *counts = (int *)(space + 2 * ncatG);
  double *rK = space, *freqK = space + ncatG;
  if (alpha == 0) {
    if (rates)
      for (h = 0; h < ls; h++)
        rates[h] = 1;
  } else {
    if (ncatG > 1) {
      DiscreteGamma(freqK, rK, alpha, alpha, ncatG, 0);
      MultiNomial(ls, ncatG, freqK, counts, space + 3 * ncatG);
      for (ir = 0, h = 0; ir < ncatG; ir++)
        for (j = 0; j < counts[ir]; j++)
          rates[h++] = rK[ir];
    } else
      for (h = 0; h < ls; h++)
        rates[h] = rndgamma(alpha) / alpha;
    if (cdf) {
      for (h = 1; h < ls; h++)
        rates[h] += rates[h - 1];
      abyx(1 / rates[ls - 1], rates, ls);
    }
  }
  return (0);
}
```