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

```asm
<ldumat>:
  endbr64
  push   %r15
  mov    %rsi,%r15
  mov    $0x8,%esi
  push   %r14
  mov    %ecx,%r14d
  push   %r13
  mov    %rdi,%r13
  mov    %edx,%edi
  push   %r12
  push   %rbp
  push   %rbx
  mov    %edx,%ebx
  mov    %ebx,%r12d
  imul   %ebx,%r12d
  sub    $0x18,%rsp
  call   <calloc@plt>
  mov    %rax,%rbp
  test   %r12d,%r12d
  je     L227
  movslq %r12d,%r12
  xor    %esi,%esi
  mov    %r15,%rdi
  lea    0x0(,%r12,8),%rdx
  call   <memset@plt>
L128:
  mov    %r14d,%eax
  mov    %ebx,%edi
  lea    -0x8(%r15,%r12,8),%r10
  imul   %r14d,%eax
  sub    %r14d,%edi
  lea    -0x1(%r14),%r12d
  cltq
  lea    -0x8(%r13,%rax,8),%r13
  test   %edi,%edi
  jle    L137
  movslq %ebx,%rax
  mov    $0xfffffffffffffff8,%rcx
  movsd  D25(%rip),%xmm6
  mov    %r10,%rdx
  shl    $0x3,%rax
  sub    %rax,%rcx
  xor    %eax,%eax
  nopl   0x0(%rax,%rax,1)
L47:
  add    $0x1,%eax
  movsd  %xmm6,(%rdx)
  add    %rcx,%rdx
  cmp    %eax,%edi
  jne    L47
  movslq %edi,%rax
  imul   %rcx,%rax
  add    %rax,%r10
L30:
  test   %r12d,%r12d
  js     L131
  lea    -0x1(%rbx),%ecx
  movslq %ebx,%rsi
  movslq %r14d,%r11
  mov    $0xfffffffffffffff8,%rax
  mov    %ecx,0xc(%rsp)
  lea    0x8(%rbp),%rcx
  shl    $0x3,%rsi
  mov    %rax,%r14
  mov    %rcx,(%rsp)
  shl    $0x3,%r11
  sub    %rsi,%rax
  add    $0x1,%r12d
  mov    %ebx,0x8(%rsp)
  sub    %r11,%r14
  sub    $0x1,%edi
  pxor   %xmm4,%xmm4
  movq   D4(%rip),%xmm5
  movsd  D25(%rip),%xmm6
  mov    %rax,%rbx
  movapd %xmm5,%xmm7
  nopl   0x0(%rax)
L182:
  movsd  0x0(%r13),%xmm3
  lea    0x8(%r10),%r8
  ucomisd %xmm4,%xmm3
  jp     L40
  je     L116
L40:
  lea    0x1(%rdi),%r9d
  lea    0x0(%r13,%r11,1),%rcx
  test   %r9d,%r9d
  jle    L59
  mov    (%rsp),%r15
  mov    %edi,%edx
  mov    %rbp,%rax
  lea    (%r15,%rdx,8),%r15
  mov    %rbp,%rdx
  nopl   0x0(%rax,%rax,1)
L95:
  movsd  (%rcx),%xmm0
  add    $0x8,%rdx
  add    %r11,%rcx
  movsd  %xmm0,-0x8(%rdx)
  cmp    %r15,%rdx
  jne    L95
  movsd  0x0(%r13),%xmm3
  movapd %xmm6,%xmm0
  lea    (%r10,%rsi,1),%rdx
  subsd  %xmm3,%xmm0
  movapd %xmm3,%xmm1
  xorpd  %xmm7,%xmm1
  movsd  %xmm0,(%r10)
  nopl   0x0(%rax)
L11:
  movsd  (%rax),%xmm0
  add    $0x8,%rax
  mulsd  %xmm1,%xmm0
  movsd  %xmm0,(%rdx)
  add    %rsi,%rdx
  cmp    %r15,%rax
  jne    L11
L61:
  cmp    %r12d,0x8(%rsp)
  jle    L63
  mov    0xc(%rsp),%eax
  sub    %r12d,%eax
  lea    0x10(%r10,%rax,8),%r15
  mov    (%rsp),%rax
  lea    (%rax,%rdi,8),%rdi
  nop
L104:
  lea    (%r8,%rsi,1),%rax
  test   %r9d,%r9d
  jle    L228
  mov    %rbp,%rdx
  mov    %rax,%rcx
  movapd %xmm4,%xmm1
  nopw   0x0(%rax,%rax,1)
L214:
  movsd  (%rdx),%xmm0
  mulsd  (%rcx),%xmm0
  add    $0x8,%rdx
  add    %rsi,%rcx
  addsd  %xmm0,%xmm1
  cmp    %rdx,%rdi
  jne    L214
  mulsd  %xmm3,%xmm1
  mov    %rbp,%rdx
  nop
L7:
  movsd  (%rdx),%xmm2
  movsd  (%rax),%xmm0
  add    $0x8,%rdx
  mulsd  %xmm1,%xmm2
  subsd  %xmm2,%xmm0
  movsd  %xmm0,(%rax)
  add    %rsi,%rax
  cmp    %rdx,%rdi
  jne    L7
L157:
  add    $0x8,%r8
  xorpd  %xmm5,%xmm1
  movsd  %xmm1,-0x8(%r8)
  cmp    %r15,%r8
  jne    L104
L63:
  add    %r14,%r13
  add    %rbx,%r10
  mov    %r9d,%edi
  sub    $0x1,%r12d
  jne    L182
L131:
  add    $0x18,%rsp
  mov    %rbp,%rdi
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  jmp    1070 <free@plt>
  nopl   0x0(%rax)
L116:
  lea    0x1(%rdi),%r9d
  movsd  %xmm6,(%r10)
  lea    (%rsi,%r10,1),%rax
  test   %r9d,%r9d
  jle    L63
  lea    0x10(%r10,%rdi,8),%rdx
  nopw   0x0(%rax,%rax,1)
L238:
  add    $0x8,%r8
  movq   $0x0,-0x8(%r8)
  movq   $0x0,(%rax)
  add    %rsi,%rax
  cmp    %rdx,%r8
  jne    L238
  add    %r14,%r13
  add    %rbx,%r10
  mov    %r9d,%edi
  sub    $0x1,%r12d
  jne    L182
  jmp    L131
L228:
  movapd %xmm3,%xmm1
  mulsd  %xmm4,%xmm1
  jmp    L157
L137:
  jne    L30
  mov    $0xfffffffffffffff8,%rax
  movslq %r14d,%rdx
  lea    -0x2(%r14),%r12d
  mov    $0x1,%edi
  shl    $0x3,%rdx
  mov    %rax,%rcx
  sub    %rdx,%rcx
  movslq %ebx,%rdx
  add    %rcx,%r13
  shl    $0x3,%rdx
  mov    D25(%rip),%rcx
  sub    %rdx,%rax
  mov    %rcx,(%r10)
  add    %rax,%r10
  jmp    L30
L59:
  movapd %xmm6,%xmm0
  subsd  %xmm3,%xmm0
  movsd  %xmm0,(%r10)
  jmp    L61
L227:
  xor    %r12d,%r12d
  jmp    L128
```

tool_calls: name: unpack_data
arguments: data_label: D25
data_type: f64
name: unpack_data
arguments: data_label: D4
data_type: f32[2]

tool: D25: 1

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

assistant: ```c
void ldumat(double *a, double *u, int m, int n) {
  double *p0, *q0, *p, *q, *w;
  int i, j, k, mm;
  double s, h;
  w = (double *)calloc(m, sizeof(double));
  for (i = 0, mm = m * m, q = u; i < mm; ++i)
    *q++ = 0.;
  p0 = a + n * n - 1;
  q0 = u + m * m - 1;
  mm = m - n;
  i = n - 1;
  for (j = 0; j < mm; ++j, q0 -= m + 1)
    *q0 = 1.;
  if (mm == 0) {
    p0 -= n + 1;
    *q0 = 1.;
    q0 -= m + 1;
    --i;
    ++mm;
  }
  for (; i >= 0; --i, ++mm, p0 -= n + 1, q0 -= m + 1) {
    if (*p0 != 0.) {
      for (j = 0, p = p0 + n, h = 1.; j < mm; p += n)
        w[j++] = *p;
      h = *p0;
      *q0 = 1. - h;
      for (j = 0, q = q0 + m; j < mm; q += m)
        *q = -h * w[j++];
      for (k = i + 1, q = q0 + 1; k < m; ++k) {
        for (j = 0, p = q + m, s = 0.; j < mm; p += m)
          s += w[j++] * *p;
        s *= h;
        for (j = 0, p = q + m; j < mm; p += m)
          *p -= s * w[j++];
        *q++ = -s;
      }
    } else {
      *q0 = 1.;
      for (j = 0, p = q0 + 1, q = q0 + m; j < mm; ++j, q += m)
        *q = *p++ = 0.;
    }
  }
  free(w);
}
```