user: 下面的汇编指令对应的C源代码是什么:

```asm
<sprstp>:
  endbr64
  push   %r15
  push   %r14
  push   %r13
  push   %r12
  push   %rbp
  push   %rbx
  sub    $0x38,%rsp
  mov    %rdi,%r15
  mov    %rsi,%rbx
  mov    %rdx,%r12
  mov    %rcx,%r13
  mov    0x8(%rsi),%rbp
  mov    %rbp,%rax
  sub    $0x2,%rax
  mov    %rax,0x10(%rsp)
  je     L161
  lea    -0x1(%rbp),%rdx
  mov    $0x1,%eax
L13:
  movss  (%r15,%rax,4),%xmm0
  movss  %xmm0,(%r12,%rax,4)
  add    $0x1,%rax
  cmp    %rdx,%rax
  jne    L13
L161:
  lea    0x0(,%rbp,8),%rax
  mov    %rax,0x8(%rsp)
  sub    $0x8,%rax
  lea    (%rbx,%rax,1),%r14
  lea    0x0(%r13,%rax,1),%rdx
  mov    (%r14),%rdi
  sub    0x8(%rbx),%rdi
  mov    %r14,%rsi
  call   <iindexx@plt>
  mov    0x8(%rbx),%r8
  mov    (%r14),%rax
  sub    $0x1,%rax
  cmp    %rax,%r8
  ja     L208
  mov    $0x0,%edi
  lea    -0x1(%rbp),%r9
  mov    $0x1,%r10d
  jmp    L207
L201:
  lea    (%rdx,%rcx,1),%rax
  shr    %rax
  cmp    %rsi,(%rbx,%rax,8)
  cmovbe %rax,%rdx
  cmova  %rax,%rcx
L225:
  mov    %rcx,%rax
  sub    %rdx,%rax
  cmp    $0x1,%rax
  ja     L201
  mov    %rdx,0x0(%r13,%r8,8)
  add    $0x1,%r8
  mov    (%r14),%rax
  sub    $0x1,%rax
  cmp    %r8,%rax
  jb     L181
L207:
  mov    %r9,%rcx
  mov    %r9,%rsi
  add    0x0(%r13,%r8,8),%rsi
  movss  (%r15,%rsi,4),%xmm0
  movss  %xmm0,(%r12,%r8,4)
  lea    0x1(%rdi),%rax
  lea    (%rbx,%rsi,8),%rdx
  mov    (%rdx),%rdi
  cmp    %rdi,%rax
  ja     L172
L43:
  mov    %r8,0x0(%r13,%rax,8)
  add    $0x1,%rax
  mov    (%rdx),%rdi
  cmp    %rax,%rdi
  jae    L43
L172:
  mov    %r10,%rdx
  jmp    L225
L208:
  mov    $0x0,%edi
L181:
  add    $0x1,%rdi
  cmp    %rdi,%rbp
  jbe    L55
  lea    0x0(%r13,%rdi,8),%rax
  mov    0x8(%rsp),%r15
  add    %r13,%r15
L146:
  mov    (%r14),%rdx
  mov    %rdx,(%rax)
  add    $0x8,%rax
  cmp    %r15,%rax
  jne    L146
L55:
  mov    0x10(%rsp),%rax
  test   %rax,%rax
  je     L106
  lea    0x8(%r13),%rbx
  mov    %rbx,0x20(%rsp)
  lea    (%rbx,%rax,8),%rax
  mov    %rax,0x28(%rsp)
  jmp    L220
L52:
  mov    %r11,0x0(%r13,%r10,8)
  movss  %xmm1,(%r12,%r10,4)
  add    $0x1,%rbp
  add    $0x8,%r15
  add    $0x4,%r14
  cmp    0x8(%rsp),%rbp
  ja     L218
L177:
  mov    (%r15),%r11
  movss  (%r14),%xmm1
  mov    %r14,%rsi
  mov    %r15,%rcx
  mov    %rbp,%rdx
L131:
  mov    %rdx,%r10
  sub    %rax,%rdx
  mov    (%rcx,%r9,1),%rdi
  cmp    %r11,%rdi
  jbe    L52
  mov    %rdi,(%rcx)
  movss  (%rsi,%r8,1),%xmm0
  movss  %xmm0,(%rsi)
  add    %r9,%rcx
  add    %r8,%rsi
  lea    (%rbx,%rdx,1),%rdi
  cmp    %rdi,%rax
  jb     L131
  mov    %rdx,%r10
  jmp    L52
L218:
  cmpq   $0x5,0x10(%rsp)
  jbe    L170
L118:
  mov    %rax,0x10(%rsp)
  movabs $0xaaaaaaaaaaaaaaab,%rcx
  mul    %rcx
  shr    %rdx
  mov    %rdx,%rax
  mov    0x18(%rsp),%rcx
  lea    (%rcx,%rdx,1),%rbp
  cmp    0x8(%rsp),%rbp
  ja     L218
  lea    0x0(%r13,%rbp,8),%r15
  lea    (%r12,%rbp,4),%r14
  mov    %rdx,%r9
  neg    %r9
  shl    $0x3,%r9
  neg    %rdx
  lea    0x0(,%rdx,4),%r8
  jmp    L177
L170:
  addq   $0x8,0x20(%rsp)
  mov    0x20(%rsp),%rax
  cmp    0x28(%rsp),%rax
  je     L106
L220:
  mov    0x20(%rsp),%rax
  mov    0x8(%rax),%rcx
  mov    (%rax),%rax
  mov    %rax,0x18(%rsp)
  mov    %rcx,%rdx
  sub    %rax,%rdx
  mov    $0x1,%eax
L20:
  lea    0x1(%rax,%rax,2),%rax
  cmp    %rax,%rdx
  jae    L20
  lea    -0x1(%rcx),%rbx
  mov    %rbx,0x8(%rsp)
  mov    $0x1,%ebx
  sub    0x18(%rsp),%rbx
  jmp    L118
L106:
  add    $0x38,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
```

assistant: ```c
void sprstp(float sa[], unsigned long ija[], float sb[], unsigned long ijb[]) {
  void iindexx(unsigned long n, long arr[], unsigned long indx[]);
  unsigned long j, jl, jm, jp, ju, k, m, n2, noff, inc, iv;
  float v;
  n2 = ija[1];
  for (j = 1; j <= n2 - 2; j++)
    sb[j] = sa[j];
  iindexx(ija[n2 - 1] - ija[1], (long *)&ija[n2 - 1], &ijb[n2 - 1]);
  jp = 0;
  for (k = ija[1]; k <= ija[n2 - 1] - 1; k++) {
    m = ijb[k] + n2 - 1;
    sb[k] = sa[m];
    for (j = jp + 1; j <= ija[m]; j++)
      ijb[j] = k;
    jp = ija[m];
    jl = 1;
    ju = n2 - 1;
    while (ju - jl > 1) {
      jm = (ju + jl) / 2;
      if (ija[jm] > m)
        ju = jm;
      else
        jl = jm;
    }
    ijb[k] = jl;
  }
  for (j = jp + 1; j < n2; j++)
    ijb[j] = ija[n2 - 1];
  for (j = 1; j <= n2 - 2; j++) {
    jl = ijb[j + 1] - ijb[j];
    noff = ijb[j] - 1;
    inc = 1;
    do {
      inc *= 3;
      inc++;
    } while (inc <= jl);
    do {
      inc /= 3;
      for (k = noff + inc + 1; k <= noff + jl; k++) {
        iv = ijb[k];
        v = sb[k];
        m = k;
        while (ijb[m - inc] > iv) {
          ijb[m] = ijb[m - inc];
          sb[m] = sb[m - inc];
          m -= inc;
          if (m - noff <= inc)
            break;
        }
        ijb[m] = iv;
        sb[m] = v;
      }
    } while (inc > 1);
  }
}
```