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

```asm
<eikonalEquationEvaluate>:
  endbr64
  push   %r15
  movslq %edi,%rax
  push   %r14
  push   %r13
  push   %r12
  push   %rbp
  push   %rbx
  sub    $0x48,%rsp
  mov    %r8,(%rsp)
  mov    0x88(%rsp),%r13
  mov    0x80(%rsp),%r8
  test   %eax,%eax
  jle    L116
  mov    %rax,0x8(%rsp)
  mov    %esi,%eax
  mov    %esi,%r15d
  mov    %esi,%r14d
  shr    %eax
  movapd %xmm0,%xmm5
  movsd  D30(%rip),%xmm6
  xor    %ebp,%ebp
  lea    -0x1(%rax),%esi
  lea    -0x1(%r15),%eax
  mov    %rdx,%r10
  mov    %rcx,%rbx
  add    $0x1,%rax
  add    $0x1,%rsi
  and    $0xfffffffe,%r14d
  xor    %r12d,%r12d
  mov    %rax,0x10(%rsp)
  xorpd  D8(%rip),%xmm5
  shl    $0x4,%rsi
  pxor   %xmm4,%xmm4
  movsd  D44(%rip),%xmm7
  nopw   0x0(%rax,%rax,1)
L3:
  movsd  (%r10,%r12,8),%xmm0
  mov    (%rsp),%rax
  movsd  %xmm0,(%rax,%r12,8)
  mov    0x90(%rsp),%rax
  movsd  %xmm6,(%r9,%r12,8)
  movq   $0x0,(%r8,%r12,8)
  movsd  %xmm5,(%rax,%r12,8)
  test   %r15d,%r15d
  jle    L90
  cmp    $0x1,%r15d
  je     L57
  movslq %ebp,%rdi
  movapd %xmm4,%xmm1
  lea    (%rbx,%rdi,8),%rax
  lea    (%rax,%rsi,1),%rdx
  xchg   %ax,%ax
L144:
  movupd (%rax),%xmm3
  add    $0x10,%rax
  mulpd  %xmm3,%xmm3
  addsd  %xmm3,%xmm1
  unpckhpd %xmm3,%xmm3
  addsd  %xmm3,%xmm1
  cmp    %rdx,%rax
  jne    L144
  mov    %r14d,%eax
  cmp    %r14d,%r15d
  je     L168
L99:
  add    %ebp,%eax
  cltq
  movsd  (%rbx,%rax,8),%xmm2
  mulsd  %xmm2,%xmm2
  addsd  %xmm2,%xmm1
L168:
  ucomisd %xmm1,%xmm4
  ja     L119
  sqrtsd %xmm1,%xmm1
  movsd  %xmm1,(%r8,%r12,8)
L33:
  lea    0x0(,%rdi,8),%rax
  addsd  %xmm7,%xmm1
  lea    0x8(%rax),%rdx
  lea    0x0(%r13,%rax,1),%rcx
  lea    (%rbx,%rdx,1),%r11
  cmp    %r11,%rcx
  je     L40
  cmp    $0x1,%r15d
  je     L40
  movapd %xmm1,%xmm2
  add    %rbx,%rax
  xor    %edx,%edx
  unpcklpd %xmm2,%xmm2
  nopl   0x0(%rax,%rax,1)
L163:
  movupd (%rax,%rdx,1),%xmm0
  divpd  %xmm2,%xmm0
  movups %xmm0,(%rcx,%rdx,1)
  add    $0x10,%rdx
  cmp    %rsi,%rdx
  jne    L163
  cmp    %r14d,%r15d
  je     L96
  lea    0x0(%rbp,%r14,1),%eax
  cltq
  movsd  (%rbx,%rax,8),%xmm2
  divsd  %xmm1,%xmm2
  movsd  %xmm2,0x0(%r13,%rax,8)
L96:
  add    $0x1,%r12
  add    %r15d,%ebp
  cmp    0x8(%rsp),%r12
  jne    L3
L116:
  add    $0x48,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
  nopl   0x0(%rax,%rax,1)
L40:
  add    0x10(%rsp),%rdi
  shl    $0x3,%rdi
  jmp    L169
  nopl   0x0(%rax,%rax,1)
L166:
  add    $0x8,%rdx
L169:
  movsd  (%rbx,%rax,1),%xmm0
  divsd  %xmm1,%xmm0
  movsd  %xmm0,0x0(%r13,%rax,1)
  mov    %rdx,%rax
  cmp    %rdx,%rdi
  jne    L166
  jmp    L96
  cs nopw 0x0(%rax,%rax,1)
L90:
  movq   $0x0,(%r8,%r12,8)
  jmp    L96
L57:
  movapd %xmm4,%xmm1
  xor    %eax,%eax
  movslq %ebp,%rdi
  jmp    L99
L119:
  movapd %xmm1,%xmm0
  mov    %r8,0x80(%rsp)
  mov    %r9,0x38(%rsp)
  mov    %r10,0x30(%rsp)
  mov    %rsi,0x20(%rsp)
  mov    %rdi,0x18(%rsp)
  movsd  %xmm5,0x28(%rsp)
  call   <sqrt@plt>
  mov    D30(%rip),%rax
  mov    0x18(%rsp),%rdi
  pxor   %xmm4,%xmm4
  mov    0x20(%rsp),%rsi
  movsd  0x28(%rsp),%xmm5
  movapd %xmm0,%xmm1
  movq   %rax,%xmm6
  mov    0x30(%rsp),%r10
  mov    0x38(%rsp),%r9
  mov    0x80(%rsp),%r8
  mov    D44(%rip),%rax
  movsd  %xmm0,(%r8,%r12,8)
  movq   %rax,%xmm7
  jmp    L33
```

tool_calls: name: extract_symbolic_data
arguments: tag: D30
type_infered: f64
name: extract_symbolic_data
arguments: tag: D8
type_infered: f32[4]
name: extract_symbolic_data
arguments: tag: D44
type_infered: byte[8]

tool: D30: 1

tool: {"D8": [0.0,-0.0,0.0,0.0]}

tool: {"D44": [0x11,0xea,0x2d,0x81,0x99,0x97,0x71,0x3d]}

assistant: ```c
void eikonalEquationEvaluate(int nPoints, int nSpace, double rhs, double *u, double *grad_u, double *m, double *dm, double *H, double *dH, double *r) {
  int i, I;
  double normGradU;
  for (i = 0; i < nPoints; i++) {
    m[i] = u[i];
    dm[i] = 1.0;
    H[i] = 0.0;
    r[i] = -rhs;
    normGradU = 0.0;
    for (I = 0; I < nSpace; I++)
      normGradU += grad_u[i * nSpace + I] * grad_u[i * nSpace + I];
    normGradU = sqrt(normGradU);
    H[i] = normGradU;
    for (I = 0; I < nSpace; I++) {
      dH[i * nSpace + I] = grad_u[i * nSpace + I] / (normGradU + 1.0e-12);
    }
  }
}
```