user: Translate the below assembly code into C programming language:

```asm
<mulmp4>:
  endbr64
  push   %rbp
  mov    %rsp,%rbp
  mov    %rdi,-0x18(%rbp)
  mov    %rsi,-0x20(%rbp)
  mov    %rdx,-0x28(%rbp)
  mov    -0x20(%rbp),%rax
  movss  (%rax),%xmm1
  mov    -0x28(%rbp),%rax
  movss  (%rax),%xmm0
  mulss  %xmm0,%xmm1
  mov    -0x20(%rbp),%rax
  add    $0x10,%rax
  movss  (%rax),%xmm2
  mov    -0x28(%rbp),%rax
  add    $0x4,%rax
  movss  (%rax),%xmm0
  mulss  %xmm2,%xmm0
  addss  %xmm0,%xmm1
  mov    -0x20(%rbp),%rax
  add    $0x20,%rax
  movss  (%rax),%xmm2
  mov    -0x28(%rbp),%rax
  add    $0x8,%rax
  movss  (%rax),%xmm0
  mulss  %xmm2,%xmm0
  addss  %xmm0,%xmm1
  mov    -0x20(%rbp),%rax
  add    $0x30,%rax
  movss  (%rax),%xmm2
  mov    -0x28(%rbp),%rax
  add    $0xc,%rax
  movss  (%rax),%xmm0
  mulss  %xmm2,%xmm0
  addss  %xmm1,%xmm0
  movss  %xmm0,-0x10(%rbp)
  mov    -0x20(%rbp),%rax
  add    $0x4,%rax
  movss  (%rax),%xmm1
  mov    -0x28(%rbp),%rax
  movss  (%rax),%xmm0
  mulss  %xmm0,%xmm1
  mov    -0x20(%rbp),%rax
  add    $0x14,%rax
  movss  (%rax),%xmm2
  mov    -0x28(%rbp),%rax
  add    $0x4,%rax
  movss  (%rax),%xmm0
  mulss  %xmm2,%xmm0
  addss  %xmm0,%xmm1
  mov    -0x20(%rbp),%rax
  add    $0x24,%rax
  movss  (%rax),%xmm2
  mov    -0x28(%rbp),%rax
  add    $0x8,%rax
  movss  (%rax),%xmm0
  mulss  %xmm2,%xmm0
  addss  %xmm0,%xmm1
  mov    -0x20(%rbp),%rax
  add    $0x34,%rax
  movss  (%rax),%xmm2
  mov    -0x28(%rbp),%rax
  add    $0xc,%rax
  movss  (%rax),%xmm0
  mulss  %xmm2,%xmm0
  addss  %xmm1,%xmm0
  movss  %xmm0,-0xc(%rbp)
  mov    -0x20(%rbp),%rax
  add    $0x8,%rax
  movss  (%rax),%xmm1
  mov    -0x28(%rbp),%rax
  movss  (%rax),%xmm0
  mulss  %xmm0,%xmm1
  mov    -0x20(%rbp),%rax
  add    $0x18,%rax
  movss  (%rax),%xmm2
  mov    -0x28(%rbp),%rax
  add    $0x4,%rax
  movss  (%rax),%xmm0
  mulss  %xmm2,%xmm0
  addss  %xmm0,%xmm1
  mov    -0x20(%rbp),%rax
  add    $0x28,%rax
  movss  (%rax),%xmm2
  mov    -0x28(%rbp),%rax
  add    $0x8,%rax
  movss  (%rax),%xmm0
  mulss  %xmm2,%xmm0
  addss  %xmm0,%xmm1
  mov    -0x20(%rbp),%rax
  add    $0x38,%rax
  movss  (%rax),%xmm2
  mov    -0x28(%rbp),%rax
  add    $0xc,%rax
  movss  (%rax),%xmm0
  mulss  %xmm2,%xmm0
  addss  %xmm1,%xmm0
  movss  %xmm0,-0x8(%rbp)
  mov    -0x20(%rbp),%rax
  add    $0xc,%rax
  movss  (%rax),%xmm1
  mov    -0x28(%rbp),%rax
  movss  (%rax),%xmm0
  mulss  %xmm0,%xmm1
  mov    -0x20(%rbp),%rax
  add    $0x1c,%rax
  movss  (%rax),%xmm2
  mov    -0x28(%rbp),%rax
  add    $0x4,%rax
  movss  (%rax),%xmm0
  mulss  %xmm2,%xmm0
  addss  %xmm0,%xmm1
  mov    -0x20(%rbp),%rax
  add    $0x2c,%rax
  movss  (%rax),%xmm2
  mov    -0x28(%rbp),%rax
  add    $0x8,%rax
  movss  (%rax),%xmm0
  mulss  %xmm2,%xmm0
  addss  %xmm0,%xmm1
  mov    -0x20(%rbp),%rax
  add    $0x3c,%rax
  movss  (%rax),%xmm2
  mov    -0x28(%rbp),%rax
  add    $0xc,%rax
  movss  (%rax),%xmm0
  mulss  %xmm2,%xmm0
  addss  %xmm1,%xmm0
  movss  %xmm0,-0x4(%rbp)
  mov    -0x18(%rbp),%rax
  movss  -0x10(%rbp),%xmm0
  movss  %xmm0,(%rax)
  mov    -0x18(%rbp),%rax
  add    $0x4,%rax
  movss  -0xc(%rbp),%xmm0
  movss  %xmm0,(%rax)
  mov    -0x18(%rbp),%rax
  add    $0x8,%rax
  movss  -0x8(%rbp),%xmm0
  movss  %xmm0,(%rax)
  mov    -0x18(%rbp),%rax
  add    $0xc,%rax
  movss  -0x4(%rbp),%xmm0
  movss  %xmm0,(%rax)
  nop
  pop    %rbp
  ret
```

assistant: ```c
void mulmp4(float *result, const float *s1, const float *s2) {
  const float x = s1[(0 + (0 << 2))] * s2[0] + s1[(0 + (1 << 2))] * s2[1] + s1[(0 + (2 << 2))] * s2[2] + s1[(0 + (3 << 2))] * s2[3];
  const float y = s1[(1 + (0 << 2))] * s2[0] + s1[(1 + (1 << 2))] * s2[1] + s1[(1 + (2 << 2))] * s2[2] + s1[(1 + (3 << 2))] * s2[3];
  const float z = s1[(2 + (0 << 2))] * s2[0] + s1[(2 + (1 << 2))] * s2[1] + s1[(2 + (2 << 2))] * s2[2] + s1[(2 + (3 << 2))] * s2[3];
  const float w = s1[(3 + (0 << 2))] * s2[0] + s1[(3 + (1 << 2))] * s2[1] + s1[(3 + (2 << 2))] * s2[2] + s1[(3 + (3 << 2))] * s2[3];
  result[0] = x;
  result[1] = y;
  result[2] = z;
  result[3] = w;
}
```