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

```asm
<dJacTx2D>:
  endbr64
  push   %r15
  push   %r14
  push   %r13
  push   %r12
  push   %rbp
  push   %rbx
  sub    $0x148,%rsp
  mov    %rdi,0x58(%rsp)
  mov    %rsi,0x60(%rsp)
  mov    %rdx,0x88(%rsp)
  mov    %fs:0x28,%rax
  mov    %rax,0x138(%rsp)
  xor    %eax,%eax
  cvttsd2si (%rcx),%edx
  cvttsd2si 0x8(%rcx),%eax
  mov    %eax,0x54(%rsp)
  lea    0x1(%rdx),%ecx
  add    $0x1,%eax
  imul   %ecx,%eax
  test   %edx,%edx
  jle    L70
  movslq %ecx,%rbx
  mov    %rbx,0x68(%rsp)
  lea    0x8(,%rbx,8),%r15
  mov    %edx,%esi
  mov    %rsi,0x80(%rsp)
  lea    -0x10(%r15),%rsi
  mov    %rsi,0x70(%rsp)
  cltq
  shl    $0x3,%rax
  mov    %rax,0x8(%rsp)
  movq   $0x0,0x78(%rsp)
  mov    %r15,%r12
  jmp    L60
L68:
  add    $0x1,%r14d
  mov    0x68(%rsp),%rbx
  add    %rbx,0x18(%rsp)
  mov    0x70(%rsp),%rbx
  add    %rbx,0x20(%rsp)
  cmp    %r14d,0x54(%rsp)
  je     L78
L53:
  mov    0x18(%rsp),%rax
  shl    $0x3,%rax
  mov    %rax,0x28(%rsp)
  mov    0x60(%rsp),%rdx
  add    %rax,%rdx
  mov    $0x0,%eax
L85:
  movsd  (%rdx),%xmm0
  movsd  %xmm0,0x90(%rsp,%rax,1)
  movsd  0x8(%rdx),%xmm3
  movsd  %xmm3,0xa0(%rsp,%rax,1)
  movsd  (%rdx,%r12,1),%xmm1
  movsd  %xmm1,0xc0(%rsp,%rax,1)
  movsd  -0x8(%rdx,%r12,1),%xmm2
  movsd  %xmm2,0xb0(%rsp,%rax,1)
  addsd  %xmm3,%xmm0
  addsd  %xmm2,%xmm0
  addsd  %xmm1,%xmm0
  mulsd  D1(%rip),%xmm0
  movsd  %xmm0,0xd0(%rsp,%rax,1)
  movq   $0x0,0xe0(%rsp,%rax,1)
  movq   $0x0,0xf0(%rsp,%rax,1)
  movq   $0x0,0x100(%rsp,%rax,1)
  movq   $0x0,0x110(%rsp,%rax,1)
  movq   $0x0,0x120(%rsp,%rax,1)
  add    0x8(%rsp),%rdx
  add    $0x8,%rax
  cmp    $0x10,%rax
  jne    L85
  mov    0x20(%rsp),%rax
  movsd  (%rax),%xmm0
  movsd  %xmm0,0x10(%rsp)
  lea    0x120(%rsp),%rbx
  lea    0xf0(%rsp),%r13
  lea    0xe0(%rsp),%rax
  lea    0xd0(%rsp),%r15
  lea    0xa0(%rsp),%rbp
  lea    0x90(%rsp),%rdi
  mov    %rbx,%r9
  mov    %r13,%r8
  mov    %rax,0x30(%rsp)
  mov    %rax,%rcx
  mov    %r15,%rdx
  mov    %rbp,%rsi
  mov    %rdi,0x38(%rsp)
  call   <dVTTriangle2D@plt>
  lea    0x110(%rsp),%rax
  lea    0xc0(%rsp),%r10
  mov    %rbx,%r9
  mov    %rax,0x40(%rsp)
  mov    %rax,%r8
  mov    %r13,%rcx
  movsd  0x10(%rsp),%xmm0
  mov    %r15,%rdx
  mov    %r10,0x48(%rsp)
  mov    %r10,%rsi
  mov    %rbp,%rdi
  call   <dVTTriangle2D@plt>
  lea    0x100(%rsp),%r13
  lea    0xb0(%rsp),%rbp
  mov    %rbx,%r9
  mov    %r13,%r8
  mov    0x40(%rsp),%rcx
  movsd  0x10(%rsp),%xmm0
  mov    %r15,%rdx
  mov    %rbp,%rsi
  mov    0x48(%rsp),%rdi
  call   <dVTTriangle2D@plt>
  mov    %rbx,%r9
  mov    0x30(%rsp),%r8
  mov    %r13,%rcx
  movsd  0x10(%rsp),%xmm0
  mov    %r15,%rdx
  mov    0x38(%rsp),%rsi
  mov    %rbp,%rdi
  call   <dVTTriangle2D@plt>
  mov    0x28(%rsp),%rax
  add    0x58(%rsp),%rax
  mov    $0x0,%edx
  lea    0xe0(%rsp),%r8
  lea    0xf0(%rsp),%rdi
  lea    0x110(%rsp),%rsi
  mov    %r13,%rcx
L30:
  movsd  0x120(%rsp,%rdx,1),%xmm0
  mulsd  D1(%rip),%xmm0
  movsd  %xmm0,0x120(%rsp,%rdx,1)
  movapd %xmm0,%xmm1
  addsd  (%rdx,%r8,1),%xmm1
  addsd  (%rax),%xmm1
  movsd  %xmm1,(%rax)
  movapd %xmm0,%xmm1
  addsd  (%rdx,%rdi,1),%xmm1
  addsd  0x8(%rax),%xmm1
  movsd  %xmm1,0x8(%rax)
  movapd %xmm0,%xmm1
  addsd  (%rdx,%rsi,1),%xmm1
  addsd  (%rax,%r12,1),%xmm1
  movsd  %xmm1,(%rax,%r12,1)
  addsd  (%rdx,%rcx,1),%xmm0
  addsd  -0x8(%rax,%r12,1),%xmm0
  movsd  %xmm0,-0x8(%rax,%r12,1)
  add    0x8(%rsp),%rax
  add    $0x8,%rdx
  cmp    $0x10,%rdx
  jne    L30
  jmp    L68
L78:
  addq   $0x1,0x78(%rsp)
  mov    0x78(%rsp),%rax
  cmp    0x80(%rsp),%rax
  je     L70
L60:
  cmpl   $0x0,0x54(%rsp)
  jle    L78
  mov    0x88(%rsp),%rsi
  mov    0x78(%rsp),%rax
  lea    (%rsi,%rax,8),%rcx
  mov    %rcx,0x20(%rsp)
  mov    %rax,0x18(%rsp)
  mov    $0x0,%r13d
  mov    %r13d,%r14d
  jmp    L53
L70:
  mov    0x138(%rsp),%rax
  sub    %fs:0x28,%rax
  jne    L4
  add    $0x148,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
L4:
  call   <__stack_chk_fail@plt>
```

tool_calls: name: parse_data
arguments: data_tag: D1
type_guessed: f64

tool: D1: 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;
      }
    }
  }
}
```