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

```asm
<l2norm>:
  endbr64
  push   %r15
  push   %r14
  push   %r13
  push   %r12
  push   %rbp
  push   %rbx
  sub    $0x18,%rsp
  mov    %esi,%r10d
  mov    %edx,%r11d
  mov    %ecx,%edx
  mov    %r9d,%ebp
  mov    0x50(%rsp),%r12d
  mov    0x58(%rsp),%rcx
  mov    0x60(%rsp),%r9
  movq   $0x0,(%r9)
  movq   $0x0,0x8(%r9)
  movq   $0x0,0x10(%r9)
  movq   $0x0,0x18(%r9)
  movq   $0x0,0x20(%r9)
  cmp    %r8d,%edx
  jg     L91
  mov    %r8d,%eax
  movslq %edx,%rsi
  imul   $0x29428,%rsi,%r8
  add    %rcx,%r8
  sub    %edx,%eax
  add    %rsi,%rax
  imul   $0x29428,%rax,%rax
  lea    0x29428(%rcx,%rax,1),%r13
  movslq %ebp,%rdx
  imul   $0xa28,%rdx,%rax
  lea    0x28(%rax),%r15
  lea    -0x3(%r11),%ecx
  lea    (%rcx,%rcx,4),%rcx
  mov    %r12d,%ebx
  sub    %ebp,%ebx
  add    %rdx,%rbx
  imul   $0xa28,%rbx,%rbx
  pxor   %xmm1,%xmm1
  movapd %xmm1,%xmm2
  movapd %xmm1,%xmm3
  movapd %xmm1,%xmm4
  movapd %xmm1,%xmm5
  lea    0x50(%rax,%rcx,8),%r14
  add    $0xa50,%rbx
  jmp    L41
L92:
  movsd  (%rax),%xmm0
  mulsd  %xmm0,%xmm0
  addsd  %xmm0,%xmm5
  movsd  0x8(%rax),%xmm0
  mulsd  %xmm0,%xmm0
  addsd  %xmm0,%xmm4
  movsd  0x10(%rax),%xmm0
  mulsd  %xmm0,%xmm0
  addsd  %xmm0,%xmm3
  movsd  0x18(%rax),%xmm0
  mulsd  %xmm0,%xmm0
  addsd  %xmm0,%xmm2
  movsd  0x20(%rax),%xmm0
  mulsd  %xmm0,%xmm0
  addsd  %xmm0,%xmm1
  add    $0x28,%rax
  cmp    %rdx,%rax
  jne    L92
L25:
  add    $0xa28,%rcx
  add    $0xa28,%rdx
  cmp    %rsi,%rcx
  je     L101
L50:
  mov    %rcx,%rax
  cmp    $0x2,%r11d
  jg     L92
  jmp    L25
L101:
  add    $0x29428,%r8
  cmp    %r13,%r8
  je     L11
L41:
  cmp    %r12d,%ebp
  jg     L101
  lea    (%r15,%r8,1),%rcx
  lea    (%r14,%r8,1),%rdx
  lea    (%rbx,%r8,1),%rsi
  jmp    L50
L91:
  pxor   %xmm1,%xmm1
  movapd %xmm1,%xmm2
  movapd %xmm1,%xmm3
  movapd %xmm1,%xmm4
  movapd %xmm1,%xmm5
L11:
  addsd  (%r9),%xmm5
  movsd  %xmm5,(%r9)
  addsd  0x8(%r9),%xmm4
  movsd  %xmm4,0x8(%r9)
  addsd  0x10(%r9),%xmm3
  movsd  %xmm3,0x10(%r9)
  addsd  0x18(%r9),%xmm2
  movsd  %xmm2,0x18(%r9)
  addsd  0x20(%r9),%xmm1
  movsd  %xmm1,0x20(%r9)
  lea    -0x2(%rdi),%eax
  sub    $0x2,%r10d
  imul   %r10d,%eax
  sub    $0x2,%r11d
  imul   %r11d,%eax
  pxor   %xmm6,%xmm6
  cvtsi2sd %eax,%xmm6
  movsd  %xmm6,0x8(%rsp)
  mov    %r9,%rbx
  lea    0x28(%r9),%rbp
L44:
  mov    %rbx,%r12
  movsd  (%rbx),%xmm0
  divsd  0x8(%rsp),%xmm0
  pxor   %xmm7,%xmm7
  ucomisd %xmm0,%xmm7
  ja     L10
  sqrtsd %xmm0,%xmm0
L59:
  movsd  %xmm0,(%r12)
  add    $0x8,%rbx
  cmp    %rbx,%rbp
  jne    L44
  add    $0x18,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
L10:
  call   <sqrt@plt>
  jmp    L59
```

assistant: ```c
void l2norm(int nx0, int ny0, int nz0, int ist, int iend, int jst, int jend, double v[64][64 / 2 * 2 + 1][64 / 2 * 2 + 1][5], double sum[5]) {
  int i, j, k, m;
  double sum0 = 0.0, sum1 = 0.0, sum2 = 0.0, sum3 = 0.0, sum4 = 0.0;
#pragma omp single
  for (m = 0; m < 5; m++) {
    sum[m] = 0.0;
  }
#pragma omp for nowait
  for (i = ist; i <= iend; i++) {
    for (j = jst; j <= jend; j++) {
      for (k = 1; k <= nz0 - 2; k++) {
        sum0 = sum0 + v[i][j][k][0] * v[i][j][k][0];
        sum1 = sum1 + v[i][j][k][1] * v[i][j][k][1];
        sum2 = sum2 + v[i][j][k][2] * v[i][j][k][2];
        sum3 = sum3 + v[i][j][k][3] * v[i][j][k][3];
        sum4 = sum4 + v[i][j][k][4] * v[i][j][k][4];
      }
    }
  }
#pragma omp critical
  {
    sum[0] += sum0;
    sum[1] += sum1;
    sum[2] += sum2;
    sum[3] += sum3;
    sum[4] += sum4;
  }
#pragma omp barrier
#pragma omp single
  for (m = 0; m < 5; m++) {
    sum[m] = sqrt(sum[m] / ((nx0 - 2) * (ny0 - 2) * (nz0 - 2)));
  }
}
```