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

```asm
<mycs>:
  endbr64
  sub    $0x138,%rsp
  movsd  0x28(%rdi),%xmm5
  movsd  0x18(%rdi),%xmm0
  mov    %rdi,%rax
  mov    %fs:0x28,%rdx
  mov    %rdx,0x128(%rsp)
  xor    %edx,%edx
  movsd  0x8(%rdi),%xmm12
  movsd  0x20(%rdi),%xmm6
  mov    %rsi,%rcx
  movsd  0xa8(%rdi),%xmm8
  movsd  %xmm0,0x48(%rsp)
  addsd  %xmm5,%xmm0
  movsd  0x98(%rdi),%xmm2
  movsd  %xmm5,0x50(%rsp)
  movsd  0xb8(%rdi),%xmm5
  movsd  0x38(%rdi),%xmm11
  movsd  %xmm8,0x58(%rsp)
  movsd  0xc8(%rdi),%xmm1
  addsd  %xmm5,%xmm8
  movapd %xmm0,%xmm4
  movsd  %xmm5,0x60(%rsp)
  movsd  0xb0(%rdi),%xmm5
  addsd  %xmm12,%xmm4
  movsd  %xmm6,0x8(%rsp)
  movsd  %xmm5,0x10(%rsp)
  movapd %xmm8,%xmm3
  addsd  %xmm2,%xmm3
  addsd  %xmm11,%xmm4
  addsd  %xmm1,%xmm3
  addsd  %xmm6,%xmm4
  addsd  %xmm5,%xmm3
  comisd %xmm3,%xmm4
  ja     L46
  movsd  D79(%rip),%xmm5
  movsd  D24(%rip),%xmm4
L56:
  movapd %xmm12,%xmm6
  movapd %xmm11,%xmm7
  movsd  0x50(%rax),%xmm3
  movsd  0x80(%rax),%xmm10
  addsd  %xmm2,%xmm6
  addsd  %xmm1,%xmm7
  movapd %xmm2,%xmm14
  movsd  D56(%rip),%xmm9
  movsd  %xmm3,0x18(%rsp)
  addsd  %xmm1,%xmm14
  movsd  %xmm10,0x20(%rsp)
  movsd  %xmm6,0x78(%rsp)
  addsd  %xmm3,%xmm6
  movsd  0x60(%rax),%xmm3
  movsd  %xmm7,0x80(%rsp)
  addsd  %xmm10,%xmm7
  movsd  0x10(%rsp),%xmm10
  movsd  %xmm3,0x28(%rsp)
  movsd  %xmm14,0x90(%rsp)
  movapd %xmm6,%xmm13
  movapd %xmm5,%xmm6
  movsd  0x50(%rsp),%xmm5
  subsd  %xmm7,%xmm13
  addsd  0x60(%rsp),%xmm5
  movsd  0x70(%rax),%xmm7
  movsd  %xmm7,0x30(%rsp)
  mulsd  %xmm9,%xmm13
  addsd  %xmm5,%xmm7
  unpcklpd %xmm13,%xmm6
  movaps %xmm6,0xa0(%rsp)
  movsd  0x48(%rsp),%xmm6
  addsd  0x58(%rsp),%xmm6
  movapd %xmm6,%xmm15
  addsd  %xmm3,%xmm15
  movapd %xmm14,%xmm3
  movsd  0x8(%rsp),%xmm14
  addsd  %xmm3,%xmm10
  movsd  0x58(%rax),%xmm3
  subsd  %xmm7,%xmm15
  movapd %xmm12,%xmm7
  movsd  %xmm3,0x38(%rsp)
  addsd  %xmm11,%xmm7
  mulsd  %xmm9,%xmm15
  movsd  %xmm7,0x88(%rsp)
  addsd  %xmm7,%xmm14
  movsd  0x48(%rax),%xmm7
  movsd  %xmm7,0x68(%rsp)
  addsd  %xmm3,%xmm7
  subsd  %xmm10,%xmm14
  movsd  0x88(%rax),%xmm10
  movapd %xmm7,%xmm3
  movsd  0x78(%rax),%xmm7
  movsd  %xmm10,0x40(%rsp)
  addsd  %xmm3,%xmm2
  mulsd  %xmm9,%xmm14
  movsd  %xmm7,0x70(%rsp)
  addsd  %xmm10,%xmm7
  addsd  %xmm12,%xmm2
  addsd  0x18(%rsp),%xmm2
  addsd  %xmm7,%xmm1
  addsd  %xmm11,%xmm1
  addsd  0x20(%rsp),%xmm1
  movapd %xmm4,%xmm11
  comisd %xmm1,%xmm2
  ja     L45
  movsd  D79(%rip),%xmm11
L45:
  movapd %xmm14,%xmm2
  movsd  0x68(%rsp),%xmm1
  addsd  0x20(%rsp),%xmm7
  unpcklpd %xmm11,%xmm2
  addsd  0x8(%rsp),%xmm0
  addsd  0x10(%rsp),%xmm8
  movaps %xmm2,0xb0(%rsp)
  movsd  0x28(%rsp),%xmm2
  addsd  0x18(%rsp),%xmm3
  addsd  %xmm1,%xmm6
  movsd  0x70(%rsp),%xmm10
  movsd  0x30(%rsp),%xmm11
  addsd  %xmm1,%xmm2
  addsd  0x38(%rsp),%xmm11
  movsd  0x28(%rsp),%xmm1
  subsd  %xmm7,%xmm3
  addsd  0x40(%rsp),%xmm11
  subsd  %xmm8,%xmm0
  addsd  %xmm10,%xmm6
  addsd  %xmm10,%xmm2
  addsd  0x38(%rsp),%xmm5
  addsd  0x40(%rsp),%xmm5
  unpcklpd %xmm3,%xmm0
  movsd  0x30(%rsp),%xmm3
  addsd  %xmm6,%xmm1
  subsd  %xmm11,%xmm2
  mulpd  D56(%rip),%xmm0
  movapd %xmm4,%xmm11
  addsd  %xmm5,%xmm3
  mulsd  %xmm9,%xmm2
  comisd %xmm3,%xmm1
  movhpd %xmm0,0x98(%rsp)
  movapd %xmm0,%xmm7
  movapd %xmm2,%xmm12
  ja     L26
  movsd  D79(%rip),%xmm11
L26:
  movq   D21(%rip),%xmm2
  movapd %xmm15,%xmm3
  movapd %xmm14,%xmm1
  xor    %edx,%edx
  andpd  %xmm2,%xmm13
  andpd  %xmm2,%xmm3
  andpd  %xmm2,%xmm1
  movapd %xmm13,%xmm8
  addsd  %xmm4,%xmm1
  addsd  %xmm4,%xmm8
  addsd  %xmm3,%xmm8
  movapd 0xa0(%rsp),%xmm3
  divsd  %xmm8,%xmm15
  movapd %xmm8,%xmm9
  movapd %xmm12,%xmm8
  andpd  %xmm2,%xmm8
  unpcklpd %xmm9,%xmm9
  addsd  %xmm8,%xmm1
  movapd 0xb0(%rsp),%xmm8
  divpd  %xmm9,%xmm3
  movapd %xmm1,%xmm9
  unpcklpd %xmm9,%xmm9
  movsd  %xmm15,0xd0(%rsp)
  divsd  %xmm1,%xmm12
  movapd %xmm7,%xmm1
  movsd  0x98(%rsp),%xmm7
  andpd  %xmm2,%xmm1
  andpd  %xmm2,%xmm7
  addsd  %xmm7,%xmm1
  addsd  %xmm4,%xmm1
  movaps %xmm3,0xc0(%rsp)
  andpd  %xmm2,%xmm3
  movapd %xmm1,%xmm4
  unpcklpd %xmm4,%xmm4
  divpd  %xmm9,%xmm8
  movsd  %xmm12,0xe8(%rsp)
  divpd  %xmm4,%xmm0
  movups %xmm8,0xd8(%rsp)
  unpckhpd %xmm8,%xmm8
  andpd  %xmm2,%xmm8
  comisd %xmm3,%xmm8
  divsd  %xmm1,%xmm11
  movaps %xmm0,0xf0(%rsp)
  movsd  %xmm11,0x100(%rsp)
  movapd %xmm11,%xmm0
  andpd  %xmm2,%xmm0
  jbe    L69
  movapd %xmm8,%xmm3
  mov    $0x1,%edx
L69:
  movsd  (%rax),%xmm8
  ucomisd %xmm3,%xmm0
  movsd  0x30(%rax),%xmm7
  addsd  %xmm6,%xmm6
  movsd  0x10(%rax),%xmm4
  movsd  0x48(%rsp),%xmm0
  addsd  %xmm5,%xmm5
  mov    $0x2,%esi
  movapd %xmm8,%xmm1
  movsd  0x40(%rax),%xmm13
  addsd  0x88(%rsp),%xmm0
  addsd  %xmm7,%xmm1
  addsd  0x50(%rsp),%xmm0
  cmova  %esi,%edx
  movsd  D68(%rip),%xmm9
  movsd  0x90(%rax),%xmm15
  addsd  %xmm4,%xmm8
  movsd  0xc0(%rax),%xmm14
  addsd  %xmm13,%xmm7
  movsd  0xa0(%rax),%xmm12
  movsd  0xd0(%rax),%xmm11
  movslq %edx,%rax
  movapd %xmm1,%xmm3
  addsd  %xmm0,%xmm0
  addsd  %xmm15,%xmm1
  mov    %rax,%rdx
  addsd  %xmm4,%xmm3
  addsd  %xmm15,%xmm8
  shl    $0x5,%rdx
  addsd  %xmm14,%xmm7
  addsd  %xmm13,%xmm4
  addsd  %xmm14,%xmm1
  addsd  %xmm13,%xmm3
  addsd  %xmm12,%xmm8
  addsd  %xmm11,%xmm7
  addsd  %xmm12,%xmm4
  addsd  %xmm6,%xmm1
  movsd  0x28(%rsp),%xmm6
  addsd  %xmm3,%xmm0
  movsd  0x8(%rsp),%xmm3
  addsd  %xmm11,%xmm4
  mulsd  %xmm9,%xmm6
  mulsd  %xmm9,%xmm3
  addsd  %xmm4,%xmm5
  movsd  0x30(%rsp),%xmm4
  mulsd  %xmm9,%xmm4
  addsd  %xmm6,%xmm1
  addsd  %xmm3,%xmm0
  movapd %xmm15,%xmm3
  addsd  %xmm14,%xmm3
  addsd  %xmm12,%xmm3
  addsd  %xmm4,%xmm5
  addsd  %xmm11,%xmm3
  subsd  %xmm5,%xmm1
  movapd %xmm3,%xmm10
  movsd  0x58(%rsp),%xmm3
  addsd  0x90(%rsp),%xmm3
  addsd  0x60(%rsp),%xmm3
  addsd  %xmm3,%xmm3
  addsd  %xmm10,%xmm3
  movsd  0x10(%rsp),%xmm10
  mulsd  %xmm9,%xmm10
  addsd  %xmm10,%xmm3
  subsd  %xmm3,%xmm0
  movsd  0x78(%rsp),%xmm3
  addsd  0x68(%rsp),%xmm3
  addsd  0x38(%rsp),%xmm3
  movapd %xmm0,%xmm4
  addsd  %xmm3,%xmm3
  andpd  %xmm2,%xmm4
  addsd  %xmm8,%xmm3
  movsd  0x18(%rsp),%xmm8
  mulsd  %xmm9,%xmm8
  addsd  %xmm8,%xmm3
  movsd  0x80(%rsp),%xmm8
  addsd  0x70(%rsp),%xmm8
  addsd  0x40(%rsp),%xmm8
  addsd  %xmm8,%xmm8
  addsd  %xmm7,%xmm8
  movsd  0x20(%rsp),%xmm7
  mulsd  %xmm9,%xmm7
  addsd  %xmm7,%xmm8
  subsd  %xmm8,%xmm3
  movapd %xmm3,%xmm5
  andpd  %xmm2,%xmm5
  addsd  %xmm5,%xmm4
  movapd %xmm1,%xmm5
  andpd  %xmm2,%xmm5
  addsd  %xmm5,%xmm4
  addsd  D36(%rip),%xmm4
  divsd  %xmm4,%xmm0
  divsd  %xmm4,%xmm3
  movsd  %xmm0,0x108(%rsp)
  andpd  %xmm2,%xmm0
  divsd  %xmm4,%xmm1
  movsd  %xmm3,0x110(%rsp)
  andpd  %xmm2,%xmm3
  maxsd  %xmm0,%xmm3
  movapd %xmm1,%xmm4
  andpd  %xmm2,%xmm4
  maxsd  %xmm3,%xmm4
  movsd  0xc0(%rsp,%rdx,1),%xmm3
  andpd  %xmm3,%xmm2
  comisd %xmm4,%xmm2
  ja     L8
  lea    (%rax,%rax,2),%rdx
  movsd  0xd0(%rsp,%rdx,8),%xmm1
L9:
  lea    (%rax,%rax,2),%rax
  movsd  %xmm1,0x10(%rcx)
  movupd 0xc0(%rsp,%rax,8),%xmm6
  movups %xmm6,(%rcx)
  mov    0x128(%rsp),%rax
  sub    %fs:0x28,%rax
  jne    L21
  add    $0x138,%rsp
  ret
  nopw   0x0(%rax,%rax,1)
L46:
  movsd  D24(%rip),%xmm4
  movapd %xmm4,%xmm5
  jmp    L56
  nopl   0x0(%rax)
L8:
  mov    $0x3,%eax
  jmp    L9
L21:
  call   <__stack_chk_fail@plt>
```

tool_calls: name: read_data
arguments: tag: D79
type_guessed: byte[8]
name: read_data
arguments: tag: D24
type_guessed: f64
name: read_data
arguments: tag: D56
type_guessed: f64[2]
name: read_data
arguments: tag: D21
type_guessed: qword
name: read_data
arguments: tag: D68
type_guessed: f64
name: read_data
arguments: tag: D36
type_guessed: byte[8]

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

tool: D24: 1

tool: {"D56": [0.5,0.5]}

tool: D21: 9223372036854775807

tool: D68: 4

tool: {"D36": [0xa0,0xc2,0xeb,0xfe,0x4b,0x48,0xb4,0x39]}

assistant: ```c
void mycs(double c[3][3][3], double mxyz[3]) {
  double m1, m2, m[4][3], t0, t1, t2;
  int cn;
  m1 = c[0][1][0] + c[0][1][2] + c[0][0][1] + c[0][2][1] + c[0][1][1];
  m2 = c[2][1][0] + c[2][1][2] + c[2][0][1] + c[2][2][1] + c[2][1][1];
  m[0][0] = m1 > m2 ? 1. : -1.;
  m1 = c[0][0][1] + c[2][0][1] + c[1][0][1];
  m2 = c[0][2][1] + c[2][2][1] + c[1][2][1];
  m[0][1] = 0.5 * (m1 - m2);
  m1 = c[0][1][0] + c[2][1][0] + c[1][1][0];
  m2 = c[0][1][2] + c[2][1][2] + c[1][1][2];
  m[0][2] = 0.5 * (m1 - m2);
  m1 = c[0][0][1] + c[0][2][1] + c[0][1][1];
  m2 = c[2][0][1] + c[2][2][1] + c[2][1][1];
  m[1][0] = 0.5 * (m1 - m2);
  m1 = c[1][0][0] + c[1][0][2] + c[2][0][1] + c[0][0][1] + c[1][0][1];
  m2 = c[1][2][0] + c[1][2][2] + c[2][2][1] + c[0][2][1] + c[1][2][1];
  m[1][1] = m1 > m2 ? 1. : -1.;
  m1 = c[1][0][0] + c[1][1][0] + c[1][2][0];
  m2 = c[1][0][2] + c[1][1][2] + c[1][2][2];
  m[1][2] = 0.5 * (m1 - m2);
  m1 = c[0][1][0] + c[0][1][2] + c[0][1][1];
  m2 = c[2][1][0] + c[2][1][2] + c[2][1][1];
  m[2][0] = 0.5 * (m1 - m2);
  m1 = c[1][0][0] + c[1][0][2] + c[1][0][1];
  m2 = c[1][2][0] + c[1][2][2] + c[1][2][1];
  m[2][1] = 0.5 * (m1 - m2);
  m1 = c[0][1][0] + c[2][1][0] + c[1][0][0] + c[1][2][0] + c[1][1][0];
  m2 = c[0][1][2] + c[2][1][2] + c[1][0][2] + c[1][2][2] + c[1][1][2];
  m[2][2] = m1 > m2 ? 1. : -1.;
  t0 = fabs(m[0][0]) + fabs(m[0][1]) + fabs(m[0][2]);
  m[0][0] /= t0;
  m[0][1] /= t0;
  m[0][2] /= t0;
  t0 = fabs(m[1][0]) + fabs(m[1][1]) + fabs(m[1][2]);
  m[1][0] /= t0;
  m[1][1] /= t0;
  m[1][2] /= t0;
  t0 = fabs(m[2][0]) + fabs(m[2][1]) + fabs(m[2][2]);
  m[2][0] /= t0;
  m[2][1] /= t0;
  m[2][2] /= t0;
  t0 = fabs(m[0][0]);
  t1 = fabs(m[1][1]);
  t2 = fabs(m[2][2]);
  cn = 0;
  if (t1 > t0) {
    t0 = t1;
    cn = 1;
  }
  if (t2 > t0)
    cn = 2;
  m1 = c[0][0][0] + c[0][2][0] + c[0][0][2] + c[0][2][2] + 2. * (c[0][0][1] + c[0][2][1] + c[0][1][0] + c[0][1][2]) + 4. * c[0][1][1];
  m2 = c[2][0][0] + c[2][2][0] + c[2][0][2] + c[2][2][2] + 2. * (c[2][0][1] + c[2][2][1] + c[2][1][0] + c[2][1][2]) + 4. * c[2][1][1];
  m[3][0] = m1 - m2;
  m1 = c[0][0][0] + c[0][0][2] + c[2][0][0] + c[2][0][2] + 2. * (c[0][0][1] + c[2][0][1] + c[1][0][0] + c[1][0][2]) + 4. * c[1][0][1];
  m2 = c[0][2][0] + c[0][2][2] + c[2][2][0] + c[2][2][2] + 2. * (c[0][2][1] + c[2][2][1] + c[1][2][0] + c[1][2][2]) + 4. * c[1][2][1];
  m[3][1] = m1 - m2;
  m1 = c[0][0][0] + c[0][2][0] + c[2][0][0] + c[2][2][0] + 2. * (c[0][1][0] + c[2][1][0] + c[1][0][0] + c[1][2][0]) + 4. * c[1][1][0];
  m2 = c[0][0][2] + c[0][2][2] + c[2][0][2] + c[2][2][2] + 2. * (c[0][1][2] + c[2][1][2] + c[1][0][2] + c[1][2][2]) + 4. * c[1][1][2];
  m[3][2] = m1 - m2;
  t0 = fabs(m[3][0]) + fabs(m[3][1]) + fabs(m[3][2]) + 1.e-30;
  m[3][0] /= t0;
  m[3][1] /= t0;
  m[3][2] /= t0;
  t0 = fabs(m[3][0]);
  t1 = fabs(m[3][1]);
  t2 = fabs(m[3][2]);
  if (t1 > t0)
    t0 = t1;
  if (t2 > t0)
    t0 = t2;
  if (fabs(m[cn][cn]) > t0)
    cn = 3;
  mxyz[0] = m[cn][0];
  mxyz[1] = m[cn][1];
  mxyz[2] = m[cn][2];
  return;
}
```