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

```asm
<test_mm_ninini_b_c_bc_i4_d>:
  endbr64
  push   %rbp
  mov    %rsp,%rbp
  push   %r14
  push   %r13
  push   %r12
  push   %rbx
  sub    $0x10,%rsp
  mov    %fs:0x28,%rax
  mov    %rax,-0x28(%rbp)
  xor    %eax,%eax
  mov    %rsp,%rax
L125:
  cmp    %rax,%rsp
  je     L16
  sub    $0x1000,%rsp
  orq    $0x0,0xff8(%rsp)
  jmp    L125
L16:
  sub    $0x790,%rsp
  orq    $0x0,0x788(%rsp)
  mov    %rsp,%r13
  mov    %rsp,%rax
L84:
  cmp    %rax,%rsp
  je     L78
  sub    $0x1000,%rsp
  orq    $0x0,0xff8(%rsp)
  jmp    L84
L78:
  sub    $0x900,%rsp
  orq    $0x0,0x8f8(%rsp)
  mov    %rsp,%rdi
  mov    %rsp,%rax
L69:
  cmp    %rax,%rsp
  je     L28
  sub    $0x1000,%rsp
  orq    $0x0,0xff8(%rsp)
  jmp    L69
L28:
  sub    $0x840,%rsp
  orq    $0x0,0x838(%rsp)
  mov    %rsp,%r12
  mov    %r13,%rsi
  mov    $0x18,%ecx
L39:
  lea    -0x17(%rcx),%eax
  mov    %rsi,%rdx
L85:
  mov    %eax,(%rdx)
  add    $0x1,%eax
  add    $0x4,%rdx
  cmp    %ecx,%eax
  jne    L85
  add    $0x5c,%rsi
  add    $0x15,%ecx
  cmp    $0x1d1,%ecx
  jne    L39
  mov    %rdi,%rsi
  mov    $0x272,%ecx
L56:
  lea    -0x271(%rcx),%eax
  mov    %rsi,%rdx
L140:
  mov    %eax,(%rdx)
  add    $0x19,%eax
  add    $0x4,%rdx
  cmp    %ecx,%eax
  jne    L140
  add    $0x64,%rsi
  add    $0x1,%ecx
  cmp    $0x289,%ecx
  jne    L56
  mov    %r12,%rbx
  lea    0x64(%r12),%rdx
  lea    0x898(%r12),%rcx
L31:
  lea    -0x64(%rdx),%rax
L65:
  movl   $0x0,(%rax)
  add    $0x4,%rax
  cmp    %rdx,%rax
  jne    L65
  add    $0x64,%rdx
  cmp    %rcx,%rdx
  jne    L31
  call   <xmp_desc_of@plt>
  mov    %eax,%r14d
  mov    %r13,%rdi
  call   <xmp_desc_of@plt>
  mov    %eax,%r13d
  mov    %r12,%rdi
  call   <xmp_desc_of@plt>
  mov    %eax,%edi
  mov    %r14d,%edx
  mov    %r13d,%esi
  call   <xmp_matmul@plt>
  mov    $0x0,%edi
  mov    $0x0,%edx
  movsd  D7(%rip),%xmm7
  movsd  D41(%rip),%xmm6
  movsd  D44(%rip),%xmm5
  movsd  D21(%rip),%xmm4
  movsd  D35(%rip),%xmm3
L2:
  pxor   %xmm2,%xmm2
  cvtsi2sd %edx,%xmm2
  mulsd  %xmm7,%xmm2
  mov    $0x0,%eax
L129:
  pxor   %xmm0,%xmm0
  cvtsi2sd %eax,%xmm0
  mulsd  %xmm6,%xmm0
  pxor   %xmm8,%xmm8
  cvtsi2sdl (%rbx,%rax,4),%xmm8
  movapd %xmm2,%xmm1
  addsd  %xmm0,%xmm1
  mulsd  %xmm5,%xmm1
  addsd  %xmm4,%xmm1
  mulsd  %xmm2,%xmm0
  mulsd  %xmm3,%xmm0
  addsd  %xmm1,%xmm0
  movapd %xmm8,%xmm1
  subsd  %xmm0,%xmm1
  pxor   %xmm0,%xmm0
  cvtsi2sd %edi,%xmm0
  addsd  %xmm1,%xmm0
  cvttsd2si %xmm0,%edi
  add    $0x1,%rax
  cmp    $0x19,%rax
  jne    L129
  add    $0x1,%edx
  add    $0x64,%rbx
  cmp    $0x15,%edx
  jne    L2
  call   <chk_int2@plt>
  mov    -0x28(%rbp),%rax
  sub    %fs:0x28,%rax
  jne    L135
  lea    -0x20(%rbp),%rsp
  pop    %rbx
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %rbp
  ret
L135:
  call   <__stack_chk_fail@plt>
```

tool_calls: name: unpack_data
arguments: data_tag: D7
type_infered: byte[8]
name: unpack_data
arguments: data_tag: D41
type_infered: byte[8]
name: unpack_data
arguments: data_tag: D44
type_infered: byte[8]
name: unpack_data
arguments: data_tag: D21
type_infered: byte[8]
name: unpack_data
arguments: data_tag: D35
type_infered: byte[8]

tool: {"D7": [0x00,0x00,0x00,0x00,0x00,0x00,0x35,0x40]}

tool: {"D41": [0x00,0x00,0x00,0x00,0x00,0x00,0x39,0x40]}

tool: {"D44": [0x00,0x00,0x00,0x00,0x00,0x40,0x71,0x40]}

tool: {"D21": [0x00,0x00,0x00,0x00,0x00,0xe4,0xb0,0x40]}

tool: {"D35": [0x00,0x00,0x00,0x00,0x00,0x00,0x37,0x40]}

assistant: ```c
void test_mm_ninini_b_c_bc_i4_d() {
  int n1 = 21, n2 = 23, n3 = 25;
  int a[n1][n2], b[n2][n3], x[n1][n3];
#pragma xmp nodes p(2, 2)
#pragma xmp template tx(0 : n1 - 1, 0 : n2 - 1)
#pragma xmp template ty(0 : n2 - 1, 0 : n3 - 1)
#pragma xmp template tz(0 : n1 - 1, 0 : n3 - 1)
#pragma xmp distribute tx(block, block) onto p
#pragma xmp distribute ty(cyclic, cyclic) onto p
#pragma xmp distribute tz(cyclic(2), cyclic(2)) onto p
#pragma xmp align a[*][j] with tx(*, j)
#pragma xmp align b[*][j] with ty(*, j)
#pragma xmp align x[*][j] with tz(*, j)
  int i, j, ierr;
  double rn1, rn2, rn3, rn4, rn5, rn6, ra, rb;
  for (i = 0; i < n1; i++) {
#pragma xmp loop(j) on tx(*, j)
    for (j = 0; j < n2; j++) {
      a[i][j] = i * n1 + j + 1;
    }
  }
  for (i = 0; i < n2; i++) {
#pragma xmp loop(j) on ty(*, j)
    for (j = 0; j < n3; j++) {
      b[i][j] = j * n3 + i + 1;
    }
  }
  for (i = 0; i < n1; i++) {
#pragma xmp loop(j) on tz(*, j)
    for (j = 0; j < n3; j++) {
      x[i][j] = 0;
    }
  }
  xmp_matmul(xmp_desc_of(x), xmp_desc_of(a), xmp_desc_of(b));
  ierr = 0;
  rn1 = n1;
  rn2 = n2;
  rn3 = n3;
  rn4 = rn2 * (rn2 + 1) * (2 * rn2 + 1) / 6.0;
  rn5 = rn2 * (rn2 + 1) / 2.0;
  for (i = 0; i < n1; i++) {
#pragma xmp loop(j) on tz(*, j)
    for (j = 0; j < n3; j++) {
      ra = i * rn1;
      rb = j * rn3;
      rn6 = rn4 + rn5 * (ra + rb) + ra * rb * rn2;
      ierr = ierr + (x[i][j] - rn6);
    }
  }
#pragma xmp task on p(1, 1 : 2)
  { chk_int2(ierr); }
}
```