user: 请将这段汇编代码片段转换成等效的C代码:

```asm
<bbob2009_compute_rotation>:
  endbr64
  push   %rbp
  mov    %rsp,%rbp
  sub    $0x1000,%rsp
  orq    $0x0,(%rsp)
  sub    $0x1000,%rsp
  orq    $0x0,(%rsp)
  sub    $0x1000,%rsp
  orq    $0x0,(%rsp)
  sub    $0xed0,%rsp
  mov    %rdi,-0x3eb8(%rbp)
  mov    %rsi,-0x3ec0(%rbp)
  mov    %rdx,-0x3ec8(%rbp)
  mov    %fs:0x28,%rax
  mov    %rax,-0x8(%rbp)
  xor    %eax,%eax
  mov    -0x3ec8(%rbp),%rax
  imul   %rax,%rax
  cmp    $0x7cf,%rax
  setbe  %al
  movzbl %al,%eax
  mov    %eax,%edi
  call   <assert@plt>
  mov    -0x3ec8(%rbp),%rax
  imul   %rax,%rax
  mov    %rax,%rcx
  mov    -0x3ec0(%rbp),%rdx
  lea    -0x3e90(%rbp),%rax
  mov    %rcx,%rsi
  mov    %rax,%rdi
  call   <bbob2009_gauss@plt>
  mov    -0x3ec8(%rbp),%rcx
  mov    -0x3ec8(%rbp),%rdx
  lea    -0x3e90(%rbp),%rsi
  mov    -0x3eb8(%rbp),%rax
  mov    %rax,%rdi
  call   <bbob2009_reshape@plt>
  movq   $0x0,-0x3ea8(%rbp)
  jmp    L75
L13:
  movq   $0x0,-0x3ea0(%rbp)
  jmp    L109
L112:
  pxor   %xmm0,%xmm0
  movsd  %xmm0,-0x3eb0(%rbp)
  movq   $0x0,-0x3e98(%rbp)
  jmp    L97
L2:
  mov    -0x3e98(%rbp),%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0x3eb8(%rbp),%rax
  add    %rdx,%rax
  mov    (%rax),%rdx
  mov    -0x3ea8(%rbp),%rax
  shl    $0x3,%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm1
  mov    -0x3e98(%rbp),%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0x3eb8(%rbp),%rax
  add    %rdx,%rax
  mov    (%rax),%rdx
  mov    -0x3ea0(%rbp),%rax
  shl    $0x3,%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  mulsd  %xmm1,%xmm0
  movsd  -0x3eb0(%rbp),%xmm1
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x3eb0(%rbp)
  addq   $0x1,-0x3e98(%rbp)
L97:
  mov    -0x3e98(%rbp),%rax
  cmp    %rax,-0x3ec8(%rbp)
  ja     L2
  movq   $0x0,-0x3e98(%rbp)
  jmp    L57
L38:
  mov    -0x3e98(%rbp),%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0x3eb8(%rbp),%rax
  add    %rdx,%rax
  mov    (%rax),%rdx
  mov    -0x3ea8(%rbp),%rax
  shl    $0x3,%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  mov    -0x3e98(%rbp),%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0x3eb8(%rbp),%rax
  add    %rdx,%rax
  mov    (%rax),%rdx
  mov    -0x3ea0(%rbp),%rax
  shl    $0x3,%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm1
  mulsd  -0x3eb0(%rbp),%xmm1
  mov    -0x3e98(%rbp),%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0x3eb8(%rbp),%rax
  add    %rdx,%rax
  mov    (%rax),%rdx
  mov    -0x3ea8(%rbp),%rax
  shl    $0x3,%rax
  add    %rdx,%rax
  subsd  %xmm1,%xmm0
  movsd  %xmm0,(%rax)
  addq   $0x1,-0x3e98(%rbp)
L57:
  mov    -0x3e98(%rbp),%rax
  cmp    %rax,-0x3ec8(%rbp)
  ja     L38
  addq   $0x1,-0x3ea0(%rbp)
L109:
  mov    -0x3ea0(%rbp),%rax
  cmp    -0x3ea8(%rbp),%rax
  jl     L112
  pxor   %xmm0,%xmm0
  movsd  %xmm0,-0x3eb0(%rbp)
  movq   $0x0,-0x3e98(%rbp)
  jmp    L35
L58:
  mov    -0x3e98(%rbp),%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0x3eb8(%rbp),%rax
  add    %rdx,%rax
  mov    (%rax),%rdx
  mov    -0x3ea8(%rbp),%rax
  shl    $0x3,%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm1
  mov    -0x3e98(%rbp),%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0x3eb8(%rbp),%rax
  add    %rdx,%rax
  mov    (%rax),%rdx
  mov    -0x3ea8(%rbp),%rax
  shl    $0x3,%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  mulsd  %xmm1,%xmm0
  movsd  -0x3eb0(%rbp),%xmm1
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x3eb0(%rbp)
  addq   $0x1,-0x3e98(%rbp)
L35:
  mov    -0x3e98(%rbp),%rax
  cmp    %rax,-0x3ec8(%rbp)
  ja     L58
  movq   $0x0,-0x3e98(%rbp)
  jmp    L5
L105:
  mov    -0x3eb0(%rbp),%rax
  movq   %rax,%xmm0
  call   <sqrt@plt>
  movapd %xmm0,%xmm1
  mov    -0x3e98(%rbp),%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0x3eb8(%rbp),%rax
  add    %rdx,%rax
  mov    (%rax),%rdx
  mov    -0x3ea8(%rbp),%rax
  shl    $0x3,%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm0
  mov    -0x3e98(%rbp),%rax
  lea    0x0(,%rax,8),%rdx
  mov    -0x3eb8(%rbp),%rax
  add    %rdx,%rax
  mov    (%rax),%rdx
  mov    -0x3ea8(%rbp),%rax
  shl    $0x3,%rax
  add    %rdx,%rax
  divsd  %xmm1,%xmm0
  movsd  %xmm0,(%rax)
  addq   $0x1,-0x3e98(%rbp)
L5:
  mov    -0x3e98(%rbp),%rax
  cmp    %rax,-0x3ec8(%rbp)
  ja     L105
  addq   $0x1,-0x3ea8(%rbp)
L75:
  mov    -0x3ea8(%rbp),%rax
  cmp    %rax,-0x3ec8(%rbp)
  ja     L13
  nop
  mov    -0x8(%rbp),%rax
  sub    %fs:0x28,%rax
  je     L94
  call   <__stack_chk_fail@plt>
L94:
  leave
  ret
```

assistant: ```c
void bbob2009_compute_rotation(double **B, long seed, size_t DIM) {
  double prod;
  double gvect[2000];
  long i, j, k;
  assert(DIM * DIM < 2000);
  bbob2009_gauss(gvect, DIM * DIM, seed);
  bbob2009_reshape(B, gvect, DIM, DIM);
  for (i = 0; i < DIM; i++) {
    for (j = 0; j < i; j++) {
      prod = 0;
      for (k = 0; k < DIM; k++)
        prod += B[k][i] * B[k][j];
      for (k = 0; k < DIM; k++)
        B[k][i] -= prod * B[k][j];
    }
    prod = 0;
    for (k = 0; k < DIM; k++)
      prod += B[k][i] * B[k][i];
    for (k = 0; k < DIM; k++)
      B[k][i] /= sqrt(prod);
  }
}
```