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

```asm
<kendl2>:
  endbr64
  push   %r15
  movslq %esi,%rax
  mov    %rdi,%r11
  movslq %edx,%rdi
  push   %r14
  push   %r13
  push   %r12
  mov    %r8,%r12
  push   %rbp
  mov    %rcx,%rbp
  mov    %eax,%ecx
  push   %rbx
  imul   %edx,%ecx
  sub    $0x18,%rsp
  mov    (%r11,%rax,8),%rax
  mov    %r9,(%rsp)
  movss  (%rax,%rdi,4),%xmm5
  cmp    $0x1,%ecx
  jle    L36
  pxor   %xmm2,%xmm2
  movslq %ecx,%r8
  xor    %ebx,%ebx
  lea    -0x1(%r8),%r13
  movaps %xmm2,%xmm0
  movaps %xmm2,%xmm4
  movss  D34(%rip),%xmm6
  nopw   0x0(%rax,%rax,1)
L46:
  mov    %rbx,%rax
  add    $0x1,%rbx
  cqto
  idiv   %rdi
  mov    %rax,%r9
  mov    0x8(%r11,%rax,8),%rax
  mov    %rdx,%r10
  movss  0x4(%rax,%rdx,4),%xmm3
  addss  %xmm3,%xmm5
  cmp    %rbx,%r8
  jle    L87
  mov    %rbx,%rcx
  jmp    L41
  nopl   0x0(%rax,%rax,1)
L17:
  addss  %xmm1,%xmm2
  addss  %xmm1,%xmm0
  jg     L25
  xorps  %xmm6,%xmm1
L25:
  addss  %xmm1,%xmm4
L57:
  add    $0x1,%rcx
  cmp    %rcx,%r8
  je     L87
L41:
  mov    %rcx,%rax
  cqto
  idiv   %rdi
  mov    %rax,%r14
  mov    %rdx,%rsi
  mov    0x8(%r11,%rax,8),%rax
  sub    %r9,%r14
  sub    %r10,%rsi
  mov    %r14,%r15
  movss  0x4(%rax,%rdx,4),%xmm1
  imul   %rsi,%r15
  mulss  %xmm3,%xmm1
  test   %r15,%r15
  jne    L17
  test   %r14,%r14
  je     L94
  addss  %xmm1,%xmm2
L94:
  test   %rsi,%rsi
  je     L57
  add    $0x1,%rcx
  addss  %xmm1,%xmm0
  cmp    %rcx,%r8
  jne    L41
L87:
  cmp    %r13,%rbx
  jne    L46
  mulss  %xmm2,%xmm0
L69:
  movss  %xmm5,0xc(%rsp)
  movss  %xmm4,0x8(%rsp)
  call   <sqrt@plt>
  movss  0x8(%rsp),%xmm4
  movss  0xc(%rsp),%xmm5
  movsd  D25(%rip),%xmm1
  divss  %xmm0,%xmm4
  cvtss2sd %xmm5,%xmm5
  movsd  D57(%rip),%xmm0
  mulsd  %xmm5,%xmm1
  mulsd  %xmm5,%xmm0
  subsd  D53(%rip),%xmm5
  addsd  D32(%rip),%xmm0
  mulsd  %xmm5,%xmm1
  divsd  %xmm1,%xmm0
  movss  %xmm4,0x0(%rbp)
  movss  %xmm4,0x8(%rsp)
  cvtsd2ss %xmm0,%xmm0
  call   <sqrt@plt>
  movss  0x8(%rsp),%xmm4
  divss  %xmm0,%xmm4
  movss  %xmm4,(%r12)
  movaps %xmm4,%xmm0
  call   <fabs@plt>
  divsd  D27(%rip),%xmm0
  cvtsd2ss %xmm0,%xmm0
  call   <erfcc@plt>
  mov    (%rsp),%rax
  movss  %xmm0,(%rax)
  add    $0x18,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
L36:
  pxor   %xmm0,%xmm0
  movaps %xmm0,%xmm4
  jmp    L69
```

tool_calls: name: disassemble_data
arguments: data_label: D34
data_type: byte[4]
name: disassemble_data
arguments: data_label: D25
data_type: f64
name: disassemble_data
arguments: data_label: D57
data_type: f64
name: disassemble_data
arguments: data_label: D53
data_type: f64
name: disassemble_data
arguments: data_label: D32
data_type: f64
name: disassemble_data
arguments: data_label: D27
data_type: f64

tool: {"D34": [0x00,0x00,0x00,0x80]}

tool: D25: 9

tool: D57: 4

tool: D53: 1

tool: D32: 10

tool: D27: 1.4142136

assistant: ```c
void kendl2(float **tab, int i, int j, float *tau, float *z, float *prob) {
  float erfcc(float x);
  long nn, mm, m2, m1, lj, li, l, kj, ki, k;
  float svar, s = 0.0, points, pairs, en2 = 0.0, en1 = 0.0;
  nn = i * j;
  points = tab[i][j];
  for (k = 0; k <= nn - 2; k++) {
    ki = (k / j);
    kj = k - j * ki;
    points += tab[ki + 1][kj + 1];
    for (l = k + 1; l <= nn - 1; l++) {
      li = l / j;
      lj = l - j * li;
      mm = (m1 = li - ki) * (m2 = lj - kj);
      pairs = tab[ki + 1][kj + 1] * tab[li + 1][lj + 1];
      if (mm) {
        en1 += pairs;
        en2 += pairs;
        s += (mm > 0 ? pairs : -pairs);
      } else {
        if (m1)
          en1 += pairs;
        if (m2)
          en2 += pairs;
      }
    }
  }
  *tau = s / sqrt(en1 * en2);
  svar = (4.0 * points + 10.0) / (9.0 * points * (points - 1.0));
  *z = (*tau) / sqrt(svar);
  *prob = erfcc(fabs(*z) / 1.4142136);
}
```