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

```asm
<projillum>:
  endbr64
  push   %r15
  push   %r14
  push   %r13
  push   %r12
  push   %rbp
  push   %rbx
  sub    $0x58,%rsp
  mov    %rdi,%r14
  mov    %esi,0x1c(%rsp)
  mov    %edx,%r15d
  movsd  %xmm0,0x38(%rsp)
  test   %edx,%edx
  jle    L184
  mov    %rdi,%rax
  lea    -0x1(%rdx),%edx
  lea    0x2(%rdi,%rdx,2),%rdx
L109:
  movw   $0xffff,(%rax)
  add    $0x2,%rax
  cmp    %rdx,%rax
  jne    L109
L184:
  movsd  0x38(%rsp),%xmm3
  mulsd  D59(%rip),%xmm3
  movq   %xmm3,%rbx
  movapd %xmm3,%xmm0
  xorpd  D137(%rip),%xmm0
  call   <sin@plt>
  movq   %xmm0,%r12
  movq   %rbx,%xmm0
  call   <cos@plt>
  movsd  %xmm0,0x20(%rsp)
  mov    $0x65,%ebp
  movsd  D95(%rip),%xmm3
  movsd  %xmm3,0x8(%rsp)
  mov    $0x0,%ebx
  movl   $0x0,0x18(%rsp)
  mov    $0x1,%r13d
  movq   %r12,%xmm3
  xorpd  D137(%rip),%xmm3
  movsd  %xmm3,0x28(%rsp)
  movsd  D103(%rip),%xmm3
  movsd  %xmm3,0x30(%rsp)
  mov    %r14,%r12
  jmp    L125
L86:
  movapd %xmm2,%xmm0
  call   <atan2@plt>
  divsd  D59(%rip),%xmm0
  movq   %xmm0,%r14
  jmp    L0
L51:
  mov    %edi,%esi
  sub    %ebx,%esi
  mov    %r8d,%eax
  mov    0x18(%rsp),%r11d
  sub    %r11d,%eax
  pxor   %xmm3,%xmm3
  cvtsi2sd %eax,%xmm3
  pxor   %xmm0,%xmm0
  cvtsi2sd %esi,%xmm0
  divsd  %xmm0,%xmm3
  sar    $0x1f,%esi
  or     $0x1,%esi
  mov    %esi,%r9d
  neg    %esi
  movslq %esi,%rsi
  add    %rsi,%rsi
  movslq %r15d,%rax
  sub    $0x1,%rax
  movslq %ebx,%rdx
  sub    %rdx,%rax
  lea    (%r12,%rax,2),%rcx
  mov    $0x0,%edx
  movsd  D82(%rip),%xmm6
  movsd  D24(%rip),%xmm5
  movsd  D133(%rip),%xmm7
  pxor   %xmm4,%xmm4
  cvtsi2sd %r11d,%xmm4
  mov    $0x1,%r10d
  movsd  0x30(%rsp),%xmm8
  jmp    L146
L3:
  addsd  %xmm4,%xmm1
  cvttsd2si %xmm1,%eax
  test   %eax,%eax
  cmove  %r10d,%eax
  mov    %ax,(%rcx)
  add    %r9d,%edx
  add    %rsi,%rcx
  lea    (%rbx,%rdx,1),%eax
  cmp    %eax,%edi
  je     L81
L146:
  pxor   %xmm0,%xmm0
  cvtsi2sd %edx,%xmm0
  mulsd  %xmm3,%xmm0
  addsd  %xmm6,%xmm0
  movapd %xmm0,%xmm1
  movapd %xmm0,%xmm2
  andpd  %xmm8,%xmm2
  ucomisd %xmm2,%xmm5
  jbe    L3
  cvttsd2si %xmm0,%rax
  pxor   %xmm2,%xmm2
  cvtsi2sd %rax,%xmm2
  movapd %xmm2,%xmm1
  cmpnlesd %xmm0,%xmm1
  andpd  %xmm7,%xmm1
  subsd  %xmm1,%xmm2
  movapd %xmm8,%xmm1
  andnpd %xmm0,%xmm1
  orpd   %xmm2,%xmm1
  jmp    L3
L81:
  movsd  0x8(%rsp),%xmm7
  addsd  D43(%rip),%xmm7
  movsd  %xmm7,0x8(%rsp)
  mov    $0x0,%r13d
  sub    $0x1,%ebp
  je     L165
  mov    %edi,%ebx
  mov    %r8d,0x18(%rsp)
L125:
  lea    0x48(%rsp),%rdi
  lea    0x40(%rsp),%rsi
  movsd  0x8(%rsp),%xmm0
  call   <sincos@plt>
  movsd  0x40(%rsp),%xmm2
  movsd  0x48(%rsp),%xmm0
  movsd  0x28(%rsp),%xmm1
  mulsd  %xmm0,%xmm1
  mulsd  0x20(%rsp),%xmm0
  movsd  %xmm0,0x10(%rsp)
  pxor   %xmm5,%xmm5
  ucomisd %xmm5,%xmm2
  setp   %al
  mov    $0x1,%edx
  cmovne %edx,%eax
  test   %al,%al
  jne    L86
  ucomisd %xmm5,%xmm1
  setp   %al
  cmovne %edx,%eax
  movq   %xmm5,%r14
  test   %al,%al
  jne    L86
L0:
  movsd  0x10(%rsp),%xmm0
  call   <asin@plt>
  pxor   %xmm1,%xmm1
  cvtsi2sd %r15d,%xmm1
  divsd  D59(%rip),%xmm0
  addsd  D23(%rip),%xmm0
  movapd %xmm1,%xmm2
  divsd  D86(%rip),%xmm2
  mulsd  %xmm2,%xmm0
  subsd  %xmm0,%xmm1
  cvttsd2si %xmm1,%edi
  pxor   %xmm0,%xmm0
  cvtsi2sdl 0x1c(%rsp),%xmm0
  divsd  D77(%rip),%xmm0
  movq   %r14,%xmm6
  mulsd  %xmm6,%xmm0
  cvttsd2si %xmm0,%r8d
  test   %r13d,%r13d
  jne    L81
  cmp    %ebx,%edi
  jne    L51
  test   %r8d,%r8d
  mov    $0x1,%edx
  cmovne %r8d,%edx
  lea    -0x1(%r15),%eax
  sub    %ebx,%eax
  cltq
  mov    %dx,(%r12,%rax,2)
  jmp    L81
L165:
  mov    %r12,%r14
  pxor   %xmm0,%xmm0
  lea    -0x1(%r15),%r9d
  mov    $0xffffffff,%ecx
  comisd 0x38(%rsp),%xmm0
  ja     L61
  mov    $0x1,%ecx
  mov    $0x0,%r9d
L61:
  mov    %r15d,%esi
  shr    $0x1f,%esi
  add    %r15d,%esi
  sar    %esi
  cmp    %esi,%r9d
  je     L29
  movslq %ecx,%r10
  lea    (%r10,%r10,1),%r8
  movslq %r9d,%rax
  lea    (%r14,%rax,2),%rdx
  mov    %r9d,%eax
L154:
  mov    %rdx,%rdi
  cmpw   $0xffff,(%rdx)
  jne    L69
  add    %ecx,%eax
  add    %r8,%rdx
  cmp    %esi,%eax
  jne    L154
L29:
  add    $0x58,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
L69:
  mov    0x1c(%rsp),%ebx
  mov    %ebx,%edx
  shr    $0x1f,%edx
  add    %ebx,%edx
  sar    %edx
  mov    %edx,%esi
  mov    %dx,(%rdi)
  cmp    %eax,%r9d
  je     L29
  neg    %ecx
  movslq %ecx,%rdx
  cltq
  neg    %r10
  lea    (%r14,%r10,2),%rcx
L127:
  mov    %si,(%rcx,%rax,2)
  add    %rdx,%rax
  cmp    %eax,%r9d
  jne    L127
  jmp    L29
```

tool_calls: name: parse_data
arguments: label: D59
data_type: byte[8]
name: parse_data
arguments: label: D137
data_type: f32[4]
name: parse_data
arguments: label: D95
data_type: byte[8]
name: parse_data
arguments: label: D103
data_type: byte[8]
name: parse_data
arguments: label: D82
data_type: f64
name: parse_data
arguments: label: D24
data_type: f32[2]
name: parse_data
arguments: label: D133
data_type: f32[2]
name: parse_data
arguments: label: D43
data_type: byte[8]
name: parse_data
arguments: label: D23
data_type: f32[2]
name: parse_data
arguments: label: D86
data_type: f64
name: parse_data
arguments: label: D77
data_type: f64

tool: {"D59": [0x39,0x9d,0x52,0xa2,0x46,0xdf,0x91,0x3f]}

tool: {"D137": [0.0,-0.0,0.0,0.0]}

tool: {"D95": [0x18,0x2d,0x44,0x54,0xfb,0x21,0xf9,0xbf]}

tool: {"D103": [0xff,0xff,0xff,0xff,0xff,0xff,0xff,0x7f]}

tool: D82: 0.5

tool: {"D24": [0.0,176.0]}

tool: {"D133": [0.0,1.875]}

tool: {"D43": [0x1a,0x27,0x17,0x92,0xbf,0x15,0xa0,0x3f]}

tool: {"D23": [0.0,3.3515625]}

tool: D86: 180

tool: D77: 360

assistant: ```c
void projillum(wtab, xdots, ydots, dec) short *wtab;
int xdots, ydots;
double dec;
{
  int i, ftf = 1, ilon, ilat, lilon = 0, lilat = 0, xt;
  double m, x, y, z, th, lon, lat, s, c;
  for (i = 0; i < ydots; i++)
    wtab[i] = -1;
  s = sin(-((dec) * (3.14159265358979323846 / 180.0)));
  c = cos(-((dec) * (3.14159265358979323846 / 180.0)));
  for (th = -(3.14159265358979323846 / 2); th <= 3.14159265358979323846 / 2 + 0.001; th += 3.14159265358979323846 / 100) {
    x = -s * sin(th);
    y = cos(th);
    z = c * sin(th);
    lon = (y == 0 && x == 0) ? 0.0 : ((atan2(y, x)) / (3.14159265358979323846 / 180.0));
    lat = ((asin(z)) / (3.14159265358979323846 / 180.0));
    ilat = ydots - (lat + 90) * (ydots / 180.0);
    ilon = lon * (xdots / 360.0);
    if (ftf) {
      lilon = ilon;
      lilat = ilat;
      ftf = 0;
    } else {
      if (lilat == ilat) {
        wtab[(ydots - 1) - ilat] = ilon == 0 ? 1 : ilon;
      } else {
        m = ((double)(ilon - lilon)) / (ilat - lilat);
        for (i = lilat; i != ilat; i += (((ilat - lilat) < 0) ? -1 : ((ilat - lilat) > 0 ? 1 : 0))) {
          xt = lilon + floor((m * (i - lilat)) + 0.5);
          wtab[(ydots - 1) - i] = xt == 0 ? 1 : xt;
        }
      }
      lilon = ilon;
      lilat = ilat;
    }
  }
  if (dec < 0.0) {
    ilat = ydots - 1;
    lilat = -1;
  } else {
    ilat = 0;
    lilat = 1;
  }
  for (i = ilat; i != ydots / 2; i += lilat) {
    if (wtab[i] != -1) {
      while (1) {
        wtab[i] = xdots / 2;
        if (i == ilat)
          break;
        i -= lilat;
      }
      break;
    }
  }
}
```