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

```asm
<ATM_Spectrum>:
  endbr64
  push   %r13
  mov    %rcx,%r13
  push   %r12
  mov    %rdx,%r12
  push   %rbp
  mov    %rsi,%rbp
  mov    %r12,%rsi
  push   %rbx
  mov    %rdi,%rbx
  sub    $0x28,%rsp
  mov    %fs:0x28,%rax
  mov    %rax,0x18(%rsp)
  xor    %eax,%eax
  lea    0x10(%rsp),%rdx
  call   <ATM_Dos@plt>
  movsd  (%rbx),%xmm1
  movsd  (%r12),%xmm0
  movsd  0x8(%rbx),%xmm4
  movsd  0x10(%rbx),%xmm3
  comisd %xmm0,%xmm1
  movsd  0x18(%rbx),%xmm5
  movsd  0x0(%rbp),%xmm2
  movsd  0x8(%rbp),%xmm10
  movsd  0x10(%rbp),%xmm9
  movsd  0x18(%rbp),%xmm7
  ja     L50
  comisd %xmm1,%xmm0
  movapd %xmm4,%xmm11
  movapd %xmm3,%xmm13
  movapd %xmm5,%xmm12
  movapd %xmm10,%xmm8
  movapd %xmm9,%xmm15
  movapd %xmm7,%xmm6
  subsd  %xmm1,%xmm11
  subsd  %xmm1,%xmm13
  subsd  %xmm1,%xmm12
  subsd  %xmm2,%xmm8
  subsd  %xmm2,%xmm15
  subsd  %xmm2,%xmm6
  jbe    L46
  comisd %xmm0,%xmm4
  jbe    L46
  divsd  %xmm11,%xmm8
  movapd %xmm15,%xmm3
  subsd  %xmm1,%xmm0
  mulsd  D29(%rip),%xmm0
  divsd  %xmm13,%xmm3
  divsd  %xmm12,%xmm6
  addsd  %xmm3,%xmm8
  addsd  %xmm6,%xmm8
  mulsd  %xmm8,%xmm0
  addsd  %xmm2,%xmm0
  mulsd  0x10(%rsp),%xmm0
  movsd  %xmm0,0x0(%r13)
  jmp    L35
  nopw   0x0(%rax,%rax,1)
L46:
  movapd %xmm5,%xmm1
  comisd %xmm4,%xmm0
  movapd %xmm7,%xmm14
  subsd  %xmm4,%xmm1
  subsd  %xmm10,%xmm14
  movapd %xmm7,%xmm10
  subsd  %xmm9,%xmm10
  movsd  %xmm1,0x8(%rsp)
  movapd %xmm5,%xmm1
  subsd  %xmm3,%xmm1
  ja     L96
L2:
  comisd %xmm3,%xmm0
  ja     L102
L13:
  comisd %xmm5,%xmm0
  jbe    L35
L50:
  movq   $0x0,0x0(%r13)
L35:
  mov    0x18(%rsp),%rax
  sub    %fs:0x28,%rax
  jne    L82
  add    $0x28,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  ret
  nopw   0x0(%rax,%rax,1)
L96:
  comisd %xmm0,%xmm3
  jbe    L2
  mulsd  %xmm11,%xmm15
  movapd %xmm3,%xmm5
  subsd  %xmm0,%xmm3
  mulsd  %xmm6,%xmm11
  subsd  %xmm4,%xmm5
  subsd  %xmm4,%xmm0
  movapd %xmm15,%xmm7
  divsd  %xmm13,%xmm7
  divsd  %xmm12,%xmm11
  addsd  %xmm8,%xmm7
  movsd  D29(%rip),%xmm8
  addsd  %xmm7,%xmm11
  mulsd  %xmm8,%xmm11
  mulsd  %xmm11,%xmm3
  divsd  %xmm5,%xmm3
  addsd  %xmm2,%xmm3
  movapd %xmm1,%xmm2
  mulsd  %xmm6,%xmm2
  mulsd  %xmm14,%xmm1
  divsd  0x8(%rsp),%xmm1
  divsd  %xmm12,%xmm2
  addsd  %xmm2,%xmm1
  addsd  %xmm10,%xmm1
  mulsd  %xmm8,%xmm1
  subsd  %xmm1,%xmm6
  mulsd  %xmm0,%xmm6
  divsd  %xmm5,%xmm6
  addsd  %xmm6,%xmm3
  mulsd  0x10(%rsp),%xmm3
  movsd  %xmm3,0x0(%r13)
  jmp    L35
  nopl   0x0(%rax,%rax,1)
L102:
  comisd %xmm0,%xmm5
  jbe    L13
  divsd  %xmm12,%xmm6
  divsd  0x8(%rsp),%xmm14
  subsd  %xmm5,%xmm0
  mulsd  D29(%rip),%xmm0
  divsd  %xmm1,%xmm10
  addsd  %xmm14,%xmm6
  addsd  %xmm10,%xmm6
  mulsd  %xmm6,%xmm0
  addsd  %xmm7,%xmm0
  mulsd  0x10(%rsp),%xmm0
  movsd  %xmm0,0x0(%r13)
  jmp    L35
L82:
  call   <__stack_chk_fail@plt>
```

tool_calls: name: read_data
arguments: data_label: D29
type_infered: byte[8]

tool: {"D29": [0x55,0x55,0x55,0x55,0x55,0x55,0xd5,0x3f]}

assistant: ```c
void ATM_Spectrum(double *et, double *at, double *e, double *spectrum) {
  double a21, a31, a32, a41;
  double e21, e31, e32, e41, e42, e43;
  double a42, a43;
  double dos;
  --at;
  ATM_Dos(et, e, &dos);
  --et;
  e21 = et[2] - et[1];
  e31 = et[3] - et[1];
  e32 = et[3] - et[2];
  e41 = et[4] - et[1];
  e42 = et[4] - et[2];
  e43 = et[4] - et[3];
  a21 = at[2] - at[1];
  a31 = at[3] - at[1];
  a32 = at[3] - at[2];
  a41 = at[4] - at[1];
  a42 = at[4] - at[2];
  a43 = at[4] - at[3];
  if (*e < et[1]) {
    *spectrum = 0.;
  } else if (*e > et[1] && *e < et[2]) {
    *spectrum = dos * (at[1] + (*e - et[1]) * .33333333333333331 * (a21 / e21 + a31 / e31 + a41 / e41));
  } else if (*e > et[2] && *e < et[3]) {
    *spectrum = dos * (at[1] + (a21 + e21 * a31 / e31 + e21 * a41 / e41) * .33333333333333331 * (et[3] - *e) / e32 + (at[4] - at[1] - (e43 * a41 / e41 + e43 * a42 / e42 + a43) * .33333333333333331) * (*e - et[2]) / e32);
  } else if (*e > et[3] && *e < et[4]) {
    *spectrum = dos * (at[4] + (*e - et[4]) * .33333333333333331 * (a41 / e41 + a42 / e42 + a43 / e43));
  } else if (*e > et[4]) {
    *spectrum = 0.;
  }
}
```