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

```asm
<lu>:
  endbr64
  push   %r15
  movsd  D12(%rip),%xmm1
  push   %r14
  mov    %edi,%r14d
  mov    %rsi,%rdi
  xor    %esi,%esi
  mov    %r14d,%eax
  push   %r13
  mov    %rdx,%r13
  imul   %r14d,%eax
  push   %r12
  push   %rbp
  push   %rbx
  movslq %eax,%r8
  test   %eax,%eax
  je     L63
  nop
L156:
  mov    %esi,%eax
  movsd  (%rdi,%rsi,8),%xmm0
  cltd
  idiv   %r14d
  movsd  %xmm0,(%rcx,%rsi,8)
  cmp    %eax,%edx
  je     L27
  movq   $0x0,0x0(%r13,%rsi,8)
  add    $0x1,%rsi
  cmp    %r8,%rsi
  jne    L156
L63:
  lea    -0x1(%r14),%eax
  mov    %eax,-0x10(%rsp)
  test   %eax,%eax
  jle    L11
  test   %r14d,%r14d
  jle    L11
  mov    %eax,-0x1c(%rsp)
  lea    0x1(%r14),%eax
  mov    %r14d,%r15d
  movq   D2(%rip),%xmm3
  mov    %eax,-0xc(%rsp)
  movslq %r14d,%rax
  mov    %rax,-0x50(%rsp)
  shl    $0x3,%rax
  mov    %rax,-0x48(%rsp)
  mov    %r14d,-0x18(%rsp)
  movl   $0x0,-0x14(%rsp)
  movl   $0x0,-0x20(%rsp)
  movq   $0x0,-0x28(%rsp)
  movl   $0x0,-0x34(%rsp)
  nopl   0x0(%rax)
L39:
  mov    -0x34(%rsp),%eax
  mov    %eax,%r10d
  add    $0x1,%eax
  mov    %eax,-0x34(%rsp)
  cmp    %eax,%r14d
  jle    L109
  movslq -0x14(%rsp),%rax
  mov    -0x28(%rsp),%r11
  movslq -0x18(%rsp),%rdx
  lea    (%rcx,%rax,8),%r12
  movslq -0x20(%rsp),%rax
  lea    (%r11,%rdx,1),%rsi
  mov    %rdx,%rdi
  mov    %rax,%rbp
  mov    %rax,%rbx
  lea    0x1(%rax,%r11,1),%rax
  shl    $0x3,%rsi
  shl    $0x3,%rax
  sub    %rdx,%rbx
  mov    %rax,-0x40(%rsp)
  lea    -0x8(%rcx,%rax,1),%r9
  mov    %r15d,%eax
  shr    %eax
  lea    -0x1(%rax),%r8d
  mov    %r15d,%eax
  and    $0xfffffffe,%eax
  add    $0x1,%r8
  add    %eax,%r10d
  mov    %eax,-0x38(%rsp)
  shl    $0x4,%r8
  lea    (%r10,%rbp,1),%eax
  mov    -0x34(%rsp),%ebp
  cltq
  lea    (%rcx,%rax,8),%r11
  mov    -0x1c(%rsp),%eax
  lea    0x8(%rcx,%rax,8),%rax
  mov    %rax,-0x30(%rsp)
  nopl   0x0(%rax,%rax,1)
L61:
  movsd  (%rcx,%rsi,1),%xmm1
  xorpd  %xmm3,%xmm1
  divsd  (%r12),%xmm1
  cmp    $0x1,%r15d
  je     L126
  cmp    -0x40(%rsp),%rsi
  je     L126
  movapd %xmm1,%xmm2
  lea    (%rcx,%rsi,1),%rdx
  xor    %eax,%eax
  unpcklpd %xmm2,%xmm2
  nopw   0x0(%rax,%rax,1)
L58:
  movupd (%r9,%rax,1),%xmm0
  movupd (%rdx,%rax,1),%xmm4
  mulpd  %xmm2,%xmm0
  addpd  %xmm4,%xmm0
  movups %xmm0,(%rdx,%rax,1)
  add    $0x10,%rax
  cmp    %r8,%rax
  jne    L58
  cmp    %r15d,-0x38(%rsp)
  je     L23
  movsd  (%r11),%xmm0
  lea    (%rdi,%r10,1),%eax
  cltq
  lea    (%rcx,%rax,8),%rax
  mulsd  %xmm1,%xmm0
  addsd  (%rax),%xmm0
  movsd  %xmm0,(%rax)
L23:
  xorpd  %xmm3,%xmm1
  add    $0x1,%ebp
  sub    -0x50(%rsp),%rbx
  add    %r14d,%edi
  movsd  %xmm1,0x0(%r13,%rsi,1)
  add    -0x48(%rsp),%rsi
  cmp    %ebp,%r14d
  jne    L61
L109:
  mov    -0xc(%rsp),%ebx
  addq   $0x1,-0x28(%rsp)
  sub    $0x1,%r15d
  add    %ebx,-0x14(%rsp)
  mov    -0x10(%rsp),%ebx
  add    %r14d,-0x20(%rsp)
  subl   $0x1,-0x1c(%rsp)
  add    %r14d,-0x18(%rsp)
  cmp    %ebx,-0x34(%rsp)
  jne    L39
L11:
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
  nopl   0x0(%rax,%rax,1)
L126:
  mov    -0x30(%rsp),%rdx
  lea    (%rcx,%rsi,1),%rax
  add    %rsi,%rdx
  nopl   0x0(%rax)
L12:
  movsd  (%rax,%rbx,8),%xmm0
  add    $0x8,%rax
  mulsd  %xmm1,%xmm0
  addsd  -0x8(%rax),%xmm0
  movsd  %xmm0,-0x8(%rax)
  cmp    %rax,%rdx
  jne    L12
  jmp    L23
L27:
  movsd  %xmm1,0x0(%r13,%rsi,8)
  add    $0x1,%rsi
  cmp    %r8,%rsi
  jne    L156
  jmp    L63
```

tool_calls: name: disassemble_data
arguments: tag: D12
type_guessed: byte[8]
name: disassemble_data
arguments: tag: D2
type_guessed: byte[8]

tool: {"D12": [0x00,0x00,0x00,0x00,0x00,0x00,0xf0,0x3f]}

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

assistant: ```c
void lu(int n, double *a, double *l, double *u) {
  int i, j, k;
  for (i = 0; i < n * n; i++) {
    u[i] = a[i];
    if ((i / n) == (i % n))
      l[i] = 1;
    else
      l[i] = 0;
  }
  for (j = 0; j < n - 1; j++) {
    for (i = j + 1; i < n; i++) {
      double m = -u[i * n + j] / u[j * n + j];
      for (k = j; k < n; k++)
        u[i * n + k] = (m * u[j * n + k]) + u[i * n + k];
      l[i * n + j] = -m;
    }
  }
}
```