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

```asm
<turbulence>:
  endbr64
  push   %r15
  mov    %edi,%eax
  push   %r14
  sub    %edx,%eax
  push   %r13
  mov    %esi,%r13d
  push   %r12
  mov    %edi,%r12d
  push   %rbp
  lea    (%r12,%r13,1),%ebp
  push   %rbx
  mov    %edx,%ebx
  lea    0x0(%rbp,%r12,1),%r15d
  mov    %ebx,%r14d
  sub    $0x28,%rsp
  movsd  %xmm0,0x18(%rsp)
  pxor   %xmm0,%xmm0
  cvtsi2sd %eax,%xmm0
  movsd  %xmm1,0x8(%rsp)
  mulsd  %xmm1,%xmm0
  divsd  D74(%rip),%xmm0
  cvttsd2si %xmm0,%edi
  call   <sin@plt>
  imul   $0xffffffe9,%ebx,%eax
  movsd  %xmm0,0x10(%rsp)
  pxor   %xmm0,%xmm0
  add    %r13d,%eax
  cvtsi2sd %eax,%xmm0
  mulsd  0x8(%rsp),%xmm0
  divsd  D74(%rip),%xmm0
  cvttsd2si %xmm0,%edi
  call   <sin@plt>
  imul   $0xffffffed,%ebx,%eax
  movsd  0x10(%rsp),%xmm2
  addsd  %xmm0,%xmm2
  pxor   %xmm0,%xmm0
  add    %ebp,%eax
  add    %r13d,%ebp
  cvtsi2sd %eax,%xmm0
  mulsd  0x8(%rsp),%xmm0
  divsd  D74(%rip),%xmm0
  cvttsd2si %xmm0,%edi
  movsd  %xmm2,0x10(%rsp)
  call   <sin@plt>
  imul   $0xffffffeb,%ebx,%eax
  movsd  0x10(%rsp),%xmm3
  addsd  %xmm0,%xmm3
  pxor   %xmm0,%xmm0
  add    %r15d,%eax
  add    %r12d,%r15d
  cvtsi2sd %eax,%xmm0
  mulsd  0x8(%rsp),%xmm0
  divsd  D74(%rip),%xmm0
  cvttsd2si %xmm0,%edi
  movsd  %xmm3,0x10(%rsp)
  call   <sin@plt>
  movsd  0x10(%rsp),%xmm4
  lea    0x0(,%rbx,8),%eax
  sub    %eax,%r14d
  imul   $0xfffffff5,%ebx,%ebx
  lea    0x0(%rbp,%r14,2),%eax
  addsd  %xmm0,%xmm4
  pxor   %xmm0,%xmm0
  add    %r13d,%ebp
  cvtsi2sd %eax,%xmm0
  mulsd  0x8(%rsp),%xmm0
  add    %r14d,%ebp
  divsd  D74(%rip),%xmm0
  cvttsd2si %xmm0,%edi
  add    %ebx,%r15d
  movsd  %xmm4,0x10(%rsp)
  call   <sin@plt>
  movsd  0x10(%rsp),%xmm5
  addsd  %xmm0,%xmm5
  pxor   %xmm0,%xmm0
  cvtsi2sd %r15d,%xmm0
  mulsd  0x8(%rsp),%xmm0
  divsd  D74(%rip),%xmm0
  cvttsd2si %xmm0,%edi
  movsd  %xmm5,0x10(%rsp)
  call   <sin@plt>
  movsd  0x10(%rsp),%xmm6
  addsd  %xmm0,%xmm6
  pxor   %xmm0,%xmm0
  cvtsi2sd %ebp,%xmm0
  mulsd  0x8(%rsp),%xmm0
  divsd  D74(%rip),%xmm0
  cvttsd2si %xmm0,%edi
  movsd  %xmm6,0x10(%rsp)
  call   <sin@plt>
  lea    0x0(,%r12,8),%eax
  movsd  0x10(%rsp),%xmm7
  sub    %r12d,%eax
  lea    (%rax,%r13,8),%eax
  addsd  %xmm0,%xmm7
  pxor   %xmm0,%xmm0
  add    %r14d,%eax
  cvtsi2sd %eax,%xmm0
  mulsd  0x8(%rsp),%xmm0
  divsd  D74(%rip),%xmm0
  cvttsd2si %xmm0,%edi
  movsd  %xmm7,0x10(%rsp)
  call   <sin@plt>
  addsd  0x10(%rsp),%xmm0
  mulsd  D37(%rip),%xmm0
  mulsd  0x18(%rsp),%xmm0
  add    $0x28,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
```

tool_calls: name: parse_data
arguments: label: D74
type: f64
name: parse_data
arguments: label: D37
type: byte[8]

tool: D74: 3.141592654

tool: {"D37": [0x00,0x00,0x00,0x00,0x00,0x00,0xc0,0x3f]}

assistant: ```c
double turbulence(int x, int y, int seed, double strength, double frequency) {
  double offset = sin((x - seed) * frequency / (3.141592654)) + sin((y - seed * 23) * frequency / (3.141592654)) + sin((x + y - seed * 19) * frequency / (3.141592654)) + sin((2 * x + y - seed * 21) * frequency / (3.141592654)) + sin((x + 2 * y - seed * 14) * frequency / (3.141592654)) + sin((3 * x + y - seed * 11) * frequency / (3.141592654)) + sin((x + 3 * y - seed * 7) * frequency / (3.141592654)) + sin((7 * x + 8 * y - seed * 7) * frequency / (3.141592654));
  offset /= 8;
  return strength * offset;
}
```