user: Please show me the C implementation that matches the assembly code listed below:

```asm
<test_mm_njnini_b_c_bc_i4_d>:
  endbr64
  push   %rbp
  mov    %rsp,%rbp
  push   %r15
  push   %r14
  push   %r13
  push   %r12
  push   %rbx
  sub    $0x168,%rsp
  mov    %fs:0x28,%rax
  mov    %rax,-0x38(%rbp)
  xor    %eax,%eax
  mov    %rsp,%rax
  mov    %rax,-0x188(%rbp)
  movl   $0x15,-0xcc(%rbp)
  movl   $0x17,-0xc8(%rbp)
  movl   $0x19,-0xc4(%rbp)
  mov    -0xc8(%rbp),%r10d
  mov    -0xcc(%rbp),%r11d
  movslq %r10d,%rax
  sub    $0x1,%rax
  mov    %rax,-0xc0(%rbp)
  movslq %r10d,%rax
  mov    %rax,%rcx
  mov    $0x0,%ebx
  movslq %r10d,%rax
  lea    0x0(,%rax,4),%rbx
  movslq %r11d,%rax
  sub    $0x1,%rax
  mov    %rax,-0xb8(%rbp)
  movslq %r10d,%rax
  mov    %rax,%r8
  mov    $0x0,%r9d
  movslq %r11d,%rax
  mov    %rax,%rsi
  mov    $0x0,%edi
  mov    %r9,%rdx
  imul   %rsi,%rdx
  mov    %rdi,%rax
  imul   %r8,%rax
  lea    (%rdx,%rax,1),%rcx
  mov    %r8,%rax
  mul    %rsi
  add    %rdx,%rcx
  mov    %rcx,%rdx
  movslq %r10d,%rax
  mov    %rax,%r14
  mov    $0x0,%r15d
  movslq %r11d,%rax
  mov    %rax,%r12
  mov    $0x0,%r13d
  mov    %r15,%rdx
  imul   %r12,%rdx
  mov    %r13,%rax
  imul   %r14,%rax
  lea    (%rdx,%rax,1),%rcx
  mov    %r14,%rax
  mul    %r12
  add    %rdx,%rcx
  mov    %rcx,%rdx
  movslq %r10d,%rdx
  movslq %r11d,%rax
  imul   %rdx,%rax
  lea    0x0(,%rax,4),%rdx
  mov    $0x10,%eax
  sub    $0x1,%rax
  add    %rdx,%rax
  mov    $0x10,%edi
  mov    $0x0,%edx
  div    %rdi
  imul   $0x10,%rax,%rdx
  mov    %rdx,%rax
  and    $0xfffffffffffff000,%rax
  mov    %rsp,%rcx
  sub    %rax,%rcx
L234:
  cmp    %rcx,%rsp
  je     L39
  sub    $0x1000,%rsp
  orq    $0x0,0xff8(%rsp)
  jmp    L234
L39:
  mov    %rdx,%rax
  and    $0xfff,%eax
  sub    %rax,%rsp
  mov    %rdx,%rax
  and    $0xfff,%eax
  test   %rax,%rax
  je     L184
  mov    %rdx,%rax
  and    $0xfff,%eax
  sub    $0x8,%rax
  add    %rsp,%rax
  orq    $0x0,(%rax)
L184:
  mov    %rsp,%rax
  add    $0x3,%rax
  shr    $0x2,%rax
  shl    $0x2,%rax
  mov    %rax,-0xb0(%rbp)
  mov    -0xc4(%rbp),%esi
  mov    -0xc8(%rbp),%edi
  movslq %esi,%rax
  sub    $0x1,%rax
  mov    %rax,-0xa8(%rbp)
  movslq %esi,%rax
  mov    %rax,-0x170(%rbp)
  movq   $0x0,-0x168(%rbp)
  movslq %esi,%rax
  lea    0x0(,%rax,4),%r8
  movslq %edi,%rax
  sub    $0x1,%rax
  mov    %rax,-0xa0(%rbp)
  movslq %esi,%rax
  mov    %rax,-0xf0(%rbp)
  movq   $0x0,-0xe8(%rbp)
  movslq %edi,%rax
  mov    %rax,-0x100(%rbp)
  movq   $0x0,-0xf8(%rbp)
  mov    -0xf0(%rbp),%r14
  mov    -0xe8(%rbp),%r15
  mov    %r15,%rdx
  mov    -0x100(%rbp),%r9
  mov    -0xf8(%rbp),%r10
  imul   %r9,%rdx
  mov    %r10,%rax
  imul   %r14,%rax
  lea    (%rdx,%rax,1),%rcx
  mov    %r14,%rax
  mul    %r9
  add    %rdx,%rcx
  mov    %rcx,%rdx
  movslq %esi,%rax
  mov    %rax,-0x110(%rbp)
  movq   $0x0,-0x108(%rbp)
  movslq %edi,%rax
  mov    %rax,-0x120(%rbp)
  movq   $0x0,-0x118(%rbp)
  mov    -0x110(%rbp),%r9
  mov    -0x108(%rbp),%r10
  mov    %r10,%rdx
  mov    -0x120(%rbp),%r11
  mov    -0x118(%rbp),%r12
  imul   %r11,%rdx
  mov    %r12,%rax
  imul   %r9,%rax
  lea    (%rdx,%rax,1),%rcx
  mov    %r9,%rax
  mul    %r11
  add    %rdx,%rcx
  mov    %rcx,%rdx
  movslq %esi,%rdx
  movslq %edi,%rax
  imul   %rdx,%rax
  lea    0x0(,%rax,4),%rdx
  mov    $0x10,%eax
  sub    $0x1,%rax
  add    %rdx,%rax
  mov    $0x10,%edi
  mov    $0x0,%edx
  div    %rdi
  imul   $0x10,%rax,%rdx
  mov    %rdx,%rax
  and    $0xfffffffffffff000,%rax
  mov    %rsp,%rcx
  sub    %rax,%rcx
L186:
  cmp    %rcx,%rsp
  je     L247
  sub    $0x1000,%rsp
  orq    $0x0,0xff8(%rsp)
  jmp    L186
L247:
  mov    %rdx,%rax
  and    $0xfff,%eax
  sub    %rax,%rsp
  mov    %rdx,%rax
  and    $0xfff,%eax
  test   %rax,%rax
  je     L71
  mov    %rdx,%rax
  and    $0xfff,%eax
  sub    $0x8,%rax
  add    %rsp,%rax
  orq    $0x0,(%rax)
L71:
  mov    %rsp,%rax
  add    $0x3,%rax
  shr    $0x2,%rax
  shl    $0x2,%rax
  mov    %rax,-0x98(%rbp)
  mov    -0xc4(%rbp),%esi
  mov    -0xcc(%rbp),%edi
  movslq %esi,%rax
  sub    $0x1,%rax
  mov    %rax,-0x90(%rbp)
  movslq %esi,%rax
  mov    %rax,-0x180(%rbp)
  movq   $0x0,-0x178(%rbp)
  movslq %esi,%rax
  lea    0x0(,%rax,4),%r13
  movslq %edi,%rax
  sub    $0x1,%rax
  mov    %rax,-0x88(%rbp)
  movslq %esi,%rax
  mov    %rax,-0x130(%rbp)
  movq   $0x0,-0x128(%rbp)
  movslq %edi,%rax
  mov    %rax,-0x140(%rbp)
  movq   $0x0,-0x138(%rbp)
  mov    -0x130(%rbp),%r9
  mov    -0x128(%rbp),%r10
  mov    %r10,%rdx
  mov    -0x140(%rbp),%r11
  mov    -0x138(%rbp),%r12
  imul   %r11,%rdx
  mov    %r12,%rax
  imul   %r9,%rax
  lea    (%rdx,%rax,1),%rcx
  mov    %r9,%rax
  mul    %r11
  add    %rdx,%rcx
  mov    %rcx,%rdx
  movslq %esi,%rax
  mov    %rax,-0x150(%rbp)
  movq   $0x0,-0x148(%rbp)
  movslq %edi,%rax
  mov    %rax,-0x160(%rbp)
  movq   $0x0,-0x158(%rbp)
  mov    -0x150(%rbp),%r11
  mov    -0x148(%rbp),%r12
  mov    %r12,%rdx
  mov    -0x160(%rbp),%r14
  mov    -0x158(%rbp),%r15
  imul   %r14,%rdx
  mov    %r15,%rax
  imul   %r11,%rax
  lea    (%rdx,%rax,1),%rcx
  mov    %r11,%rax
  mul    %r14
  add    %rdx,%rcx
  mov    %rcx,%rdx
  movslq %esi,%rdx
  movslq %edi,%rax
  imul   %rdx,%rax
  lea    0x0(,%rax,4),%rdx
  mov    $0x10,%eax
  sub    $0x1,%rax
  add    %rdx,%rax
  mov    $0x10,%edi
  mov    $0x0,%edx
  div    %rdi
  imul   $0x10,%rax,%rax
  mov    %rax,%rcx
  and    $0xfffffffffffff000,%rcx
  mov    %rsp,%rdx
  sub    %rcx,%rdx
L35:
  cmp    %rdx,%rsp
  je     L210
  sub    $0x1000,%rsp
  orq    $0x0,0xff8(%rsp)
  jmp    L35
L210:
  mov    %rax,%rdx
  and    $0xfff,%edx
  sub    %rdx,%rsp
  mov    %rax,%rdx
  and    $0xfff,%edx
  test   %rdx,%rdx
  je     L251
  and    $0xfff,%eax
  sub    $0x8,%rax
  add    %rsp,%rax
  orq    $0x0,(%rax)
L251:
  mov    %rsp,%rax
  add    $0x3,%rax
  shr    $0x2,%rax
  shl    $0x2,%rax
  mov    %rax,-0x80(%rbp)
  movl   $0x0,-0xd8(%rbp)
  jmp    L10
L63:
  movl   $0x0,-0xd4(%rbp)
  jmp    L88
L4:
  mov    -0xd8(%rbp),%eax
  imul   -0xcc(%rbp),%eax
  mov    %eax,%edx
  mov    -0xd4(%rbp),%eax
  add    %edx,%eax
  mov    %rbx,%rsi
  shr    $0x2,%rsi
  lea    0x1(%rax),%ecx
  mov    -0xb0(%rbp),%rax
  mov    -0xd4(%rbp),%edx
  movslq %edx,%rdi
  mov    -0xd8(%rbp),%edx
  movslq %edx,%rdx
  imul   %rsi,%rdx
  add    %rdi,%rdx
  mov    %ecx,(%rax,%rdx,4)
  addl   $0x1,-0xd4(%rbp)
L88:
  mov    -0xd4(%rbp),%eax
  cmp    -0xc8(%rbp),%eax
  jl     L4
  addl   $0x1,-0xd8(%rbp)
L10:
  mov    -0xd8(%rbp),%eax
  cmp    -0xcc(%rbp),%eax
  jl     L63
  movl   $0x0,-0xd8(%rbp)
  jmp    L146
L76:
  movl   $0x0,-0xd4(%rbp)
  jmp    L69
L196:
  mov    -0xd4(%rbp),%eax
  imul   -0xc4(%rbp),%eax
  mov    %eax,%edx
  mov    -0xd8(%rbp),%eax
  add    %edx,%eax
  mov    %r8,%rsi
  shr    $0x2,%rsi
  lea    0x1(%rax),%ecx
  mov    -0x98(%rbp),%rax
  mov    -0xd4(%rbp),%edx
  movslq %edx,%rdi
  mov    -0xd8(%rbp),%edx
  movslq %edx,%rdx
  imul   %rsi,%rdx
  add    %rdi,%rdx
  mov    %ecx,(%rax,%rdx,4)
  addl   $0x1,-0xd4(%rbp)
L69:
  mov    -0xd4(%rbp),%eax
  cmp    -0xc4(%rbp),%eax
  jl     L196
  addl   $0x1,-0xd8(%rbp)
L146:
  mov    -0xd8(%rbp),%eax
  cmp    -0xc8(%rbp),%eax
  jl     L76
  movl   $0x0,-0xd8(%rbp)
  jmp    L171
L90:
  movl   $0x0,-0xd4(%rbp)
  jmp    L61
L96:
  mov    %r13,%rcx
  shr    $0x2,%rcx
  mov    -0x80(%rbp),%rax
  mov    -0xd4(%rbp),%edx
  movslq %edx,%rsi
  mov    -0xd8(%rbp),%edx
  movslq %edx,%rdx
  imul   %rcx,%rdx
  add    %rsi,%rdx
  movl   $0x0,(%rax,%rdx,4)
  addl   $0x1,-0xd4(%rbp)
L61:
  mov    -0xd4(%rbp),%eax
  cmp    -0xc4(%rbp),%eax
  jl     L96
  addl   $0x1,-0xd8(%rbp)
L171:
  mov    -0xd8(%rbp),%eax
  cmp    -0xcc(%rbp),%eax
  jl     L90
  mov    -0x98(%rbp),%rax
  mov    %rax,%rdi
  call   <xmp_desc_of@plt>
  mov    %eax,%r12d
  mov    -0xb0(%rbp),%rax
  mov    %rax,%rdi
  call   <xmp_desc_of@plt>
  mov    %eax,%ebx
  mov    -0x80(%rbp),%rax
  mov    %rax,%rdi
  call   <xmp_desc_of@plt>
  mov    %r12d,%edx
  mov    %ebx,%esi
  mov    %eax,%edi
  call   <xmp_matmul@plt>
  movl   $0x0,-0xd0(%rbp)
  pxor   %xmm0,%xmm0
  cvtsi2sdl -0xcc(%rbp),%xmm0
  movsd  %xmm0,-0x78(%rbp)
  pxor   %xmm0,%xmm0
  cvtsi2sdl -0xc8(%rbp),%xmm0
  movsd  %xmm0,-0x70(%rbp)
  pxor   %xmm0,%xmm0
  cvtsi2sdl -0xc4(%rbp),%xmm0
  movsd  %xmm0,-0x68(%rbp)
  movsd  -0x70(%rbp),%xmm1
  movsd  D44(%rip),%xmm0
  addsd  %xmm1,%xmm0
  movapd %xmm0,%xmm1
  mulsd  -0x70(%rbp),%xmm1
  movsd  -0x70(%rbp),%xmm0
  movapd %xmm0,%xmm2
  addsd  %xmm0,%xmm2
  movsd  D44(%rip),%xmm0
  addsd  %xmm2,%xmm0
  mulsd  %xmm1,%xmm0
  movsd  D18(%rip),%xmm1
  divsd  %xmm1,%xmm0
  movsd  %xmm0,-0x60(%rbp)
  movsd  -0x70(%rbp),%xmm1
  movsd  D44(%rip),%xmm0
  addsd  %xmm1,%xmm0
  mulsd  -0x70(%rbp),%xmm0
  movsd  D1(%rip),%xmm1
  divsd  %xmm1,%xmm0
  movsd  %xmm0,-0x58(%rbp)
  movl   $0x0,-0xd8(%rbp)
  jmp    L192
L257:
  movl   $0x0,-0xd4(%rbp)
  jmp    L15
L216:
  pxor   %xmm0,%xmm0
  cvtsi2sdl -0xd8(%rbp),%xmm0
  movsd  -0x78(%rbp),%xmm1
  mulsd  %xmm1,%xmm0
  movsd  %xmm0,-0x50(%rbp)
  pxor   %xmm0,%xmm0
  cvtsi2sdl -0xd4(%rbp),%xmm0
  movsd  -0x68(%rbp),%xmm1
  mulsd  %xmm1,%xmm0
  movsd  %xmm0,-0x48(%rbp)
  movsd  -0x50(%rbp),%xmm0
  addsd  -0x48(%rbp),%xmm0
  mulsd  -0x58(%rbp),%xmm0
  movapd %xmm0,%xmm1
  addsd  -0x60(%rbp),%xmm1
  movsd  -0x50(%rbp),%xmm0
  mulsd  -0x48(%rbp),%xmm0
  mulsd  -0x70(%rbp),%xmm0
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x40(%rbp)
  pxor   %xmm1,%xmm1
  cvtsi2sdl -0xd0(%rbp),%xmm1
  mov    %r13,%rcx
  shr    $0x2,%rcx
  mov    -0x80(%rbp),%rax
  mov    -0xd4(%rbp),%edx
  movslq %edx,%rsi
  mov    -0xd8(%rbp),%edx
  movslq %edx,%rdx
  imul   %rcx,%rdx
  add    %rsi,%rdx
  mov    (%rax,%rdx,4),%eax
  pxor   %xmm0,%xmm0
  cvtsi2sd %eax,%xmm0
  subsd  -0x40(%rbp),%xmm0
  addsd  %xmm1,%xmm0
  cvttsd2si %xmm0,%eax
  mov    %eax,-0xd0(%rbp)
  addl   $0x1,-0xd4(%rbp)
L15:
  mov    -0xd4(%rbp),%eax
  cmp    -0xc4(%rbp),%eax
  jl     L216
  addl   $0x1,-0xd8(%rbp)
L192:
  mov    -0xd8(%rbp),%eax
  cmp    -0xcc(%rbp),%eax
  jl     L257
  mov    -0xd0(%rbp),%eax
  mov    %eax,%edi
  call   <chk_int2@plt>
  mov    -0x188(%rbp),%rsp
  nop
  mov    -0x38(%rbp),%rax
  sub    %fs:0x28,%rax
  je     L60
  call   <__stack_chk_fail@plt>
L60:
  lea    -0x28(%rbp),%rsp
  pop    %rbx
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  pop    %rbp
  ret
```

tool_calls: name: retrieve_data
arguments: data_label: D44
type_infered: byte[8]
name: retrieve_data
arguments: data_label: D18
type_infered: f64
name: retrieve_data
arguments: data_label: D1
type_infered: f64

tool: {"D44": [0x00,0x00,0x00,0x00,0x00,0x00,0xf0,0x3f]}

tool: D18: 6

tool: D1: 2

assistant: ```c
void test_mm_njnini_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][*] with tx(i, *)
#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;
#pragma xmp loop(i) on tx(i, *)
  for (i = 0; i < n1; i++) {
    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); }
}
```