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

```asm
<invert>:
  endbr64
  push   %r15
  push   %r14
  xor    %r14d,%r14d
  push   %r13
  mov    %r14,%r15
  mov    %rdi,%r13
  push   %r12
  mov    %rdx,%r12
  push   %rbp
  push   %rbx
  sub    $0x1d8,%rsp
  mov    %rsi,0x20(%rsp)
  lea    0x30(%rsp),%rbp
  lea    0x58(%rsp),%rbx
  mov    %rdx,0x28(%rsp)
  mov    %fs:0x28,%rax
  mov    %rax,0x1c8(%rsp)
  xor    %eax,%eax
  mov    %rbp,0x8(%rsp)
  nopl   (%rax)
L94:
  mov    %r14,%rsi
  neg    %rsi
  shl    $0x2,%rsi
  lea    0x28(%rsi),%rax
  cmp    $0x8,%eax
  jae    L99
  test   $0x4,%al
  jne    L194
  test   %eax,%eax
  je     L92
  movzbl 0x0(%r13),%eax
  mov    %al,0x0(%rbp)
L92:
  mov    %r15,%rcx
  mov    0x8(%rsp),%r8
  neg    %rcx
  test   %r14,%r14
  je     L124
  nop
L144:
  movss  0x0(%rbp,%rcx,4),%xmm0
  movss  (%r8),%xmm2
  mov    %rbp,%rax
  mulss  %xmm0,%xmm2
  jmp    L197
  nopl   0x0(%rax)
L10:
  movss  (%rax,%rcx,4),%xmm0
L197:
  mulss  %xmm2,%xmm0
  movss  (%rax),%xmm1
  add    $0x4,%rax
  subss  %xmm0,%xmm1
  movss  %xmm1,-0x4(%rax)
  cmp    %rbx,%rax
  jne    L10
  add    $0x2c,%r8
  add    $0xa,%rcx
  jne    L144
  movss  0x0(%rbp),%xmm2
  mov    %rsi,0x18(%rsp)
  movaps %xmm2,%xmm0
  movss  %xmm2,0x14(%rsp)
  call   <fabs@plt>
  movss  0x14(%rsp),%xmm2
  mov    0x18(%rsp),%rsi
  movss  D59(%rip),%xmm4
  comiss %xmm0,%xmm4
  ja     L98
  mov    0x20(%rsp),%rax
  mov    0x8(%rsp),%rcx
  movss  (%rax,%r14,4),%xmm1
  mov    0x28(%rsp),%rax
  sub    %rsi,%rcx
  movss  %xmm1,(%r12)
  nopl   0x0(%rax,%rax,1)
L32:
  movss  (%rax),%xmm0
  mulss  (%rcx),%xmm0
  add    $0x4,%rax
  add    $0x28,%rcx
  subss  %xmm0,%xmm1
  movss  %xmm1,(%r12)
  cmp    %r12,%rax
  jne    L32
L55:
  movss  D49(%rip),%xmm0
  movss  D48(%rip),%xmm3
  divss  %xmm2,%xmm0
  movss  %xmm0,0x0(%rbp)
  mulss  %xmm1,%xmm0
  comiss %xmm0,%xmm3
  jb     L69
  comiss D35(%rip),%xmm0
  jae    L161
  movss  D35(%rip),%xmm0
L161:
  add    $0x1,%r14
  add    $0x2c,%rbp
  add    $0x2c,%r13
  add    $0x4,%r12
  movss  %xmm0,-0x4(%r12)
  add    $0xa,%r15
  add    $0x28,%rbx
  cmp    $0xa,%r14
  jne    L94
L60:
  mov    0x1c8(%rsp),%rax
  sub    %fs:0x28,%rax
  jne    L266
  add    $0x1d8,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
L99:
  mov    0x0(%r13),%rcx
  lea    0x8(%rbp),%r9
  mov    %r13,%r11
  and    $0xfffffffffffffff8,%r9
  mov    %rcx,0x0(%rbp)
  mov    %eax,%ecx
  mov    -0x8(%r13,%rcx,1),%r8
  mov    %r8,-0x8(%rbp,%rcx,1)
  mov    %rbp,%rcx
  sub    %r9,%rcx
  add    %ecx,%eax
  sub    %rcx,%r11
  and    $0xfffffff8,%eax
  cmp    $0x8,%eax
  jb     L92
  and    $0xfffffff8,%eax
  xor    %ecx,%ecx
L185:
  mov    %ecx,%edx
  add    $0x8,%ecx
  mov    (%r11,%rdx,1),%rdi
  mov    %rdi,(%r9,%rdx,1)
  cmp    %eax,%ecx
  jb     L185
  mov    %r15,%rcx
  mov    0x8(%rsp),%r8
  neg    %rcx
  test   %r14,%r14
  jne    L144
L124:
  movss  0x0(%rbp),%xmm2
  mov    %rsi,0x18(%rsp)
  movaps %xmm2,%xmm0
  movss  %xmm2,0x14(%rsp)
  call   <fabs@plt>
  movss  0x14(%rsp),%xmm2
  mov    0x18(%rsp),%rsi
  movss  D59(%rip),%xmm5
  comiss %xmm0,%xmm5
  ja     L98
  mov    0x20(%rsp),%rax
  movss  (%rax),%xmm1
  movss  %xmm1,(%r12)
  jmp    L55
L69:
  movss  D48(%rip),%xmm0
  jmp    L161
L194:
  mov    0x0(%r13),%ecx
  mov    %eax,%eax
  mov    %ecx,0x0(%rbp)
  mov    -0x4(%r13,%rax,1),%ecx
  mov    %ecx,-0x4(%rbp,%rax,1)
  jmp    L92
L98:
  mov    %r14,%r13
  mov    $0xa,%eax
  mov    %rsi,%r15
  xor    %esi,%esi
  sub    %r13,%rax
  mov    0x28(%rsp),%r13
  shl    $0x2,%rax
  sub    %r15,%r13
  cmp    $0x8,%eax
  jb     L205
  lea    0x8(%r13),%rcx
  mov    %eax,%edx
  movq   $0x0,0x0(%r13)
  movq   $0x0,-0x8(%r13,%rdx,1)
  and    $0xfffffffffffffff8,%rcx
  sub    %rcx,%r13
  add    %r13d,%eax
  and    $0xfffffff8,%eax
  cmp    $0x8,%eax
  jb     L60
  and    $0xfffffff8,%eax
  xor    %edx,%edx
L44:
  mov    %edx,%edi
  add    $0x8,%edx
  mov    %rsi,(%rcx,%rdi,1)
  cmp    %eax,%edx
  jb     L44
  jmp    L60
L205:
  test   $0x4,%al
  jne    L147
  test   %eax,%eax
  je     L60
  movb   $0x0,0x0(%r13)
  jmp    L60
L147:
  mov    %eax,%eax
  movl   $0x0,0x0(%r13)
  movl   $0x0,-0x4(%r13,%rax,1)
  jmp    L60
L266:
  call   <__stack_chk_fail@plt>
```

tool_calls: name: read_data
arguments: data_tag: D59
type_guessed: byte[4]
name: read_data
arguments: data_tag: D49
type_guessed: byte[4]
name: read_data
arguments: data_tag: D48
type_guessed: byte[4]
name: read_data
arguments: data_tag: D35
type_guessed: byte[4]

tool: {"D59": [0xff,0xe6,0xdb,0x2e]}

tool: {"D49": [0x00,0x00,0x80,0x3f]}

tool: {"D48": [0x77,0xbe,0x7f,0x3f]}

tool: {"D35": [0x77,0xbe,0x7f,0xbf]}

assistant: ```c
void invert(float *phi, float *psi, float *rc) {
  float save;
  long i, j, k;
  float v[100];
  for (j = 0; j < 10; ++j) {
    for (i = j; i < 10; ++i) {
      v[i + j * 10] = phi[i + j * 10];
    }
    for (k = 0; k < j; ++k) {
      save = v[j + k * 10] * v[k + k * 10];
      for (i = j; i < 10; ++i) {
        v[i + j * 10] -= v[i + k * 10] * save;
      }
    }
    if ((fabs(v[j + j * 10])) < 1e-10f) {
      goto L100;
    }
    rc[j] = psi[j];
    for (k = 0; k < j; ++k) {
      rc[j] -= rc[k] * v[j + k * 10];
    }
    v[j + j * 10] = 1.f / v[j + j * 10];
    rc[j] *= v[j + j * 10];
    rc[j] = ((((rc[j]) <= (.999f) ? (rc[j]) : (.999f))) >= (-.999f) ? (((rc[j]) <= (.999f) ? (rc[j]) : (.999f))) : (-.999f));
  }
  return;
L100:
  for (i = j; i < 10; ++i) {
    rc[i] = 0.f;
  }
}
```