user: Demonstrate the C code that generates the assembly output below:

```asm
<sample_H>:
  endbr64
  push   %rbp
  mov    %rsp,%rbp
  push   %rbx
  sub    $0x48,%rsp
  mov    %edi,-0x24(%rbp)
  mov    %esi,-0x28(%rbp)
  mov    %rdx,-0x30(%rbp)
  mov    %rcx,-0x38(%rbp)
  mov    %r8,-0x40(%rbp)
  movsd  %xmm0,-0x48(%rbp)
  mov    %r9,-0x50(%rbp)
  movl   $0x0,-0x1c(%rbp)
  jmp    L76
L22:
  movl   $0x0,-0x18(%rbp)
  jmp    L41
L5:
  mov    -0x18(%rbp),%eax
  mov    %eax,-0x14(%rbp)
  jmp    L62
L21:
  mov    -0x14(%rbp),%eax
  cmp    -0x18(%rbp),%eax
  jne    L59
  mov    -0x40(%rbp),%rax
  add    $0x8,%rax
  movsd  (%rax),%xmm0
  mov    -0x40(%rbp),%rax
  mov    (%rax),%rax
  mov    -0x1c(%rbp),%edx
  movslq %edx,%rdx
  lea    0x0(,%rdx,8),%rcx
  mov    0x10(%rbp),%rdx
  add    %rcx,%rdx
  mov    (%rdx),%rcx
  mov    -0x18(%rbp),%edx
  movslq %edx,%rdx
  shl    $0x3,%rdx
  add    %rcx,%rdx
  mov    (%rdx),%rcx
  mov    -0x14(%rbp),%edx
  movslq %edx,%rdx
  shl    $0x3,%rdx
  lea    (%rcx,%rdx,1),%rbx
  movapd %xmm0,%xmm1
  movq   %rax,%xmm0
  call   <rgamma@plt>
  movq   %xmm0,%rax
  mov    %rax,(%rbx)
  mov    -0x1c(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x50(%rbp),%rax
  add    %rdx,%rax
  mov    (%rax),%rdx
  mov    -0x18(%rbp),%eax
  cltq
  shl    $0x3,%rax
  add    %rdx,%rax
  mov    (%rax),%rdx
  mov    -0x18(%rbp),%eax
  cltq
  shl    $0x2,%rax
  add    %rdx,%rax
  movl   $0xffffffff,(%rax)
  jmp    L32
L59:
  mov    -0x1c(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x50(%rbp),%rax
  add    %rdx,%rax
  mov    (%rax),%rdx
  mov    -0x18(%rbp),%eax
  cltq
  shl    $0x3,%rax
  add    %rdx,%rax
  mov    (%rax),%rdx
  mov    -0x14(%rbp),%eax
  cltq
  shl    $0x2,%rax
  lea    (%rdx,%rax,1),%rbx
  mov    -0x48(%rbp),%rax
  movq   %rax,%xmm0
  call   <rbernoulli@plt>
  mov    %eax,(%rbx)
  mov    -0x1c(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x50(%rbp),%rax
  add    %rdx,%rax
  mov    (%rax),%rdx
  mov    -0x18(%rbp),%eax
  cltq
  shl    $0x3,%rax
  add    %rdx,%rax
  mov    (%rax),%rdx
  mov    -0x14(%rbp),%eax
  cltq
  shl    $0x2,%rax
  add    %rdx,%rax
  mov    (%rax),%eax
  cmp    $0x1,%eax
  jne    L91
  mov    -0x38(%rbp),%rax
  add    $0x8,%rax
  movsd  (%rax),%xmm0
  mov    -0x38(%rbp),%rax
  mov    (%rax),%rax
  mov    -0x1c(%rbp),%edx
  movslq %edx,%rdx
  lea    0x0(,%rdx,8),%rcx
  mov    0x10(%rbp),%rdx
  add    %rcx,%rdx
  mov    (%rdx),%rcx
  mov    -0x18(%rbp),%edx
  movslq %edx,%rdx
  shl    $0x3,%rdx
  add    %rcx,%rdx
  mov    (%rdx),%rcx
  mov    -0x14(%rbp),%edx
  movslq %edx,%rdx
  shl    $0x3,%rdx
  lea    (%rcx,%rdx,1),%rbx
  movapd %xmm0,%xmm1
  movq   %rax,%xmm0
  call   <rgamma@plt>
  movq   %xmm0,%rax
  mov    %rax,(%rbx)
  jmp    L49
L91:
  mov    -0x30(%rbp),%rax
  add    $0x8,%rax
  movsd  (%rax),%xmm0
  mov    -0x30(%rbp),%rax
  mov    (%rax),%rax
  mov    -0x1c(%rbp),%edx
  movslq %edx,%rdx
  lea    0x0(,%rdx,8),%rcx
  mov    0x10(%rbp),%rdx
  add    %rcx,%rdx
  mov    (%rdx),%rcx
  mov    -0x18(%rbp),%edx
  movslq %edx,%rdx
  shl    $0x3,%rdx
  add    %rcx,%rdx
  mov    (%rdx),%rcx
  mov    -0x14(%rbp),%edx
  movslq %edx,%rdx
  shl    $0x3,%rdx
  lea    (%rcx,%rdx,1),%rbx
  movapd %xmm0,%xmm1
  movq   %rax,%xmm0
  call   <rgamma@plt>
  movq   %xmm0,%rax
  mov    %rax,(%rbx)
L49:
  mov    -0x1c(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x50(%rbp),%rax
  add    %rdx,%rax
  mov    (%rax),%rdx
  mov    -0x18(%rbp),%eax
  cltq
  shl    $0x3,%rax
  add    %rdx,%rax
  mov    (%rax),%rdx
  mov    -0x14(%rbp),%eax
  cltq
  shl    $0x2,%rax
  add    %rdx,%rax
  mov    -0x1c(%rbp),%edx
  movslq %edx,%rdx
  lea    0x0(,%rdx,8),%rcx
  mov    -0x50(%rbp),%rdx
  add    %rcx,%rdx
  mov    (%rdx),%rcx
  mov    -0x14(%rbp),%edx
  movslq %edx,%rdx
  shl    $0x3,%rdx
  add    %rcx,%rdx
  mov    (%rdx),%rcx
  mov    -0x18(%rbp),%edx
  movslq %edx,%rdx
  shl    $0x2,%rdx
  add    %rcx,%rdx
  mov    (%rax),%eax
  mov    %eax,(%rdx)
  mov    -0x1c(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    0x10(%rbp),%rax
  add    %rdx,%rax
  mov    (%rax),%rdx
  mov    -0x18(%rbp),%eax
  cltq
  shl    $0x3,%rax
  add    %rdx,%rax
  mov    (%rax),%rdx
  mov    -0x14(%rbp),%eax
  cltq
  shl    $0x3,%rax
  add    %rax,%rdx
  mov    -0x1c(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rcx
  mov    0x10(%rbp),%rax
  add    %rcx,%rax
  mov    (%rax),%rcx
  mov    -0x14(%rbp),%eax
  cltq
  shl    $0x3,%rax
  add    %rcx,%rax
  mov    (%rax),%rcx
  mov    -0x18(%rbp),%eax
  cltq
  shl    $0x3,%rax
  add    %rcx,%rax
  movsd  (%rdx),%xmm0
  movsd  %xmm0,(%rax)
L32:
  addl   $0x1,-0x14(%rbp)
L62:
  mov    -0x14(%rbp),%eax
  cmp    -0x28(%rbp),%eax
  jl     L21
  addl   $0x1,-0x18(%rbp)
L41:
  mov    -0x18(%rbp),%eax
  cmp    -0x28(%rbp),%eax
  jl     L5
  addl   $0x1,-0x1c(%rbp)
L76:
  mov    -0x1c(%rbp),%eax
  cmp    -0x24(%rbp),%eax
  jl     L22
  nop
  nop
  mov    -0x8(%rbp),%rbx
  leave
  ret
```

assistant: ```c
void sample_H(int D, int p, double *par_0, double *par_1, double *par_2, double par_p_t, int ***L, double ***Lambda) {
  for (int d = 0; d < D; d++) {
    for (int i = 0; i < p; i++) {
      for (int j = i; j < p; j++) {
        if (j == i) {
          Lambda[d][i][j] = rgamma(par_2[0], par_2[1]);
          L[d][i][i] = -1;
        } else {
          L[d][i][j] = rbernoulli(par_p_t);
          if (L[d][i][j] == 1) {
            Lambda[d][i][j] = rgamma(par_1[0], par_1[1]);
          } else {
            Lambda[d][i][j] = rgamma(par_0[0], par_0[1]);
          }
          L[d][j][i] = L[d][i][j];
          Lambda[d][j][i] = Lambda[d][i][j];
        }
      }
    }
  }
}
```