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


x86 assembly (clang, -O0):
<SetUp3>:
push   %rbp
mov    %rsp,%rbp
sub    $0x80,%rsp
mov    0x18(%rbp),%rax
mov    0x10(%rbp),%rax
mov    %edi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %rcx,-0x20(%rbp)
mov    %r8,-0x28(%rbp)
mov    %r9,-0x30(%rbp)
movl   $0x0,-0x38(%rbp)
cmpl   $0x6,-0x38(%rbp)
jge    13d2 <SetUp3+0x2c2>
movl   $0x0,-0x34(%rbp)
cmpl   $0x3,-0x34(%rbp)
jge    11c0 <SetUp3+0xb0>
mov    -0x20(%rbp),%rax
movslq -0x38(%rbp),%rcx
imul   $0x18,%rcx,%rcx
add    %rcx,%rax
movslq -0x34(%rbp),%rcx
movsd  0xe91(%rip),%xmm0
divsd  (%rax,%rcx,8),%xmm0
movslq -0x34(%rbp),%rax
movsd  %xmm0,-0x60(%rbp,%rax,8)
mov    -0x28(%rbp),%rax
movslq -0x38(%rbp),%rcx
imul   $0x18,%rcx,%rcx
add    %rcx,%rax
movslq -0x34(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movslq -0x34(%rbp),%rax
mulsd  -0x60(%rbp,%rax,8),%xmm0
movslq -0x34(%rbp),%rax
movsd  %xmm0,-0x80(%rbp,%rax,8)
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    1152 <SetUp3+0x42>
mov    -0x10(%rbp),%rax
movslq -0x38(%rbp),%rcx
mov    (%rax,%rcx,8),%eax
mov    %eax,-0x3c(%rbp)
mov    -0x3c(%rbp),%eax
mov    -0x10(%rbp),%rcx
movslq -0x38(%rbp),%rdx
cmp    0x4(%rcx,%rdx,8),%eax
jg     13bf <SetUp3+0x2af>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x48(%rbp)
movl   $0x0,-0x40(%rbp)
mov    -0x40(%rbp),%eax
mov    -0x10(%rbp),%rcx
movslq -0x38(%rbp),%rdx
cmp    (%rcx,%rdx,8),%eax
jge    1236 <SetUp3+0x126>
mov    -0x30(%rbp),%rax
movslq -0x3c(%rbp),%rcx
shl    $0x10,%rcx
add    %rcx,%rax
movslq -0x40(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
addsd  -0x48(%rbp),%xmm0
movsd  %xmm0,-0x48(%rbp)
mov    -0x40(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x40(%rbp)
jmp    11f2 <SetUp3+0xe2>
mov    -0x10(%rbp),%rax
movslq -0x38(%rbp),%rcx
mov    0x4(%rax,%rcx,8),%eax
add    $0x1,%eax
mov    %eax,-0x40(%rbp)
mov    -0x40(%rbp),%eax
cmp    -0x8(%rbp),%eax
jge    1284 <SetUp3+0x174>
mov    -0x30(%rbp),%rax
movslq -0x3c(%rbp),%rcx
shl    $0x10,%rcx
add    %rcx,%rax
movslq -0x40(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
addsd  -0x48(%rbp),%xmm0
movsd  %xmm0,-0x48(%rbp)
mov    -0x40(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x40(%rbp)
jmp    1248 <SetUp3+0x138>
movsd  -0x48(%rbp),%xmm0
mov    -0x18(%rbp),%rax
movslq -0x3c(%rbp),%rcx
subsd  (%rax,%rcx,8),%xmm0
movaps 0xd73(%rip),%xmm1
pand   %xmm1,%xmm0
movsd  0xd57(%rip),%xmm1
mulsd  -0x48(%rbp),%xmm1
ucomisd %xmm1,%xmm0
jbe    12d2 <SetUp3+0x1c2>
lea    0xd61(%rip),%rdi
mov    $0x0,%al
call   1030 <printf@plt>
movl   $0x0,-0x4(%rbp)
jmp    13d9 <SetUp3+0x2c9>
movsd  -0x48(%rbp),%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
movsd  %xmm0,-0x48(%rbp)
mov    -0x10(%rbp),%rax
movslq -0x38(%rbp),%rcx
mov    (%rax,%rcx,8),%eax
mov    %eax,-0x40(%rbp)
mov    -0x40(%rbp),%eax
mov    -0x10(%rbp),%rcx
movslq -0x38(%rbp),%rdx
cmp    0x4(%rcx,%rdx,8),%eax
jg     133f <SetUp3+0x22f>
mov    -0x30(%rbp),%rax
movslq -0x3c(%rbp),%rcx
shl    $0x10,%rcx
add    %rcx,%rax
movslq -0x40(%rbp),%rcx
xorps  %xmm0,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x40(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x40(%rbp)
jmp    1301 <SetUp3+0x1f1>
movl   $0x0,-0x34(%rbp)
cmpl   $0x3,-0x34(%rbp)
jge    13ac <SetUp3+0x29c>
movslq -0x34(%rbp),%rax
movsd  -0x60(%rbp,%rax,8),%xmm0
mulsd  -0x48(%rbp),%xmm0
mov    0x10(%rbp),%rax
movslq -0x34(%rbp),%rcx
shl    $0x10,%rcx
add    %rcx,%rax
movslq -0x3c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
movslq -0x34(%rbp),%rax
movsd  -0x80(%rbp,%rax,8),%xmm0
mulsd  -0x48(%rbp),%xmm0
mov    0x18(%rbp),%rax
movslq -0x34(%rbp),%rcx
shl    $0x10,%rcx
add    %rcx,%rax
movslq -0x3c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    1346 <SetUp3+0x236>
jmp    13b1 <SetUp3+0x2a1>
mov    -0x3c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x3c(%rbp)
jmp    11ce <SetUp3+0xbe>
jmp    13c4 <SetUp3+0x2b4>
mov    -0x38(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x38(%rbp)
jmp    1141 <SetUp3+0x31>
movl   $0x1,-0x4(%rbp)
mov    -0x4(%rbp),%eax
add    $0x80,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<SetUp3>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0xd8,%rsp
mov    %r8,0x50(%rsp)
mov    %rcx,0x48(%rsp)
mov    %rdx,0xb8(%rsp)
mov    %rsi,0x58(%rsp)
lea    -0x2(%rdi),%eax
mov    %eax,0x8(%rsp)
mov    %r9,0x28(%rsp)
lea    0x18(%r9),%rax
mov    %rax,0x40(%rsp)
mov    %rdi,0x30(%rsp)
mov    %edi,%esi
xor    %ecx,%ecx
mov    %rsi,0x68(%rsp)
jmp    1183 <SetUp3+0x63>
nopl   (%rax)
mov    0x18(%rsp),%rcx
add    $0x1,%rcx
cmp    $0x6,%rcx
je     1503 <SetUp3+0x3e3>
mov    0x58(%rsp),%rax
mov    (%rax,%rcx,8),%edx
mov    %rcx,0x18(%rsp)
mov    0x4(%rax,%rcx,8),%edi
mov    %rdx,0x38(%rsp)
cmp    %edx,%edi
jl     1170 <SetUp3+0x50>
mov    0x18(%rsp),%rax
lea    0x0(,%rax,8),%rax
lea    (%rax,%rax,2),%rax
movsd  0xe4a(%rip),%xmm0
movapd %xmm0,%xmm1
mov    0x48(%rsp),%rdx
divsd  (%rdx,%rax,1),%xmm1
mov    0x50(%rsp),%rcx
movsd  (%rcx,%rax,1),%xmm2
movsd  %xmm1,0xb0(%rsp)
mulsd  %xmm1,%xmm2
movsd  %xmm2,0xa8(%rsp)
movapd %xmm0,%xmm1
divsd  0x8(%rdx,%rax,1),%xmm1
movsd  0x8(%rcx,%rax,1),%xmm2
movsd  %xmm1,0xa0(%rsp)
mulsd  %xmm1,%xmm2
movsd  %xmm2,0x98(%rsp)
divsd  0x10(%rdx,%rax,1),%xmm0
movsd  0x10(%rcx,%rax,1),%xmm1
movsd  %xmm0,0x90(%rsp)
mulsd  %xmm0,%xmm1
movsd  %xmm1,0x88(%rsp)
mov    0x38(%rsp),%r13
movslq %r13d,%rbx
lea    0x1(%rdi),%eax
mov    %eax,0x14(%rsp)
movslq %eax,%rcx
mov    %rbx,%rax
shl    $0xd,%rax
add    %rbx,%rax
mov    %rax,0x80(%rsp)
mov    %edi,%eax
sub    %r13d,%eax
lea    0x8(,%rax,8),%rax
mov    %rax,0x78(%rsp)
mov    %rcx,0x20(%rsp)
mov    %ecx,%eax
sub    %r13d,%eax
mov    %eax,0x10(%rsp)
mov    %r13,%rax
add    $0xffffffffffffffff,%rax
mov    %rax,0x60(%rsp)
mov    0x8(%rsp),%eax
sub    %edi,%eax
mov    %eax,0xc(%rsp)
not    %edi
add    0x30(%rsp),%edi
mov    %r13d,%r8d
and    $0x3,%r8d
mov    %r13d,%r15d
and    $0xfffffffc,%r15d
and    $0x3,%edi
mov    %rbx,%r14
shl    $0x10,%r14
mov    0x40(%rsp),%rax
lea    (%rax,%r14,1),%rbp
add    0x28(%rsp),%r14
xor    %r12d,%r12d
mov    %rdi,0x70(%rsp)
nopw   0x0(%rax,%rax,1)
cmpl   $0x0,0x38(%rsp)
jle    12e0 <SetUp3+0x1c0>
cmpq   $0x3,0x60(%rsp)
jae    12f0 <SetUp3+0x1d0>
xorpd  %xmm2,%xmm2
xor    %eax,%eax
jmp    1321 <SetUp3+0x201>
nopw   0x0(%rax,%rax,1)
xorpd  %xmm2,%xmm2
jmp    133e <SetUp3+0x21e>
cs nopw 0x0(%rax,%rax,1)
xorpd  %xmm2,%xmm2
xor    %eax,%eax
cs nopw 0x0(%rax,%rax,1)
addsd  -0x18(%rbp,%rax,8),%xmm2
addsd  -0x10(%rbp,%rax,8),%xmm2
addsd  -0x8(%rbp,%rax,8),%xmm2
addsd  0x0(%rbp,%rax,8),%xmm2
add    $0x4,%rax
cmp    %rax,%r15
jne    1300 <SetUp3+0x1e0>
test   %r8,%r8
je     133e <SetUp3+0x21e>
lea    (%r14,%rax,8),%rax
xor    %ecx,%ecx
nopl   0x0(%rax)
addsd  (%rax,%rcx,8),%xmm2
add    $0x1,%rcx
cmp    %rcx,%r8
jne    1330 <SetUp3+0x210>
mov    0x14(%rsp),%eax
cmp    0x30(%rsp),%eax
jge    13a3 <SetUp3+0x283>
test   %edi,%edi
je     14ce <SetUp3+0x3ae>
mov    %edi,%ecx
mov    0x20(%rsp),%rax
nopw   0x0(%rax,%rax,1)
addsd  (%r14,%rax,8),%xmm2
add    $0x1,%rax
add    $0xffffffff,%ecx
jne    1360 <SetUp3+0x240>
cmpl   $0x3,0xc(%rsp)
jb     13a3 <SetUp3+0x283>
cs nopw 0x0(%rax,%rax,1)
addsd  (%r14,%rax,8),%xmm2
addsd  0x8(%r14,%rax,8),%xmm2
addsd  0x10(%r14,%rax,8),%xmm2
addsd  0x18(%r14,%rax,8),%xmm2
add    $0x4,%rax
cmp    %eax,%esi
jne    1380 <SetUp3+0x260>
movapd %xmm2,%xmm0
mov    0xb8(%rsp),%rax
subsd  (%rax,%rbx,8),%xmm0
andpd  0xc54(%rip),%xmm0
movapd %xmm2,%xmm1
mulsd  0xc40(%rip),%xmm1
ucomisd %xmm1,%xmm0
ja     14e3 <SetUp3+0x3c3>
mov    %r12,%rax
shl    $0x10,%rax
add    0x28(%rsp),%rax
mov    0x80(%rsp),%rcx
lea    (%rax,%rcx,8),%rdi
xorpd  0xc2e(%rip),%xmm2
movapd %xmm2,0xc0(%rsp)
xor    %esi,%esi
mov    0x78(%rsp),%rdx
mov    %r8,%r13
call   1040 <memset@plt>
mov    %r13,%r8
mov    0x70(%rsp),%rdi
mov    0x68(%rsp),%rsi
movsd  0xb0(%rsp),%xmm0
movapd 0xc0(%rsp),%xmm1
mulsd  %xmm1,%xmm0
mov    0x110(%rsp),%rcx
movsd  %xmm0,(%rcx,%rbx,8)
movsd  0xa8(%rsp),%xmm0
mulsd  %xmm1,%xmm0
mov    0x118(%rsp),%rax
movsd  %xmm0,(%rax,%rbx,8)
movsd  0xa0(%rsp),%xmm0
mulsd  %xmm1,%xmm0
movsd  %xmm0,0x10000(%rcx,%rbx,8)
movsd  0x98(%rsp),%xmm0
mulsd  %xmm1,%xmm0
movsd  %xmm0,0x10000(%rax,%rbx,8)
movsd  0x90(%rsp),%xmm0
mulsd  %xmm1,%xmm0
movsd  %xmm0,0x20000(%rcx,%rbx,8)
mulsd  0x88(%rsp),%xmm1
movsd  %xmm1,0x20000(%rax,%rbx,8)
add    $0x1,%rbx
add    $0x1,%r12
add    $0x10000,%rbp
add    $0x10000,%r14
cmp    %r12d,0x10(%rsp)
jne    12c0 <SetUp3+0x1a0>
jmp    1170 <SetUp3+0x50>
mov    0x20(%rsp),%rax
cmpl   $0x3,0xc(%rsp)
jae    1380 <SetUp3+0x260>
jmp    13a3 <SetUp3+0x283>
lea    0xb46(%rip),%rdi
call   1030 <puts@plt>
xor    %eax,%eax
add    $0xd8,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret
mov    $0x1,%eax
jmp    14f1 <SetUp3+0x3d1>

LLVM IR (-O0):
@.str = [42 x i8] c"Total form factor is too far from unity.\0A\00"
define i32 @SetUp3(i32 %0, [2 x i32]* %1, double* %2, [3 x double]* %3, [3 x double]* %4, [8192 x double]* %5, [8192 x double]* %6, [8192 x double]* %7) {
%9 = alloca i32
%10 = alloca i32
%11 = alloca [2 x i32]*
%12 = alloca double*
%13 = alloca [3 x double]*
%14 = alloca [3 x double]*
%15 = alloca [8192 x double]*
%16 = alloca [8192 x double]*
%17 = alloca [8192 x double]*
%18 = alloca i32
%19 = alloca i32
%20 = alloca i32
%21 = alloca i32
%22 = alloca double
%23 = alloca [3 x double]
%24 = alloca [3 x double]
store i32 %0, i32* %10
store [2 x i32]* %1, [2 x i32]** %11
store double* %2, double** %12
store [3 x double]* %3, [3 x double]** %13
store [3 x double]* %4, [3 x double]** %14
store [8192 x double]* %5, [8192 x double]** %15
store [8192 x double]* %6, [8192 x double]** %16
store [8192 x double]* %7, [8192 x double]** %17
store i32 0, i32* %19
br label %25
25: ; preds = %212, %8
%26 = load i32, i32* %19
%27 = icmp slt i32 %26, 6
br i1 %27, label %28, label %215
28: ; preds = %25
store i32 0, i32* %18
br label %29
29: ; preds = %61, %28
%30 = load i32, i32* %18
%31 = icmp slt i32 %30, 3
br i1 %31, label %32, label %64
32: ; preds = %29
%33 = load [3 x double]*, [3 x double]** %13
%34 = load i32, i32* %19
%35 = sext i32 %34 to i64
%36 = getelementptr inbounds [3 x double], [3 x double]* %33, i64 %35
%37 = load i32, i32* %18
%38 = sext i32 %37 to i64
%39 = getelementptr inbounds [3 x double], [3 x double]* %36, i64 0, i64 %38
%40 = load double, double* %39
%41 = fdiv double 1.000000e+00, %40
%42 = load i32, i32* %18
%43 = sext i32 %42 to i64
%44 = getelementptr inbounds [3 x double], [3 x double]* %23, i64 0, i64 %43
store double %41, double* %44
%45 = load [3 x double]*, [3 x double]** %14
%46 = load i32, i32* %19
%47 = sext i32 %46 to i64
%48 = getelementptr inbounds [3 x double], [3 x double]* %45, i64 %47
%49 = load i32, i32* %18
%50 = sext i32 %49 to i64
%51 = getelementptr inbounds [3 x double], [3 x double]* %48, i64 0, i64 %50
%52 = load double, double* %51
%53 = load i32, i32* %18
%54 = sext i32 %53 to i64
%55 = getelementptr inbounds [3 x double], [3 x double]* %23, i64 0, i64 %54
%56 = load double, double* %55
%57 = fmul double %52, %56
%58 = load i32, i32* %18
%59 = sext i32 %58 to i64
%60 = getelementptr inbounds [3 x double], [3 x double]* %24, i64 0, i64 %59
store double %57, double* %60
br label %61
61: ; preds = %32
%62 = load i32, i32* %18
%63 = add nsw i32 %62, 1
store i32 %63, i32* %18
br label %29
64: ; preds = %29
%65 = load [2 x i32]*, [2 x i32]** %11
%66 = load i32, i32* %19
%67 = sext i32 %66 to i64
%68 = getelementptr inbounds [2 x i32], [2 x i32]* %65, i64 %67
%69 = getelementptr inbounds [2 x i32], [2 x i32]* %68, i64 0, i64 0
%70 = load i32, i32* %69
store i32 %70, i32* %20
br label %71
71: ; preds = %208, %64
%72 = load i32, i32* %20
%73 = load [2 x i32]*, [2 x i32]** %11
%74 = load i32, i32* %19
%75 = sext i32 %74 to i64
%76 = getelementptr inbounds [2 x i32], [2 x i32]* %73, i64 %75
%77 = getelementptr inbounds [2 x i32], [2 x i32]* %76, i64 0, i64 1
%78 = load i32, i32* %77
%79 = icmp sle i32 %72, %78
br i1 %79, label %80, label %211
80: ; preds = %71
store double 0.000000e+00, double* %22
store i32 0, i32* %21
br label %81
81: ; preds = %101, %80
%82 = load i32, i32* %21
%83 = load [2 x i32]*, [2 x i32]** %11
%84 = load i32, i32* %19
%85 = sext i32 %84 to i64
%86 = getelementptr inbounds [2 x i32], [2 x i32]* %83, i64 %85
%87 = getelementptr inbounds [2 x i32], [2 x i32]* %86, i64 0, i64 0
%88 = load i32, i32* %87
%89 = icmp slt i32 %82, %88
br i1 %89, label %90, label %104
90: ; preds = %81
%91 = load [8192 x double]*, [8192 x double]** %15
%92 = load i32, i32* %20
%93 = sext i32 %92 to i64
%94 = getelementptr inbounds [8192 x double], [8192 x double]* %91, i64 %93
%95 = load i32, i32* %21
%96 = sext i32 %95 to i64
%97 = getelementptr inbounds [8192 x double], [8192 x double]* %94, i64 0, i64 %96
%98 = load double, double* %97
%99 = load double, double* %22
%100 = fadd double %99, %98
store double %100, double* %22
br label %101
101: ; preds = %90
%102 = load i32, i32* %21
%103 = add nsw i32 %102, 1
store i32 %103, i32* %21
br label %81
104: ; preds = %81
%105 = load [2 x i32]*, [2 x i32]** %11
%106 = load i32, i32* %19
%107 = sext i32 %106 to i64
%108 = getelementptr inbounds [2 x i32], [2 x i32]* %105, i64 %107
%109 = getelementptr inbounds [2 x i32], [2 x i32]* %108, i64 0, i64 1
%110 = load i32, i32* %109
%111 = add nsw i32 %110, 1
store i32 %111, i32* %21
br label %112
112: ; preds = %127, %104
%113 = load i32, i32* %21
%114 = load i32, i32* %10
%115 = icmp slt i32 %113, %114
br i1 %115, label %116, label %130
116: ; preds = %112
%117 = load [8192 x double]*, [8192 x double]** %15
%118 = load i32, i32* %20
%119 = sext i32 %118 to i64
%120 = getelementptr inbounds [8192 x double], [8192 x double]* %117, i64 %119
%121 = load i32, i32* %21
%122 = sext i32 %121 to i64
%123 = getelementptr inbounds [8192 x double], [8192 x double]* %120, i64 0, i64 %122
%124 = load double, double* %123
%125 = load double, double* %22
%126 = fadd double %125, %124
store double %126, double* %22
br label %127
127: ; preds = %116
%128 = load i32, i32* %21
%129 = add nsw i32 %128, 1
store i32 %129, i32* %21
br label %112
130: ; preds = %112
%131 = load double, double* %22
%132 = load double*, double** %12
%133 = load i32, i32* %20
%134 = sext i32 %133 to i64
%135 = getelementptr inbounds double, double* %132, i64 %134
%136 = load double, double* %135
%137 = fsub double %131, %136
%138 = call double @llvm.fabs.f64(double %137)
%139 = load double, double* %22
%140 = fmul double 5.000000e-10, %139
%141 = fcmp ogt double %138, %140
br i1 %141, label %142, label %144
142: ; preds = %130
%143 = call i32 (i8*, ...) @printf(i8* getelementptr inbounds ([42 x i8], [42 x i8]* @.str, i64 0, i64 0))
store i32 0, i32* %9
br label %216
144: ; preds = %130
%145 = load double, double* %22
%146 = fneg double %145
store double %146, double* %22
%147 = load [2 x i32]*, [2 x i32]** %11
%148 = load i32, i32* %19
%149 = sext i32 %148 to i64
%150 = getelementptr inbounds [2 x i32], [2 x i32]* %147, i64 %149
%151 = getelementptr inbounds [2 x i32], [2 x i32]* %150, i64 0, i64 0
%152 = load i32, i32* %151
store i32 %152, i32* %21
br label %153
153: ; preds = %170, %144
%154 = load i32, i32* %21
%155 = load [2 x i32]*, [2 x i32]** %11
%156 = load i32, i32* %19
%157 = sext i32 %156 to i64
%158 = getelementptr inbounds [2 x i32], [2 x i32]* %155, i64 %157
%159 = getelementptr inbounds [2 x i32], [2 x i32]* %158, i64 0, i64 1
%160 = load i32, i32* %159
%161 = icmp sle i32 %154, %160
br i1 %161, label %162, label %173
162: ; preds = %153
%163 = load [8192 x double]*, [8192 x double]** %15
%164 = load i32, i32* %20
%165 = sext i32 %164 to i64
%166 = getelementptr inbounds [8192 x double], [8192 x double]* %163, i64 %165
%167 = load i32, i32* %21
%168 = sext i32 %167 to i64
%169 = getelementptr inbounds [8192 x double], [8192 x double]* %166, i64 0, i64 %168
store double 0.000000e+00, double* %169
br label %170
170: ; preds = %162
%171 = load i32, i32* %21
%172 = add nsw i32 %171, 1
store i32 %172, i32* %21
br label %153
173: ; preds = %153
store i32 0, i32* %18
br label %174
174: ; preds = %204, %173
%175 = load i32, i32* %18
%176 = icmp slt i32 %175, 3
br i1 %176, label %177, label %207
177: ; preds = %174
%178 = load i32, i32* %18
%179 = sext i32 %178 to i64
%180 = getelementptr inbounds [3 x double], [3 x double]* %23, i64 0, i64 %179
%181 = load double, double* %180
%182 = load double, double* %22
%183 = fmul double %181, %182
%184 = load [8192 x double]*, [8192 x double]** %16
%185 = load i32, i32* %18
%186 = sext i32 %185 to i64
%187 = getelementptr inbounds [8192 x double], [8192 x double]* %184, i64 %186
%188 = load i32, i32* %20
%189 = sext i32 %188 to i64
%190 = getelementptr inbounds [8192 x double], [8192 x double]* %187, i64 0, i64 %189
store double %183, double* %190
%191 = load i32, i32* %18
%192 = sext i32 %191 to i64
%193 = getelementptr inbounds [3 x double], [3 x double]* %24, i64 0, i64 %192
%194 = load double, double* %193
%195 = load double, double* %22
%196 = fmul double %194, %195
%197 = load [8192 x double]*, [8192 x double]** %17
%198 = load i32, i32* %18
%199 = sext i32 %198 to i64
%200 = getelementptr inbounds [8192 x double], [8192 x double]* %197, i64 %199
%201 = load i32, i32* %20
%202 = sext i32 %201 to i64
%203 = getelementptr inbounds [8192 x double], [8192 x double]* %200, i64 0, i64 %202
store double %196, double* %203
br label %204
204: ; preds = %177
%205 = load i32, i32* %18
%206 = add nsw i32 %205, 1
store i32 %206, i32* %18
br label %174
207: ; preds = %174
br label %208
208: ; preds = %207
%209 = load i32, i32* %20
%210 = add nsw i32 %209, 1
store i32 %210, i32* %20
br label %71
211: ; preds = %71
br label %212
212: ; preds = %211
%213 = load i32, i32* %19
%214 = add nsw i32 %213, 1
store i32 %214, i32* %19
br label %25
215: ; preds = %25
store i32 1, i32* %9
br label %216
216: ; preds = %215, %142
%217 = load i32, i32* %9
ret i32 %217
}

LLVM IR (-O2):
@str = [41 x i8] c"Total form factor is too far from unity.\00"
define i32 @SetUp3(i32 %0, [2 x i32]* nocapture %1, double* nocapture %2, [3 x double]* nocapture %3, [3 x double]* nocapture %4, [8192 x double]* nocapture %5, [8192 x double]* nocapture writeonly %6, [8192 x double]* nocapture writeonly %7) {
%9 = add i32 %0, -2
br label %10
10: ; preds = %8, %167
%11 = phi i64 [ 0, %8 ], [ %168, %167 ]
%12 = getelementptr inbounds [3 x double], [3 x double]* %3, i64 %11, i64 0
%13 = load double, double* %12
%14 = fdiv double 1.000000e+00, %13
%15 = getelementptr inbounds [3 x double], [3 x double]* %4, i64 %11, i64 0
%16 = load double, double* %15
%17 = fmul double %14, %16
%18 = getelementptr inbounds [3 x double], [3 x double]* %3, i64 %11, i64 1
%19 = load double, double* %18
%20 = fdiv double 1.000000e+00, %19
%21 = getelementptr inbounds [3 x double], [3 x double]* %4, i64 %11, i64 1
%22 = load double, double* %21
%23 = fmul double %20, %22
%24 = getelementptr inbounds [3 x double], [3 x double]* %3, i64 %11, i64 2
%25 = load double, double* %24
%26 = fdiv double 1.000000e+00, %25
%27 = getelementptr inbounds [3 x double], [3 x double]* %4, i64 %11, i64 2
%28 = load double, double* %27
%29 = fmul double %26, %28
%30 = getelementptr inbounds [2 x i32], [2 x i32]* %1, i64 %11, i64 0
%31 = load i32, i32* %30
%32 = getelementptr inbounds [2 x i32], [2 x i32]* %1, i64 %11, i64 1
%33 = load i32, i32* %32
%34 = icmp slt i32 %33, %31
br i1 %34, label %167, label %35
35: ; preds = %10
%36 = icmp sgt i32 %31, 0
%37 = add i32 %33, 1
%38 = icmp slt i32 %37, %0
%39 = sext i32 %37 to i64
%40 = sext i32 %31 to i64
%41 = mul nsw i64 %40, 8193
%42 = sub i32 %33, %31
%43 = zext i32 %42 to i64
%44 = shl nuw nsw i64 %43, 3
%45 = add nuw nsw i64 %44, 8
%46 = sub i32 %37, %31
%47 = zext i32 %31 to i64
%48 = add nsw i64 %47, -1
%49 = xor i32 %33, -1
%50 = add i32 %49, %0
%51 = sub i32 %9, %33
%52 = and i64 %47, 3
%53 = icmp ult i64 %48, 3
%54 = and i64 %47, 4294967292
%55 = icmp eq i64 %52, 0
%56 = and i32 %50, 3
%57 = icmp eq i32 %56, 0
%58 = icmp ult i32 %51, 3
br label %59
59: ; preds = %35, %149
%60 = phi i64 [ %40, %35 ], [ %163, %149 ]
%61 = phi i64 [ 0, %35 ], [ %164, %149 ]
%62 = getelementptr [8192 x double], [8192 x double]* %5, i64 %61, i64 %41
%63 = bitcast double* %62 to i8*
br i1 %36, label %64, label %79
64: ; preds = %59
br i1 %53, label %65, label %96
65: ; preds = %96, %64
%66 = phi double [ undef, %64 ], [ %114, %96 ]
%67 = phi i64 [ 0, %64 ], [ %115, %96 ]
%68 = phi double [ 0.000000e+00, %64 ], [ %114, %96 ]
br i1 %55, label %79, label %69
69: ; preds = %65, %69
%70 = phi i64 [ %76, %69 ], [ %67, %65 ]
%71 = phi double [ %75, %69 ], [ %68, %65 ]
%72 = phi i64 [ %77, %69 ], [ 0, %65 ]
%73 = getelementptr inbounds [8192 x double], [8192 x double]* %5, i64 %60, i64 %70
%74 = load double, double* %73
%75 = fadd double %71, %74
%76 = add nuw nsw i64 %70, 1
%77 = add i64 %72, 1
%78 = icmp eq i64 %77, %52
br i1 %78, label %79, label %69
79: ; preds = %65, %69, %59
%80 = phi double [ 0.000000e+00, %59 ], [ %66, %65 ], [ %75, %69 ]
br i1 %38, label %81, label %139
81: ; preds = %79
br i1 %57, label %92, label %82
82: ; preds = %81, %82
%83 = phi i64 [ %89, %82 ], [ %39, %81 ]
%84 = phi double [ %88, %82 ], [ %80, %81 ]
%85 = phi i32 [ %90, %82 ], [ 0, %81 ]
%86 = getelementptr inbounds [8192 x double], [8192 x double]* %5, i64 %60, i64 %83
%87 = load double, double* %86
%88 = fadd double %84, %87
%89 = add nsw i64 %83, 1
%90 = add i32 %85, 1
%91 = icmp eq i32 %90, %56
br i1 %91, label %92, label %82
92: ; preds = %82, %81
%93 = phi double [ undef, %81 ], [ %88, %82 ]
%94 = phi i64 [ %39, %81 ], [ %89, %82 ]
%95 = phi double [ %80, %81 ], [ %88, %82 ]
br i1 %58, label %139, label %118
96: ; preds = %64, %96
%97 = phi i64 [ %115, %96 ], [ 0, %64 ]
%98 = phi double [ %114, %96 ], [ 0.000000e+00, %64 ]
%99 = phi i64 [ %116, %96 ], [ 0, %64 ]
%100 = getelementptr inbounds [8192 x double], [8192 x double]* %5, i64 %60, i64 %97
%101 = load double, double* %100
%102 = fadd double %98, %101
%103 = or i64 %97, 1
%104 = getelementptr inbounds [8192 x double], [8192 x double]* %5, i64 %60, i64 %103
%105 = load double, double* %104
%106 = fadd double %102, %105
%107 = or i64 %97, 2
%108 = getelementptr inbounds [8192 x double], [8192 x double]* %5, i64 %60, i64 %107
%109 = load double, double* %108
%110 = fadd double %106, %109
%111 = or i64 %97, 3
%112 = getelementptr inbounds [8192 x double], [8192 x double]* %5, i64 %60, i64 %111
%113 = load double, double* %112
%114 = fadd double %110, %113
%115 = add nuw nsw i64 %97, 4
%116 = add i64 %99, 4
%117 = icmp eq i64 %116, %54
br i1 %117, label %65, label %96
118: ; preds = %92, %118
%119 = phi i64 [ %136, %118 ], [ %94, %92 ]
%120 = phi double [ %135, %118 ], [ %95, %92 ]
%121 = getelementptr inbounds [8192 x double], [8192 x double]* %5, i64 %60, i64 %119
%122 = load double, double* %121
%123 = fadd double %120, %122
%124 = add nsw i64 %119, 1
%125 = getelementptr inbounds [8192 x double], [8192 x double]* %5, i64 %60, i64 %124
%126 = load double, double* %125
%127 = fadd double %123, %126
%128 = add nsw i64 %119, 2
%129 = getelementptr inbounds [8192 x double], [8192 x double]* %5, i64 %60, i64 %128
%130 = load double, double* %129
%131 = fadd double %127, %130
%132 = add nsw i64 %119, 3
%133 = getelementptr inbounds [8192 x double], [8192 x double]* %5, i64 %60, i64 %132
%134 = load double, double* %133
%135 = fadd double %131, %134
%136 = add nsw i64 %119, 4
%137 = trunc i64 %136 to i32
%138 = icmp eq i32 %137, %0
br i1 %138, label %139, label %118
139: ; preds = %92, %118, %79
%140 = phi double [ %80, %79 ], [ %93, %92 ], [ %135, %118 ]
%141 = getelementptr inbounds double, double* %2, i64 %60
%142 = load double, double* %141
%143 = fsub double %140, %142
%144 = tail call double @llvm.fabs.f64(double %143)
%145 = fmul double %140, 5.000000e-10
%146 = fcmp ogt double %144, %145
br i1 %146, label %147, label %149
147: ; preds = %139
%148 = tail call i32 @puts(i8* dereferenceable(1) getelementptr inbounds ([41 x i8], [41 x i8]* @str, i64 0, i64 0))
br label %170
149: ; preds = %139
%150 = fneg double %140
%151 = fmul double %14, %150
%152 = getelementptr inbounds [8192 x double], [8192 x double]* %6, i64 0, i64 %60
store double %151, double* %152
%153 = fmul double %17, %150
%154 = getelementptr inbounds [8192 x double], [8192 x double]* %7, i64 0, i64 %60
store double %153, double* %154
%155 = fmul double %20, %150
%156 = getelementptr inbounds [8192 x double], [8192 x double]* %6, i64 1, i64 %60
store double %155, double* %156
%157 = fmul double %23, %150
%158 = getelementptr inbounds [8192 x double], [8192 x double]* %7, i64 1, i64 %60
store double %157, double* %158
%159 = fmul double %26, %150
%160 = getelementptr inbounds [8192 x double], [8192 x double]* %6, i64 2, i64 %60
store double %159, double* %160
%161 = fmul double %29, %150
%162 = getelementptr inbounds [8192 x double], [8192 x double]* %7, i64 2, i64 %60
store double %161, double* %162
%163 = add nsw i64 %60, 1
%164 = add nuw nsw i64 %61, 1
%165 = trunc i64 %164 to i32
%166 = icmp eq i32 %46, %165
br i1 %166, label %167, label %59
167: ; preds = %149, %10
%168 = add nuw nsw i64 %11, 1
%169 = icmp eq i64 %168, 6
br i1 %169, label %170, label %10
170: ; preds = %167, %147
%171 = phi i32 [ 0, %147 ], [ 1, %167 ]
ret i32 %171
}