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

```asm
<matrix_product>:
  endbr64
  push   %rbp
  mov    %rsp,%rbp
  sub    $0x40,%rsp
  mov    %rdi,-0x28(%rbp)
  mov    %rsi,-0x30(%rbp)
  mov    %rdx,-0x38(%rbp)
  mov    %ecx,-0x3c(%rbp)
  mov    %r8d,-0x40(%rbp)
  mov    %fs:0x28,%rax
  mov    %rax,-0x8(%rbp)
  xor    %eax,%eax
  cmpl   $0x4,-0x3c(%rbp)
  jne    L2
  cmpl   $0x4,-0x40(%rbp)
  jne    L2
  movl   $0x2,-0xc(%rbp)
  mov    -0xc(%rbp),%eax
  mov    %eax,-0x10(%rbp)
  mov    -0x10(%rbp),%eax
  mov    %eax,-0x14(%rbp)
  mov    -0x14(%rbp),%eax
  mov    %eax,-0x18(%rbp)
  jmp    L123
L2:
  cmpl   $0x9,-0x3c(%rbp)
  jne    L1
  cmpl   $0x9,-0x40(%rbp)
  jne    L1
  movl   $0x3,-0xc(%rbp)
  mov    -0xc(%rbp),%eax
  mov    %eax,-0x10(%rbp)
  mov    -0x10(%rbp),%eax
  mov    %eax,-0x14(%rbp)
  mov    -0x14(%rbp),%eax
  mov    %eax,-0x18(%rbp)
  jmp    L123
L1:
  cmpl   $0x10,-0x3c(%rbp)
  jne    L127
  cmpl   $0x10,-0x40(%rbp)
  jne    L127
  movl   $0x4,-0xc(%rbp)
  mov    -0xc(%rbp),%eax
  mov    %eax,-0x10(%rbp)
  mov    -0x10(%rbp),%eax
  mov    %eax,-0x14(%rbp)
  mov    -0x14(%rbp),%eax
  mov    %eax,-0x18(%rbp)
  jmp    L123
L127:
  cmpl   $0x9,-0x3c(%rbp)
  jne    L5
  cmpl   $0x3,-0x40(%rbp)
  jne    L5
  movl   $0x3,-0x10(%rbp)
  mov    -0x10(%rbp),%eax
  mov    %eax,-0x14(%rbp)
  mov    -0x14(%rbp),%eax
  mov    %eax,-0x18(%rbp)
  movl   $0x1,-0xc(%rbp)
  jmp    L123
L5:
  cmpl   $0x4,-0x3c(%rbp)
  jne    L94
  cmpl   $0x2,-0x40(%rbp)
  jne    L94
  movl   $0x2,-0x10(%rbp)
  mov    -0x10(%rbp),%eax
  mov    %eax,-0x14(%rbp)
  mov    -0x14(%rbp),%eax
  mov    %eax,-0x18(%rbp)
  movl   $0x1,-0xc(%rbp)
  jmp    L123
L94:
  cmpl   $0x1,-0x3c(%rbp)
  jne    L157
  cmpl   $0x1,-0x40(%rbp)
  jne    L157
  movl   $0x1,-0xc(%rbp)
  mov    -0xc(%rbp),%eax
  mov    %eax,-0x10(%rbp)
  mov    -0x10(%rbp),%eax
  mov    %eax,-0x14(%rbp)
  mov    -0x14(%rbp),%eax
  mov    %eax,-0x18(%rbp)
  jmp    L123
L157:
  cmpl   $0x6,-0x3c(%rbp)
  jne    L199
  cmpl   $0x2,-0x40(%rbp)
  jne    L199
  mov    -0x30(%rbp),%rax
  movss  (%rax),%xmm1
  mov    -0x38(%rbp),%rax
  movss  (%rax),%xmm0
  mulss  %xmm0,%xmm1
  mov    -0x30(%rbp),%rax
  add    $0x4,%rax
  movss  (%rax),%xmm2
  mov    -0x38(%rbp),%rax
  add    $0x4,%rax
  movss  (%rax),%xmm0
  mulss  %xmm2,%xmm0
  addss  %xmm0,%xmm1
  mov    -0x30(%rbp),%rax
  add    $0x8,%rax
  movss  (%rax),%xmm0
  addss  %xmm1,%xmm0
  mov    -0x28(%rbp),%rax
  movss  %xmm0,(%rax)
  mov    -0x30(%rbp),%rax
  add    $0xc,%rax
  movss  (%rax),%xmm1
  mov    -0x38(%rbp),%rax
  movss  (%rax),%xmm0
  mulss  %xmm0,%xmm1
  mov    -0x30(%rbp),%rax
  add    $0x10,%rax
  movss  (%rax),%xmm2
  mov    -0x38(%rbp),%rax
  add    $0x4,%rax
  movss  (%rax),%xmm0
  mulss  %xmm2,%xmm0
  addss  %xmm0,%xmm1
  mov    -0x30(%rbp),%rax
  add    $0x14,%rax
  movss  (%rax),%xmm0
  mov    -0x28(%rbp),%rax
  add    $0x4,%rax
  addss  %xmm1,%xmm0
  movss  %xmm0,(%rax)
  mov    $0x2,%eax
  jmp    L125
L199:
  mov    -0x40(%rbp),%edx
  mov    -0x3c(%rbp),%eax
  mov    %eax,%esi
  lea    D8(%rip),%rax
  mov    %rax,%rdi
  call   <fail@plt>
L123:
  mov    -0x14(%rbp),%eax
  cmp    -0x10(%rbp),%eax
  sete   %al
  movzbl %al,%eax
  mov    %eax,%edi
  call   <assert@plt>
  mov    -0x14(%rbp),%r9d
  mov    -0x18(%rbp),%r8d
  mov    -0x30(%rbp),%rcx
  lea    -0x1c(%rbp),%rdx
  lea    -0x20(%rbp),%rsi
  mov    -0x28(%rbp),%rax
  sub    $0x8,%rsp
  mov    -0xc(%rbp),%edi
  push   %rdi
  mov    -0x10(%rbp),%edi
  push   %rdi
  push   -0x38(%rbp)
  mov    %rax,%rdi
  call   <matrix_product_clean@plt>
  add    $0x20,%rsp
  mov    -0x20(%rbp),%edx
  mov    -0x1c(%rbp),%eax
  imul   %edx,%eax
L125:
  mov    -0x8(%rbp),%rdx
  sub    %fs:0x28,%rdx
  je     L47
  call   <__stack_chk_fail@plt>
L47:
  leave
  ret
```

tool_calls: name: extract_symbolic_data
arguments: data_label: D8
type_infered: 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;
}
```