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

```asm
<invband_compute>:
  endbr64
  push   %r15
  push   %r14
  push   %r13
  push   %r12
  push   %rbp
  push   %rbx
  mov    %ecx,-0x4(%rsp)
  test   %edx,%edx
  jle    L121
  mov    %rdi,%r9
  mov    %rsi,%r13
  mov    %edx,%r11d
  mov    %ecx,%ebx
  movslq %edx,%rsi
  shl    $0x3,%rsi
  mov    %edx,%r10d
  mov    %rsi,%r8
  mov    $0x0,%edi
  movsd  D5(%rip),%xmm3
  lea    0x1(%rcx),%ecx
  jmp    L123
L144:
  add    %rsi,%rax
L132:
  add    $0x1,%edx
  cmp    %ecx,%edx
  je     L55
  movsd  0x0(%r13,%rax,1),%xmm0
  divsd  %xmm1,%xmm0
  movsd  %xmm0,0x0(%r13,%rax,1)
  test   %edx,%edx
  jle    L144
  movq   $0x0,(%r9,%rax,1)
  jmp    L144
L55:
  add    $0x1,%rdi
  add    $0x8,%r8
  cmp    %r10,%rdi
  je     L84
L123:
  movsd  0x0(%r13,%rdi,8),%xmm1
  movapd %xmm1,%xmm2
  mulsd  %xmm1,%xmm2
  movapd %xmm3,%xmm0
  divsd  %xmm2,%xmm0
  movsd  %xmm0,(%r9,%rdi,8)
  test   %ebx,%ebx
  js     L55
  movsd  0x0(%r13,%rdi,8),%xmm0
  divsd  %xmm1,%xmm0
  movsd  %xmm0,0x0(%r13,%rdi,8)
  mov    %r8,%rax
  mov    $0x0,%edx
  jmp    L132
L54:
  mov    %edx,%eax
  neg    %eax
  lea    (%rdx,%rbp,1),%edi
L47:
  imul   %r11d,%eax
  add    %edi,%eax
  cltq
  movsd  (%r9,%rax,8),%xmm1
  mulsd  (%rcx),%xmm1
  movsd  (%rsi),%xmm0
  subsd  %xmm1,%xmm0
  movsd  %xmm0,(%rsi)
  add    $0x1,%edx
  add    %r12,%rcx
  cmp    %r10d,%edx
  je     L93
L101:
  lea    (%rdx,%rbx,1),%eax
  cmp    %eax,%r8d
  jge    L54
  mov    %edx,%eax
  mov    %r15d,%edi
  jmp    L47
L93:
  lea    0x1(%r8),%eax
  add    $0x1,%ebp
  sub    $0x1,%r10d
  add    %r11d,%r14d
  cmp    -0x14(%rsp),%eax
  je     L10
  mov    %eax,%r8d
L140:
  movslq %r14d,%rax
  lea    (%r9,%rax,8),%rsi
  mov    %r8d,%ebx
  mov    $0x1,%edx
  sub    %r8d,%edx
  mov    -0x28(%rsp),%rax
  lea    0x0(%r13,%rax,1),%rcx
  mov    %ebp,%r15d
  jmp    L101
L10:
  mov    -0x10(%rsp),%r15
  mov    -0x28(%rsp),%rax
  mov    $0x1,%edx
L167:
  movsd  0x0(%r13,%rax,1),%xmm1
  mulsd  (%r9,%rax,1),%xmm1
  movsd  (%r9,%r15,8),%xmm0
  subsd  %xmm1,%xmm0
  movsd  %xmm0,(%r9,%r15,8)
  mov    %edx,%ecx
  add    $0x1,%edx
  add    -0x20(%rsp),%rax
  cmp    %ecx,%r8d
  jne    L167
L105:
  sub    $0x1,%r15
  subq   $0x8,-0x28(%rsp)
  test   %r15d,%r15d
  js     L121
L149:
  mov    -0x8(%rsp),%eax
  sub    %r15d,%eax
  mov    -0x4(%rsp),%ebx
  cmp    %ebx,%eax
  cmovg  %ebx,%eax
  test   %eax,%eax
  jle    L105
  lea    0x1(%r15),%ebp
  mov    %eax,%r10d
  lea    (%r11,%r15,1),%r14d
  add    $0x1,%eax
  mov    %eax,-0x14(%rsp)
  mov    $0x1,%r8d
  mov    -0x20(%rsp),%r12
  mov    %r15,-0x10(%rsp)
  jmp    L140
L84:
  lea    -0x1(%r11),%eax
  mov    %eax,-0x8(%rsp)
  movslq %eax,%r15
  movslq %r11d,%rax
  lea    (%r15,%rax,1),%rdx
  lea    0x0(,%rdx,8),%rbx
  mov    %rbx,-0x28(%rsp)
  shl    $0x3,%rax
  mov    %rax,-0x20(%rsp)
  jmp    L149
L121:
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
```

tool_calls: name: read_data
arguments: data_tag: D5
data_type: f64

tool: D5: 1

assistant: ```c
void invband_compute(double *dataptr, double *L, int n, int m) {
  int i, j, k;
  int idx, idy;
  double *data, *odata;
  double diag;
  data = dataptr;
  for (i = 0; i < n; i++) {
    diag = L[i];
    data[i] = 1.0 / (diag * diag);
    for (j = 0; j <= m; j++) {
      L[j * n + i] /= diag;
      if (j > 0) {
        data[j * n + i] = 0;
      }
    }
  }
  for (i = n - 1; i >= 0; i--) {
    for (j = 1; j <= (m < n - 1 - i ? m : n - 1 - i); j++) {
      for (k = 1; k <= (n - 1 - i < m ? n - 1 - i : m); k++) {
        idx = (j < k ? k - j : j - k);
        idy = (j < k ? i + j : i + k);
        data[j * n + i] -= L[k * n + i] * data[idx * n + idy];
      }
    }
    for (k = 1; k <= (n - 1 - i < m ? n - 1 - i : m); k++) {
      data[i] -= L[k * n + i] * data[k * n + i];
    }
  }
  return;
}
```