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

```asm
<SetUp3>:
  endbr64
  push   %rbp
  mov    %rsp,%rbp
  sub    $0xa0,%rsp
  mov    %edi,-0x64(%rbp)
  mov    %rsi,-0x70(%rbp)
  mov    %rdx,-0x78(%rbp)
  mov    %rcx,-0x80(%rbp)
  mov    %r8,-0x88(%rbp)
  mov    %r9,-0x90(%rbp)
  mov    0x10(%rbp),%rax
  mov    %rax,-0x98(%rbp)
  mov    0x18(%rbp),%rax
  mov    %rax,-0xa0(%rbp)
  mov    %fs:0x28,%rax
  mov    %rax,-0x8(%rbp)
  xor    %eax,%eax
  movl   $0x0,-0x54(%rbp)
  jmp    L114
L75:
  movl   $0x0,-0x58(%rbp)
  jmp    L148
L21:
  mov    -0x54(%rbp),%eax
  movslq %eax,%rdx
  mov    %rdx,%rax
  add    %rax,%rax
  add    %rdx,%rax
  shl    $0x3,%rax
  mov    %rax,%rdx
  mov    -0x80(%rbp),%rax
  add    %rax,%rdx
  mov    -0x58(%rbp),%eax
  cltq
  movsd  (%rdx,%rax,8),%xmm1
  movsd  D21(%rip),%xmm0
  divsd  %xmm1,%xmm0
  mov    -0x58(%rbp),%eax
  cltq
  movsd  %xmm0,-0x40(%rbp,%rax,8)
  mov    -0x54(%rbp),%eax
  movslq %eax,%rdx
  mov    %rdx,%rax
  add    %rax,%rax
  add    %rdx,%rax
  shl    $0x3,%rax
  mov    %rax,%rdx
  mov    -0x88(%rbp),%rax
  add    %rax,%rdx
  mov    -0x58(%rbp),%eax
  cltq
  movsd  (%rdx,%rax,8),%xmm1
  mov    -0x58(%rbp),%eax
  cltq
  movsd  -0x40(%rbp,%rax,8),%xmm0
  mulsd  %xmm1,%xmm0
  mov    -0x58(%rbp),%eax
  cltq
  movsd  %xmm0,-0x20(%rbp,%rax,8)
  addl   $0x1,-0x58(%rbp)
L148:
  cmpl   $0x2,-0x58(%rbp)
  jle    L21
  mov    -0x54(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x70(%rbp),%rax
  add    %rdx,%rax
  mov    (%rax),%eax
  mov    %eax,-0x50(%rbp)
  jmp    L81
L20:
  pxor   %xmm0,%xmm0
  movsd  %xmm0,-0x48(%rbp)
  movl   $0x0,-0x4c(%rbp)
  jmp    L107
L72:
  mov    -0x50(%rbp),%eax
  cltq
  shl    $0x10,%rax
  mov    %rax,%rdx
  mov    -0x90(%rbp),%rax
  add    %rax,%rdx
  mov    -0x4c(%rbp),%eax
  cltq
  movsd  (%rdx,%rax,8),%xmm0
  movsd  -0x48(%rbp),%xmm1
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x48(%rbp)
  addl   $0x1,-0x4c(%rbp)
L107:
  mov    -0x54(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x70(%rbp),%rax
  add    %rdx,%rax
  mov    (%rax),%eax
  cmp    %eax,-0x4c(%rbp)
  jl     L72
  mov    -0x54(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x70(%rbp),%rax
  add    %rdx,%rax
  mov    0x4(%rax),%eax
  add    $0x1,%eax
  mov    %eax,-0x4c(%rbp)
  jmp    L106
L17:
  mov    -0x50(%rbp),%eax
  cltq
  shl    $0x10,%rax
  mov    %rax,%rdx
  mov    -0x90(%rbp),%rax
  add    %rax,%rdx
  mov    -0x4c(%rbp),%eax
  cltq
  movsd  (%rdx,%rax,8),%xmm0
  movsd  -0x48(%rbp),%xmm1
  addsd  %xmm1,%xmm0
  movsd  %xmm0,-0x48(%rbp)
  addl   $0x1,-0x4c(%rbp)
L106:
  mov    -0x4c(%rbp),%eax
  cmp    -0x64(%rbp),%eax
  jl     L17
  mov    -0x50(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x78(%rbp),%rax
  add    %rdx,%rax
  movsd  (%rax),%xmm1
  movsd  -0x48(%rbp),%xmm0
  subsd  %xmm1,%xmm0
  movq   D39(%rip),%xmm1
  andpd  %xmm1,%xmm0
  movsd  -0x48(%rbp),%xmm2
  movsd  D6(%rip),%xmm1
  mulsd  %xmm2,%xmm1
  comisd %xmm1,%xmm0
  jbe    L123
  lea    D7(%rip),%rax
  mov    %rax,%rdi
  call   <puts@plt>
  mov    $0x0,%eax
  jmp    L54
L123:
  movsd  -0x48(%rbp),%xmm0
  movq   D32(%rip),%xmm1
  xorpd  %xmm1,%xmm0
  movsd  %xmm0,-0x48(%rbp)
  mov    -0x54(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x70(%rbp),%rax
  add    %rdx,%rax
  mov    (%rax),%eax
  mov    %eax,-0x4c(%rbp)
  jmp    L90
L105:
  mov    -0x50(%rbp),%eax
  cltq
  shl    $0x10,%rax
  mov    %rax,%rdx
  mov    -0x90(%rbp),%rax
  add    %rax,%rdx
  mov    -0x4c(%rbp),%eax
  cltq
  pxor   %xmm0,%xmm0
  movsd  %xmm0,(%rdx,%rax,8)
  addl   $0x1,-0x4c(%rbp)
L90:
  mov    -0x54(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x70(%rbp),%rax
  add    %rdx,%rax
  mov    0x4(%rax),%eax
  cmp    %eax,-0x4c(%rbp)
  jle    L105
  movl   $0x0,-0x58(%rbp)
  jmp    L165
L126:
  mov    -0x58(%rbp),%eax
  cltq
  movsd  -0x40(%rbp,%rax,8),%xmm0
  mov    -0x58(%rbp),%eax
  cltq
  shl    $0x10,%rax
  mov    %rax,%rdx
  mov    -0x98(%rbp),%rax
  add    %rax,%rdx
  mulsd  -0x48(%rbp),%xmm0
  mov    -0x50(%rbp),%eax
  cltq
  movsd  %xmm0,(%rdx,%rax,8)
  mov    -0x58(%rbp),%eax
  cltq
  movsd  -0x20(%rbp,%rax,8),%xmm0
  mov    -0x58(%rbp),%eax
  cltq
  shl    $0x10,%rax
  mov    %rax,%rdx
  mov    -0xa0(%rbp),%rax
  add    %rax,%rdx
  mulsd  -0x48(%rbp),%xmm0
  mov    -0x50(%rbp),%eax
  cltq
  movsd  %xmm0,(%rdx,%rax,8)
  addl   $0x1,-0x58(%rbp)
L165:
  cmpl   $0x2,-0x58(%rbp)
  jle    L126
  addl   $0x1,-0x50(%rbp)
L81:
  mov    -0x54(%rbp),%eax
  cltq
  lea    0x0(,%rax,8),%rdx
  mov    -0x70(%rbp),%rax
  add    %rdx,%rax
  mov    0x4(%rax),%eax
  cmp    %eax,-0x50(%rbp)
  jle    L20
  addl   $0x1,-0x54(%rbp)
L114:
  cmpl   $0x5,-0x54(%rbp)
  jle    L75
  mov    $0x1,%eax
L54:
  mov    -0x8(%rbp),%rdx
  sub    %fs:0x28,%rdx
  je     L153
  call   <__stack_chk_fail@plt>
L153:
  leave
  ret
```

tool_calls: name: unpack_data
arguments: data_label: D21
type_guessed: f64
name: unpack_data
arguments: data_label: D39
type_guessed: qword
name: unpack_data
arguments: data_label: D6
type_guessed: byte[8]
name: unpack_data
arguments: data_label: D7
type_guessed: string
name: unpack_data
arguments: data_label: D32
type_guessed: f32[2]

tool: D21: 1

tool: D39: 9223372036854775807

tool: {"D6": [0x95,0xd6,0x26,0xe8,0x0b,0x2e,0x01,0x3e]}

tool: D7: Total form factor is too far from unity.

tool: {"D32": [0.0,-0.0]}

assistant: ```c
SetUp3(npatch, loop, area, rho, emiss, coeff, diag, rhs) int npatch, loop[][2];
double area[], rho[][3], emiss[][3], coeff[][8192], diag[][8192], rhs[][8192];
{
  int j, iface, ipatch, jpatch;
  double tmp1, vtmp1[3], vtmp2[3];
  for (iface = 0; iface < 6; iface++) {
    for (j = 0; j < 3; j++) {
      vtmp1[j] = 1.0 / rho[iface][j];
      vtmp2[j] = emiss[iface][j] * vtmp1[j];
    }
    for (ipatch = loop[iface][0]; ipatch <= loop[iface][1]; ipatch++) {
      tmp1 = 0.0;
      for (jpatch = 0; jpatch < loop[iface][0]; jpatch++) {
        tmp1 += coeff[ipatch][jpatch];
      }
      for (jpatch = loop[iface][1] + 1; jpatch < npatch; jpatch++) {
        tmp1 += coeff[ipatch][jpatch];
      }
      if (fabs(tmp1 - area[ipatch]) > (0.5e-9 * tmp1)) {
        printf("Total form factor is too far from unity.\n");
        return ((1 == 0));
      }
      tmp1 = -tmp1;
      for (jpatch = loop[iface][0]; jpatch <= loop[iface][1]; jpatch++) {
        coeff[ipatch][jpatch] = 0.0;
      }
      for (j = 0; j < 3; j++) {
        diag[j][ipatch] = vtmp1[j] * tmp1;
        rhs[j][ipatch] = vtmp2[j] * tmp1;
      }
    }
  }
  return ((!(1 == 0)));
}
```