user: What is the c source code of the assembly code below?

```asm
<DerivePrincipalStresses>:
  endbr64
  push   %r15
  mov    $0x3,%ecx
  mov    %rsi,%r15
  mov    $0x3,%esi
  push   %r14
  push   %r13
  mov    %rdi,%r13
  mov    $0x1,%edi
  push   %r12
  push   %rbp
  push   %rbx
  mov    %rdx,%rbx
  mov    $0x1,%edx
  sub    $0x8,%rsp
  call   <dmatrix@plt>
  mov    $0x3,%ecx
  mov    $0x1,%edx
  mov    $0x3,%esi
  mov    $0x1,%edi
  mov    %rax,%r14
  call   <dmatrix@plt>
  mov    $0x3,%esi
  mov    $0x1,%edi
  mov    %rax,%r12
  call   <dvector@plt>
  mov    0x10(%r14),%rdx
  mov    0x18(%r14),%rcx
  mov    %r14,%rdi
  mov    %rax,%rbp
  movsd  0x0(%r13),%xmm0
  mov    0x8(%r14),%rax
  mov    %rbp,%rsi
  movsd  %xmm0,0x8(%rax)
  movsd  0x20(%r13),%xmm0
  movsd  %xmm0,0x10(%rdx)
  movsd  0x40(%r13),%xmm0
  movsd  %xmm0,0x18(%rcx)
  movsd  0x8(%r13),%xmm0
  movsd  %xmm0,0x10(%rax)
  movsd  0x10(%r13),%xmm0
  movsd  %xmm0,0x18(%rax)
  movsd  0x28(%r13),%xmm0
  movsd  %xmm0,0x18(%rdx)
  movsd  0x10(%rax),%xmm0
  movsd  %xmm0,0x8(%rdx)
  movsd  0x18(%rax),%xmm0
  movsd  %xmm0,0x8(%rcx)
  movsd  0x18(%rdx),%xmm0
  mov    %r12,%rdx
  movsd  %xmm0,0x10(%rcx)
  call   <jacobi@plt>
  movsd  0x8(%rbp),%xmm1
  mov    %r12,%rsi
  mov    %rbp,%rdi
  movq   D11(%rip),%xmm0
  xorpd  %xmm0,%xmm1
  movsd  %xmm1,0x8(%rbp)
  movsd  0x10(%rbp),%xmm1
  xorpd  %xmm0,%xmm1
  movsd  %xmm1,0x10(%rbp)
  movsd  0x18(%rbp),%xmm1
  xorpd  %xmm0,%xmm1
  movsd  %xmm1,0x18(%rbp)
  call   <eigsrt@plt>
  movq   D11(%rip),%xmm0
  movsd  0x8(%rbp),%xmm1
  mov    %rbx,%rdx
  mov    0x8(%r12),%rdi
  mov    0x10(%r12),%rsi
  mov    $0x8,%eax
  xorpd  %xmm0,%xmm1
  mov    0x18(%r12),%rcx
  movsd  %xmm1,0x8(%rbp)
  movsd  0x10(%rbp),%xmm1
  xorpd  %xmm0,%xmm1
  movsd  %xmm1,0x10(%rbp)
  movsd  0x18(%rbp),%xmm1
  xorpd  %xmm0,%xmm1
  movsd  %xmm1,0x18(%rbp)
L5:
  movsd  0x0(%rbp,%rax,1),%xmm0
  add    $0x18,%rdx
  movsd  %xmm0,-0x8(%r15,%rax,1)
  movsd  (%rdi,%rax,1),%xmm0
  movsd  %xmm0,-0x18(%rdx)
  movsd  (%rsi,%rax,1),%xmm0
  movsd  %xmm0,-0x10(%rdx)
  movsd  (%rcx,%rax,1),%xmm0
  add    $0x8,%rax
  movsd  %xmm0,-0x8(%rdx)
  cmp    $0x20,%rax
  jne    L5
  mov    %r14,%rdi
  mov    $0x3,%r8d
  mov    $0x1,%ecx
  mov    $0x3,%edx
  mov    $0x1,%esi
  call   <free_dmatrix@plt>
  mov    %r12,%rdi
  mov    $0x3,%r8d
  mov    $0x1,%ecx
  mov    $0x3,%edx
  mov    $0x1,%esi
  call   <free_dmatrix@plt>
  mov    %rbp,%rdi
  mov    $0x3,%edx
  mov    $0x1,%esi
  call   <free_dvector@plt>
  add    $0x8,%rsp
  mov    $0x1,%eax
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
```

tool_calls: name: disassemble_data
arguments: tag: D11
datatype: byte[8]

tool: {"D11": [0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x80]}

assistant: ```c
int DerivePrincipalStresses(double stressTensor[3][3], double sp[3], double cn[3][3]) {
  double **a, **v, *d;
  int i, j;
  a = dmatrix(1, 3, 1, 3);
  v = dmatrix(1, 3, 1, 3);
  d = dvector(1, 3);
  a[1][1] = stressTensor[0][0];
  a[2][2] = stressTensor[1][1];
  a[3][3] = stressTensor[2][2];
  a[1][2] = stressTensor[0][1];
  a[1][3] = stressTensor[0][2];
  a[2][3] = stressTensor[1][2];
  a[2][1] = a[1][2];
  a[3][1] = a[1][3];
  a[3][2] = a[2][3];
  jacobi(a, d, v);
  d[1] *= -1.0f;
  d[2] *= -1.0f;
  d[3] *= -1.0f;
  eigsrt(d, v);
  d[1] *= -1.0f;
  d[2] *= -1.0f;
  d[3] *= -1.0f;
  for (i = 0; i < 3; i++) {
    sp[i] = d[i + 1];
    for (j = 0; j < 3; j++) {
      cn[i][j] = v[j + 1][i + 1];
    }
  }
  free_dmatrix(a, 1, 3, 1, 3);
  free_dmatrix(v, 1, 3, 1, 3);
  free_dvector(d, 1, 3);
  return (1);
}
```