user: 请将这段汇编代码片段转换成等效的C代码：

```asm
<dJacTx2D>:
  endbr64
  push   %r15
  push   %r14
  push   %r13
  push   %r12
  push   %rbp
  push   %rbx
  sub    $0x148,%rsp
  mov    %rdx,0x88(%rsp)
  cvttsd2si (%rcx),%edx
  mov    %rdi,0x58(%rsp)
  mov    %rsi,0x60(%rsp)
  mov    %fs:0x28,%rax
  mov    %rax,0x138(%rsp)
  xor    %eax,%eax
  cvttsd2si 0x8(%rcx),%eax
  lea    0x1(%rdx),%ecx
  mov    %eax,0x54(%rsp)
  add    $0x1,%eax
  imul   %ecx,%eax
  test   %edx,%edx
  jle    L23
  movslq %ecx,%rdi
  movslq %edx,%rdx
  cltq
  movsd  D5(%rip),%xmm2
  lea    0x8(,%rdi,8),%r11
  mov    %rdi,0x68(%rsp)
  lea    0x0(,%rax,8),%r10
  mov    %rdx,0x80(%rsp)
  lea    -0x10(%r11),%rdx
  mov    %r11,%rbp
  mov    %rdx,0x70(%rsp)
  movq   $0x0,0x78(%rsp)
  nopw   0x0(%rax,%rax,1)
L50:
  mov    0x54(%rsp),%eax
  test   %eax,%eax
  jle    L17
  mov    0x78(%rsp),%rax
  mov    0x88(%rsp),%rdx
  movl   $0x0,0x50(%rsp)
  lea    0xd0(%rsp),%r13
  lea    0x120(%rsp),%r12
  mov    %r13,%rbx
  lea    0x90(%rsp),%r11
  lea    (%rdx,%rax,8),%rdi
  mov    %rax,0x28(%rsp)
  lea    0xa0(%rsp),%rax
  mov    %r12,%r13
  mov    %rax,0x18(%rsp)
  lea    0xe0(%rsp),%rdx
  lea    0xc0(%rsp),%rax
  mov    %rax,0x10(%rsp)
  lea    0xb0(%rsp),%rax
  lea    0x110(%rsp),%r15
  mov    %rdi,0x20(%rsp)
  lea    0x100(%rsp),%r14
  mov    %rdx,(%rsp)
  mov    %rax,0x8(%rsp)
  lea    0xf0(%rsp),%rax
  mov    %rax,%r12
  nopl   0x0(%rax,%rax,1)
L54:
  mov    0x28(%rsp),%rax
  mov    0x60(%rsp),%rcx
  mov    0x8(%rsp),%rdx
  mov    0x10(%rsp),%rsi
  shl    $0x3,%rax
  mov    0x18(%rsp),%rdi
  mov    %rax,0x30(%rsp)
  add    %rax,%rcx
  xor    %eax,%eax
L42:
  movsd  (%rcx),%xmm4
  movsd  0x8(%rcx),%xmm0
  movsd  -0x8(%rcx,%rbp,1),%xmm3
  movsd  (%rcx,%rbp,1),%xmm1
  add    %r10,%rcx
  movsd  %xmm0,(%rdi,%rax,1)
  addsd  %xmm4,%xmm0
  movsd  %xmm4,(%r11,%rax,1)
  movsd  %xmm1,(%rsi,%rax,1)
  addsd  %xmm3,%xmm0
  movsd  %xmm3,(%rdx,%rax,1)
  addsd  %xmm1,%xmm0
  mulsd  %xmm2,%xmm0
  movsd  %xmm0,(%rbx,%rax,1)
  add    $0x8,%rax
  cmp    $0x10,%rax
  jne    L42
  mov    (%rsp),%rax
  mov    0x20(%rsp),%rsi
  pxor   %xmm0,%xmm0
  mov    %r11,%rdi
  movaps %xmm0,0x0(%r13)
  mov    %r13,%r9
  mov    %r12,%r8
  mov    %rbx,%rdx
  movaps %xmm0,(%r15)
  mov    %rax,%rcx
  movaps %xmm0,(%r14)
  movaps %xmm0,(%r12)
  movaps %xmm0,(%rax)
  movsd  (%rsi),%xmm0
  mov    0x18(%rsp),%rsi
  mov    %r10,0x48(%rsp)
  mov    %r11,0x40(%rsp)
  movsd  %xmm0,0x38(%rsp)
  call   <dVTTriangle2D@plt>
  mov    0x10(%rsp),%rsi
  mov    %r13,%r9
  mov    %r15,%r8
  movsd  0x38(%rsp),%xmm0
  mov    0x18(%rsp),%rdi
  mov    %r12,%rcx
  mov    %rbx,%rdx
  call   <dVTTriangle2D@plt>
  mov    0x8(%rsp),%rsi
  mov    %r13,%r9
  mov    %r14,%r8
  movsd  0x38(%rsp),%xmm0
  mov    0x10(%rsp),%rdi
  mov    %r15,%rcx
  mov    %rbx,%rdx
  call   <dVTTriangle2D@plt>
  mov    0x40(%rsp),%r11
  mov    %r14,%rcx
  mov    %rbx,%rdx
  movsd  0x38(%rsp),%xmm0
  mov    0x8(%rsp),%rdi
  mov    %r13,%r9
  mov    (%rsp),%r8
  mov    %r11,%rsi
  mov    %r11,0x38(%rsp)
  call   <dVTTriangle2D@plt>
  mov    0x48(%rsp),%r10
  mov    (%rsp),%rdx
  xor    %ecx,%ecx
  mov    D5(%rip),%rdi
  mov    0x38(%rsp),%r11
  mov    0x30(%rsp),%rax
  add    0x58(%rsp),%rax
  movq   %rdi,%xmm2
L61:
  movsd  0x0(%r13,%rcx,1),%xmm0
  movsd  (%rdx,%rcx,1),%xmm1
  mulsd  %xmm2,%xmm0
  addsd  %xmm0,%xmm1
  addsd  (%rax),%xmm1
  movsd  %xmm0,0x0(%r13,%rcx,1)
  movsd  %xmm1,(%rax)
  movsd  (%r12,%rcx,1),%xmm1
  addsd  %xmm0,%xmm1
  addsd  0x8(%rax),%xmm1
  movsd  %xmm1,0x8(%rax)
  movsd  (%r15,%rcx,1),%xmm1
  addsd  %xmm0,%xmm1
  addsd  (%r14,%rcx,1),%xmm0
  add    $0x8,%rcx
  addsd  (%rax,%rbp,1),%xmm1
  addsd  -0x8(%rbp,%rax,1),%xmm0
  movsd  %xmm1,(%rax,%rbp,1)
  movsd  %xmm0,-0x8(%rax,%rbp,1)
  add    %r10,%rax
  cmp    $0x10,%rcx
  jne    L61
  addl   $0x1,0x50(%rsp)
  mov    0x68(%rsp),%rdi
  mov    0x70(%rsp),%rdx
  add    %rdi,0x28(%rsp)
  add    %rdx,0x20(%rsp)
  mov    0x50(%rsp),%eax
  cmp    %eax,0x54(%rsp)
  jne    L54
L17:
  addq   $0x1,0x78(%rsp)
  mov    0x78(%rsp),%rax
  cmp    %rax,0x80(%rsp)
  jne    L50
L23:
  mov    0x138(%rsp),%rax
  sub    %fs:0x28,%rax
  jne    L2
  add    $0x148,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
L2:
  call   <__stack_chk_fail@plt>
```

tool_calls: name: disassemble_data
arguments: tag: D5
type_guessed: f64

tool: D5: 0.25

assistant: ```c
void dJacTx2D(double *prod, const double *yc, const double *x, const double *m, int n) {
  double yA[2], yB[2], yC[2], yD[2], yM[2], vol;
  double rA[2], rB[2], rC[2], rD[2], rM[2];
  double xp;
  int i1, i2, r, id, idn, comp;
  const int m1 = m[0], m2 = m[1];
  const int nn = (m1 + 1) * (m2 + 1);
  const int inc_x = 1;
  const int inc_y = m1 + 1;
  const double *ptr;
  double *wptr;
#pragma omp parallel for default(shared) private(ptr, wptr, i1, i2, r, id, idn, comp, xp, yA, yB, yC, yD, yM, rA, rB, rC, rD, rM)
  for (i1 = 0; i1 < m1; i1++) {
    for (i2 = 0; i2 < m2; i2++) {
      idn = i1 + i2 * inc_y;
      id = i1 + i2 * m1;
      comp = 0;
      for (r = 0; r < 2; r++) {
        ptr = yc + idn + comp;
        yA[r] = *ptr;
        ptr = ptr + inc_x;
        yB[r] = *ptr;
        ptr = ptr + inc_y;
        yD[r] = *ptr;
        ptr = ptr - inc_x;
        yC[r] = *ptr;
        yM[r] = 0.25 * (yA[r] + yB[r] + yC[r] + yD[r]);
        rA[r] = 0.0;
        rB[r] = 0.0;
        rC[r] = 0.0;
        rD[r] = 0.0;
        rM[r] = 0.0;
        comp += nn;
      }
      xp = x[id];
      dVTTriangle2D(yA, yB, yM, xp, rA, rB, rM);
      dVTTriangle2D(yB, yD, yM, xp, rB, rD, rM);
      dVTTriangle2D(yD, yC, yM, xp, rD, rC, rM);
      dVTTriangle2D(yC, yA, yM, xp, rC, rA, rM);
      comp = 0;
      for (r = 0; r < 2; r++) {
        wptr = prod + idn + comp;
        rM[r] *= .25;
#pragma omp atomic
        *wptr += rA[r] + rM[r];
        wptr += inc_x;
#pragma omp atomic
        *wptr += rB[r] + rM[r];
        wptr += inc_y;
#pragma omp atomic
        *wptr += rD[r] + rM[r];
        wptr -= inc_x;
#pragma omp atomic
        *wptr += rC[r] + rM[r];
        comp += nn;
      }
    }
  }
}
```