user: 请显示与下面列出的汇编代码相匹配的C语言实现：

```asm
<paw_test_triang_matrix_inverse>:
  endbr64
  push   %r13
  mov    $0x3,%esi
  mov    $0x3,%edi
  push   %r12
  push   %rbp
  push   %rbx
  sub    $0x8,%rsp
  call   <paw_alloc_2d_array@plt>
  mov    $0x3,%esi
  mov    $0x3,%edi
  mov    %rax,%rbx
  call   <paw_alloc_2d_array@plt>
  mov    $0x3,%esi
  mov    $0x3,%edi
  mov    %rax,%r12
  call   <paw_alloc_2d_array@plt>
  movapd D28(%rip),%xmm0
  mov    D63(%rip),%rdx
  mov    %r12,%rcx
  mov    %rax,%rbp
  mov    (%rbx),%rax
  mov    D20(%rip),%rsi
  lea    D45(%rip),%rdi
  lea    0x18(%rbp),%r13
  mov    %rdx,0x10(%rax)
  mov    D48(%rip),%rdx
  movups %xmm0,(%rax)
  mov    0x8(%rbx),%rax
  movapd D66(%rip),%xmm0
  mov    %rdx,0x10(%rax)
  mov    %rbx,%rdx
  movups %xmm0,(%rax)
  mov    0x10(%rbx),%rax
  pxor   %xmm0,%xmm0
  mov    %rsi,0x10(%rax)
  mov    $0x3,%esi
  movups %xmm0,(%rax)
  call   <paw_triang_matrix_inverse@plt>
  mov    (%r12),%rdi
  mov    (%rbx),%rcx
  pxor   %xmm0,%xmm0
  mov    0x0(%rbp),%rax
  mov    0x8(%r12),%rsi
  mov    0x10(%r12),%rdx
  lea    D21(%rip),%r12
  movq   $0x0,(%rax)
  movsd  (%rdi),%xmm2
  mulsd  (%rcx),%xmm2
  addsd  %xmm0,%xmm2
  movsd  %xmm2,(%rax)
  movsd  0x8(%rcx),%xmm1
  mulsd  (%rsi),%xmm1
  addsd  %xmm2,%xmm1
  movsd  %xmm1,(%rax)
  movsd  0x10(%rcx),%xmm2
  mulsd  (%rdx),%xmm2
  movq   $0x0,0x8(%rax)
  addsd  %xmm1,%xmm2
  movsd  %xmm2,(%rax)
  movsd  (%rcx),%xmm2
  mulsd  0x8(%rdi),%xmm2
  addsd  %xmm0,%xmm2
  movsd  %xmm2,0x8(%rax)
  movsd  0x8(%rcx),%xmm1
  mulsd  0x8(%rsi),%xmm1
  addsd  %xmm2,%xmm1
  movsd  %xmm1,0x8(%rax)
  movsd  0x10(%rcx),%xmm2
  mulsd  0x8(%rdx),%xmm2
  movq   $0x0,0x10(%rax)
  addsd  %xmm1,%xmm2
  movsd  %xmm2,0x8(%rax)
  movsd  (%rcx),%xmm2
  mulsd  0x10(%rdi),%xmm2
  addsd  %xmm0,%xmm2
  movsd  %xmm2,0x10(%rax)
  movsd  0x10(%rsi),%xmm1
  mulsd  0x8(%rcx),%xmm1
  addsd  %xmm2,%xmm1
  movsd  %xmm1,0x10(%rax)
  movsd  0x10(%rcx),%xmm2
  mulsd  0x10(%rdx),%xmm2
  mov    0x8(%rbx),%rcx
  addsd  %xmm1,%xmm2
  movsd  %xmm2,0x10(%rax)
  mov    0x8(%rbp),%rax
  movq   $0x0,(%rax)
  movsd  (%rcx),%xmm2
  mulsd  (%rdi),%xmm2
  addsd  %xmm0,%xmm2
  movsd  %xmm2,(%rax)
  movsd  0x8(%rcx),%xmm1
  mulsd  (%rsi),%xmm1
  addsd  %xmm2,%xmm1
  movsd  %xmm1,(%rax)
  movsd  0x10(%rcx),%xmm2
  mulsd  (%rdx),%xmm2
  movq   $0x0,0x8(%rax)
  addsd  %xmm1,%xmm2
  movsd  %xmm2,(%rax)
  movsd  (%rcx),%xmm2
  mulsd  0x8(%rdi),%xmm2
  addsd  %xmm0,%xmm2
  movsd  %xmm2,0x8(%rax)
  movsd  0x8(%rcx),%xmm1
  mulsd  0x8(%rsi),%xmm1
  addsd  %xmm2,%xmm1
  movsd  %xmm1,0x8(%rax)
  movsd  0x10(%rcx),%xmm2
  mulsd  0x8(%rdx),%xmm2
  movq   $0x0,0x10(%rax)
  addsd  %xmm1,%xmm2
  movsd  %xmm2,0x8(%rax)
  movsd  (%rcx),%xmm2
  mulsd  0x10(%rdi),%xmm2
  addsd  %xmm0,%xmm2
  movsd  %xmm2,0x10(%rax)
  movsd  0x8(%rcx),%xmm1
  mulsd  0x10(%rsi),%xmm1
  addsd  %xmm2,%xmm1
  movsd  %xmm1,0x10(%rax)
  movsd  0x10(%rcx),%xmm2
  mulsd  0x10(%rdx),%xmm2
  mov    0x10(%rbx),%rcx
  addsd  %xmm1,%xmm2
  movsd  %xmm2,0x10(%rax)
  mov    0x10(%rbp),%rax
  movq   $0x0,(%rax)
  movsd  (%rcx),%xmm2
  mulsd  (%rdi),%xmm2
  addsd  %xmm0,%xmm2
  movsd  %xmm2,(%rax)
  movsd  0x8(%rcx),%xmm1
  mulsd  (%rsi),%xmm1
  addsd  %xmm2,%xmm1
  movsd  %xmm1,(%rax)
  movsd  0x10(%rcx),%xmm2
  mulsd  (%rdx),%xmm2
  movq   $0x0,0x8(%rax)
  addsd  %xmm1,%xmm2
  movsd  %xmm2,(%rax)
  movsd  (%rcx),%xmm2
  mulsd  0x8(%rdi),%xmm2
  addsd  %xmm0,%xmm2
  movsd  %xmm2,0x8(%rax)
  movsd  0x8(%rcx),%xmm1
  mulsd  0x8(%rsi),%xmm1
  addsd  %xmm2,%xmm1
  movsd  %xmm1,0x8(%rax)
  movsd  0x10(%rcx),%xmm2
  mulsd  0x8(%rdx),%xmm2
  movq   $0x0,0x10(%rax)
  addsd  %xmm1,%xmm2
  movsd  %xmm2,0x8(%rax)
  movsd  (%rcx),%xmm1
  mulsd  0x10(%rdi),%xmm1
  addsd  %xmm0,%xmm1
  movsd  %xmm1,0x10(%rax)
  movsd  0x8(%rcx),%xmm0
  mulsd  0x10(%rsi),%xmm0
  addsd  %xmm1,%xmm0
  movsd  %xmm0,0x10(%rax)
  movsd  0x10(%rcx),%xmm1
  mulsd  0x10(%rdx),%xmm1
  addsd  %xmm0,%xmm1
  movsd  %xmm1,0x10(%rax)
L15:
  xor    %ebx,%ebx
L14:
  mov    0x0(%rbp),%rax
  mov    %r12,%rdi
  movsd  (%rax,%rbx,1),%xmm0
  mov    $0x1,%eax
  add    $0x8,%rbx
  call   <printf@plt>
  cmp    $0x18,%rbx
  jne    L14
  mov    $0xa,%edi
  add    $0x8,%rbp
  call   <putchar@plt>
  cmp    %rbp,%r13
  jne    L15
  add    $0x8,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  ret
```

tool_calls: name: extract_symbolic_data
arguments: data_label: D28
datatype: f64[2]
name: extract_symbolic_data
arguments: data_label: D63
datatype: f64
name: extract_symbolic_data
arguments: data_label: D20
datatype: f64
name: extract_symbolic_data
arguments: data_label: D45
datatype: string
name: extract_symbolic_data
arguments: data_label: D48
datatype: f64
name: extract_symbolic_data
arguments: data_label: D66
datatype: f64[2]
name: extract_symbolic_data
arguments: data_label: D21
datatype: string

tool: {"D28": [1.0,2.0]}

tool: D63: 3

tool: D20: 4

tool: D45: u

tool: D48: 6.3

tool: {"D66": [0.0,3.0]}

tool: D21: %f	

assistant: ```c
void paw_test_triang_matrix_inverse() {
  int i, j, k;
  int n;
  double **a;
  double **b;
  double **c;
  n = 3;
  a = paw_alloc_2d_array(n, n);
  b = paw_alloc_2d_array(n, n);
  c = paw_alloc_2d_array(n, n);
  a[0][0] = 1.0;
  a[0][1] = 2.0;
  a[0][2] = 3.0;
  a[1][0] = 0.0;
  a[1][1] = 3.0;
  a[1][2] = 6.3;
  a[2][0] = 0.0;
  a[2][1] = 0.0;
  a[2][2] = 4.0;
  paw_triang_matrix_inverse("u", n, a, b);
  for (i = 0; i < n; i++) {
    for (j = 0; j < n; j++) {
      c[i][j] = 0.0;
      for (k = 0; k < n; k++) {
        c[i][j] = c[i][j] + a[i][k] * b[k][j];
      }
    }
  }
  for (i = 0; i < n; i++) {
    for (j = 0; j < n; j++) {
      printf("%f\t", c[i][j]);
    }
    printf("\n");
  }
}
```