user: Demonstrate the C code that generates the assembly output below:

```asm
<linear_GRP_solver_LAG>:
  endbr64
  push   %r15
  push   %r14
  push   %r13
  push   %r12
  push   %rbp
  push   %rbx
  sub    $0xb8,%rsp
  mov    %rdi,%r13
  mov    %rsi,%r14
  movapd %xmm0,%xmm13
  movsd  %xmm0,0x18(%rsp)
  movsd  %xmm1,0x8(%rsp)
  movsd  %xmm2,0x78(%rsp)
  movsd  %xmm3,0x80(%rsp)
  movsd  %xmm4,0x28(%rsp)
  movsd  %xmm5,0x10(%rsp)
  movsd  %xmm6,0x60(%rsp)
  movsd  %xmm7,0x68(%rsp)
  mov    %fs:0x28,%rax
  mov    %rax,0xa8(%rsp)
  xor    %eax,%eax
  movsd  D4(%rip),%xmm0
  movapd %xmm0,%xmm7
  addsd  0x110(%rsp),%xmm7
  movsd  %xmm7,0x70(%rsp)
  movsd  0x110(%rsp),%xmm1
  subsd  %xmm0,%xmm1
  divsd  %xmm7,%xmm1
  movsd  %xmm1,0x30(%rsp)
  movsd  0x110(%rsp),%xmm5
  mulsd  0xf0(%rsp),%xmm5
  movapd %xmm5,%xmm0
  divsd  %xmm13,%xmm0
  pxor   %xmm1,%xmm1
  ucomisd %xmm0,%xmm1
  ja     L145
  sqrtsd %xmm0,%xmm0
  movsd  %xmm0,0x40(%rsp)
L95:
  movsd  0x110(%rsp),%xmm6
  mulsd  0xf8(%rsp),%xmm6
  divsd  0x8(%rsp),%xmm6
  pxor   %xmm0,%xmm0
  ucomisd %xmm6,%xmm0
  ja     L71
  sqrtsd %xmm6,%xmm6
L77:
  movsd  0x18(%rsp),%xmm4
  movsd  0x40(%rsp),%xmm5
  mulsd  %xmm5,%xmm4
  movsd  %xmm4,0x48(%rsp)
  movsd  0x8(%rsp),%xmm3
  mulsd  %xmm6,%xmm3
  movsd  %xmm3,0x50(%rsp)
  lea    0xa0(%rsp),%rdx
  lea    0x98(%rsp),%rsi
  lea    0x90(%rsp),%rdi
  mov    $0x1f4,%ecx
  movsd  0x118(%rsp),%xmm7
  movsd  0xf8(%rsp),%xmm4
  movsd  0xf0(%rsp),%xmm3
  movsd  0x10(%rsp),%xmm2
  movsd  0x28(%rsp),%xmm1
  movsd  0x110(%rsp),%xmm0
  call   <Riemann_solver_exact@plt>
  mov    0x98(%rsp),%r12
  movq   %r12,%xmm2
  comisd 0xf0(%rsp),%xmm2
  jbe    L119
  movsd  0x30(%rsp),%xmm7
  movapd %xmm7,%xmm0
  mulsd  0xf0(%rsp),%xmm0
  addsd  %xmm2,%xmm0
  mulsd  0x18(%rsp),%xmm0
  movq   %r12,%xmm1
  mulsd  %xmm7,%xmm1
  addsd  0xf0(%rsp),%xmm1
  divsd  %xmm1,%xmm0
  movq   %xmm0,%rbx
L46:
  movq   %r12,%xmm4
  comisd 0xf8(%rsp),%xmm4
  jbe    L74
  movsd  0x30(%rsp),%xmm2
  movapd %xmm2,%xmm1
  mulsd  0xf8(%rsp),%xmm1
  addsd  %xmm4,%xmm1
  mulsd  0x8(%rsp),%xmm1
  movq   %r12,%xmm0
  mulsd  %xmm2,%xmm0
  addsd  0xf8(%rsp),%xmm0
  divsd  %xmm0,%xmm1
  movq   %xmm1,%rbp
L22:
  mov    0x90(%rsp),%r15
  mov    %r15,0x8(%r14)
  mov    %r12,0x10(%r14)
  mov    %rbx,(%r14)
  mov    %rbp,0x18(%r14)
  movq   %r12,%xmm4
  mulsd  0x110(%rsp),%xmm4
  movq   %xmm4,%r14
  movapd %xmm4,%xmm0
  movq   %rbx,%xmm6
  divsd  %xmm6,%xmm0
  pxor   %xmm1,%xmm1
  ucomisd %xmm0,%xmm1
  ja     L168
  sqrtsd %xmm0,%xmm0
  movsd  %xmm0,0x18(%rsp)
L109:
  movq   %r14,%xmm0
  movq   %rbp,%xmm7
  divsd  %xmm7,%xmm0
  pxor   %xmm1,%xmm1
  ucomisd %xmm0,%xmm1
  ja     L129
  sqrtsd %xmm0,%xmm0
  movsd  %xmm0,0x20(%rsp)
L31:
  movq   %rbp,%xmm4
  mulsd  0x20(%rsp),%xmm4
  movsd  %xmm4,0x58(%rsp)
  movq   %rbx,%xmm7
  mulsd  0x18(%rsp),%xmm7
  movsd  %xmm7,0x38(%rsp)
  movsd  0x28(%rsp),%xmm0
  subsd  0x10(%rsp),%xmm0
  movsd  0xf0(%rsp),%xmm1
  subsd  0xf8(%rsp),%xmm1
  mulsd  %xmm0,%xmm0
  mulsd  %xmm1,%xmm1
  addsd  %xmm1,%xmm0
  pxor   %xmm1,%xmm1
  ucomisd %xmm0,%xmm1
  ja     L167
  sqrtsd %xmm0,%xmm0
L106:
  movsd  0x118(%rsp),%xmm1
  mulsd  D72(%rip),%xmm1
  comisd %xmm0,%xmm1
  jbe    L162
  movsd  D4(%rip),%xmm2
  movapd %xmm2,%xmm6
  divsd  0x38(%rsp),%xmm6
  movsd  %xmm6,0x28(%rsp)
  movq   D12(%rip),%xmm1
  movapd %xmm6,%xmm0
  xorpd  %xmm1,%xmm0
  movsd  0x48(%rsp),%xmm3
  mulsd  %xmm3,%xmm0
  mulsd  0x60(%rsp),%xmm3
  addsd  0x100(%rsp),%xmm3
  mulsd  %xmm3,%xmm0
  divsd  0x58(%rsp),%xmm2
  xorpd  %xmm2,%xmm1
  movsd  0x50(%rsp),%xmm3
  mulsd  %xmm3,%xmm1
  mulsd  0x68(%rsp),%xmm3
  subsd  0x108(%rsp),%xmm3
  mulsd  %xmm3,%xmm1
  movsd  D109(%rip),%xmm11
L138:
  movsd  0x28(%rsp),%xmm7
  movapd %xmm7,%xmm5
  mulsd  %xmm11,%xmm5
  movapd %xmm0,%xmm3
  mulsd  %xmm2,%xmm3
  mulsd  %xmm1,%xmm7
  subsd  %xmm7,%xmm3
  movapd %xmm2,%xmm4
  subsd  %xmm5,%xmm4
  divsd  %xmm4,%xmm3
  movsd  %xmm3,0x8(%r13)
  mulsd  %xmm11,%xmm0
  subsd  %xmm1,%xmm0
  subsd  %xmm2,%xmm5
  divsd  %xmm5,%xmm0
  movsd  %xmm0,0x10(%r13)
  movsd  D4(%rip),%xmm1
  movapd %xmm1,%xmm2
  movsd  0x18(%rsp),%xmm6
  divsd  %xmm6,%xmm2
  divsd  %xmm6,%xmm2
  mulsd  %xmm0,%xmm2
  movsd  %xmm2,0x0(%r13)
  movsd  0x20(%rsp),%xmm5
  divsd  %xmm5,%xmm1
  divsd  %xmm5,%xmm1
  mulsd  %xmm0,%xmm1
  movsd  %xmm1,0x18(%r13)
  mov    0xa8(%rsp),%rax
  sub    %fs:0x28,%rax
  jne    L137
  add    $0xb8,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
L145:
  call   <sqrt@plt>
  movsd  %xmm0,0x40(%rsp)
  jmp    L95
L71:
  movapd %xmm6,%xmm0
  call   <sqrt@plt>
  movapd %xmm0,%xmm6
  jmp    L77
L119:
  movq   %r12,%xmm0
  divsd  0xf0(%rsp),%xmm0
  movsd  D4(%rip),%xmm1
  divsd  0x110(%rsp),%xmm1
  call   <pow@plt>
  mulsd  0x18(%rsp),%xmm0
  movq   %xmm0,%rbx
  jmp    L46
L74:
  movq   %r12,%xmm0
  divsd  0xf8(%rsp),%xmm0
  movsd  D4(%rip),%xmm1
  divsd  0x110(%rsp),%xmm1
  call   <pow@plt>
  mulsd  0x8(%rsp),%xmm0
  movq   %xmm0,%rbp
  jmp    L22
L168:
  call   <sqrt@plt>
  movsd  %xmm0,0x18(%rsp)
  jmp    L109
L129:
  call   <sqrt@plt>
  movsd  %xmm0,0x20(%rsp)
  jmp    L31
L167:
  call   <sqrt@plt>
  jmp    L106
L162:
  movsd  D4(%rip),%xmm5
  movsd  0x38(%rsp),%xmm7
  divsd  %xmm7,%xmm5
  movsd  %xmm5,0x28(%rsp)
  movsd  0x110(%rsp),%xmm0
  mulsd  D58(%rip),%xmm0
  subsd  D4(%rip),%xmm0
  movsd  %xmm0,0x88(%rsp)
  movapd %xmm0,%xmm1
  mulsd  D97(%rip),%xmm1
  movapd %xmm7,%xmm0
  divsd  0x48(%rsp),%xmm0
  divsd  0x70(%rsp),%xmm1
  call   <pow@plt>
  movsd  0x40(%rsp),%xmm2
  mulsd  %xmm2,%xmm2
  movapd %xmm2,%xmm1
  mulsd  0x78(%rsp),%xmm1
  subsd  0x100(%rsp),%xmm1
  movsd  D4(%rip),%xmm2
  movsd  0x48(%rsp),%xmm6
  divsd  %xmm6,%xmm2
  divsd  0x88(%rsp),%xmm2
  mulsd  %xmm2,%xmm1
  subsd  D4(%rip),%xmm0
  mulsd  %xmm0,%xmm1
  movsd  0x100(%rsp),%xmm5
  divsd  %xmm6,%xmm5
  movapd %xmm5,%xmm0
  addsd  0x60(%rsp),%xmm0
  addsd  %xmm0,%xmm1
  movsd  %xmm1,0x40(%rsp)
  mulsd  0x38(%rsp),%xmm6
  movapd %xmm6,%xmm0
  pxor   %xmm1,%xmm1
  ucomisd %xmm6,%xmm1
  ja     L59
  sqrtsd %xmm0,%xmm0
L47:
  movq   D12(%rip),%xmm5
  xorpd  %xmm5,%xmm0
  mulsd  0x40(%rsp),%xmm0
  movq   %r12,%xmm7
  subsd  0xf8(%rsp),%xmm7
  movq   %r15,%xmm9
  subsd  0x10(%rsp),%xmm9
  movapd %xmm7,%xmm10
  divsd  %xmm9,%xmm10
  movsd  D97(%rip),%xmm8
  movapd %xmm7,%xmm12
  mulsd  %xmm8,%xmm12
  movsd  0x30(%rsp),%xmm14
  movapd %xmm14,%xmm3
  mulsd  0xf8(%rsp),%xmm3
  movq   %r12,%xmm6
  addsd  %xmm6,%xmm3
  movapd %xmm12,%xmm1
  divsd  %xmm3,%xmm1
  movsd  D93(%rip),%xmm6
  movapd %xmm6,%xmm11
  subsd  %xmm1,%xmm11
  movapd %xmm10,%xmm2
  xorpd  %xmm5,%xmm2
  movsd  0x58(%rsp),%xmm4
  divsd  %xmm4,%xmm2
  divsd  %xmm4,%xmm2
  movsd  D4(%rip),%xmm4
  movapd %xmm11,%xmm1
  subsd  %xmm4,%xmm1
  divsd  %xmm10,%xmm1
  subsd  %xmm1,%xmm2
  movapd %xmm14,%xmm1
  mulsd  %xmm7,%xmm1
  mulsd  %xmm8,%xmm1
  divsd  %xmm3,%xmm1
  addsd  %xmm4,%xmm1
  movsd  0x50(%rsp),%xmm4
  mulsd  %xmm4,%xmm4
  divsd  %xmm10,%xmm4
  mulsd  %xmm4,%xmm1
  movsd  0x8(%rsp),%xmm15
  mulsd  %xmm15,%xmm9
  mulsd  %xmm8,%xmm9
  addsd  %xmm10,%xmm9
  addsd  %xmm9,%xmm1
  mulsd  0x68(%rsp),%xmm1
  movapd %xmm14,%xmm4
  mulsd  %xmm8,%xmm4
  mulsd  %xmm7,%xmm4
  divsd  %xmm3,%xmm4
  movapd %xmm4,%xmm3
  addsd  %xmm6,%xmm3
  mulsd  0x108(%rsp),%xmm3
  subsd  %xmm3,%xmm1
  xorpd  %xmm5,%xmm12
  divsd  %xmm15,%xmm12
  mulsd  0x80(%rsp),%xmm12
  addsd  %xmm12,%xmm1
  jmp    L138
L59:
  call   <sqrt@plt>
  jmp    L47
L137:
  call   <__stack_chk_fail@plt>
```

tool_calls: name: extract_symbolic_data
arguments: label: D4
type_guessed: f64
name: extract_symbolic_data
arguments: label: D72
type_guessed: byte[8]
name: extract_symbolic_data
arguments: label: D12
type_guessed: byte[8]
name: extract_symbolic_data
arguments: label: D109
type_guessed: byte[8]
name: extract_symbolic_data
arguments: label: D58
type_guessed: f64
name: extract_symbolic_data
arguments: label: D97
type_guessed: f64
name: extract_symbolic_data
arguments: label: D93
type_guessed: f64

tool: D4: 1

tool: {"D72": [0x00,0x00,0x00,0x00,0x00,0x40,0x8f,0x40]}

tool: {"D12": [0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x80]}

tool: {"D109": [0x00,0x00,0x00,0x00,0x00,0x00,0xf0,0xbf]}

tool: D58: 3

tool: D97: 0.5

tool: D93: 2

assistant: ```c
void linear_GRP_solver_LAG(double *direvative, double *mid, double rho_L, double rho_R, double s_rho_L, double s_rho_R, double u_L, double u_R, double s_u_L, double s_u_R, double p_L, double p_R, double s_p_L, double s_p_R, double gamma, double eps) {
  double dist;
  int CRW[2];
  double c_L, c_R, g_L, g_R;
  double rho_star, u_star, p_star, rho_star_L, rho_star_R, c_star_L, c_star_R, g_star_L, g_star_R, beta_star;
  double a_L, b_L, d_L, a_R, b_R, d_R, L_rho, L_u, L_p, sigma;
  double zeta = (gamma - 1.0) / (gamma + 1.0);
  c_L = sqrt(gamma * p_L / rho_L);
  c_R = sqrt(gamma * p_R / rho_R);
  g_L = rho_L * c_L;
  g_R = rho_R * c_R;
  Riemann_solver_exact(&u_star, &p_star, gamma, u_L, u_R, p_L, p_R, c_L, c_R, CRW, eps, 500);
  if (p_star > p_L)
    rho_star_L = rho_L * (p_star + zeta * p_L) / (p_L + zeta * p_star);
  else
    rho_star_L = rho_L * pow(p_star / p_L, 1.0 / gamma);
  if (p_star > p_R)
    rho_star_R = rho_R * (p_star + zeta * p_R) / (p_R + zeta * p_star);
  else
    rho_star_R = rho_R * pow(p_star / p_R, 1.0 / gamma);
  mid[1] = u_star;
  mid[2] = p_star;
  mid[0] = rho_star_L;
  mid[3] = rho_star_R;
  c_star_L = sqrt(gamma * p_star / rho_star_L);
  c_star_R = sqrt(gamma * p_star / rho_star_R);
  g_star_R = rho_star_R * c_star_R;
  g_star_L = rho_star_L * c_star_L;
  beta_star = g_star_L / g_L;
  dist = sqrt((u_L - u_R) * (u_L - u_R) + (p_L - p_R) * (p_L - p_R));
  if ((dist < eps * 1000)) {
    a_L = 1.0;
    b_L = 1.0 / g_star_L;
    d_L = -1.0 / g_star_L * g_L * (g_L * s_u_L + s_p_L);
    a_R = -1.0;
    b_R = 1.0 / g_star_R;
    d_R = -1.0 / g_star_R * g_R * (g_R * s_u_R - s_p_R);
  } else {
    a_L = 1.0;
    b_L = 1.0 / g_star_L;
    d_L = (s_u_L + s_p_L / g_L) + 1.0 / g_L / (3.0 * gamma - 1.0) * (c_L * c_L * s_rho_L - s_p_L) * (pow(beta_star, (3.0 * gamma - 1.0) / 2.0 / (gamma + 1.0)) - 1.0);
    d_L = (-1.0) * sqrt(g_L * g_star_L) * d_L;
    sigma = (p_star - p_R) / (u_star - u_R);
    a_R = 2.0 - 0.5 * (p_star - p_R) / (p_star + zeta * p_R);
    b_R = -1.0 * sigma / g_star_R / g_star_R - (a_R - 1) / sigma;
    L_rho = -1.0 * (p_star - p_R) / 2.0 / rho_R;
    L_u = sigma + rho_R * (u_star - u_R) / 2.0 + g_R * g_R / sigma * (1.0 + zeta * (p_star - p_R) / 2.0 / (p_star + zeta * p_R));
    L_p = -1.0 * (2 + zeta / 2.0 * (p_star - p_R) / (p_star + zeta * p_R));
    d_R = L_u * s_u_R + L_p * s_p_R + L_rho * s_rho_R;
  }
  direvative[1] = (d_L * b_R - d_R * b_L) / (a_L * b_R - a_R * b_L);
  direvative[2] = (d_L * a_R - d_R * a_L) / (b_L * a_R - b_R * a_L);
  direvative[0] = 1.0 / c_star_L / c_star_L * direvative[2];
  direvative[3] = 1.0 / c_star_R / c_star_R * direvative[2];
}
```