void init4d(int dim1, int dim2, int dim3, int dim4, double *a) {
  double shift1, shift2, shift3;
  if (dim1 < 100 && dim2 < 100 && dim3 < 100 && dim4 < 100)
    shift1 = 0.001;
  else
    shift1 = 0.0001;
  shift2 = shift1 * shift1;
  shift3 = shift2 * shift1;
  for (int i = 0; i < dim1; i++)
    for (int j = 0; j < dim2; j++)
      for (int k = 0; k < dim3; k++)
        for (int l = 0; l < dim4; l++)
          a[i * dim2 * dim3 * dim4 + j * dim3 * dim4 + k * dim4 + l] =
              (double)(i + j * shift1 + k * shift2 + l * shift3);
  return;
}


x86 assembly (clang, -O0):
<init4d>:
push   %rbp
mov    %rsp,%rbp
mov    %edi,-0x4(%rbp)
mov    %esi,-0x8(%rbp)
mov    %edx,-0xc(%rbp)
mov    %ecx,-0x10(%rbp)
mov    %r8,-0x18(%rbp)
cmpl   $0x64,-0x4(%rbp)
jge    114e <init4d+0x4e>
cmpl   $0x64,-0x8(%rbp)
jge    114e <init4d+0x4e>
cmpl   $0x64,-0xc(%rbp)
jge    114e <init4d+0x4e>
cmpl   $0x64,-0x10(%rbp)
jge    114e <init4d+0x4e>
movsd  0xec4(%rip),%xmm0
movsd  %xmm0,-0x20(%rbp)
jmp    115b <init4d+0x5b>
movsd  0xeaa(%rip),%xmm0
movsd  %xmm0,-0x20(%rbp)
movsd  -0x20(%rbp),%xmm0
mulsd  -0x20(%rbp),%xmm0
movsd  %xmm0,-0x28(%rbp)
movsd  -0x28(%rbp),%xmm0
mulsd  -0x20(%rbp),%xmm0
movsd  %xmm0,-0x30(%rbp)
movl   $0x0,-0x34(%rbp)
mov    -0x34(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    127b <init4d+0x17b>
movl   $0x0,-0x38(%rbp)
mov    -0x38(%rbp),%eax
cmp    -0x8(%rbp),%eax
jge    1268 <init4d+0x168>
movl   $0x0,-0x3c(%rbp)
mov    -0x3c(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    1255 <init4d+0x155>
movl   $0x0,-0x40(%rbp)
mov    -0x40(%rbp),%eax
cmp    -0x10(%rbp),%eax
jge    1242 <init4d+0x142>
cvtsi2sdl -0x34(%rbp),%xmm1
cvtsi2sdl -0x38(%rbp),%xmm0
movsd  -0x20(%rbp),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
cvtsi2sdl -0x3c(%rbp),%xmm1
movsd  -0x28(%rbp),%xmm2
mulsd  %xmm2,%xmm1
addsd  %xmm0,%xmm1
cvtsi2sdl -0x40(%rbp),%xmm0
movsd  -0x30(%rbp),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x18(%rbp),%rax
mov    -0x34(%rbp),%ecx
imul   -0x8(%rbp),%ecx
imul   -0xc(%rbp),%ecx
imul   -0x10(%rbp),%ecx
mov    -0x38(%rbp),%edx
imul   -0xc(%rbp),%edx
imul   -0x10(%rbp),%edx
add    %edx,%ecx
mov    -0x3c(%rbp),%edx
imul   -0x10(%rbp),%edx
add    %edx,%ecx
add    -0x40(%rbp),%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x40(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x40(%rbp)
jmp    11b9 <init4d+0xb9>
jmp    1247 <init4d+0x147>
mov    -0x3c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x3c(%rbp)
jmp    11a6 <init4d+0xa6>
jmp    125a <init4d+0x15a>
mov    -0x38(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x38(%rbp)
jmp    1193 <init4d+0x93>
jmp    126d <init4d+0x16d>
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    1180 <init4d+0x80>
pop    %rbp
ret

x86 assembly (clang, -O2):
<init4d>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
mov    %edi,-0x4c(%rsp)
cmp    $0x64,%edi
setl   %al
mov    %esi,-0x58(%rsp)
cmp    $0x64,%esi
setl   %bl
and    %al,%bl
mov    %edx,-0x50(%rsp)
cmp    $0x64,%edx
setl   %al
and    %bl,%al
cmp    $0x64,%ecx
setl   %dl
test   %dl,%al
jne    1165 <init4d+0x65>
movsd  0xeca(%rip),%xmm0
movsd  %xmm0,-0x30(%rsp)
je     1175 <init4d+0x75>
movsd  0xec1(%rip),%xmm13
je     1180 <init4d+0x80>
movsd  0xec7(%rip),%xmm2
cmpl   $0x0,-0x4c(%rsp)
jg     1193 <init4d+0x93>
jmp    1452 <init4d+0x352>
movsd  0xe93(%rip),%xmm0
movsd  %xmm0,-0x30(%rsp)
jne    1146 <init4d+0x46>
movsd  0xe9a(%rip),%xmm13
jne    1151 <init4d+0x51>
movsd  0xea0(%rip),%xmm2
cmpl   $0x0,-0x4c(%rsp)
jle    1452 <init4d+0x352>
mov    -0x50(%rsp),%edx
mov    %edx,%r9d
mov    %ecx,%eax
mov    %eax,%r10d
and    $0xfffffffc,%r10d
lea    -0x4(%r10),%rsi
mov    %rsi,-0x8(%rsp)
mov    %rsi,%rdi
shr    $0x2,%rdi
add    $0x1,%rdi
movapd %xmm2,%xmm3
unpcklpd %xmm2,%xmm3
mov    %rdi,%rsi
mov    %rdi,-0x10(%rsp)
and    $0xfffffffffffffffe,%rdi
mov    %rdi,-0x18(%rsp)
lea    0x30(%r8),%rsi
mov    %rsi,-0x20(%rsp)
mov    %ecx,%esi
imul   %edx,%esi
mov    %esi,-0x34(%rsp)
imul   -0x58(%rsp),%esi
mov    %esi,-0x38(%rsp)
xorpd  %xmm0,%xmm0
movsd  %xmm0,-0x40(%rsp)
xor    %edx,%edx
mov    %rdx,-0x48(%rsp)
movdqa 0xe3b(%rip),%xmm8
movdqa 0xe43(%rip),%xmm6
movdqa 0xe4a(%rip),%xmm15
movdqa 0xe51(%rip),%xmm9
movdqa 0xe59(%rip),%xmm4
movl   $0x0,-0x54(%rsp)
jmp    1277 <init4d+0x177>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
movsd  -0x40(%rsp),%xmm0
addsd  0xde2(%rip),%xmm0
movsd  %xmm0,-0x40(%rsp)
mov    -0x54(%rsp),%esi
add    $0x1,%esi
mov    -0x48(%rsp),%rdx
add    -0x38(%rsp),%edx
mov    %rdx,-0x48(%rsp)
mov    %esi,-0x54(%rsp)
cmp    -0x4c(%rsp),%esi
je     1452 <init4d+0x352>
cmpl   $0x0,-0x58(%rsp)
jle    1240 <init4d+0x140>
mov    -0x54(%rsp),%edx
imul   -0x58(%rsp),%edx
mov    %rdx,-0x28(%rsp)
xor    %r13d,%r13d
xorpd  %xmm12,%xmm12
mov    -0x48(%rsp),%rdx
mov    %edx,%ebp
jmp    12b8 <init4d+0x1b8>
nopl   (%rax)
addsd  0xd87(%rip),%xmm12
add    $0x1,%r13d
add    -0x34(%rsp),%ebp
cmp    -0x58(%rsp),%r13d
je     1240 <init4d+0x140>
cmpl   $0x0,-0x50(%rsp)
jle    12a0 <init4d+0x1a0>
movapd %xmm12,%xmm14
mulsd  -0x30(%rsp),%xmm14
addsd  -0x40(%rsp),%xmm14
mov    -0x28(%rsp),%rdx
lea    (%rdx,%r13,1),%edi
imul   -0x50(%rsp),%edi
mov    %ebp,%r12d
xor    %esi,%esi
jmp    12fc <init4d+0x1fc>
nopw   0x0(%rax,%rax,1)
add    $0x1,%rsi
add    %ecx,%r12d
cmp    %r9,%rsi
je     12a0 <init4d+0x1a0>
test   %ecx,%ecx
jle    12f0 <init4d+0x1f0>
xorps  %xmm5,%xmm5
cvtsi2sd %esi,%xmm5
movslq %r12d,%r15
mulsd  %xmm13,%xmm5
addsd  %xmm14,%xmm5
cmp    $0x4,%ecx
jae    1320 <init4d+0x220>
xor    %edx,%edx
jmp    1415 <init4d+0x315>
movapd %xmm5,%xmm1
unpcklpd %xmm5,%xmm1
cmpq   $0x0,-0x8(%rsp)
je     1442 <init4d+0x342>
mov    -0x20(%rsp),%rdx
lea    (%rdx,%r15,8),%r14
mov    -0x18(%rsp),%r11
xor    %edx,%edx
movdqa %xmm8,%xmm11
nopl   0x0(%rax)
movdqa %xmm11,%xmm7
cvtdq2pd %xmm11,%xmm10
paddd  %xmm6,%xmm7
cvtdq2pd %xmm7,%xmm7
mulpd  %xmm3,%xmm10
addpd  %xmm1,%xmm10
mulpd  %xmm3,%xmm7
addpd  %xmm1,%xmm7
movupd %xmm10,-0x30(%r14,%rdx,8)
movupd %xmm7,-0x20(%r14,%rdx,8)
movdqa %xmm11,%xmm7
paddd  %xmm15,%xmm7
movdqa %xmm11,%xmm0
cvtdq2pd %xmm7,%xmm7
paddd  %xmm9,%xmm0
cvtdq2pd %xmm0,%xmm0
mulpd  %xmm3,%xmm7
addpd  %xmm1,%xmm7
mulpd  %xmm3,%xmm0
addpd  %xmm1,%xmm0
movupd %xmm7,-0x10(%r14,%rdx,8)
movupd %xmm0,(%r14,%rdx,8)
add    $0x8,%rdx
paddd  %xmm4,%xmm11
add    $0xfffffffffffffffe,%r11
jne    1350 <init4d+0x250>
testb  $0x1,-0x10(%rsp)
je     1409 <init4d+0x309>
lea    (%rdi,%rsi,1),%ebx
cvtdq2pd %xmm11,%xmm0
imul   %ecx,%ebx
paddd  %xmm6,%xmm11
cvtdq2pd %xmm11,%xmm7
movslq %ebx,%rbx
mulpd  %xmm3,%xmm0
addpd  %xmm1,%xmm0
mulpd  %xmm3,%xmm7
addpd  %xmm1,%xmm7
add    %rbx,%rdx
movupd %xmm0,(%r8,%rdx,8)
movupd %xmm7,0x10(%r8,%rdx,8)
mov    %r10,%rdx
cmp    %rax,%r10
je     12f0 <init4d+0x1f0>
lea    (%r8,%r15,8),%rbx
nopl   0x0(%rax)
xorps  %xmm0,%xmm0
cvtsi2sd %edx,%xmm0
mulsd  %xmm2,%xmm0
addsd  %xmm5,%xmm0
movsd  %xmm0,(%rbx,%rdx,8)
add    $0x1,%rdx
cmp    %rdx,%rax
jne    1420 <init4d+0x320>
jmp    12f0 <init4d+0x1f0>
xor    %edx,%edx
movdqa %xmm8,%xmm11
testb  $0x1,-0x10(%rsp)
jne    13d1 <init4d+0x2d1>
jmp    1409 <init4d+0x309>
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

LLVM IR (-O0):
define void @init4d(i32 %0, i32 %1, i32 %2, i32 %3, double* %4) {
%6 = alloca i32
%7 = alloca i32
%8 = alloca i32
%9 = alloca i32
%10 = alloca double*
%11 = alloca double
%12 = alloca double
%13 = alloca double
%14 = alloca i32
%15 = alloca i32
%16 = alloca i32
%17 = alloca i32
store i32 %0, i32* %6
store i32 %1, i32* %7
store i32 %2, i32* %8
store i32 %3, i32* %9
store double* %4, double** %10
%18 = load i32, i32* %6
%19 = icmp slt i32 %18, 100
br i1 %19, label %20, label %30
20: ; preds = %5
%21 = load i32, i32* %7
%22 = icmp slt i32 %21, 100
br i1 %22, label %23, label %30
23: ; preds = %20
%24 = load i32, i32* %8
%25 = icmp slt i32 %24, 100
br i1 %25, label %26, label %30
26: ; preds = %23
%27 = load i32, i32* %9
%28 = icmp slt i32 %27, 100
br i1 %28, label %29, label %30
29: ; preds = %26
store double 1.000000e-03, double* %11
br label %31
30: ; preds = %26, %23, %20, %5
store double 1.000000e-04, double* %11
br label %31
31: ; preds = %30, %29
%32 = load double, double* %11
%33 = load double, double* %11
%34 = fmul double %32, %33
store double %34, double* %12
%35 = load double, double* %12
%36 = load double, double* %11
%37 = fmul double %35, %36
store double %37, double* %13
store i32 0, i32* %14
br label %38
38: ; preds = %106, %31
%39 = load i32, i32* %14
%40 = load i32, i32* %6
%41 = icmp slt i32 %39, %40
br i1 %41, label %42, label %109
42: ; preds = %38
store i32 0, i32* %15
br label %43
43: ; preds = %102, %42
%44 = load i32, i32* %15
%45 = load i32, i32* %7
%46 = icmp slt i32 %44, %45
br i1 %46, label %47, label %105
47: ; preds = %43
store i32 0, i32* %16
br label %48
48: ; preds = %98, %47
%49 = load i32, i32* %16
%50 = load i32, i32* %8
%51 = icmp slt i32 %49, %50
br i1 %51, label %52, label %101
52: ; preds = %48
store i32 0, i32* %17
br label %53
53: ; preds = %94, %52
%54 = load i32, i32* %17
%55 = load i32, i32* %9
%56 = icmp slt i32 %54, %55
br i1 %56, label %57, label %97
57: ; preds = %53
%58 = load i32, i32* %14
%59 = sitofp i32 %58 to double
%60 = load i32, i32* %15
%61 = sitofp i32 %60 to double
%62 = load double, double* %11
%63 = call double @llvm.fmuladd.f64(double %61, double %62, double %59)
%64 = load i32, i32* %16
%65 = sitofp i32 %64 to double
%66 = load double, double* %12
%67 = call double @llvm.fmuladd.f64(double %65, double %66, double %63)
%68 = load i32, i32* %17
%69 = sitofp i32 %68 to double
%70 = load double, double* %13
%71 = call double @llvm.fmuladd.f64(double %69, double %70, double %67)
%72 = load double*, double** %10
%73 = load i32, i32* %14
%74 = load i32, i32* %7
%75 = mul nsw i32 %73, %74
%76 = load i32, i32* %8
%77 = mul nsw i32 %75, %76
%78 = load i32, i32* %9
%79 = mul nsw i32 %77, %78
%80 = load i32, i32* %15
%81 = load i32, i32* %8
%82 = mul nsw i32 %80, %81
%83 = load i32, i32* %9
%84 = mul nsw i32 %82, %83
%85 = add nsw i32 %79, %84
%86 = load i32, i32* %16
%87 = load i32, i32* %9
%88 = mul nsw i32 %86, %87
%89 = add nsw i32 %85, %88
%90 = load i32, i32* %17
%91 = add nsw i32 %89, %90
%92 = sext i32 %91 to i64
%93 = getelementptr inbounds double, double* %72, i64 %92
store double %71, double* %93
br label %94
94: ; preds = %57
%95 = load i32, i32* %17
%96 = add nsw i32 %95, 1
store i32 %96, i32* %17
br label %53
97: ; preds = %53
br label %98
98: ; preds = %97
%99 = load i32, i32* %16
%100 = add nsw i32 %99, 1
store i32 %100, i32* %16
br label %48
101: ; preds = %48
br label %102
102: ; preds = %101
%103 = load i32, i32* %15
%104 = add nsw i32 %103, 1
store i32 %104, i32* %15
br label %43
105: ; preds = %43
br label %106
106: ; preds = %105
%107 = load i32, i32* %14
%108 = add nsw i32 %107, 1
store i32 %108, i32* %14
br label %38
109: ; preds = %38
ret void
}

LLVM IR (-O2):
define void @init4d(i32 %0, i32 %1, i32 %2, i32 %3, double* nocapture writeonly %4) {
%6 = icmp slt i32 %0, 100
%7 = icmp slt i32 %1, 100
%8 = and i1 %6, %7
%9 = icmp slt i32 %2, 100
%10 = and i1 %8, %9
%11 = icmp slt i32 %3, 100
%12 = and i1 %10, %11
%13 = select i1 %12, double 1.000000e-03, double 1.000000e-04
%14 = select i1 %12, double 0x3EB0C6F7A0B5ED8D, double 1.000000e-08
%15 = select i1 %12, double 1.000000e-09, double 0x3D719799812DEA11
%16 = icmp sgt i32 %0, 0
br i1 %16, label %17, label %43
17: ; preds = %5
%18 = icmp sgt i32 %1, 0
%19 = icmp sgt i32 %2, 0
%20 = icmp sgt i32 %3, 0
%21 = zext i32 %2 to i64
%22 = zext i32 %3 to i64
%23 = and i64 %22, 4294967292
%24 = add nsw i64 %23, -4
%25 = lshr exact i64 %24, 2
%26 = add nuw nsw i64 %25, 1
%27 = icmp ult i32 %3, 4
%28 = and i64 %22, 4294967292
%29 = insertelement <2 x double> poison, double %15, i64 0
%30 = shufflevector <2 x double> %29, <2 x double> poison, <2 x i32> zeroinitializer
%31 = insertelement <2 x double> poison, double %15, i64 0
%32 = shufflevector <2 x double> %31, <2 x double> poison, <2 x i32> zeroinitializer
%33 = and i64 %26, 1
%34 = icmp eq i64 %24, 0
%35 = and i64 %26, 9223372036854775806
%36 = icmp eq i64 %33, 0
%37 = icmp eq i64 %28, %22
br label %38
38: ; preds = %17, %51
%39 = phi i32 [ 0, %17 ], [ %52, %51 ]
br i1 %18, label %40, label %51
40: ; preds = %38
%41 = sitofp i32 %39 to double
%42 = mul nsw i32 %39, %1
br label %44
43: ; preds = %51, %5
ret void
44: ; preds = %40, %115
%45 = phi i32 [ 0, %40 ], [ %116, %115 ]
br i1 %19, label %46, label %115
46: ; preds = %44
%47 = sitofp i32 %45 to double
%48 = tail call double @llvm.fmuladd.f64(double %47, double %13, double %41)
%49 = add i32 %45, %42
%50 = mul i32 %49, %2
br label %54
51: ; preds = %115, %38
%52 = add nuw nsw i32 %39, 1
%53 = icmp eq i32 %52, %0
br i1 %53, label %43, label %38
54: ; preds = %46, %118
%55 = phi i64 [ 0, %46 ], [ %119, %118 ]
br i1 %20, label %56, label %118
56: ; preds = %54
%57 = trunc i64 %55 to i32
%58 = sitofp i32 %57 to double
%59 = tail call double @llvm.fmuladd.f64(double %58, double %14, double %48)
%60 = add i32 %50, %57
%61 = mul i32 %60, %3
%62 = sext i32 %61 to i64
br i1 %27, label %113, label %63
63: ; preds = %56
%64 = insertelement <2 x double> poison, double %59, i64 0
%65 = shufflevector <2 x double> %64, <2 x double> poison, <2 x i32> zeroinitializer
%66 = insertelement <2 x double> poison, double %59, i64 0
%67 = shufflevector <2 x double> %66, <2 x double> poison, <2 x i32> zeroinitializer
br i1 %34, label %98, label %68
68: ; preds = %63, %68
%69 = phi i64 [ %94, %68 ], [ 0, %63 ]
%70 = phi <2 x i32> [ %95, %68 ], [ <i32 0, i32 1>, %63 ]
%71 = phi i64 [ %96, %68 ], [ 0, %63 ]
%72 = add <2 x i32> %70, <i32 2, i32 2>
%73 = sitofp <2 x i32> %70 to <2 x double>
%74 = sitofp <2 x i32> %72 to <2 x double>
%75 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %73, <2 x double> %30, <2 x double> %65)
%76 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %74, <2 x double> %32, <2 x double> %67)
%77 = add nsw i64 %69, %62
%78 = getelementptr inbounds double, double* %4, i64 %77
%79 = bitcast double* %78 to <2 x double>*
store <2 x double> %75, <2 x double>* %79
%80 = getelementptr inbounds double, double* %78, i64 2
%81 = bitcast double* %80 to <2 x double>*
store <2 x double> %76, <2 x double>* %81
%82 = or i64 %69, 4
%83 = add <2 x i32> %70, <i32 4, i32 4>
%84 = add <2 x i32> %70, <i32 6, i32 6>
%85 = sitofp <2 x i32> %83 to <2 x double>
%86 = sitofp <2 x i32> %84 to <2 x double>
%87 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %85, <2 x double> %30, <2 x double> %65)
%88 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %86, <2 x double> %32, <2 x double> %67)
%89 = add nsw i64 %82, %62
%90 = getelementptr inbounds double, double* %4, i64 %89
%91 = bitcast double* %90 to <2 x double>*
store <2 x double> %87, <2 x double>* %91
%92 = getelementptr inbounds double, double* %90, i64 2
%93 = bitcast double* %92 to <2 x double>*
store <2 x double> %88, <2 x double>* %93
%94 = add nuw i64 %69, 8
%95 = add <2 x i32> %70, <i32 8, i32 8>
%96 = add i64 %71, 2
%97 = icmp eq i64 %96, %35
br i1 %97, label %98, label %68
98: ; preds = %68, %63
%99 = phi i64 [ 0, %63 ], [ %94, %68 ]
%100 = phi <2 x i32> [ <i32 0, i32 1>, %63 ], [ %95, %68 ]
br i1 %36, label %112, label %101
101: ; preds = %98
%102 = add <2 x i32> %100, <i32 2, i32 2>
%103 = sitofp <2 x i32> %100 to <2 x double>
%104 = sitofp <2 x i32> %102 to <2 x double>
%105 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %103, <2 x double> %30, <2 x double> %65)
%106 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %104, <2 x double> %32, <2 x double> %67)
%107 = add nsw i64 %99, %62
%108 = getelementptr inbounds double, double* %4, i64 %107
%109 = bitcast double* %108 to <2 x double>*
store <2 x double> %105, <2 x double>* %109
%110 = getelementptr inbounds double, double* %108, i64 2
%111 = bitcast double* %110 to <2 x double>*
store <2 x double> %106, <2 x double>* %111
br label %112
112: ; preds = %98, %101
br i1 %37, label %118, label %113
113: ; preds = %56, %112
%114 = phi i64 [ 0, %56 ], [ %28, %112 ]
br label %121
115: ; preds = %118, %44
%116 = add nuw nsw i32 %45, 1
%117 = icmp eq i32 %116, %1
br i1 %117, label %51, label %44
118: ; preds = %121, %112, %54
%119 = add nuw nsw i64 %55, 1
%120 = icmp eq i64 %119, %21
br i1 %120, label %115, label %54
121: ; preds = %113, %121
%122 = phi i64 [ %128, %121 ], [ %114, %113 ]
%123 = trunc i64 %122 to i32
%124 = sitofp i32 %123 to double
%125 = tail call double @llvm.fmuladd.f64(double %124, double %15, double %59)
%126 = add nsw i64 %122, %62
%127 = getelementptr inbounds double, double* %4, i64 %126
store double %125, double* %127
%128 = add nuw nsw i64 %122, 1
%129 = icmp eq i64 %128, %22
br i1 %129, label %118, label %121
}