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

```asm
<compute_flow_contrib_rm>:
  endbr64
  push   %r15
  mov    %rdx,%rax
  push   %r14
  push   %r13
  movslq %edi,%r13
  push   %r12
  lea    0x0(,%r13,4),%r12
  push   %rbp
  lea    (%rsi,%r12,1),%r15
  push   %rbx
  mov    %r9,%rbx
  sub    $0x88,%rsp
  mov    (%r15),%edi
  mov    %rcx,(%rsp)
  mov    0xc0(%rsp),%r9
  mov    0xc8(%rsp),%r14
  mov    %fs:0x28,%rdx
  mov    %rdx,0x78(%rsp)
  xor    %edx,%edx
  movl   $0x0,0x4c(%rsp)
  movq   $0x0,0x50(%rsp)
  movq   $0x0,0x58(%rsp)
  movq   $0x0,0x60(%rsp)
  movq   $0x0,0x68(%rsp)
  test   %edi,%edi
  jg     L75
L68:
  pxor   %xmm0,%xmm0
  movapd %xmm0,%xmm3
  movapd %xmm0,%xmm1
  movapd %xmm0,%xmm2
L43:
  mov    (%rsp),%rax
  movsd  %xmm0,(%rax)
  movsd  %xmm3,0x8(%rax)
  movsd  %xmm1,0x10(%rax)
  movsd  %xmm2,0x18(%rax)
  mov    0x78(%rsp),%rax
  sub    %fs:0x28,%rax
  jne    L62
  add    $0x88,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
  xchg   %ax,%ax
L75:
  mov    %rsi,%r10
  addl   $0x1,(%rax,%r12,1)
  lea    0x4c(%rsp),%rsi
  mov    %r8,%rbp
  lea    0x50(%rsp),%rdi
  sub    $0x8,%rsp
  mov    %rbp,%rcx
  mov    %r10,%rdx
  push   %r14
  mov    %rbx,%r8
  mov    %r9,0x18(%rsp)
  call   <compute_evt_plane_avg@plt>
  mov    (%r15),%ecx
  pop    %rax
  pop    %rdx
  test   %ecx,%ecx
  jle    L68
  mov    0x8(%rsp),%r9
  shl    $0x3,%r13
  add    %r14,%r12
  movq   $0x0,0x10(%rsp)
  add    %r13,%rbp
  add    %r13,%rbx
  movq   $0x0,0x8(%rsp)
  lea    (%r9,%r13,1),%r14
  movl   $0x0,0x1c(%rsp)
  xor    %r13d,%r13d
  jmp    L19
L24:
  add    $0x1,%r13d
  add    $0x64,%r12
  add    $0xc8,%r14
  add    $0xc8,%rbp
  add    $0xc8,%rbx
  cmp    %r13d,(%r15)
  jle    L10
L19:
  pxor   %xmm0,%xmm0
  cvtsi2sdl (%r12),%xmm0
  call   <abs@plt>
  pxor   %xmm5,%xmm5
  comisd %xmm5,%xmm0
  jbe    L24
  movsd  (%r14),%xmm0
  call   <abs@plt>
  movsd  D49(%rip),%xmm6
  comisd %xmm0,%xmm6
  jb     L24
  movsd  (%rbx),%xmm0
  movsd  0x0(%rbp),%xmm7
  addl   $0x1,0x1c(%rsp)
  addsd  %xmm0,%xmm0
  movsd  %xmm7,0x20(%rsp)
  cvttsd2si %xmm0,%edi
  call   <cos@plt>
  movsd  %xmm0,0x30(%rsp)
  movsd  (%rbx),%xmm0
  addsd  %xmm0,%xmm0
  call   <sin@plt>
  movsd  0x0(%rbp),%xmm3
  movsd  0x30(%rsp),%xmm1
  mulsd  0x20(%rsp),%xmm1
  mov    0x4c(%rsp),%edi
  mulsd  %xmm0,%xmm3
  movsd  0x50(%rsp),%xmm0
  movsd  %xmm3,0x28(%rsp)
  call   <subMean@plt>
  movsd  0x28(%rsp),%xmm3
  mov    0x4c(%rsp),%edi
  movsd  %xmm0,0x20(%rsp)
  movsd  0x58(%rsp),%xmm0
  movapd %xmm3,%xmm1
  call   <subMean@plt>
  movsd  0x20(%rsp),%xmm2
  movapd %xmm2,%xmm1
  call   <atan2@plt>
  movsd  0x0(%rbp),%xmm7
  movsd  %xmm0,0x38(%rsp)
  movsd  D58(%rip),%xmm0
  mulsd  (%rbx),%xmm0
  movsd  %xmm7,0x20(%rsp)
  cvttsd2si %xmm0,%edi
  call   <cos@plt>
  movsd  %xmm0,0x30(%rsp)
  movsd  D58(%rip),%xmm0
  mulsd  (%rbx),%xmm0
  call   <sin@plt>
  movsd  0x0(%rbp),%xmm4
  movsd  0x30(%rsp),%xmm1
  mulsd  0x20(%rsp),%xmm1
  mov    0x4c(%rsp),%edi
  mulsd  %xmm0,%xmm4
  movsd  0x60(%rsp),%xmm0
  movsd  %xmm4,0x28(%rsp)
  call   <subMean@plt>
  movsd  0x28(%rsp),%xmm4
  mov    0x4c(%rsp),%edi
  movsd  %xmm0,0x20(%rsp)
  movsd  0x68(%rsp),%xmm0
  movapd %xmm4,%xmm1
  call   <subMean@plt>
  movsd  0x20(%rsp),%xmm2
  movapd %xmm2,%xmm1
  call   <atan2@plt>
  movsd  0x38(%rsp),%xmm3
  movsd  (%rbx),%xmm1
  movsd  %xmm0,0x20(%rsp)
  movsd  D11(%rip),%xmm0
  addsd  %xmm1,%xmm1
  mulsd  %xmm3,%xmm0
  subsd  %xmm0,%xmm1
  cvttsd2si %xmm1,%edi
  call   <cos@plt>
  movsd  D58(%rip),%xmm1
  mulsd  (%rbx),%xmm1
  movsd  0x20(%rsp),%xmm2
  divsd  D58(%rip),%xmm2
  addsd  0x8(%rsp),%xmm0
  movsd  %xmm0,0x8(%rsp)
  subsd  %xmm2,%xmm1
  cvttsd2si %xmm1,%edi
  call   <cos@plt>
  addsd  0x10(%rsp),%xmm0
  movsd  %xmm0,0x10(%rsp)
  jmp    L24
  nopl   0x0(%rax)
L10:
  mov    0x1c(%rsp),%eax
  test   %eax,%eax
  je     L68
  pxor   %xmm2,%xmm2
  movsd  0x8(%rsp),%xmm0
  movsd  0x10(%rsp),%xmm1
  cvtsi2sd %eax,%xmm2
  divsd  %xmm2,%xmm0
  divsd  %xmm2,%xmm1
  movapd %xmm0,%xmm3
  mulsd  %xmm0,%xmm3
  movapd %xmm1,%xmm2
  mulsd  %xmm1,%xmm2
  jmp    L43
L62:
  call   <__stack_chk_fail@plt>
```

tool_calls: name: disassemble_data
arguments: data_tag: D49
type_guessed: f64
name: disassemble_data
arguments: data_tag: D58
type_guessed: f64
name: disassemble_data
arguments: data_tag: D11
type_guessed: f32[2]

tool: D49: 1

tool: D58: 3

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

assistant: ```c
void compute_flow_contrib_rm(int ibin, int *ptCount, int *evCount, double *EventBinFlowContrib, double ptArray[10000][25], double phiArray[10000][25], double rapArray[10000][25], int chArray[10000][25]) {
  int nch = 0;
  int i;
  int nPartsEvtPlane = 0;
  double rapCutFlow = 1.0;
  double evtPlanesAvg[4] = {0.0, 0.0, 0.0, 0.0};
  double evtPlanes[2] = {0.0, 0.0};
  double meanCosVn = 0.0;
  double meanSinVn = 0.0;
  double cosVn = 0.0;
  double sinVn = 0.0;
  double v2Sum = 0.0;
  double v3Sum = 0.0;
  if (ptCount[ibin] > 0) {
    evCount[ibin]++;
    compute_evt_plane_avg(evtPlanesAvg, &nPartsEvtPlane, ptCount, ptArray, phiArray, rapArray, chArray);
    for (i = 0; i < ptCount[ibin]; i++) {
      if (abs(chArray[i][ibin]) > 0 && abs(rapArray[i][ibin]) <= rapCutFlow) {
        nch++;
        cosVn = ptArray[i][ibin] * cos(2.0 * phiArray[i][ibin]);
        sinVn = ptArray[i][ibin] * sin(2.0 * phiArray[i][ibin]);
        meanCosVn = subMean(nPartsEvtPlane, evtPlanesAvg[0], cosVn);
        meanSinVn = subMean(nPartsEvtPlane, evtPlanesAvg[1], sinVn);
        evtPlanes[0] = atan2(meanSinVn, meanCosVn) / 2.0;
        cosVn = ptArray[i][ibin] * cos(3.0 * phiArray[i][ibin]);
        sinVn = ptArray[i][ibin] * sin(3.0 * phiArray[i][ibin]);
        meanCosVn = subMean(nPartsEvtPlane, evtPlanesAvg[2], cosVn);
        meanSinVn = subMean(nPartsEvtPlane, evtPlanesAvg[3], sinVn);
        evtPlanes[1] = atan2(meanSinVn, meanCosVn) / 3.0;
        v2Sum += cos(2 * phiArray[i][ibin] - evtPlanes[0]);
        v3Sum += cos(3 * phiArray[i][ibin] - evtPlanes[1]);
      }
    }
  }
  if (nch > 0) {
    EventBinFlowContrib[0] = v2Sum / (double)nch;
    EventBinFlowContrib[1] = pow(v2Sum / (double)nch, 2.0);
    EventBinFlowContrib[2] = v3Sum / (double)nch;
    EventBinFlowContrib[3] = pow(v3Sum / (double)nch, 2.0);
  } else {
    EventBinFlowContrib[0] = 0.0;
    EventBinFlowContrib[1] = 0.0;
    EventBinFlowContrib[2] = 0.0;
    EventBinFlowContrib[3] = 0.0;
  }
}
```