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

```asm
<gauss_eliminate>:
  endbr64
  push   %r15
  push   %r14
  push   %r13
  push   %r12
  mov    %esi,%r12d
  mov    $0x8,%esi
  push   %rbp
  push   %rbx
  mov    %rdi,%rbx
  mov    %r12d,%edi
  sub    $0x28,%rsp
  call   <calloc@plt>
  mov    %rax,%r11
  test   %r12d,%r12d
  jle    L188
  lea    0x1(%r12),%r14d
  movslq %r12d,%rbp
  mov    %r11,0x18(%rsp)
  xor    %r8d,%r8d
  mov    %r14d,%r15d
  mov    %r14d,%r9d
  mov    %rbp,0x8(%rsp)
  lea    -0x2(%r12),%r10d
  and    $0xfffffffe,%r15d
  shr    %r9d
  mov    %r10d,0x14(%rsp)
  lea    0x8(,%rbp,8),%r13
  movslq %r15d,%rdi
  shl    $0x4,%r9
  add    $0x2,%r10
  shl    $0x3,%rdi
  nopl   0x0(%rax)
L87:
  lea    0x0(,%r8,8),%r11
  mov    $0x1,%ecx
  cmp    $0x1,%r12d
  jne    L128
  jmp    L168
L93:
  add    $0x1,%rcx
  cmp    %r10,%rcx
  je     L168
  nopl   (%rax)
L128:
  cmp    %r8d,%ecx
  jle    L93
  mov    (%rbx,%rcx,8),%rax
  mov    (%rbx,%r8,8),%rsi
  movsd  (%rax,%r11,1),%xmm0
  lea    0x8(%rsi),%rdx
  movapd %xmm0,%xmm2
  divsd  (%rsi,%r11,1),%xmm2
  cmp    %rdx,%rax
  je     L16
  movapd %xmm2,%xmm3
  xor    %edx,%edx
  unpcklpd %xmm3,%xmm3
  nopl   0x0(%rax)
L146:
  movupd (%rsi,%rdx,1),%xmm1
  movupd (%rax,%rdx,1),%xmm0
  mulpd  %xmm3,%xmm1
  subpd  %xmm1,%xmm0
  movups %xmm0,(%rax,%rdx,1)
  add    $0x10,%rdx
  cmp    %r9,%rdx
  jne    L146
  cmp    %r15d,%r14d
  je     L93
  mulsd  (%rsi,%rdi,1),%xmm2
  add    %rdi,%rax
  add    $0x1,%rcx
  movsd  (%rax),%xmm1
  subsd  %xmm2,%xmm1
  movsd  %xmm1,(%rax)
  cmp    %r10,%rcx
  jne    L128
L168:
  add    $0x1,%r8
  cmp    %rbp,%r8
  jne    L87
  mov    0x18(%rsp),%r11
L62:
  mov    0x8(%rsp),%rdi
  lea    0x0(,%rdi,8),%rbp
  mov    -0x8(%rbx,%rbp,1),%rax
  movsd  (%rax,%rdi,8),%xmm0
  divsd  -0x8(%rax,%rbp,1),%xmm0
  mov    0x14(%rsp),%eax
  movsd  %xmm0,-0x8(%r11,%rbp,1)
  test   %eax,%eax
  js     L110
  movslq %eax,%r8
  mov    $0x2,%r10d
  pxor   %xmm3,%xmm3
  mov    $0x1,%r14d
  lea    (%r11,%r8,8),%rdi
  nopw   0x0(%rax,%rax,1)
L122:
  cmp    %r12d,%r8d
  mov    %r14d,%r15d
  mov    (%rbx,%r8,8),%rsi
  mov    %r8d,%edx
  cmovl  %r10d,%r15d
  lea    0x0(,%r8,8),%r13
  cmp    %r8d,%r12d
  jle    L185
  mov    %r15d,%ecx
  lea    (%rsi,%r13,1),%r9
  xor    %eax,%eax
  shr    %ecx
  movapd %xmm3,%xmm1
  shl    $0x4,%rcx
  nopl   (%rax)
L97:
  movupd (%rdi,%rax,1),%xmm0
  movupd (%r9,%rax,1),%xmm4
  add    $0x10,%rax
  mulpd  %xmm4,%xmm0
  addsd  %xmm0,%xmm1
  unpckhpd %xmm0,%xmm0
  addsd  %xmm0,%xmm1
  cmp    %rax,%rcx
  jne    L97
  mov    %r15d,%eax
  and    $0xfffffffe,%eax
  add    %eax,%edx
  cmp    %eax,%r15d
  je     L100
L124:
  movslq %edx,%rax
  movsd  (%rsi,%rax,8),%xmm0
  mulsd  (%r11,%rax,8),%xmm0
  addsd  %xmm0,%xmm1
L100:
  movsd  (%rsi,%rbp,1),%xmm0
  sub    $0x1,%r8
  add    $0x1,%r10d
  sub    $0x8,%rdi
  subsd  %xmm1,%xmm0
  divsd  (%rsi,%r13,1),%xmm0
  movsd  %xmm0,0x8(%rdi)
  test   %r8d,%r8d
  jns    L122
L110:
  add    $0x28,%rsp
  mov    %r11,%rax
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
  xchg   %ax,%ax
L16:
  xor    %edx,%edx
  nopw   0x0(%rax,%rax,1)
L11:
  movsd  (%rsi,%rdx,1),%xmm1
  movsd  (%rax,%rdx,1),%xmm0
  mulsd  %xmm2,%xmm1
  subsd  %xmm1,%xmm0
  movsd  %xmm0,(%rax,%rdx,1)
  add    $0x8,%rdx
  cmp    %r13,%rdx
  jne    L11
  jmp    L93
L185:
  movapd %xmm3,%xmm1
  jmp    L124
L188:
  movslq %r12d,%rax
  mov    %rax,0x8(%rsp)
  lea    -0x2(%r12),%eax
  mov    %eax,0x14(%rsp)
  jmp    L62
```

assistant: ```c
double *gauss_eliminate(double **aug, int n) {
  int i, j, k;
  double c, sum = 0.0;
  double *res = calloc(n, sizeof(double));
  for (j = 0; j < n; j++) {
    for (i = 0; i < n; i++) {
      if (i > j) {
        c = aug[i][j] / aug[j][j];
        for (k = 0; k < n + 1; k++)
          aug[i][k] = aug[i][k] - c * aug[j][k];
      }
    }
  }
  res[n - 1] = aug[n - 1][n] / aug[n - 1][n - 1];
  for (i = n - 2; i >= 0; i--) {
    sum = 0;
    for (j = i; j < n; j++)
      sum = sum + aug[i][j] * res[j];
    res[i] = (aug[i][n] - sum) / aug[i][i];
  }
  return res;
}
```