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

```asm
<test_mm_aanj_b_c_bc_i4_d>:
  endbr64
  push   %rbp
  mov    %rsp,%rbp
  push   %r13
  push   %r12
  push   %rbx
  sub    $0x18,%rsp
  mov    %fs:0x28,%rax
  mov    %rax,-0x28(%rbp)
  mov    %rsp,%rax
  cmp    %rax,%rsp
  je     L16
L79:
  sub    $0x1000,%rsp
  orq    $0x0,0xff8(%rsp)
  cmp    %rax,%rsp
  jne    L79
L16:
  sub    $0x790,%rsp
  orq    $0x0,0x788(%rsp)
  mov    %rsp,%rax
  mov    %rsp,%r12
  cmp    %rax,%rsp
  je     L27
L48:
  sub    $0x1000,%rsp
  orq    $0x0,0xff8(%rsp)
  cmp    %rax,%rsp
  jne    L48
L27:
  sub    $0x900,%rsp
  orq    $0x0,0x8f8(%rsp)
  mov    %rsp,%rax
  mov    %rsp,%r8
  cmp    %rax,%rsp
  je     L26
L8:
  sub    $0x1000,%rsp
  orq    $0x0,0xff8(%rsp)
  cmp    %rax,%rsp
  jne    L8
L26:
  sub    $0x840,%rsp
  orq    $0x0,0x838(%rsp)
  movdqa D74(%rip),%xmm6
  movdqa D97(%rip),%xmm5
  movdqa D7(%rip),%xmm4
  movdqa D64(%rip),%xmm3
  movdqa D42(%rip),%xmm2
  lea    0x3(%rsp),%rdx
  mov    %r12,%rax
  mov    %rdx,%rsi
  mov    %rdx,%rbx
  mov    $0x15,%edx
  shr    $0x2,%rsi
  and    $0xfffffffffffffffc,%rbx
  xchg   %ax,%ax
L108:
  lea    -0x15(%rdx),%edi
  lea    0x1(%rdx),%ecx
  mov    %edx,0x50(%rax)
  add    $0x5c,%rax
  movd   %edi,%xmm0
  mov    %ecx,-0x8(%rax)
  lea    0x2(%rdx),%ecx
  add    $0x15,%edx
  pshufd $0x0,%xmm0,%xmm0
  mov    %ecx,-0x4(%rax)
  movdqa %xmm0,%xmm1
  paddd  %xmm6,%xmm1
  movups %xmm1,-0x5c(%rax)
  movdqa %xmm0,%xmm1
  paddd  %xmm5,%xmm1
  movups %xmm1,-0x4c(%rax)
  movdqa %xmm0,%xmm1
  paddd  %xmm4,%xmm1
  movups %xmm1,-0x3c(%rax)
  movdqa %xmm0,%xmm1
  paddd  %xmm2,%xmm0
  paddd  %xmm3,%xmm1
  movups %xmm0,-0x1c(%rax)
  movups %xmm1,-0x2c(%rax)
  cmp    $0x1ce,%edx
  jne    L108
  movdqa D9(%rip),%xmm7
  movdqa D44(%rip),%xmm6
  mov    %r8,%rax
  xor    %edx,%edx
  movdqa D123(%rip),%xmm5
  movdqa D119(%rip),%xmm4
  movdqa D12(%rip),%xmm3
  movdqa D137(%rip),%xmm2
  nopw   0x0(%rax,%rax,1)
L87:
  movd   %edx,%xmm0
  lea    0x259(%rdx),%ecx
  add    $0x1,%edx
  add    $0x64,%rax
  pshufd $0x0,%xmm0,%xmm0
  mov    %ecx,-0x4(%rax)
  movdqa %xmm0,%xmm1
  paddd  %xmm7,%xmm1
  movups %xmm1,-0x64(%rax)
  movdqa %xmm0,%xmm1
  paddd  %xmm6,%xmm1
  movups %xmm1,-0x54(%rax)
  movdqa %xmm0,%xmm1
  paddd  %xmm5,%xmm1
  movups %xmm1,-0x44(%rax)
  movdqa %xmm0,%xmm1
  paddd  %xmm4,%xmm1
  movups %xmm1,-0x34(%rax)
  movdqa %xmm0,%xmm1
  paddd  %xmm2,%xmm0
  paddd  %xmm3,%xmm1
  movups %xmm0,-0x14(%rax)
  movups %xmm1,-0x24(%rax)
  cmp    $0x17,%edx
  jne    L87
  lea    0x8(%rbx),%rdi
  mov    %rbx,%rcx
  xor    %eax,%eax
  movq   $0x0,0x0(,%rsi,4)
  movq   $0x0,0x82c(%rbx)
  and    $0xfffffffffffffff8,%rdi
  sub    %rdi,%rcx
  add    $0x834,%ecx
  shr    $0x3,%ecx
  rep stos %rax,%es:(%rdi)
  mov    %r8,%rdi
  call   <xmp_desc_of@plt>
  mov    %r12,%rdi
  mov    %eax,%r13d
  call   <xmp_desc_of@plt>
  mov    %rbx,%rdi
  mov    %eax,%r12d
  call   <xmp_desc_of@plt>
  mov    %r13d,%edx
  mov    %r12d,%esi
  mov    %eax,%edi
  call   <xmp_matmul@plt>
  movsd  D69(%rip),%xmm7
  mov    %rbx,%rdx
  movsd  D45(%rip),%xmm6
  movsd  D89(%rip),%xmm5
  movsd  D51(%rip),%xmm4
  xor    %edi,%edi
  xor    %ecx,%ecx
  movsd  D116(%rip),%xmm8
L3:
  pxor   %xmm3,%xmm3
  xor    %eax,%eax
  cvtsi2sd %ecx,%xmm3
  mulsd  %xmm8,%xmm3
  nop
L11:
  pxor   %xmm0,%xmm0
  movapd %xmm3,%xmm1
  pxor   %xmm2,%xmm2
  cvtsi2sd %eax,%xmm0
  cvtsi2sdl (%rdx,%rax,4),%xmm2
  add    $0x1,%rax
  mulsd  %xmm7,%xmm0
  addsd  %xmm0,%xmm1
  mulsd  %xmm3,%xmm0
  mulsd  %xmm6,%xmm1
  mulsd  %xmm4,%xmm0
  addsd  %xmm5,%xmm1
  addsd  %xmm1,%xmm0
  subsd  %xmm0,%xmm2
  pxor   %xmm0,%xmm0
  cvtsi2sd %edi,%xmm0
  movapd %xmm2,%xmm1
  addsd  %xmm0,%xmm1
  cvttsd2si %xmm1,%edi
  cmp    $0x19,%rax
  jne    L11
  add    $0x1,%ecx
  add    $0x64,%rdx
  cmp    $0x15,%ecx
  jne    L3
  call   <chk_int2@plt>
  mov    -0x28(%rbp),%rax
  sub    %fs:0x28,%rax
  jne    L57
  lea    -0x18(%rbp),%rsp
  pop    %rbx
  pop    %r12
  pop    %r13
  pop    %rbp
  ret
L57:
  call   <__stack_chk_fail@plt>
```

tool_calls: name: read_data
arguments: data_tag: D74
datatype: byte[16]
name: read_data
arguments: data_tag: D97
datatype: byte[16]
name: read_data
arguments: data_tag: D7
datatype: byte[16]
name: read_data
arguments: data_tag: D64
datatype: byte[16]
name: read_data
arguments: data_tag: D42
datatype: byte[16]
name: read_data
arguments: data_tag: D9
datatype: byte[16]
name: read_data
arguments: data_tag: D44
datatype: byte[16]
name: read_data
arguments: data_tag: D123
datatype: byte[16]
name: read_data
arguments: data_tag: D119
datatype: byte[16]
name: read_data
arguments: data_tag: D12
datatype: byte[16]
name: read_data
arguments: data_tag: D137
datatype: byte[16]
name: read_data
arguments: data_tag: D69
datatype: byte[8]
name: read_data
arguments: data_tag: D45
datatype: byte[8]
name: read_data
arguments: data_tag: D89
datatype: byte[8]
name: read_data
arguments: data_tag: D51
datatype: byte[8]
name: read_data
arguments: data_tag: D116
datatype: byte[8]

tool: {"D74": [0x01,0x00,0x00,0x00,0x02,0x00,0x00,0x00,0x03,0x00,0x00,0x00,0x04,0x00,0x00,0x00]}

tool: {"D97": [0x05,0x00,0x00,0x00,0x06,0x00,0x00,0x00,0x07,0x00,0x00,0x00,0x08,0x00,0x00,0x00]}

tool: {"D7": [0x09,0x00,0x00,0x00,0x0a,0x00,0x00,0x00,0x0b,0x00,0x00,0x00,0x0c,0x00,0x00,0x00]}

tool: {"D64": [0x0d,0x00,0x00,0x00,0x0e,0x00,0x00,0x00,0x0f,0x00,0x00,0x00,0x10,0x00,0x00,0x00]}

tool: {"D42": [0x11,0x00,0x00,0x00,0x12,0x00,0x00,0x00,0x13,0x00,0x00,0x00,0x14,0x00,0x00,0x00]}

tool: {"D9": [0x01,0x00,0x00,0x00,0x1a,0x00,0x00,0x00,0x33,0x00,0x00,0x00,0x4c,0x00,0x00,0x00]}

tool: {"D44": [0x65,0x00,0x00,0x00,0x7e,0x00,0x00,0x00,0x97,0x00,0x00,0x00,0xb0,0x00,0x00,0x00]}

tool: {"D123": [0xc9,0x00,0x00,0x00,0xe2,0x00,0x00,0x00,0xfb,0x00,0x00,0x00,0x14,0x01,0x00,0x00]}

tool: {"D119": [0x2d,0x01,0x00,0x00,0x46,0x01,0x00,0x00,0x5f,0x01,0x00,0x00,0x78,0x01,0x00,0x00]}

tool: {"D12": [0x91,0x01,0x00,0x00,0xaa,0x01,0x00,0x00,0xc3,0x01,0x00,0x00,0xdc,0x01,0x00,0x00]}

tool: {"D137": [0xf5,0x01,0x00,0x00,0x0e,0x02,0x00,0x00,0x27,0x02,0x00,0x00,0x40,0x02,0x00,0x00]}

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

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

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

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

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

assistant: ```c
void test_mm_aanj_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[i][j] with tx(i, j)
#pragma xmp align b[i][j] with ty(i, j)
#pragma xmp align x[i][*] with tz(i, *)
  int i, j, ierr;
  double rn1, rn2, rn3, rn4, rn5, rn6, ra, rb;
#pragma xmp loop(i, j) on tx(i, j)
  for (i = 0; i < n1; i++) {
    for (j = 0; j < n2; j++) {
      a[i][j] = i * n1 + j + 1;
    }
  }
#pragma xmp loop(i, j) on ty(i, j)
  for (i = 0; i < n2; i++) {
    for (j = 0; j < n3; j++) {
      b[i][j] = j * n3 + i + 1;
    }
  }
#pragma xmp loop(i) on tz(i, *)
  for (i = 0; i < n1; i++) {
    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;
#pragma xmp loop(i) on tz(i, *)
  for (i = 0; i < n1; i++) {
    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 : 2, 1)
  { chk_int2(ierr); }
}
```