user: Demonstrate the C code that generates the assembly output below:

```asm
<rand_gamma>:
  endbr64
  push   %rbx
  movapd %xmm0,%xmm7
  sub    $0x30,%rsp
  movsd  %xmm0,0x28(%rsp)
  movsd  D133(%rip),%xmm0
  comisd %xmm7,%xmm0
  ja     L63
  movsd  D19(%rip),%xmm7
  comisd 0x28(%rsp),%xmm7
  jae    L32
  movsd  D6(%rip),%xmm0
  comisd 0x28(%rsp),%xmm0
  ja     L122
  movsd  0x28(%rsp),%xmm7
  movapd %xmm7,%xmm6
  subsd  D19(%rip),%xmm6
  mulsd  D52(%rip),%xmm7
  movsd  %xmm6,0x10(%rsp)
  movapd %xmm7,%xmm0
  subsd  D120(%rip),%xmm0
  movq   %xmm0,%rbx
  xchg   %ax,%ax
L104:
  xor    %eax,%eax
  call   <rand_uniopen@plt>
  xor    %eax,%eax
  movsd  %xmm0,0x8(%rsp)
  call   <rand_uniopen@plt>
  movsd  0x8(%rsp),%xmm1
  pxor   %xmm5,%xmm5
  movsd  D19(%rip),%xmm2
  movapd %xmm0,%xmm3
  movq   %rbx,%xmm0
  subsd  %xmm1,%xmm2
  mulsd  %xmm1,%xmm2
  divsd  %xmm2,%xmm0
  ucomisd %xmm0,%xmm5
  ja     L78
  sqrtsd %xmm0,%xmm0
  subsd  D63(%rip),%xmm1
  movsd  0x10(%rsp),%xmm4
  mulsd  %xmm0,%xmm1
  addsd  %xmm1,%xmm4
  comisd %xmm5,%xmm4
  jb     L104
L107:
  movsd  D46(%rip),%xmm0
  movsd  %xmm1,0x8(%rsp)
  mulsd  %xmm2,%xmm0
  mulsd  %xmm2,%xmm0
  mulsd  %xmm2,%xmm0
  movapd %xmm1,%xmm2
  addsd  %xmm1,%xmm2
  mulsd  %xmm1,%xmm2
  mulsd  %xmm3,%xmm0
  divsd  %xmm4,%xmm2
  mulsd  %xmm3,%xmm0
  movsd  D19(%rip),%xmm3
  subsd  %xmm2,%xmm3
  comisd %xmm0,%xmm3
  jae    L85
  movsd  %xmm4,0x20(%rsp)
  call   <log@plt>
  movsd  0x20(%rsp),%xmm4
  movsd  %xmm0,0x18(%rsp)
  movapd %xmm4,%xmm7
  divsd  0x10(%rsp),%xmm7
  movapd %xmm7,%xmm0
  call   <log@plt>
  mulsd  0x10(%rsp),%xmm0
  movsd  0x8(%rsp),%xmm1
  movsd  0x20(%rsp),%xmm4
  subsd  %xmm1,%xmm0
  addsd  %xmm0,%xmm0
  comisd 0x18(%rsp),%xmm0
  jb     L104
L85:
  movsd  D9(%rip),%xmm0
  comisd %xmm4,%xmm0
  jbe    L39
  movsd  0x28(%rsp),%xmm7
  comisd %xmm4,%xmm7
  jbe    L39
  comisd %xmm7,%xmm0
  jbe    L138
  nopl   0x0(%rax)
L63:
  movsd  0x28(%rsp),%xmm4
L39:
  add    $0x30,%rsp
  movapd %xmm4,%xmm0
  pop    %rbx
  ret
L32:
  xor    %eax,%eax
  call   <rand_uniopen@plt>
  movsd  %xmm0,0x10(%rsp)
  movsd  D19(%rip),%xmm0
  addsd  0x28(%rsp),%xmm0
  call   <rand_gamma>
  movsd  D19(%rip),%xmm6
  divsd  0x28(%rsp),%xmm6
  cvttsd2si %xmm6,%edi
  movsd  0x10(%rsp),%xmm1
  movsd  %xmm0,0x8(%rsp)
  movapd %xmm1,%xmm0
  call   <pow@plt>
  movsd  0x8(%rsp),%xmm4
  mulsd  %xmm0,%xmm4
  jmp    L85
L138:
  movapd %xmm0,%xmm4
  add    $0x30,%rsp
  movapd %xmm4,%xmm0
  pop    %rbx
  ret
L122:
  xor    %eax,%eax
  call   <rand_exp@plt>
  movapd %xmm0,%xmm4
  jmp    L85
L78:
  movsd  %xmm2,0x20(%rsp)
  movsd  %xmm3,0x18(%rsp)
  movsd  %xmm1,0x8(%rsp)
  call   <sqrt@plt>
  movsd  0x8(%rsp),%xmm1
  pxor   %xmm7,%xmm7
  subsd  D63(%rip),%xmm1
  movsd  0x10(%rsp),%xmm4
  movsd  0x18(%rsp),%xmm3
  movsd  0x20(%rsp),%xmm2
  mulsd  %xmm0,%xmm1
  addsd  %xmm1,%xmm4
  comisd %xmm7,%xmm4
  jb     L104
  jmp    L107
```

tool_calls: name: disassemble_data
arguments: data_tag: D133
data_type: f64
name: disassemble_data
arguments: data_tag: D19
data_type: byte[8]
name: disassemble_data
arguments: data_tag: D6
data_type: f64
name: disassemble_data
arguments: data_tag: D52
data_type: byte[8]
name: disassemble_data
arguments: data_tag: D120
data_type: f64
name: disassemble_data
arguments: data_tag: D63
data_type: f64
name: disassemble_data
arguments: data_tag: D46
data_type: byte[8]
name: disassemble_data
arguments: data_tag: D9
data_type: byte[8]

tool: D133: 1e-05

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

tool: D6: 1.00001

tool: {"D52": [0x00,0x00,0x00,0x00,0x00,0x00,0x08,0x40]}

tool: D120: 0.75

tool: D63: 0.5

tool: {"D46": [0x00,0x00,0x00,0x00,0x00,0x00,0x50,0x40]}

tool: {"D9": [0xa0,0xc2,0xeb,0xfe,0x4b,0x48,0xb4,0x39]}

assistant: ```c
double rand_gamma(double a) {
  double b, c, X, Y, Z, U, V, W;
  if (a < 0.00001) {
    X = a;
  } else if (a <= 1) {
    U = rand_uniopen();
    X = rand_gamma(1 + a) * pow(U, 1 / a);
  } else if (a < 1.00001) {
    X = rand_exp();
  } else {
    b = a - 1;
    c = 3 * a - 0.75;
    for (;;) {
      U = rand_uniopen();
      V = rand_uniopen();
      W = U * (1 - U);
      Y = sqrt(c / W) * (U - 0.5);
      X = b + Y;
      if (X >= 0) {
        Z = 64 * W * W * W * V * V;
        if (Z <= 1 - 2 * Y * Y / X || log(Z) <= 2 * (b * log(X / b) - Y))
          break;
      }
    }
  }
  return X < 1e-30 && X < a ? (a < 1e-30 ? a : 1e-30) : X;
}
```