user: What would the C source look like for the assembly code provided here?

```asm
<matmatkij>:
  endbr64
  push   %r15
  push   %r14
  push   %r13
  push   %r12
  push   %rbp
  push   %rbx
  movslq %edi,%rbx
  mov    %rdx,-0x30(%rsp)
  mov    %rcx,-0x48(%rsp)
  mov    %rbx,-0x40(%rsp)
  test   %r8d,%r8d
  jle    L61
  mov    %rsi,%rax
  mov    %r8d,%ebp
  mov    %r8d,%r13d
  mov    %rdx,%rsi
  mov    %rax,-0x10(%rsp)
  lea    0x4(%rdx),%rax
  mov    %ebp,%r14d
  lea    0x0(,%rbx,4),%r12
  mov    %rax,-0x38(%rsp)
  lea    -0x1(%r8),%eax
  and    $0xfffffffc,%r14d
  xor    %r15d,%r15d
  mov    %eax,-0x20(%rsp)
  movslq %r8d,%rax
  movq   $0x0,-0x18(%rsp)
  shl    $0x2,%rax
  mov    %rax,-0x28(%rsp)
  mov    %r8d,%eax
  shr    $0x2,%eax
  lea    -0x1(%rax),%r8d
  mov    %ebp,%eax
  and    $0x3,%eax
  add    $0x1,%r8
  mov    %eax,-0x1c(%rsp)
  shl    $0x4,%r8
  nopl   (%rax)
L34:
  mov    -0x48(%rsp),%rax
  mov    -0x10(%rsp),%rcx
  xor    %r10d,%r10d
  xor    %r11d,%r11d
L132:
  mov    -0x28(%rsp),%rbx
  lea    0x4(%rcx),%rdi
  lea    (%rbx,%rax,1),%rdx
  cmp    %rdx,%rcx
  setae  %dl
  cmp    %rdi,%rax
  setae  %dil
  or     %edi,%edx
  mov    %rax,%rdi
  sub    -0x38(%rsp),%rdi
  cmp    $0x8,%rdi
  seta   %dil
  test   %dil,%dl
  je     L83
  cmp    $0x1,%ebp
  je     L83
  cmpl   $0x2,-0x20(%rsp)
  jbe    L70
  movss  (%rcx),%xmm1
  xor    %edx,%edx
  shufps $0x0,%xmm1,%xmm1
  nopl   0x0(%rax,%rax,1)
L44:
  movups (%rsi,%rdx,1),%xmm0
  movups (%rax,%rdx,1),%xmm2
  mulps  %xmm1,%xmm0
  addps  %xmm2,%xmm0
  movups %xmm0,(%rax,%rdx,1)
  add    $0x10,%rdx
  cmp    %r8,%rdx
  jne    L44
  cmp    %ebp,%r14d
  je     L115
  mov    -0x1c(%rsp),%edi
  mov    %edi,%ebx
  cmp    $0x1,%edi
  je     L116
  mov    %r14d,-0x4c(%rsp)
  mov    %r14d,%edi
L11:
  mov    -0x48(%rsp),%rdx
  lea    (%rdi,%r10,1),%r9
  add    %r15,%rdi
  movss  (%rcx),%xmm1
  lea    (%rdx,%r9,4),%r9
  mov    -0x30(%rsp),%rdx
  shufps $0xe0,%xmm1,%xmm1
  movq   (%rdx,%rdi,4),%xmm0
  mov    -0x4c(%rsp),%edx
  mov    %ebx,%edi
  and    $0xfffffffe,%edi
  mulps  %xmm1,%xmm0
  movq   (%r9),%xmm1
  add    %edi,%edx
  addps  %xmm0,%xmm1
  movlps %xmm1,(%r9)
  cmp    %ebx,%edi
  je     L115
L50:
  movslq %edx,%rdx
  movss  (%rcx),%xmm0
  mulss  (%rsi,%rdx,4),%xmm0
  addss  (%rax,%rdx,4),%xmm0
  movss  %xmm0,(%rax,%rdx,4)
L115:
  add    $0x1,%r11d
  add    %r12,%rax
  add    %r12,%rcx
  add    -0x40(%rsp),%r10
  cmp    %r11d,%ebp
  jne    L132
  addq   $0x1,-0x18(%rsp)
  mov    -0x18(%rsp),%rax
  add    %r12,%rsi
  addq   $0x4,-0x10(%rsp)
  add    -0x40(%rsp),%r15
  add    %r12,-0x38(%rsp)
  cmp    %r13,%rax
  jne    L34
L61:
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
  nop
L83:
  xor    %edx,%edx
  nopw   0x0(%rax,%rax,1)
L110:
  movss  (%rsi,%rdx,4),%xmm0
  mulss  (%rcx),%xmm0
  addss  (%rax,%rdx,4),%xmm0
  movss  %xmm0,(%rax,%rdx,4)
  add    $0x1,%rdx
  cmp    %r13,%rdx
  jne    L110
  jmp    L115
L70:
  movl   $0x0,-0x4c(%rsp)
  mov    %ebp,%ebx
  xor    %edi,%edi
  jmp    L11
L116:
  mov    %r14d,%edx
  jmp    L50
```

assistant: ```c
void matmatkij(int LD, float A[][LD], float B[][LD], float C[][LD], int N) {
  int i = 0, j = 0, k = 0;
  for (k = 0; k < N; k++)
    for (i = 0; i < N; i++)
      for (j = 0; j < N; j++)
        C[i][j] += A[i][k] * B[k][j];
}
```