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

```asm
<trsrch>:
  endbr64
  push   %r15
  push   %r14
  push   %r13
  push   %r12
  push   %rbp
  push   %rbx
  sub    $0x98,%rsp
  mov    %rdi,0x30(%rsp)
  mov    %rsi,0x38(%rsp)
  mov    %rdx,%r12
  mov    %rcx,%r15
  mov    %r8d,%ebx
  mov    %r9,%rbp
  movsd  %xmm0,0x20(%rsp)
  mov    0xd8(%rsp),%rax
  mov    %rax,0x70(%rsp)
  mov    0xe0(%rsp),%rax
  mov    %rax,0x78(%rsp)
  movsd  %xmm1,0x40(%rsp)
  mov    0xe8(%rsp),%r14
  movsd  %xmm2,0x48(%rsp)
  mov    0xf0(%rsp),%rax
  mov    %rax,0x50(%rsp)
  mov    %fs:0x28,%rax
  mov    %rax,0x88(%rsp)
  xor    %eax,%eax
  mov    D13(%rip),%rax
  mov    %rax,0x38(%r14)
  movslq %r8d,%rax
  mov    %rax,0x10(%rsp)
  lea    0x0(,%r8,8),%r13d
  mov    %r13d,%edi
  call   <malloc@plt>
  mov    %rax,0x68(%rsp)
  mov    %r13d,%edi
  call   <malloc@plt>
  mov    %rax,0x8(%rsp)
  mov    %r13d,%edi
  call   <malloc@plt>
  mov    %rax,0x18(%rsp)
  movq   $0x0,0x80(%rsp)
  test   %ebx,%ebx
  jle    L144
  mov    %ebx,%edx
  mov    $0x0,%eax
  pxor   %xmm3,%xmm3
L158:
  movsd  0x0(%rbp,%rax,8),%xmm2
  movsd  (%r15,%rax,8),%xmm1
  movapd %xmm2,%xmm0
  mulsd  %xmm1,%xmm0
  mulsd  %xmm2,%xmm0
  mulsd  %xmm1,%xmm0
  addsd  %xmm0,%xmm3
  add    $0x1,%rax
  cmp    %rdx,%rax
  jne    L158
  pxor   %xmm0,%xmm0
  ucomisd %xmm3,%xmm0
  ja     L193
L60:
  sqrtsd %xmm3,%xmm3
  movsd  %xmm3,0x28(%rsp)
  cmpl   $0x1,0xd0(%rsp)
  je     L104
L185:
  movsd  D18(%rip),%xmm0
  ucomisd (%r14),%xmm0
  jp     L18
  jne    L18
L104:
  movq   $0x0,0x8(%r14)
  test   %ebx,%ebx
  jle    L52
L154:
  mov    %ebx,%eax
  mov    0x68(%rsp),%rdi
  lea    (%rdi,%rax,8),%rsi
  mov    %rdi,%rax
  mov    %r12,%rcx
  sub    %rdi,%rcx
  mov    %rbp,%rdx
  sub    %rdi,%rdx
L177:
  movsd  (%rcx,%rax,1),%xmm0
  divsd  (%rdx,%rax,1),%xmm0
  movsd  %xmm0,(%rax)
  add    $0x8,%rax
  cmp    %rsi,%rax
  jne    L177
  mov    %ebx,%esi
  mov    0x68(%rsp),%rdi
  call   <l2norm@plt>
  mulsd  %xmm0,%xmm0
  movsd  %xmm0,0x58(%rsp)
  mov    %ebx,%r8d
  mov    0x10(%rsp),%rcx
  shl    $0x3,%rcx
  mov    $0x0,%edi
  mov    $0x0,%esi
  pxor   %xmm6,%xmm6
  movsd  %xmm6,0x10(%rsp)
  mov    0x70(%rsp),%r9
  jmp    L15
L144:
  pxor   %xmm3,%xmm3
  jmp    L60
L193:
  movapd %xmm3,%xmm0
  call   <sqrt@plt>
  movsd  %xmm0,0x28(%rsp)
  cmpl   $0x1,0xd0(%rsp)
  jne    L185
  movq   $0x0,0x8(%r14)
  jmp    L154
L63:
  movslq %edi,%rax
  add    %rsi,%rax
  lea    (%r9,%rax,8),%rdx
  mov    %rsi,%rax
  pxor   %xmm2,%xmm2
L137:
  movsd  0x0(%rbp,%rax,8),%xmm1
  movsd  (%rdx),%xmm0
  mulsd  (%r12,%rax,8),%xmm0
  mulsd  %xmm1,%xmm1
  divsd  %xmm1,%xmm0
  addsd  %xmm0,%xmm2
  add    $0x1,%rax
  add    %rcx,%rdx
  cmp    %eax,%ebx
  jg     L137
L110:
  mulsd  %xmm2,%xmm2
  addsd  0x10(%rsp),%xmm2
  movsd  %xmm2,0x10(%rsp)
  add    $0x1,%rsi
  add    %ebx,%edi
  cmp    %r8,%rsi
  je     L79
L15:
  pxor   %xmm2,%xmm2
  cmp    %esi,%ebx
  jg     L63
  jmp    L110
L16:
  movsd  0x58(%rsp),%xmm0
  call   <sqrt@plt>
  jmp    L14
L52:
  mov    %ebx,%esi
  mov    0x68(%rsp),%rdi
  call   <l2norm@plt>
  mulsd  %xmm0,%xmm0
  movsd  %xmm0,0x58(%rsp)
  pxor   %xmm7,%xmm7
  movsd  %xmm7,0x10(%rsp)
L79:
  pxor   %xmm0,%xmm0
  movsd  0x58(%rsp),%xmm5
  ucomisd %xmm5,%xmm0
  ja     L16
  sqrtsd %xmm5,%xmm5
  movapd %xmm5,%xmm0
L14:
  mulsd  0x58(%rsp),%xmm0
  divsd  0x10(%rsp),%xmm0
  movsd  0x20(%rsp),%xmm1
  minsd  %xmm0,%xmm1
  movsd  %xmm1,(%r14)
L18:
  mov    $0x4,%r13d
  lea    0x80(%rsp),%rax
  mov    %rax,0x58(%rsp)
  mov    %r15,0x60(%rsp)
  mov    %ebx,0x10(%rsp)
  mov    0x70(%rsp),%rbx
  mov    0x78(%rsp),%r15
L109:
  push   0x8(%rsp)
  push   %r14
  movsd  0x58(%rsp),%xmm1
  movsd  0x38(%rsp),%xmm0
  mov    %r15,%r9
  mov    %rbx,%r8
  mov    %rbp,%rcx
  mov    0x20(%rsp),%edx
  mov    0x70(%rsp),%rsi
  mov    %r12,%rdi
  call   <trstep@plt>
  movsd  (%r14),%xmm0
  movsd  %xmm0,0x10(%r14)
  push   0x60(%rsp)
  push   0x70(%rsp)
  push   0x38(%rsp)
  push   %r14
  push   $0x0
  push   %r15
  push   %rbx
  push   %r13
  movsd  0x90(%rsp),%xmm1
  movsd  0x70(%rsp),%xmm0
  mov    %rbp,%r9
  mov    0x60(%rsp),%r8d
  mov    0x58(%rsp),%rcx
  mov    %r12,%rdx
  mov    0x88(%rsp),%rsi
  mov    0x80(%rsp),%rdi
  call   <trupdate@plt>
  mov    %eax,%r13d
  add    $0x50,%rsp
  cmp    $0x1,%eax
  jle    L44
  cmp    $0xf,%eax
  jne    L109
L44:
  mov    0x68(%rsp),%rdi
  call   <free@plt>
  mov    0x8(%rsp),%rdi
  call   <free@plt>
  mov    0x18(%rsp),%rdi
  call   <free@plt>
  mov    0x88(%rsp),%rax
  sub    %fs:0x28,%rax
  jne    L72
  mov    %r13d,%eax
  add    $0x98,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
L72:
  call   <__stack_chk_fail@plt>
```

tool_calls: name: unpack_data
arguments: tag: D13
type_infered: f64
name: unpack_data
arguments: tag: D18
type_infered: f32[2]

tool: D13: 1

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

assistant: ```c
int trsrch(double (*funcpt)(double *, int), double *xi, double *jac, double *sN, int N, double *dx, double maxstep, int iter, double *L, double *hess, double stol, double *ioval, double eps, double *x) {
  int retval, i, j, method;
  double alpha, beta;
  double nlen, t1, funcprev;
  double *temp, *step, *xprev;
  retval = 4;
  nlen = 0.0;
  ioval[7] = 1.0;
  method = 0;
  temp = (double *)malloc(sizeof(double) * N);
  step = (double *)malloc(sizeof(double) * N);
  xprev = (double *)malloc(sizeof(double) * N);
  funcprev = 0.0;
  for (i = 0; i < N; ++i) {
    nlen += dx[i] * sN[i] * dx[i] * sN[i];
  }
  nlen = sqrt(nlen);
  if (iter == 1 || ioval[0] == -1) {
    ioval[1] = 0.;
    for (i = 0; i < N; ++i) {
      temp[i] = jac[i] / dx[i];
    }
    alpha = l2norm(temp, N);
    alpha = alpha * alpha;
    beta = 0;
    for (i = 0; i < N; ++i) {
      t1 = 0.;
      for (j = i; j < N; ++j) {
        t1 += L[j * N + i] * jac[j] / (dx[j] * dx[j]);
      }
      beta += t1 * t1;
    }
    ioval[0] = alpha * sqrt(alpha) / beta;
    if (ioval[0] > maxstep) {
      ioval[0] = maxstep;
    }
  }
  while (retval >= 2 && retval != 15) {
    trstep(jac, sN, N, dx, L, hess, nlen, ioval, eps, step);
    ioval[2] = ioval[0];
    retval = trupdate(funcpt, xi, jac, step, N, dx, maxstep, retval, L, hess, stol, method, ioval, xprev, &funcprev, x);
  }
  free(temp);
  free(step);
  free(xprev);
  return retval;
}
```