int PolygonD_Clip(int innumpoints, const double *inpoints, double planenormalx,
                  double planenormaly, double planenormalz, double planedist,
                  double epsilon, int outfrontmaxpoints,
                  double *outfrontpoints) {
  int i, frontcount = 0;
  const double *n, *p;
  double frac, pdist, ndist;
  if (innumpoints < 1)
    return 0;
  n = inpoints;
  ndist = n[0] * planenormalx + n[1] * planenormaly + n[2] * planenormalz -
          planedist;
  for (i = 0; i < innumpoints; i++) {
    p = n;
    pdist = ndist;
    n = inpoints + ((i + 1) < innumpoints ? (i + 1) : 0) * 3;
    ndist = n[0] * planenormalx + n[1] * planenormaly + n[2] * planenormalz -
            planedist;
    if (pdist >= -epsilon) {
      if (frontcount < outfrontmaxpoints) {
        *outfrontpoints++ = p[0];
        *outfrontpoints++ = p[1];
        *outfrontpoints++ = p[2];
      }
      frontcount++;
    }
    if ((pdist > epsilon && ndist < -epsilon) ||
        (pdist < -epsilon && ndist > epsilon)) {
      frac = pdist / (pdist - ndist);
      if (frontcount < outfrontmaxpoints) {
        *outfrontpoints++ = p[0] + frac * (n[0] - p[0]);
        *outfrontpoints++ = p[1] + frac * (n[1] - p[1]);
        *outfrontpoints++ = p[2] + frac * (n[2] - p[2]);
      }
      frontcount++;
    }
  }
  return frontcount;
}


x86 assembly (clang, -O0):
<PolygonD_Clip>:
push   %rbp
mov    %rsp,%rbp
push   %rax
mov    %edi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
movsd  %xmm0,-0x18(%rbp)
movsd  %xmm1,-0x20(%rbp)
movsd  %xmm2,-0x28(%rbp)
movsd  %xmm3,-0x30(%rbp)
movsd  %xmm4,-0x38(%rbp)
mov    %edx,-0x3c(%rbp)
mov    %rcx,-0x48(%rbp)
movl   $0x0,-0x50(%rbp)
cmpl   $0x1,-0x8(%rbp)
jge    1149 <PolygonD_Clip+0x49>
movl   $0x0,-0x4(%rbp)
jmp    1449 <PolygonD_Clip+0x349>
mov    -0x10(%rbp),%rax
mov    %rax,-0x58(%rbp)
mov    -0x58(%rbp),%rax
movsd  (%rax),%xmm1
movsd  -0x18(%rbp),%xmm2
mov    -0x58(%rbp),%rax
movsd  0x8(%rax),%xmm0
mulsd  -0x20(%rbp),%xmm0
mulsd  %xmm2,%xmm1
addsd  %xmm0,%xmm1
mov    -0x58(%rbp),%rax
movsd  0x10(%rax),%xmm0
movsd  -0x28(%rbp),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
subsd  -0x30(%rbp),%xmm0
movsd  %xmm0,-0x78(%rbp)
movl   $0x0,-0x4c(%rbp)
mov    -0x4c(%rbp),%eax
cmp    -0x8(%rbp),%eax
jge    1443 <PolygonD_Clip+0x343>
mov    -0x58(%rbp),%rax
mov    %rax,-0x60(%rbp)
movsd  -0x78(%rbp),%xmm0
movsd  %xmm0,-0x70(%rbp)
mov    -0x10(%rbp),%rax
mov    %rax,-0x80(%rbp)
mov    -0x4c(%rbp),%eax
add    $0x1,%eax
cmp    -0x8(%rbp),%eax
jge    11e1 <PolygonD_Clip+0xe1>
mov    -0x4c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x84(%rbp)
jmp    11ee <PolygonD_Clip+0xee>
xor    %eax,%eax
mov    %eax,-0x84(%rbp)
jmp    11ee <PolygonD_Clip+0xee>
mov    -0x80(%rbp),%rax
mov    -0x84(%rbp),%ecx
imul   $0x3,%ecx,%ecx
movslq %ecx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    %rax,-0x58(%rbp)
mov    -0x58(%rbp),%rax
movsd  (%rax),%xmm1
movsd  -0x18(%rbp),%xmm2
mov    -0x58(%rbp),%rax
movsd  0x8(%rax),%xmm0
mulsd  -0x20(%rbp),%xmm0
mulsd  %xmm2,%xmm1
addsd  %xmm0,%xmm1
mov    -0x58(%rbp),%rax
movsd  0x10(%rax),%xmm0
movsd  -0x28(%rbp),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
subsd  -0x30(%rbp),%xmm0
movsd  %xmm0,-0x78(%rbp)
movsd  -0x70(%rbp),%xmm0
movsd  -0x38(%rbp),%xmm1
movq   %xmm1,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm1
ucomisd %xmm1,%xmm0
jb     12df <PolygonD_Clip+0x1df>
mov    -0x50(%rbp),%eax
cmp    -0x3c(%rbp),%eax
jge    12d6 <PolygonD_Clip+0x1d6>
mov    -0x60(%rbp),%rax
movsd  (%rax),%xmm0
mov    -0x48(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x48(%rbp)
movsd  %xmm0,(%rax)
mov    -0x60(%rbp),%rax
movsd  0x8(%rax),%xmm0
mov    -0x48(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x48(%rbp)
movsd  %xmm0,(%rax)
mov    -0x60(%rbp),%rax
movsd  0x10(%rax),%xmm0
mov    -0x48(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x48(%rbp)
movsd  %xmm0,(%rax)
mov    -0x50(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x50(%rbp)
movsd  -0x70(%rbp),%xmm0
ucomisd -0x38(%rbp),%xmm0
jbe    131a <PolygonD_Clip+0x21a>
movsd  -0x78(%rbp),%xmm1
movsd  -0x38(%rbp),%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
ucomisd %xmm1,%xmm0
ja     1355 <PolygonD_Clip+0x255>
movsd  -0x70(%rbp),%xmm1
movsd  -0x38(%rbp),%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
ucomisd %xmm1,%xmm0
jbe    1430 <PolygonD_Clip+0x330>
movsd  -0x78(%rbp),%xmm0
ucomisd -0x38(%rbp),%xmm0
jbe    1430 <PolygonD_Clip+0x330>
movsd  -0x70(%rbp),%xmm0
movsd  -0x70(%rbp),%xmm1
subsd  -0x78(%rbp),%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x68(%rbp)
mov    -0x50(%rbp),%eax
cmp    -0x3c(%rbp),%eax
jge    1427 <PolygonD_Clip+0x327>
mov    -0x60(%rbp),%rax
movsd  (%rax),%xmm1
movsd  -0x68(%rbp),%xmm0
mov    -0x58(%rbp),%rax
movsd  (%rax),%xmm2
mov    -0x60(%rbp),%rax
subsd  (%rax),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x48(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x48(%rbp)
movsd  %xmm0,(%rax)
mov    -0x60(%rbp),%rax
movsd  0x8(%rax),%xmm1
movsd  -0x68(%rbp),%xmm0
mov    -0x58(%rbp),%rax
movsd  0x8(%rax),%xmm2
mov    -0x60(%rbp),%rax
subsd  0x8(%rax),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x48(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x48(%rbp)
movsd  %xmm0,(%rax)
mov    -0x60(%rbp),%rax
movsd  0x10(%rax),%xmm1
movsd  -0x68(%rbp),%xmm0
mov    -0x58(%rbp),%rax
movsd  0x10(%rax),%xmm2
mov    -0x60(%rbp),%rax
subsd  0x10(%rax),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x48(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x48(%rbp)
movsd  %xmm0,(%rax)
mov    -0x50(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x50(%rbp)
jmp    1435 <PolygonD_Clip+0x335>
mov    -0x4c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x4c(%rbp)
jmp    119b <PolygonD_Clip+0x9b>
mov    -0x50(%rbp),%eax
mov    %eax,-0x4(%rbp)
mov    -0x4(%rbp),%eax
add    $0x8,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<PolygonD_Clip>:
push   %rbp
push   %r14
push   %rbx
test   %edi,%edi
jle    1291 <PolygonD_Clip+0x191>
movsd  0x8(%rsi),%xmm5
mulsd  %xmm1,%xmm5
movsd  (%rsi),%xmm7
mulsd  %xmm0,%xmm7
addsd  %xmm5,%xmm7
movsd  0x10(%rsi),%xmm10
mulsd  %xmm2,%xmm10
addsd  %xmm7,%xmm10
subsd  %xmm3,%xmm10
movapd 0xec1(%rip),%xmm8
xorpd  %xmm4,%xmm8
mov    %edi,%r8d
neg    %r8d
xor    %eax,%eax
mov    $0x1,%r10d
mov    $0x3,%r9d
mov    %rsi,%r11
jmp    1180 <PolygonD_Clip+0x80>
nopl   (%rax)
add    $0x1,%eax
lea    (%r8,%r10,1),%ebp
add    $0x1,%ebp
mov    %r10d,%ebx
add    $0x1,%ebx
add    $0x3,%r9d
mov    %ebx,%r10d
cmp    $0x1,%ebp
je     1293 <PolygonD_Clip+0x193>
mov    %r11,%r14
movapd %xmm10,%xmm9
cmp    %edi,%r10d
mov    $0x0,%ebx
cmovl  %r9d,%ebx
movsd  0x8(%rsi,%rbx,8),%xmm5
mulsd  %xmm1,%xmm5
movsd  (%rsi,%rbx,8),%xmm7
mulsd  %xmm0,%xmm7
movsd  0x10(%rsi,%rbx,8),%xmm10
mulsd  %xmm2,%xmm10
ucomisd %xmm8,%xmm9
jb     11e4 <PolygonD_Clip+0xe4>
cmp    %edx,%eax
jge    11e1 <PolygonD_Clip+0xe1>
movsd  (%r14),%xmm6
movsd  %xmm6,(%rcx)
movsd  0x8(%r14),%xmm6
movsd  %xmm6,0x8(%rcx)
movsd  0x10(%r14),%xmm6
movsd  %xmm6,0x10(%rcx)
add    $0x18,%rcx
add    $0x1,%eax
lea    (%rsi,%rbx,8),%r11
addsd  %xmm5,%xmm7
addsd  %xmm7,%xmm10
subsd  %xmm3,%xmm10
ucomisd %xmm4,%xmm9
jbe    1204 <PolygonD_Clip+0x104>
ucomisd %xmm10,%xmm8
ja     121a <PolygonD_Clip+0x11a>
ucomisd %xmm9,%xmm8
jbe    1163 <PolygonD_Clip+0x63>
ucomisd %xmm4,%xmm10
jbe    1163 <PolygonD_Clip+0x63>
cmp    %edx,%eax
jge    1160 <PolygonD_Clip+0x60>
movapd %xmm9,%xmm5
subsd  %xmm10,%xmm5
divsd  %xmm5,%xmm9
movsd  (%r14),%xmm5
movsd  (%r11),%xmm6
subsd  %xmm5,%xmm6
mulsd  %xmm9,%xmm6
addsd  %xmm5,%xmm6
movsd  %xmm6,(%rcx)
movsd  0x8(%r14),%xmm5
movsd  0x8(%r11),%xmm6
subsd  %xmm5,%xmm6
mulsd  %xmm9,%xmm6
addsd  %xmm5,%xmm6
movsd  %xmm6,0x8(%rcx)
movsd  0x10(%r14),%xmm5
movsd  0x10(%r11),%xmm6
subsd  %xmm5,%xmm6
mulsd  %xmm9,%xmm6
addsd  %xmm5,%xmm6
movsd  %xmm6,0x10(%rcx)
add    $0x18,%rcx
jmp    1160 <PolygonD_Clip+0x60>
xor    %eax,%eax
pop    %rbx
pop    %r14
pop    %rbp
ret

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

LLVM IR (-O2):
define i32 @PolygonD_Clip(i32 %0, double* nocapture %1, double %2, double %3, double %4, double %5, double %6, i32 %7, double* nocapture writeonly %8) {
%10 = icmp slt i32 %0, 1
br i1 %10, label %97, label %11
11: ; preds = %9
%12 = getelementptr inbounds double, double* %1, i64 2
%13 = load double, double* %12
%14 = load double, double* %1
%15 = getelementptr inbounds double, double* %1, i64 1
%16 = load double, double* %15
%17 = fmul double %16, %3
%18 = tail call double @llvm.fmuladd.f64(double %14, double %2, double %17)
%19 = tail call double @llvm.fmuladd.f64(double %13, double %4, double %18)
%20 = fsub double %19, %5
%21 = fneg double %6
br label %22
22: ; preds = %11, %93
%23 = phi double [ %20, %11 ], [ %42, %93 ]
%24 = phi double* [ %8, %11 ], [ %95, %93 ]
%25 = phi i32 [ 0, %11 ], [ %28, %93 ]
%26 = phi double* [ %1, %11 ], [ %33, %93 ]
%27 = phi i32 [ 0, %11 ], [ %94, %93 ]
%28 = add nuw nsw i32 %25, 1
%29 = icmp slt i32 %28, %0
%30 = mul nsw i32 %28, 3
%31 = select i1 %29, i32 %30, i32 0
%32 = zext i32 %31 to i64
%33 = getelementptr inbounds double, double* %1, i64 %32
%34 = load double, double* %33
%35 = getelementptr inbounds double, double* %33, i64 1
%36 = load double, double* %35
%37 = fmul double %36, %3
%38 = tail call double @llvm.fmuladd.f64(double %34, double %2, double %37)
%39 = getelementptr inbounds double, double* %33, i64 2
%40 = load double, double* %39
%41 = tail call double @llvm.fmuladd.f64(double %40, double %4, double %38)
%42 = fsub double %41, %5
%43 = fcmp ult double %23, %21
br i1 %43, label %58, label %44
44: ; preds = %22
%45 = icmp slt i32 %27, %7
br i1 %45, label %46, label %55
46: ; preds = %44
%47 = load double, double* %26
%48 = getelementptr inbounds double, double* %24, i64 1
store double %47, double* %24
%49 = getelementptr inbounds double, double* %26, i64 1
%50 = load double, double* %49
%51 = getelementptr inbounds double, double* %24, i64 2
store double %50, double* %48
%52 = getelementptr inbounds double, double* %26, i64 2
%53 = load double, double* %52
%54 = getelementptr inbounds double, double* %24, i64 3
store double %53, double* %51
br label %55
55: ; preds = %46, %44
%56 = phi double* [ %54, %46 ], [ %24, %44 ]
%57 = add nsw i32 %27, 1
br label %58
58: ; preds = %55, %22
%59 = phi i32 [ %57, %55 ], [ %27, %22 ]
%60 = phi double* [ %56, %55 ], [ %24, %22 ]
%61 = fcmp ogt double %23, %6
%62 = fcmp olt double %42, %21
%63 = select i1 %61, i1 %62, i1 false
br i1 %63, label %68, label %64
64: ; preds = %58
%65 = fcmp olt double %23, %21
%66 = fcmp ogt double %42, %6
%67 = select i1 %65, i1 %66, i1 false
br i1 %67, label %68, label %93
68: ; preds = %64, %58
%69 = fsub double %23, %42
%70 = fdiv double %23, %69
%71 = icmp slt i32 %59, %7
br i1 %71, label %72, label %90
72: ; preds = %68
%73 = load double, double* %26
%74 = load double, double* %33
%75 = fsub double %74, %73
%76 = tail call double @llvm.fmuladd.f64(double %70, double %75, double %73)
%77 = getelementptr inbounds double, double* %60, i64 1
store double %76, double* %60
%78 = getelementptr inbounds double, double* %26, i64 1
%79 = load double, double* %78
%80 = load double, double* %35
%81 = fsub double %80, %79
%82 = tail call double @llvm.fmuladd.f64(double %70, double %81, double %79)
%83 = getelementptr inbounds double, double* %60, i64 2
store double %82, double* %77
%84 = getelementptr inbounds double, double* %26, i64 2
%85 = load double, double* %84
%86 = load double, double* %39
%87 = fsub double %86, %85
%88 = tail call double @llvm.fmuladd.f64(double %70, double %87, double %85)
%89 = getelementptr inbounds double, double* %60, i64 3
store double %88, double* %83
br label %90
90: ; preds = %72, %68
%91 = phi double* [ %89, %72 ], [ %60, %68 ]
%92 = add nsw i32 %59, 1
br label %93
93: ; preds = %64, %90
%94 = phi i32 [ %92, %90 ], [ %59, %64 ]
%95 = phi double* [ %91, %90 ], [ %60, %64 ]
%96 = icmp eq i32 %28, %0
br i1 %96, label %97, label %22
97: ; preds = %93, %9
%98 = phi i32 [ 0, %9 ], [ %94, %93 ]
ret i32 %98
}