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

```asm
<dmat_scale_add_2x2>:
  endbr64
  push   %r15
  push   %r14
  push   %r13
  push   %r12
  push   %rbp
  push   %rbx
  mov    %rdi,%rbx
  mov    %rdx,%rax
  mov    %rcx,%rdx
  ucomisd D7(%rip),%xmm0
  jp     L60
  mov    %rdi,%r14
  jne    L60
  cmpl   $0x0,0x68(%rsp)
  jle    L108
  mov    %rdi,%r12
  mov    0x48(%rsp),%rbp
  mov    0x50(%rsp),%rbx
  mov    %rax,%rcx
  mov    0x38(%rsp),%r11
  mov    0x58(%rsp),%r10
  mov    0x40(%rsp),%rdi
  mov    0x60(%rsp),%r13
  mov    0x68(%rsp),%eax
  lea    -0x1(%rax),%eax
  lea    0x8(%r14,%rax,8),%r15
  mov    0x68(%rsp),%eax
  lea    0x0(,%rax,8),%r14
L8:
  mov    $0x0,%eax
  mov    %r15,-0x8(%rsp)
L105:
  mov    (%r12),%r15
  movsd  (%r15,%rax,1),%xmm0
  mov    (%r8),%r15
  addsd  (%r15,%rax,1),%xmm0
  mov    0x0(%rbp),%r15
  movsd  %xmm0,(%r15,%rax,1)
  mov    (%rsi),%r15
  movsd  (%r15,%rax,1),%xmm0
  mov    (%r9),%r15
  addsd  (%r15,%rax,1),%xmm0
  mov    (%rbx),%r15
  movsd  %xmm0,(%r15,%rax,1)
  mov    (%rcx),%r15
  movsd  (%r15,%rax,1),%xmm0
  mov    (%r11),%r15
  addsd  (%r15,%rax,1),%xmm0
  mov    (%r10),%r15
  movsd  %xmm0,(%r15,%rax,1)
  mov    (%rdx),%r15
  movsd  (%r15,%rax,1),%xmm0
  mov    (%rdi),%r15
  addsd  (%r15,%rax,1),%xmm0
  mov    0x0(%r13),%r15
  movsd  %xmm0,(%r15,%rax,1)
  add    $0x8,%rax
  cmp    %r14,%rax
  jne    L105
  mov    -0x8(%rsp),%r15
  add    $0x8,%r12
  add    $0x8,%r8
  add    $0x8,%rbp
  add    $0x8,%rsi
  add    $0x8,%r9
  add    $0x8,%rbx
  add    $0x8,%rcx
  add    $0x8,%r11
  add    $0x8,%r10
  add    $0x8,%rdx
  add    $0x8,%rdi
  add    $0x8,%r13
  cmp    %r15,%r12
  jne    L8
L108:
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
L60:
  ucomisd D4(%rip),%xmm0
  jp     L123
  jne    L123
  cmpl   $0x0,0x68(%rsp)
  jle    L108
  mov    %rbx,%r13
  mov    %rbx,%r14
  mov    0x48(%rsp),%rbp
  mov    0x50(%rsp),%rbx
  mov    %rax,%rcx
  mov    0x38(%rsp),%r11
  mov    0x58(%rsp),%r10
  mov    0x40(%rsp),%rdi
  mov    0x60(%rsp),%r15
  mov    0x68(%rsp),%eax
  lea    -0x1(%rax),%eax
  lea    0x8(%r13,%rax,8),%r13
  mov    0x68(%rsp),%eax
  lea    0x0(,%rax,8),%r12
  mov    %r8,%rax
  mov    %r12,%r8
  mov    %rax,%r12
L24:
  mov    $0x0,%eax
  mov    %r13,-0x8(%rsp)
  mov    %r14,%r13
L104:
  mov    0x0(%r13),%r14
  movsd  (%r14,%rax,1),%xmm0
  mov    (%r12),%r14
  subsd  (%r14,%rax,1),%xmm0
  mov    0x0(%rbp),%r14
  movsd  %xmm0,(%r14,%rax,1)
  mov    (%rsi),%r14
  movsd  (%r14,%rax,1),%xmm0
  mov    (%r9),%r14
  subsd  (%r14,%rax,1),%xmm0
  mov    (%rbx),%r14
  movsd  %xmm0,(%r14,%rax,1)
  mov    (%rcx),%r14
  movsd  (%r14,%rax,1),%xmm0
  mov    (%r11),%r14
  subsd  (%r14,%rax,1),%xmm0
  mov    (%r10),%r14
  movsd  %xmm0,(%r14,%rax,1)
  mov    (%rdx),%r14
  movsd  (%r14,%rax,1),%xmm0
  mov    (%rdi),%r14
  subsd  (%r14,%rax,1),%xmm0
  mov    (%r15),%r14
  movsd  %xmm0,(%r14,%rax,1)
  add    $0x8,%rax
  cmp    %r8,%rax
  jne    L104
  mov    -0x8(%rsp),%rax
  mov    %r13,%r14
  mov    %rax,%r13
  add    $0x8,%r14
  add    $0x8,%r12
  add    $0x8,%rbp
  add    $0x8,%rsi
  add    $0x8,%r9
  add    $0x8,%rbx
  add    $0x8,%rcx
  add    $0x8,%r11
  add    $0x8,%r10
  add    $0x8,%rdx
  add    $0x8,%rdi
  add    $0x8,%r15
  cmp    %rax,%r14
  jne    L24
  jmp    L108
L123:
  cmpl   $0x0,0x68(%rsp)
  jle    L108
  mov    %rbx,%rdi
  mov    %rbx,%r10
  mov    0x48(%rsp),%r14
  mov    0x50(%rsp),%r13
  mov    %rax,%rcx
  mov    0x38(%rsp),%r12
  mov    0x58(%rsp),%rbp
  mov    0x40(%rsp),%rbx
  mov    0x60(%rsp),%r11
  mov    0x68(%rsp),%eax
  lea    -0x1(%rax),%eax
  lea    0x8(%rdi,%rax,8),%rax
  mov    0x68(%rsp),%r15d
  shl    $0x3,%r15
  mov    %rax,%rdi
L59:
  mov    $0x0,%eax
  mov    %rdi,-0x8(%rsp)
L97:
  mov    (%r8),%rdi
  movapd %xmm0,%xmm1
  mulsd  (%rdi,%rax,1),%xmm1
  mov    (%r10),%rdi
  addsd  (%rdi,%rax,1),%xmm1
  mov    (%r14),%rdi
  movsd  %xmm1,(%rdi,%rax,1)
  mov    (%r9),%rdi
  movapd %xmm0,%xmm1
  mulsd  (%rdi,%rax,1),%xmm1
  mov    (%rsi),%rdi
  addsd  (%rdi,%rax,1),%xmm1
  mov    0x0(%r13),%rdi
  movsd  %xmm1,(%rdi,%rax,1)
  mov    (%r12),%rdi
  movapd %xmm0,%xmm1
  mulsd  (%rdi,%rax,1),%xmm1
  mov    (%rcx),%rdi
  addsd  (%rdi,%rax,1),%xmm1
  mov    0x0(%rbp),%rdi
  movsd  %xmm1,(%rdi,%rax,1)
  mov    (%rbx),%rdi
  movapd %xmm0,%xmm1
  mulsd  (%rdi,%rax,1),%xmm1
  mov    (%rdx),%rdi
  addsd  (%rdi,%rax,1),%xmm1
  mov    (%r11),%rdi
  movsd  %xmm1,(%rdi,%rax,1)
  add    $0x8,%rax
  cmp    %r15,%rax
  jne    L97
  mov    -0x8(%rsp),%rdi
  add    $0x8,%r10
  add    $0x8,%r8
  add    $0x8,%r14
  add    $0x8,%rsi
  add    $0x8,%r9
  add    $0x8,%r13
  add    $0x8,%rcx
  add    $0x8,%r12
  add    $0x8,%rbp
  add    $0x8,%rdx
  add    $0x8,%rbx
  add    $0x8,%r11
  cmp    %rdi,%r10
  jne    L59
  jmp    L108
```

tool_calls: name: unpack_data
arguments: tag: D7
type_infered: f64
name: unpack_data
arguments: tag: D4
type_infered: byte[8]

tool: D7: 1

tool: {"D4": [0x00,0x00,0x00,0x00,0x00,0x00,0xf0,0xbf]}

assistant: ```c
void dmat_scale_add_2x2(double **a11, double **a12, double **a21, double **a22, double **b11, double **b12, double **b21, double **b22, double **P11, double **P12, double **P21, double **P22, double alpha, int n) {
  int i;
  int j;
  if (alpha == 1.) {
    for (i = 0; i < n; ++i) {
      for (j = 0; j < n; ++j) {
        P11[i][j] = a11[i][j] + b11[i][j];
        P12[i][j] = a12[i][j] + b12[i][j];
        P21[i][j] = a21[i][j] + b21[i][j];
        P22[i][j] = a22[i][j] + b22[i][j];
      }
    }
  } else if (alpha == -1.) {
    for (i = 0; i < n; ++i) {
      for (j = 0; j < n; ++j) {
        P11[i][j] = a11[i][j] - b11[i][j];
        P12[i][j] = a12[i][j] - b12[i][j];
        P21[i][j] = a21[i][j] - b21[i][j];
        P22[i][j] = a22[i][j] - b22[i][j];
      }
    }
  } else {
    for (i = 0; i < n; ++i) {
      for (j = 0; j < n; ++j) {
        P11[i][j] = a11[i][j] + b11[i][j] * alpha;
        P12[i][j] = a12[i][j] + b12[i][j] * alpha;
        P21[i][j] = a21[i][j] + b21[i][j] * alpha;
        P22[i][j] = a22[i][j] + b22[i][j] * alpha;
      }
    }
  }
}
```