user: Converte this assembly snippet into its equivalent C code:

```asm
<compute>:
  endbr64
  push   %rbp
  mov    %rsp,%rbp
  mov    %rdi,-0x18(%rbp)
  mov    %rsi,-0x20(%rbp)
  mov    %rdx,-0x28(%rbp)
  mov    %rcx,-0x30(%rbp)
  mov    %r8d,-0x34(%rbp)
  movl   $0x1,-0x8(%rbp)
  jmp    L21
L31:
  mov    -0x8(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x20(%rbp),%rax
  add    %rdx,%rax
  mov    (%rax),%rax
  add    $0x8,%rax
  mov    (%rax),%rdx
  mov    -0x8(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rcx
  mov    -0x18(%rbp),%rax
  add    %rcx,%rax
  mov    (%rax),%rax
  lea    0x8(%rax),%rcx
  mov    %rdx,%rax
  add    %rax,%rax
  add    %rdx,%rax
  mov    %rax,(%rcx)
  mov    -0x8(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x18(%rbp),%rax
  add    %rdx,%rax
  mov    (%rax),%rax
  add    $0x8,%rax
  mov    (%rax),%rcx
  mov    -0x8(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x30(%rbp),%rax
  add    %rdx,%rax
  mov    (%rax),%rax
  add    $0x8,%rax
  mov    (%rax),%rdx
  mov    -0x8(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rsi
  mov    -0x28(%rbp),%rax
  add    %rsi,%rax
  mov    (%rax),%rax
  imul   %rcx,%rdx
  mov    %rdx,(%rax)
  movl   $0x2,-0x4(%rbp)
  jmp    L28
L6:
  mov    -0x8(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x28(%rbp),%rax
  add    %rdx,%rax
  mov    (%rax),%rdx
  mov    -0x4(%rbp),%eax
  cltq
  shl    $0x3,%rax
  sub    $0x8,%rax
  add    %rdx,%rax
  mov    (%rax),%rcx
  mov    -0x8(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x28(%rbp),%rax
  add    %rdx,%rax
  mov    (%rax),%rdx
  mov    -0x4(%rbp),%eax
  cltq
  shl    $0x3,%rax
  sub    $0x8,%rax
  add    %rdx,%rax
  mov    (%rax),%rax
  mov    -0x8(%rbp),%edx
  movslq %edx,%rdx
  lea    0x0(,%rdx,8),%rsi
  mov    -0x20(%rbp),%rdx
  add    %rsi,%rdx
  mov    (%rdx),%rsi
  mov    -0x4(%rbp),%edx
  movslq %edx,%rdx
  shl    $0x3,%rdx
  add    %rsi,%rdx
  imul   %rcx,%rax
  mov    %rax,(%rdx)
  mov    -0x8(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x20(%rbp),%rax
  add    %rdx,%rax
  mov    (%rax),%rdx
  mov    -0x4(%rbp),%eax
  cltq
  shl    $0x3,%rax
  add    %rdx,%rax
  mov    (%rax),%rdx
  mov    -0x8(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rcx
  mov    -0x18(%rbp),%rax
  add    %rcx,%rax
  mov    (%rax),%rcx
  mov    -0x4(%rbp),%eax
  cltq
  shl    $0x3,%rax
  add    %rax,%rcx
  mov    %rdx,%rax
  add    %rax,%rax
  add    %rdx,%rax
  mov    %rax,(%rcx)
  mov    -0x8(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x18(%rbp),%rax
  add    %rdx,%rax
  mov    (%rax),%rdx
  mov    -0x4(%rbp),%eax
  cltq
  shl    $0x3,%rax
  add    %rdx,%rax
  mov    (%rax),%rcx
  mov    -0x8(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x30(%rbp),%rax
  add    %rdx,%rax
  mov    (%rax),%rdx
  mov    -0x4(%rbp),%eax
  cltq
  shl    $0x3,%rax
  add    %rdx,%rax
  mov    (%rax),%rax
  mov    -0x8(%rbp),%edx
  movslq %edx,%rdx
  lea    0x0(,%rdx,8),%rsi
  mov    -0x28(%rbp),%rdx
  add    %rsi,%rdx
  mov    (%rdx),%rsi
  mov    -0x4(%rbp),%edx
  movslq %edx,%rdx
  shl    $0x3,%rdx
  sub    $0x8,%rdx
  add    %rsi,%rdx
  imul   %rcx,%rax
  mov    %rax,(%rdx)
  addl   $0x1,-0x4(%rbp)
L28:
  mov    -0x4(%rbp),%eax
  cmp    -0x34(%rbp),%eax
  jl     L6
  mov    -0x8(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x28(%rbp),%rax
  add    %rdx,%rax
  mov    (%rax),%rdx
  mov    -0x34(%rbp),%eax
  cltq
  shl    $0x3,%rax
  sub    $0x8,%rax
  add    %rdx,%rax
  mov    (%rax),%rcx
  mov    -0x8(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x28(%rbp),%rax
  add    %rdx,%rax
  mov    (%rax),%rdx
  mov    -0x34(%rbp),%eax
  cltq
  shl    $0x3,%rax
  sub    $0x8,%rax
  add    %rdx,%rax
  mov    (%rax),%rax
  mov    -0x8(%rbp),%edx
  movslq %edx,%rdx
  lea    0x0(,%rdx,8),%rsi
  mov    -0x20(%rbp),%rdx
  add    %rsi,%rdx
  mov    (%rdx),%rsi
  mov    -0x34(%rbp),%edx
  movslq %edx,%rdx
  shl    $0x3,%rdx
  add    %rsi,%rdx
  imul   %rcx,%rax
  mov    %rax,(%rdx)
  addl   $0x1,-0x8(%rbp)
L21:
  mov    -0x8(%rbp),%eax
  cmp    -0x34(%rbp),%eax
  jl     L31
  nop
  nop
  pop    %rbp
  ret
```

assistant: ```c
void compute(unsigned long **a, unsigned long **b, unsigned long **c, unsigned long **d, int N) {
  int i, j;
#pragma omp parallel for
  for (i = 1; i < N; i++) {
    a[i][1] = 3 * b[i][1];
    c[i][0] = a[i][1] * d[i][1];
#pragma omp parallel for
    for (j = 2; j < N; j++) {
      b[i][j] = c[i][j - 1] * c[i][j - 1];
      a[i][j] = 3 * b[i][j];
      c[i][j - 1] = a[i][j] * d[i][j];
    }
    b[i][N] = c[i][N - 1] * c[i][N - 1];
  }
}
```