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

```asm
<form_flat_clusters_maxclust_monocrit>:
  endbr64
  push   %rbp
  mov    %rsp,%rbp
  sub    $0x90,%rsp
  mov    %rdi,-0x78(%rbp)
  mov    %rsi,-0x80(%rbp)
  mov    %rdx,-0x88(%rbp)
  mov    %ecx,-0x8c(%rbp)
  mov    %r8d,-0x90(%rbp)
  movsd  D4(%rip),%xmm0
  movsd  %xmm0,-0x40(%rbp)
  pxor   %xmm0,%xmm0
  movsd  %xmm0,-0x38(%rbp)
  movl   $0x1,-0x58(%rbp)
  pxor   %xmm0,%xmm0
  cvtsi2sdl -0x8c(%rbp),%xmm0
  movsd  D16(%rip),%xmm2
  movapd %xmm0,%xmm1
  divsd  %xmm2,%xmm1
  mov    -0x8c(%rbp),%eax
  cltq
  shr    $0x3,%rax
  test   %rax,%rax
  js     L238
  pxor   %xmm0,%xmm0
  cvtsi2sd %rax,%xmm0
  jmp    L69
L238:
  mov    %rax,%rdx
  shr    %rdx
  and    $0x1,%eax
  or     %rax,%rdx
  pxor   %xmm0,%xmm0
  cvtsi2sd %rdx,%xmm0
  addsd  %xmm0,%xmm0
L69:
  ucomisd %xmm0,%xmm1
  jp     L222
  ucomisd %xmm0,%xmm1
  jne    L222
  mov    -0x8c(%rbp),%eax
  cltq
  shr    $0x3,%rax
  jmp    L132
L222:
  mov    -0x8c(%rbp),%eax
  cltq
  shr    $0x3,%rax
  add    $0x1,%eax
L132:
  mov    %eax,-0x54(%rbp)
  movl   $0x0,-0x64(%rbp)
  mov    -0x8c(%rbp),%eax
  cltq
  shl    $0x3,%rax
  lea    -0x10(%rax),%rdx
  mov    -0x80(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  movsd  %xmm0,-0x38(%rbp)
  mov    -0x8c(%rbp),%eax
  cltq
  shl    $0x2,%rax
  mov    %rax,%rdi
  call   <malloc@plt>
  mov    %rax,-0x30(%rbp)
  mov    -0x54(%rbp),%eax
  cltq
  mov    %rax,%rdi
  call   <malloc@plt>
  mov    %rax,-0x28(%rbp)
  mov    -0x54(%rbp),%eax
  cltq
  mov    %rax,%rdi
  call   <malloc@plt>
  mov    %rax,-0x20(%rbp)
  mov    -0x8c(%rbp),%eax
  lea    -0x1(%rax),%edx
  mov    -0x64(%rbp),%eax
  cltq
  lea    0x0(,%rax,4),%rcx
  mov    -0x30(%rbp),%rax
  add    %rcx,%rax
  add    %edx,%edx
  mov    %edx,(%rax)
  mov    -0x54(%rbp),%eax
  cltq
  mov    -0x28(%rbp),%rdx
  mov    %rdx,%rcx
  mov    %rax,%rdx
  mov    $0x0,%esi
  mov    %rcx,%rdi
  call   <memset@plt>
  mov    -0x54(%rbp),%eax
  cltq
  mov    -0x20(%rbp),%rdx
  mov    %rdx,%rcx
  mov    %rax,%rdx
  mov    $0x0,%esi
  mov    %rcx,%rdi
  call   <memset@plt>
  movl   $0x0,-0x60(%rbp)
  mov    -0x8c(%rbp),%eax
  sub    $0x2,%eax
  mov    %eax,-0x5c(%rbp)
  jmp    L126
L75:
  mov    -0x5c(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x80(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  movsd  %xmm0,-0x18(%rbp)
  movsd  -0x18(%rbp),%xmm0
  comisd -0x38(%rbp),%xmm0
  ja     L195
  movl   $0x0,-0x64(%rbp)
  mov    -0x8c(%rbp),%eax
  lea    -0x1(%rax),%edx
  mov    -0x64(%rbp),%eax
  cltq
  lea    0x0(,%rax,4),%rcx
  mov    -0x30(%rbp),%rax
  add    %rcx,%rax
  add    %edx,%edx
  mov    %edx,(%rax)
  mov    -0x54(%rbp),%eax
  cltq
  mov    -0x28(%rbp),%rdx
  mov    %rdx,%rcx
  mov    %rax,%rdx
  mov    $0x0,%esi
  mov    %rcx,%rdi
  call   <memset@plt>
  mov    -0x54(%rbp),%eax
  cltq
  mov    -0x20(%rbp),%rdx
  mov    %rdx,%rcx
  mov    %rax,%rdx
  mov    $0x0,%esi
  mov    %rcx,%rdi
  call   <memset@plt>
  movl   $0x0,-0x60(%rbp)
  movl   $0xffffffff,-0x50(%rbp)
  jmp    L97
L251:
  mov    -0x64(%rbp),%eax
  cltq
  lea    0x0(,%rax,4),%rdx
  mov    -0x30(%rbp),%rax
  add    %rdx,%rax
  mov    (%rax),%eax
  mov    %eax,-0x4c(%rbp)
  mov    -0x4c(%rbp),%eax
  sub    -0x8c(%rbp),%eax
  shl    $0x2,%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x78(%rbp),%rax
  add    %rdx,%rax
  mov    %rax,-0x10(%rbp)
  mov    -0x10(%rbp),%rax
  movsd  (%rax),%xmm0
  cvttsd2si %xmm0,%eax
  mov    %eax,-0x48(%rbp)
  mov    -0x10(%rbp),%rax
  add    $0x8,%rax
  movsd  (%rax),%xmm0
  cvttsd2si %xmm0,%eax
  mov    %eax,-0x44(%rbp)
  mov    -0x4c(%rbp),%eax
  sub    -0x8c(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x80(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  movsd  %xmm0,-0x8(%rbp)
  movsd  -0x18(%rbp),%xmm0
  comisd -0x8(%rbp),%xmm0
  jb     L44
  addl   $0x1,-0x60(%rbp)
  subl   $0x1,-0x64(%rbp)
  mov    -0x4c(%rbp),%eax
  sub    -0x8c(%rbp),%eax
  cltq
  shr    $0x3,%rax
  mov    %rax,%rdx
  mov    -0x28(%rbp),%rax
  add    %rdx,%rax
  movzbl (%rax),%eax
  mov    %eax,%edi
  mov    -0x4c(%rbp),%eax
  sub    -0x8c(%rbp),%eax
  not    %eax
  and    $0x7,%eax
  mov    $0x1,%esi
  mov    %eax,%ecx
  shl    %cl,%esi
  mov    %esi,%eax
  or     %eax,%edi
  mov    %edi,%ecx
  mov    -0x28(%rbp),%rax
  add    %rdx,%rax
  mov    %ecx,%edx
  mov    %dl,(%rax)
  mov    -0x4c(%rbp),%eax
  sub    -0x8c(%rbp),%eax
  cltq
  shr    $0x3,%rax
  mov    %rax,%rdx
  mov    -0x20(%rbp),%rax
  add    %rdx,%rax
  movzbl (%rax),%eax
  mov    %eax,%edi
  mov    -0x4c(%rbp),%eax
  sub    -0x8c(%rbp),%eax
  not    %eax
  and    $0x7,%eax
  mov    $0x1,%esi
  mov    %eax,%ecx
  shl    %cl,%esi
  mov    %esi,%eax
  or     %eax,%edi
  mov    %edi,%ecx
  mov    -0x20(%rbp),%rax
  add    %rdx,%rax
  mov    %ecx,%edx
  mov    %dl,(%rax)
  jmp    L97
L44:
  mov    -0x4c(%rbp),%eax
  sub    -0x8c(%rbp),%eax
  cltq
  shr    $0x3,%rax
  mov    %rax,%rdx
  mov    -0x28(%rbp),%rax
  add    %rdx,%rax
  movzbl (%rax),%eax
  movzbl %al,%edx
  mov    -0x4c(%rbp),%eax
  sub    -0x8c(%rbp),%eax
  not    %eax
  and    $0x7,%eax
  mov    %eax,%ecx
  sar    %cl,%edx
  mov    %edx,%eax
  and    $0x1,%eax
  test   %eax,%eax
  jne    L142
  mov    -0x4c(%rbp),%eax
  sub    -0x8c(%rbp),%eax
  cltq
  shr    $0x3,%rax
  mov    %rax,%rdx
  mov    -0x28(%rbp),%rax
  add    %rdx,%rax
  movzbl (%rax),%eax
  mov    %eax,%edi
  mov    -0x4c(%rbp),%eax
  sub    -0x8c(%rbp),%eax
  not    %eax
  and    $0x7,%eax
  mov    $0x1,%esi
  mov    %eax,%ecx
  shl    %cl,%esi
  mov    %esi,%eax
  or     %eax,%edi
  mov    %edi,%ecx
  mov    -0x28(%rbp),%rax
  add    %rdx,%rax
  mov    %ecx,%edx
  mov    %dl,(%rax)
  mov    -0x48(%rbp),%eax
  cmp    -0x8c(%rbp),%eax
  jl     L22
  mov    -0x64(%rbp),%eax
  cltq
  add    $0x1,%rax
  lea    0x0(,%rax,4),%rdx
  mov    -0x30(%rbp),%rax
  add    %rax,%rdx
  mov    -0x48(%rbp),%eax
  mov    %eax,(%rdx)
  addl   $0x1,-0x64(%rbp)
  jmp    L97
L22:
  mov    -0x48(%rbp),%eax
  cmp    -0x8c(%rbp),%eax
  jge    L142
  addl   $0x1,-0x60(%rbp)
L142:
  mov    -0x4c(%rbp),%eax
  sub    -0x8c(%rbp),%eax
  cltq
  shr    $0x3,%rax
  mov    %rax,%rdx
  mov    -0x20(%rbp),%rax
  add    %rdx,%rax
  movzbl (%rax),%eax
  movzbl %al,%edx
  mov    -0x4c(%rbp),%eax
  sub    -0x8c(%rbp),%eax
  not    %eax
  and    $0x7,%eax
  mov    %eax,%ecx
  sar    %cl,%edx
  mov    %edx,%eax
  and    $0x1,%eax
  test   %eax,%eax
  jne    L96
  mov    -0x44(%rbp),%eax
  cmp    -0x8c(%rbp),%eax
  jl     L211
  mov    -0x4c(%rbp),%eax
  sub    -0x8c(%rbp),%eax
  cltq
  shr    $0x3,%rax
  mov    %rax,%rdx
  mov    -0x20(%rbp),%rax
  add    %rdx,%rax
  movzbl (%rax),%eax
  mov    %eax,%edi
  mov    -0x4c(%rbp),%eax
  sub    -0x8c(%rbp),%eax
  not    %eax
  and    $0x7,%eax
  mov    $0x1,%esi
  mov    %eax,%ecx
  shl    %cl,%esi
  mov    %esi,%eax
  or     %eax,%edi
  mov    %edi,%ecx
  mov    -0x20(%rbp),%rax
  add    %rdx,%rax
  mov    %ecx,%edx
  mov    %dl,(%rax)
  mov    -0x64(%rbp),%eax
  cltq
  add    $0x1,%rax
  lea    0x0(,%rax,4),%rdx
  mov    -0x30(%rbp),%rax
  add    %rax,%rdx
  mov    -0x44(%rbp),%eax
  mov    %eax,(%rdx)
  addl   $0x1,-0x64(%rbp)
  jmp    L97
L211:
  mov    -0x44(%rbp),%eax
  cmp    -0x8c(%rbp),%eax
  jge    L96
  addl   $0x1,-0x60(%rbp)
L96:
  subl   $0x1,-0x64(%rbp)
L97:
  cmpl   $0x0,-0x64(%rbp)
  jns    L251
  movsd  -0x18(%rbp),%xmm0
  comisd -0x40(%rbp),%xmm0
  jbe    L209
  mov    -0x60(%rbp),%eax
  cmp    -0x90(%rbp),%eax
  jle    L209
  movsd  -0x18(%rbp),%xmm0
  movsd  %xmm0,-0x40(%rbp)
  jmp    L231
L209:
  movsd  -0x38(%rbp),%xmm0
  comisd -0x18(%rbp),%xmm0
  jbe    L231
  mov    -0x60(%rbp),%eax
  cmp    -0x90(%rbp),%eax
  jg     L231
  movsd  -0x18(%rbp),%xmm0
  movsd  %xmm0,-0x38(%rbp)
  mov    -0x60(%rbp),%eax
  mov    %eax,-0x58(%rbp)
  jmp    L231
L195:
  nop
L231:
  subl   $0x1,-0x5c(%rbp)
L126:
  cmpl   $0x0,-0x5c(%rbp)
  jns    L75
  mov    -0x8c(%rbp),%ecx
  mov    -0x38(%rbp),%rdi
  mov    -0x88(%rbp),%rdx
  mov    -0x80(%rbp),%rsi
  mov    -0x78(%rbp),%rax
  movq   %rdi,%xmm0
  mov    %rax,%rdi
  mov    $0x1,%eax
  call   <form_flat_clusters_from_monotonic_criterion@plt>
  mov    -0x30(%rbp),%rax
  mov    %rax,%rdi
  call   <free@plt>
  mov    -0x28(%rbp),%rax
  mov    %rax,%rdi
  call   <free@plt>
  mov    -0x20(%rbp),%rax
  mov    %rax,%rdi
  call   <free@plt>
  nop
  leave
  ret
```

tool_calls: name: parse_data
arguments: data_label: D4
type_guessed: f32[2]
name: parse_data
arguments: data_label: D16
type_guessed: f32[2]

tool: {"D4": [0.0,-1.875]}

tool: {"D16": [0.0,2.5]}

assistant: ```c
void form_flat_clusters_maxclust_monocrit(const double *Z, const double *mono_crit, int *T, int n, int mc) {
  int *curNode;
  int ndid, lid, rid, k, nc, g, ms;
  unsigned char *lvisited, *rvisited;
  const double *Zrow;
  double thresh, maxmono_crit;
  double max_illegal = -1.0;
  double min_legal = 0.0;
  int min_legal_nc = 1;
  const int bff = (((((double)(n)) / (double)(sizeof(unsigned char) * 8)) == ((double)(((n)) / ((sizeof(unsigned char) * 8)))) ? (((n)) / ((sizeof(unsigned char) * 8))) : (((n)) / ((sizeof(unsigned char) * 8)) + 1)));
  k = 0;
  min_legal = mono_crit[n - 2];
  curNode = (int *)malloc(n * sizeof(int));
  lvisited = (unsigned char *)malloc(bff);
  rvisited = (unsigned char *)malloc(bff);
  curNode[k] = (n * 2) - 2;
  bzero(lvisited, bff);
  bzero(rvisited, bff);
  nc = 0;
  ;
  for (g = n - 2; g >= 0; g--) {
    thresh = mono_crit[g];
    if (thresh > min_legal) {
      continue;
    }
    k = 0;
    curNode[k] = (n * 2) - 2;
    bzero(lvisited, bff);
    bzero(rvisited, bff);
    nc = 0;
    ms = -1;
    while (k >= 0) {
      ndid = curNode[k];
      Zrow = Z + ((ndid - n) * 4);
      lid = (int)Zrow[0];
      rid = (int)Zrow[1];
      maxmono_crit = mono_crit[ndid - n];
      if (maxmono_crit <= thresh) {
        nc++;
        k--;
        ((lvisited)[(ndid - n) / (sizeof(unsigned char) * 8)] |= ((0x1) << (((sizeof(unsigned char) * 8) - 1) - ((ndid - n) % (sizeof(unsigned char) * 8)))));
        ((rvisited)[(ndid - n) / (sizeof(unsigned char) * 8)] |= ((0x1) << (((sizeof(unsigned char) * 8) - 1) - ((ndid - n) % (sizeof(unsigned char) * 8)))));
        continue;
      }
      if (!(((lvisited)[(ndid - n) / (sizeof(unsigned char) * 8)] >> (((sizeof(unsigned char) * 8) - 1) - ((ndid - n) % (sizeof(unsigned char) * 8)))) & 0x1)) {
        ((lvisited)[(ndid - n) / (sizeof(unsigned char) * 8)] |= ((0x1) << (((sizeof(unsigned char) * 8) - 1) - ((ndid - n) % (sizeof(unsigned char) * 8)))));
        if (lid >= n) {
          curNode[k + 1] = lid;
          k++;
          continue;
        } else if (lid < n) {
          nc++;
        }
      }
      if (!(((rvisited)[(ndid - n) / (sizeof(unsigned char) * 8)] >> (((sizeof(unsigned char) * 8) - 1) - ((ndid - n) % (sizeof(unsigned char) * 8)))) & 0x1)) {
        if (rid >= n) {
          ((rvisited)[(ndid - n) / (sizeof(unsigned char) * 8)] |= ((0x1) << (((sizeof(unsigned char) * 8) - 1) - ((ndid - n) % (sizeof(unsigned char) * 8)))));
          curNode[k + 1] = rid;
          k++;
          continue;
        } else if (rid < n) {
          nc++;
        }
      }
      k--;
    }
    if (thresh > max_illegal && nc > mc) {
      ;
      max_illegal = thresh;
      continue;
    }
    if (thresh < min_legal && nc <= mc) {
      min_legal = thresh;
      min_legal_nc = nc;
      ;
    }
  }
  form_flat_clusters_from_monotonic_criterion(Z, mono_crit, T, min_legal, n);
  free(curNode);
  free(lvisited);
  free(rvisited);
}
```