user: 请将这段汇编代码片段转换成等效的C代码：

```asm
<slevmar_trans_mat_mat_mult>:
  endbr64
  push   %rbp
  mov    %rsp,%rbp
  push   %r15
  push   %r14
  push   %r13
  push   %r12
  push   %rbx
  mov    %rdi,-0x40(%rbp)
  mov    %rsi,-0x48(%rbp)
  mov    %edx,-0x4c(%rbp)
  mov    %ecx,-0x50(%rbp)
  movl   $0x20,-0x2c(%rbp)
  mov    $0x0,%r13d
  jmp    L30
L10:
  mov    $0x0,%ebx
  jmp    L23
L164:
  mov    %ebx,%eax
  imul   -0x50(%rbp),%eax
  cltq
  lea    0x0(,%rax,4),%rdx
  mov    -0x48(%rbp),%rax
  add    %rdx,%rax
  mov    %rax,%rcx
  cmp    %ebx,%r13d
  mov    %ebx,%r12d
  cmovge %r13d,%r12d
  jmp    L153
L80:
  movslq %r12d,%rax
  shl    $0x2,%rax
  add    %rcx,%rax
  pxor   %xmm0,%xmm0
  movss  %xmm0,(%rax)
  add    $0x1,%r12d
L153:
  mov    -0x2c(%rbp),%eax
  lea    0x0(%r13,%rax,1),%edx
  mov    -0x50(%rbp),%eax
  cmp    %eax,%edx
  cmovle %edx,%eax
  cmp    %eax,%r12d
  jl     L80
  add    $0x1,%ebx
L23:
  cmp    -0x50(%rbp),%ebx
  jl     L164
  mov    $0x0,%r15d
  jmp    L72
L61:
  mov    $0x0,%ebx
  jmp    L106
L51:
  mov    %ebx,%eax
  imul   -0x50(%rbp),%eax
  cltq
  lea    0x0(,%rax,4),%rdx
  mov    -0x48(%rbp),%rax
  add    %rdx,%rax
  mov    %rax,%rcx
  cmp    %ebx,%r13d
  mov    %ebx,%r12d
  cmovge %r13d,%r12d
  jmp    L122
L152:
  pxor   %xmm2,%xmm2
  movaps %xmm2,%xmm5
  mov    %r15d,%r14d
  jmp    L101
L157:
  mov    %r14d,%eax
  imul   -0x50(%rbp),%eax
  cltq
  lea    0x0(,%rax,4),%rdx
  mov    -0x40(%rbp),%rax
  lea    (%rdx,%rax,1),%rsi
  movslq %ebx,%rax
  shl    $0x2,%rax
  add    %rsi,%rax
  movss  (%rax),%xmm1
  movslq %r12d,%rax
  shl    $0x2,%rax
  add    %rsi,%rax
  movss  (%rax),%xmm0
  mulss  %xmm1,%xmm0
  movaps %xmm5,%xmm3
  addss  %xmm3,%xmm0
  movaps %xmm0,%xmm5
  add    $0x1,%r14d
L101:
  mov    -0x2c(%rbp),%eax
  lea    (%r15,%rax,1),%edx
  mov    -0x4c(%rbp),%eax
  cmp    %eax,%edx
  cmovle %edx,%eax
  cmp    %eax,%r14d
  jl     L157
  movslq %r12d,%rax
  shl    $0x2,%rax
  mov    %rcx,%rdi
  add    %rdi,%rax
  movss  (%rax),%xmm0
  movslq %r12d,%rax
  shl    $0x2,%rax
  add    %rdi,%rax
  movaps %xmm5,%xmm4
  addss  %xmm4,%xmm0
  movss  %xmm0,(%rax)
  add    $0x1,%r12d
L122:
  mov    -0x2c(%rbp),%eax
  lea    0x0(%r13,%rax,1),%edx
  mov    -0x50(%rbp),%eax
  cmp    %eax,%edx
  cmovle %edx,%eax
  cmp    %eax,%r12d
  jl     L152
  add    $0x1,%ebx
L106:
  cmp    -0x50(%rbp),%ebx
  jl     L51
  mov    -0x2c(%rbp),%eax
  add    %r15d,%eax
  mov    %eax,%r15d
L72:
  cmp    -0x4c(%rbp),%r15d
  jl     L61
  mov    -0x2c(%rbp),%eax
  add    %r13d,%eax
  mov    %eax,%r13d
L30:
  cmp    -0x50(%rbp),%r13d
  jl     L10
  mov    $0x0,%ebx
  jmp    L177
L108:
  mov    $0x0,%r12d
  jmp    L171
L43:
  mov    %r12d,%eax
  imul   -0x50(%rbp),%eax
  add    %ebx,%eax
  cltq
  lea    0x0(,%rax,4),%rdx
  mov    -0x48(%rbp),%rax
  add    %rax,%rdx
  mov    %ebx,%eax
  imul   -0x50(%rbp),%eax
  add    %r12d,%eax
  cltq
  lea    0x0(,%rax,4),%rcx
  mov    -0x48(%rbp),%rax
  add    %rcx,%rax
  movss  (%rdx),%xmm0
  movss  %xmm0,(%rax)
  add    $0x1,%r12d
L171:
  cmp    %ebx,%r12d
  jl     L43
  add    $0x1,%ebx
L177:
  cmp    -0x50(%rbp),%ebx
  jl     L108
  nop
  nop
  pop    %rbx
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  pop    %rbp
  ret
```

assistant: ```c
void slevmar_trans_mat_mat_mult(float *a, float *b, int n, int m) {
  register int i, j, k, jj, kk;
  register float sum, *bim, *akm;
  const int bsize = 32;
  for (jj = 0; jj < m; jj += bsize) {
    for (i = 0; i < m; ++i) {
      bim = b + i * m;
      for (j = (((jj) >= (i)) ? (jj) : (i)); j < (((jj + bsize) <= (m)) ? (jj + bsize) : (m)); ++j)
        bim[j] = 0.0;
    }
    for (kk = 0; kk < n; kk += bsize) {
      for (i = 0; i < m; ++i) {
        bim = b + i * m;
        for (j = (((jj) >= (i)) ? (jj) : (i)); j < (((jj + bsize) <= (m)) ? (jj + bsize) : (m)); ++j) {
          sum = 0.0;
          for (k = kk; k < (((kk + bsize) <= (n)) ? (kk + bsize) : (n)); ++k) {
            akm = a + k * m;
            sum += akm[i] * akm[j];
          }
          bim[j] += sum;
        }
      }
    }
  }
  for (i = 0; i < m; ++i)
    for (j = 0; j < i; ++j)
      b[i * m + j] = b[j * m + i];
}
```