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

```asm
<bitrv2>:
  endbr64
  push   %r15
  mov    %edi,%r8d
  sar    %edi
  push   %r14
  sar    $0x2,%r8d
  push   %r13
  push   %r12
  push   %rbp
  mov    %rdx,%rbp
  push   %rbx
  sub    $0x70,%rsp
  mov    %edi,-0x48(%rsp)
  cmp    $0x8,%r8d
  jle    L105
  mov    $0x1,%eax
  nopl   0x0(%rax)
L35:
  sar    $0x2,%r8d
  add    %eax,%eax
  cmp    $0x8,%r8d
  jg     L35
  mov    %eax,(%rsp)
  lea    0x0(,%rax,4),%r14d
  je     L28
  test   %eax,%eax
  jle    L111
L121:
  movslq (%rsp),%rax
  xor    %r15d,%r15d
  xor    %ebx,%ebx
  movslq %r14d,%rdi
  mov    %r15d,-0x78(%rsp)
  lea    (%rsi,%rax,4),%r12
  movslq -0x48(%rsp),%rsi
  mov    $0x2,%eax
  sub    %rsi,%rax
  mov    %rsi,-0x70(%rsp)
  mov    %rsi,%r10
  mov    %rax,-0x68(%rsp)
  lea    0x20(%rbp),%rax
  mov    %rax,-0x60(%rsp)
  nopl   (%rax)
L24:
  mov    (%r12,%rbx,1),%eax
  addl   $0x1,-0x78(%rsp)
  add    %ebx,%eax
  lea    0x2(%rax),%edx
  add    %r10d,%eax
  movslq %edx,%rcx
  add    %r14d,%edx
  lea    0x0(%rbp,%rcx,8),%rsi
  movslq %eax,%rcx
  add    %r14d,%eax
  movslq %edx,%rdx
  lea    0x0(%rbp,%rcx,8),%rcx
  movupd (%rsi),%xmm0
  cltq
  lea    0x0(%rbp,%rdx,8),%rdx
  movupd (%rcx),%xmm5
  lea    0x0(%rbp,%rax,8),%rax
  movups %xmm5,(%rsi)
  movups %xmm0,(%rcx)
  movupd (%rax),%xmm6
  movupd (%rdx),%xmm0
  movups %xmm6,(%rdx)
  movups %xmm0,(%rax)
  mov    -0x78(%rsp),%eax
  cmp    (%rsp),%eax
  jge    L111
  mov    0x4(%r12,%rbx,1),%r8d
  lea    0x0(,%rax,4),%r13d
  mov    -0x60(%rsp),%r11
  mov    %r12,%r9
  movslq %r8d,%rax
  add    %r10d,%r8d
  lea    (%rdi,%rax,1),%rdx
  lea    0x0(%rbp,%rax,8),%rsi
  movslq %r8d,%r8
  add    -0x68(%rsp),%r8
  lea    0x0(%rbp,%rdx,8),%rcx
  mov    -0x70(%rsp),%rdx
  lea    0x0(%rbp,%r8,8),%r8
  add    %rax,%rdx
  add    %rbx,%rax
  lea    0x0(%rbp,%rdx,8),%rdx
  lea    (%r11,%rax,8),%r11
  nopl   0x0(%rax)
L67:
  mov    (%r9),%eax
  movupd (%rsi),%xmm0
  add    $0x20,%rsi
  add    $0x4,%r9
  add    $0x20,%rcx
  add    $0x20,%r8
  add    %r13d,%eax
  movslq %eax,%r15
  add    %r14d,%eax
  lea    0x0(%rbp,%r15,8),%r15
  movupd (%r15),%xmm1
  movups %xmm1,-0x20(%rsi)
  movups %xmm0,(%r15)
  movslq %eax,%r15
  movupd -0x20(%rcx),%xmm0
  add    $0x2,%eax
  lea    0x0(%rbp,%r15,8),%r15
  movupd (%r15),%xmm2
  movups %xmm2,-0x20(%rcx)
  movups %xmm0,(%r15)
  movslq %eax,%r15
  movupd (%rdx,%rdi,8),%xmm0
  sub    %r14d,%eax
  lea    0x0(%rbp,%r15,8),%r15
  movupd (%r15),%xmm3
  movups %xmm3,(%rdx,%rdi,8)
  movups %xmm0,(%r15)
  movslq %eax,%r15
  movupd (%rdx),%xmm0
  add    %r10d,%eax
  lea    0x0(%rbp,%r15,8),%r15
  movupd (%r15),%xmm4
  movups %xmm4,(%rdx)
  movups %xmm0,(%r15)
  movslq %eax,%r15
  movupd 0x10(%rdx),%xmm0
  add    %r14d,%eax
  lea    0x0(%rbp,%r15,8),%r15
  movupd (%r15),%xmm5
  movups %xmm5,0x10(%rdx)
  movups %xmm0,(%r15)
  movslq %eax,%r15
  movupd 0x10(%rdx,%rdi,8),%xmm0
  sub    $0x2,%eax
  lea    0x0(%rbp,%r15,8),%r15
  movupd (%r15),%xmm6
  movups %xmm6,0x10(%rdx,%rdi,8)
  add    $0x20,%rdx
  movups %xmm0,(%r15)
  movslq %eax,%r15
  sub    %r14d,%eax
  movupd -0x10(%rcx),%xmm0
  lea    0x0(%rbp,%r15,8),%r15
  cltq
  movupd (%r15),%xmm7
  lea    0x0(%rbp,%rax,8),%rax
  movups %xmm7,-0x10(%rcx)
  movups %xmm0,(%r15)
  movupd -0x20(%r8),%xmm0
  movupd (%rax),%xmm1
  movups %xmm1,-0x20(%r8)
  movups %xmm0,(%rax)
  cmp    %rsi,%r11
  jne    L67
  add    $0x4,%rbx
  jmp    L24
  nopw   0x0(%rax,%rax,1)
L111:
  add    $0x70,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
L28:
  mov    (%rsp),%eax
  test   %eax,%eax
  jle    L111
  lea    0x0(,%rax,8),%r8d
  lea    -0x2(%r14,%r14,1),%eax
  mov    %eax,0x4(%rsp)
L47:
  movq   $0x1,-0x20(%rsp)
  mov    -0x48(%rsp),%edi
  movslq %r14d,%rbx
  mov    %rbx,0x10(%rsp)
  lea    0x2(%rdi),%eax
  mov    %eax,0x8(%rsp)
  mov    %edi,%eax
  sub    %r14d,%eax
  mov    %eax,0xc(%rsp)
  movslq (%rsp),%rax
  lea    (%rsi,%rax,4),%rax
  movslq %edi,%rsi
  mov    %rbx,%rdi
  mov    %rax,-0x18(%rsp)
  mov    %rbx,%rax
  lea    (%rbx,%rbx,1),%rbx
  neg    %rax
  mov    %rbx,0x18(%rsp)
  shl    $0x3,%rax
  mov    %rsi,0x20(%rsp)
  mov    %rax,0x28(%rsp)
  lea    0x0(,%rsi,8),%rax
  mov    %rax,-0x10(%rsp)
  mov    %rbx,%rax
  add    %rdi,%rax
  mov    %rax,0x38(%rsp)
  shl    $0x3,%rax
  mov    %rax,0x30(%rsp)
  mov    %rbx,%rax
  add    %rsi,%rax
  mov    %rax,0x40(%rsp)
  add    %rdi,%rax
  mov    %rax,-0x8(%rsp)
  mov    %rsi,%rax
  sub    %rbx,%rax
  sub    %rdi,%rax
  mov    %rax,0x48(%rsp)
  mov    $0x2,%eax
  sub    %rdi,%rax
  mov    %rax,0x50(%rsp)
  lea    0x2(%rdi,%rbx,1),%rax
  sub    %rsi,%rbx
  mov    %rax,0x58(%rsp)
  mov    %rbx,%rax
  sub    %rdi,%rax
  mov    %r14d,%edi
  mov    %rax,0x60(%rsp)
  cs nopw 0x0(%rax,%rax,1)
L46:
  mov    -0x20(%rsp),%rbx
  mov    -0x18(%rsp),%r10
  mov    -0x48(%rsp),%r11d
  mov    -0x4(%r10,%rbx,4),%eax
  lea    (%rbx,%rbx,1),%esi
  mov    %esi,-0x44(%rsp)
  lea    -0x2(%rsi,%rax,1),%eax
  add    %eax,%eax
  lea    0x2(%rax),%edx
  add    %r11d,%eax
  movslq %edx,%rcx
  add    %edi,%edx
  lea    0x0(%rbp,%rcx,8),%rsi
  movslq %eax,%rcx
  add    %r8d,%eax
  lea    0x0(%rbp,%rcx,8),%rcx
  movupd (%rsi),%xmm0
  movupd (%rcx),%xmm7
  movups %xmm7,(%rsi)
  movups %xmm0,(%rcx)
  movslq %edx,%rcx
  add    %edi,%edx
  lea    0x0(%rbp,%rcx,8),%rsi
  movslq %eax,%rcx
  sub    %edi,%eax
  lea    0x0(%rbp,%rcx,8),%rcx
  movupd (%rsi),%xmm0
  movupd (%rcx),%xmm7
  movups %xmm7,(%rsi)
  movups %xmm0,(%rcx)
  movslq %edx,%rcx
  sub    $0x2,%edx
  lea    0x0(%rbp,%rcx,8),%rsi
  movslq %eax,%rcx
  lea    0x0(%rbp,%rcx,8),%rcx
  movupd (%rsi),%xmm0
  movupd (%rcx),%xmm5
  movups %xmm5,(%rsi)
  movups %xmm0,(%rcx)
  movslq %edx,%rcx
  add    0x8(%rsp),%edx
  lea    0x0(%rbp,%rcx,8),%rsi
  mov    %eax,%ecx
  add    $0x2,%eax
  sub    %r11d,%ecx
  movupd (%rsi),%xmm0
  movslq %ecx,%rcx
  lea    0x0(%rbp,%rcx,8),%rcx
  movupd (%rcx),%xmm6
  movups %xmm6,(%rsi)
  movups %xmm0,(%rcx)
  movslq %edx,%rcx
  sub    0xc(%rsp),%edx
  lea    0x0(%rbp,%rcx,8),%rsi
  movslq %eax,%rcx
  add    0x4(%rsp),%eax
  movslq %edx,%rdx
  lea    0x0(%rbp,%rcx,8),%rcx
  movupd (%rsi),%xmm0
  cltq
  lea    0x0(%rbp,%rdx,8),%rdx
  movupd (%rcx),%xmm7
  lea    0x0(%rbp,%rax,8),%rax
  movups %xmm7,(%rsi)
  movups %xmm0,(%rcx)
  movupd (%rdx),%xmm0
  movupd (%rax),%xmm5
  movups %xmm5,(%rdx)
  movups %xmm0,(%rax)
  cmp    %ebx,(%rsp)
  jle    L111
  mov    (%r10,%rbx,4),%esi
  mov    0x38(%rsp),%r9
  lea    (%rsi,%rsi,1),%eax
  mov    %esi,-0x50(%rsp)
  mov    0x10(%rsp),%rsi
  movslq %eax,%rdx
  lea    (%rsi,%rdx,1),%rcx
  lea    0x0(%rbp,%rcx,8),%rcx
  mov    %rcx,-0x70(%rsp)
  mov    0x18(%rsp),%rcx
  add    %rdx,%rcx
  lea    0x0(%rbp,%rcx,8),%r14
  lea    (%r9,%rdx,1),%rcx
  lea    0x0(%rbp,%rcx,8),%r13
  mov    -0x8(%rsp),%rcx
  add    %rdx,%rcx
  add    0x40(%rsp),%rdx
  lea    0x0(%rbp,%rcx,8),%r15
  lea    0x0(%rbp,%rdx,8),%r12
  mov    %r15,-0x68(%rsp)
  lea    (%rax,%rdi,2),%r15d
  lea    (%rdi,%r15,1),%eax
  mov    %r15d,-0x58(%rsp)
  movslq %eax,%rdx
  add    %r11d,%eax
  add    0x48(%rsp),%rdx
  lea    0x0(%rbp,%rdx,8),%rbx
  movslq %eax,%rdx
  sub    %edi,%eax
  sub    %r9,%rdx
  lea    0x0(%rbp,%rdx,8),%r11
  movslq %eax,%rdx
  sub    %edi,%eax
  add    0x50(%rsp),%rdx
  lea    0x0(%rbp,%rdx,8),%r10
  movslq %eax,%rdx
  sub    %edi,%eax
  add    %rsi,%rdx
  movslq %eax,%r15
  lea    0x2(%rdi,%rax,1),%eax
  lea    0x0(%rbp,%rdx,8),%r9
  mov    0x58(%rsp),%rdx
  add    %r15,%rdx
  lea    0x0(%rbp,%rdx,8),%rdx
  mov    %rdx,-0x60(%rsp)
  movslq %eax,%rdx
  add    0x60(%rsp),%rdx
  add    %edi,%eax
  lea    0x0(%rbp,%rdx,8),%rdx
  movslq %eax,%rcx
  add    %edi,%eax
  mov    %rdx,-0x78(%rsp)
  mov    0x20(%rsp),%rdx
  cltq
  sub    -0x8(%rsp),%rax
  neg    %rdx
  sub    %rsi,%rdx
  add    %rcx,%rdx
  lea    0x0(%rbp,%rax,8),%rcx
  movslq -0x50(%rsp),%rax
  lea    0x0(%rbp,%rdx,8),%rsi
  mov    %rcx,-0x50(%rsp)
  shl    $0x4,%rax
  lea    0x0(%rbp,%rax,1),%rdx
  movslq -0x58(%rsp),%rax
  mov    %rdx,-0x40(%rsp)
  mov    0x28(%rsp),%rdx
  add    -0x10(%rsp),%rdx
  lea    (%rdx,%rax,8),%rax
  lea    0x0(%rbp,%rax,1),%rdx
  mov    %rdx,-0x28(%rsp)
  mov    0x30(%rsp),%rdx
  sub    -0x10(%rsp),%rdx
  lea    (%rdx,%r15,8),%rax
  lea    0x0(%rbp,%rax,1),%rdx
  mov    -0x20(%rsp),%rax
  mov    %rdx,-0x30(%rsp)
  xor    %edx,%edx
  shl    $0x5,%rax
  mov    %rax,-0x38(%rsp)
  mov    -0x18(%rsp),%rax
  mov    %rax,-0x58(%rsp)
  nopl   0x0(%rax)
L119:
  mov    -0x44(%rsp),%eax
  add    $0x20,%r14
  add    $0x20,%r13
  add    $0x20,%r12
  add    $0x20,%rbx
  add    $0x20,%r11
  add    $0x20,%r10
  add    $0x20,%r9
  mov    -0x58(%rsp),%rcx
  add    $0x20,%rsi
  add    (%rcx),%eax
  mov    -0x40(%rsp),%rcx
  add    %eax,%eax
  movslq %eax,%r15
  movupd (%rcx,%rdx,1),%xmm0
  add    %r8d,%eax
  lea    0x0(%rbp,%r15,8),%r15
  movupd (%r15),%xmm2
  movups %xmm2,(%rcx,%rdx,1)
  mov    -0x70(%rsp),%rcx
  movups %xmm0,(%r15)
  movslq %eax,%r15
  sub    %edi,%eax
  lea    0x0(%rbp,%r15,8),%r15
  movupd (%rcx),%xmm0
  movupd (%r15),%xmm3
  movups %xmm3,(%rcx)
  mov    -0x68(%rsp),%rcx
  movups %xmm0,(%r15)
  movslq %eax,%r15
  movupd -0x20(%r14),%xmm0
  add    %r8d,%eax
  lea    0x0(%rbp,%r15,8),%r15
  movupd (%r15),%xmm4
  movups %xmm4,-0x20(%r14)
  movups %xmm0,(%r15)
  movslq %eax,%r15
  movupd -0x20(%r13),%xmm0
  add    $0x2,%eax
  lea    0x0(%rbp,%r15,8),%r15
  movupd (%r15),%xmm5
  movups %xmm5,-0x20(%r13)
  movups %xmm0,(%r15)
  movslq %eax,%r15
  movupd (%rcx),%xmm0
  sub    %r8d,%eax
  lea    0x0(%rbp,%r15,8),%r15
  movupd (%r15),%xmm6
  movups %xmm6,(%rcx)
  mov    -0x28(%rsp),%rcx
  movups %xmm0,(%r15)
  movslq %eax,%r15
  movupd -0x20(%r12),%xmm0
  add    %edi,%eax
  lea    0x0(%rbp,%r15,8),%r15
  movupd (%r15),%xmm7
  movups %xmm7,-0x20(%r12)
  movups %xmm0,(%r15)
  movslq %eax,%r15
  movupd (%rcx,%rdx,1),%xmm0
  sub    %r8d,%eax
  lea    0x0(%rbp,%r15,8),%r15
  movupd (%r15),%xmm2
  movups %xmm2,(%rcx,%rdx,1)
  mov    -0x60(%rsp),%rcx
  movups %xmm0,(%r15)
  movslq %eax,%r15
  movupd -0x20(%rbx),%xmm0
  add    -0x48(%rsp),%eax
  lea    0x0(%rbp,%r15,8),%r15
  movupd (%r15),%xmm3
  movups %xmm3,-0x20(%rbx)
  movups %xmm0,(%r15)
  movslq %eax,%r15
  movupd -0x10(%r11),%xmm0
  add    %r8d,%eax
  lea    0x0(%rbp,%r15,8),%r15
  movupd (%r15),%xmm4
  movups %xmm4,-0x10(%r11)
  movups %xmm0,(%r15)
  movslq %eax,%r15
  movupd -0x20(%r10),%xmm0
  sub    %edi,%eax
  lea    0x0(%rbp,%r15,8),%r15
  movupd (%r15),%xmm5
  movups %xmm5,-0x20(%r10)
  movups %xmm0,(%r15)
  movslq %eax,%r15
  movupd -0x10(%r9),%xmm0
  add    %r8d,%eax
  lea    0x0(%rbp,%r15,8),%r15
  movupd (%r15),%xmm6
  movups %xmm6,-0x10(%r9)
  movups %xmm0,(%r15)
  movslq %eax,%r15
  movupd (%rcx),%xmm0
  sub    $0x2,%eax
  lea    0x0(%rbp,%r15,8),%r15
  movupd (%r15),%xmm7
  movups %xmm7,(%rcx)
  mov    -0x30(%rsp),%rcx
  movups %xmm0,(%r15)
  movslq %eax,%r15
  sub    %r8d,%eax
  lea    0x0(%rbp,%r15,8),%r15
  movupd 0x10(%rcx,%rdx,1),%xmm0
  movupd (%r15),%xmm1
  movups %xmm1,0x10(%rcx,%rdx,1)
  mov    -0x78(%rsp),%rcx
  add    $0x20,%rdx
  movups %xmm0,(%r15)
  movslq %eax,%r15
  add    %edi,%eax
  lea    0x0(%rbp,%r15,8),%r15
  movupd (%rcx),%xmm0
  movupd (%r15),%xmm2
  addq   $0x4,-0x58(%rsp)
  addq   $0x20,-0x70(%rsp)
  movups %xmm2,(%rcx)
  mov    -0x50(%rsp),%rcx
  movups %xmm0,(%r15)
  movslq %eax,%r15
  sub    %r8d,%eax
  movupd -0x20(%rsi),%xmm0
  lea    0x0(%rbp,%r15,8),%r15
  cltq
  add    $0x20,%rcx
  addq   $0x20,-0x68(%rsp)
  movupd (%r15),%xmm3
  lea    0x0(%rbp,%rax,8),%rax
  addq   $0x20,-0x60(%rsp)
  addq   $0x20,-0x78(%rsp)
  movups %xmm3,-0x20(%rsi)
  movups %xmm0,(%r15)
  movupd (%rax),%xmm4
  movupd -0x20(%rcx),%xmm0
  movups %xmm4,-0x20(%rcx)
  mov    %rcx,-0x50(%rsp)
  movups %xmm0,(%rax)
  cmp    -0x38(%rsp),%rdx
  jne    L119
  addq   $0x1,-0x20(%rsp)
  jmp    L46
L105:
  je     L135
  movl   $0x1,(%rsp)
  mov    $0x4,%r14d
  jmp    L121
L135:
  movl   $0x6,0x4(%rsp)
  mov    $0x4,%r14d
  movl   $0x1,(%rsp)
  jmp    L47
```

assistant: ```c
void bitrv2(int n, int *ip, double *a) {
  int j, j1, k, k1, l, m, nh, nm;
  double xr, xi, yr, yi;
  m = 1;
  for (l = n >> 2; l > 8; l >>= 2) {
    m <<= 1;
  }
  nh = n >> 1;
  nm = 4 * m;
  if (l == 8) {
    for (k = 0; k < m; k++) {
      for (j = 0; j < k; j++) {
        j1 = 4 * j + 2 * ip[m + k];
        k1 = 4 * k + 2 * ip[m + j];
        xr = a[j1];
        xi = a[j1 + 1];
        yr = a[k1];
        yi = a[k1 + 1];
        a[j1] = yr;
        a[j1 + 1] = yi;
        a[k1] = xr;
        a[k1 + 1] = xi;
        j1 += nm;
        k1 += 2 * nm;
        xr = a[j1];
        xi = a[j1 + 1];
        yr = a[k1];
        yi = a[k1 + 1];
        a[j1] = yr;
        a[j1 + 1] = yi;
        a[k1] = xr;
        a[k1 + 1] = xi;
        j1 += nm;
        k1 -= nm;
        xr = a[j1];
        xi = a[j1 + 1];
        yr = a[k1];
        yi = a[k1 + 1];
        a[j1] = yr;
        a[j1 + 1] = yi;
        a[k1] = xr;
        a[k1 + 1] = xi;
        j1 += nm;
        k1 += 2 * nm;
        xr = a[j1];
        xi = a[j1 + 1];
        yr = a[k1];
        yi = a[k1 + 1];
        a[j1] = yr;
        a[j1 + 1] = yi;
        a[k1] = xr;
        a[k1 + 1] = xi;
        j1 += nh;
        k1 += 2;
        xr = a[j1];
        xi = a[j1 + 1];
        yr = a[k1];
        yi = a[k1 + 1];
        a[j1] = yr;
        a[j1 + 1] = yi;
        a[k1] = xr;
        a[k1 + 1] = xi;
        j1 -= nm;
        k1 -= 2 * nm;
        xr = a[j1];
        xi = a[j1 + 1];
        yr = a[k1];
        yi = a[k1 + 1];
        a[j1] = yr;
        a[j1 + 1] = yi;
        a[k1] = xr;
        a[k1 + 1] = xi;
        j1 -= nm;
        k1 += nm;
        xr = a[j1];
        xi = a[j1 + 1];
        yr = a[k1];
        yi = a[k1 + 1];
        a[j1] = yr;
        a[j1 + 1] = yi;
        a[k1] = xr;
        a[k1 + 1] = xi;
        j1 -= nm;
        k1 -= 2 * nm;
        xr = a[j1];
        xi = a[j1 + 1];
        yr = a[k1];
        yi = a[k1 + 1];
        a[j1] = yr;
        a[j1 + 1] = yi;
        a[k1] = xr;
        a[k1 + 1] = xi;
        j1 += 2;
        k1 += nh;
        xr = a[j1];
        xi = a[j1 + 1];
        yr = a[k1];
        yi = a[k1 + 1];
        a[j1] = yr;
        a[j1 + 1] = yi;
        a[k1] = xr;
        a[k1 + 1] = xi;
        j1 += nm;
        k1 += 2 * nm;
        xr = a[j1];
        xi = a[j1 + 1];
        yr = a[k1];
        yi = a[k1 + 1];
        a[j1] = yr;
        a[j1 + 1] = yi;
        a[k1] = xr;
        a[k1 + 1] = xi;
        j1 += nm;
        k1 -= nm;
        xr = a[j1];
        xi = a[j1 + 1];
        yr = a[k1];
        yi = a[k1 + 1];
        a[j1] = yr;
        a[j1 + 1] = yi;
        a[k1] = xr;
        a[k1 + 1] = xi;
        j1 += nm;
        k1 += 2 * nm;
        xr = a[j1];
        xi = a[j1 + 1];
        yr = a[k1];
        yi = a[k1 + 1];
        a[j1] = yr;
        a[j1 + 1] = yi;
        a[k1] = xr;
        a[k1 + 1] = xi;
        j1 -= nh;
        k1 -= 2;
        xr = a[j1];
        xi = a[j1 + 1];
        yr = a[k1];
        yi = a[k1 + 1];
        a[j1] = yr;
        a[j1 + 1] = yi;
        a[k1] = xr;
        a[k1 + 1] = xi;
        j1 -= nm;
        k1 -= 2 * nm;
        xr = a[j1];
        xi = a[j1 + 1];
        yr = a[k1];
        yi = a[k1 + 1];
        a[j1] = yr;
        a[j1 + 1] = yi;
        a[k1] = xr;
        a[k1 + 1] = xi;
        j1 -= nm;
        k1 += nm;
        xr = a[j1];
        xi = a[j1 + 1];
        yr = a[k1];
        yi = a[k1 + 1];
        a[j1] = yr;
        a[j1 + 1] = yi;
        a[k1] = xr;
        a[k1 + 1] = xi;
        j1 -= nm;
        k1 -= 2 * nm;
        xr = a[j1];
        xi = a[j1 + 1];
        yr = a[k1];
        yi = a[k1 + 1];
        a[j1] = yr;
        a[j1 + 1] = yi;
        a[k1] = xr;
        a[k1 + 1] = xi;
      }
      k1 = 4 * k + 2 * ip[m + k];
      j1 = k1 + 2;
      k1 += nh;
      xr = a[j1];
      xi = a[j1 + 1];
      yr = a[k1];
      yi = a[k1 + 1];
      a[j1] = yr;
      a[j1 + 1] = yi;
      a[k1] = xr;
      a[k1 + 1] = xi;
      j1 += nm;
      k1 += 2 * nm;
      xr = a[j1];
      xi = a[j1 + 1];
      yr = a[k1];
      yi = a[k1 + 1];
      a[j1] = yr;
      a[j1 + 1] = yi;
      a[k1] = xr;
      a[k1 + 1] = xi;
      j1 += nm;
      k1 -= nm;
      xr = a[j1];
      xi = a[j1 + 1];
      yr = a[k1];
      yi = a[k1 + 1];
      a[j1] = yr;
      a[j1 + 1] = yi;
      a[k1] = xr;
      a[k1 + 1] = xi;
      j1 -= 2;
      k1 -= nh;
      xr = a[j1];
      xi = a[j1 + 1];
      yr = a[k1];
      yi = a[k1 + 1];
      a[j1] = yr;
      a[j1 + 1] = yi;
      a[k1] = xr;
      a[k1 + 1] = xi;
      j1 += nh + 2;
      k1 += nh + 2;
      xr = a[j1];
      xi = a[j1 + 1];
      yr = a[k1];
      yi = a[k1 + 1];
      a[j1] = yr;
      a[j1 + 1] = yi;
      a[k1] = xr;
      a[k1 + 1] = xi;
      j1 -= nh - nm;
      k1 += 2 * nm - 2;
      xr = a[j1];
      xi = a[j1 + 1];
      yr = a[k1];
      yi = a[k1 + 1];
      a[j1] = yr;
      a[j1 + 1] = yi;
      a[k1] = xr;
      a[k1 + 1] = xi;
    }
  } else {
    for (k = 0; k < m; k++) {
      for (j = 0; j < k; j++) {
        j1 = 4 * j + ip[m + k];
        k1 = 4 * k + ip[m + j];
        xr = a[j1];
        xi = a[j1 + 1];
        yr = a[k1];
        yi = a[k1 + 1];
        a[j1] = yr;
        a[j1 + 1] = yi;
        a[k1] = xr;
        a[k1 + 1] = xi;
        j1 += nm;
        k1 += nm;
        xr = a[j1];
        xi = a[j1 + 1];
        yr = a[k1];
        yi = a[k1 + 1];
        a[j1] = yr;
        a[j1 + 1] = yi;
        a[k1] = xr;
        a[k1 + 1] = xi;
        j1 += nh;
        k1 += 2;
        xr = a[j1];
        xi = a[j1 + 1];
        yr = a[k1];
        yi = a[k1 + 1];
        a[j1] = yr;
        a[j1 + 1] = yi;
        a[k1] = xr;
        a[k1 + 1] = xi;
        j1 -= nm;
        k1 -= nm;
        xr = a[j1];
        xi = a[j1 + 1];
        yr = a[k1];
        yi = a[k1 + 1];
        a[j1] = yr;
        a[j1 + 1] = yi;
        a[k1] = xr;
        a[k1 + 1] = xi;
        j1 += 2;
        k1 += nh;
        xr = a[j1];
        xi = a[j1 + 1];
        yr = a[k1];
        yi = a[k1 + 1];
        a[j1] = yr;
        a[j1 + 1] = yi;
        a[k1] = xr;
        a[k1 + 1] = xi;
        j1 += nm;
        k1 += nm;
        xr = a[j1];
        xi = a[j1 + 1];
        yr = a[k1];
        yi = a[k1 + 1];
        a[j1] = yr;
        a[j1 + 1] = yi;
        a[k1] = xr;
        a[k1 + 1] = xi;
        j1 -= nh;
        k1 -= 2;
        xr = a[j1];
        xi = a[j1 + 1];
        yr = a[k1];
        yi = a[k1 + 1];
        a[j1] = yr;
        a[j1 + 1] = yi;
        a[k1] = xr;
        a[k1 + 1] = xi;
        j1 -= nm;
        k1 -= nm;
        xr = a[j1];
        xi = a[j1 + 1];
        yr = a[k1];
        yi = a[k1 + 1];
        a[j1] = yr;
        a[j1 + 1] = yi;
        a[k1] = xr;
        a[k1 + 1] = xi;
      }
      k1 = 4 * k + ip[m + k];
      j1 = k1 + 2;
      k1 += nh;
      xr = a[j1];
      xi = a[j1 + 1];
      yr = a[k1];
      yi = a[k1 + 1];
      a[j1] = yr;
      a[j1 + 1] = yi;
      a[k1] = xr;
      a[k1 + 1] = xi;
      j1 += nm;
      k1 += nm;
      xr = a[j1];
      xi = a[j1 + 1];
      yr = a[k1];
      yi = a[k1 + 1];
      a[j1] = yr;
      a[j1 + 1] = yi;
      a[k1] = xr;
      a[k1 + 1] = xi;
    }
  }
}
```