user: What would the C source look like for the assembly code provided here?

```asm
<matupd_>:
  endbr64
  push   %r15
  mov    %rdx,%r10
  lea    0x1(%rsi),%edx
  push   %r14
  movslq %edx,%rdx
  push   %r13
  push   %r12
  push   %rbp
  mov    %edi,%ebp
  lea    0x0(,%rdx,8),%rdi
  push   %rbx
  mov    %esi,%ebx
  mov    %r9,%rsi
  lea    0x1(%rbp),%edx
  sub    %rdi,%rsi
  movslq %edx,%rdx
  shl    $0x3,%rdx
  sub    $0x68,%rsp
  sub    %rdx,%rcx
  mov    %rsi,0x20(%rsp)
  mov    %r8,%rsi
  mov    0xb8(%rsp),%eax
  sub    %rdi,%rsi
  mov    %rdi,0x8(%rsp)
  mov    0xa0(%rsp),%r14
  mov    %rsi,0x28(%rsp)
  mov    %r10,%rsi
  mov    0xb0(%rsp),%r13
  sub    %rdx,%rsi
  mov    %rcx,0x30(%rsp)
  mov    0xc0(%rsp),%r12
  mov    %rsi,0x38(%rsp)
  movsd  %xmm1,0x40(%rsp)
  movsd  %xmm2,0x50(%rsp)
  movsd  %xmm3,0x48(%rsp)
  cmp    %eax,%ebx
  jl     L78
  mov    0xc8(%rsp),%rdi
  mov    %eax,(%r12)
  movsd  %xmm0,0x8(%rsp)
  add    (%rdi),%eax
  mov    %ebp,%edi
  sub    $0x2,%eax
  cltd
  idiv   %ebx
  add    $0x1,%edx
  mov    %edx,0x0(%r13)
  imul   %ebp,%edx
  movslq %edx,%rdx
  lea    0x8(%rsi,%rdx,8),%rdx
  mov    %r14,%rsi
  call   <dcopy_@plt>
  mov    0x0(%r13),%eax
  mov    0x30(%rsp),%rcx
  mov    %ebp,%edi
  mov    0xa8(%rsp),%rsi
  imul   %ebp,%eax
  cltq
  lea    0x8(%rcx,%rax,8),%rdx
  call   <dcopy_@plt>
  mov    0xd0(%rsp),%rax
  movsd  0x8(%rsp),%xmm0
  divsd  0x40(%rsp),%xmm0
  movsd  %xmm0,(%rax)
L22:
  mov    0xc8(%rsp),%rax
  mov    (%r12),%edx
  mov    (%rax),%eax
  mov    %edx,0x18(%rsp)
  cmp    $0x1,%edx
  jle    L51
  mov    %r12,0xc0(%rsp)
  mov    $0x1,%r15d
  mov    %ebp,%r12d
  mov    %eax,%r13d
  mov    %ebx,%ebp
  jmp    L93
  nopl   0x0(%rax,%rax,1)
L69:
  mov    0xc0(%rsp),%rax
  mov    (%rax),%edx
L93:
  mov    %r12d,%ecx
  mov    0x28(%rsp),%rax
  add    %ebp,%edx
  mov    %r14,%rsi
  imul   %r13d,%ecx
  movslq %edx,%rdx
  mov    %r12d,%edi
  add    %ebx,%ebp
  movslq %ecx,%rcx
  lea    0x8(,%rcx,8),%r8
  lea    (%rax,%rdx,8),%rcx
  mov    0x30(%rsp),%rax
  mov    %rcx,0x8(%rsp)
  lea    (%rax,%r8,1),%rdx
  mov    %r8,0x10(%rsp)
  call   <ddot_@plt>
  mov    0x8(%rsp),%rcx
  mov    %r14,%rdx
  mov    %r12d,%edi
  mov    0xc0(%rsp),%rax
  mov    0x10(%rsp),%r8
  movsd  %xmm0,(%rcx)
  mov    (%rax),%ecx
  mov    0x20(%rsp),%rax
  imul   %ebx,%ecx
  add    %r15d,%ecx
  add    $0x1,%r15d
  movslq %ecx,%rcx
  lea    (%rax,%rcx,8),%rcx
  mov    0x38(%rsp),%rax
  mov    %rcx,0x8(%rsp)
  lea    (%rax,%r8,1),%rsi
  call   <ddot_@plt>
  mov    %r13d,%eax
  mov    0x8(%rsp),%rcx
  cltd
  idiv   %ebx
  movsd  %xmm0,(%rcx)
  lea    0x1(%rdx),%r13d
  cmp    %r15d,0x18(%rsp)
  jne    L69
  mov    0xc0(%rsp),%r12
  mov    (%r12),%eax
  mov    %eax,0x18(%rsp)
L51:
  mov    0x18(%rsp),%eax
  movsd  0x50(%rsp),%xmm5
  imul   %eax,%ebx
  add    %eax,%ebx
  ucomisd D0(%rip),%xmm5
  mov    0x20(%rsp),%rax
  movslq %ebx,%rbx
  lea    (%rax,%rbx,8),%rax
  jp     L26
  jne    L26
L39:
  movsd  0x48(%rsp),%xmm6
  movsd  0x40(%rsp),%xmm7
  movsd  %xmm6,(%rax)
  mov    0x28(%rsp),%rax
  movsd  %xmm7,(%rax,%rbx,8)
  add    $0x68,%rsp
  xor    %eax,%eax
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
  nopl   0x0(%rax)
L78:
  mov    0x0(%r13),%eax
  mov    %ebp,%edi
  mov    %r9,0x10(%rsp)
  mov    %r8,%r15
  movsd  %xmm0,0x18(%rsp)
  cltd
  idiv   %ebx
  mov    0xc8(%rsp),%rax
  add    $0x1,%edx
  mov    %edx,0x0(%r13)
  mov    (%rax),%eax
  cltd
  idiv   %ebx
  mov    0xc8(%rsp),%rax
  add    $0x1,%edx
  mov    %edx,(%rax)
  mov    0x0(%r13),%eax
  imul   %ebp,%eax
  cltq
  lea    0x8(%rsi,%rax,8),%rdx
  mov    %r14,%rsi
  call   <dcopy_@plt>
  mov    0x0(%r13),%eax
  mov    0x30(%rsp),%rdi
  mov    0xa8(%rsp),%rsi
  imul   %ebp,%eax
  cltq
  lea    0x8(%rdi,%rax,8),%rdx
  mov    %ebp,%edi
  call   <dcopy_@plt>
  movsd  0x18(%rsp),%xmm0
  divsd  0x40(%rsp),%xmm0
  mov    0xd0(%rsp),%rax
  mov    0x10(%rsp),%r9
  movsd  %xmm0,(%rax)
  mov    (%r12),%eax
  cmp    $0x1,%eax
  mov    %eax,0x18(%rsp)
  jle    L51
  mov    0x8(%rsp),%rax
  mov    $0x1,%ecx
  mov    %ebp,0x10(%rsp)
  mov    %ebx,0x5c(%rsp)
  mov    %ecx,%ebx
  lea    -0x8(%rax),%rdi
  mov    %r14,0xa0(%rsp)
  mov    %rax,%rbp
  mov    %r15,%r14
  mov    %rdi,0x8(%rsp)
  mov    %r9,%r15
  nopl   0x0(%rax)
L40:
  mov    %ebx,%r13d
  lea    0x0(%rbp,%r15,1),%rsi
  mov    %r15,%rdx
  add    $0x1,%ebx
  mov    %r13d,%edi
  call   <dcopy_@plt>
  mov    (%r12),%edi
  mov    %r14,%rdx
  add    %rbp,%r14
  mov    %r14,%rsi
  sub    %r13d,%edi
  call   <dcopy_@plt>
  add    0x8(%rsp),%r15
  cmp    %ebx,0x18(%rsp)
  jne    L40
  mov    0x10(%rsp),%ebp
  mov    0x5c(%rsp),%ebx
  mov    0xa0(%rsp),%r14
  jmp    L22
  nopl   0x0(%rax,%rax,1)
L26:
  movsd  0x50(%rsp),%xmm0
  mulsd  %xmm0,%xmm0
  mulsd  0x48(%rsp),%xmm0
  movsd  %xmm0,0x48(%rsp)
  jmp    L39
```

tool_calls: name: disassemble_data
arguments: label: D0
type: f64

tool: D0: 1

assistant: ```c
int matupd_(int n, int m, double *ws, double *wy, double *sy, double *ss, double *d, double *r, int *itail, int iupdat, int *col, int *head, double *theta, double rr, double dr, double stp, double dtd) {
  int ws_offset, wy_offset, sy_offset, ss_offset, i1, i2;
  int j, pointr;
  --r;
  --d;
  ss_offset = 1 + m;
  ss -= ss_offset;
  sy_offset = 1 + m;
  sy -= sy_offset;
  wy_offset = 1 + n;
  wy -= wy_offset;
  ws_offset = 1 + n;
  ws -= ws_offset;
  if (iupdat <= m) {
    *col = iupdat;
    *itail = (*head + iupdat - 2) % m + 1;
  } else {
    *itail = *itail % m + 1;
    *head = *head % m + 1;
  }
  dcopy_(n, &d[1], &ws[*itail * n + 1]);
  dcopy_(n, &r[1], &wy[*itail * n + 1]);
  *theta = rr / dr;
  if (iupdat > m) {
    i1 = *col - 1;
    for (j = 1; j <= i1; ++j) {
      dcopy_(j, &ss[(j + 1) * m + 2], &ss[j * m + 1]);
      i2 = *col - j;
      dcopy_(i2, &sy[j + 1 + (j + 1) * m], &sy[j + j * m]);
    }
  }
  pointr = *head;
  i1 = *col - 1;
  for (j = 1; j <= i1; ++j) {
    sy[*col + j * m] = ddot_(n, &d[1], &wy[pointr * n + 1]);
    ss[j + *col * m] = ddot_(n, &ws[pointr * n + 1], &d[1]);
    pointr = pointr % m + 1;
  }
  if (stp == 1.) {
    ss[*col + *col * m] = dtd;
  } else {
    ss[*col + *col * m] = stp * stp * dtd;
  }
  sy[*col + *col * m] = dr;
  return 0;
}
```