user: 请显示与下面列出的汇编代码相匹配的C语言实现:

```asm
<gaussj>:
  endbr64
  push   %r15
  push   %r14
  push   %r13
  push   %r12
  mov    %esi,%r12d
  push   %rbp
  push   %rbx
  sub    $0x78,%rsp
  mov    %rdi,0x20(%rsp)
  mov    $0x1,%edi
  mov    %rdx,0x30(%rsp)
  mov    %ecx,0x2c(%rsp)
  call   <ivector@plt>
  mov    %r12d,%esi
  mov    $0x1,%edi
  mov    %rax,0x48(%rsp)
  call   <ivector@plt>
  mov    %r12d,%esi
  mov    $0x1,%edi
  mov    %rax,0x50(%rsp)
  call   <ivector@plt>
  mov    %rax,%rbp
  test   %r12d,%r12d
  jle    L224
  mov    %r12d,%edx
  lea    -0x1(%r12),%r15d
  lea    0x4(%rax),%rdi
  xor    %esi,%esi
  shl    $0x2,%rdx
  mov    %r15,0x68(%rsp)
  call   <memset@plt>
  mov    0x2c(%rsp),%eax
  movq   $0x1,0x18(%rsp)
  movl   $0x1,0x28(%rsp)
  lea    -0x1(%rax),%ebx
  movl   $0x1,0x14(%rsp)
  lea    0x0(,%rbx,4),%rax
  add    $0x2,%rbx
  mov    %rax,0x60(%rsp)
  mov    %r15,%rax
  lea    0x2(%r15),%r15
  lea    0x8(,%rax,4),%rax
  mov    %r15,%r14
  mov    %rax,0x58(%rsp)
  nopl   0x0(%rax)
L56:
  mov    $0x1,%r15d
  mov    %rbx,0x40(%rsp)
  mov    0x20(%rsp),%r13
  movl   $0x0,(%rsp)
  mov    %r15,%rbx
  mov    %r14,0x38(%rsp)
  jmp    L167
  nopl   0x0(%rax)
L238:
  add    $0x1,%rbx
  cmp    %ebx,%r12d
  jl     L282
L167:
  cmpl   $0x1,0x0(%rbp,%rbx,4)
  mov    %ebx,0x8(%rsp)
  je     L238
  mov    $0x1,%r14d
  jmp    L88
  nopl   0x0(%rax)
L111:
  cmp    $0x1,%eax
  jle    L169
  lea    D0(%rip),%rdi
  call   <NRerror@plt>
L169:
  add    $0x1,%r14
  cmp    %r14d,%r12d
  jl     L238
L88:
  mov    0x0(%rbp,%r14,4),%eax
  lea    0x0(,%r14,4),%r15
  test   %eax,%eax
  jne    L111
  mov    0x0(%r13,%rbx,8),%rax
  movss  (%rax,%r14,4),%xmm0
  call   <fabs@plt>
  comiss (%rsp),%xmm0
  jb     L169
  mov    0x0(%r13,%rbx,8),%rax
  movss  (%rax,%r15,1),%xmm0
  call   <fabs@plt>
  mov    0x8(%rsp),%eax
  mov    %r14d,0x14(%rsp)
  movss  %xmm0,(%rsp)
  mov    %eax,0x28(%rsp)
  jmp    L169
  nopl   (%rax)
L282:
  movslq 0x14(%rsp),%r8
  mov    0x20(%rsp),%rsi
  mov    0x38(%rsp),%r14
  mov    0x40(%rsp),%rbx
  lea    0x0(,%r8,8),%rdx
  lea    0x0(,%r8,4),%r15
  lea    (%rsi,%rdx,1),%r9
  addl   $0x1,0x0(%rbp,%r15,1)
  mov    (%r9),%rsi
  cmp    0x28(%rsp),%r8d
  jne    L102
L116:
  mov    0x18(%rsp),%rax
  mov    0x50(%rsp),%rdi
  pxor   %xmm4,%xmm4
  mov    0x28(%rsp),%ecx
  mov    %ecx,(%rdi,%rax,4)
  mov    0x48(%rsp),%rdi
  mov    0x14(%rsp),%ecx
  mov    %ecx,(%rdi,%rax,4)
  lea    (%rsi,%r15,1),%rax
  movss  (%rax),%xmm1
  ucomiss %xmm4,%xmm1
  jp     L295
  je     L178
L295:
  movss  D7(%rip),%xmm0
  mov    0x58(%rsp),%rdi
  movl   $0x3f800000,(%rax)
  lea    0x4(%rsi),%rax
  divss  %xmm1,%xmm0
  add    %rsi,%rdi
  cs nopw 0x0(%rax,%rax,1)
L275:
  movss  (%rax),%xmm1
  add    $0x4,%rax
  mulss  %xmm0,%xmm1
  movss  %xmm1,-0x4(%rax)
  cmp    %rax,%rdi
  jne    L275
  mov    0x2c(%rsp),%eax
  test   %eax,%eax
  jle    L281
  mov    0x30(%rsp),%rax
  mov    0x60(%rsp),%rcx
  mov    (%rax,%r8,8),%rdi
  lea    0x4(%rdi),%rax
  lea    0x8(%rdi,%rcx,1),%rdi
  nopl   (%rax)
L241:
  movss  (%rax),%xmm1
  add    $0x4,%rax
  mulss  %xmm0,%xmm1
  movss  %xmm1,-0x4(%rax)
  cmp    %rax,%rdi
  jne    L241
L281:
  mov    0x14(%rsp),%r8d
  mov    0x20(%rsp),%r10
  mov    $0x1,%edi
  mov    0x30(%rsp),%r9
  mov    0x2c(%rsp),%r11d
  nopw   0x0(%rax,%rax,1)
L13:
  cmp    %edi,%r8d
  je     L71
  mov    (%r10,%rdi,8),%rcx
  lea    (%rcx,%r15,1),%rax
  movss  (%rax),%xmm1
  movl   $0x0,(%rax)
  mov    $0x1,%eax
  nopl   0x0(%rax)
L90:
  movss  (%rsi,%rax,4),%xmm2
  movss  (%rcx,%rax,4),%xmm0
  mulss  %xmm1,%xmm2
  subss  %xmm2,%xmm0
  movss  %xmm0,(%rcx,%rax,4)
  add    $0x1,%rax
  cmp    %rax,%r14
  jne    L90
  test   %r11d,%r11d
  jle    L71
  mov    (%r9,%rdi,8),%rcx
  mov    (%r9,%rdx,1),%r13
  mov    $0x1,%eax
  nopw   0x0(%rax,%rax,1)
L84:
  movss  0x0(%r13,%rax,4),%xmm2
  movss  (%rcx,%rax,4),%xmm0
  mulss  %xmm1,%xmm2
  subss  %xmm2,%xmm0
  movss  %xmm0,(%rcx,%rax,4)
  add    $0x1,%rax
  cmp    %rax,%rbx
  jne    L84
L71:
  add    $0x1,%rdi
  cmp    %rdi,%r14
  jne    L13
  addq   $0x1,0x18(%rsp)
  mov    0x18(%rsp),%rax
  cmp    %eax,%r12d
  jge    L56
  mov    0x20(%rsp),%rbx
  mov    0x68(%rsp),%rax
  movslq %r12d,%r8
  mov    0x48(%rsp),%r10
  mov    0x50(%rsp),%r11
  lea    0x10(%rbx,%rax,8),%r9
  jmp    L260
  nopl   0x0(%rax)
L218:
  sub    $0x1,%r8
  test   %r8d,%r8d
  jle    L224
L260:
  movslq (%r11,%r8,4),%rdx
  movslq (%r10,%r8,4),%rax
  cmp    %eax,%edx
  je     L218
  lea    0x0(,%rdx,4),%rdi
  lea    0x0(,%rax,4),%rsi
  lea    0x8(%rbx),%rdx
  nopl   0x0(%rax)
L172:
  mov    (%rdx),%rax
  add    $0x8,%rdx
  lea    (%rax,%rdi,1),%rcx
  add    %rsi,%rax
  movss  (%rcx),%xmm0
  movss  (%rax),%xmm1
  movss  %xmm1,(%rcx)
  movss  %xmm0,(%rax)
  cmp    %r9,%rdx
  jne    L172
  sub    $0x1,%r8
  test   %r8d,%r8d
  jg     L260
L224:
  mov    %r12d,%edx
  mov    %rbp,%rdi
  mov    $0x1,%esi
  call   <free_ivector@plt>
  mov    0x50(%rsp),%rdi
  mov    %r12d,%edx
  mov    $0x1,%esi
  call   <free_ivector@plt>
  mov    0x48(%rsp),%rdi
  add    $0x78,%rsp
  mov    %r12d,%edx
  pop    %rbx
  mov    $0x1,%esi
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  jmp    10c0 <free_ivector@plt>
  xchg   %ax,%ax
L102:
  mov    0x20(%rsp),%rax
  movslq 0x28(%rsp),%r10
  mov    (%rax,%r10,8),%rdi
  mov    $0x1,%eax
  nopl   0x0(%rax,%rax,1)
L154:
  movss  (%rdi,%rax,4),%xmm0
  movss  (%rsi,%rax,4),%xmm1
  movss  %xmm1,(%rdi,%rax,4)
  movss  %xmm0,(%rsi,%rax,4)
  add    $0x1,%rax
  cmp    %rax,%r14
  jne    L154
  mov    0x2c(%rsp),%ecx
  test   %ecx,%ecx
  jle    L116
  mov    0x30(%rsp),%rax
  mov    (%rax,%r10,8),%r10
  mov    (%rax,%r8,8),%rdi
  mov    $0x1,%eax
  nopl   0x0(%rax,%rax,1)
L109:
  movss  (%r10,%rax,4),%xmm0
  movss  (%rdi,%rax,4),%xmm1
  movss  %xmm1,(%r10,%rax,4)
  movss  %xmm0,(%rdi,%rax,4)
  add    $0x1,%rax
  cmp    %rax,%rbx
  jne    L109
  jmp    L116
L178:
  lea    D11(%rip),%rdi
  mov    %r9,0x38(%rsp)
  mov    %rdx,0x8(%rsp)
  mov    %r8,(%rsp)
  call   <NRerror@plt>
  mov    0x38(%rsp),%r9
  mov    0x8(%rsp),%rdx
  mov    (%rsp),%r8
  mov    (%r9),%rsi
  lea    (%rsi,%r15,1),%rax
  movss  (%rax),%xmm1
  jmp    L295
```

tool_calls: name: extract_symbolic_data
arguments: label: D0
type: string
name: extract_symbolic_data
arguments: label: D7
type: f32
name: extract_symbolic_data
arguments: label: D11
type: string

tool: D0: gaussj: Singular Matrix-1

tool: D7: 1

tool: D11: gaussj: Singular Matrix-2

assistant: ```c
void gaussj(float **A, int n, float **B, int m) {
  int *indxc, *indxr, *ipiv;
  int i, icol = 1, irow = 1, j, k, l, ll;
  float big, dum, pivinv, temp;
  indxc = ivector(1, n);
  indxr = ivector(1, n);
  ipiv = ivector(1, n);
  for (j = 1; j <= n; j++)
    ipiv[j] = 0;
  for (i = 1; i <= n; i++) {
    big = 0.0;
    for (j = 1; j <= n; j++)
      if (ipiv[j] != 1)
        for (k = 1; k <= n; k++) {
          if (ipiv[k] == 0) {
            if (fabs(A[j][k]) >= big) {
              big = fabs(A[j][k]);
              irow = j;
              icol = k;
            }
          } else if (ipiv[k] > 1)
            NRerror("gaussj: Singular Matrix-1");
        }
    ++(ipiv[icol]);
    if (irow != icol) {
      for (l = 1; l <= n; l++) {
        temp = (A[irow][l]);
        (A[irow][l]) = (A[icol][l]);
        (A[icol][l]) = temp;
      }
      for (l = 1; l <= m; l++) {
        temp = (B[irow][l]);
        (B[irow][l]) = (B[icol][l]);
        (B[icol][l]) = temp;
      }
    }
    indxr[i] = irow;
    indxc[i] = icol;
    if (A[icol][icol] == 0.0)
      NRerror("gaussj: Singular Matrix-2");
    pivinv = 1.0 / A[icol][icol];
    A[icol][icol] = 1.0;
    for (l = 1; l <= n; l++)
      A[icol][l] *= pivinv;
    for (l = 1; l <= m; l++)
      B[icol][l] *= pivinv;
    for (ll = 1; ll <= n; ll++)
      if (ll != icol) {
        dum = A[ll][icol];
        A[ll][icol] = 0.0;
        for (l = 1; l <= n; l++)
          A[ll][l] -= A[icol][l] * dum;
        for (l = 1; l <= m; l++)
          B[ll][l] -= B[icol][l] * dum;
      }
  }
  for (l = n; l >= 1; l--) {
    if (indxr[l] != indxc[l])
      for (k = 1; k <= n; k++) {
        temp = (A[k][indxr[l]]);
        (A[k][indxr[l]]) = (A[k][indxc[l]]);
        (A[k][indxc[l]]) = temp;
      };
  }
  free_ivector(ipiv, 1, n);
  free_ivector(indxr, 1, n);
  free_ivector(indxc, 1, n);
}
```