user: 请显示与下面列出的汇编代码相匹配的C语言实现：

```asm
<decompose>:
  endbr64
  push   %r15
  movapd %xmm0,%xmm2
  push   %r14
  push   %r13
  push   %r12
  push   %rbp
  push   %rbx
  sub    $0x38,%rsp
  mov    %rdx,0x20(%rsp)
  movsd  %xmm1,0x18(%rsp)
  test   %edi,%edi
  jle    L2
  lea    -0x1(%rdi),%eax
  mov    %edi,%r13d
  mov    %rsi,%rbx
  mov    %rdx,%r15
  mov    %eax,0x2c(%rsp)
  lea    0x1(%rdi),%r12d
  mov    $0x1,%r14d
  jmp    L49
  nopl   (%rax)
L62:
  divsd  %xmm1,%xmm2
  mov    %r15,%rdi
  add    $0x1,%r14
  add    $0x8,%r15
  movsd  %xmm1,0x8(%rsp)
  movapd %xmm2,%xmm0
  call   <modf@plt>
  movsd  0x8(%rsp),%xmm1
  movapd %xmm1,%xmm2
  mulsd  %xmm0,%xmm2
  cmp    %r14,%r12
  je     L59
L49:
  movsd  -0x8(%rbx,%r14,8),%xmm1
  movapd %xmm2,%xmm0
  movsd  %xmm2,0x8(%rsp)
  movsd  %xmm1,0x10(%rsp)
  call   <fmod@plt>
  pxor   %xmm5,%xmm5
  movsd  0x8(%rsp),%xmm2
  movsd  0x10(%rsp),%xmm1
  movapd %xmm0,%xmm3
  addsd  %xmm0,%xmm3
  comisd %xmm3,%xmm5
  jbe    L57
  xorpd  D32(%rip),%xmm3
L57:
  comisd %xmm1,%xmm3
  jbe    L39
  pxor   %xmm5,%xmm5
  movapd %xmm1,%xmm3
  xorpd  D32(%rip),%xmm3
  movapd %xmm5,%xmm4
  cmpltsd %xmm0,%xmm4
  andpd  %xmm4,%xmm3
  andnpd %xmm1,%xmm4
  orpd   %xmm4,%xmm3
  addsd  %xmm3,%xmm0
L39:
  pxor   %xmm6,%xmm6
  comisd %xmm0,%xmm6
  jbe    L24
  xorpd  D32(%rip),%xmm0
L24:
  movsd  0x18(%rsp),%xmm7
  mov    %r15,%rbp
  mov    %r14d,%ecx
  mov    %r14d,%edx
  comisd %xmm0,%xmm7
  jb     L62
  pxor   %xmm6,%xmm6
  movapd %xmm1,%xmm0
  comisd %xmm2,%xmm6
  ja     L13
  mulsd  D54(%rip),%xmm0
L47:
  addsd  %xmm2,%xmm0
  mov    %r15,%rdi
  mov    %edx,0x10(%rsp)
  mov    %ecx,0x8(%rsp)
  divsd  %xmm1,%xmm0
  call   <modf@plt>
  movslq 0x10(%rsp),%rdx
  mov    0x8(%rsp),%ecx
  cmp    %edx,%r13d
  jle    L6
  mov    0x2c(%rsp),%r15d
  mov    0x20(%rsp),%rax
  add    $0x38,%rsp
  xor    %esi,%esi
  pop    %rbx
  pop    %rbp
  sub    %ecx,%r15d
  pop    %r12
  lea    (%rax,%rdx,8),%rdi
  pop    %r13
  lea    0x8(,%r15,8),%r8
  pop    %r14
  pop    %r15
  mov    %r8,%rdx
  jmp    1090 <memset@plt>
L2:
  mov    %rdx,%rbp
  sub    $0x8,%rbp
  nopl   0x0(%rax)
L59:
  addsd  0x0(%rbp),%xmm2
  movsd  %xmm2,0x0(%rbp)
L6:
  add    $0x38,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
  nopl   0x0(%rax)
L13:
  xorpd  D32(%rip),%xmm0
  mulsd  D54(%rip),%xmm0
  jmp    L47
```

tool_calls: name: unpack_data
arguments: tag: D32
type_infered: f32[4]
name: unpack_data
arguments: tag: D54
type_infered: f32[2]

tool: {"D32": [0.0,-0.0,0.0,0.0]}

tool: {"D54": [0.0,1.75]}

assistant: ```c
void decompose(value, uncer, nbasis, basis, count) double value;
double uncer;
int nbasis;
double *basis;
double *count;
{
  int i;
  for (i = 0; i < nbasis; i++) {
    double r = fmod(value, basis[i]);
    if (((2 * r) < 0 ? -(2 * r) : (2 * r)) > basis[i])
      r += r > 0 ? -basis[i] : basis[i];
    if (((r) < 0 ? -(r) : (r)) <= uncer) {
      double half = value < 0 ? -basis[i] / 2 : basis[i] / 2;
      modf((value + half) / basis[i], count + i);
      break;
    }
    value = basis[i] * modf(value / basis[i], count + i);
  }
  if (i >= nbasis) {
    count[--i] += value;
  } else {
    for (i++; i < nbasis; i++)
      count[i] = 0;
  }
}
```