user: 请从给定的汇编指令推导出C源代码：

```asm
<gaussj>:
  endbr64
  push   %r15
  push   %r14
  push   %r13
  push   %r12
  push   %rbp
  push   %rbx
  sub    $0x338,%rsp
  mov    %rdx,0x10(%rsp)
  mov    %fs:0x28,%rax
  mov    %rax,0x328(%rsp)
  xor    %eax,%eax
  mov    $0x1,%eax
  cmp    $0x40,%esi
  jg     L171
  mov    %rdi,%rbx
  mov    %ecx,%r11d
  test   %esi,%esi
  jle    L86
  lea    0x20(%rsp),%rax
  lea    -0x1(%rsi),%edx
  lea    0x24(%rsp,%rdx,4),%rdx
L146:
  movl   $0x0,(%rax)
  add    $0x4,%rax
  cmp    %rdx,%rax
  jne    L146
  mov    %esi,%edi
  movslq %esi,%r8
  shl    $0x3,%r8
  movq   $0x0,0x8(%rsp)
  movl   $0x0,0x18(%rsp)
  mov    $0x0,%r15d
  lea    0x20(%rsp),%rcx
  movq   D23(%rip),%xmm3
L396:
  mov    $0x0,%ebp
  pxor   %xmm1,%xmm1
  jmp    L96
L354:
  cmp    $0x1,%r9d
  jg     L326
L176:
  add    $0x1,%rax
  add    %esi,%edx
  cmp    %rdi,%rax
  je     L117
L270:
  mov    (%rcx,%rax,4),%r9d
  test   %r9d,%r9d
  jne    L354
  movslq %edx,%r9
  movsd  (%rbx,%r9,8),%xmm0
  andpd  %xmm3,%xmm0
  comisd %xmm1,%xmm0
  jb     L176
  movapd %xmm0,%xmm1
  mov    %r12d,0x18(%rsp)
  mov    %eax,%r15d
  jmp    L176
L117:
  add    $0x1,%rbp
  cmp    %rdi,%rbp
  je     L106
L96:
  mov    %ebp,%edx
  mov    %ebp,%r12d
  cmpl   $0x1,(%rcx,%rbp,4)
  je     L117
  mov    $0x0,%eax
  jmp    L270
L106:
  movslq %r15d,%rax
  addl   $0x1,0x20(%rsp,%rax,4)
  mov    0x18(%rsp),%eax
  cmp    %eax,%r15d
  jne    L318
L91:
  mov    0x8(%rsp),%rax
  mov    0x18(%rsp),%edx
  mov    %edx,0x220(%rsp,%rax,4)
  mov    %r15d,0x120(%rsp,%rax,4)
  mov    %r15d,%eax
  imul   %esi,%eax
  mov    %eax,0x1c(%rsp)
  add    %r15d,%eax
  cltq
  lea    (%rbx,%rax,8),%rax
  movsd  (%rax),%xmm1
  ucomisd D1(%rip),%xmm1
  jp     L151
  je     L306
L151:
  movsd  D15(%rip),%xmm0
  divsd  %xmm1,%xmm0
  movsd  D15(%rip),%xmm4
  movsd  %xmm4,(%rax)
  movslq %r15d,%r9
  lea    0x0(,%r9,8),%rax
  lea    (%rbx,%rax,1),%r10
  mov    $0x0,%edx
L141:
  movapd %xmm0,%xmm1
  mulsd  (%r10),%xmm1
  movsd  %xmm1,(%r10)
  mov    %edx,%ebp
  add    $0x1,%edx
  add    %r8,%r10
  cmp    %edx,%esi
  jne    L141
  test   %r11d,%r11d
  jle    L173
  add    0x10(%rsp),%rax
  mov    $0x0,%edx
L262:
  movapd %xmm0,%xmm1
  mulsd  (%rax),%xmm1
  movsd  %xmm1,(%rax)
  add    $0x1,%edx
  add    %r8,%rax
  cmp    %edx,%r11d
  jne    L262
L173:
  mov    %rbx,%r13
  mov    0x10(%rsp),%r14
  mov    $0x0,%r12d
  jmp    L375
L318:
  movslq %eax,%rbp
  lea    0x0(,%rbp,8),%rax
  lea    (%rbx,%rax,1),%rdx
  mov    $0x0,%r10d
  movslq %r15d,%r9
  sub    %rbp,%r9
L394:
  movsd  (%rdx),%xmm0
  movsd  (%rdx,%r9,8),%xmm1
  movsd  %xmm1,(%rdx)
  movsd  %xmm0,(%rdx,%r9,8)
  add    $0x1,%r10d
  add    %r8,%rdx
  cmp    %r10d,%esi
  jne    L394
  test   %r11d,%r11d
  jle    L91
  add    0x10(%rsp),%rax
  mov    $0x0,%edx
L186:
  movsd  (%rax),%xmm0
  movsd  (%rax,%r9,8),%xmm1
  movsd  %xmm1,(%rax)
  movsd  %xmm0,(%rax,%r9,8)
  add    $0x1,%edx
  add    %r8,%rax
  cmp    %edx,%r11d
  jne    L186
  jmp    L91
L217:
  lea    0x1(%r12),%eax
  add    $0x8,%r13
  sub    $0x1,%r9
  add    $0x8,%r14
  cmp    %r12d,%ebp
  je     L187
  mov    %eax,%r12d
L375:
  cmp    %r12d,%r15d
  je     L217
  mov    0x1c(%rsp),%eax
  add    %r12d,%eax
  cltq
  lea    (%rbx,%rax,8),%rax
  movsd  (%rax),%xmm1
  movq   $0x0,(%rax)
  mov    %r13,%rax
  mov    $0x0,%edx
L208:
  movapd %xmm1,%xmm2
  mulsd  (%rax,%r9,8),%xmm2
  movsd  (%rax),%xmm0
  subsd  %xmm2,%xmm0
  movsd  %xmm0,(%rax)
  mov    %edx,%r10d
  add    $0x1,%edx
  add    %r8,%rax
  cmp    %r10d,%ebp
  jne    L208
  test   %r11d,%r11d
  jle    L217
  mov    %r14,%rax
  mov    $0x0,%edx
L13:
  movapd %xmm1,%xmm2
  mulsd  (%rax,%r9,8),%xmm2
  movsd  (%rax),%xmm0
  subsd  %xmm2,%xmm0
  movsd  %xmm0,(%rax)
  add    $0x1,%edx
  add    %r8,%rax
  cmp    %edx,%r11d
  jne    L13
  jmp    L217
L187:
  addq   $0x1,0x8(%rsp)
  mov    0x8(%rsp),%rax
  cmp    %rdi,%rax
  jne    L396
  movslq %ebp,%rbp
  lea    -0x1(%rsi),%r9d
  lea    0x8(%rbx),%r8
  jmp    L130
L303:
  sub    $0x1,%rbp
  test   %ebp,%ebp
  js     L379
L130:
  mov    0x220(%rsp,%rbp,4),%ecx
  mov    0x120(%rsp,%rbp,4),%eax
  cmp    %eax,%ecx
  je     L303
  test   %esi,%esi
  jle    L303
  imul   %esi,%eax
  mov    %eax,%edx
  imul   %esi,%ecx
  movslq %ecx,%rcx
  lea    (%rbx,%rcx,8),%rax
  lea    (%r9,%rcx,1),%rdi
  lea    (%r8,%rdi,8),%rdi
  movslq %edx,%rdx
  sub    %rcx,%rdx
L397:
  movsd  (%rax),%xmm0
  movsd  (%rax,%rdx,8),%xmm1
  movsd  %xmm1,(%rax)
  movsd  %xmm0,(%rax,%rdx,8)
  add    $0x8,%rax
  cmp    %rax,%rdi
  jne    L397
  jmp    L303
L379:
  mov    $0x0,%eax
  jmp    L171
L86:
  mov    $0x0,%eax
  jmp    L171
L326:
  mov    $0x1,%eax
L171:
  mov    0x328(%rsp),%rdx
  sub    %fs:0x28,%rdx
  jne    L18
  add    $0x338,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
L306:
  mov    $0x1,%eax
  jmp    L171
L18:
  call   <__stack_chk_fail@plt>
```

tool_calls: name: parse_data
arguments: data_label: D23
data_type: qword
name: parse_data
arguments: data_label: D1
data_type: f64
name: parse_data
arguments: data_label: D15
data_type: f64

tool: D23: 9223372036854775807

tool: D1: 0

tool: D15: 1

assistant: ```c
int gaussj(double *a, int n, double *b, int m) {
  int icol = 0, ipiv[64], irow = 0, i, j, k, l, indxc[64], indxr[64], ll;
  double pivinv, big, dum;
  if (n > 64)
    return (1);
  for (j = 0; j < n; ++j)
    ipiv[j] = 0;
  for (i = 0; i < n; ++i) {
    big = 0.0;
    for (j = 0; j < n; ++j)
      if (ipiv[j] != 1)
        for (k = 0; k < n; ++k) {
          if (ipiv[k] == 0) {
            if (fabs(a[j + k * n]) >= big) {
              big = fabs(a[j + k * n]);
              irow = j;
              icol = k;
            }
          } else if (ipiv[k] > 1)
            return (1);
        }
    ++ipiv[icol];
    if (irow != icol) {
      for (l = 0; l < n; ++l) {
        dum = a[irow + l * n];
        a[irow + l * n] = a[icol + l * n];
        a[icol + l * n] = dum;
      }
      for (l = 0; l < m; ++l) {
        dum = b[irow + l * n];
        b[irow + l * n] = b[icol + l * n];
        b[icol + l * n] = dum;
      }
    }
    indxr[i] = irow;
    indxc[i] = icol;
    if (a[icol + icol * n] == 0.0)
      return (1);
    pivinv = 1.0 / a[icol + icol * n];
    a[icol + icol * n] = 1.0;
    for (l = 0; l < n; ++l)
      a[icol + l * n] *= pivinv;
    for (l = 0; l < m; ++l)
      b[icol + l * n] *= pivinv;
    for (ll = 0; ll < n; ++ll)
      if (ll != icol) {
        dum = a[ll + icol * n];
        a[ll + icol * n] = 0.0;
        for (l = 0; l < n; ++l)
          a[ll + l * n] -= a[icol + l * n] * dum;
        for (l = 0; l < m; ++l)
          b[ll + l * n] -= b[icol + l * n] * dum;
      }
  }
  for (l = n - 1; l >= 0; --l)
    if (indxr[l] != indxc[l])
      for (k = 0; k < n; ++k) {
        dum = a[k + indxr[l] * n];
        a[k + indxr[l] * n] = a[k + indxc[l] * n];
        a[k + indxc[l] * n] = dum;
      }
  return (0);
}
```