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

```asm
<find_dadjFast>:
  endbr64
  push   %rbp
  mov    %esi,%eax
  mov    %edi,%r11d
  mov    %rcx,%rsi
  push   %rbx
  mov    (%rsi,%r11,8),%rbp
  mov    %eax,%ecx
  pxor   %xmm1,%xmm1
  mov    (%r8,%r11,8),%r9
  pxor   %xmm8,%xmm8
  movapd %xmm0,%xmm3
  movsd  D31(%rip),%xmm7
  lea    0x0(,%rcx,4),%r10
  cvtss2sd 0x0(%rbp,%rcx,4),%xmm1
  movapd %xmm1,%xmm2
  cvtss2sd (%r9,%rcx,4),%xmm8
  subsd  %xmm8,%xmm2
  comisd %xmm7,%xmm2
  ja     L63
  movsd  D104(%rip),%xmm7
  maxsd  %xmm2,%xmm7
L63:
  movsd  D8(%rip),%xmm4
  sub    $0x1,%edx
  mulsd  %xmm3,%xmm4
  cmp    $0x1,%edi
  je     L59
L101:
  pxor   %xmm1,%xmm1
  movapd %xmm4,%xmm0
  comisd %xmm4,%xmm1
  jbe    L134
L7:
  test   %eax,%eax
  jne    L73
L156:
  pxor   %xmm9,%xmm9
  movapd %xmm9,%xmm5
  cmp    %edx,%eax
  jae    L50
L137:
  pxor   %xmm6,%xmm6
  comisd %xmm4,%xmm5
  cvtss2sd 0x4(%r9,%r10,1),%xmm6
  subsd  %xmm8,%xmm6
  movapd %xmm6,%xmm8
  mulsd  %xmm6,%xmm8
  ja     L8
L135:
  movapd %xmm4,%xmm10
  xorpd  D0(%rip),%xmm10
  maxsd  %xmm5,%xmm10
  movapd %xmm10,%xmm5
L157:
  comisd %xmm4,%xmm6
  ja     L55
  xorpd  D0(%rip),%xmm4
  maxsd  %xmm6,%xmm4
L55:
  mulsd  %xmm3,%xmm3
  movsd  D21(%rip),%xmm6
  pop    %rbx
  mulsd  %xmm2,%xmm2
  pop    %rbp
  mulsd  %xmm3,%xmm6
  mulsd  D63(%rip),%xmm3
  addsd  %xmm6,%xmm9
  addsd  %xmm6,%xmm8
  addsd  %xmm1,%xmm3
  divsd  %xmm3,%xmm1
  movsd  D109(%rip),%xmm3
  mulsd  %xmm1,%xmm0
  movapd %xmm2,%xmm1
  addsd  %xmm6,%xmm1
  divsd  %xmm1,%xmm2
  movapd %xmm6,%xmm1
  mulsd  %xmm3,%xmm0
  divsd  %xmm9,%xmm1
  mulsd  %xmm7,%xmm2
  mulsd  %xmm3,%xmm2
  addsd  %xmm2,%xmm0
  divsd  %xmm8,%xmm6
  mulsd  %xmm5,%xmm1
  mulsd  %xmm3,%xmm1
  addsd  %xmm1,%xmm0
  movapd %xmm6,%xmm1
  mulsd  %xmm4,%xmm1
  mulsd  %xmm3,%xmm1
  addsd  %xmm1,%xmm0
  ret
  nopl   0x0(%rax)
L134:
  movapd %xmm4,%xmm5
  xorpd  D0(%rip),%xmm0
  cmpltsd %xmm1,%xmm5
  andpd  %xmm5,%xmm0
  andnpd %xmm1,%xmm5
  orpd   %xmm5,%xmm0
  test   %eax,%eax
  je     L156
  jmp    L73
  nopw   0x0(%rax,%rax,1)
L59:
  cmp    %edx,%eax
  jae    L101
  test   %eax,%eax
  je     L101
  pxor   %xmm5,%xmm5
  pxor   %xmm0,%xmm0
  mov    (%rsi),%rbx
  mov    (%r8),%r11
  cvtss2sd -0x4(%r9,%r10,1),%xmm5
  cvtss2sd -0x4(%rbp,%r10,1),%xmm0
  subsd  %xmm5,%xmm0
  mov    0x10(%rsi),%rdi
  mov    0x10(%r8),%rsi
  addsd  %xmm1,%xmm0
  pxor   %xmm1,%xmm1
  cvtss2sd 0x4(%rbp,%r10,1),%xmm1
  subsd  %xmm8,%xmm0
  addsd  %xmm1,%xmm0
  pxor   %xmm1,%xmm1
  cvtss2sd 0x4(%r9,%r10,1),%xmm1
  subsd  %xmm1,%xmm0
  pxor   %xmm1,%xmm1
  cvtss2sd -0x4(%rbx,%r10,1),%xmm1
  addsd  %xmm1,%xmm0
  pxor   %xmm1,%xmm1
  cvtss2sd -0x4(%r11,%r10,1),%xmm1
  subsd  %xmm1,%xmm0
  pxor   %xmm1,%xmm1
  cvtss2sd (%rbx,%rcx,4),%xmm1
  addsd  %xmm1,%xmm0
  pxor   %xmm1,%xmm1
  cvtss2sd (%r11,%rcx,4),%xmm1
  subsd  %xmm1,%xmm0
  pxor   %xmm1,%xmm1
  cvtss2sd 0x4(%rbx,%r10,1),%xmm1
  addsd  %xmm1,%xmm0
  pxor   %xmm1,%xmm1
  cvtss2sd 0x4(%r11,%r10,1),%xmm1
  subsd  %xmm1,%xmm0
  pxor   %xmm1,%xmm1
  cvtss2sd -0x4(%rdi,%r10,1),%xmm1
  addsd  %xmm1,%xmm0
  pxor   %xmm1,%xmm1
  cvtss2sd -0x4(%rsi,%r10,1),%xmm1
  subsd  %xmm1,%xmm0
  pxor   %xmm1,%xmm1
  cvtss2sd (%rdi,%rcx,4),%xmm1
  addsd  %xmm1,%xmm0
  pxor   %xmm1,%xmm1
  cvtss2sd (%rsi,%rcx,4),%xmm1
  subsd  %xmm1,%xmm0
  pxor   %xmm1,%xmm1
  cvtss2sd 0x4(%rdi,%r10,1),%xmm1
  addsd  %xmm1,%xmm0
  pxor   %xmm1,%xmm1
  cvtss2sd 0x4(%rsi,%r10,1),%xmm1
  subsd  %xmm1,%xmm0
  divsd  D79(%rip),%xmm0
  comisd %xmm4,%xmm0
  movapd %xmm0,%xmm1
  mulsd  %xmm0,%xmm1
  ja     L158
  movapd %xmm4,%xmm5
  xorpd  D0(%rip),%xmm5
  comisd %xmm0,%xmm5
  jbe    L7
  movapd %xmm5,%xmm0
  nopl   0x0(%rax)
L73:
  pxor   %xmm5,%xmm5
  cvtss2sd -0x4(%r9,%r10,1),%xmm5
L79:
  subsd  %xmm8,%xmm5
  movapd %xmm5,%xmm9
  mulsd  %xmm5,%xmm9
  cmp    %edx,%eax
  jb     L137
L50:
  comisd %xmm4,%xmm5
  pxor   %xmm8,%xmm8
  movapd %xmm8,%xmm6
  jbe    L135
L8:
  movapd %xmm4,%xmm5
  jmp    L157
  nopw   0x0(%rax,%rax,1)
L158:
  movapd %xmm4,%xmm0
  jmp    L79
```

tool_calls: name: extract_symbolic_data
arguments: data_label: D31
type_guessed: f32[2]
name: extract_symbolic_data
arguments: data_label: D104
type_guessed: f32[2]
name: extract_symbolic_data
arguments: data_label: D8
type_guessed: f64
name: extract_symbolic_data
arguments: data_label: D0
type_guessed: f32[4]
name: extract_symbolic_data
arguments: data_label: D21
type_guessed: f64
name: extract_symbolic_data
arguments: data_label: D63
type_guessed: f64
name: extract_symbolic_data
arguments: data_label: D109
type_guessed: f32[2]
name: extract_symbolic_data
arguments: data_label: D79
type_guessed: f64

tool: {"D31": [0.0,1.875]}

tool: {"D104": [0.0,-1.875]}

tool: D8: 0.33

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

tool: D21: 4

tool: D63: 18

tool: {"D109": [0.0,1.625]}

tool: D79: 9

assistant: ```c
double find_dadjFast(const unsigned i, const unsigned j, const unsigned nRows, const float *obsloc[3], const float *rszloc[3], const double readNoiseAmp) {
  const double mval = (double)*(rszloc[i] + j);
  const double dval0 = (double)*(obsloc[i] + j) - mval;
  double dval0u = dval0;
  if (dval0u > 1)
    dval0u = 1;
  else if (dval0u < -1)
    dval0u = -1;
  double dval9 = 0.0;
  if (i == 1 && j < nRows - 1 && j > 0) {
    dval9 = (double)*(obsloc[i] + j - 1) - (double)*(rszloc[i] + j - 1) + (double)*(obsloc[i] + j) - (double)*(rszloc[i] + j) + (double)*(obsloc[i] + j + 1) - (double)*(rszloc[i] + j + 1) + (double)*(obsloc[i - 1] + j - 1) - (double)*(rszloc[i - 1] + j - 1) + (double)*(obsloc[i - 1] + j) - (double)*(rszloc[i - 1] + j) + (double)*(obsloc[i - 1] + j + 1) - (double)*(rszloc[i - 1] + j + 1) + (double)*(obsloc[i + 1] + j - 1) - (double)*(rszloc[i + 1] + j - 1) + (double)*(obsloc[i + 1] + j) - (double)*(rszloc[i + 1] + j) + (double)*(obsloc[i + 1] + j + 1) - (double)*(rszloc[i + 1] + j + 1);
    dval9 = dval9 / 9.0;
  }
  const double readNoiseAmpFraction = 0.33;
  double dval9u = dval9;
  if (dval9u > readNoiseAmp * readNoiseAmpFraction)
    dval9u = readNoiseAmp * readNoiseAmpFraction;
  else if (dval9u < readNoiseAmp * -readNoiseAmpFraction)
    dval9u = readNoiseAmp * -readNoiseAmpFraction;
  const double dmod1 = j > 0 ? (double)*(rszloc[i] + j - 1) - mval : 0;
  const double dmod2 = j < nRows - 1 ? (double)*(rszloc[i] + j + 1) - mval : 0;
  double dmod1u = dmod1;
  if (dmod1u > readNoiseAmp * readNoiseAmpFraction)
    dmod1u = readNoiseAmp * readNoiseAmpFraction;
  else if (dmod1u < readNoiseAmp * -readNoiseAmpFraction)
    dmod1u = readNoiseAmp * -readNoiseAmpFraction;
  double dmod2u = dmod2;
  if (dmod2u > readNoiseAmp * readNoiseAmpFraction)
    dmod2u = readNoiseAmp * readNoiseAmpFraction;
  else if (dmod2u < readNoiseAmp * -readNoiseAmpFraction)
    dmod2u = readNoiseAmp * -readNoiseAmpFraction;
  const double readNoiseAmp2 = readNoiseAmp * readNoiseAmp;
  const double w0 = dval0 * dval0 / (dval0 * dval0 + 4.0 * readNoiseAmp2);
  const double w9 = dval9 * dval9 / (dval9 * dval9 + 18.0 * readNoiseAmp2);
  const double w1 = 4 * readNoiseAmp2 / (dmod1 * dmod1 + 4.0 * readNoiseAmp2);
  const double w2 = 4 * readNoiseAmp2 / (dmod2 * dmod2 + 4.0 * readNoiseAmp2);
  return dval0u * w0 * 0.25f + dval9u * w9 * 0.25f + dmod1u * w1 * 0.25f + dmod2u * w2 * 0.25f;
}
```