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

```asm
<ldumat>:
  endbr64
  push   %r15
  push   %r14
  push   %r13
  push   %r12
  push   %rbp
  push   %rbx
  sub    $0x28,%rsp
  mov    %rdi,%r13
  mov    %rsi,%r12
  mov    %edx,%ebx
  mov    %ecx,%ebp
  mov    $0x8,%esi
  mov    %edx,%edi
  call   <calloc@plt>
  mov    %rax,%r9
  mov    %ebx,%ecx
  imul   %ebx,%ecx
  test   %ecx,%ecx
  jle    L216
  mov    %ecx,%eax
  lea    (%r12,%rax,8),%rdx
  mov    %r12,%rax
L86:
  add    $0x8,%rax
  movq   $0x0,-0x8(%rax)
  cmp    %rdx,%rax
  jne    L86
L216:
  mov    %ebp,%eax
  imul   %ebp,%eax
  cltq
  lea    -0x8(%r13,%rax,8),%r14
  movslq %ecx,%rcx
  lea    -0x8(%r12,%rcx,8),%r11
  mov    %ebx,%r10d
  sub    %ebp,%r10d
  lea    -0x1(%rbp),%r13d
  test   %r10d,%r10d
  jle    L163
  movslq %ebx,%rax
  shl    $0x3,%rax
  mov    $0xfffffffffffffff8,%rcx
  sub    %rax,%rcx
  mov    %r11,%rdx
  mov    $0x0,%eax
  movsd  D2(%rip),%xmm0
L128:
  movsd  %xmm0,(%rdx)
  add    $0x1,%eax
  add    %rcx,%rdx
  cmp    %eax,%r10d
  jne    L128
  mov    %r10d,%eax
  imul   %rcx,%rax
  add    %rax,%r11
L74:
  test   %r13d,%r13d
  js     L224
  movslq %ebx,%rsi
  shl    $0x3,%rsi
  movslq %ebp,%rbp
  shl    $0x3,%rbp
  mov    $0xfffffffffffffff8,%rax
  mov    %rax,%rcx
  sub    %rbp,%rcx
  mov    %rcx,0x8(%rsp)
  sub    %rsi,%rax
  mov    %rax,0x10(%rsp)
  add    $0x1,%r13d
  lea    -0x1(%r10),%r15d
  pxor   %xmm6,%xmm6
  lea    -0x1(%rbx),%eax
  mov    %eax,0x1c(%rsp)
  movapd %xmm6,%xmm10
  movq   D40(%rip),%xmm4
  movsd  D2(%rip),%xmm9
  movapd %xmm4,%xmm8
  movapd %xmm6,%xmm5
  movapd %xmm9,%xmm7
  jmp    L132
L163:
  jne    L74
  movslq %ebp,%rdx
  shl    $0x3,%rdx
  mov    $0xfffffffffffffff8,%rax
  mov    %rax,%rcx
  sub    %rdx,%rcx
  add    %rcx,%r14
  mov    D2(%rip),%rcx
  mov    %rcx,(%r11)
  movslq %ebx,%rdx
  shl    $0x3,%rdx
  sub    %rdx,%rax
  add    %rax,%r11
  lea    -0x2(%rbp),%r13d
  mov    $0x1,%r10d
  jmp    L74
L243:
  lea    (%r14,%rbp,1),%rcx
  test   %r10d,%r10d
  jle    L192
  mov    %r9,%rax
  mov    %r15d,%edx
  lea    0x8(%r9,%rdx,8),%rdi
  mov    %r9,%rdx
L59:
  movsd  (%rcx),%xmm0
  movsd  %xmm0,(%rdx)
  add    %rbp,%rcx
  add    $0x8,%rdx
  cmp    %rdi,%rdx
  jne    L59
  movsd  (%r14),%xmm3
  movapd %xmm9,%xmm0
  subsd  %xmm3,%xmm0
  movsd  %xmm0,(%r11)
  lea    (%r11,%rsi,1),%rdx
  movapd %xmm3,%xmm1
  xorpd  %xmm8,%xmm1
L143:
  movapd %xmm1,%xmm0
  mulsd  (%rax),%xmm0
  movsd  %xmm0,(%rdx)
  add    %rsi,%rdx
  add    $0x8,%rax
  cmp    %rdi,%rax
  jne    L143
L256:
  lea    0x8(%r11),%r8
  cmp    %r13d,%ebx
  jle    L248
  mov    0x1c(%rsp),%eax
  sub    %r13d,%eax
  lea    0x10(%r11,%rax,8),%r12
  mov    %r15d,%eax
  lea    0x8(%r9,%rax,8),%rdi
  jmp    L236
L0:
  movapd %xmm1,%xmm2
  mulsd  (%rdx),%xmm2
  movsd  (%rax),%xmm0
  subsd  %xmm2,%xmm0
  movsd  %xmm0,(%rax)
  add    %rsi,%rax
  add    $0x8,%rdx
  cmp    %rdi,%rdx
  jne    L0
L26:
  add    $0x8,%r8
  xorpd  %xmm4,%xmm1
  movsd  %xmm1,-0x8(%r8)
  cmp    %r8,%r12
  je     L248
L236:
  lea    (%r8,%rsi,1),%rax
  test   %r10d,%r10d
  jle    L169
  mov    %r9,%rdx
  mov    %rax,%rcx
  pxor   %xmm1,%xmm1
L255:
  movsd  (%rdx),%xmm0
  mulsd  (%rcx),%xmm0
  addsd  %xmm0,%xmm1
  add    %rsi,%rcx
  add    $0x8,%rdx
  cmp    %rdi,%rdx
  jne    L255
  mulsd  %xmm3,%xmm1
  mov    %r9,%rdx
  jmp    L0
L248:
  add    $0x1,%r10d
  add    0x8(%rsp),%r14
  add    0x10(%rsp),%r11
  add    $0x1,%r15d
  sub    $0x1,%r13d
  je     L224
L132:
  movsd  (%r14),%xmm3
  ucomisd %xmm6,%xmm3
  jp     L243
  ucomisd %xmm5,%xmm3
  jne    L243
  movsd  %xmm7,(%r11)
  lea    0x8(%r11),%rax
  lea    (%rsi,%r11,1),%rdx
  test   %r10d,%r10d
  jle    L248
  mov    %r15d,%ecx
  lea    0x10(%r11,%rcx,8),%rcx
L178:
  add    $0x8,%rax
  movq   $0x0,-0x8(%rax)
  movq   $0x0,(%rdx)
  add    %rsi,%rdx
  cmp    %rcx,%rax
  jne    L178
  jmp    L248
L224:
  mov    %r9,%rdi
  call   <free@plt>
  add    $0x28,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
L169:
  movapd %xmm3,%xmm1
  mulsd  %xmm10,%xmm1
  jmp    L26
L192:
  movsd  D2(%rip),%xmm0
  subsd  %xmm3,%xmm0
  movsd  %xmm0,(%r11)
  jmp    L256
```

tool_calls: name: parse_data
arguments: label: D2
data_type: f64
name: parse_data
arguments: label: D40
data_type: f32[2]

tool: D2: 1

tool: {"D40": [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);
}
```