user: Please show me the C implementation that matches the assembly code listed below:

```asm
<main>:
  endbr64
  push   %r14
  xor    %eax,%eax
  push   %rbp
  push   %rbx
  mov    %rsi,%rbx
  add    $0xffffffffffffff80,%rsp
  mov    0x8(%rsi),%rdi
  call   <atof@plt>
  mov    0x10(%rbx),%rdi
  mov    %eax,%ebp
  xor    %eax,%eax
  call   <atof@plt>
  pxor   %xmm1,%xmm1
  pxor   %xmm6,%xmm6
  cvtsi2sd %ebp,%xmm1
  divsd  D19(%rip),%xmm1
  addsd  D110(%rip),%xmm1
  cvtsi2sd %eax,%xmm6
  movapd %xmm1,%xmm0
  subsd  D110(%rip),%xmm0
  movq   %xmm6,%rbx
  pxor   %xmm6,%xmm6
  divsd  %xmm1,%xmm0
  ucomisd %xmm0,%xmm6
  ja     L11
  sqrtsd %xmm0,%xmm0
L7:
  mulsd  D196(%rip),%xmm0
  lea    D231(%rip),%rbp
  call   <abs@plt>
  pxor   %xmm8,%xmm8
  movq   %rbx,%xmm0
  lea    0x78(%rsp),%rdi
  cvtsi2sd %eax,%xmm8
  lea    0x70(%rsp),%rsi
  mov    $0x64,%ebx
  mulsd  D223(%rip),%xmm0
  divsd  D38(%rip),%xmm0
  movsd  %xmm8,(%rsp)
  call   <sincos@plt>
  movsd  (%rsp),%xmm8
  movsd  0x78(%rsp),%xmm6
  pxor   %xmm3,%xmm3
  movsd  D239(%rip),%xmm7
  movsd  %xmm3,(%rsp)
  movapd %xmm3,%xmm2
  movapd %xmm3,%xmm5
  mulsd  %xmm8,%xmm6
  movq   %xmm3,%r14
  mulsd  0x70(%rsp),%xmm8
  nopl   0x0(%rax)
L21:
  movapd %xmm7,%xmm12
  movapd %xmm2,%xmm11
  movsd  %xmm5,0x58(%rsp)
  movsd  (%rsp),%xmm5
  mulsd  %xmm7,%xmm12
  mov    D178(%rip),%rax
  movsd  %xmm8,0x68(%rsp)
  mulsd  %xmm2,%xmm11
  movapd %xmm5,%xmm4
  movsd  %xmm6,0x60(%rsp)
  mulsd  %xmm5,%xmm4
  movq   %rax,%xmm1
  movsd  %xmm7,0x20(%rsp)
  movsd  %xmm2,0x10(%rsp)
  movapd %xmm12,%xmm0
  movsd  %xmm12,0x50(%rsp)
  addsd  %xmm11,%xmm0
  movsd  %xmm11,0x48(%rsp)
  movsd  %xmm4,0x40(%rsp)
  addsd  %xmm4,%xmm0
  call   <pow@plt>
  movsd  0x20(%rsp),%xmm7
  movsd  (%rsp),%xmm3
  mov    %rbp,%rdi
  movsd  0x10(%rsp),%xmm2
  movsd  %xmm0,0x8(%rsp)
  mov    $0x4,%eax
  movq   %r14,%xmm0
  movapd %xmm7,%xmm1
  movsd  %xmm7,0x10(%rsp)
  movsd  %xmm2,0x18(%rsp)
  call   <printf@plt>
  mov    D243(%rip),%rax
  movsd  0x8(%rsp),%xmm13
  movq   %rax,%xmm1
  movapd %xmm13,%xmm0
  call   <pow@plt>
  movsd  0x10(%rsp),%xmm7
  movsd  D104(%rip),%xmm10
  mov    D243(%rip),%rax
  movsd  0x8(%rsp),%xmm13
  mulsd  %xmm7,%xmm10
  movsd  %xmm7,0x38(%rsp)
  mulsd  (%rsp),%xmm10
  movq   %rax,%xmm1
  mulsd  D17(%rip),%xmm10
  divsd  %xmm0,%xmm10
  movapd %xmm13,%xmm0
  movsd  %xmm10,0x30(%rsp)
  call   <pow@plt>
  movsd  0x18(%rsp),%xmm2
  movsd  D104(%rip),%xmm3
  mov    D243(%rip),%rax
  movsd  0x8(%rsp),%xmm13
  mulsd  %xmm2,%xmm3
  movsd  %xmm2,0x28(%rsp)
  mulsd  (%rsp),%xmm3
  movq   %rax,%xmm1
  mulsd  D17(%rip),%xmm3
  divsd  %xmm0,%xmm3
  movapd %xmm13,%xmm0
  movsd  %xmm3,0x20(%rsp)
  call   <pow@plt>
  movsd  0x40(%rsp),%xmm4
  movsd  0x50(%rsp),%xmm12
  movapd %xmm0,%xmm1
  movsd  0x48(%rsp),%xmm11
  movsd  0x58(%rsp),%xmm5
  addsd  %xmm4,%xmm4
  movsd  0x60(%rsp),%xmm6
  mov    D178(%rip),%rax
  movsd  %xmm5,0x10(%rsp)
  movsd  %xmm6,0x8(%rsp)
  movapd %xmm4,%xmm0
  subsd  %xmm12,%xmm0
  subsd  %xmm11,%xmm0
  mulsd  D17(%rip),%xmm0
  movsd  D110(%rip),%xmm11
  divsd  %xmm1,%xmm0
  movapd %xmm5,%xmm1
  mulsd  %xmm5,%xmm1
  movsd  %xmm0,0x18(%rsp)
  movapd %xmm6,%xmm0
  mulsd  %xmm6,%xmm0
  addsd  %xmm1,%xmm0
  addsd  %xmm1,%xmm0
  divsd  D84(%rip),%xmm0
  subsd  %xmm0,%xmm11
  movq   %rax,%xmm1
  movapd %xmm11,%xmm0
  call   <pow@plt>
  movsd  0x10(%rsp),%xmm5
  movsd  0x68(%rsp),%xmm8
  movsd  D110(%rip),%xmm12
  movsd  0x8(%rsp),%xmm6
  movsd  0x30(%rsp),%xmm10
  movsd  0x20(%rsp),%xmm3
  movapd %xmm5,%xmm1
  movapd %xmm5,%xmm11
  divsd  %xmm0,%xmm12
  movsd  0x18(%rsp),%xmm4
  movapd %xmm8,%xmm0
  mov    D240(%rip),%rax
  mulsd  %xmm3,%xmm1
  movsd  0x38(%rsp),%xmm7
  movsd  0x28(%rsp),%xmm2
  mulsd  %xmm10,%xmm0
  movq   %rax,%xmm9
  mulsd  %xmm6,%xmm10
  mulsd  %xmm4,%xmm11
  mulsd  %xmm8,%xmm3
  subsd  %xmm10,%xmm1
  mulsd  D179(%rip),%xmm1
  subsd  %xmm0,%xmm11
  movapd %xmm4,%xmm0
  movsd  (%rsp),%xmm4
  mulsd  %xmm6,%xmm0
  mulsd  D93(%rip),%xmm11
  subsd  %xmm3,%xmm0
  mulsd  D179(%rip),%xmm0
  movsd  D110(%rip),%xmm3
  divsd  %xmm9,%xmm1
  divsd  %xmm9,%xmm11
  mulsd  %xmm12,%xmm1
  addsd  %xmm1,%xmm8
  addsd  %xmm8,%xmm4
  movsd  %xmm4,(%rsp)
  movq   %r14,%xmm4
  divsd  %xmm9,%xmm0
  addsd  %xmm4,%xmm3
  movq   %xmm3,%r14
  mulsd  %xmm12,%xmm11
  addsd  %xmm11,%xmm6
  addsd  %xmm6,%xmm2
  mulsd  %xmm12,%xmm0
  addsd  %xmm0,%xmm5
  addsd  %xmm5,%xmm7
  sub    $0x1,%ebx
  jne    L21
  sub    $0xffffffffffffff80,%rsp
  xor    %eax,%eax
  pop    %rbx
  pop    %rbp
  pop    %r14
  ret
L11:
  call   <sqrt@plt>
  jmp    L7
```

tool_calls: name: extract_symbolic_data
arguments: label: D19
data_type: byte[8]
name: extract_symbolic_data
arguments: label: D110
data_type: f64
name: extract_symbolic_data
arguments: label: D196
data_type: f64
name: extract_symbolic_data
arguments: label: D231
data_type: string
name: extract_symbolic_data
arguments: label: D223
data_type: f64
name: extract_symbolic_data
arguments: label: D38
data_type: f32[2]
name: extract_symbolic_data
arguments: label: D239
data_type: byte[8]
name: extract_symbolic_data
arguments: label: D178
data_type: f64
name: extract_symbolic_data
arguments: label: D243
data_type: f32[2]
name: extract_symbolic_data
arguments: label: D104
data_type: f64
name: extract_symbolic_data
arguments: label: D17
data_type: byte[8]
name: extract_symbolic_data
arguments: label: D84
data_type: byte[8]
name: extract_symbolic_data
arguments: label: D240
data_type: byte[8]
name: extract_symbolic_data
arguments: label: D179
data_type: byte[8]
name: extract_symbolic_data
arguments: label: D93
data_type: byte[8]

tool: {"D19": [0x41,0x8c,0xbf,0x5a,0x2f,0xa9,0xe4,0x3d]}

tool: D110: 1

tool: D196: 299792458

tool: D231: %g,%g,%g,%g


tool: D223: 3.14159265

tool: {"D38": [0.0,3.6015625]}

tool: {"D239": [0x9a,0x99,0x99,0x99,0x19,0xea,0xc8,0x40]}

tool: D178: 0.5

tool: {"D243": [0.0,2.3125]}

tool: D104: 3

tool: {"D17": [0xd5,0xb7,0xef,0xf1,0x6d,0xb1,0x5d,0xc1]}

tool: {"D84": [0x56,0xd7,0xcd,0xea,0xd4,0xf4,0x73,0x43]}

tool: {"D240": [0x58,0x35,0x71,0x18,0x99,0x90,0x60,0x3a]}

tool: {"D179": [0x2d,0x90,0x13,0x8c,0xd9,0xa4,0x07,0x3c]}

tool: {"D93": [0x2d,0x90,0x13,0x8c,0xd9,0xa4,0x07,0xbc]}

assistant: ```c
int main(int argc, char **argv) {
  double k = atof(argv[1]);
  double alp = atof(argv[2]);
  double p[3];
  double b0 = 3.0e-5;
  double rt = 6378.1;
  double B[3];
  double v[3];
  double q = 1.602176e-19;
  double c = 299792458.0;
  double m = 1.672621e-27;
  double c2 = pow(c, 2.0);
  double g = (k / (m * c2)) + 1;
  double v0 = abs(c * sqrt((g - 1) / g));
  p[0] = 2 * rt;
  p[1] = 0;
  p[2] = 0;
  v[0] = 0;
  v[1] = v0 * sin(alp * 3.14159265 / 180);
  v[2] = v0 * cos(alp * 3.14159265 / 180);
  double ax = 0.0;
  double ay = 0.0;
  double az = 0.0;
  double t = 0.0;
  double maxt = 100.0;
  double deltat = 1.0;
  void f(double x, double y, double z) {
    double r = pow(pow(x, 2.) + pow(y, 2.) + pow(z, 2.), 0.5);
    double vec[3] = {3.0 * x * z, 3.0 * y * z, 2 * pow(z, 2.) - pow(x, 2.) - pow(y, 2.)};
    printf("%g,%g,%g,%g\n", t, p[0], p[1], p[2]);
    B[0] = -b0 * pow(rt, 3.) * vec[0] / pow(r, 5.);
    B[1] = -b0 * pow(rt, 3.) * vec[1] / pow(r, 5.);
    B[2] = -b0 * pow(rt, 3.) * vec[2] / pow(r, 5.);
    double v2 = pow(v[0], 2.) + pow(v[1], 2.) + pow(v[0], 2.);
    g = 1.0 / pow(1.0 - (v2 / c2), 0.5);
    ax = q * (v[1] * B[2] - v[2] * B[1]) / m * g;
    ay = -q * (v[0] * B[2] - v[2] * B[0]) / m * g;
    az = q * (v[0] * B[1] - v[1] * B[0]) / m * g;
  }
  while (t < maxt) {
    f(p[0], p[1], p[2]);
    v[0] += ax;
    v[1] += ay;
    v[2] += az;
    p[0] += v[0];
    p[1] += v[1];
    p[2] += v[2];
    t += deltat;
  }
  return 0;
}
```