user: 这段汇编代码用C语言表示是什么样的:

```asm
<gridGenRotBounds>:
  endbr64
  push   %r15
  push   %r14
  push   %r13
  push   %r12
  push   %rbp
  push   %rbx
  sub    $0x68,%rsp
  mov    %edi,%ebp
  mov    %esi,%ebx
  mov    %edx,%r12d
  mov    %rcx,0x50(%rsp)
  mov    %r8,0x40(%rsp)
  mov    %r9,0x58(%rsp)
  call   <gridInqXpole@plt>
  movsd  %xmm0,0x28(%rsp)
  mov    %ebp,%edi
  call   <gridInqYpole@plt>
  movsd  %xmm0,0x8(%rsp)
  mov    %ebp,%edi
  call   <gridInqAngle@plt>
  movsd  %xmm0,0x10(%rsp)
  movslq %r12d,%rax
  mov    %rax,0x48(%rsp)
  test   %rax,%rax
  jle    L47
  movq   $0x0,0x38(%rsp)
  movslq %ebx,%r12
  jmp    L20
L55:
  mov    0x38(%rsp),%rax
  shl    $0x4,%rax
  mov    0x40(%rsp),%rsi
  movsd  0x8(%rsi,%rax,1),%xmm6
  movsd  %xmm6,0x20(%rsp)
  movsd  (%rsi,%rax,1),%xmm7
  movsd  %xmm7,0x18(%rsp)
L39:
  test   %r12,%r12
  jle    L29
  mov    0x50(%rsp),%r14
  mov    0x38(%rsp),%rbx
  imul   %r12,%rbx
  shl    $0x5,%rbx
  mov    0x58(%rsp),%rax
  lea    (%rax,%rbx,1),%rbp
  add    0xa0(%rsp),%rbx
  mov    $0x0,%r15d
L22:
  movsd  (%r14),%xmm5
  movsd  %xmm5,0x30(%rsp)
  mov    0x8(%r14),%r13
  movsd  0x10(%rsp),%xmm4
  movsd  0x28(%rsp),%xmm3
  movsd  0x8(%rsp),%xmm2
  movapd %xmm5,%xmm1
  movsd  0x18(%rsp),%xmm0
  call   <lamrot_to_lam@plt>
  movsd  %xmm0,0x0(%rbp)
  movsd  0x10(%rsp),%xmm4
  movsd  0x28(%rsp),%xmm3
  movsd  0x8(%rsp),%xmm2
  movq   %r13,%xmm1
  movsd  0x18(%rsp),%xmm0
  call   <lamrot_to_lam@plt>
  movsd  %xmm0,0x8(%rbp)
  movsd  0x10(%rsp),%xmm4
  movsd  0x28(%rsp),%xmm3
  movsd  0x8(%rsp),%xmm2
  movq   %r13,%xmm1
  movsd  0x20(%rsp),%xmm0
  call   <lamrot_to_lam@plt>
  movsd  %xmm0,0x10(%rbp)
  movsd  0x10(%rsp),%xmm4
  movsd  0x28(%rsp),%xmm3
  movsd  0x8(%rsp),%xmm2
  movsd  0x30(%rsp),%xmm1
  movsd  0x20(%rsp),%xmm0
  call   <lamrot_to_lam@plt>
  movsd  %xmm0,0x18(%rbp)
  movsd  0x10(%rsp),%xmm3
  movsd  0x8(%rsp),%xmm2
  movsd  0x30(%rsp),%xmm1
  movsd  0x18(%rsp),%xmm0
  call   <phirot_to_phi@plt>
  movsd  %xmm0,(%rbx)
  movsd  0x10(%rsp),%xmm3
  movsd  0x8(%rsp),%xmm2
  movq   %r13,%xmm1
  movsd  0x18(%rsp),%xmm0
  call   <phirot_to_phi@plt>
  movsd  %xmm0,0x8(%rbx)
  movsd  0x10(%rsp),%xmm3
  movsd  0x8(%rsp),%xmm2
  movq   %r13,%xmm1
  movsd  0x20(%rsp),%xmm0
  call   <phirot_to_phi@plt>
  movsd  %xmm0,0x10(%rbx)
  movsd  0x10(%rsp),%xmm3
  movsd  0x8(%rsp),%xmm2
  movsd  0x30(%rsp),%xmm1
  movsd  0x20(%rsp),%xmm0
  call   <phirot_to_phi@plt>
  movsd  %xmm0,0x18(%rbx)
  add    $0x1,%r15
  add    $0x10,%r14
  add    $0x20,%rbp
  add    $0x20,%rbx
  cmp    %r12,%r15
  jne    L22
L29:
  addq   $0x1,0x38(%rsp)
  mov    0x38(%rsp),%rax
  cmp    0x48(%rsp),%rax
  jge    L47
L20:
  mov    0x40(%rsp),%rdx
  movsd  (%rdx),%xmm0
  comisd 0x8(%rdx),%xmm0
  jbe    L55
  mov    0x38(%rsp),%rax
  shl    $0x4,%rax
  movsd  (%rdx,%rax,1),%xmm6
  movsd  %xmm6,0x20(%rsp)
  movsd  0x8(%rdx,%rax,1),%xmm7
  movsd  %xmm7,0x18(%rsp)
  jmp    L39
L47:
  add    $0x68,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
```

assistant: ```c
void gridGenRotBounds(int gridID, int nx, int ny, double *xbounds, double *ybounds, double *xbounds2D, double *ybounds2D) {
  long i, j, index;
  double minlon, maxlon;
  double minlat, maxlat;
  double xpole, ypole, angle;
  xpole = gridInqXpole(gridID);
  ypole = gridInqYpole(gridID);
  angle = gridInqAngle(gridID);
  for (j = 0; j < ny; j++) {
    if (ybounds[0] > ybounds[1]) {
      maxlat = ybounds[2 * j];
      minlat = ybounds[2 * j + 1];
    } else {
      maxlat = ybounds[2 * j + 1];
      minlat = ybounds[2 * j];
    }
    for (i = 0; i < nx; i++) {
      minlon = xbounds[2 * i];
      maxlon = xbounds[2 * i + 1];
      index = j * 4 * nx + 4 * i;
      xbounds2D[index + 0] = lamrot_to_lam(minlat, minlon, ypole, xpole, angle);
      xbounds2D[index + 1] = lamrot_to_lam(minlat, maxlon, ypole, xpole, angle);
      xbounds2D[index + 2] = lamrot_to_lam(maxlat, maxlon, ypole, xpole, angle);
      xbounds2D[index + 3] = lamrot_to_lam(maxlat, minlon, ypole, xpole, angle);
      ybounds2D[index + 0] = phirot_to_phi(minlat, minlon, ypole, angle);
      ybounds2D[index + 1] = phirot_to_phi(minlat, maxlon, ypole, angle);
      ybounds2D[index + 2] = phirot_to_phi(maxlat, maxlon, ypole, angle);
      ybounds2D[index + 3] = phirot_to_phi(maxlat, minlon, ypole, angle);
    }
  }
}
```