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

```asm
<fitLine>:
  endbr64
  push   %r15
  xor    %r15d,%r15d
  push   %r14
  push   %r13
  push   %r12
  push   %rbp
  push   %rbx
  sub    $0xc8,%rsp
  mov    %rdx,0x68(%rsp)
  mov    %rcx,0x70(%rsp)
  mov    %r8d,0x58(%rsp)
  mov    %r9d,0x1c(%rsp)
  mov    %fs:0x28,%rax
  mov    %rax,0xb8(%rsp)
  xor    %eax,%eax
  test   %esi,%esi
  jle    L92
  mov    %rdi,%r12
  mov    %esi,%ebx
  cmp    $0x1,%esi
  je     L142
  pxor   %xmm3,%xmm3
  lea    -0x1(%rsi),%r15d
  movsd  D27(%rip),%xmm1
  movslq %esi,%r13
  cvtsd2ss %xmm0,%xmm3
  pxor   %xmm0,%xmm0
  movss  %xmm3,0x8(%rsp)
  lea    0x0(,%r13,4),%rbp
  cvtsi2sd %r15d,%xmm0
  mov    %rbp,%rdi
  divsd  %xmm0,%xmm1
  movsd  %xmm1,0x60(%rsp)
  call   <malloc@plt>
  mov    %rbp,%rdi
  mov    %rax,%r14
  call   <malloc@plt>
  mov    %r13,%rdi
  mov    $0x1,%esi
  mov    %rax,%rbp
  call   <calloc@plt>
  cmp    $0x2,%r15d
  movss  0x8(%rsp),%xmm3
  mov    %rax,%r13
  jbe    L15
  mov    %ebx,%edx
  movsd  0x60(%rsp),%xmm5
  movdqa D73(%rip),%xmm2
  mov    %rbp,%rax
  shr    $0x2,%edx
  movdqa D41(%rip),%xmm6
  movapd D29(%rip),%xmm4
  shl    $0x4,%rdx
  unpcklpd %xmm5,%xmm5
  add    %rbp,%rdx
  nopl   0x0(%rax)
L174:
  movdqa %xmm2,%xmm0
  add    $0x10,%rax
  paddd  %xmm6,%xmm2
  cvtdq2pd %xmm0,%xmm1
  mulpd  %xmm5,%xmm1
  pshufd $0xee,%xmm0,%xmm0
  cvtdq2pd %xmm0,%xmm0
  mulpd  %xmm5,%xmm0
  addpd  %xmm4,%xmm1
  addpd  %xmm4,%xmm0
  cvtpd2ps %xmm1,%xmm1
  cvtpd2ps %xmm0,%xmm0
  movlhps %xmm0,%xmm1
  movups %xmm1,-0x10(%rax)
  cmp    %rax,%rdx
  jne    L174
  mov    %ebx,%eax
  pxor   %xmm4,%xmm4
  and    $0xfffffffc,%eax
  test   $0x3,%bl
  je     L164
  pxor   %xmm0,%xmm0
  movsd  D69(%rip),%xmm1
  movslq %eax,%rcx
  pxor   %xmm4,%xmm4
  cvtsi2sd %eax,%xmm0
  mulsd  0x60(%rsp),%xmm0
  add    $0x1,%eax
  lea    0x0(,%rcx,4),%rdx
  subsd  %xmm1,%xmm0
  cvtsd2ss %xmm0,%xmm0
  movss  %xmm0,0x0(%rbp,%rcx,4)
  cmp    %eax,%ebx
  jle    L164
L95:
  pxor   %xmm0,%xmm0
  movsd  0x60(%rsp),%xmm7
  cvtsi2sd %eax,%xmm0
  add    $0x1,%eax
  mulsd  %xmm7,%xmm0
  subsd  %xmm1,%xmm0
  cvtsd2ss %xmm0,%xmm0
  movss  %xmm0,0x4(%rbp,%rdx,1)
  cmp    %ebx,%eax
  jge    L164
  pxor   %xmm0,%xmm0
  cvtsi2sd %eax,%xmm0
  mulsd  %xmm7,%xmm0
  subsd  %xmm1,%xmm0
  cvtsd2ss %xmm0,%xmm0
  movss  %xmm0,0x8(%rbp,%rdx,1)
L164:
  mov    %ebx,%edx
  movapd %xmm4,%xmm2
  movapd %xmm4,%xmm5
  xor    %eax,%eax
  movq   $0x0,0xb0(%rsp)
  cs nopw 0x0(%rax,%rax,1)
L28:
  pxor   %xmm0,%xmm0
  pxor   %xmm1,%xmm1
  cvtss2sd 0x0(%rbp,%rax,4),%xmm0
  movapd %xmm0,%xmm6
  cvtss2sd (%r12,%rax,4),%xmm1
  add    $0x1,%rax
  mulsd  %xmm0,%xmm6
  addsd  %xmm1,%xmm2
  mulsd  %xmm1,%xmm0
  addsd  %xmm6,%xmm5
  addsd  %xmm0,%xmm4
  cmp    %rdx,%rax
  jne    L28
  movsd  %xmm4,0xa0(%rsp)
  divsd  %xmm5,%xmm4
  pxor   %xmm0,%xmm0
  mov    %ebx,%r15d
  cvtsi2sd %ebx,%xmm0
  mov    $0x5,%eax
  movsd  %xmm2,0xa8(%rsp)
  movsd  %xmm5,0x98(%rsp)
  divsd  %xmm0,%xmm2
  mulsd  D43(%rip),%xmm0
  cvttsd2si %xmm0,%esi
  movapd %xmm4,%xmm1
  movsd  %xmm4,0x78(%rsp)
  cmp    %eax,%esi
  cmovge %esi,%eax
  mov    %eax,0x5c(%rsp)
  mov    0x1c(%rsp),%eax
  test   %eax,%eax
  jle    L84
  lea    0x90(%rsp),%rax
  xor    %r11d,%r11d
  pxor   %xmm7,%xmm7
  mov    %ebx,%r9d
  mov    %rax,0x20(%rsp)
  lea    0x88(%rsp),%rax
  cvtss2sd %xmm3,%xmm7
  mov    %r11d,%r15d
  mov    %rax,0x50(%rsp)
  lea    0x98(%rsp),%rax
  mov    %rax,0x28(%rsp)
  lea    0xa8(%rsp),%rax
  mov    %rax,0x38(%rsp)
  lea    0xb0(%rsp),%rax
  mov    %rax,0x40(%rsp)
  lea    0xa0(%rsp),%rax
  mov    %rax,0x30(%rsp)
  movsd  %xmm7,0x48(%rsp)
  jmp    L119
  nopl   (%rax)
L19:
  movsd  0xb0(%rsp),%xmm0
  movsd  0x98(%rsp),%xmm5
  movsd  0xa0(%rsp),%xmm1
  movsd  0xa8(%rsp),%xmm2
  movapd %xmm0,%xmm4
  divsd  %xmm5,%xmm4
  movapd %xmm1,%xmm3
  mulsd  %xmm4,%xmm3
  mulsd  %xmm0,%xmm4
  subsd  %xmm3,%xmm2
  pxor   %xmm3,%xmm3
  cvtsi2sd %eax,%xmm3
  subsd  %xmm4,%xmm3
  divsd  %xmm3,%xmm2
  mulsd  %xmm2,%xmm0
  subsd  %xmm0,%xmm1
  divsd  %xmm5,%xmm1
  cmp    %eax,0x18(%rsp)
  jle    L123
  cmp    %eax,0x5c(%rsp)
  jg     L123
  add    $0x1,%r15d
  cmp    %r15d,0x1c(%rsp)
  je     L123
L119:
  movapd %xmm2,%xmm0
  mov    %ebx,%ecx
  mov    %rbp,%rdx
  mov    %r14,%rsi
  mov    %r12,%rdi
  mov    $0x2,%eax
  mov    %r9d,0x18(%rsp)
  movsd  %xmm1,0x10(%rsp)
  movsd  %xmm2,0x8(%rsp)
  call   <flattenData@plt>
  mov    0x20(%rsp),%r8
  mov    %ebx,%edx
  mov    %r13,%rsi
  mov    0x50(%rsp),%rcx
  mov    %r14,%rdi
  xor    %eax,%eax
  call   <computeSigma@plt>
  mov    0x58(%rsp),%eax
  mov    %ebx,%r8d
  mov    %r13,%rcx
  movsd  0x48(%rsp),%xmm0
  mov    %rbp,%rdx
  mov    %r14,%rsi
  mov    %r12,%rdi
  mulsd  0x90(%rsp),%xmm0
  push   %rax
  mov    $0x1,%eax
  push   0x40(%rsp)
  push   0x50(%rsp)
  push   0x48(%rsp)
  mov    0x48(%rsp),%r9
  call   <rejectPixels@plt>
  add    $0x20,%rsp
  movsd  0x8(%rsp),%xmm2
  movsd  0x10(%rsp),%xmm1
  test   %eax,%eax
  mov    %eax,%r9d
  jg     L19
L123:
  mov    %r9d,%r15d
L84:
  subsd  %xmm1,%xmm2
  mulsd  0x60(%rsp),%xmm1
  mov    0x68(%rsp),%rax
  cvtsd2ss %xmm2,%xmm2
  movss  %xmm2,(%rax)
  cvtsd2ss %xmm1,%xmm1
  comiss D73(%rip),%xmm1
  movaps %xmm1,%xmm0
  jae    L76
  xorps  D26(%rip),%xmm0
L76:
  movsd  D22(%rip),%xmm2
  cvtss2sd %xmm0,%xmm0
  comisd %xmm0,%xmm2
  ja     L152
  mov    0x70(%rsp),%rax
  movss  %xmm1,(%rax)
L56:
  mov    %r14,%rdi
  call   <free@plt>
  mov    %rbp,%rdi
  call   <free@plt>
  mov    %r13,%rdi
  call   <free@plt>
L92:
  mov    0xb8(%rsp),%rax
  sub    %fs:0x28,%rax
  jne    L66
  add    $0xc8,%rsp
  mov    %r15d,%eax
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
  nopw   0x0(%rax,%rax,1)
L142:
  movss  0x4(%rdi),%xmm0
  mov    $0x1,%r15d
  movss  %xmm0,(%rdx)
  movl   $0x0,(%rcx)
  jmp    L92
  nopw   0x0(%rax,%rax,1)
L152:
  movsd  0x60(%rsp),%xmm0
  mulsd  0x78(%rsp),%xmm0
  mov    0x70(%rsp),%rax
  cvtsd2ss %xmm0,%xmm0
  movss  %xmm0,(%rax)
  jmp    L56
  xchg   %ax,%ax
L15:
  movsd  0x60(%rsp),%xmm0
  pxor   %xmm4,%xmm4
  mov    $0x1,%eax
  xor    %edx,%edx
  movsd  D69(%rip),%xmm1
  mulsd  %xmm4,%xmm0
  subsd  %xmm1,%xmm0
  cvtsd2ss %xmm0,%xmm0
  movss  %xmm0,0x0(%rbp)
  jmp    L95
L66:
  call   <__stack_chk_fail@plt>
```

tool_calls: name: extract_symbolic_data
arguments: label: D27
type_guessed: f64
name: extract_symbolic_data
arguments: label: D73
type_guessed: f32[4]
name: extract_symbolic_data
arguments: label: D41
type_guessed: byte[16]
name: extract_symbolic_data
arguments: label: D29
type_guessed: f32[4]
name: extract_symbolic_data
arguments: label: D69
type_guessed: f64
name: extract_symbolic_data
arguments: label: D43
type_guessed: f64
name: extract_symbolic_data
arguments: label: D26
type_guessed: byte[16]
name: extract_symbolic_data
arguments: label: D22
type_guessed: f64

tool: D27: 2

tool: {"D73": [0.0,1.401298464324817e-45,2.802596928649634e-45,4.203895392974451e-45]}

tool: {"D41": [0x04,0x00,0x00,0x00,0x04,0x00,0x00,0x00,0x04,0x00,0x00,0x00,0x04,0x00,0x00,0x00]}

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

tool: D69: 1

tool: D43: 0.5

tool: {"D26": [0x00,0x00,0x00,0x80,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00]}

tool: D22: 0.001

assistant: ```c
int fitLine(data, npix, zstart, zslope, krej, ngrow, maxiter)
float *data;
int npix;
float *zstart;
float *zslope;
float krej;
int ngrow;
int maxiter;
{
  int i, ngoodpix, last_ngoodpix, minpix, niter;
  double xscale, z0, dz, o_dz, x, z, mean, sigma, threshold;
  double sumxsqr, sumxz, sumz, sumx, rowrat;
  float *flat, *normx;
  char *badpix;
  if (npix <= 0)
    return (0);
  else if (npix == 1) {
    *zstart = data[1];
    *zslope = 0.0;
    return (1);
  } else
    xscale = 2.0 / (npix - 1);
  flat = (float *)malloc(npix * sizeof(float));
  normx = (float *)malloc(npix * sizeof(float));
  badpix = (char *)calloc(npix, sizeof(char));
  for (i = 0; i < npix; i++)
    normx[i] = i * xscale - 1.0;
  sumxsqr = 0;
  sumxz = 0;
  sumx = 0;
  sumz = 0;
  for (i = 0; i < npix; i++) {
    x = normx[i];
    z = data[i];
    sumxsqr = sumxsqr + (x * x);
    sumxz = sumxz + z * x;
    sumz = sumz + z;
  }
  z0 = sumz / npix;
  dz = o_dz = sumxz / sumxsqr;
  ngoodpix = npix;
  minpix = ((5) > ((int)(npix * 0.5)) ? (5) : ((int)(npix * 0.5)));
  for (niter = 0; niter < maxiter; niter++) {
    last_ngoodpix = ngoodpix;
    flattenData(data, flat, normx, npix, z0, dz);
    ngoodpix = computeSigma(flat, badpix, npix, &mean, &sigma);
    threshold = sigma * krej;
    ngoodpix = rejectPixels(data, flat, normx, badpix, npix, &sumxsqr, &sumxz, &sumx, &sumz, threshold, ngrow);
    if (ngoodpix > 0) {
      rowrat = sumx / sumxsqr;
      z0 = (sumz - rowrat * sumxz) / (ngoodpix - rowrat * sumx);
      dz = (sumxz - z0 * sumx) / sumxsqr;
    }
    if (ngoodpix >= last_ngoodpix || ngoodpix < minpix)
      break;
  }
  *zstart = z0 - dz;
  *zslope = dz * xscale;
  if (((*zslope) >= 0 ? (*zslope) : -(*zslope)) < 0.001)
    *zslope = o_dz * xscale;
  free((float *)flat);
  free((float *)normx);
  free((char *)badpix);
  return (ngoodpix);
}
```