user: What is the c source code of the assembly code below?

```asm
<ROE_2D_solver>:
  endbr64
  push   %r14
  push   %rbp
  push   %rbx
  sub    $0x140,%rsp
  mov    %rdi,%rbx
  movapd %xmm1,%xmm15
  movsd  %xmm1,0x50(%rsp)
  movapd %xmm2,%xmm12
  movsd  %xmm2,0x8(%rsp)
  movapd %xmm3,%xmm13
  movsd  %xmm3,0x28(%rsp)
  movapd %xmm4,%xmm11
  movsd  %xmm4,0x30(%rsp)
  movapd %xmm5,%xmm10
  movsd  %xmm5,0x58(%rsp)
  movapd %xmm6,%xmm14
  movsd  %xmm6,0x60(%rsp)
  movapd %xmm7,%xmm9
  movsd  %xmm7,0x68(%rsp)
  mov    %rsi,%rbp
  mov    %fs:0x28,%rax
  mov    %rax,0x138(%rsp)
  xor    %eax,%eax
  movapd %xmm0,%xmm5
  subsd  D89(%rip),%xmm5
  movsd  %xmm5,0x40(%rsp)
  divsd  %xmm5,%xmm0
  movapd %xmm3,%xmm2
  mulsd  %xmm3,%xmm2
  movapd %xmm2,%xmm1
  movapd %xmm4,%xmm3
  mulsd  %xmm4,%xmm3
  addsd  %xmm3,%xmm1
  movsd  D97(%rip),%xmm4
  mulsd  %xmm4,%xmm1
  movapd %xmm15,%xmm2
  mulsd  %xmm0,%xmm2
  divsd  %xmm12,%xmm2
  addsd  %xmm2,%xmm1
  movapd %xmm1,%xmm8
  movsd  %xmm1,0x48(%rsp)
  movsd  0x168(%rsp),%xmm1
  mulsd  %xmm1,%xmm1
  movsd  0x170(%rsp),%xmm2
  mulsd  %xmm2,%xmm2
  addsd  %xmm2,%xmm1
  mulsd  %xmm4,%xmm1
  mulsd  %xmm7,%xmm0
  divsd  0x160(%rsp),%xmm0
  addsd  %xmm0,%xmm1
  movsd  %xmm1,0x10(%rsp)
  movapd %xmm12,%xmm0
  mulsd  %xmm13,%xmm0
  movsd  0x160(%rsp),%xmm6
  mulsd  0x168(%rsp),%xmm6
  movapd %xmm12,%xmm3
  mulsd  %xmm11,%xmm3
  movsd  0x160(%rsp),%xmm5
  mulsd  0x170(%rsp),%xmm5
  movapd %xmm0,%xmm2
  addsd  %xmm6,%xmm2
  movapd %xmm2,%xmm1
  mulsd  %xmm4,%xmm1
  mulsd  %xmm10,%xmm1
  movapd %xmm3,%xmm2
  addsd  %xmm5,%xmm2
  mulsd  %xmm4,%xmm2
  mulsd  %xmm14,%xmm2
  addsd  %xmm2,%xmm1
  movsd  %xmm1,(%rdi)
  movapd %xmm11,%xmm7
  mulsd  %xmm0,%xmm7
  movapd %xmm6,%xmm1
  mulsd  0x170(%rsp),%xmm1
  addsd  %xmm1,%xmm7
  mulsd  %xmm4,%xmm7
  mulsd  %xmm0,%xmm13
  movapd %xmm13,%xmm1
  addsd  %xmm15,%xmm1
  movapd %xmm6,%xmm2
  mulsd  0x168(%rsp),%xmm2
  addsd  %xmm2,%xmm1
  addsd  %xmm9,%xmm1
  mulsd  %xmm4,%xmm1
  mulsd  %xmm10,%xmm1
  movapd %xmm14,%xmm2
  mulsd  %xmm7,%xmm2
  addsd  %xmm2,%xmm1
  movsd  %xmm1,0x8(%rdi)
  mulsd  %xmm3,%xmm11
  movapd %xmm11,%xmm2
  addsd  %xmm15,%xmm2
  movapd %xmm5,%xmm1
  mulsd  0x170(%rsp),%xmm1
  addsd  %xmm2,%xmm1
  addsd  %xmm9,%xmm1
  mulsd  %xmm4,%xmm1
  mulsd  %xmm14,%xmm1
  mulsd  %xmm10,%xmm7
  addsd  %xmm7,%xmm1
  movsd  %xmm1,0x10(%rdi)
  mulsd  %xmm8,%xmm0
  movsd  0x10(%rsp),%xmm7
  mulsd  %xmm7,%xmm6
  addsd  %xmm6,%xmm0
  mulsd  %xmm4,%xmm0
  mulsd  %xmm10,%xmm0
  mulsd  %xmm8,%xmm3
  mulsd  %xmm7,%xmm5
  addsd  %xmm5,%xmm3
  mulsd  %xmm4,%xmm3
  mulsd  %xmm14,%xmm3
  addsd  %xmm3,%xmm0
  movsd  %xmm0,0x18(%rdi)
  movapd %xmm12,%xmm0
  mulsd  0x160(%rsp),%xmm0
  pxor   %xmm1,%xmm1
  ucomisd %xmm0,%xmm1
  ja     L10
  sqrtsd %xmm0,%xmm0
  movsd  %xmm0,0x38(%rsp)
L278:
  pxor   %xmm0,%xmm0
  movsd  0x8(%rsp),%xmm3
  ucomisd %xmm3,%xmm0
  ja     L86
  sqrtsd %xmm3,%xmm3
  movapd %xmm3,%xmm0
L25:
  mulsd  0x28(%rsp),%xmm0
  movsd  %xmm0,0x18(%rsp)
  pxor   %xmm0,%xmm0
  ucomisd 0x160(%rsp),%xmm0
  ja     L15
  movsd  0x160(%rsp),%xmm0
  sqrtsd %xmm0,%xmm0
L143:
  mulsd  0x168(%rsp),%xmm0
  addsd  0x18(%rsp),%xmm0
  movq   %xmm0,%r14
  pxor   %xmm0,%xmm0
  movsd  0x8(%rsp),%xmm6
  ucomisd %xmm6,%xmm0
  ja     L119
  sqrtsd %xmm6,%xmm6
  movsd  %xmm6,0x18(%rsp)
L120:
  pxor   %xmm0,%xmm0
  ucomisd 0x160(%rsp),%xmm0
  ja     L31
  movsd  0x160(%rsp),%xmm0
  sqrtsd %xmm0,%xmm0
L129:
  addsd  0x18(%rsp),%xmm0
  movq   %r14,%xmm7
  divsd  %xmm0,%xmm7
  movsd  %xmm7,0x18(%rsp)
  pxor   %xmm0,%xmm0
  movsd  0x8(%rsp),%xmm7
  ucomisd %xmm7,%xmm0
  ja     L201
  sqrtsd %xmm7,%xmm7
  movapd %xmm7,%xmm0
L191:
  mulsd  0x30(%rsp),%xmm0
  movsd  %xmm0,0x20(%rsp)
  pxor   %xmm0,%xmm0
  ucomisd 0x160(%rsp),%xmm0
  ja     L149
  movsd  0x160(%rsp),%xmm0
  sqrtsd %xmm0,%xmm0
L92:
  mulsd  0x170(%rsp),%xmm0
  addsd  0x20(%rsp),%xmm0
  movq   %xmm0,%r14
  pxor   %xmm0,%xmm0
  movsd  0x8(%rsp),%xmm3
  ucomisd %xmm3,%xmm0
  ja     L64
  sqrtsd %xmm3,%xmm3
  movsd  %xmm3,0x20(%rsp)
L293:
  pxor   %xmm0,%xmm0
  ucomisd 0x160(%rsp),%xmm0
  ja     L171
  movsd  0x160(%rsp),%xmm0
  sqrtsd %xmm0,%xmm0
L98:
  addsd  0x20(%rsp),%xmm0
  movq   %r14,%xmm7
  divsd  %xmm0,%xmm7
  movsd  %xmm7,0x20(%rsp)
  pxor   %xmm0,%xmm0
  movsd  0x8(%rsp),%xmm4
  ucomisd %xmm4,%xmm0
  ja     L260
  sqrtsd %xmm4,%xmm4
  movapd %xmm4,%xmm0
L124:
  mulsd  0x48(%rsp),%xmm0
  movsd  %xmm0,0x48(%rsp)
  pxor   %xmm0,%xmm0
  ucomisd 0x160(%rsp),%xmm0
  ja     L301
  movsd  0x160(%rsp),%xmm0
  sqrtsd %xmm0,%xmm0
L35:
  mulsd  0x10(%rsp),%xmm0
  addsd  0x48(%rsp),%xmm0
  movq   %xmm0,%r14
  pxor   %xmm0,%xmm0
  movsd  0x8(%rsp),%xmm5
  ucomisd %xmm5,%xmm0
  ja     L49
  sqrtsd %xmm5,%xmm5
  movsd  %xmm5,0x10(%rsp)
L166:
  pxor   %xmm0,%xmm0
  ucomisd 0x160(%rsp),%xmm0
  ja     L302
  movsd  0x160(%rsp),%xmm0
  sqrtsd %xmm0,%xmm0
L114:
  addsd  0x10(%rsp),%xmm0
  movq   %r14,%xmm5
  divsd  %xmm0,%xmm5
  movsd  %xmm5,0x10(%rsp)
  movsd  D97(%rip),%xmm0
  movsd  0x18(%rsp),%xmm6
  movapd %xmm6,%xmm2
  mulsd  %xmm0,%xmm2
  mulsd  %xmm6,%xmm2
  subsd  %xmm2,%xmm5
  movapd %xmm5,%xmm1
  movsd  0x20(%rsp),%xmm7
  mulsd  %xmm7,%xmm0
  mulsd  %xmm7,%xmm0
  subsd  %xmm0,%xmm1
  mulsd  0x40(%rsp),%xmm1
  pxor   %xmm0,%xmm0
  ucomisd %xmm1,%xmm0
  ja     L174
  sqrtsd %xmm1,%xmm1
L265:
  movsd  0x58(%rsp),%xmm9
  movapd %xmm9,%xmm3
  movsd  0x18(%rsp),%xmm5
  mulsd  %xmm5,%xmm3
  movsd  0x60(%rsp),%xmm8
  movapd %xmm8,%xmm0
  movsd  0x20(%rsp),%xmm10
  mulsd  %xmm10,%xmm0
  addsd  %xmm0,%xmm3
  movsd  D89(%rip),%xmm0
  movsd  %xmm0,0xb0(%rsp)
  movsd  %xmm0,0xb8(%rsp)
  movsd  %xmm0,0xc0(%rsp)
  movq   $0x0,0xc8(%rsp)
  movapd %xmm9,%xmm0
  mulsd  %xmm1,%xmm0
  movapd %xmm5,%xmm2
  subsd  %xmm0,%xmm2
  movsd  %xmm2,0xd0(%rsp)
  movsd  %xmm5,0xd8(%rsp)
  addsd  %xmm5,%xmm0
  movsd  %xmm0,0xe0(%rsp)
  movapd %xmm8,%xmm0
  xorpd  D58(%rip),%xmm0
  movsd  %xmm0,0xe8(%rsp)
  movapd %xmm8,%xmm0
  mulsd  %xmm1,%xmm0
  movapd %xmm10,%xmm2
  subsd  %xmm0,%xmm2
  movsd  %xmm2,0xf0(%rsp)
  movsd  %xmm10,0xf8(%rsp)
  movapd %xmm10,%xmm4
  addsd  %xmm10,%xmm0
  movsd  %xmm0,0x100(%rsp)
  movsd  %xmm9,0x108(%rsp)
  movapd %xmm3,%xmm2
  mulsd  %xmm1,%xmm2
  movsd  0x10(%rsp),%xmm12
  movapd %xmm12,%xmm0
  subsd  %xmm2,%xmm0
  movsd  %xmm0,0x110(%rsp)
  movapd %xmm5,%xmm0
  movapd %xmm5,%xmm6
  mulsd  %xmm5,%xmm0
  mulsd  %xmm10,%xmm4
  addsd  %xmm4,%xmm0
  movsd  D97(%rip),%xmm5
  mulsd  %xmm5,%xmm0
  movsd  %xmm0,0x118(%rsp)
  addsd  %xmm12,%xmm2
  movsd  %xmm2,0x120(%rsp)
  mulsd  %xmm8,%xmm6
  movapd %xmm10,%xmm0
  mulsd  %xmm9,%xmm0
  subsd  %xmm6,%xmm0
  movsd  %xmm0,0x128(%rsp)
  movsd  0x68(%rsp),%xmm2
  subsd  0x50(%rsp),%xmm2
  movapd %xmm9,%xmm4
  mulsd  0x168(%rsp),%xmm4
  movapd %xmm8,%xmm0
  mulsd  0x170(%rsp),%xmm0
  addsd  %xmm0,%xmm4
  movsd  0x28(%rsp),%xmm13
  movapd %xmm13,%xmm0
  mulsd  %xmm9,%xmm0
  movsd  0x30(%rsp),%xmm15
  movapd %xmm15,%xmm6
  mulsd  %xmm8,%xmm6
  addsd  %xmm6,%xmm0
  subsd  %xmm0,%xmm4
  movsd  0x38(%rsp),%xmm11
  movapd %xmm11,%xmm0
  mulsd  %xmm1,%xmm0
  mulsd  %xmm0,%xmm4
  movapd %xmm1,%xmm6
  mulsd  %xmm1,%xmm6
  movapd %xmm2,%xmm0
  subsd  %xmm4,%xmm0
  mulsd  %xmm5,%xmm0
  divsd  %xmm6,%xmm0
  movsd  %xmm0,0x90(%rsp)
  movsd  0x160(%rsp),%xmm0
  subsd  0x8(%rsp),%xmm0
  movapd %xmm2,%xmm7
  divsd  %xmm6,%xmm7
  subsd  %xmm7,%xmm0
  movsd  %xmm0,0x98(%rsp)
  addsd  %xmm4,%xmm2
  mulsd  %xmm5,%xmm2
  divsd  %xmm6,%xmm2
  movsd  %xmm2,0xa0(%rsp)
  movapd %xmm8,%xmm2
  mulsd  0x168(%rsp),%xmm2
  movapd %xmm9,%xmm0
  mulsd  0x170(%rsp),%xmm0
  subsd  %xmm2,%xmm0
  movapd %xmm13,%xmm4
  mulsd  %xmm8,%xmm4
  movapd %xmm15,%xmm2
  mulsd  %xmm9,%xmm2
  subsd  %xmm4,%xmm2
  subsd  %xmm2,%xmm0
  mulsd  %xmm11,%xmm0
  movsd  %xmm0,0xa8(%rsp)
  movapd %xmm3,%xmm2
  subsd  %xmm1,%xmm2
  movq   D22(%rip),%xmm0
  movapd %xmm2,%xmm4
  andpd  %xmm0,%xmm4
  movsd  %xmm4,0x70(%rsp)
  movapd %xmm3,%xmm5
  andpd  %xmm0,%xmm5
  movsd  %xmm5,0x78(%rsp)
  addsd  %xmm1,%xmm3
  andpd  %xmm3,%xmm0
  movsd  %xmm0,0x80(%rsp)
  movsd  %xmm5,0x88(%rsp)
  movsd  0x178(%rsp),%xmm5
  comisd %xmm4,%xmm5
  jbe    L249
  mulsd  %xmm2,%xmm2
  movapd %xmm5,%xmm1
  mulsd  %xmm5,%xmm1
  addsd  %xmm1,%xmm2
  movsd  D97(%rip),%xmm1
  divsd  %xmm5,%xmm1
  mulsd  %xmm1,%xmm2
  movsd  %xmm2,0x70(%rsp)
L249:
  movsd  0x178(%rsp),%xmm6
  comisd %xmm0,%xmm6
  jbe    L181
  mulsd  %xmm3,%xmm3
  movapd %xmm6,%xmm0
  mulsd  %xmm6,%xmm0
  addsd  %xmm0,%xmm3
  movsd  D97(%rip),%xmm0
  divsd  %xmm6,%xmm0
  mulsd  %xmm0,%xmm3
  movsd  %xmm3,0x80(%rsp)
L181:
  movq   $0x0,0x0(%rbp)
  lea    0xb0(%rsp),%rcx
  mov    $0x0,%edx
  movsd  D84(%rip),%xmm2
  lea    0x90(%rsp),%rsi
L314:
  movsd  0x0(%rbp),%xmm0
  maxsd  0x70(%rsp,%rdx,8),%xmm0
  movsd  %xmm0,0x0(%rbp)
  movsd  (%rbx,%rdx,8),%xmm1
  mov    $0x0,%eax
L93:
  movapd %xmm2,%xmm0
  mulsd  0x70(%rsp,%rax,1),%xmm0
  mulsd  (%rax,%rsi,1),%xmm0
  mulsd  (%rcx,%rax,1),%xmm0
  addsd  %xmm0,%xmm1
  add    $0x8,%rax
  cmp    $0x20,%rax
  jne    L93
  movsd  %xmm1,(%rbx,%rdx,8)
  add    $0x1,%rdx
  add    $0x20,%rcx
  cmp    $0x4,%rdx
  jne    L314
  mov    0x138(%rsp),%rax
  sub    %fs:0x28,%rax
  jne    L18
  add    $0x140,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r14
  ret
L10:
  call   <sqrt@plt>
  movsd  %xmm0,0x38(%rsp)
  jmp    L278
L86:
  movsd  0x8(%rsp),%xmm0
  call   <sqrt@plt>
  jmp    L25
L15:
  movsd  0x160(%rsp),%xmm0
  call   <sqrt@plt>
  jmp    L143
L119:
  movsd  0x8(%rsp),%xmm0
  call   <sqrt@plt>
  movsd  %xmm0,0x18(%rsp)
  jmp    L120
L31:
  movsd  0x160(%rsp),%xmm0
  call   <sqrt@plt>
  jmp    L129
L201:
  movsd  0x8(%rsp),%xmm0
  call   <sqrt@plt>
  jmp    L191
L149:
  movsd  0x160(%rsp),%xmm0
  call   <sqrt@plt>
  jmp    L92
L64:
  movsd  0x8(%rsp),%xmm0
  call   <sqrt@plt>
  movsd  %xmm0,0x20(%rsp)
  jmp    L293
L171:
  movsd  0x160(%rsp),%xmm0
  call   <sqrt@plt>
  jmp    L98
L260:
  movsd  0x8(%rsp),%xmm0
  call   <sqrt@plt>
  jmp    L124
L301:
  movsd  0x160(%rsp),%xmm0
  call   <sqrt@plt>
  jmp    L35
L49:
  movsd  0x8(%rsp),%xmm0
  call   <sqrt@plt>
  movsd  %xmm0,0x10(%rsp)
  jmp    L166
L302:
  movsd  0x160(%rsp),%xmm0
  call   <sqrt@plt>
  jmp    L114
L174:
  movapd %xmm1,%xmm0
  call   <sqrt@plt>
  movapd %xmm0,%xmm1
  jmp    L265
L18:
  call   <__stack_chk_fail@plt>
```

tool_calls: name: unpack_data
arguments: data_tag: D89
type_guessed: f64
name: unpack_data
arguments: data_tag: D97
type_guessed: f64
name: unpack_data
arguments: data_tag: D58
type_guessed: f32[4]
name: unpack_data
arguments: data_tag: D22
type_guessed: qword
name: unpack_data
arguments: data_tag: D84
type_guessed: f32[2]

tool: D89: 1

tool: D97: 0.5

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

tool: D22: 9223372036854775807

tool: {"D84": [0.0,-1.75]}

assistant: ```c
void ROE_2D_solver(double *F, double gamma, double P_L, double RHO_L, double U_L, double V_L, double n_x, double n_y, double P_R, double RHO_R, double U_R, double V_R, double *lambda_max, double delta) {
  double H_L, H_R;
  H_L = gamma / (gamma - 1.0) * P_L / RHO_L + 0.5 * (U_L * U_L + V_L * V_L);
  H_R = gamma / (gamma - 1.0) * P_R / RHO_R + 0.5 * (U_R * U_R + V_R * V_R);
  F[0] = 0.5 * (RHO_L * U_L + RHO_R * U_R) * n_x + 0.5 * (RHO_L * V_L + RHO_R * V_R) * n_y;
  F[1] = 0.5 * (RHO_L * U_L * U_L + P_L + RHO_R * U_R * U_R + P_R) * n_x + 0.5 * (RHO_L * U_L * V_L + RHO_R * U_R * V_R) * n_y;
  F[2] = 0.5 * (RHO_L * U_L * V_L + RHO_R * U_R * V_R) * n_x + 0.5 * (RHO_L * V_L * V_L + P_L + RHO_R * V_R * V_R + P_R) * n_y;
  F[3] = 0.5 * (RHO_L * U_L * H_L + RHO_R * U_R * H_R) * n_x + 0.5 * (RHO_L * V_L * H_L + RHO_R * V_R * H_R) * n_y;
  double RHO_S, U_S, V_S, H_S, C_S;
  RHO_S = sqrt(RHO_L * RHO_R);
  U_S = (U_L * sqrt(RHO_L) + U_R * sqrt(RHO_R)) / (sqrt(RHO_L) + sqrt(RHO_R));
  V_S = (V_L * sqrt(RHO_L) + V_R * sqrt(RHO_R)) / (sqrt(RHO_L) + sqrt(RHO_R));
  H_S = (H_L * sqrt(RHO_L) + H_R * sqrt(RHO_R)) / (sqrt(RHO_L) + sqrt(RHO_R));
  C_S = sqrt((gamma - 1.0) * (H_S - 0.5 * U_S * U_S - 0.5 * V_S * V_S));
  double qn_S, qt_S;
  qn_S = U_S * n_x + V_S * n_y;
  qt_S = -U_S * n_y + V_S * n_x;
  double qn_R, qt_R;
  qn_R = U_R * n_x + V_R * n_y;
  qt_R = -U_R * n_y + V_R * n_x;
  double qn_L, qt_L;
  qn_L = U_L * n_x + V_L * n_y;
  qt_L = -U_L * n_y + V_L * n_x;
  double R[4][4];
  double lambda[4], W[4];
  R[0][0] = 1.0;
  R[0][1] = 1.0;
  R[0][2] = 1.0;
  R[0][3] = 0.0;
  R[1][0] = U_S - C_S * n_x;
  R[1][1] = U_S;
  R[1][2] = U_S + C_S * n_x;
  R[1][3] = -n_y;
  R[2][0] = V_S - C_S * n_y;
  R[2][1] = V_S;
  R[2][2] = V_S + C_S * n_y;
  R[2][3] = n_x;
  R[3][0] = H_S - qn_S * C_S;
  R[3][1] = 0.5 * (U_S * U_S + V_S * V_S);
  R[3][2] = H_S + qn_S * C_S;
  R[3][3] = qt_S;
  int i, j;
  W[0] = 0.5 * ((P_R - P_L) - RHO_S * C_S * (qn_R - qn_L)) / (C_S * C_S);
  W[1] = (RHO_R - RHO_L) - (P_R - P_L) / (C_S * C_S);
  W[2] = 0.5 * ((P_R - P_L) + RHO_S * C_S * (qn_R - qn_L)) / (C_S * C_S);
  W[3] = RHO_S * (qt_R - qt_L);
  lambda[0] = fabs(qn_S - C_S);
  lambda[1] = fabs(qn_S);
  lambda[2] = fabs(qn_S + C_S);
  lambda[3] = fabs(qn_S);
  double delta_1 = 0.01;
  double delta_2 = 0.2;
  if (lambda[0] < delta)
    lambda[0] = 0.5 / delta * (lambda[0] * lambda[0] + delta * delta);
  if (lambda[2] < delta)
    lambda[2] = 0.5 / delta * (lambda[2] * lambda[2] + delta * delta);
  *lambda_max = 0;
  for (i = 0; i < 4; i++) {
    *lambda_max = (*lambda_max > lambda[i] ? *lambda_max : lambda[i]);
    for (j = 0; j < 4; j++) {
      F[i] += -0.5 * lambda[j] * W[j] * R[i][j];
    }
  }
}
```