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

```asm
<SetUp3>:
  endbr64
  push   %r15
  lea    0x90(%rcx),%rax
  push   %r14
  push   %r13
  push   %r12
  push   %rbp
  push   %rbx
  sub    $0xb8,%rsp
  mov    0xf0(%rsp),%rbp
  mov    0xf8(%rsp),%r14
  mov    %edi,0x14(%rsp)
  mov    %rdx,0x18(%rsp)
  mov    %r9,0xa0(%rsp)
  mov    %rbp,%r15
  mov    %r14,%rbp
  mov    %rcx,0x88(%rsp)
  mov    %r8,0x90(%rsp)
  mov    %rsi,0x98(%rsp)
  mov    %rax,0xa8(%rsp)
L103:
  mov    0x88(%rsp),%rax
  mov    0x90(%rsp),%rdi
  movsd  D30(%rip),%xmm4
  movsd  D30(%rip),%xmm5
  divsd  (%rax),%xmm4
  divsd  0x8(%rax),%xmm5
  movsd  D30(%rip),%xmm6
  divsd  0x10(%rax),%xmm6
  movsd  %xmm4,0x20(%rsp)
  mulsd  (%rdi),%xmm4
  mov    0x98(%rsp),%rax
  movsd  %xmm5,0x38(%rsp)
  mulsd  0x8(%rdi),%xmm5
  movsd  %xmm6,0x48(%rsp)
  mulsd  0x10(%rdi),%xmm6
  mov    0x4(%rax),%r11d
  mov    (%rax),%r13d
  mov    %r11d,0x2c(%rsp)
  movsd  %xmm4,0x30(%rsp)
  movsd  %xmm5,0x40(%rsp)
  movsd  %xmm6,0x50(%rsp)
  cmp    %r11d,%r13d
  jg     L42
  movslq %r13d,%r12
  lea    0x1(%r11),%r9d
  mov    0x14(%rsp),%r10d
  mov    0xa0(%rsp),%rcx
  mov    %r12,%rax
  movslq %r9d,%rdx
  mov    %r12,%rbx
  mov    %r9d,0x28(%rsp)
  shl    $0xd,%rax
  shl    $0x10,%rbx
  add    %rdx,%rax
  mov    %r10d,%edx
  add    %rcx,%rbx
  lea    (%rcx,%rax,8),%r8
  mov    %r13d,%eax
  sub    %r11d,%edx
  mov    $0x1,%ecx
  shr    %eax
  lea    -0x1(%rdx),%esi
  lea    -0x1(%rax),%edi
  mov    %r13d,%eax
  and    $0xfffffffe,%eax
  mov    %eax,0x6c(%rsp)
  lea    0x2(%r11),%eax
  cmp    %eax,%r10d
  cmovge %esi,%ecx
  cmp    $0x2,%edx
  sete   %dl
  cmp    %eax,%r10d
  mov    %rbp,%r10
  setl   %al
  mov    %ecx,%esi
  mov    %ecx,0x74(%rsp)
  shr    %ecx
  sub    $0x1,%ecx
  and    $0xfffffffe,%esi
  mov    %r13d,%ebp
  mov    %rcx,0x78(%rsp)
  mov    %r9d,%ecx
  add    %esi,%ecx
  or     %al,%dl
  mov    %r11d,%eax
  mov    %esi,0x70(%rsp)
  cmovne %r9d,%ecx
  sub    %r13d,%eax
  add    $0x1,%rdi
  mov    %dl,0x6b(%rsp)
  lea    0x8(,%rax,8),%rax
  shl    $0x4,%rdi
  mov    %r8,%r13
  mov    %rax,0x60(%rsp)
  lea    0x0(,%r12,8),%rax
  mov    %ecx,%r14d
  mov    %rax,0x58(%rsp)
  mov    %rdi,0x80(%rsp)
  nopl   0x0(%rax)
L100:
  pxor   %xmm0,%xmm0
  test   %ebp,%ebp
  jle    L16
  cmp    $0x1,%ebp
  je     L60
  mov    0x80(%rsp),%rax
  pxor   %xmm0,%xmm0
  lea    (%rax,%rbx,1),%rdx
  mov    %rbx,%rax
  nopl   0x0(%rax)
L27:
  movsd  (%rax),%xmm1
  add    $0x10,%rax
  addsd  %xmm1,%xmm0
  movsd  -0x8(%rax),%xmm1
  addsd  %xmm1,%xmm0
  cmp    %rax,%rdx
  jne    L27
  movslq 0x6c(%rsp),%rax
  cmp    %ebp,%eax
  je     L16
L4:
  addsd  (%rbx,%rax,8),%xmm0
L16:
  mov    0x14(%rsp),%ecx
  cmp    %ecx,0x28(%rsp)
  jge    L122
  cmpb   $0x0,0x6b(%rsp)
  jne    L25
  mov    0x78(%rsp),%rax
  lea    0x1(%rax),%rdx
  mov    %r13,%rax
  shl    $0x4,%rdx
  add    %r13,%rdx
  nopl   0x0(%rax)
L26:
  movsd  (%rax),%xmm1
  add    $0x10,%rax
  addsd  %xmm1,%xmm0
  movsd  -0x8(%rax),%xmm1
  addsd  %xmm1,%xmm0
  cmp    %rdx,%rax
  jne    L26
  mov    0x74(%rsp),%esi
  cmp    %esi,0x70(%rsp)
  je     L122
L25:
  movslq %r14d,%rax
  addsd  (%rbx,%rax,8),%xmm0
L122:
  movsd  D41(%rip),%xmm2
  mov    0x18(%rsp),%rax
  movapd %xmm0,%xmm1
  mulsd  %xmm0,%xmm2
  subsd  (%rax,%r12,8),%xmm1
  andpd  D16(%rip),%xmm1
  comisd %xmm2,%xmm1
  ja     L128
  mov    0x58(%rsp),%rax
  mov    0x60(%rsp),%rdx
  xor    %esi,%esi
  add    $0x10000,%r13
  xorpd  D54(%rip),%xmm0
  mov    %r10,0xf8(%rsp)
  lea    (%rbx,%rax,1),%rdi
  add    $0x10000,%rbx
  movsd  %xmm0,0x8(%rsp)
  call   <memset@plt>
  movsd  0x8(%rsp),%xmm0
  movsd  0x20(%rsp),%xmm1
  mov    0xf8(%rsp),%r10
  mulsd  %xmm0,%xmm1
  movsd  %xmm1,(%r15,%r12,8)
  movsd  0x30(%rsp),%xmm1
  mulsd  %xmm0,%xmm1
  movsd  %xmm1,(%r10,%r12,8)
  movsd  0x38(%rsp),%xmm1
  mulsd  %xmm0,%xmm1
  movsd  %xmm1,0x10000(%r15,%r12,8)
  movsd  0x40(%rsp),%xmm1
  mulsd  %xmm0,%xmm1
  movsd  %xmm1,0x10000(%r10,%r12,8)
  movsd  0x48(%rsp),%xmm1
  mulsd  %xmm0,%xmm1
  mulsd  0x50(%rsp),%xmm0
  movsd  %xmm1,0x20000(%r15,%r12,8)
  movsd  %xmm0,0x20000(%r10,%r12,8)
  add    $0x1,%r12
  cmp    %r12d,0x2c(%rsp)
  jge    L100
  mov    %r10,%rbp
L42:
  addq   $0x18,0x88(%rsp)
  mov    0x88(%rsp),%rax
  addq   $0x18,0x90(%rsp)
  addq   $0x8,0x98(%rsp)
  cmp    0xa8(%rsp),%rax
  jne    L103
  mov    $0x1,%eax
  jmp    L18
L60:
  pxor   %xmm0,%xmm0
  xor    %eax,%eax
  jmp    L4
L128:
  lea    D42(%rip),%rdi
  call   <puts@plt>
  xor    %eax,%eax
L18:
  add    $0xb8,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
```

tool_calls: name: parse_data
arguments: tag: D30
datatype: f64
name: parse_data
arguments: tag: D41
datatype: byte[8]
name: parse_data
arguments: tag: D16
datatype: qword[2]
name: parse_data
arguments: tag: D54
datatype: f32[4]
name: parse_data
arguments: tag: D42
datatype: string

tool: D30: 1

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

tool: {"D16": [0x7fffffffffffffff,0x0000000000000000]}

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

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

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)));
}
```