user: What is the c source code of the assembly code below?

```asm
<bitset>:
  endbr64
  push   %r15
  push   %r14
  push   %r13
  push   %r12
  push   %rbp
  push   %rbx
  sub    $0x838,%rsp
  mov    %rdi,(%rsp)
  mov    %rsi,%r15
  mov    %edx,%r14d
  mov    %rcx,0x18(%rsp)
  mov    %r8d,0x8(%rsp)
  mov    %fs:0x28,%rax
  mov    %rax,0x828(%rsp)
  xor    %eax,%eax
  mov    (%rsi),%ebp
  cmp    $0x1,%edx
  jle    L108
  lea    0x4(%rsi),%rdx
  lea    -0x2(%r14),%eax
  lea    0x8(%rsi,%rax,4),%rcx
  mov    %ebp,%r13d
L91:
  mov    (%rdx),%eax
  cmp    %eax,%r13d
  cmovl  %eax,%r13d
  cmp    %eax,%ebp
  cmovg  %eax,%ebp
  add    $0x4,%rdx
  cmp    %rcx,%rdx
  jne    L91
  cmp    %ebp,%r13d
  jge    L116
  mov    %ebp,%eax
  mov    %r13d,%ebp
  mov    %eax,%r13d
L116:
  test   %ebp,%ebp
  mov    $0x0,%eax
  cmovs  %eax,%ebp
  lea    -0x1(%r9),%eax
  cmp    %r9d,%r13d
  cmovge %eax,%r13d
  cmp    %ebp,%r13d
  jl     L34
  mov    0x8(%rsp),%eax
  imul   %ebp,%eax
  mov    %eax,0xc(%rsp)
  mov    D7(%rip),%r12d
  mov    %r13d,0x10(%rsp)
  jmp    L41
L108:
  mov    %ebp,%r13d
  jmp    L116
L6:
  add    $0x4,%rbx
  cmp    %r13,%rbx
  je     L105
L114:
  mov    0x4(%rbx),%eax
  add    (%rbx),%eax
  mov    %eax,%ecx
  shr    $0x1f,%ecx
  add    %eax,%ecx
  sar    %ecx
  mov    %ebp,%r8d
  mov    %r14d,%edx
  mov    %r15,%rsi
  mov    (%rsp),%rdi
  mov    $0x0,%eax
  call   <inside@plt>
  test   %eax,%eax
  je     L6
  mov    (%rbx),%ecx
  test   %ecx,%ecx
  mov    $0x0,%eax
  cmovs  %eax,%ecx
  mov    0x4(%rbx),%edx
  mov    0x8(%rsp),%edi
  lea    -0x1(%rdi),%eax
  cmp    %edi,%edx
  cmovge %eax,%edx
  cmp    %ecx,%edx
  jl     L6
  movslq %ecx,%rsi
  movslq 0x14(%rsp),%rax
  add    %rax,%rsi
  mov    0x18(%rsp),%rdi
  lea    (%rdi,%rsi,4),%rax
  sub    %ecx,%edx
  add    %rsi,%rdx
  lea    0x4(%rdi,%rdx,4),%rdx
L111:
  mov    %r12d,(%rax)
  add    $0x4,%rax
  cmp    %rdx,%rax
  jne    L111
  jmp    L6
L105:
  add    $0x1,%ebp
  mov    0x8(%rsp),%esi
  add    %esi,0xc(%rsp)
  cmp    %ebp,0x10(%rsp)
  jl     L34
L41:
  lea    0x20(%rsp),%r8
  mov    %ebp,%ecx
  mov    %r14d,%edx
  mov    %r15,%rsi
  mov    (%rsp),%rdi
  mov    $0x0,%eax
  call   <xxsect@plt>
  cmp    $0x1,%eax
  jle    L105
  mov    0xc(%rsp),%esi
  mov    %esi,0x14(%rsp)
  lea    0x20(%rsp),%rbx
  lea    -0x2(%rax),%eax
  lea    0x24(%rsp,%rax,4),%r13
  jmp    L114
L34:
  mov    0x828(%rsp),%rax
  sub    %fs:0x28,%rax
  jne    L49
  add    $0x838,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
L49:
  call   <__stack_chk_fail@plt>
```

tool_calls: name: disassemble_data
arguments: tag: D7
datatype: f32

tool: D7: 1

assistant: ```c
void bitset(xverts, yverts, nverts, mask, npix, nrast) float mask[];
int xverts[], yverts[], nverts;
int npix, nrast;
{
  int xsects[512];
  int ifrast, ilrast;
  int ifpix, ilpix;
  int nxsect = 0;
  int i, j, k, n, ix, itemp;
  ifrast = ilrast = yverts[0];
  for (i = 1; i < nverts; i++) {
    ilrast = ((yverts[i]) > (ilrast) ? (yverts[i]) : (ilrast));
    ifrast = ((yverts[i]) < (ifrast) ? (yverts[i]) : (ifrast));
  }
  if (ilrast < ifrast)
    (itemp = (ifrast), (ifrast) = (ilrast), (ilrast) = itemp);
  ifrast = ((0) > (ifrast) ? (0) : (ifrast));
  ilrast = (((nrast - 1)) < (ilrast) ? ((nrast - 1)) : (ilrast));
  for (j = ifrast; j <= ilrast; j++)
    if ((nxsect = xxsect(xverts, yverts, nverts, j, xsects)) > 1)
      for (i = 1; i < nxsect; i++) {
        n = (j * npix);
        ix = (xsects[i - 1] + xsects[i]) / 2;
        if (inside(xverts, yverts, nverts, ix, j)) {
          ifpix = ((0) > (xsects[i - 1]) ? (0) : (xsects[i - 1]));
          ilpix = (((npix - 1)) < (xsects[i]) ? ((npix - 1)) : (xsects[i]));
          for (k = ifpix; k <= ilpix; k++)
            mask[n + k] = 1.0;
        }
      }
}
```