user: 请从给定的汇编指令推导出C源代码：

```asm
<colorMap_rainbow>:
  endbr64
  push   %r15
  push   %r14
  push   %r13
  push   %r12
  mov    $0x3,%r12d
  push   %rbp
  push   %rbx
  mov    %rdi,%rbx
  sub    $0x18,%rsp
  mulsd  D87(%rip),%xmm0
  addsd  D76(%rip),%xmm0
  movapd %xmm0,%xmm5
  subsd  D59(%rip),%xmm5
  movq   %xmm0,%rbp
  movapd %xmm5,%xmm0
  movq   %xmm5,%r15
  call   <fabs@plt>
  movq   %rbp,%xmm1
  subsd  D88(%rip),%xmm1
  mov    %eax,%r13d
  movapd %xmm1,%xmm0
  movsd  %xmm1,0x8(%rsp)
  call   <fabs@plt>
  pxor   %xmm0,%xmm0
  mov    %eax,%r8d
  mov    %r12d,%eax
  sub    %r13d,%eax
  sub    %r8d,%eax
  cmp    $0xffffffff,%eax
  jl     L7
  movq   %r15,%xmm0
  call   <fabs@plt>
  movsd  0x8(%rsp),%xmm1
  mov    %eax,%r13d
  movapd %xmm1,%xmm0
  sub    %r13d,%r12d
  call   <fabs@plt>
  pxor   %xmm0,%xmm0
  sub    %eax,%r12d
  mov    %r12d,%eax
  shr    $0x1f,%eax
  add    %r12d,%eax
  sar    %eax
  cvtsi2sd %eax,%xmm0
L7:
  movq   %rbp,%xmm6
  movsd  %xmm0,(%rbx)
  subsd  D39(%rip),%xmm6
  mov    $0x4,%r12d
  movapd %xmm6,%xmm0
  movq   %xmm6,%r14
  call   <fabs@plt>
  movq   %r15,%xmm0
  mov    %eax,%r13d
  call   <fabs@plt>
  pxor   %xmm0,%xmm0
  mov    %eax,%r8d
  mov    %r12d,%eax
  sub    %r13d,%eax
  sub    %r8d,%eax
  cmp    $0xffffffff,%eax
  jl     L19
  movq   %r14,%xmm0
  call   <fabs@plt>
  movq   %r15,%xmm0
  mov    %eax,%r13d
  call   <fabs@plt>
  sub    %r13d,%r12d
  pxor   %xmm0,%xmm0
  sub    %eax,%r12d
  mov    %r12d,%eax
  shr    $0x1f,%eax
  add    %r12d,%eax
  sar    %eax
  cvtsi2sd %eax,%xmm0
L19:
  movq   %rbp,%xmm1
  movsd  %xmm0,0x8(%rbx)
  mov    $0x3,%ebp
  subsd  D19(%rip),%xmm1
  movapd %xmm1,%xmm0
  movsd  %xmm1,0x8(%rsp)
  call   <fabs@plt>
  movq   %r14,%xmm0
  mov    %eax,%r12d
  call   <fabs@plt>
  mov    %eax,%r8d
  mov    %ebp,%eax
  sub    %r12d,%eax
  sub    %r8d,%eax
  cmp    $0xffffffff,%eax
  jge    L21
  movupd (%rbx),%xmm4
  pxor   %xmm1,%xmm1
  add    $0x10,%rbx
  movapd D19(%rip),%xmm3
  movq   $0x0,(%rbx)
  movapd %xmm4,%xmm0
  cmpltpd %xmm4,%xmm1
  cmpltpd %xmm3,%xmm0
  movdqa %xmm0,%xmm2
  pand   %xmm1,%xmm0
  andpd  %xmm0,%xmm4
  andnpd %xmm3,%xmm0
  pandn  %xmm1,%xmm2
  orpd   %xmm4,%xmm0
  andpd  %xmm2,%xmm3
  andnpd %xmm0,%xmm2
  pxor   %xmm0,%xmm0
  orpd   %xmm3,%xmm2
  andpd  %xmm2,%xmm1
  movups %xmm1,-0x10(%rbx)
L10:
  movsd  %xmm0,(%rbx)
  add    $0x18,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
  nopl   0x0(%rax)
L21:
  movsd  0x8(%rsp),%xmm1
  add    $0x10,%rbx
  movapd %xmm1,%xmm0
  call   <fabs@plt>
  movq   %r14,%xmm0
  mov    %eax,%r12d
  call   <fabs@plt>
  pxor   %xmm0,%xmm0
  movupd -0x10(%rbx),%xmm4
  movapd D19(%rip),%xmm3
  mov    %eax,%r8d
  mov    %ebp,%eax
  pxor   %xmm1,%xmm1
  mov    D19(%rip),%rcx
  sub    %r12d,%eax
  cmpltpd %xmm4,%xmm1
  sub    %r8d,%eax
  mov    %eax,%edx
  shr    $0x1f,%edx
  add    %eax,%edx
  sar    %edx
  cvtsi2sd %edx,%xmm0
  movsd  %xmm0,(%rbx)
  movapd %xmm4,%xmm0
  cmpltpd %xmm3,%xmm0
  movdqa %xmm0,%xmm2
  pand   %xmm1,%xmm0
  andpd  %xmm0,%xmm4
  andnpd %xmm3,%xmm0
  pandn  %xmm1,%xmm2
  orpd   %xmm4,%xmm0
  andpd  %xmm2,%xmm3
  andnpd %xmm0,%xmm2
  movq   %rcx,%xmm0
  orpd   %xmm3,%xmm2
  andpd  %xmm2,%xmm1
  movups %xmm1,-0x10(%rbx)
  cmp    $0x1,%eax
  jg     L10
  pxor   %xmm0,%xmm0
  jmp    L10
```

tool_calls: name: extract_symbolic_data
arguments: data_label: D87
type_guessed: byte[8]
name: extract_symbolic_data
arguments: data_label: D76
type_guessed: f64
name: extract_symbolic_data
arguments: data_label: D59
type_guessed: f32[2]
name: extract_symbolic_data
arguments: data_label: D88
type_guessed: f32[2]
name: extract_symbolic_data
arguments: data_label: D39
type_guessed: f64
name: extract_symbolic_data
arguments: data_label: D19
type_guessed: f64[2]

tool: {"D87": [0x33,0x33,0x33,0x33,0x33,0x33,0x17,0x40]}

tool: D76: 0.1

tool: {"D59": [0.0,2.25]}

tool: {"D88": [0.0,2.3125]}

tool: D39: 2

tool: {"D19": [1.0,1.0]}

assistant: ```c
void colorMap_rainbow(GLfloat fin, GLfloat *fcol) {
  int i;
  GLfloat xmin, xmax;
  GLfloat dx = 0.10;
  GLfloat ga = (6.0 - 2.0 * dx) * fin + dx;
  fcol[0] = ((0) > ((3 - fabs(ga - 4) - fabs(ga - 5)) / 2) ? (0) : ((3 - fabs(ga - 4) - fabs(ga - 5)) / 2));
  fcol[1] = ((0) > ((4 - fabs(ga - 2) - fabs(ga - 4)) / 2) ? (0) : ((4 - fabs(ga - 2) - fabs(ga - 4)) / 2));
  fcol[2] = ((0) > ((3 - fabs(ga - 1) - fabs(ga - 2)) / 2) ? (0) : ((3 - fabs(ga - 1) - fabs(ga - 2)) / 2));
  for (i = 0; i < 3; i++) {
    xmin = ((fcol[i]) > (0.0) ? (fcol[i]) : (0.0));
    fcol[i] = ((xmin) < (1.0) ? (xmin) : (1.0));
  }
}
```