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

```asm
<gcproj>:
  endbr64
  push   %r15
  push   %r14
  mov    %rdx,%r14
  push   %r13
  mov    %rcx,%r13
  push   %r12
  mov    %rdi,%r12
  push   %rbp
  mov    %r9,%rbp
  push   %rbx
  mov    %rsi,%rbx
  sub    $0x48,%rsp
  mov    0x98(%rsp),%edx
  mov    0x88(%rsp),%r15
  mov    0x90(%rsp),%rax
  test   %edx,%edx
  jne    L5
  movss  (%r8),%xmm1
  movss  (%r12),%xmm0
  lea    0x38(%rsp),%rdi
  lea    0x30(%rsp),%rsi
  mov    0x80(%rsp),%rax
  mov    %rsi,0x18(%rsp)
  divss  %xmm1,%xmm0
  mov    %rdi,0x10(%rsp)
  movss  %xmm1,0x20(%rsp)
  cvtss2sd %xmm0,%xmm0
  subsd  (%rax),%xmm0
  call   <sincos@plt>
  movss  0x20(%rsp),%xmm1
  movss  (%rbx),%xmm0
  movsd  0x38(%rsp),%xmm3
  movsd  0x30(%rsp),%xmm7
  divss  %xmm1,%xmm0
  mov    0x18(%rsp),%rsi
  mov    0x10(%rsp),%rdi
  movsd  %xmm3,0x8(%rsp)
  movsd  %xmm7,(%rsp)
  cvtss2sd %xmm0,%xmm0
  subsd  0x0(%rbp),%xmm0
  call   <sincos@plt>
  movsd  0x8(%rsp),%xmm3
  movsd  0x38(%rsp),%xmm2
  pxor   %xmm4,%xmm4
  movsd  D26(%rip),%xmm1
  movsd  0x30(%rsp),%xmm5
  movapd %xmm3,%xmm0
  mulsd  %xmm3,%xmm0
  mulsd  %xmm2,%xmm0
  mulsd  %xmm2,%xmm0
  addsd  %xmm1,%xmm0
  ucomisd %xmm0,%xmm4
  ja     L90
  sqrtsd %xmm0,%xmm0
L23:
  mulsd  (%rsp),%xmm2
  mulsd  %xmm5,%xmm3
  mulsd  %xmm0,%xmm2
  mulsd  %xmm0,%xmm3
  movapd %xmm2,%xmm5
  mulsd  %xmm2,%xmm5
  subsd  %xmm5,%xmm1
  movapd %xmm1,%xmm0
  movapd %xmm3,%xmm1
  mulsd  %xmm3,%xmm1
  subsd  %xmm1,%xmm0
  ucomisd %xmm0,%xmm4
  ja     L14
  sqrtsd %xmm0,%xmm0
L40:
  movsd  (%r15),%xmm1
  movsd  0x8(%r15),%xmm5
  movsd  0x20(%r15),%xmm6
  mulsd  %xmm3,%xmm5
  mulsd  %xmm2,%xmm1
  mulsd  %xmm3,%xmm6
  mulsd  0x38(%r15),%xmm3
  addsd  %xmm5,%xmm1
  movsd  0x10(%r15),%xmm5
  mulsd  %xmm0,%xmm5
  addsd  %xmm5,%xmm1
  movsd  0x18(%r15),%xmm5
  mulsd  %xmm2,%xmm5
  mulsd  0x30(%r15),%xmm2
  addsd  %xmm6,%xmm5
  movsd  0x28(%r15),%xmm6
  addsd  %xmm3,%xmm2
  mulsd  %xmm0,%xmm6
  mulsd  0x40(%r15),%xmm0
  addsd  %xmm6,%xmm5
  addsd  %xmm0,%xmm2
  ucomisd %xmm4,%xmm2
  jp     L95
  jne    L95
  movl   $0x0,0x0(%r13)
L115:
  ucomisd %xmm4,%xmm1
  jp     L140
  pxor   %xmm0,%xmm0
  je     L85
L140:
  divsd  %xmm1,%xmm5
  movsd  %xmm4,0x8(%rsp)
  movsd  %xmm1,(%rsp)
  movapd %xmm5,%xmm0
  call   <atan@plt>
  movsd  0x8(%rsp),%xmm4
  movsd  (%rsp),%xmm1
  divsd  D91(%rip),%xmm0
  cvtsd2ss %xmm0,%xmm0
L85:
  comisd %xmm1,%xmm4
  jbe    L89
  subss  D24(%rip),%xmm0
L89:
  movss  D55(%rip),%xmm1
  movss  D52(%rip),%xmm2
  movss  %xmm0,(%r14)
  comiss %xmm0,%xmm1
  jbe    L195
  nopw   0x0(%rax,%rax,1)
L124:
  addss  %xmm2,%xmm0
  comiss %xmm0,%xmm1
  ja     L124
  movss  %xmm0,(%r14)
  add    $0x48,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
  nopl   (%rax)
L5:
  pxor   %xmm0,%xmm0
  lea    0x38(%rsp),%r15
  lea    0x30(%rsp),%rsi
  mov    %rax,0x90(%rsp)
  cvtss2sd (%r14),%xmm0
  mulsd  D91(%rip),%xmm0
  mov    %r15,%rdi
  mov    %r8,0x18(%rsp)
  mov    %rsi,0x10(%rsp)
  call   <sincos@plt>
  movsd  D82(%rip),%xmm0
  pxor   %xmm1,%xmm1
  movsd  0x30(%rsp),%xmm2
  cvtss2sd 0x0(%r13),%xmm1
  movsd  0x38(%rsp),%xmm7
  mov    0x10(%rsp),%rsi
  mov    %r15,%rdi
  subsd  %xmm1,%xmm0
  mulsd  D91(%rip),%xmm0
  movsd  %xmm2,0x8(%rsp)
  movsd  %xmm7,(%rsp)
  call   <sincos@plt>
  movsd  0x38(%rsp),%xmm0
  mov    0x18(%rsp),%r8
  movsd  0x8(%rsp),%xmm2
  mov    0x90(%rsp),%rax
  movsd  0x30(%rsp),%xmm1
  mulsd  %xmm0,%xmm2
  movsd  0x8(%rax),%xmm4
  mulsd  (%rsp),%xmm0
  movapd %xmm2,%xmm3
  mulsd  %xmm0,%xmm4
  mulsd  (%rax),%xmm3
  mulsd  0x20(%rax),%xmm0
  mulsd  0x18(%rax),%xmm2
  addsd  %xmm4,%xmm3
  movsd  0x10(%rax),%xmm4
  mulsd  %xmm1,%xmm4
  addsd  %xmm0,%xmm2
  mulsd  0x28(%rax),%xmm1
  addsd  %xmm4,%xmm3
  addsd  %xmm1,%xmm2
  movsd  D26(%rip),%xmm1
  movapd %xmm1,%xmm0
  movapd %xmm2,%xmm4
  mulsd  %xmm2,%xmm4
  subsd  %xmm4,%xmm0
  pxor   %xmm4,%xmm4
  ucomisd %xmm0,%xmm4
  ja     L17
  sqrtsd %xmm0,%xmm0
L80:
  movapd %xmm3,%xmm5
  mulsd  %xmm3,%xmm3
  divsd  %xmm0,%xmm5
  subsd  %xmm3,%xmm1
  ucomisd %xmm1,%xmm4
  ja     L41
  sqrtsd %xmm1,%xmm1
L121:
  pxor   %xmm7,%xmm7
  mov    %r8,0x8(%rsp)
  cvtss2sd (%r8),%xmm7
  movsd  %xmm5,0x10(%rsp)
  movsd  %xmm7,(%rsp)
  divsd  %xmm1,%xmm2
  movapd %xmm2,%xmm0
  call   <asin@plt>
  movsd  0x10(%rsp),%xmm5
  mov    0x8(%rsp),%r8
  pxor   %xmm6,%xmm6
  mov    0x80(%rsp),%rax
  addsd  (%rax),%xmm0
  mulsd  (%rsp),%xmm0
  cvtsd2ss %xmm0,%xmm0
  movss  %xmm0,(%r12)
  movapd %xmm5,%xmm0
  cvtss2sd (%r8),%xmm6
  movsd  %xmm6,(%rsp)
  call   <asin@plt>
  addsd  0x0(%rbp),%xmm0
  mulsd  (%rsp),%xmm0
  cvtsd2ss %xmm0,%xmm0
  movss  %xmm0,(%rbx)
L195:
  add    $0x48,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
  nopl   0x0(%rax)
L95:
  movapd %xmm5,%xmm3
  movapd %xmm1,%xmm0
  mulsd  %xmm5,%xmm3
  mulsd  %xmm1,%xmm0
  addsd  %xmm3,%xmm0
  ucomisd %xmm0,%xmm4
  ja     L103
  sqrtsd %xmm0,%xmm0
L78:
  divsd  %xmm2,%xmm0
  movsd  %xmm4,0x18(%rsp)
  movsd  %xmm5,0x10(%rsp)
  movsd  %xmm1,0x8(%rsp)
  movsd  %xmm2,(%rsp)
  call   <atan@plt>
  movsd  (%rsp),%xmm2
  movsd  D82(%rip),%xmm3
  divsd  D91(%rip),%xmm0
  movsd  0x18(%rsp),%xmm4
  movsd  0x8(%rsp),%xmm1
  subsd  %xmm0,%xmm3
  movsd  0x10(%rsp),%xmm5
  comisd %xmm2,%xmm4
  cvtsd2ss %xmm3,%xmm3
  jbe    L167
  subss  D24(%rip),%xmm3
L167:
  movss  %xmm3,0x0(%r13)
  jmp    L115
L90:
  movsd  %xmm4,0x28(%rsp)
  movsd  %xmm1,0x20(%rsp)
  movsd  %xmm3,0x18(%rsp)
  movsd  %xmm2,0x10(%rsp)
  movsd  %xmm5,0x8(%rsp)
  call   <sqrt@plt>
  movsd  0x28(%rsp),%xmm4
  movsd  0x20(%rsp),%xmm1
  movsd  0x18(%rsp),%xmm3
  movsd  0x10(%rsp),%xmm2
  movsd  0x8(%rsp),%xmm5
  jmp    L23
L17:
  mov    %r8,0x10(%rsp)
  movsd  %xmm4,0x20(%rsp)
  movsd  %xmm1,0x18(%rsp)
  movsd  %xmm2,0x8(%rsp)
  movsd  %xmm3,(%rsp)
  call   <sqrt@plt>
  movsd  0x20(%rsp),%xmm4
  mov    0x10(%rsp),%r8
  movsd  0x18(%rsp),%xmm1
  movsd  0x8(%rsp),%xmm2
  movsd  (%rsp),%xmm3
  jmp    L80
L14:
  movsd  %xmm4,0x10(%rsp)
  movsd  %xmm3,0x8(%rsp)
  movsd  %xmm2,(%rsp)
  call   <sqrt@plt>
  movsd  0x10(%rsp),%xmm4
  movsd  0x8(%rsp),%xmm3
  movsd  (%rsp),%xmm2
  jmp    L40
L41:
  movapd %xmm1,%xmm0
  mov    %r8,0x10(%rsp)
  movsd  %xmm5,0x8(%rsp)
  movsd  %xmm2,(%rsp)
  call   <sqrt@plt>
  mov    0x10(%rsp),%r8
  movsd  0x8(%rsp),%xmm5
  movsd  (%rsp),%xmm2
  movapd %xmm0,%xmm1
  jmp    L121
L103:
  movsd  %xmm4,0x18(%rsp)
  movsd  %xmm2,0x10(%rsp)
  movsd  %xmm5,0x8(%rsp)
  movsd  %xmm1,(%rsp)
  call   <sqrt@plt>
  movsd  0x18(%rsp),%xmm4
  movsd  0x10(%rsp),%xmm2
  movsd  0x8(%rsp),%xmm5
  movsd  (%rsp),%xmm1
  jmp    L78
```

tool_calls: name: extract_symbolic_data
arguments: tag: D26
data_type: f64
name: extract_symbolic_data
arguments: tag: D91
data_type: f64
name: extract_symbolic_data
arguments: tag: D24
data_type: f32
name: extract_symbolic_data
arguments: tag: D55
data_type: byte[4]
name: extract_symbolic_data
arguments: tag: D52
data_type: f32
name: extract_symbolic_data
arguments: tag: D82
data_type: f64

tool: D26: 1

tool: D91: 0.017453292

tool: D24: 180

tool: {"D55": [0x00,0x00,0x34,0xc3]}

tool: D52: 360

tool: D82: 90

assistant: ```c
void gcproj(float *xf, float *yf, float *rlon, float *rlat, float *ref_rad, double *g0, double *b0, double *amat, double *ainv, int gflag) {
  double xp, yp, zp;
  double xg, yg, zg;
  double arg;
  double cosG, sinG;
  double cosB, sinB;
  double rperd = 0.017453292;
  if (gflag == 0) {
    arg = (*xf) / (*ref_rad) - (*b0);
    cosB = cos(arg);
    sinB = sin(arg);
    arg = (*yf) / (*ref_rad) - (*g0);
    cosG = cos(arg);
    sinG = sin(arg);
    arg = sqrt(1.0 + sinB * sinB * sinG * sinG);
    xp = sinG * cosB * arg;
    yp = sinB * cosG * arg;
    zp = sqrt(1.0 - xp * xp - yp * yp);
    xg = xp * amat[0] + yp * amat[1] + zp * amat[2];
    yg = xp * amat[3] + yp * amat[4] + zp * amat[5];
    zg = xp * amat[6] + yp * amat[7] + zp * amat[8];
    if (zg != (double)(0.0)) {
      arg = sqrt(xg * xg + yg * yg) / zg;
      (*rlat) = 90.0 - atan(arg) / rperd;
      if (zg < (double)(0.0))
        (*rlat) = (*rlat) - 180.0;
    } else
      (*rlat) = 0.0;
    if (xg != (double)(0.0)) {
      arg = yg / xg;
      (*rlon) = atan(arg) / rperd;
    } else
      (*rlon) = 0.0;
    if (xg < (double)(0.0))
      (*rlon) = (*rlon) - 180.0;
    while ((*rlon) < (double)(-180.0))
      (*rlon) = (*rlon) + 360.0;
  } else {
    arg = (*rlon) * rperd;
    cosG = cos(arg);
    sinG = sin(arg);
    arg = (90.0 - (*rlat)) * rperd;
    cosB = cos(arg);
    sinB = sin(arg);
    xg = sinB * cosG;
    yg = sinB * sinG;
    zg = cosB;
    xp = xg * ainv[0] + yg * ainv[1] + zg * ainv[2];
    yp = xg * ainv[3] + yg * ainv[4] + zg * ainv[5];
    zp = xg * ainv[6] + yg * ainv[7] + zg * ainv[8];
    sinG = xp / sqrt(1.0 - yp * yp);
    sinB = yp / sqrt(1.0 - xp * xp);
    *xf = (double)(*ref_rad) * (asin(sinB) + (*b0));
    *yf = (double)(*ref_rad) * (asin(sinG) + (*g0));
  }
}
```