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

```asm
<matrix_product>:
  endbr64
  push   %r15
  push   %r14
  push   %r13
  push   %r12
  push   %rbp
  push   %rbx
  sub    $0x28,%rsp
  mov    %rdi,%r14
  mov    %rsi,%rbx
  mov    %rdx,%rbp
  mov    %fs:0x28,%rax
  mov    %rax,0x18(%rsp)
  xor    %eax,%eax
  cmp    $0x4,%ecx
  sete   %al
  cmp    $0x4,%r8d
  jne    L31
  test   %al,%al
  jne    L85
L31:
  cmp    $0x9,%ecx
  sete   %dl
  cmp    $0x9,%r8d
  jne    L138
  test   %dl,%dl
  je     L138
  mov    $0x3,%r15d
  mov    $0x3,%r13d
  mov    $0x3,%r12d
  movl   $0x3,0xc(%rsp)
L29:
  cmp    %r13d,%r12d
  sete   %dil
  movzbl %dil,%edi
  call   <assert@plt>
  lea    0x14(%rsp),%rdx
  lea    0x10(%rsp),%rsi
  sub    $0x8,%rsp
  push   %r15
  push   %r13
  push   %rbp
  mov    %r12d,%r9d
  mov    0x2c(%rsp),%r8d
  mov    %rbx,%rcx
  mov    %r14,%rdi
  call   <matrix_product_clean@plt>
  mov    0x30(%rsp),%eax
  imul   0x34(%rsp),%eax
  add    $0x20,%rsp
L56:
  mov    0x18(%rsp),%rdx
  sub    %fs:0x28,%rdx
  jne    L98
  add    $0x28,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
L138:
  cmp    $0x10,%ecx
  jne    L154
  cmp    $0x10,%r8d
  jne    L154
  mov    $0x4,%r15d
  mov    $0x4,%r13d
  mov    $0x4,%r12d
  movl   $0x4,0xc(%rsp)
  jmp    L29
L154:
  cmp    $0x3,%r8d
  jne    L148
  test   %dl,%dl
  je     L148
  mov    $0x1,%r15d
  mov    $0x3,%r13d
  mov    $0x3,%r12d
  movl   $0x3,0xc(%rsp)
  jmp    L29
L148:
  cmp    $0x2,%r8d
  sete   %dl
  test   %al,%al
  je     L57
  test   %dl,%dl
  je     L57
  mov    $0x1,%r15d
  mov    $0x2,%r13d
  mov    $0x2,%r12d
  movl   $0x2,0xc(%rsp)
  jmp    L29
L57:
  cmp    $0x1,%ecx
  jne    L61
  cmp    $0x1,%r8d
  jne    L61
  mov    $0x1,%r15d
  mov    $0x1,%r13d
  mov    $0x1,%r12d
  movl   $0x1,0xc(%rsp)
  jmp    L29
L61:
  cmp    $0x6,%ecx
  jne    L39
  test   %dl,%dl
  je     L39
  movss  (%rbx),%xmm0
  mulss  0x0(%rbp),%xmm0
  movss  0x4(%rbx),%xmm1
  mulss  0x4(%rbp),%xmm1
  addss  %xmm1,%xmm0
  addss  0x8(%rbx),%xmm0
  movss  %xmm0,(%r14)
  movss  0xc(%rbx),%xmm0
  mulss  0x0(%rbp),%xmm0
  movss  0x10(%rbx),%xmm1
  mulss  0x4(%rbp),%xmm1
  addss  %xmm1,%xmm0
  addss  0x14(%rbx),%xmm0
  movss  %xmm0,0x4(%r14)
  mov    $0x2,%eax
  jmp    L56
L39:
  mov    %r8d,%edx
  mov    %ecx,%esi
  lea    D8(%rip),%rdi
  call   <fail@plt>
  jmp    L29
L85:
  mov    $0x2,%r15d
  mov    $0x2,%r13d
  mov    $0x2,%r12d
  movl   $0x2,0xc(%rsp)
  jmp    L29
L98:
  call   <__stack_chk_fail@plt>
```

tool_calls: name: retrieve_data
arguments: data_tag: D8
type: string

tool: D8: bad matrix product (%d %d)

assistant: ```c
int matrix_product(float *ab, float *a, float *b, int na, int nb) {
  int a_nrows, a_ncols, b_nrows, b_ncols;
  if (na == 4 && nb == 4) {
    a_nrows = a_ncols = b_nrows = b_ncols = 2;
  } else if (na == 9 && nb == 9) {
    a_nrows = a_ncols = b_nrows = b_ncols = 3;
  } else if (na == 16 && nb == 16) {
    a_nrows = a_ncols = b_nrows = b_ncols = 4;
  } else if (na == 9 && nb == 3) {
    a_nrows = a_ncols = b_nrows = 3;
    b_ncols = 1;
  } else if (na == 4 && nb == 2) {
    a_nrows = a_ncols = b_nrows = 2;
    b_ncols = 1;
  } else if (na == 1 && nb == 1) {
    a_nrows = a_ncols = b_nrows = b_ncols = 1;
  } else if (na == 6 && nb == 2) {
    ab[0] = a[0] * b[0] + a[1] * b[1] + a[2];
    ab[1] = a[3] * b[0] + a[4] * b[1] + a[5];
    return 2;
  } else
    fail("bad matrix product (%d %d)", na, nb);
  assert(a_ncols == b_nrows);
  int ab_nrows, ab_ncols;
  matrix_product_clean(ab, &ab_nrows, &ab_ncols, a, a_nrows, a_ncols, b, b_nrows, b_ncols);
  return ab_nrows * ab_ncols;
}
```