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

```asm
<double_exp_matrix>:
  endbr64
  push   %r15
  push   %r14
  push   %r13
  push   %r12
  push   %rbp
  push   %rbx
  sub    $0xa8,%rsp
  mov    %edi,%ebp
  mov    %ecx,%r14d
  mov    %r8,%rbx
  mov    %fs:0x28,%rax
  mov    %rax,0x98(%rsp)
  xor    %eax,%eax
  test   %ecx,%ecx
  je     L38
  mov    %rdx,%r13
  test   %edi,%edi
  je     L308
  mov    %rdx,%rdi
  lea    -0x1(%rbp),%eax
  shl    $0x5,%rax
  lea    0x20(%rdx,%rax,1),%r9
  jmp    L146
L38:
  test   %edi,%edi
  je     L227
  mov    %r8,%rcx
  movsd  D8(%rip),%xmm0
  jmp    L288
L143:
  movsd  %xmm0,(%rdx,%rax,8)
L121:
  add    $0x1,%rax
  cmp    %eax,%esi
  jbe    L25
L259:
  cmp    %eax,%r14d
  je     L143
  movq   $0x0,(%rdx,%rax,8)
  jmp    L121
L25:
  add    $0x1,%r14d
  add    $0x20,%rcx
  cmp    %r14d,%ebp
  je     L227
L288:
  mov    %rcx,%rdx
  mov    $0x0,%eax
  test   %esi,%esi
  jne    L259
  jmp    L25
L108:
  mov    %rdi,%rcx
  mov    %r8,%rdx
  mov    $0x0,%eax
L72:
  movsd  (%rcx,%rax,8),%xmm0
  movsd  %xmm0,(%rdx,%rax,8)
  add    $0x1,%rax
  cmp    %eax,%esi
  ja     L72
L207:
  add    $0x20,%rdi
  add    $0x20,%r8
  cmp    %r9,%rdi
  je     L308
L146:
  test   %esi,%esi
  jne    L108
  jmp    L207
L308:
  cmp    $0x1,%r14d
  je     L227
  lea    0x10(%rsp),%r10
  lea    -0x1(%rbp),%ecx
  mov    %rcx,%rax
  shl    $0x5,%rax
  lea    0x30(%rsp,%rax,1),%r12
  lea    -0x1(%rsi),%eax
  lea    0x0(,%rax,8),%rdx
  mov    %rbx,%rdi
  lea    0x8(%rbx,%rdx,1),%rbx
  mov    %rbx,0x8(%rsp)
  lea    (%rax,%rcx,4),%rcx
  lea    0x28(%rdi,%rcx,8),%rbx
  lea    0x8(%r13,%rdx,1),%r9
  mov    $0x1,%r15d
  not    %rax
  lea    0x0(,%rax,8),%r11
  mov    %ebp,(%rsp)
  mov    %r10,%rbp
  mov    %r14d,0x4(%rsp)
  mov    %rdi,%r14
  jmp    L21
L97:
  movsd  (%rcx,%rax,8),%xmm0
  movsd  %xmm0,(%rdx,%rax,8)
  add    $0x1,%rax
  cmp    %eax,%esi
  ja     L97
L331:
  add    $0x20,%rdi
  add    $0x20,%rdx
  cmp    %r12,%rdx
  je     L113
L217:
  mov    %rdi,%rcx
  mov    $0x0,%eax
  test   %esi,%esi
  jne    L97
  jmp    L331
L249:
  movq   $0x0,(%rax)
  add    $0x8,%rax
  cmp    %rdx,%rax
  jne    L249
L30:
  add    $0x20,%rdx
  cmp    %rbx,%rdx
  je     L18
L335:
  lea    (%r11,%rdx,1),%rax
  test   %esi,%esi
  jne    L249
  jmp    L30
L113:
  mov    0x8(%rsp),%rdx
  jmp    L335
L221:
  add    $0x8,%r8
  add    $0x8,%rcx
  cmp    %r8,%r9
  je     L12
L203:
  mov    %r8,%rdx
  mov    $0x0,%eax
L267:
  movsd  (%rdi,%rax,8),%xmm0
  mulsd  (%rdx),%xmm0
  addsd  (%rcx),%xmm0
  movsd  %xmm0,(%rcx)
  add    $0x1,%rax
  add    $0x20,%rdx
  cmp    %eax,%esi
  ja     L267
  jmp    L221
L12:
  add    $0x20,%rdi
  add    $0x20,%r10
  cmp    %r12,%rdi
  je     L200
L144:
  mov    %r10,%rcx
  mov    %r13,%r8
  test   %esi,%esi
  jne    L203
  jmp    L12
L18:
  mov    %r14,%r10
  mov    %rbp,%rdi
  jmp    L144
L200:
  add    $0x1,%r15d
  cmp    %r15d,0x4(%rsp)
  jbe    L227
L21:
  mov    %rbp,%rdx
  mov    %r14,%rdi
  cmpl   $0x0,(%rsp)
  jne    L217
  jmp    L200
L227:
  mov    0x98(%rsp),%rax
  sub    %fs:0x28,%rax
  jne    L82
  add    $0xa8,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
L82:
  call   <__stack_chk_fail@plt>
```

tool_calls: name: parse_data
arguments: tag: D8
datatype: f64

tool: D8: 1

assistant: ```c
void double_exp_matrix(unsigned int lines, unsigned int columns, double m1[4][4], unsigned int expNumber, double result[4][4]) {
  unsigned int i, j, k, l;
  double m2[4][4];
  if (expNumber == 0) {
    for (i = 0; i < lines; i++) {
      for (j = 0; j < columns; j++) {
        if (i == j) {
          result[i][j] = 1.0;
        } else {
          result[i][j] = 0.0;
        }
      }
    }
    return;
  }
  for (i = 0; i < lines; i++)
    for (j = 0; j < columns; j++)
      result[i][j] = m1[i][j];
  if (expNumber == 1) {
    return;
  }
  for (l = 1; l < expNumber; l++) {
    for (i = 0; i < lines; i++)
      for (j = 0; j < columns; j++)
        m2[i][j] = result[i][j];
    for (i = 0; i < lines; i++)
      for (j = 0; j < columns; j++)
        result[i][j] = 0;
    for (i = 0; i < lines; i++) {
      for (j = 0; j < columns; j++) {
        for (k = 0; k < columns; k++) {
          result[i][j] = result[i][j] + (m2[i][k] * m1[k][j]);
        }
      }
    }
  }
}
```