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

```asm
<ApproxDistance>:
  endbr64
  push   %r12
  push   %rbp
  sub    $0x58,%rsp
  movsd  D76(%rip),%xmm4
  lea    0x48(%rsp),%rbp
  lea    0x40(%rsp),%r12
  movsd  %xmm1,0x30(%rsp)
  mulsd  %xmm4,%xmm2
  mov    %r12,%rsi
  mov    %rbp,%rdi
  movsd  %xmm3,0x8(%rsp)
  mulsd  %xmm4,%xmm0
  movsd  %xmm2,0x10(%rsp)
  movsd  %xmm0,0x18(%rsp)
  addsd  %xmm2,%xmm0
  mulsd  D47(%rip),%xmm0
  call   <sincos@plt>
  movsd  0x10(%rsp),%xmm2
  mov    %r12,%rsi
  movsd  0x18(%rsp),%xmm3
  movsd  0x40(%rsp),%xmm5
  movsd  D47(%rip),%xmm0
  mov    %rbp,%rdi
  subsd  %xmm2,%xmm3
  movsd  0x48(%rsp),%xmm10
  movsd  %xmm5,0x28(%rsp)
  movsd  %xmm10,0x20(%rsp)
  mulsd  %xmm3,%xmm0
  call   <sincos@plt>
  movsd  D15(%rip),%xmm2
  mov    %r12,%rsi
  movsd  0x30(%rsp),%xmm1
  movsd  0x40(%rsp),%xmm7
  movsd  0x48(%rsp),%xmm4
  mov    %rbp,%rdi
  mulsd  %xmm2,%xmm1
  movsd  %xmm7,0x18(%rsp)
  movsd  %xmm4,0x10(%rsp)
  movapd %xmm1,%xmm0
  movsd  0x8(%rsp),%xmm1
  mulsd  %xmm2,%xmm1
  subsd  %xmm1,%xmm0
  mulsd  D47(%rip),%xmm0
  call   <sincos@plt>
  movsd  0x10(%rsp),%xmm4
  pxor   %xmm3,%xmm3
  movsd  0x28(%rsp),%xmm5
  movsd  0x40(%rsp),%xmm2
  movsd  0x48(%rsp),%xmm1
  movapd %xmm4,%xmm6
  movapd %xmm5,%xmm0
  movsd  0x18(%rsp),%xmm7
  movsd  0x20(%rsp),%xmm10
  mulsd  %xmm4,%xmm6
  mulsd  %xmm5,%xmm0
  movapd %xmm7,%xmm8
  mulsd  %xmm7,%xmm8
  mulsd  %xmm2,%xmm6
  mulsd  %xmm1,%xmm0
  mulsd  %xmm2,%xmm8
  mulsd  %xmm2,%xmm6
  mulsd  %xmm1,%xmm0
  mulsd  %xmm2,%xmm8
  addsd  %xmm0,%xmm6
  movapd %xmm10,%xmm0
  mulsd  %xmm10,%xmm0
  mulsd  %xmm1,%xmm0
  mulsd  %xmm1,%xmm0
  addsd  %xmm0,%xmm8
  ucomisd %xmm8,%xmm3
  ja     L45
  movapd %xmm8,%xmm1
  sqrtsd %xmm1,%xmm1
L32:
  pxor   %xmm3,%xmm3
  ucomisd %xmm6,%xmm3
  ja     L22
  movapd %xmm6,%xmm0
  sqrtsd %xmm0,%xmm0
L23:
  movsd  %xmm4,0x30(%rsp)
  movsd  %xmm7,0x28(%rsp)
  movsd  %xmm10,0x20(%rsp)
  movsd  %xmm5,0x18(%rsp)
  movsd  %xmm8,0x10(%rsp)
  movsd  %xmm6,0x8(%rsp)
  call   <atan2@plt>
  movsd  0x8(%rsp),%xmm6
  pxor   %xmm5,%xmm5
  movsd  0x10(%rsp),%xmm8
  movsd  0x28(%rsp),%xmm7
  movsd  0x20(%rsp),%xmm10
  movapd %xmm0,%xmm3
  movapd %xmm6,%xmm1
  movsd  0x30(%rsp),%xmm4
  mulsd  %xmm8,%xmm1
  ucomisd %xmm1,%xmm5
  movsd  0x18(%rsp),%xmm5
  ja     L10
  sqrtsd %xmm1,%xmm1
L41:
  divsd  %xmm3,%xmm1
  addsd  %xmm8,%xmm8
  addsd  %xmm6,%xmm6
  movsd  D18(%rip),%xmm9
  addsd  %xmm3,%xmm3
  mulsd  D23(%rip),%xmm3
  mulsd  D17(%rip),%xmm1
  movapd %xmm1,%xmm2
  addsd  %xmm9,%xmm1
  subsd  %xmm9,%xmm2
  divsd  %xmm8,%xmm2
  movsd  D66(%rip),%xmm8
  add    $0x58,%rsp
  pop    %rbp
  pop    %r12
  mulsd  %xmm8,%xmm2
  mulsd  %xmm10,%xmm2
  mulsd  %xmm10,%xmm2
  movapd %xmm2,%xmm0
  movapd %xmm1,%xmm2
  divsd  %xmm6,%xmm2
  mulsd  %xmm7,%xmm0
  mulsd  %xmm7,%xmm0
  addsd  %xmm9,%xmm0
  mulsd  %xmm8,%xmm2
  movapd %xmm2,%xmm1
  mulsd  %xmm5,%xmm1
  mulsd  %xmm5,%xmm1
  mulsd  %xmm4,%xmm1
  mulsd  %xmm4,%xmm1
  subsd  %xmm1,%xmm0
  mulsd  %xmm3,%xmm0
  ret
L45:
  movapd %xmm8,%xmm0
  movsd  %xmm4,0x30(%rsp)
  movsd  %xmm7,0x28(%rsp)
  movsd  %xmm10,0x20(%rsp)
  movsd  %xmm5,0x18(%rsp)
  movsd  %xmm6,0x10(%rsp)
  movsd  %xmm8,0x8(%rsp)
  call   <sqrt@plt>
  movsd  0x30(%rsp),%xmm4
  movsd  0x28(%rsp),%xmm7
  movsd  0x18(%rsp),%xmm5
  movsd  0x20(%rsp),%xmm10
  movapd %xmm0,%xmm1
  movsd  0x10(%rsp),%xmm6
  movsd  0x8(%rsp),%xmm8
  jmp    L32
L10:
  movsd  %xmm0,0x18(%rsp)
  movapd %xmm1,%xmm0
  movsd  %xmm4,0x38(%rsp)
  movsd  %xmm7,0x30(%rsp)
  movsd  %xmm10,0x28(%rsp)
  movsd  %xmm5,0x20(%rsp)
  movsd  %xmm8,0x10(%rsp)
  movsd  %xmm6,0x8(%rsp)
  call   <sqrt@plt>
  movsd  0x38(%rsp),%xmm4
  movsd  0x30(%rsp),%xmm7
  movsd  0x20(%rsp),%xmm5
  movsd  0x28(%rsp),%xmm10
  movapd %xmm0,%xmm1
  movsd  0x18(%rsp),%xmm3
  movsd  0x10(%rsp),%xmm8
  movsd  0x8(%rsp),%xmm6
  jmp    L41
L22:
  movapd %xmm6,%xmm0
  movsd  %xmm4,0x38(%rsp)
  movsd  %xmm7,0x30(%rsp)
  movsd  %xmm10,0x28(%rsp)
  movsd  %xmm5,0x20(%rsp)
  movsd  %xmm8,0x18(%rsp)
  movsd  %xmm1,0x10(%rsp)
  movsd  %xmm6,0x8(%rsp)
  call   <sqrt@plt>
  movsd  0x38(%rsp),%xmm4
  movsd  0x30(%rsp),%xmm7
  movsd  0x28(%rsp),%xmm10
  movsd  0x20(%rsp),%xmm5
  movsd  0x18(%rsp),%xmm8
  movsd  0x10(%rsp),%xmm1
  movsd  0x8(%rsp),%xmm6
  jmp    L23
```

tool_calls: name: unpack_data
arguments: data_label: D76
type_infered: f64
name: unpack_data
arguments: data_label: D47
type_infered: byte[8]
name: unpack_data
arguments: data_label: D15
type_infered: byte[8]
name: unpack_data
arguments: data_label: D18
type_infered: f64
name: unpack_data
arguments: data_label: D23
type_infered: f64
name: unpack_data
arguments: data_label: D17
type_infered: byte[8]
name: unpack_data
arguments: data_label: D66
type_infered: byte[8]

tool: D76: 0.01745329252

tool: {"D47": [0x00,0x00,0x00,0x00,0x00,0x00,0xe0,0x3f]}

tool: {"D15": [0x11,0xdd,0x52,0xa2,0x46,0xdf,0x91,0xbf]}

tool: D18: 1

tool: D23: 6378.135

tool: {"D17": [0x00,0x00,0x00,0x00,0x00,0x00,0x08,0x40]}

tool: {"D66": [0x76,0x38,0x67,0xc3,0x49,0x77,0x6b,0x3f]}

assistant: ```c
double ApproxDistance(double lat1, double lon1, double lat2, double lon2) {
  double F, G, L, sing, cosl, cosf, sinl, sinf, cosg, S, C, W, R, H1, H2, D;
  lat1 = 0.01745329252 * lat1;
  lon1 = -0.01745329252 * lon1;
  lat2 = 0.01745329252 * lat2;
  lon2 = -0.01745329252 * lon2;
  F = (lat1 + lat2) / 2.0;
  G = (lat1 - lat2) / 2.0;
  L = (lon1 - lon2) / 2.0;
  sing = sin(G);
  cosl = cos(L);
  cosf = cos(F);
  sinl = sin(L);
  sinf = sin(F);
  cosg = cos(G);
  S = sing * sing * cosl * cosl + cosf * cosf * sinl * sinl;
  C = cosg * cosg * cosl * cosl + sinf * sinf * sinl * sinl;
  W = atan2(sqrt(S), sqrt(C));
  R = sqrt((S * C)) / W;
  H1 = (3 * R - 1.0) / (2.0 * C);
  H2 = (3 * R + 1.0) / (2.0 * S);
  D = 2 * W * 6378.135;
  return (D * (1 + 1.0 / 298.26 * H1 * sinf * sinf * cosg * cosg - 1.0 / 298.26 * H2 * cosf * cosf * sing * sing));
}
```