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

```asm
<imc_rotate_WCS>:
  endbr64
  push   %rbp
  mov    %rsp,%rbp
  push   %rbx
  sub    $0x58,%rsp
  mov    %rdi,-0x48(%rbp)
  movss  %xmm0,-0x4c(%rbp)
  movss  %xmm1,-0x50(%rbp)
  movsd  %xmm2,-0x58(%rbp)
  movsd  %xmm3,-0x60(%rbp)
  mov    %fs:0x28,%rax
  mov    %rax,-0x18(%rbp)
  xor    %eax,%eax
  lea    -0x40(%rbp),%rdx
  mov    -0x48(%rbp),%rax
  lea    D89(%rip),%rcx
  mov    %rcx,%rsi
  mov    %rax,%rdi
  call   <imget_fits_float_value@plt>
  cmp    $0x1,%eax
  jne    L50
  lea    -0x3c(%rbp),%rdx
  mov    -0x48(%rbp),%rax
  lea    D258(%rip),%rcx
  mov    %rcx,%rsi
  mov    %rax,%rdi
  call   <imget_fits_float_value@plt>
  cmp    $0x1,%eax
  jne    L50
  movss  -0x40(%rbp),%xmm0
  cvtss2sd %xmm0,%xmm0
  mulsd  -0x58(%rbp),%xmm0
  movss  -0x3c(%rbp),%xmm1
  cvtss2sd %xmm1,%xmm1
  movapd %xmm1,%xmm2
  mulsd  -0x60(%rbp),%xmm2
  movapd %xmm0,%xmm1
  subsd  %xmm2,%xmm1
  pxor   %xmm0,%xmm0
  cvtss2sd -0x4c(%rbp),%xmm0
  addsd  %xmm1,%xmm0
  pxor   %xmm4,%xmm4
  cvtsd2ss %xmm0,%xmm4
  movd   %xmm4,%edx
  mov    -0x48(%rbp),%rax
  lea    D398(%rip),%rcx
  movd   %edx,%xmm0
  lea    D191(%rip),%rdx
  lea    D89(%rip),%rsi
  mov    %rax,%rdi
  call   <imupdate_fitsf@plt>
  movss  -0x40(%rbp),%xmm0
  cvtss2sd %xmm0,%xmm0
  movapd %xmm0,%xmm1
  mulsd  -0x60(%rbp),%xmm1
  movss  -0x3c(%rbp),%xmm0
  cvtss2sd %xmm0,%xmm0
  mulsd  -0x58(%rbp),%xmm0
  addsd  %xmm0,%xmm1
  pxor   %xmm0,%xmm0
  cvtss2sd -0x50(%rbp),%xmm0
  addsd  %xmm1,%xmm0
  pxor   %xmm5,%xmm5
  cvtsd2ss %xmm0,%xmm5
  movd   %xmm5,%edx
  mov    -0x48(%rbp),%rax
  lea    D398(%rip),%rcx
  movd   %edx,%xmm0
  lea    D191(%rip),%rdx
  lea    D258(%rip),%rsi
  mov    %rax,%rdi
  call   <imupdate_fitsf@plt>
  lea    -0x38(%rbp),%rdx
  mov    -0x48(%rbp),%rax
  lea    D311(%rip),%rcx
  mov    %rcx,%rsi
  mov    %rax,%rdi
  call   <imget_fits_float_value@plt>
  mov    %eax,%ebx
  lea    -0x34(%rbp),%rdx
  mov    -0x48(%rbp),%rax
  lea    D202(%rip),%rcx
  mov    %rcx,%rsi
  mov    %rax,%rdi
  call   <imget_fits_float_value@plt>
  add    %eax,%ebx
  lea    -0x30(%rbp),%rdx
  mov    -0x48(%rbp),%rax
  lea    D108(%rip),%rcx
  mov    %rcx,%rsi
  mov    %rax,%rdi
  call   <imget_fits_float_value@plt>
  add    %eax,%ebx
  lea    -0x2c(%rbp),%rdx
  mov    -0x48(%rbp),%rax
  lea    D390(%rip),%rcx
  mov    %rcx,%rsi
  mov    %rax,%rdi
  call   <imget_fits_float_value@plt>
  add    %ebx,%eax
  test   %eax,%eax
  je     L2
  movss  -0x38(%rbp),%xmm0
  cvtss2sd %xmm0,%xmm0
  mulsd  -0x58(%rbp),%xmm0
  movss  -0x30(%rbp),%xmm1
  cvtss2sd %xmm1,%xmm1
  mulsd  -0x60(%rbp),%xmm1
  subsd  %xmm1,%xmm0
  pxor   %xmm6,%xmm6
  cvtsd2ss %xmm0,%xmm6
  movd   %xmm6,%edx
  mov    -0x48(%rbp),%rax
  lea    D398(%rip),%rcx
  movd   %edx,%xmm0
  lea    D210(%rip),%rdx
  lea    D311(%rip),%rsi
  mov    %rax,%rdi
  call   <imupdate_fitsf@plt>
  movss  -0x34(%rbp),%xmm0
  cvtss2sd %xmm0,%xmm0
  mulsd  -0x58(%rbp),%xmm0
  movss  -0x2c(%rbp),%xmm1
  cvtss2sd %xmm1,%xmm1
  mulsd  -0x60(%rbp),%xmm1
  subsd  %xmm1,%xmm0
  pxor   %xmm7,%xmm7
  cvtsd2ss %xmm0,%xmm7
  movd   %xmm7,%edx
  mov    -0x48(%rbp),%rax
  lea    D398(%rip),%rcx
  movd   %edx,%xmm0
  lea    D210(%rip),%rdx
  lea    D202(%rip),%rsi
  mov    %rax,%rdi
  call   <imupdate_fitsf@plt>
  movss  -0x38(%rbp),%xmm0
  cvtss2sd %xmm0,%xmm0
  movapd %xmm0,%xmm1
  mulsd  -0x60(%rbp),%xmm1
  movss  -0x30(%rbp),%xmm0
  cvtss2sd %xmm0,%xmm0
  mulsd  -0x58(%rbp),%xmm0
  addsd  %xmm1,%xmm0
  pxor   %xmm3,%xmm3
  cvtsd2ss %xmm0,%xmm3
  movd   %xmm3,%edx
  mov    -0x48(%rbp),%rax
  lea    D398(%rip),%rcx
  movd   %edx,%xmm0
  lea    D210(%rip),%rdx
  lea    D108(%rip),%rsi
  mov    %rax,%rdi
  call   <imupdate_fitsf@plt>
  movss  -0x34(%rbp),%xmm0
  cvtss2sd %xmm0,%xmm0
  movapd %xmm0,%xmm1
  mulsd  -0x60(%rbp),%xmm1
  movss  -0x2c(%rbp),%xmm0
  cvtss2sd %xmm0,%xmm0
  mulsd  -0x58(%rbp),%xmm0
  addsd  %xmm1,%xmm0
  pxor   %xmm4,%xmm4
  cvtsd2ss %xmm0,%xmm4
  movd   %xmm4,%edx
  mov    -0x48(%rbp),%rax
  lea    D398(%rip),%rcx
  movd   %edx,%xmm0
  lea    D210(%rip),%rdx
  lea    D390(%rip),%rsi
  mov    %rax,%rdi
  call   <imupdate_fitsf@plt>
L2:
  lea    -0x28(%rbp),%rdx
  mov    -0x48(%rbp),%rax
  lea    D410(%rip),%rcx
  mov    %rcx,%rsi
  mov    %rax,%rdi
  call   <imget_fits_float_value@plt>
  cmp    $0x1,%eax
  jne    L48
  lea    -0x24(%rbp),%rdx
  mov    -0x48(%rbp),%rax
  lea    D56(%rip),%rcx
  mov    %rcx,%rsi
  mov    %rax,%rdi
  call   <imget_fits_float_value@plt>
  cmp    $0x1,%eax
  jne    L48
  lea    -0x20(%rbp),%rdx
  mov    -0x48(%rbp),%rax
  lea    D54(%rip),%rcx
  mov    %rcx,%rsi
  mov    %rax,%rdi
  call   <imget_fits_float_value@plt>
  cmp    $0x1,%eax
  jne    L48
  lea    -0x1c(%rbp),%rdx
  mov    -0x48(%rbp),%rax
  lea    D313(%rip),%rcx
  mov    %rcx,%rsi
  mov    %rax,%rdi
  call   <imget_fits_float_value@plt>
  cmp    $0x1,%eax
  jne    L48
  movss  -0x28(%rbp),%xmm0
  cvtss2sd %xmm0,%xmm0
  movapd %xmm0,%xmm1
  mulsd  -0x58(%rbp),%xmm1
  movss  -0x20(%rbp),%xmm0
  cvtss2sd %xmm0,%xmm0
  mulsd  -0x60(%rbp),%xmm0
  addsd  %xmm1,%xmm0
  pxor   %xmm5,%xmm5
  cvtsd2ss %xmm0,%xmm5
  movd   %xmm5,%edx
  mov    -0x48(%rbp),%rax
  lea    D398(%rip),%rcx
  movd   %edx,%xmm0
  lea    D210(%rip),%rdx
  lea    D410(%rip),%rsi
  mov    %rax,%rdi
  call   <imupdate_fitsf@plt>
  movss  -0x24(%rbp),%xmm0
  cvtss2sd %xmm0,%xmm0
  movapd %xmm0,%xmm1
  mulsd  -0x58(%rbp),%xmm1
  movss  -0x1c(%rbp),%xmm0
  cvtss2sd %xmm0,%xmm0
  mulsd  -0x60(%rbp),%xmm0
  addsd  %xmm1,%xmm0
  pxor   %xmm6,%xmm6
  cvtsd2ss %xmm0,%xmm6
  movd   %xmm6,%edx
  mov    -0x48(%rbp),%rax
  lea    D398(%rip),%rcx
  movd   %edx,%xmm0
  lea    D210(%rip),%rdx
  lea    D56(%rip),%rsi
  mov    %rax,%rdi
  call   <imupdate_fitsf@plt>
  movss  -0x28(%rbp),%xmm0
  movss  D358(%rip),%xmm1
  xorps  %xmm1,%xmm0
  cvtss2sd %xmm0,%xmm0
  movapd %xmm0,%xmm1
  mulsd  -0x60(%rbp),%xmm1
  movss  -0x20(%rbp),%xmm0
  cvtss2sd %xmm0,%xmm0
  mulsd  -0x58(%rbp),%xmm0
  addsd  %xmm1,%xmm0
  pxor   %xmm7,%xmm7
  cvtsd2ss %xmm0,%xmm7
  movd   %xmm7,%edx
  mov    -0x48(%rbp),%rax
  lea    D398(%rip),%rcx
  movd   %edx,%xmm0
  lea    D210(%rip),%rdx
  lea    D54(%rip),%rsi
  mov    %rax,%rdi
  call   <imupdate_fitsf@plt>
  movss  -0x24(%rbp),%xmm0
  movss  D358(%rip),%xmm1
  xorps  %xmm1,%xmm0
  cvtss2sd %xmm0,%xmm0
  movapd %xmm0,%xmm1
  mulsd  -0x60(%rbp),%xmm1
  movss  -0x1c(%rbp),%xmm0
  cvtss2sd %xmm0,%xmm0
  mulsd  -0x58(%rbp),%xmm0
  addsd  %xmm1,%xmm0
  pxor   %xmm3,%xmm3
  cvtsd2ss %xmm0,%xmm3
  movd   %xmm3,%edx
  mov    -0x48(%rbp),%rax
  lea    D398(%rip),%rcx
  movd   %edx,%xmm0
  lea    D210(%rip),%rdx
  lea    D313(%rip),%rsi
  mov    %rax,%rdi
  call   <imupdate_fitsf@plt>
L48:
  movsd  -0x58(%rbp),%xmm0
  mov    -0x60(%rbp),%rdx
  mov    -0x48(%rbp),%rax
  movapd %xmm0,%xmm1
  movq   %rdx,%xmm0
  lea    D175(%rip),%rdx
  mov    %rdx,%rsi
  mov    %rax,%rdi
  call   <imaddhistf@plt>
  mov    $0x1,%eax
  jmp    L31
L50:
  mov    $0x0,%eax
L31:
  mov    -0x18(%rbp),%rdx
  sub    %fs:0x28,%rdx
  je     L53
  call   <__stack_chk_fail@plt>
L53:
  mov    -0x8(%rbp),%rbx
  leave
  ret
```

tool_calls: name: unpack_data
arguments: data_label: D89
datatype: string
name: unpack_data
arguments: data_label: D258
datatype: string
name: unpack_data
arguments: data_label: D398
datatype: string
name: unpack_data
arguments: data_label: D191
datatype: string
name: unpack_data
arguments: data_label: D311
datatype: string
name: unpack_data
arguments: data_label: D202
datatype: string
name: unpack_data
arguments: data_label: D108
datatype: string
name: unpack_data
arguments: data_label: D390
datatype: string
name: unpack_data
arguments: data_label: D210
datatype: string
name: unpack_data
arguments: data_label: D410
datatype: string
name: unpack_data
arguments: data_label: D56
datatype: string
name: unpack_data
arguments: data_label: D54
datatype: string
name: unpack_data
arguments: data_label: D313
datatype: string
name: unpack_data
arguments: data_label: D358
datatype: byte[4]
name: unpack_data
arguments: data_label: D175
datatype: string

tool: D89: CRPIX1

tool: D258: CRPIX2

tool: 

tool: D191: %20.1f / %-47.47s

tool: D311: CD1_1

tool: D202: CD1_2

tool: D108: CD2_1

tool: D390: CD2_2

tool: D210: %20.8f / %-47.47s

tool: D410: PC001001

tool: D56: PC001002

tool: D54: PC002001

tool: D313: PC002002

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

tool: D175: WCS is rotated sin=%f cos=%f

assistant: ```c
int imc_rotate_WCS(struct icom *icom, float dx, float dy, double cos_rot, double sin_rot) {
  float crpix1, crpix2;
  float cd11, cd12, cd21, cd22;
  float pc11, pc12, pc21, pc22;
  if (imget_fits_float_value(icom, "CRPIX1", &crpix1) == 1 && imget_fits_float_value(icom, "CRPIX2", &crpix2) == 1) {
    imupdate_fitsf(icom, "CRPIX1", "%20.1f / %-47.47s", crpix1 * cos_rot - crpix2 * sin_rot + dx, "");
    imupdate_fitsf(icom, "CRPIX2", "%20.1f / %-47.47s", crpix1 * sin_rot + crpix2 * cos_rot + dy, "");
    if ((imget_fits_float_value(icom, "CD1_1", &cd11) + imget_fits_float_value(icom, "CD1_2", &cd12) + imget_fits_float_value(icom, "CD2_1", &cd21) + imget_fits_float_value(icom, "CD2_2", &cd22)) != 0) {
      imupdate_fitsf(icom, "CD1_1", "%20.8f / %-47.47s", (cd11 * cos_rot - cd21 * sin_rot), "");
      imupdate_fitsf(icom, "CD1_2", "%20.8f / %-47.47s", (cd12 * cos_rot - cd22 * sin_rot), "");
      imupdate_fitsf(icom, "CD2_1", "%20.8f / %-47.47s", (cd11 * sin_rot + cd21 * cos_rot), "");
      imupdate_fitsf(icom, "CD2_2", "%20.8f / %-47.47s", (cd12 * sin_rot + cd22 * cos_rot), "");
    }
    if (imget_fits_float_value(icom, "PC001001", &pc11) == 1 && imget_fits_float_value(icom, "PC001002", &pc12) == 1 && imget_fits_float_value(icom, "PC002001", &pc21) == 1 && imget_fits_float_value(icom, "PC002002", &pc22) == 1) {
      imupdate_fitsf(icom, "PC001001", "%20.8f / %-47.47s", (pc11 * cos_rot + pc21 * sin_rot), "");
      imupdate_fitsf(icom, "PC001002", "%20.8f / %-47.47s", (pc12 * cos_rot + pc22 * sin_rot), "");
      imupdate_fitsf(icom, "PC002001", "%20.8f / %-47.47s", (-pc11 * sin_rot + pc21 * cos_rot), "");
      imupdate_fitsf(icom, "PC002002", "%20.8f / %-47.47s", (-pc12 * sin_rot + pc22 * cos_rot), "");
    }
    imaddhistf(icom, "WCS is rotated sin=%f cos=%f", sin_rot, cos_rot);
    return 1;
  }
  return 0;
}
```