void Average_block(float *ima, int x, int y, int z, int neighborhoodsize,
                   float *average, double weight, int sx, int sy, int sz) {
  int x_pos, y_pos, z_pos;
  int is_outside;
  int a, b, c, ns, sxy, index;
  extern int rician;
  ns = 2 * neighborhoodsize + 1;
  sxy = sx * sy;
  for (c = 0; c < ns; c++) {
    for (b = 0; b < ns; b++) {
      for (a = 0; a < ns; a++) {
        is_outside = 0;
        x_pos = x + a - neighborhoodsize;
        y_pos = y + b - neighborhoodsize;
        z_pos = z + c - neighborhoodsize;
        if ((z_pos < 0) || (z_pos > sz - 1))
          is_outside = 1;
        if ((y_pos < 0) || (y_pos > sy - 1))
          is_outside = 1;
        if ((x_pos < 0) || (x_pos > sx - 1))
          is_outside = 1;
        index = a + ns * b + ns * ns * c;
        if (rician) {
          if (is_outside)
            average[index] += ima[z * (sxy) + (y * sx) + x] *
                              ima[z * (sxy) + (y * sx) + x] * (float)weight;
          else
            average[index] += ima[z_pos * (sxy) + (y_pos * sx) + x_pos] *
                              ima[z_pos * (sxy) + (y_pos * sx) + x_pos] *
                              (float)weight;
        } else {
          if (is_outside)
            average[index] += ima[z * (sxy) + (y * sx) + x] * (float)weight;
          else
            average[index] +=
                ima[z_pos * (sxy) + (y_pos * sx) + x_pos] * (float)weight;
        }
      }
    }
  }
}


x86 assembly (clang, -O0):
<Average_block>:
push   %rbp
mov    %rsp,%rbp
mov    0x20(%rbp),%eax
mov    0x18(%rbp),%eax
mov    0x10(%rbp),%eax
mov    %rdi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
mov    %edx,-0x10(%rbp)
mov    %ecx,-0x14(%rbp)
mov    %r8d,-0x18(%rbp)
mov    %r9,-0x20(%rbp)
movsd  %xmm0,-0x28(%rbp)
mov    -0x18(%rbp),%eax
shl    $0x1,%eax
add    $0x1,%eax
mov    %eax,-0x48(%rbp)
mov    0x10(%rbp),%eax
imul   0x18(%rbp),%eax
mov    %eax,-0x4c(%rbp)
movl   $0x0,-0x44(%rbp)
mov    -0x44(%rbp),%eax
cmp    -0x48(%rbp),%eax
jge    13da <Average_block+0x2da>
movl   $0x0,-0x40(%rbp)
mov    -0x40(%rbp),%eax
cmp    -0x48(%rbp),%eax
jge    13c7 <Average_block+0x2c7>
movl   $0x0,-0x3c(%rbp)
mov    -0x3c(%rbp),%eax
cmp    -0x48(%rbp),%eax
jge    13b4 <Average_block+0x2b4>
movl   $0x0,-0x38(%rbp)
mov    -0xc(%rbp),%eax
add    -0x3c(%rbp),%eax
sub    -0x18(%rbp),%eax
mov    %eax,-0x2c(%rbp)
mov    -0x10(%rbp),%eax
add    -0x40(%rbp),%eax
sub    -0x18(%rbp),%eax
mov    %eax,-0x30(%rbp)
mov    -0x14(%rbp),%eax
add    -0x44(%rbp),%eax
sub    -0x18(%rbp),%eax
mov    %eax,-0x34(%rbp)
cmpl   $0x0,-0x34(%rbp)
jl     11bc <Average_block+0xbc>
mov    -0x34(%rbp),%eax
mov    0x20(%rbp),%ecx
sub    $0x1,%ecx
cmp    %ecx,%eax
jle    11c3 <Average_block+0xc3>
movl   $0x1,-0x38(%rbp)
cmpl   $0x0,-0x30(%rbp)
jl     11de <Average_block+0xde>
mov    -0x30(%rbp),%eax
mov    0x18(%rbp),%ecx
sub    $0x1,%ecx
cmp    %ecx,%eax
jle    11e5 <Average_block+0xe5>
movl   $0x1,-0x38(%rbp)
cmpl   $0x0,-0x2c(%rbp)
jl     1200 <Average_block+0x100>
mov    -0x2c(%rbp),%eax
mov    0x10(%rbp),%ecx
sub    $0x1,%ecx
cmp    %ecx,%eax
jle    1207 <Average_block+0x107>
movl   $0x1,-0x38(%rbp)
mov    -0x3c(%rbp),%eax
mov    -0x48(%rbp),%ecx
imul   -0x40(%rbp),%ecx
add    %ecx,%eax
mov    -0x48(%rbp),%ecx
imul   -0x48(%rbp),%ecx
imul   -0x44(%rbp),%ecx
add    %ecx,%eax
mov    %eax,-0x50(%rbp)
mov    0x2db6(%rip),%rax
cmpl   $0x0,(%rax)
je     1309 <Average_block+0x209>
cmpl   $0x0,-0x38(%rbp)
je     12a3 <Average_block+0x1a3>
mov    -0x8(%rbp),%rax
mov    -0x14(%rbp),%ecx
imul   -0x4c(%rbp),%ecx
mov    -0x10(%rbp),%edx
imul   0x10(%rbp),%edx
add    %edx,%ecx
add    -0xc(%rbp),%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x14(%rbp),%ecx
imul   -0x4c(%rbp),%ecx
mov    -0x10(%rbp),%edx
imul   0x10(%rbp),%edx
add    %edx,%ecx
add    -0xc(%rbp),%ecx
movslq %ecx,%rcx
mulss  (%rax,%rcx,4),%xmm0
movsd  -0x28(%rbp),%xmm1
cvtsd2ss %xmm1,%xmm2
mov    -0x20(%rbp),%rax
movslq -0x50(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,(%rax,%rcx,4)
jmp    1304 <Average_block+0x204>
mov    -0x8(%rbp),%rax
mov    -0x34(%rbp),%ecx
imul   -0x4c(%rbp),%ecx
mov    -0x30(%rbp),%edx
imul   0x10(%rbp),%edx
add    %edx,%ecx
add    -0x2c(%rbp),%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x34(%rbp),%ecx
imul   -0x4c(%rbp),%ecx
mov    -0x30(%rbp),%edx
imul   0x10(%rbp),%edx
add    %edx,%ecx
add    -0x2c(%rbp),%ecx
movslq %ecx,%rcx
mulss  (%rax,%rcx,4),%xmm0
movsd  -0x28(%rbp),%xmm1
cvtsd2ss %xmm1,%xmm2
mov    -0x20(%rbp),%rax
movslq -0x50(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,(%rax,%rcx,4)
jmp    13a1 <Average_block+0x2a1>
cmpl   $0x0,-0x38(%rbp)
je     135a <Average_block+0x25a>
mov    -0x8(%rbp),%rax
mov    -0x14(%rbp),%ecx
imul   -0x4c(%rbp),%ecx
mov    -0x10(%rbp),%edx
imul   0x10(%rbp),%edx
add    %edx,%ecx
add    -0xc(%rbp),%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm0
movsd  -0x28(%rbp),%xmm1
cvtsd2ss %xmm1,%xmm2
mov    -0x20(%rbp),%rax
movslq -0x50(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,(%rax,%rcx,4)
jmp    139c <Average_block+0x29c>
mov    -0x8(%rbp),%rax
mov    -0x34(%rbp),%ecx
imul   -0x4c(%rbp),%ecx
mov    -0x30(%rbp),%edx
imul   0x10(%rbp),%edx
add    %edx,%ecx
add    -0x2c(%rbp),%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm0
movsd  -0x28(%rbp),%xmm1
cvtsd2ss %xmm1,%xmm2
mov    -0x20(%rbp),%rax
movslq -0x50(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,(%rax,%rcx,4)
jmp    13a1 <Average_block+0x2a1>
jmp    13a6 <Average_block+0x2a6>
mov    -0x3c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x3c(%rbp)
jmp    116a <Average_block+0x6a>
jmp    13b9 <Average_block+0x2b9>
mov    -0x40(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x40(%rbp)
jmp    1157 <Average_block+0x57>
jmp    13cc <Average_block+0x2cc>
mov    -0x44(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x44(%rbp)
jmp    1144 <Average_block+0x44>
pop    %rbp
ret

x86 assembly (clang, -O2):
<Average_block>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
lea    (%r8,%r8,1),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rsp)
test   %eax,%eax
jle    12a9 <Average_block+0x1a9>
mov    %edx,%r11d
mov    0x40(%rsp),%r10d
mov    0x38(%rsp),%ebx
mov    %r10d,%eax
imul   %ebx,%eax
sub    %r8d,%edx
mov    %edx,-0x10(%rsp)
cvtsd2ss %xmm0,%xmm0
mov    %eax,-0x24(%rsp)
imul   %ecx,%eax
sub    %r8d,%ecx
mov    %ebx,%edx
imul   %r11d,%edx
add    %esi,%edx
add    %eax,%edx
movslq %edx,%rbp
mov    -0x2c(%rsp),%eax
mov    %eax,%edx
mov    %rdx,-0x8(%rsp)
imul   %eax,%eax
mov    %eax,-0x28(%rsp)
mov    %ecx,-0x20(%rsp)
imul   %ecx,%r10d
add    %r11d,%r10d
sub    %r8d,%r10d
imul   %ebx,%r10d
sub    %r8d,%esi
mov    %esi,-0xc(%rsp)
xor    %eax,%eax
xor    %edx,%edx
jmp    11b3 <Average_block+0xb3>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
mov    -0x1c(%rsp),%edx
add    $0x1,%edx
mov    -0x18(%rsp),%eax
add    -0x28(%rsp),%eax
mov    -0x14(%rsp),%r10d
add    -0x24(%rsp),%r10d
cmp    -0x2c(%rsp),%edx
je     12a9 <Average_block+0x1a9>
mov    -0x20(%rsp),%ecx
mov    %edx,-0x1c(%rsp)
add    %edx,%ecx
sets   %dl
cmp    0x48(%rsp),%ecx
setge  %r14b
or     %dl,%r14b
mov    0x2e0e(%rip),%rcx
mov    (%rcx),%ecx
mov    %r10d,-0x14(%rsp)
mov    %r10d,%r13d
mov    %eax,-0x18(%rsp)
mov    %eax,%edx
xor    %r10d,%r10d
jmp    1202 <Average_block+0x102>
nopw   0x0(%rax,%rax,1)
add    $0x1,%r10d
mov    -0x2c(%rsp),%eax
add    %eax,%edx
add    %ebx,%r13d
cmp    %eax,%r10d
je     1190 <Average_block+0x90>
mov    -0x10(%rsp),%esi
add    %r10d,%esi
sets   %al
cmp    0x40(%rsp),%esi
setge  %r8b
or     %al,%r8b
mov    -0x8(%rsp),%r12
mov    -0xc(%rsp),%eax
mov    %eax,%r15d
mov    %edx,%esi
jmp    125f <Average_block+0x15f>
nopw   0x0(%rax,%rax,1)
lea    (%r15,%r13,1),%eax
cltq
movss  (%rdi,%rax,4),%xmm1
mulss  %xmm1,%xmm1
mulss  %xmm0,%xmm1
movslq %esi,%rsi
addss  (%r9,%rsi,4),%xmm1
movss  %xmm1,(%r9,%rsi,4)
add    $0x1,%esi
add    $0x1,%r15d
add    $0xffffffffffffffff,%r12
je     11f0 <Average_block+0xf0>
test   %r15d,%r15d
sets   %al
cmp    %ebx,%r15d
setge  %r11b
or     %al,%r11b
or     %r8b,%r11b
or     %r14b,%r11b
test   %ecx,%ecx
je     1290 <Average_block+0x190>
test   %r11b,%r11b
je     1230 <Average_block+0x130>
movss  (%rdi,%rbp,4),%xmm1
jmp    123b <Average_block+0x13b>
cs nopw 0x0(%rax,%rax,1)
nop
test   %r11b,%r11b
je     129c <Average_block+0x19c>
movss  (%rdi,%rbp,4),%xmm1
jmp    123f <Average_block+0x13f>
lea    (%r15,%r13,1),%eax
cltq
movss  (%rdi,%rax,4),%xmm1
jmp    123f <Average_block+0x13f>
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

LLVM IR (-O2):
@rician = external global i32
define void @Average_block(float* nocapture %0, i32 %1, i32 %2, i32 %3, i32 %4, float* nocapture %5, double %6, i32 %7, i32 %8, i32 %9) {
%11 = shl nsw i32 %4, 1
%12 = or i32 %11, 1
%13 = mul nsw i32 %8, %7
%14 = icmp sgt i32 %12, 0
br i1 %14, label %15, label %89
15: ; preds = %10
%16 = sub i32 %1, %4
%17 = sub i32 %2, %4
%18 = sub i32 %3, %4
%19 = fptrunc double %6 to float
%20 = mul nsw i32 %13, %3
%21 = mul nsw i32 %7, %2
%22 = add i32 %21, %1
%23 = add i32 %22, %20
%24 = sext i32 %23 to i64
%25 = getelementptr inbounds float, float* %0, i64 %24
%26 = zext i32 %12 to i64
br label %27
27: ; preds = %86, %15
%28 = phi i32 [ 0, %15 ], [ %87, %86 ]
%29 = add i32 %18, %28
%30 = icmp slt i32 %29, 0
%31 = icmp sge i32 %29, %9
%32 = or i1 %30, %31
%33 = mul i32 %28, %12
%34 = load i32, i32* @rician
%35 = icmp eq i32 %34, 0
%36 = mul nsw i32 %29, %13
br label %37
37: ; preds = %83, %27
%38 = phi i32 [ 0, %27 ], [ %84, %83 ]
%39 = add i32 %17, %38
%40 = icmp slt i32 %39, 0
%41 = icmp sge i32 %39, %8
%42 = or i1 %40, %41
%43 = add i32 %38, %33
%44 = mul i32 %43, %12
%45 = mul nsw i32 %39, %7
%46 = add nsw i32 %45, %36
br label %47
47: ; preds = %37, %75
%48 = phi i64 [ 0, %37 ], [ %81, %75 ]
%49 = trunc i64 %48 to i32
%50 = add i32 %16, %49
%51 = icmp slt i32 %50, 0
%52 = icmp sge i32 %50, %7
%53 = or i1 %51, %52
%54 = select i1 %53, i1 true, i1 %42
%55 = select i1 %54, i1 true, i1 %32
%56 = add i32 %44, %49
br i1 %35, label %67, label %57
57: ; preds = %47
br i1 %55, label %58, label %61
58: ; preds = %57
%59 = load float, float* %25
%60 = fmul float %59, %59
br label %75
61: ; preds = %57
%62 = add nsw i32 %46, %50
%63 = sext i32 %62 to i64
%64 = getelementptr inbounds float, float* %0, i64 %63
%65 = load float, float* %64
%66 = fmul float %65, %65
br label %75
67: ; preds = %47
br i1 %55, label %68, label %70
68: ; preds = %67
%69 = load float, float* %25
br label %75
70: ; preds = %67
%71 = add nsw i32 %46, %50
%72 = sext i32 %71 to i64
%73 = getelementptr inbounds float, float* %0, i64 %72
%74 = load float, float* %73
br label %75
75: ; preds = %61, %58, %70, %68
%76 = phi float [ %66, %61 ], [ %60, %58 ], [ %74, %70 ], [ %69, %68 ]
%77 = sext i32 %56 to i64
%78 = getelementptr inbounds float, float* %5, i64 %77
%79 = load float, float* %78
%80 = tail call float @llvm.fmuladd.f32(float %76, float %19, float %79)
store float %80, float* %78
%81 = add nuw nsw i64 %48, 1
%82 = icmp eq i64 %81, %26
br i1 %82, label %83, label %47
83: ; preds = %75
%84 = add nuw nsw i32 %38, 1
%85 = icmp eq i32 %84, %12
br i1 %85, label %86, label %37
86: ; preds = %83
%87 = add nuw nsw i32 %28, 1
%88 = icmp eq i32 %87, %12
br i1 %88, label %89, label %27
89: ; preds = %86, %10
ret void
}