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

```asm
<pp_sin>:
  endbr64
  push   %rbp
  mov    %rsp,%rbp
  movsd  %xmm0,-0x38(%rbp)
  movsd  -0x38(%rbp),%xmm0
  movsd  %xmm0,-0x20(%rbp)
  movl   $0x1,-0x24(%rbp)
  jmp    L41
L48:
  movsd  -0x20(%rbp),%xmm1
  movsd  D98(%rip),%xmm0
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x20(%rbp)
L41:
  movsd  D54(%rip),%xmm0
  comisd -0x20(%rbp),%xmm0
  ja     L48
  jmp    L56
L49:
  movsd  -0x20(%rbp),%xmm0
  movsd  D98(%rip),%xmm1
  subsd  %xmm1,%xmm0
  movsd  %xmm0,-0x20(%rbp)
L56:
  movsd  -0x20(%rbp),%xmm0
  comisd D114(%rip),%xmm0
  ja     L49
  movsd  -0x20(%rbp),%xmm0
  comisd D38(%rip),%xmm0
  jbe    L33
  movsd  -0x20(%rbp),%xmm0
  movsd  D28(%rip),%xmm1
  subsd  %xmm1,%xmm0
  movsd  %xmm0,-0x20(%rbp)
  movl   $0xffffffff,-0x24(%rbp)
L33:
  movsd  -0x20(%rbp),%xmm0
  movsd  D99(%rip),%xmm1
  divsd  %xmm1,%xmm0
  movsd  %xmm0,-0x18(%rbp)
  movsd  -0x18(%rbp),%xmm0
  mulsd  %xmm0,%xmm0
  movsd  %xmm0,-0x10(%rbp)
  movsd  -0x10(%rbp),%xmm1
  movsd  D108(%rip),%xmm0
  mulsd  %xmm1,%xmm0
  movsd  D71(%rip),%xmm1
  subsd  %xmm1,%xmm0
  movapd %xmm0,%xmm1
  mulsd  -0x10(%rbp),%xmm1
  movsd  D25(%rip),%xmm0
  addsd  %xmm1,%xmm0
  mulsd  -0x10(%rbp),%xmm0
  movsd  D49(%rip),%xmm1
  subsd  %xmm1,%xmm0
  movapd %xmm0,%xmm1
  mulsd  -0x10(%rbp),%xmm1
  movsd  D109(%rip),%xmm0
  addsd  %xmm1,%xmm0
  movsd  -0x18(%rbp),%xmm1
  mulsd  %xmm1,%xmm0
  movsd  %xmm0,-0x8(%rbp)
  cmpl   $0x0,-0x24(%rbp)
  jns    L34
  movsd  -0x8(%rbp),%xmm0
  movq   D21(%rip),%xmm1
  xorpd  %xmm1,%xmm0
  jmp    L14
L34:
  movsd  -0x8(%rbp),%xmm0
L14:
  movq   %xmm0,%rax
  movq   %rax,%xmm0
  pop    %rbp
  ret
```

tool_calls: name: disassemble_data
arguments: tag: D98
datatype: f64
name: disassemble_data
arguments: tag: D54
datatype: byte[8]
name: disassemble_data
arguments: tag: D114
datatype: f64[2]
name: disassemble_data
arguments: tag: D38
datatype: f64[2]
name: disassemble_data
arguments: tag: D28
datatype: f64
name: disassemble_data
arguments: tag: D99
datatype: f64
name: disassemble_data
arguments: tag: D108
datatype: byte[8]
name: disassemble_data
arguments: tag: D71
datatype: byte[8]
name: disassemble_data
arguments: tag: D25
datatype: byte[8]
name: disassemble_data
arguments: tag: D49
datatype: byte[8]
name: disassemble_data
arguments: tag: D109
datatype: f64
name: disassemble_data
arguments: tag: D21
datatype: byte[8]

tool: D98: 6.283185307179586

tool: {"D54": [0x18,0x2d,0x44,0x54,0xfb,0x21,0xf9,0xbf]}

tool: {"D114": [4.71238898038469,1.5707963267948966]}

tool: {"D38": [1.5707963267948966,3.141592653589793]}

tool: D28: 3.141592653589793

tool: D99: 1.57079632679

tool: {"D108": [0xaf,0x90,0xa8,0x48,0xf7,0xda,0x23,0x3f]}

tool: {"D71": [0x0e,0x6d,0x5b,0x64,0xcc,0x24,0x73,0x3f]}

tool: {"D25": [0xa5,0x19,0x56,0xf6,0x8a,0x66,0xb4,0x3f]}

tool: {"D49": [0xcd,0xf9,0xac,0x16,0xbc,0xab,0xe4,0x3f]}

tool: D109: 1.57079631847

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

assistant: ```c
double pp_sin(double x) {
  double xi, y, q, q2;
  int sign;
  xi = x;
  sign = 1;
  while (xi < -1.57079632679489661923)
    xi += 6.28318530717958647692;
  while (xi > 4.71238898038468985769)
    xi -= 6.28318530717958647692;
  if (xi > 1.57079632679489661923) {
    xi -= 3.141592653589793238462643;
    sign = -1;
  }
  q = xi / 1.57079632679;
  q2 = q * q;
  y = ((((.00015148419 * q2 - .00467376557) * q2 + .07968967928) * q2 - .64596371106) * q2 + 1.57079631847) * q;
  return (sign < 0 ? -y : y);
}
```