float fd_sweep(grid, field, nx, ny, nz, dx, dy, dz) float (*grid)[];
float (*field)[];
int nx, ny, nz;
float dx, dy, dz;
{
  int i, j, k;
  float delta, deltamax;
  float flux;
  float flx, fly, flz, dv;
  int izp, iyp, ixp, inxy;
  flx = dy * dz / dx;
  fly = dx * dz / dy;
  flz = dx * dy / dz;
  dv = 2. * (dx * dy + dy * dz + dx * dz);
  dv = 1. / dv;
  dv = dv / (4. * 3.1415926535897932384626433);
  inxy = nx * ny;
  deltamax = 0.;
  for (k = 1; k < nz - 1; k++) {
    izp = k * inxy;
    for (j = 1; j < ny - 1; j++) {
      iyp = j * nx;
      for (i = 1; i < nx - 1; i++) {
        ixp = i + iyp + izp;
        flux = ((*field)[ixp] - (*field)[ixp - 1]) * flx;
        flux += ((*field)[ixp] - (*field)[ixp + 1]) * flx;
        flux += ((*field)[ixp] - (*field)[ixp - nx]) * fly;
        flux += ((*field)[ixp] - (*field)[ixp + nx]) * fly;
        flux += ((*field)[ixp] - (*field)[ixp - inxy]) * flz;
        flux += ((*field)[ixp] - (*field)[ixp + inxy]) * flz;
        delta = (*grid)[ixp] - flux * dv;
        if (fabs(delta) > deltamax)
          deltamax = fabs(delta);
        (*field)[ixp] += delta * 0.5;
      }
    }
  }
  return deltamax;
}


x86 assembly (clang, -O0):
<fd_sweep>:
push   %rbp
mov    %rsp,%rbp
movsd  %xmm2,-0x68(%rbp)
movaps %xmm0,%xmm2
movsd  -0x68(%rbp),%xmm0
cvtsd2ss %xmm2,%xmm2
cvtsd2ss %xmm1,%xmm1
cvtsd2ss %xmm0,%xmm0
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %edx,-0x14(%rbp)
mov    %ecx,-0x18(%rbp)
mov    %r8d,-0x1c(%rbp)
movss  %xmm2,-0x20(%rbp)
movss  %xmm1,-0x24(%rbp)
movss  %xmm0,-0x28(%rbp)
movss  -0x24(%rbp),%xmm0
mulss  -0x28(%rbp),%xmm0
divss  -0x20(%rbp),%xmm0
movss  %xmm0,-0x44(%rbp)
movss  -0x20(%rbp),%xmm0
mulss  -0x28(%rbp),%xmm0
divss  -0x24(%rbp),%xmm0
movss  %xmm0,-0x48(%rbp)
movss  -0x20(%rbp),%xmm0
mulss  -0x24(%rbp),%xmm0
divss  -0x28(%rbp),%xmm0
movss  %xmm0,-0x4c(%rbp)
movss  -0x20(%rbp),%xmm1
movss  -0x24(%rbp),%xmm2
movss  -0x24(%rbp),%xmm0
mulss  -0x28(%rbp),%xmm0
mulss  %xmm2,%xmm1
addss  %xmm0,%xmm1
movss  -0x20(%rbp),%xmm0
movss  -0x28(%rbp),%xmm2
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
cvtss2sd %xmm0,%xmm1
movsd  0xe5c(%rip),%xmm0
mulsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,-0x50(%rbp)
movss  -0x50(%rbp),%xmm0
cvtss2sd %xmm0,%xmm1
movsd  0xe36(%rip),%xmm0
divsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,-0x50(%rbp)
movss  -0x50(%rbp),%xmm0
cvtss2sd %xmm0,%xmm0
movsd  0xe10(%rip),%xmm1
divsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,-0x50(%rbp)
mov    -0x14(%rbp),%eax
imul   -0x18(%rbp),%eax
mov    %eax,-0x60(%rbp)
xorps  %xmm0,%xmm0
movss  %xmm0,-0x3c(%rbp)
movl   $0x1,-0x34(%rbp)
mov    -0x34(%rbp),%eax
mov    -0x1c(%rbp),%ecx
sub    $0x1,%ecx
cmp    %ecx,%eax
jge    148a <fd_sweep+0x38a>
mov    -0x34(%rbp),%eax
imul   -0x60(%rbp),%eax
mov    %eax,-0x54(%rbp)
movl   $0x1,-0x30(%rbp)
mov    -0x30(%rbp),%eax
mov    -0x18(%rbp),%ecx
sub    $0x1,%ecx
cmp    %ecx,%eax
jge    1477 <fd_sweep+0x377>
mov    -0x30(%rbp),%eax
imul   -0x14(%rbp),%eax
mov    %eax,-0x58(%rbp)
movl   $0x1,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
mov    -0x14(%rbp),%ecx
sub    $0x1,%ecx
cmp    %ecx,%eax
jge    1464 <fd_sweep+0x364>
mov    -0x2c(%rbp),%eax
add    -0x58(%rbp),%eax
add    -0x54(%rbp),%eax
mov    %eax,-0x5c(%rbp)
mov    -0x10(%rbp),%rax
movslq -0x5c(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x5c(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
subss  (%rax,%rcx,4),%xmm0
mulss  -0x44(%rbp),%xmm0
movss  %xmm0,-0x40(%rbp)
mov    -0x10(%rbp),%rax
movslq -0x5c(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x5c(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
subss  (%rax,%rcx,4),%xmm0
movss  -0x44(%rbp),%xmm2
movss  -0x40(%rbp),%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,-0x40(%rbp)
mov    -0x10(%rbp),%rax
movslq -0x5c(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x5c(%rbp),%ecx
sub    -0x14(%rbp),%ecx
movslq %ecx,%rcx
subss  (%rax,%rcx,4),%xmm0
movss  -0x48(%rbp),%xmm2
movss  -0x40(%rbp),%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,-0x40(%rbp)
mov    -0x10(%rbp),%rax
movslq -0x5c(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x5c(%rbp),%ecx
add    -0x14(%rbp),%ecx
movslq %ecx,%rcx
subss  (%rax,%rcx,4),%xmm0
movss  -0x48(%rbp),%xmm2
movss  -0x40(%rbp),%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,-0x40(%rbp)
mov    -0x10(%rbp),%rax
movslq -0x5c(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x5c(%rbp),%ecx
sub    -0x60(%rbp),%ecx
movslq %ecx,%rcx
subss  (%rax,%rcx,4),%xmm0
movss  -0x4c(%rbp),%xmm2
movss  -0x40(%rbp),%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,-0x40(%rbp)
mov    -0x10(%rbp),%rax
movslq -0x5c(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x5c(%rbp),%ecx
add    -0x60(%rbp),%ecx
movslq %ecx,%rcx
subss  (%rax,%rcx,4),%xmm0
movss  -0x4c(%rbp),%xmm2
movss  -0x40(%rbp),%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,-0x40(%rbp)
mov    -0x8(%rbp),%rax
movslq -0x5c(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm1
movss  -0x40(%rbp),%xmm0
movss  -0x50(%rbp),%xmm2
movd   %xmm0,%eax
xor    $0x80000000,%eax
movd   %eax,%xmm0
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,-0x38(%rbp)
movss  -0x38(%rbp),%xmm0
cvtss2sd %xmm0,%xmm0
movaps 0xc31(%rip),%xmm1
pand   %xmm1,%xmm0
movss  -0x3c(%rbp),%xmm1
cvtss2sd %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jbe    1423 <fd_sweep+0x323>
movss  -0x38(%rbp),%xmm0
cvtss2sd %xmm0,%xmm0
movaps 0xc0a(%rip),%xmm1
pand   %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,-0x3c(%rbp)
movss  -0x38(%rbp),%xmm0
cvtss2sd %xmm0,%xmm0
mov    -0x10(%rbp),%rax
movslq -0x5c(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm1
cvtss2sd %xmm1,%xmm1
movsd  0xbd3(%rip),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,(%rax,%rcx,4)
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    125a <fd_sweep+0x15a>
jmp    1469 <fd_sweep+0x369>
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
jmp    1238 <fd_sweep+0x138>
jmp    147c <fd_sweep+0x37c>
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    1216 <fd_sweep+0x116>
movss  -0x3c(%rbp),%xmm0
pop    %rbp
ret

x86 assembly (clang, -O2):
<fd_sweep>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
movaps %xmm0,%xmm3
mov    %rcx,-0x40(%rsp)
mov    %rdx,-0x38(%rsp)
xorps  %xmm0,%xmm0
mov    %r8d,-0x48(%rsp)
cmp    $0x3,%r8d
jl     138f <fd_sweep+0x28f>
xorps  %xmm0,%xmm0
cvtsd2ss %xmm3,%xmm0
cvtsd2ss %xmm1,%xmm4
cvtsd2ss %xmm2,%xmm5
movaps %xmm4,%xmm3
mulss  %xmm5,%xmm3
movaps %xmm3,%xmm1
divss  %xmm0,%xmm1
movaps %xmm0,%xmm6
mulss  %xmm5,%xmm6
mulss  %xmm4,%xmm0
movaps %xmm0,%xmm2
addss  %xmm3,%xmm0
addss  %xmm6,%xmm0
movaps %xmm6,%xmm3
divss  %xmm4,%xmm3
addss  %xmm0,%xmm0
movss  0xe91(%rip),%xmm4
divss  %xmm0,%xmm4
xorps  %xmm0,%xmm0
cvtss2sd %xmm4,%xmm0
divsd  0xe86(%rip),%xmm0
divss  %xmm5,%xmm2
cvtsd2ss %xmm0,%xmm10
mov    -0x40(%rsp),%rcx
mov    %ecx,%r11d
mov    -0x38(%rsp),%rbp
imul   %ebp,%r11d
addl   $0xffffffff,-0x48(%rsp)
lea    -0x1(%rcx),%edx
lea    -0x1(%rbp),%eax
lea    0x1(%rcx),%r9d
imul   %ebp,%r9d
mov    %r11d,%ebx
mov    %rbx,-0x30(%rsp)
mov    %ebp,%r10d
add    $0x1,%r11d
lea    0x2(%rcx),%r15d
imul   %ebp,%r15d
add    $0x1,%r15d
lea    0x1(%rbp),%r12d
lea    (%rcx,%rcx,1),%r13d
add    $0x1,%r13d
imul   %ebp,%r13d
add    $0x1,%r13d
add    $0xffffffffffffffff,%rax
xorps  %xmm0,%xmm0
mov    $0x1,%ebp
movaps 0xe2e(%rip),%xmm8
movsd  0xe15(%rip),%xmm9
jmp    123e <fd_sweep+0x13e>
nopl   (%rax)
mov    -0x44(%rsp),%ebp
add    $0x1,%ebp
mov    -0x10(%rsp),%r9
mov    -0x30(%rsp),%rcx
add    %rcx,%r9
mov    -0x8(%rsp),%r11
add    %rcx,%r11
mov    -0x18(%rsp),%r15
add    %rcx,%r15
mov    -0x20(%rsp),%r12
add    %rcx,%r12
mov    -0x28(%rsp),%r13
add    %rcx,%r13
cmp    -0x48(%rsp),%ebp
je     138f <fd_sweep+0x28f>
cmpl   $0x3,-0x40(%rsp)
mov    %r11,-0x8(%rsp)
mov    %r9,-0x10(%rsp)
mov    %r15,-0x18(%rsp)
mov    %r12,-0x20(%rsp)
mov    %r13,-0x28(%rsp)
mov    %ebp,-0x44(%rsp)
jl     1200 <fd_sweep+0x100>
mov    $0x1,%r14d
jmp    128c <fd_sweep+0x18c>
nopw   0x0(%rax,%rax,1)
add    $0x1,%r14d
add    %r10,%r9
add    %r10,%r11
add    %r10,%r15
add    %r10,%r12
add    %r10,%r13
cmp    %edx,%r14d
je     1200 <fd_sweep+0x100>
cmpl   $0x3,-0x38(%rsp)
jl     1270 <fd_sweep+0x170>
xor    %ecx,%ecx
jmp    12c6 <fd_sweep+0x1c6>
nopw   0x0(%rax,%rax,1)
xorps  %xmm4,%xmm4
cvtss2sd %xmm7,%xmm4
mulsd  %xmm9,%xmm5
addsd  %xmm4,%xmm5
xorps  %xmm4,%xmm4
cvtsd2ss %xmm5,%xmm4
movss  %xmm4,(%rsi,%r8,4)
add    $0x1,%rcx
cmp    %rcx,%rax
je     1270 <fd_sweep+0x170>
lea    (%r9,%rcx,1),%rbp
lea    0x1(%rbp),%ebx
movslq %ebx,%r8
movss  (%rsi,%r8,4),%xmm7
movslq %ebp,%rbx
movaps %xmm7,%xmm5
subss  (%rsi,%rbx,4),%xmm5
mulss  %xmm1,%xmm5
add    $0x2,%ebx
movslq %ebx,%rbx
movaps %xmm7,%xmm6
subss  (%rsi,%rbx,4),%xmm6
mulss  %xmm1,%xmm6
addss  %xmm5,%xmm6
lea    (%r11,%rcx,1),%ebx
movslq %ebx,%rbx
movaps %xmm7,%xmm5
subss  (%rsi,%rbx,4),%xmm5
mulss  %xmm3,%xmm5
addss  %xmm6,%xmm5
lea    (%r15,%rcx,1),%ebx
movslq %ebx,%rbx
movaps %xmm7,%xmm6
subss  (%rsi,%rbx,4),%xmm6
mulss  %xmm3,%xmm6
addss  %xmm5,%xmm6
lea    (%r12,%rcx,1),%ebx
movslq %ebx,%rbx
movaps %xmm7,%xmm5
subss  (%rsi,%rbx,4),%xmm5
mulss  %xmm2,%xmm5
addss  %xmm6,%xmm5
lea    0x0(%r13,%rcx,1),%ebx
movslq %ebx,%rbx
movaps %xmm7,%xmm6
subss  (%rsi,%rbx,4),%xmm6
mulss  %xmm2,%xmm6
addss  %xmm5,%xmm6
mulss  %xmm10,%xmm6
movss  (%rdi,%r8,4),%xmm5
subss  %xmm6,%xmm5
cvtss2sd %xmm5,%xmm5
movaps %xmm5,%xmm6
andps  %xmm8,%xmm6
xorps  %xmm4,%xmm4
cvtss2sd %xmm0,%xmm4
ucomisd %xmm4,%xmm6
jbe    12a0 <fd_sweep+0x1a0>
xorps  %xmm0,%xmm0
cvtsd2ss %xmm6,%xmm0
jmp    12a0 <fd_sweep+0x1a0>
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

LLVM IR (-O2):
define float @fd_sweep([0 x float]* nocapture %0, [0 x float]* nocapture %1, i32 %2, i32 %3, i32 %4, double %5, double %6, double %7) {
%9 = fptrunc double %5 to float
%10 = fptrunc double %6 to float
%11 = fptrunc double %7 to float
%12 = fmul float %10, %11
%13 = fdiv float %12, %9
%14 = fmul float %9, %11
%15 = fdiv float %14, %10
%16 = fmul float %9, %10
%17 = fdiv float %16, %11
%18 = tail call float @llvm.fmuladd.f32(float %9, float %10, float %12)
%19 = tail call float @llvm.fmuladd.f32(float %9, float %11, float %18)
%20 = fmul float %19, 2.000000e+00
%21 = fdiv float 1.000000e+00, %20
%22 = fpext float %21 to double
%23 = fdiv double %22, 0x402921FB54442D18
%24 = fptrunc double %23 to float
%25 = mul nsw i32 %3, %2
%26 = add i32 %4, -1
%27 = icmp sgt i32 %4, 2
br i1 %27, label %28, label %111
28: ; preds = %8
%29 = add i32 %3, -1
%30 = icmp sgt i32 %3, 2
%31 = add i32 %2, -1
%32 = icmp sgt i32 %2, 2
%33 = zext i32 %31 to i64
br label %34
34: ; preds = %28, %107
%35 = phi i32 [ 1, %28 ], [ %109, %107 ]
%36 = phi float [ 0.000000e+00, %28 ], [ %108, %107 ]
%37 = mul nsw i32 %35, %25
br i1 %30, label %38, label %107
38: ; preds = %34, %103
%39 = phi i32 [ %105, %103 ], [ 1, %34 ]
%40 = phi float [ %104, %103 ], [ %36, %34 ]
br i1 %32, label %41, label %103
41: ; preds = %38
%42 = mul nsw i32 %39, %2
%43 = add i32 %42, %37
br label %44
44: ; preds = %41, %44
%45 = phi i64 [ 1, %41 ], [ %101, %44 ]
%46 = phi float [ %40, %41 ], [ %97, %44 ]
%47 = trunc i64 %45 to i32
%48 = add i32 %43, %47
%49 = sext i32 %48 to i64
%50 = getelementptr inbounds [0 x float], [0 x float]* %1, i64 0, i64 %49
%51 = load float, float* %50
%52 = add nsw i32 %48, -1
%53 = sext i32 %52 to i64
%54 = getelementptr inbounds [0 x float], [0 x float]* %1, i64 0, i64 %53
%55 = load float, float* %54
%56 = fsub float %51, %55
%57 = fmul float %13, %56
%58 = add nsw i32 %48, 1
%59 = sext i32 %58 to i64
%60 = getelementptr inbounds [0 x float], [0 x float]* %1, i64 0, i64 %59
%61 = load float, float* %60
%62 = fsub float %51, %61
%63 = tail call float @llvm.fmuladd.f32(float %62, float %13, float %57)
%64 = sub nsw i32 %48, %2
%65 = sext i32 %64 to i64
%66 = getelementptr inbounds [0 x float], [0 x float]* %1, i64 0, i64 %65
%67 = load float, float* %66
%68 = fsub float %51, %67
%69 = tail call float @llvm.fmuladd.f32(float %68, float %15, float %63)
%70 = add nsw i32 %48, %2
%71 = sext i32 %70 to i64
%72 = getelementptr inbounds [0 x float], [0 x float]* %1, i64 0, i64 %71
%73 = load float, float* %72
%74 = fsub float %51, %73
%75 = tail call float @llvm.fmuladd.f32(float %74, float %15, float %69)
%76 = sub nsw i32 %48, %25
%77 = sext i32 %76 to i64
%78 = getelementptr inbounds [0 x float], [0 x float]* %1, i64 0, i64 %77
%79 = load float, float* %78
%80 = fsub float %51, %79
%81 = tail call float @llvm.fmuladd.f32(float %80, float %17, float %75)
%82 = add nsw i32 %48, %25
%83 = sext i32 %82 to i64
%84 = getelementptr inbounds [0 x float], [0 x float]* %1, i64 0, i64 %83
%85 = load float, float* %84
%86 = fsub float %51, %85
%87 = tail call float @llvm.fmuladd.f32(float %86, float %17, float %81)
%88 = getelementptr inbounds [0 x float], [0 x float]* %0, i64 0, i64 %49
%89 = load float, float* %88
%90 = fneg float %87
%91 = tail call float @llvm.fmuladd.f32(float %90, float %24, float %89)
%92 = fpext float %91 to double
%93 = tail call double @llvm.fabs.f64(double %92)
%94 = fpext float %46 to double
%95 = fcmp ogt double %93, %94
%96 = fptrunc double %93 to float
%97 = select i1 %95, float %96, float %46
%98 = fpext float %51 to double
%99 = tail call double @llvm.fmuladd.f64(double %92, double 5.000000e-01, double %98)
%100 = fptrunc double %99 to float
store float %100, float* %50
%101 = add nuw nsw i64 %45, 1
%102 = icmp eq i64 %101, %33
br i1 %102, label %103, label %44
103: ; preds = %44, %38
%104 = phi float [ %40, %38 ], [ %97, %44 ]
%105 = add nuw nsw i32 %39, 1
%106 = icmp eq i32 %105, %29
br i1 %106, label %107, label %38
107: ; preds = %103, %34
%108 = phi float [ %36, %34 ], [ %104, %103 ]
%109 = add nuw nsw i32 %35, 1
%110 = icmp eq i32 %109, %26
br i1 %110, label %111, label %34
111: ; preds = %107, %8
%112 = phi float [ 0.000000e+00, %8 ], [ %108, %107 ]
ret float %112
}