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

```asm
<yscale_up_rgba>:
  endbr64
  push   %rbp
  mov    %rsp,%rbp
  push   %rbx
  sub    $0x58,%rsp
  mov    %rdi,-0x48(%rbp)
  mov    %esi,-0x4c(%rbp)
  mov    %rdx,-0x58(%rbp)
  mov    %rcx,-0x60(%rbp)
  mov    %fs:0x28,%rax
  mov    %rax,-0x18(%rbp)
  xor    %eax,%eax
  movl   $0x0,-0x3c(%rbp)
  jmp    L89
L67:
  movl   $0x0,-0x38(%rbp)
  jmp    L47
L80:
  mov    -0x58(%rbp),%rax
  movss  (%rax),%xmm1
  mov    -0x48(%rbp),%rax
  mov    (%rax),%rdx
  mov    -0x3c(%rbp),%ecx
  mov    -0x38(%rbp),%eax
  add    %ecx,%eax
  cltq
  shl    $0x2,%rax
  add    %rdx,%rax
  movss  (%rax),%xmm0
  mulss  %xmm0,%xmm1
  mov    -0x58(%rbp),%rax
  add    $0x4,%rax
  movss  (%rax),%xmm2
  mov    -0x48(%rbp),%rax
  add    $0x8,%rax
  mov    (%rax),%rdx
  mov    -0x3c(%rbp),%ecx
  mov    -0x38(%rbp),%eax
  add    %ecx,%eax
  cltq
  shl    $0x2,%rax
  add    %rdx,%rax
  movss  (%rax),%xmm0
  mulss  %xmm2,%xmm0
  addss  %xmm0,%xmm1
  mov    -0x58(%rbp),%rax
  add    $0x8,%rax
  movss  (%rax),%xmm2
  mov    -0x48(%rbp),%rax
  add    $0x10,%rax
  mov    (%rax),%rdx
  mov    -0x3c(%rbp),%ecx
  mov    -0x38(%rbp),%eax
  add    %ecx,%eax
  cltq
  shl    $0x2,%rax
  add    %rdx,%rax
  movss  (%rax),%xmm0
  mulss  %xmm2,%xmm0
  addss  %xmm0,%xmm1
  mov    -0x58(%rbp),%rax
  add    $0xc,%rax
  movss  (%rax),%xmm2
  mov    -0x48(%rbp),%rax
  add    $0x18,%rax
  mov    (%rax),%rdx
  mov    -0x3c(%rbp),%ecx
  mov    -0x38(%rbp),%eax
  add    %ecx,%eax
  cltq
  shl    $0x2,%rax
  add    %rdx,%rax
  movss  (%rax),%xmm0
  mulss  %xmm2,%xmm0
  addss  %xmm1,%xmm0
  mov    -0x38(%rbp),%eax
  cltq
  movss  %xmm0,-0x30(%rbp,%rax,4)
  addl   $0x1,-0x38(%rbp)
L47:
  cmpl   $0x3,-0x38(%rbp)
  jle    L80
  mov    -0x24(%rbp),%eax
  movd   %eax,%xmm0
  call   <clampf@plt>
  movd   %xmm0,%eax
  mov    %eax,-0x34(%rbp)
  movl   $0x0,-0x38(%rbp)
  jmp    L85
L39:
  pxor   %xmm0,%xmm0
  ucomiss -0x34(%rbp),%xmm0
  jp     L40
  pxor   %xmm0,%xmm0
  ucomiss -0x34(%rbp),%xmm0
  je     L7
L40:
  movss  D15(%rip),%xmm0
  ucomiss -0x34(%rbp),%xmm0
  jp     L60
  movss  D15(%rip),%xmm0
  ucomiss -0x34(%rbp),%xmm0
  je     L7
L60:
  mov    -0x38(%rbp),%eax
  cltq
  movss  -0x30(%rbp,%rax,4),%xmm0
  divss  -0x34(%rbp),%xmm0
  mov    -0x38(%rbp),%eax
  cltq
  movss  %xmm0,-0x30(%rbp,%rax,4)
  mov    -0x38(%rbp),%eax
  cltq
  mov    -0x30(%rbp,%rax,4),%eax
  movd   %eax,%xmm0
  call   <clampf@plt>
  movd   %xmm0,%eax
  mov    -0x38(%rbp),%edx
  movslq %edx,%rdx
  mov    %eax,-0x30(%rbp,%rdx,4)
L7:
  mov    -0x38(%rbp),%eax
  cltq
  mov    -0x30(%rbp,%rax,4),%eax
  mov    -0x3c(%rbp),%ecx
  mov    -0x38(%rbp),%edx
  add    %ecx,%edx
  movslq %edx,%rcx
  mov    -0x60(%rbp),%rdx
  lea    (%rcx,%rdx,1),%rbx
  movd   %eax,%xmm0
  call   <linear_sample_to_srgb@plt>
  mov    %al,(%rbx)
  addl   $0x1,-0x38(%rbp)
L85:
  cmpl   $0x2,-0x38(%rbp)
  jle    L39
  movss  -0x34(%rbp),%xmm1
  movss  D6(%rip),%xmm0
  mulss  %xmm0,%xmm1
  movd   %xmm1,%eax
  mov    -0x3c(%rbp),%edx
  movslq %edx,%rdx
  lea    0x3(%rdx),%rcx
  mov    -0x60(%rbp),%rdx
  lea    (%rcx,%rdx,1),%rbx
  movd   %eax,%xmm0
  call   <f2i@plt>
  mov    %al,(%rbx)
  addl   $0x4,-0x3c(%rbp)
L89:
  mov    -0x3c(%rbp),%eax
  cmp    -0x4c(%rbp),%eax
  jl     L67
  nop
  mov    -0x18(%rbp),%rax
  sub    %fs:0x28,%rax
  je     L25
  call   <__stack_chk_fail@plt>
L25:
  mov    -0x8(%rbp),%rbx
  leave
  ret
```

tool_calls: name: unpack_data
arguments: data_label: D15
type: byte[4]
name: unpack_data
arguments: data_label: D6
type: byte[4]

tool: {"D15": [0x00,0x00,0x80,0x3f]}

tool: {"D6": [0x00,0x00,0x7f,0x43]}

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