user: Translate the below assembly code into C programming language:

```asm
<initialize_s>:
  endbr64
  push   %r15
  push   %r14
  push   %r13
  push   %r12
  push   %rbp
  push   %rbx
  sub    $0x38,%rsp
  mov    %esi,0x2c(%rsp)
  mov    %rcx,0x18(%rsp)
  test   %esi,%esi
  jle    L21
  mov    %rdi,%r8
  mov    %rdx,%r11
  mov    0x18(%rsp),%rbp
  mov    %rdi,%r15
  mov    %rdx,0x20(%rsp)
  mov    0x2c(%rsp),%eax
  lea    -0x1(%rax),%eax
  lea    0x8(%rdi,%rax,8),%r12
  jmp    L50
L131:
  mov    (%r8),%rcx
  mov    (%rcx,%rdi,1),%rcx
  movl   $0x0,(%rcx,%rdx,1)
  add    $0x1,%eax
  mov    (%rbx),%rcx
  mov    (%rcx,%rsi,1),%ecx
  add    $0x4,%rdx
  cmp    %eax,%ecx
  jge    L131
L41:
  mov    (%r8),%rax
  add    $0x1,%ecx
  movslq %ecx,%rcx
  mov    (%rax,%rdi,1),%rax
  movl   $0x0,-0x4(%rax,%rcx,4)
  add    $0x1,%r9d
  add    $0x4,%rsi
  add    $0x8,%rdi
  cmp    %r9d,(%r10)
  jl     L112
L23:
  mov    %r11,%rbx
  mov    (%r11),%rax
  mov    (%rax,%rsi,1),%ecx
  mov    $0x0,%edx
  mov    $0x1,%eax
  test   %ecx,%ecx
  jg     L131
  jmp    L41
L112:
  add    $0x4,%rbp
  add    $0x8,%r8
  add    $0x8,%r11
  cmp    %r12,%r8
  je     L6
L50:
  mov    %rbp,%r10
  mov    $0x0,%edi
  mov    $0x0,%esi
  mov    $0x1,%r9d
  cmpl   $0x0,0x0(%rbp)
  jg     L23
  jmp    L112
L6:
  movq   $0x1,0x10(%rsp)
  jmp    L169
L88:
  mov    (%r15),%rax
  mov    (%rax,%r12,1),%rax
  lea    (%rax,%rbx,4),%rbp
  mov    $0x2,%esi
  mov    $0x1,%edi
  mov    $0x0,%eax
  call   <sample_MultSym@plt>
  mov    %eax,0x0(%rbp)
  add    $0x1,%rbx
  mov    (%r14),%rdx
  lea    0x1(%rbx),%eax
  cmp    %eax,(%rdx,%r13,1)
  jge    L88
L72:
  addq   $0x1,0x8(%rsp)
  mov    0x8(%rsp),%rax
  add    $0x1,%eax
  mov    0x18(%rsp),%rdi
  mov    0x10(%rsp),%rsi
  cmp    -0x4(%rdi,%rsi,4),%eax
  jg     L31
L38:
  mov    0x20(%rsp),%rax
  mov    %rax,%r14
  mov    0x8(%rsp),%rsi
  lea    0x0(,%rsi,4),%r13
  mov    (%rax),%rax
  lea    0x0(,%rsi,8),%r12
  mov    $0x0,%ebx
  cmpl   $0x0,(%rax,%rsi,4)
  jg     L88
  jmp    L72
L31:
  addq   $0x1,0x10(%rsp)
  mov    0x10(%rsp),%rax
  addq   $0x8,0x20(%rsp)
  add    $0x8,%r15
  cmp    %eax,0x2c(%rsp)
  jl     L21
L169:
  movq   $0x0,0x8(%rsp)
  mov    0x18(%rsp),%rax
  mov    0x10(%rsp),%rdi
  cmpl   $0x0,-0x4(%rax,%rdi,4)
  jg     L38
  jmp    L31
L21:
  add    $0x38,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
```

assistant: ```c
void initialize_s(int ***s, int D, int **M, int *N) {
  int d_10;
  int n_16;
  int m_22;
  int dvv_loop_var_1;
  int dvv_loop_var_2;
  int dvv_loop_var_3;
  for (dvv_loop_var_1 = 1; dvv_loop_var_1 <= D; dvv_loop_var_1++) {
    for (dvv_loop_var_2 = 1; dvv_loop_var_2 <= N[dvv_loop_var_1 - 1]; dvv_loop_var_2++) {
      for (dvv_loop_var_3 = 1; dvv_loop_var_3 <= M[dvv_loop_var_1 - 1][dvv_loop_var_2 - 1]; dvv_loop_var_3++) {
        s[dvv_loop_var_1 - 1][dvv_loop_var_2 - 1][dvv_loop_var_3 - 1] = 0;
      }
      s[dvv_loop_var_1 - 1][dvv_loop_var_2 - 1][(M[dvv_loop_var_1 - 1][dvv_loop_var_2 - 1]) + (1) - 1] = (0) * (((1) + (M[dvv_loop_var_1 - 1][dvv_loop_var_2 - 1])) - (1));
    }
  }
  for (d_10 = 1; d_10 <= D; d_10++) {
    for (n_16 = 1; n_16 <= N[d_10 - 1]; n_16++) {
      for (m_22 = 1; m_22 <= M[d_10 - 1][n_16 - 1]; m_22++) {
        s[d_10 - 1][n_16 - 1][m_22 - 1] = sample_MultSym(1, 2);
      }
    }
  }
}
```