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

```asm
<ROE_2D_solver>:
  endbr64
  push   %rbp
  movapd %xmm1,%xmm13
  movapd %xmm0,%xmm1
  movapd %xmm4,%xmm8
  push   %rbx
  mulsd  %xmm4,%xmm8
  movapd %xmm7,%xmm15
  mov    %rdi,%rbx
  mov    %rsi,%rbp
  sub    $0x148,%rsp
  subsd  D145(%rip),%xmm1
  movsd  D37(%rip),%xmm10
  mov    %fs:0x28,%rax
  mov    %rax,0x138(%rsp)
  xor    %eax,%eax
  movsd  %xmm15,(%rsp)
  movsd  0x160(%rsp),%xmm7
  movsd  0x168(%rsp),%xmm11
  movsd  0x170(%rsp),%xmm12
  movsd  %xmm13,0x8(%rsp)
  divsd  %xmm1,%xmm0
  movsd  %xmm1,0x18(%rsp)
  movapd %xmm3,%xmm1
  mulsd  %xmm3,%xmm1
  mulsd  %xmm7,%xmm11
  mulsd  %xmm7,%xmm12
  addsd  %xmm8,%xmm1
  movapd %xmm13,%xmm8
  mulsd  %xmm10,%xmm1
  mulsd  %xmm0,%xmm8
  mulsd  %xmm15,%xmm0
  divsd  %xmm2,%xmm8
  divsd  %xmm7,%xmm0
  addsd  %xmm8,%xmm1
  movsd  0x170(%rsp),%xmm8
  mulsd  %xmm8,%xmm8
  movapd %xmm1,%xmm14
  movsd  0x168(%rsp),%xmm1
  mulsd  %xmm1,%xmm1
  addsd  %xmm8,%xmm1
  movapd %xmm2,%xmm8
  mulsd  %xmm4,%xmm8
  mulsd  %xmm10,%xmm1
  movapd %xmm8,%xmm9
  addsd  %xmm12,%xmm9
  addsd  %xmm0,%xmm1
  mulsd  %xmm10,%xmm9
  movapd %xmm2,%xmm0
  mulsd  %xmm3,%xmm0
  movapd %xmm1,%xmm15
  mulsd  %xmm6,%xmm9
  movapd %xmm0,%xmm1
  addsd  %xmm11,%xmm1
  mulsd  %xmm10,%xmm1
  mulsd  %xmm5,%xmm1
  addsd  %xmm9,%xmm1
  movapd %xmm0,%xmm9
  mulsd  %xmm4,%xmm9
  movsd  %xmm1,(%rdi)
  movsd  0x170(%rsp),%xmm1
  mulsd  %xmm11,%xmm1
  addsd  %xmm1,%xmm9
  movapd %xmm0,%xmm1
  mulsd  %xmm3,%xmm1
  mulsd  %xmm14,%xmm0
  mulsd  %xmm10,%xmm9
  addsd  %xmm13,%xmm1
  movsd  0x168(%rsp),%xmm13
  mulsd  %xmm11,%xmm13
  mulsd  %xmm15,%xmm11
  addsd  %xmm13,%xmm1
  addsd  (%rsp),%xmm1
  movapd %xmm9,%xmm13
  mulsd  %xmm6,%xmm13
  addsd  %xmm11,%xmm0
  mulsd  %xmm5,%xmm9
  mulsd  %xmm10,%xmm1
  mulsd  %xmm10,%xmm0
  mulsd  %xmm5,%xmm1
  mulsd  %xmm5,%xmm0
  addsd  %xmm13,%xmm1
  movsd  0x170(%rsp),%xmm13
  mulsd  %xmm12,%xmm13
  movsd  %xmm1,0x8(%rdi)
  mulsd  %xmm15,%xmm12
  movapd %xmm8,%xmm1
  mulsd  %xmm14,%xmm8
  mulsd  %xmm4,%xmm1
  addsd  0x8(%rsp),%xmm1
  addsd  %xmm12,%xmm8
  addsd  %xmm13,%xmm1
  addsd  (%rsp),%xmm1
  mulsd  %xmm10,%xmm8
  mulsd  %xmm10,%xmm1
  mulsd  %xmm6,%xmm8
  mulsd  %xmm6,%xmm1
  addsd  %xmm8,%xmm0
  addsd  %xmm9,%xmm1
  movsd  %xmm0,0x18(%rdi)
  movapd %xmm2,%xmm0
  mulsd  %xmm7,%xmm0
  movsd  %xmm1,0x10(%rdi)
  pxor   %xmm1,%xmm1
  ucomisd %xmm0,%xmm1
  ja     L326
  sqrtsd %xmm0,%xmm0
  movsd  %xmm0,0x10(%rsp)
L109:
  ucomisd %xmm2,%xmm1
  ja     L282
  movapd %xmm2,%xmm8
  sqrtsd %xmm8,%xmm8
L43:
  mulsd  %xmm3,%xmm8
  ucomisd %xmm7,%xmm1
  ja     L117
  movapd %xmm7,%xmm9
  sqrtsd %xmm9,%xmm9
L222:
  mulsd  0x168(%rsp),%xmm9
  ucomisd %xmm2,%xmm1
  addsd  %xmm8,%xmm9
  ja     L286
  movapd %xmm2,%xmm8
  sqrtsd %xmm8,%xmm8
L164:
  ucomisd %xmm7,%xmm1
  ja     L348
  movapd %xmm7,%xmm0
  sqrtsd %xmm0,%xmm0
L319:
  addsd  %xmm8,%xmm0
  ucomisd %xmm2,%xmm1
  divsd  %xmm0,%xmm9
  ja     L281
  movapd %xmm2,%xmm11
  sqrtsd %xmm11,%xmm11
L90:
  mulsd  %xmm4,%xmm11
  ucomisd %xmm7,%xmm1
  ja     L55
  movapd %xmm7,%xmm8
  sqrtsd %xmm8,%xmm8
L85:
  mulsd  0x170(%rsp),%xmm8
  ucomisd %xmm2,%xmm1
  addsd  %xmm11,%xmm8
  ja     L349
  movapd %xmm2,%xmm11
  sqrtsd %xmm11,%xmm11
L14:
  ucomisd %xmm7,%xmm1
  ja     L106
  movapd %xmm7,%xmm0
  sqrtsd %xmm0,%xmm0
L182:
  addsd  %xmm11,%xmm0
  ucomisd %xmm2,%xmm1
  divsd  %xmm0,%xmm8
  ja     L113
  movapd %xmm2,%xmm11
  sqrtsd %xmm11,%xmm11
L288:
  mulsd  %xmm14,%xmm11
  ucomisd %xmm7,%xmm1
  ja     L252
  movapd %xmm7,%xmm0
  sqrtsd %xmm0,%xmm0
L33:
  mulsd  %xmm0,%xmm15
  ucomisd %xmm2,%xmm1
  addsd  %xmm11,%xmm15
  ja     L244
  movapd %xmm2,%xmm11
  sqrtsd %xmm11,%xmm11
L169:
  ucomisd %xmm7,%xmm1
  ja     L100
  movapd %xmm7,%xmm0
  sqrtsd %xmm0,%xmm0
L206:
  addsd  %xmm11,%xmm0
  movapd %xmm9,%xmm11
  mulsd  %xmm10,%xmm11
  divsd  %xmm0,%xmm15
  mulsd  %xmm9,%xmm11
  movapd %xmm15,%xmm0
  movapd %xmm15,%xmm12
  subsd  %xmm11,%xmm0
  movapd %xmm8,%xmm11
  mulsd  %xmm10,%xmm11
  mulsd  %xmm8,%xmm11
  subsd  %xmm11,%xmm0
  mulsd  0x18(%rsp),%xmm0
  ucomisd %xmm0,%xmm1
  ja     L287
  sqrtsd %xmm0,%xmm0
L96:
  movapd %xmm5,%xmm11
  movapd %xmm6,%xmm13
  movapd %xmm9,%xmm14
  movsd  %xmm9,0xd8(%rsp)
  mulsd  %xmm8,%xmm13
  movapd %xmm8,%xmm15
  movsd  %xmm8,0xf8(%rsp)
  subsd  %xmm2,%xmm7
  mulsd  %xmm9,%xmm11
  mov    D145(%rip),%rax
  movsd  %xmm5,0x108(%rsp)
  mulsd  %xmm8,%xmm15
  movq   $0x0,0xc8(%rsp)
  mov    %rax,0xb8(%rsp)
  mov    %rax,0xc0(%rsp)
  addsd  %xmm13,%xmm11
  movapd %xmm5,%xmm13
  mulsd  %xmm0,%xmm13
  subsd  %xmm13,%xmm14
  addsd  %xmm9,%xmm13
  movsd  %xmm13,0xe0(%rsp)
  movapd %xmm6,%xmm13
  xorpd  D78(%rip),%xmm13
  movsd  %xmm14,0xd0(%rsp)
  movapd %xmm8,%xmm14
  movsd  %xmm13,0xe8(%rsp)
  movapd %xmm6,%xmm13
  mulsd  %xmm0,%xmm13
  subsd  %xmm13,%xmm14
  addsd  %xmm8,%xmm13
  mulsd  %xmm5,%xmm8
  movsd  %xmm14,0xf0(%rsp)
  movapd %xmm11,%xmm14
  movsd  %xmm13,0x100(%rsp)
  mulsd  %xmm0,%xmm14
  movapd %xmm12,%xmm13
  subsd  %xmm14,%xmm13
  addsd  %xmm12,%xmm14
  movsd  0x170(%rsp),%xmm12
  mulsd  %xmm6,%xmm12
  movsd  %xmm13,0x110(%rsp)
  movapd %xmm9,%xmm13
  movsd  %xmm14,0x120(%rsp)
  mulsd  %xmm9,%xmm13
  movsd  0x10(%rsp),%xmm14
  mulsd  %xmm6,%xmm9
  addsd  %xmm15,%xmm13
  subsd  %xmm9,%xmm8
  movsd  0x168(%rsp),%xmm9
  mulsd  %xmm10,%xmm13
  movsd  %xmm8,0x128(%rsp)
  mulsd  %xmm5,%xmm9
  movsd  (%rsp),%xmm8
  subsd  0x8(%rsp),%xmm8
  movapd %xmm8,%xmm2
  movsd  %xmm13,0x118(%rsp)
  movapd %xmm3,%xmm13
  addsd  %xmm12,%xmm9
  mulsd  %xmm5,%xmm13
  movapd %xmm4,%xmm12
  mulsd  %xmm6,%xmm12
  mulsd  %xmm6,%xmm3
  mulsd  %xmm5,%xmm4
  addsd  %xmm13,%xmm12
  movapd %xmm8,%xmm13
  subsd  %xmm3,%xmm4
  movapd %xmm11,%xmm3
  subsd  %xmm12,%xmm9
  movapd %xmm14,%xmm12
  subsd  %xmm0,%xmm3
  mulsd  %xmm0,%xmm12
  mulsd  %xmm12,%xmm9
  movapd %xmm0,%xmm12
  mulsd  %xmm0,%xmm12
  addsd  %xmm9,%xmm8
  subsd  %xmm9,%xmm13
  divsd  %xmm12,%xmm2
  mulsd  %xmm10,%xmm8
  mulsd  %xmm10,%xmm13
  divsd  %xmm12,%xmm8
  subsd  %xmm2,%xmm7
  movsd  0x170(%rsp),%xmm2
  mulsd  %xmm5,%xmm2
  movapd %xmm3,%xmm5
  movsd  %xmm7,0x98(%rsp)
  movsd  0x168(%rsp),%xmm7
  mulsd  %xmm6,%xmm7
  movsd  0x178(%rsp),%xmm6
  subsd  %xmm7,%xmm2
  divsd  %xmm12,%xmm13
  subsd  %xmm4,%xmm2
  movapd %xmm11,%xmm4
  addsd  %xmm0,%xmm11
  mulsd  %xmm14,%xmm2
  movsd  %xmm8,0xa0(%rsp)
  movsd  %xmm2,0xa8(%rsp)
  movq   D169(%rip),%xmm2
  andpd  %xmm2,%xmm5
  andpd  %xmm2,%xmm4
  andpd  %xmm11,%xmm2
  comisd %xmm5,%xmm6
  movsd  %xmm4,0x78(%rsp)
  movsd  %xmm2,0x80(%rsp)
  movsd  %xmm4,0x88(%rsp)
  jbe    L303
  mulsd  %xmm3,%xmm3
  movapd %xmm10,%xmm5
  movapd %xmm6,%xmm0
  mulsd  %xmm6,%xmm0
  divsd  %xmm6,%xmm5
  addsd  %xmm0,%xmm3
  mulsd  %xmm3,%xmm5
L303:
  movsd  0x178(%rsp),%xmm6
  comisd %xmm2,%xmm6
  jbe    L324
  divsd  %xmm6,%xmm10
  movapd %xmm6,%xmm0
  mulsd  %xmm11,%xmm11
  mulsd  %xmm6,%xmm0
  addsd  %xmm0,%xmm11
  mulsd  %xmm10,%xmm11
  movsd  %xmm11,0x80(%rsp)
L324:
  mov    D145(%rip),%rax
  movsd  D142(%rip),%xmm4
  lea    0xb0(%rsp),%rcx
  xor    %edx,%edx
  movapd %xmm5,%xmm2
  movq   %rax,%xmm0
L255:
  maxsd  %xmm2,%xmm1
  movapd %xmm13,%xmm3
  xor    %eax,%eax
  movsd  %xmm1,0x0(%rbp)
  movsd  (%rbx,%rdx,1),%xmm2
  movapd %xmm5,%xmm1
L259:
  mulsd  %xmm4,%xmm1
  add    $0x8,%rax
  mulsd  %xmm3,%xmm1
  mulsd  %xmm0,%xmm1
  addsd  %xmm1,%xmm2
  cmp    $0x20,%rax
  je     L137
  movsd  0x70(%rsp,%rax,1),%xmm1
  movsd  (%rcx,%rax,1),%xmm0
  movsd  0x90(%rsp,%rax,1),%xmm3
  jmp    L259
  nop
L137:
  movsd  %xmm2,(%rbx,%rdx,1)
  add    $0x8,%rdx
  add    $0x20,%rcx
  cmp    $0x20,%rdx
  je     L132
  movsd  0x0(%rbp),%xmm1
  movsd  0x70(%rsp,%rdx,1),%xmm2
  movsd  (%rcx),%xmm0
  jmp    L255
  nopl   0x0(%rax)
L132:
  mov    0x138(%rsp),%rax
  sub    %fs:0x28,%rax
  jne    L124
  add    $0x148,%rsp
  pop    %rbx
  pop    %rbp
  ret
L326:
  movsd  %xmm7,0x160(%rsp)
  movsd  %xmm6,0x50(%rsp)
  movsd  %xmm5,0x48(%rsp)
  movsd  %xmm4,0x40(%rsp)
  movsd  %xmm3,0x38(%rsp)
  movsd  %xmm2,0x30(%rsp)
  movsd  %xmm15,0x28(%rsp)
  movsd  %xmm14,0x20(%rsp)
  call   <sqrt@plt>
  mov    D37(%rip),%rax
  pxor   %xmm1,%xmm1
  movsd  0x160(%rsp),%xmm7
  movsd  0x50(%rsp),%xmm6
  movsd  0x48(%rsp),%xmm5
  movsd  %xmm0,0x10(%rsp)
  movsd  0x40(%rsp),%xmm4
  movsd  0x38(%rsp),%xmm3
  movq   %rax,%xmm10
  movsd  0x30(%rsp),%xmm2
  movsd  0x28(%rsp),%xmm15
  movsd  0x20(%rsp),%xmm14
  jmp    L109
L124:
  call   <__stack_chk_fail@plt>
L287:
  movsd  %xmm7,0x160(%rsp)
  movsd  %xmm6,0x50(%rsp)
  movsd  %xmm5,0x48(%rsp)
  movsd  %xmm4,0x40(%rsp)
  movsd  %xmm3,0x38(%rsp)
  movsd  %xmm2,0x30(%rsp)
  movsd  %xmm8,0x20(%rsp)
  movsd  %xmm9,0x18(%rsp)
  movsd  %xmm15,0x28(%rsp)
  call   <sqrt@plt>
  mov    D37(%rip),%rax
  pxor   %xmm1,%xmm1
  movsd  0x160(%rsp),%xmm7
  movsd  0x50(%rsp),%xmm6
  movsd  0x48(%rsp),%xmm5
  movsd  0x40(%rsp),%xmm4
  movsd  0x38(%rsp),%xmm3
  movq   %rax,%xmm10
  movsd  0x30(%rsp),%xmm2
  movsd  0x28(%rsp),%xmm12
  movsd  0x20(%rsp),%xmm8
  movsd  0x18(%rsp),%xmm9
  jmp    L96
L100:
  movapd %xmm7,%xmm0
  movsd  %xmm6,0x60(%rsp)
  movsd  %xmm5,0x58(%rsp)
  movsd  %xmm4,0x50(%rsp)
  movsd  %xmm3,0x48(%rsp)
  movsd  %xmm2,0x40(%rsp)
  movsd  %xmm8,0x38(%rsp)
  movsd  %xmm9,0x30(%rsp)
  movsd  %xmm11,0x28(%rsp)
  movsd  %xmm15,0x20(%rsp)
  movsd  %xmm7,0x160(%rsp)
  call   <sqrt@plt>
  mov    D37(%rip),%rax
  pxor   %xmm1,%xmm1
  movsd  0x60(%rsp),%xmm6
  movsd  0x58(%rsp),%xmm5
  movsd  0x50(%rsp),%xmm4
  movsd  0x48(%rsp),%xmm3
  movsd  0x40(%rsp),%xmm2
  movq   %rax,%xmm10
  movsd  0x38(%rsp),%xmm8
  movsd  0x30(%rsp),%xmm9
  movsd  0x28(%rsp),%xmm11
  movsd  0x20(%rsp),%xmm15
  movsd  0x160(%rsp),%xmm7
  jmp    L206
L244:
  movapd %xmm2,%xmm0
  movsd  %xmm6,0x58(%rsp)
  movsd  %xmm7,0x160(%rsp)
  movsd  %xmm5,0x50(%rsp)
  movsd  %xmm4,0x48(%rsp)
  movsd  %xmm3,0x40(%rsp)
  movsd  %xmm8,0x38(%rsp)
  movsd  %xmm9,0x30(%rsp)
  movsd  %xmm15,0x28(%rsp)
  movsd  %xmm2,0x20(%rsp)
  call   <sqrt@plt>
  mov    D37(%rip),%rax
  pxor   %xmm1,%xmm1
  movsd  0x160(%rsp),%xmm7
  movsd  0x58(%rsp),%xmm6
  movsd  0x50(%rsp),%xmm5
  movapd %xmm0,%xmm11
  movsd  0x48(%rsp),%xmm4
  movsd  0x40(%rsp),%xmm3
  movq   %rax,%xmm10
  movsd  0x38(%rsp),%xmm8
  movsd  0x30(%rsp),%xmm9
  movsd  0x28(%rsp),%xmm15
  movsd  0x20(%rsp),%xmm2
  jmp    L169
L252:
  movapd %xmm7,%xmm0
  movsd  %xmm6,0x60(%rsp)
  movsd  %xmm5,0x58(%rsp)
  movsd  %xmm4,0x50(%rsp)
  movsd  %xmm3,0x48(%rsp)
  movsd  %xmm2,0x40(%rsp)
  movsd  %xmm8,0x38(%rsp)
  movsd  %xmm9,0x30(%rsp)
  movsd  %xmm15,0x28(%rsp)
  movsd  %xmm11,0x20(%rsp)
  movsd  %xmm7,0x160(%rsp)
  call   <sqrt@plt>
  mov    D37(%rip),%rax
  pxor   %xmm1,%xmm1
  movsd  0x60(%rsp),%xmm6
  movsd  0x58(%rsp),%xmm5
  movsd  0x50(%rsp),%xmm4
  movsd  0x48(%rsp),%xmm3
  movsd  0x40(%rsp),%xmm2
  movq   %rax,%xmm10
  movsd  0x38(%rsp),%xmm8
  movsd  0x30(%rsp),%xmm9
  movsd  0x28(%rsp),%xmm15
  movsd  0x20(%rsp),%xmm11
  movsd  0x160(%rsp),%xmm7
  jmp    L33
L113:
  movapd %xmm2,%xmm0
  movsd  %xmm6,0x60(%rsp)
  movsd  %xmm7,0x160(%rsp)
  movsd  %xmm5,0x58(%rsp)
  movsd  %xmm4,0x50(%rsp)
  movsd  %xmm3,0x48(%rsp)
  movsd  %xmm8,0x40(%rsp)
  movsd  %xmm9,0x38(%rsp)
  movsd  %xmm15,0x30(%rsp)
  movsd  %xmm14,0x28(%rsp)
  movsd  %xmm2,0x20(%rsp)
  call   <sqrt@plt>
  mov    D37(%rip),%rax
  pxor   %xmm1,%xmm1
  movsd  0x160(%rsp),%xmm7
  movsd  0x60(%rsp),%xmm6
  movsd  0x58(%rsp),%xmm5
  movapd %xmm0,%xmm11
  movsd  0x50(%rsp),%xmm4
  movsd  0x48(%rsp),%xmm3
  movq   %rax,%xmm10
  movsd  0x40(%rsp),%xmm8
  movsd  0x38(%rsp),%xmm9
  movsd  0x30(%rsp),%xmm15
  movsd  0x28(%rsp),%xmm14
  movsd  0x20(%rsp),%xmm2
  jmp    L288
L106:
  movapd %xmm7,%xmm0
  movsd  %xmm6,0x68(%rsp)
  movsd  %xmm5,0x60(%rsp)
  movsd  %xmm4,0x58(%rsp)
  movsd  %xmm3,0x50(%rsp)
  movsd  %xmm2,0x48(%rsp)
  movsd  %xmm9,0x40(%rsp)
  movsd  %xmm15,0x38(%rsp)
  movsd  %xmm14,0x30(%rsp)
  movsd  %xmm11,0x28(%rsp)
  movsd  %xmm8,0x20(%rsp)
  movsd  %xmm7,0x160(%rsp)
  call   <sqrt@plt>
  mov    D37(%rip),%rax
  pxor   %xmm1,%xmm1
  movsd  0x68(%rsp),%xmm6
  movsd  0x60(%rsp),%xmm5
  movsd  0x58(%rsp),%xmm4
  movsd  0x50(%rsp),%xmm3
  movsd  0x48(%rsp),%xmm2
  movq   %rax,%xmm10
  movsd  0x40(%rsp),%xmm9
  movsd  0x38(%rsp),%xmm15
  movsd  0x30(%rsp),%xmm14
  movsd  0x28(%rsp),%xmm11
  movsd  0x160(%rsp),%xmm7
  movsd  0x20(%rsp),%xmm8
  jmp    L182
L349:
  movapd %xmm2,%xmm0
  movsd  %xmm6,0x60(%rsp)
  movsd  %xmm7,0x160(%rsp)
  movsd  %xmm5,0x58(%rsp)
  movsd  %xmm4,0x50(%rsp)
  movsd  %xmm3,0x48(%rsp)
  movsd  %xmm9,0x40(%rsp)
  movsd  %xmm15,0x38(%rsp)
  movsd  %xmm14,0x30(%rsp)
  movsd  %xmm8,0x28(%rsp)
  movsd  %xmm2,0x20(%rsp)
  call   <sqrt@plt>
  mov    D37(%rip),%rax
  pxor   %xmm1,%xmm1
  movsd  0x160(%rsp),%xmm7
  movsd  0x60(%rsp),%xmm6
  movsd  0x58(%rsp),%xmm5
  movapd %xmm0,%xmm11
  movsd  0x50(%rsp),%xmm4
  movsd  0x48(%rsp),%xmm3
  movq   %rax,%xmm10
  movsd  0x40(%rsp),%xmm9
  movsd  0x38(%rsp),%xmm15
  movsd  0x30(%rsp),%xmm14
  movsd  0x28(%rsp),%xmm8
  movsd  0x20(%rsp),%xmm2
  jmp    L14
L55:
  movapd %xmm7,%xmm0
  movsd  %xmm6,0x60(%rsp)
  movsd  %xmm5,0x58(%rsp)
  movsd  %xmm4,0x50(%rsp)
  movsd  %xmm3,0x48(%rsp)
  movsd  %xmm2,0x40(%rsp)
  movsd  %xmm9,0x38(%rsp)
  movsd  %xmm15,0x30(%rsp)
  movsd  %xmm14,0x28(%rsp)
  movsd  %xmm11,0x20(%rsp)
  movsd  %xmm7,0x160(%rsp)
  call   <sqrt@plt>
  mov    D37(%rip),%rax
  pxor   %xmm1,%xmm1
  movsd  0x60(%rsp),%xmm6
  movsd  0x58(%rsp),%xmm5
  movsd  0x50(%rsp),%xmm4
  movapd %xmm0,%xmm8
  movsd  0x48(%rsp),%xmm3
  movsd  0x40(%rsp),%xmm2
  movq   %rax,%xmm10
  movsd  0x38(%rsp),%xmm9
  movsd  0x30(%rsp),%xmm15
  movsd  0x28(%rsp),%xmm14
  movsd  0x20(%rsp),%xmm11
  movsd  0x160(%rsp),%xmm7
  jmp    L85
L281:
  movapd %xmm2,%xmm0
  movsd  %xmm6,0x58(%rsp)
  movsd  %xmm7,0x160(%rsp)
  movsd  %xmm5,0x50(%rsp)
  movsd  %xmm4,0x48(%rsp)
  movsd  %xmm3,0x40(%rsp)
  movsd  %xmm9,0x38(%rsp)
  movsd  %xmm15,0x30(%rsp)
  movsd  %xmm14,0x28(%rsp)
  movsd  %xmm2,0x20(%rsp)
  call   <sqrt@plt>
  mov    D37(%rip),%rax
  pxor   %xmm1,%xmm1
  movsd  0x160(%rsp),%xmm7
  movsd  0x58(%rsp),%xmm6
  movsd  0x50(%rsp),%xmm5
  movapd %xmm0,%xmm11
  movsd  0x48(%rsp),%xmm4
  movsd  0x40(%rsp),%xmm3
  movq   %rax,%xmm10
  movsd  0x38(%rsp),%xmm9
  movsd  0x30(%rsp),%xmm15
  movsd  0x28(%rsp),%xmm14
  movsd  0x20(%rsp),%xmm2
  jmp    L90
L348:
  movapd %xmm7,%xmm0
  movsd  %xmm6,0x60(%rsp)
  movsd  %xmm5,0x58(%rsp)
  movsd  %xmm4,0x50(%rsp)
  movsd  %xmm3,0x48(%rsp)
  movsd  %xmm2,0x40(%rsp)
  movsd  %xmm15,0x38(%rsp)
  movsd  %xmm14,0x30(%rsp)
  movsd  %xmm8,0x28(%rsp)
  movsd  %xmm9,0x20(%rsp)
  movsd  %xmm7,0x160(%rsp)
  call   <sqrt@plt>
  mov    D37(%rip),%rax
  pxor   %xmm1,%xmm1
  movsd  0x60(%rsp),%xmm6
  movsd  0x58(%rsp),%xmm5
  movsd  0x50(%rsp),%xmm4
  movsd  0x48(%rsp),%xmm3
  movsd  0x40(%rsp),%xmm2
  movq   %rax,%xmm10
  movsd  0x38(%rsp),%xmm15
  movsd  0x30(%rsp),%xmm14
  movsd  0x28(%rsp),%xmm8
  movsd  0x20(%rsp),%xmm9
  movsd  0x160(%rsp),%xmm7
  jmp    L319
L286:
  movapd %xmm2,%xmm0
  movsd  %xmm6,0x58(%rsp)
  movsd  %xmm7,0x160(%rsp)
  movsd  %xmm5,0x50(%rsp)
  movsd  %xmm4,0x48(%rsp)
  movsd  %xmm3,0x40(%rsp)
  movsd  %xmm15,0x38(%rsp)
  movsd  %xmm14,0x30(%rsp)
  movsd  %xmm9,0x28(%rsp)
  movsd  %xmm2,0x20(%rsp)
  call   <sqrt@plt>
  mov    D37(%rip),%rax
  pxor   %xmm1,%xmm1
  movsd  0x160(%rsp),%xmm7
  movsd  0x58(%rsp),%xmm6
  movsd  0x50(%rsp),%xmm5
  movapd %xmm0,%xmm8
  movsd  0x48(%rsp),%xmm4
  movsd  0x40(%rsp),%xmm3
  movq   %rax,%xmm10
  movsd  0x38(%rsp),%xmm15
  movsd  0x30(%rsp),%xmm14
  movsd  0x28(%rsp),%xmm9
  movsd  0x20(%rsp),%xmm2
  jmp    L164
L117:
  movapd %xmm7,%xmm0
  movsd  %xmm6,0x58(%rsp)
  movsd  %xmm5,0x50(%rsp)
  movsd  %xmm4,0x48(%rsp)
  movsd  %xmm3,0x40(%rsp)
  movsd  %xmm2,0x38(%rsp)
  movsd  %xmm15,0x30(%rsp)
  movsd  %xmm14,0x28(%rsp)
  movsd  %xmm8,0x20(%rsp)
  movsd  %xmm7,0x160(%rsp)
  call   <sqrt@plt>
  mov    D37(%rip),%rax
  pxor   %xmm1,%xmm1
  movsd  0x58(%rsp),%xmm6
  movsd  0x50(%rsp),%xmm5
  movsd  0x48(%rsp),%xmm4
  movapd %xmm0,%xmm9
  movsd  0x40(%rsp),%xmm3
  movsd  0x38(%rsp),%xmm2
  movq   %rax,%xmm10
  movsd  0x30(%rsp),%xmm15
  movsd  0x28(%rsp),%xmm14
  movsd  0x160(%rsp),%xmm7
  movsd  0x20(%rsp),%xmm8
  jmp    L222
L282:
  movapd %xmm2,%xmm0
  movsd  %xmm6,0x50(%rsp)
  movsd  %xmm7,0x160(%rsp)
  movsd  %xmm5,0x48(%rsp)
  movsd  %xmm4,0x40(%rsp)
  movsd  %xmm3,0x38(%rsp)
  movsd  %xmm15,0x30(%rsp)
  movsd  %xmm14,0x28(%rsp)
  movsd  %xmm2,0x20(%rsp)
  call   <sqrt@plt>
  mov    D37(%rip),%rax
  pxor   %xmm1,%xmm1
  movsd  0x160(%rsp),%xmm7
  movsd  0x50(%rsp),%xmm6
  movsd  0x48(%rsp),%xmm5
  movapd %xmm0,%xmm8
  movsd  0x40(%rsp),%xmm4
  movsd  0x38(%rsp),%xmm3
  movq   %rax,%xmm10
  movsd  0x30(%rsp),%xmm15
  movsd  0x28(%rsp),%xmm14
  movsd  0x20(%rsp),%xmm2
  jmp    L43
```

tool_calls: name: extract_symbolic_data
arguments: data_tag: D145
data_type: f64
name: extract_symbolic_data
arguments: data_tag: D37
data_type: f64
name: extract_symbolic_data
arguments: data_tag: D78
data_type: f32[4]
name: extract_symbolic_data
arguments: data_tag: D169
data_type: qword
name: extract_symbolic_data
arguments: data_tag: D142
data_type: f32[2]

tool: D145: 1

tool: D37: 0.5

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

tool: D169: 9223372036854775807

tool: {"D142": [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];
    }
  }
}
```