user: What would the C source look like for the assembly code provided here:

```asm
<find_dadjFast>:
  endbr64
  push   %rbp
  push   %rbx
  mov    %edi,%r9d
  mov    %esi,%eax
  mov    %rcx,%rbp
  movapd %xmm0,%xmm2
  mov    %edi,%r10d
  mov    (%r8,%r10,8),%rsi
  mov    %eax,%ecx
  lea    0x0(,%rcx,4),%rdi
  pxor   %xmm3,%xmm3
  cvtss2sd (%rsi,%rcx,4),%xmm3
  mov    0x0(%rbp,%r10,8),%rbx
  pxor   %xmm5,%xmm5
  cvtss2sd (%rbx,%rcx,4),%xmm5
  movapd %xmm5,%xmm0
  subsd  %xmm3,%xmm0
  movsd  D59(%rip),%xmm4
  comisd %xmm4,%xmm0
  jbe    L62
L158:
  cmp    $0x1,%r9d
  je     L130
L40:
  movapd %xmm2,%xmm6
  mulsd  D37(%rip),%xmm6
  pxor   %xmm1,%xmm1
  comisd %xmm6,%xmm1
  ja     L16
  movapd %xmm6,%xmm10
  xorpd  D40(%rip),%xmm10
  comisd %xmm6,%xmm1
  ja     L95
  movapd %xmm1,%xmm10
L13:
  movapd %xmm10,%xmm1
  jmp    L95
L62:
  movsd  D95(%rip),%xmm4
  movapd %xmm0,%xmm1
  cmpltsd %xmm4,%xmm1
  andpd  %xmm1,%xmm4
  andnpd %xmm0,%xmm1
  orpd   %xmm1,%xmm4
  jmp    L158
L130:
  lea    -0x1(%rdx),%r9d
  cmp    %eax,%r9d
  jbe    L40
  test   %eax,%eax
  je     L40
  mov    0x0(%rbp),%r11
  mov    (%r8),%r10
  mov    0x10(%rbp),%r9
  mov    0x10(%r8),%r8
  pxor   %xmm1,%xmm1
  cvtss2sd -0x4(%rbx,%rdi,1),%xmm1
  pxor   %xmm6,%xmm6
  cvtss2sd -0x4(%rsi,%rdi,1),%xmm6
  subsd  %xmm6,%xmm1
  addsd  %xmm5,%xmm1
  subsd  %xmm3,%xmm1
  pxor   %xmm5,%xmm5
  cvtss2sd 0x4(%rbx,%rdi,1),%xmm5
  addsd  %xmm5,%xmm1
  pxor   %xmm5,%xmm5
  cvtss2sd 0x4(%rsi,%rdi,1),%xmm5
  subsd  %xmm5,%xmm1
  pxor   %xmm5,%xmm5
  cvtss2sd -0x4(%r11,%rdi,1),%xmm5
  addsd  %xmm5,%xmm1
  pxor   %xmm5,%xmm5
  cvtss2sd -0x4(%r10,%rdi,1),%xmm5
  subsd  %xmm5,%xmm1
  pxor   %xmm5,%xmm5
  cvtss2sd (%r11,%rcx,4),%xmm5
  addsd  %xmm5,%xmm1
  pxor   %xmm5,%xmm5
  cvtss2sd (%r10,%rcx,4),%xmm5
  subsd  %xmm5,%xmm1
  pxor   %xmm5,%xmm5
  cvtss2sd 0x4(%r11,%rdi,1),%xmm5
  addsd  %xmm5,%xmm1
  pxor   %xmm5,%xmm5
  cvtss2sd 0x4(%r10,%rdi,1),%xmm5
  subsd  %xmm5,%xmm1
  pxor   %xmm5,%xmm5
  cvtss2sd -0x4(%r9,%rdi,1),%xmm5
  addsd  %xmm5,%xmm1
  pxor   %xmm5,%xmm5
  cvtss2sd -0x4(%r8,%rdi,1),%xmm5
  subsd  %xmm5,%xmm1
  pxor   %xmm5,%xmm5
  cvtss2sd (%r9,%rcx,4),%xmm5
  addsd  %xmm5,%xmm1
  pxor   %xmm5,%xmm5
  cvtss2sd (%r8,%rcx,4),%xmm5
  subsd  %xmm5,%xmm1
  pxor   %xmm5,%xmm5
  cvtss2sd 0x4(%r9,%rdi,1),%xmm5
  addsd  %xmm5,%xmm1
  pxor   %xmm5,%xmm5
  cvtss2sd 0x4(%r8,%rdi,1),%xmm5
  subsd  %xmm5,%xmm1
  divsd  D19(%rip),%xmm1
  movapd %xmm2,%xmm6
  mulsd  D37(%rip),%xmm6
  comisd %xmm6,%xmm1
  ja     L2
  movapd %xmm6,%xmm10
  xorpd  D40(%rip),%xmm10
  comisd %xmm1,%xmm10
  ja     L28
  movapd %xmm1,%xmm10
  jmp    L13
L16:
  movapd %xmm6,%xmm10
  pxor   %xmm1,%xmm1
L95:
  pxor   %xmm7,%xmm7
  test   %eax,%eax
  je     L132
L28:
  pxor   %xmm7,%xmm7
  cvtss2sd -0x4(%rsi,%rdi,1),%xmm7
  subsd  %xmm3,%xmm7
L132:
  sub    $0x1,%edx
  pxor   %xmm5,%xmm5
  cmp    %eax,%edx
  jbe    L58
  pxor   %xmm5,%xmm5
  cvtss2sd 0x4(%rsi,%rdi,1),%xmm5
  subsd  %xmm3,%xmm5
L58:
  movapd %xmm6,%xmm9
  comisd %xmm6,%xmm7
  ja     L148
  xorpd  D40(%rip),%xmm9
  maxsd  %xmm7,%xmm9
L148:
  comisd %xmm6,%xmm5
  ja     L157
  xorpd  D40(%rip),%xmm6
  maxsd  %xmm5,%xmm6
L157:
  mulsd  %xmm2,%xmm2
  mulsd  %xmm0,%xmm0
  movapd %xmm0,%xmm8
  movapd %xmm2,%xmm3
  mulsd  D76(%rip),%xmm3
  mulsd  %xmm1,%xmm1
  mulsd  D42(%rip),%xmm2
  addsd  %xmm1,%xmm2
  divsd  %xmm2,%xmm1
  movapd %xmm1,%xmm0
  mulsd  %xmm10,%xmm0
  movsd  D107(%rip),%xmm2
  mulsd  %xmm2,%xmm0
  movapd %xmm8,%xmm10
  addsd  %xmm3,%xmm10
  movapd %xmm8,%xmm1
  divsd  %xmm10,%xmm1
  mulsd  %xmm4,%xmm1
  mulsd  %xmm2,%xmm1
  addsd  %xmm1,%xmm0
  mulsd  %xmm7,%xmm7
  addsd  %xmm3,%xmm7
  movapd %xmm3,%xmm1
  divsd  %xmm7,%xmm1
  mulsd  %xmm9,%xmm1
  mulsd  %xmm2,%xmm1
  addsd  %xmm1,%xmm0
  mulsd  %xmm5,%xmm5
  addsd  %xmm3,%xmm5
  divsd  %xmm5,%xmm3
  mulsd  %xmm6,%xmm3
  mulsd  %xmm2,%xmm3
  addsd  %xmm3,%xmm0
  pop    %rbx
  pop    %rbp
  ret
L2:
  movapd %xmm6,%xmm10
  jmp    L28
```

tool_calls: name: disassemble_data
arguments: tag: D59
datatype: f32[2]
name: disassemble_data
arguments: tag: D37
datatype: f64
name: disassemble_data
arguments: tag: D40
datatype: f32[4]
name: disassemble_data
arguments: tag: D95
datatype: f32[2]
name: disassemble_data
arguments: tag: D19
datatype: f64
name: disassemble_data
arguments: tag: D76
datatype: f64
name: disassemble_data
arguments: tag: D42
datatype: f64
name: disassemble_data
arguments: tag: D107
datatype: f32[2]

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

tool: D37: 0.33

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

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

tool: D19: 9

tool: D76: 4

tool: D42: 18

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

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;
}
```