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

```asm
<strip_scale_rgba>:
  endbr64
  test   %edx,%edx
  jle    L77
  push   %r15
  lea    -0x1(%rdx),%eax
  mov    %r8,%r15
  pxor   %xmm6,%xmm6
  push   %r14
  shr    $0x2,%eax
  mov    %rdi,%r14
  push   %r13
  lea    0x1(%rax),%eax
  xor    %r13d,%r13d
  push   %r12
  lea    (%rcx,%rax,4),%rax
  movslq %esi,%r12
  push   %rbp
  mov    %esi,%ebp
  push   %rbx
  mov    %rcx,%rbx
  sub    $0x28,%rsp
  mov    %rax,0x18(%rsp)
  nopl   0x0(%rax)
L38:
  test   %ebp,%ebp
  jle    L65
  lea    0x4(%r13),%r8
  lea    0x8(%r13),%rsi
  movapd %xmm6,%xmm5
  xor    %eax,%eax
  lea    0xc(%r13),%rcx
  movapd %xmm6,%xmm4
  movapd %xmm6,%xmm3
  movapd %xmm6,%xmm2
  xchg   %ax,%ax
L16:
  mov    (%r14,%rax,8),%rdx
  movss  (%r15,%rax,4),%xmm0
  add    $0x1,%rax
  movss  (%rdx,%r13,1),%xmm1
  mulss  %xmm0,%xmm1
  cvtss2sd %xmm1,%xmm1
  addsd  %xmm1,%xmm3
  movss  (%rdx,%r8,1),%xmm1
  mulss  %xmm0,%xmm1
  cvtss2sd %xmm1,%xmm1
  addsd  %xmm1,%xmm2
  movss  (%rdx,%rsi,1),%xmm1
  mulss  %xmm0,%xmm1
  mulss  (%rdx,%rcx,1),%xmm0
  cvtss2sd %xmm1,%xmm1
  addsd  %xmm1,%xmm4
  cvtss2sd %xmm0,%xmm0
  addsd  %xmm0,%xmm5
  cmp    %rax,%r12
  jne    L16
L48:
  movapd %xmm5,%xmm0
  movsd  %xmm4,0x10(%rsp)
  movsd  %xmm2,0x8(%rsp)
  movsd  %xmm3,(%rsp)
  call   <clampf@plt>
  pxor   %xmm6,%xmm6
  movsd  (%rsp),%xmm3
  movsd  0x8(%rsp),%xmm2
  ucomisd %xmm6,%xmm0
  movsd  0x10(%rsp),%xmm4
  movapd %xmm0,%xmm1
  jp     L50
  je     L12
L50:
  divsd  %xmm1,%xmm3
  movsd  %xmm1,(%rsp)
  divsd  %xmm1,%xmm2
  movapd %xmm3,%xmm0
  divsd  %xmm1,%xmm4
  movsd  %xmm2,0x10(%rsp)
  movsd  %xmm4,0x8(%rsp)
  call   <linear_sample_to_srgb@plt>
  movsd  0x10(%rsp),%xmm2
  mov    %al,(%rbx)
  movapd %xmm2,%xmm0
  call   <linear_sample_to_srgb@plt>
  movsd  0x8(%rsp),%xmm4
  mov    %al,0x1(%rbx)
  movapd %xmm4,%xmm0
  call   <linear_sample_to_srgb@plt>
  movsd  D19(%rip),%xmm0
  movsd  (%rsp),%xmm1
  mov    %al,0x2(%rbx)
L53:
  mulsd  %xmm1,%xmm0
  add    $0x4,%rbx
  add    $0x10,%r13
  call   <round@plt>
  pxor   %xmm6,%xmm6
  mov    %al,-0x1(%rbx)
  cmp    0x18(%rsp),%rbx
  jne    L38
  add    $0x28,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
  nop
L12:
  movsd  %xmm0,0x10(%rsp)
  movapd %xmm3,%xmm0
  movsd  %xmm4,0x8(%rsp)
  movsd  %xmm2,(%rsp)
  call   <linear_sample_to_srgb@plt>
  movsd  (%rsp),%xmm2
  mov    %al,(%rbx)
  movapd %xmm2,%xmm0
  call   <linear_sample_to_srgb@plt>
  movsd  0x8(%rsp),%xmm4
  mov    %al,0x1(%rbx)
  movapd %xmm4,%xmm0
  call   <linear_sample_to_srgb@plt>
  movsd  D19(%rip),%xmm0
  movsd  0x10(%rsp),%xmm1
  mov    %al,0x2(%rbx)
  jmp    L53
  nopw   0x0(%rax,%rax,1)
L65:
  movapd %xmm6,%xmm5
  movapd %xmm6,%xmm4
  movapd %xmm6,%xmm3
  movapd %xmm6,%xmm2
  jmp    L48
L77:
  ret
```

tool_calls: name: retrieve_data
arguments: data_tag: D19
data_type: byte[8]

tool: {"D19": [0x00,0x00,0x00,0x00,0x00,0xe0,0x6f,0x40]}

assistant: ```c
void strip_scale_rgba(float **in, int strip_height, int len, unsigned char *out, float *coeffs) {
  int i, j;
  double sum[4], alpha;
  for (i = 0; i < len; i += 4) {
    sum[0] = sum[1] = sum[2] = sum[3] = 0;
    for (j = 0; j < strip_height; j++) {
      sum[0] += coeffs[j] * in[j][i];
      sum[1] += coeffs[j] * in[j][i + 1];
      sum[2] += coeffs[j] * in[j][i + 2];
      sum[3] += coeffs[j] * in[j][i + 3];
    }
    alpha = clampf(sum[3]);
    if (alpha != 0) {
      sum[0] /= alpha;
      sum[1] /= alpha;
      sum[2] /= alpha;
    }
    out[0] = linear_sample_to_srgb(sum[0]);
    out[1] = linear_sample_to_srgb(sum[1]);
    out[2] = linear_sample_to_srgb(sum[2]);
    out[3] = round(alpha * 255.0f);
    out += 4;
  }
}
```