user: 请显示与下面列出的汇编代码相匹配的C语言实现:

```asm
<srecmult>:
  endbr64
  push   %r15
  push   %r14
  push   %r13
  push   %r12
  push   %rbp
  push   %rbx
  sub    $0x98,%rsp
  mov    %rdi,(%rsp)
  mov    %rsi,0x8(%rsp)
  mov    %rdx,0x10(%rsp)
  mov    %ecx,%r13d
  mov    %r9d,%r15d
  mov    0xd8(%rsp),%r14d
  lea    (%rcx,%r8,1),%eax
  add    %r9d,%eax
  cmp    $0xc0,%eax
  jg     L45
  test   %ecx,%ecx
  jle    L92
  movslq %r14d,%rsi
  shl    $0x3,%rsi
  mov    $0x0,%r14d
  mov    $0x0,%r12d
  mov    $0x0,%ebp
  mov    %r9d,%r9d
  lea    -0x1(%r8),%eax
  mov    %rax,0x18(%rsp)
  lea    0x8(%rdi),%rax
  mov    %rax,0x20(%rsp)
  jmp    L28
L75:
  mov    0x8(%rsp),%r10
  movslq %r14d,%rax
  mov    (%rsp),%rbx
  lea    (%rbx,%rax,8),%rbx
  add    0x18(%rsp),%rax
  mov    0x20(%rsp),%rdi
  lea    (%rdi,%rax,8),%rcx
  movslq %r12d,%rax
  mov    0x10(%rsp),%rdi
  lea    (%rdi,%rax,8),%r11
  mov    $0x0,%edi
L55:
  test   %r8d,%r8d
  jle    L48
  mov    %r10,%rdx
  mov    %rbx,%rax
  pxor   %xmm1,%xmm1
L85:
  movsd  (%rax),%xmm0
  mulsd  (%rdx),%xmm0
  addsd  %xmm0,%xmm1
  add    $0x8,%rax
  add    %rsi,%rdx
  cmp    %rcx,%rax
  jne    L85
L33:
  movsd  %xmm1,(%r11,%rdi,8)
  add    $0x1,%rdi
  add    $0x8,%r10
  cmp    %r9,%rdi
  jne    L55
L84:
  add    $0x1,%ebp
  add    0xe0(%rsp),%r12d
  add    0xd0(%rsp),%r14d
  cmp    %ebp,%r13d
  je     L92
L28:
  test   %r15d,%r15d
  jg     L75
  jmp    L84
L48:
  pxor   %xmm1,%xmm1
  jmp    L33
L45:
  mov    %ecx,%ebx
  shr    $0x1f,%ebx
  add    %ecx,%ebx
  sar    %ebx
  mov    %r8d,%ebp
  shr    $0x1f,%ebp
  add    %r8d,%ebp
  sar    %ebp
  mov    %r9d,%r12d
  shr    $0x1f,%r12d
  add    %r9d,%r12d
  sar    %r12d
  movslq %ebp,%rax
  mov    (%rsp),%rdi
  lea    (%rdi,%rax,8),%rcx
  mov    %rcx,0x58(%rsp)
  mov    %ebx,%edx
  imul   0xd0(%rsp),%edx
  movslq %edx,%rdx
  lea    (%rdi,%rdx,8),%rsi
  mov    %rsi,0x60(%rsp)
  add    %rdx,%rax
  lea    (%rdi,%rax,8),%rax
  mov    %rax,0x28(%rsp)
  movslq %r12d,%rdx
  lea    0x0(,%rdx,8),%rcx
  mov    0x8(%rsp),%rsi
  lea    (%rsi,%rcx,1),%r15
  mov    %r15,0x68(%rsp)
  mov    %ebp,%eax
  imul   %r14d,%eax
  cltq
  lea    (%rsi,%rax,8),%r11
  mov    %r11,0x70(%rsp)
  add    %rdx,%rax
  lea    (%rsi,%rax,8),%r15
  mov    %r15,0x30(%rsp)
  mov    0x10(%rsp),%rdi
  add    %rdi,%rcx
  mov    %rcx,0x78(%rsp)
  mov    %ebx,%eax
  imul   0xe0(%rsp),%eax
  cltq
  lea    (%rdi,%rax,8),%r11
  mov    %r11,0x80(%rsp)
  add    %rax,%rdx
  lea    (%rdi,%rdx,8),%r15
  mov    %r15,0x88(%rsp)
  lea    0x0(,%r12,8),%r13d
  mov    %ebx,%edx
  imul   %r13d,%edx
  mov    %edx,%r15d
  mov    %edx,%edi
  call   <malloc@plt>
  mov    %rax,0x18(%rsp)
  mov    %r15d,%edi
  call   <malloc@plt>
  mov    %rax,0x38(%rsp)
  mov    %r15d,%edi
  call   <malloc@plt>
  mov    %rax,0x40(%rsp)
  mov    %ebx,%edx
  imul   %ebp,%edx
  lea    0x0(,%rdx,8),%edx
  mov    %edx,0x20(%rsp)
  mov    %edx,%edi
  call   <malloc@plt>
  mov    %rax,0x48(%rsp)
  imul   %ebp,%r13d
  mov    %r13d,%edi
  call   <malloc@plt>
  mov    %rax,0x50(%rsp)
  push   %rbp
  mov    0xd8(%rsp),%r9d
  push   %r9
  mov    0xe0(%rsp),%r9d
  mov    %ebp,%r8d
  mov    %ebx,%ecx
  mov    0x58(%rsp),%rdx
  mov    0x38(%rsp),%rsi
  mov    0x10(%rsp),%rdi
  call   <madd_stride@plt>
  push   %r12
  push   %r14
  mov    %r14d,%r9d
  mov    %r12d,%r8d
  mov    %ebp,%ecx
  mov    0x70(%rsp),%rdx
  mov    0x50(%rsp),%rsi
  mov    0x28(%rsp),%rdi
  call   <madd_stride@plt>
  add    $0x18,%rsp
  push   %r12
  push   %r12
  push   %rbp
  mov    %r12d,%r9d
  mov    %ebp,%r8d
  mov    %ebx,%ecx
  mov    0x38(%rsp),%rdx
  mov    0x70(%rsp),%rsi
  mov    0x68(%rsp),%rdi
  call   <srecmult>
  add    $0x20,%rsp
  mov    0x48(%rsp),%rdi
  call   <free@plt>
  mov    0x50(%rsp),%rdi
  call   <free@plt>
  mov    0x20(%rsp),%edi
  call   <malloc@plt>
  push   %rbp
  mov    0xd8(%rsp),%r9d
  push   %r9
  mov    0xe0(%rsp),%r9d
  mov    %ebp,%r8d
  mov    %ebx,%ecx
  mov    %rax,0x58(%rsp)
  mov    %rax,%rdx
  mov    0x38(%rsp),%rsi
  mov    0x70(%rsp),%rdi
  call   <madd_stride@plt>
  mov    0xf0(%rsp),%r9d
  mov    %r9d,(%rsp)
  push   %r14
  push   %rbp
  mov    %r12d,%r9d
  mov    %ebp,%r8d
  mov    %ebx,%ecx
  mov    0xa0(%rsp),%rdx
  mov    0x28(%rsp),%rsi
  mov    0x68(%rsp),%rdi
  call   <srecmult>
  add    $0x20,%rsp
  mov    0x48(%rsp),%rdi
  call   <free@plt>
  mov    %r13d,%edi
  call   <malloc@plt>
  push   %r12
  push   %r14
  mov    %r14d,%r9d
  mov    %r12d,%r8d
  mov    %ebp,%ecx
  mov    %rax,0x58(%rsp)
  mov    %rax,%rdx
  mov    0x40(%rsp),%rsi
  mov    0x78(%rsp),%rdi
  call   <msub_stride@plt>
  mov    %r12d,(%rsp)
  push   %r12
  mov    0xe8(%rsp),%r8d
  push   %r8
  mov    %r12d,%r9d
  mov    %ebp,%r8d
  mov    %ebx,%ecx
  mov    0x58(%rsp),%rdx
  mov    0x68(%rsp),%rsi
  mov    0x20(%rsp),%rdi
  call   <srecmult>
  add    $0x20,%rsp
  mov    0x48(%rsp),%rdi
  call   <free@plt>
  mov    %r13d,%edi
  call   <malloc@plt>
  push   %r12
  push   %r14
  mov    %r14d,%r9d
  mov    %r12d,%r8d
  mov    %ebp,%ecx
  mov    %rax,0x58(%rsp)
  mov    %rax,%rdx
  mov    0x18(%rsp),%rsi
  mov    0x80(%rsp),%rdi
  call   <msub_stride@plt>
  mov    %r12d,(%rsp)
  push   %r12
  mov    0xe8(%rsp),%r9d
  push   %r9
  mov    %r12d,%r9d
  mov    %ebp,%r8d
  mov    %ebx,%ecx
  mov    0x60(%rsp),%rdx
  mov    0x68(%rsp),%rsi
  mov    0x48(%rsp),%rdi
  call   <srecmult>
  add    $0x20,%rsp
  mov    0x48(%rsp),%rdi
  call   <free@plt>
  mov    0x20(%rsp),%edi
  call   <malloc@plt>
  push   %rbp
  mov    0xd8(%rsp),%r8d
  push   %r8
  mov    0xe0(%rsp),%r9d
  mov    %ebp,%r8d
  mov    %ebx,%ecx
  mov    %rax,0x58(%rsp)
  mov    %rax,%rdx
  mov    0x68(%rsp),%rsi
  mov    0x10(%rsp),%rdi
  call   <madd_stride@plt>
  mov    0xf0(%rsp),%r9d
  mov    %r9d,(%rsp)
  push   %r14
  push   %rbp
  mov    %r12d,%r9d
  mov    %ebp,%r8d
  mov    %ebx,%ecx
  mov    0x98(%rsp),%rdx
  mov    0x50(%rsp),%rsi
  mov    0x68(%rsp),%rdi
  call   <srecmult>
  add    $0x20,%rsp
  mov    0x48(%rsp),%rdi
  call   <free@plt>
  mov    0x20(%rsp),%edi
  call   <malloc@plt>
  mov    %rax,0x48(%rsp)
  mov    %r13d,%edi
  call   <malloc@plt>
  mov    %rax,0x50(%rsp)
  push   %rbp
  mov    0xd8(%rsp),%ecx
  push   %rcx
  mov    0xe0(%rsp),%r9d
  mov    %ebp,%r8d
  mov    %ebx,%ecx
  mov    0x58(%rsp),%rdx
  mov    0x10(%rsp),%rsi
  mov    0x70(%rsp),%rdi
  call   <msub_stride@plt>
  push   %r12
  push   %r14
  mov    %r14d,%r9d
  mov    %r12d,%r8d
  mov    %ebp,%ecx
  mov    0x70(%rsp),%rdx
  mov    0x88(%rsp),%rsi
  mov    0x28(%rsp),%rdi
  call   <madd_stride@plt>
  add    $0x18,%rsp
  mov    0xe8(%rsp),%eax
  push   %rax
  push   %r12
  push   %rbp
  mov    %r12d,%r9d
  mov    %ebp,%r8d
  mov    %ebx,%ecx
  mov    0xa8(%rsp),%rdx
  mov    0x70(%rsp),%rsi
  mov    0x68(%rsp),%rdi
  call   <srecmult>
  add    $0x20,%rsp
  mov    0x48(%rsp),%rdi
  call   <free@plt>
  mov    0x50(%rsp),%rdi
  call   <free@plt>
  mov    0x20(%rsp),%edi
  call   <malloc@plt>
  mov    %rax,(%rsp)
  mov    %r13d,%edi
  call   <malloc@plt>
  mov    %rax,%r13
  push   %rbp
  mov    0xd8(%rsp),%ecx
  push   %rcx
  mov    0xe0(%rsp),%r9d
  mov    %ebp,%r8d
  mov    %ebx,%ecx
  mov    0x10(%rsp),%rdx
  mov    0x38(%rsp),%rsi
  mov    0x68(%rsp),%rdi
  call   <msub_stride@plt>
  push   %r12
  push   %r14
  mov    %r14d,%r9d
  mov    %r12d,%r8d
  mov    %ebp,%ecx
  mov    %r13,%rdx
  mov    0x50(%rsp),%rsi
  mov    0x90(%rsp),%rdi
  call   <madd_stride@plt>
  add    $0x18,%rsp
  mov    0xe8(%rsp),%eax
  push   %rax
  push   %r12
  push   %rbp
  mov    %r12d,%r9d
  mov    %ebp,%r8d
  mov    %ebx,%ecx
  mov    0x30(%rsp),%r14
  mov    %r14,%rdx
  mov    %r13,%rsi
  mov    0x20(%rsp),%rbp
  mov    %rbp,%rdi
  call   <srecmult>
  add    $0x20,%rsp
  mov    %rbp,%rdi
  call   <free@plt>
  mov    %r13,%rdi
  call   <free@plt>
  mov    %r15d,%edi
  call   <malloc@plt>
  mov    %rax,%rbp
  mov    0xe0(%rsp),%eax
  push   %rax
  mov    0xe8(%rsp),%eax
  push   %rax
  mov    %r12d,%r9d
  mov    %r12d,%r8d
  mov    %ebx,%ecx
  mov    %r14,%rdx
  mov    %r14,%rsi
  mov    0x28(%rsp),%rdi
  call   <madd_stride@plt>
  push   %r12
  mov    0xf8(%rsp),%eax
  push   %rax
  mov    %r12d,%r9d
  mov    %r12d,%r8d
  mov    %ebx,%ecx
  mov    %rbp,%rdx
  mov    0x98(%rsp),%r13
  mov    %r13,%rsi
  mov    0x60(%rsp),%rdi
  call   <msub_stride@plt>
  add    $0x20,%rsp
  mov    0xe0(%rsp),%eax
  push   %rax
  push   %r12
  mov    0xf0(%rsp),%r9d
  mov    %r12d,%r8d
  mov    %ebx,%ecx
  mov    %r14,%rdx
  mov    %rbp,%rsi
  mov    %r14,%rdi
  call   <madd_stride@plt>
  mov    %rbp,%rdi
  call   <free@plt>
  mov    %r15d,%edi
  call   <malloc@plt>
  mov    %rax,%rbp
  mov    0xf0(%rsp),%eax
  push   %rax
  mov    0xf8(%rsp),%eax
  push   %rax
  mov    %r12d,%r9d
  mov    %r12d,%r8d
  mov    %ebx,%ecx
  mov    0xa8(%rsp),%r15
  mov    %r15,%rdx
  mov    %r15,%rsi
  mov    0x38(%rsp),%rdi
  call   <madd_stride@plt>
  add    $0x20,%rsp
  push   %r12
  mov    0xe8(%rsp),%eax
  push   %rax
  mov    %r12d,%r9d
  mov    %r12d,%r8d
  mov    %ebx,%ecx
  mov    %rbp,%rdx
  mov    0x90(%rsp),%r14
  mov    %r14,%rsi
  mov    0x48(%rsp),%rdi
  call   <msub_stride@plt>
  mov    0xf0(%rsp),%eax
  push   %rax
  push   %r12
  mov    0x100(%rsp),%r9d
  mov    %r12d,%r8d
  mov    %ebx,%ecx
  mov    %r15,%rdx
  mov    %rbp,%rsi
  mov    %r15,%rdi
  call   <madd_stride@plt>
  add    $0x20,%rsp
  mov    %rbp,%rdi
  call   <free@plt>
  mov    0xe0(%rsp),%eax
  push   %rax
  mov    0xe8(%rsp),%eax
  push   %rax
  mov    %r12d,%r9d
  mov    %r12d,%r8d
  mov    %ebx,%ecx
  mov    %r13,%rdx
  mov    %r13,%rsi
  mov    0x48(%rsp),%r15
  mov    %r15,%rdi
  call   <madd_stride@plt>
  mov    0xf0(%rsp),%eax
  push   %rax
  mov    0xf8(%rsp),%eax
  push   %rax
  mov    %r12d,%r9d
  mov    %r12d,%r8d
  mov    %ebx,%ecx
  mov    %r14,%rdx
  mov    %r14,%rsi
  mov    0x60(%rsp),%rbx
  mov    %rbx,%rdi
  call   <madd_stride@plt>
  add    $0x20,%rsp
  mov    0x18(%rsp),%rdi
  call   <free@plt>
  mov    %r15,%rdi
  call   <free@plt>
  mov    %rbx,%rdi
  call   <free@plt>
L92:
  add    $0x98,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
```

assistant: ```c
void srecmult(double *A, double *B, double *C, int m, int n, int p, int sA, int sB, int sC) {
  register int i, j, k;
  int u, v, t;
  double sum;
  double *A1, *B1;
  double *a11, *a12, *a21, *a22;
  double *b11, *b12, *b21, *b22;
  double *c11, *c12, *c21, *c22;
  double *m1, *m2, *m3, *m4, *m5, *m6, *m7;
  int sm1, sm2, sm3, sm4, sm5, sm6, sm7;
  int sA1, sB1;
  if (m + n + p <= 192) {
    for (i = 0; i < m; ++i) {
      for (j = 0; j < p; ++j) {
        v = i * sA;
        u = i * sC;
        t = j + u;
        sum = 0.;
        for (k = 0; k < n; ++k) {
          sum += A[k + v] * B[j + k * sB];
        }
        C[t] = sum;
      }
    }
  } else {
    m /= 2;
    n /= 2;
    p /= 2;
    a11 = A;
    a12 = A + n;
    a21 = A + m * sA;
    a22 = A + n + m * sA;
    b11 = B;
    b12 = B + p;
    b21 = B + n * sB;
    b22 = B + p + n * sB;
    c11 = C;
    c12 = C + p;
    c21 = C + m * sC;
    c22 = C + p + m * sC;
    m1 = (double *)malloc(sizeof(double) * m * p);
    sm1 = p;
    m3 = (double *)malloc(sizeof(double) * m * p);
    sm3 = p;
    m4 = (double *)malloc(sizeof(double) * m * p);
    sm4 = p;
    m2 = c21;
    sm2 = sC;
    m5 = c12;
    sm5 = sC;
    m6 = c22;
    sm6 = sC;
    m7 = c11;
    sm7 = sC;
    sA1 = n;
    sB1 = p;
    A1 = (double *)malloc(sizeof(double) * m * n);
    B1 = (double *)malloc(sizeof(double) * n * p);
    madd_stride(a11, a22, A1, m, n, sA, sA, sA1);
    madd_stride(b11, b22, B1, n, p, sB, sB, sB1);
    srecmult(A1, B1, m1, m, n, p, sA1, sB1, sm1);
    free(A1);
    free(B1);
    A1 = (double *)malloc(sizeof(double) * m * n);
    madd_stride(a21, a22, A1, m, n, sA, sA, sA1);
    srecmult(A1, b11, m2, m, n, p, sA1, sB, sm2);
    free(A1);
    B1 = (double *)malloc(sizeof(double) * n * p);
    msub_stride(b12, b22, B1, n, p, sB, sB, sB1);
    srecmult(a11, B1, m3, m, n, p, sA, sB1, sm3);
    free(B1);
    B1 = (double *)malloc(sizeof(double) * n * p);
    msub_stride(b21, b11, B1, n, p, sB, sB, sB1);
    srecmult(a22, B1, m4, m, n, p, sA, sB1, sm4);
    free(B1);
    A1 = (double *)malloc(sizeof(double) * m * n);
    madd_stride(a11, a12, A1, m, n, sA, sA, sA1);
    srecmult(A1, b22, m5, m, n, p, sA1, sB, sm5);
    free(A1);
    A1 = (double *)malloc(sizeof(double) * m * n);
    B1 = (double *)malloc(sizeof(double) * n * p);
    msub_stride(a21, a11, A1, m, n, sA, sA, sA1);
    madd_stride(b11, b12, B1, n, p, sB, sB, sB1);
    srecmult(A1, B1, m6, m, n, p, sA1, sB1, sm6);
    free(A1);
    free(B1);
    A1 = (double *)malloc(sizeof(double) * m * n);
    B1 = (double *)malloc(sizeof(double) * n * p);
    msub_stride(a12, a22, A1, m, n, sA, sA, sA1);
    madd_stride(b21, b22, B1, n, p, sB, sB, sB1);
    srecmult(A1, B1, m7, m, n, p, sA1, sB1, sm7);
    free(A1);
    free(B1);
    A1 = (double *)malloc(sizeof(double) * m * p);
    sA1 = p;
    madd_stride(m1, m7, m7, m, p, sm1, sm7, sm7);
    msub_stride(m4, m5, A1, m, p, sm4, sm5, sA1);
    madd_stride(m7, A1, m7, m, p, sm7, sA1, sm7);
    free(A1);
    A1 = (double *)malloc(sizeof(double) * m * p);
    sA1 = p;
    madd_stride(m1, m6, m6, m, p, sm1, sm6, sm6);
    msub_stride(m3, m2, A1, m, p, sm3, sm2, sA1);
    madd_stride(m6, A1, m6, m, p, sm6, sA1, sm6);
    free(A1);
    madd_stride(m3, m5, m5, m, p, sm3, sm5, sm5);
    madd_stride(m4, m2, m2, m, p, sm4, sm2, sm2);
    free(m1);
    free(m3);
    free(m4);
  }
}
```