user: What would the C source look like for the assembly code provided here?

```asm
<l2projectVector>:
  endbr64
  push   %rbp
  mov    %rsp,%rbp
  push   %r15
  push   %r14
  mov    %edx,%r14d
  push   %r13
  mov    %esi,%r13d
  push   %r12
  mov    %r8,%r12
  push   %rbx
  sub    $0x48,%rsp
  mov    %edi,-0x50(%rbp)
  mov    %rcx,-0x70(%rbp)
  mov    %rsp,%rsi
  mov    %fs:0x28,%rax
  mov    %rax,-0x38(%rbp)
  xor    %eax,%eax
  movslq %edx,%rax
  lea    0xf(,%rax,8),%rdx
  mov    %rdx,%rax
  and    $0xfffffffffffff000,%rdx
  sub    %rdx,%rsi
  and    $0xfffffffffffffff0,%rax
  cmp    %rsi,%rsp
  je     L133
L76:
  sub    $0x1000,%rsp
  orq    $0x0,0xff8(%rsp)
  cmp    %rsi,%rsp
  jne    L76
L133:
  and    $0xfff,%eax
  sub    %rax,%rsp
  test   %rax,%rax
  jne    L95
L67:
  mov    -0x50(%rbp),%eax
  mov    %rsp,%rbx
  test   %eax,%eax
  jle    L12
  mov    %r14d,%eax
  lea    -0x1(%r14),%r15d
  movl   $0x0,-0x4c(%rbp)
  imul   %r13d,%eax
  movl   $0x0,-0x48(%rbp)
  movl   $0x0,-0x44(%rbp)
  mov    %eax,-0x5c(%rbp)
  lea    0x8(,%r15,8),%rax
  mov    %rax,-0x68(%rbp)
  mov    %r14d,%eax
  mov    %r13d,%r14d
  mov    %eax,%r13d
  nopl   0x0(%rax)
L137:
  test   %r13d,%r13d
  jle    L111
  mov    -0x68(%rbp),%rdx
  xor    %esi,%esi
  mov    %rbx,%rdi
  call   <memset@plt>
L111:
  test   %r14d,%r14d
  jle    L131
  movslq -0x4c(%rbp),%rdx
  mov    -0x70(%rbp),%rsi
  lea    -0x1(%r14),%eax
  pxor   %xmm1,%xmm1
  mov    -0x48(%rbp),%r10d
  add    %rdx,%rax
  lea    (%rsi,%rdx,8),%rdi
  lea    0x8(%rsi,%rax,8),%r11
  nopl   0x0(%rax)
L156:
  movsd  (%rdi),%xmm2
  addsd  %xmm2,%xmm1
  test   %r13d,%r13d
  jle    L6
  movslq %r10d,%rax
  lea    (%r12,%rax,8),%rcx
  xor    %eax,%eax
  xchg   %ax,%ax
L159:
  movsd  (%rcx,%rax,8),%xmm0
  mov    %rax,%rdx
  mulsd  %xmm2,%xmm0
  addsd  (%rbx,%rax,8),%xmm0
  movsd  %xmm0,(%rbx,%rax,8)
  add    $0x1,%rax
  cmp    %rdx,%r15
  jne    L159
L6:
  add    $0x8,%rdi
  add    %r13d,%r10d
  cmp    %r11,%rdi
  jne    L156
  pxor   %xmm3,%xmm3
  movsd  %xmm1,-0x58(%rbp)
  comisd %xmm3,%xmm1
  seta   %dil
  movzbl %dil,%edi
  call   <assert@plt>
  mov    -0x48(%rbp),%r10d
  movsd  -0x58(%rbp),%xmm1
  xor    %edi,%edi
L138:
  test   %r13d,%r13d
  jle    L9
  movslq %r10d,%rax
  lea    (%r12,%rax,8),%rcx
  xor    %eax,%eax
  xchg   %ax,%ax
L84:
  movsd  (%rbx,%rax,8),%xmm0
  mov    %rax,%rdx
  divsd  %xmm1,%xmm0
  movsd  %xmm0,(%rcx,%rax,8)
  add    $0x1,%rax
  cmp    %r15,%rdx
  jne    L84
L9:
  add    $0x1,%edi
  add    %r13d,%r10d
  cmp    %edi,%r14d
  jne    L138
L46:
  addl   $0x1,-0x44(%rbp)
  mov    -0x5c(%rbp),%edi
  add    %r14d,-0x4c(%rbp)
  mov    -0x44(%rbp),%eax
  add    %edi,-0x48(%rbp)
  cmp    %eax,-0x50(%rbp)
  jne    L137
L12:
  mov    -0x38(%rbp),%rax
  sub    %fs:0x28,%rax
  jne    L136
  lea    -0x28(%rbp),%rsp
  pop    %rbx
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  pop    %rbp
  ret
L131:
  xor    %edi,%edi
  call   <assert@plt>
  jmp    L46
L95:
  orq    $0x0,-0x8(%rsp,%rax,1)
  jmp    L67
L136:
  call   <__stack_chk_fail@plt>
```

assistant: ```c
void l2projectVector(const int nSimplices, const int nPointsPerSimplex, const int nSpace, double *dV, double *r) {
  int eN, k, I;
  double vol;
  double ravg[nSpace];
  for (eN = 0; eN < nSimplices; eN++) {
    vol = 0.0;
    for (I = 0; I < nSpace; I++)
      ravg[I] = 0.0;
    for (k = 0; k < nPointsPerSimplex; k++) {
      vol += dV[eN * nPointsPerSimplex + k];
      for (I = 0; I < nSpace; I++)
        ravg[I] += r[eN * nPointsPerSimplex * nSpace + k * nSpace + I] * dV[eN * nPointsPerSimplex + k];
    }
    assert(vol > 0.0);
    for (k = 0; k < nPointsPerSimplex; k++) {
      for (I = 0; I < nSpace; I++)
        r[eN * nPointsPerSimplex * nSpace + k * nSpace + I] = ravg[I] / vol;
    }
  }
}
```