void calculateSubgridError_tauRes(int nElements_global,
                                  int nQuadraturePoints_element,
                                  int nDOF_trial_element, double *tau,
                                  double *pdeResidual, double *dpdeResidual,
                                  double *subgridError, double *dsubgridError) {
  int eN, k, j;
  for (eN = 0; eN < nElements_global; eN++)
    for (k = 0; k < nQuadraturePoints_element; k++) {
      subgridError[eN * nQuadraturePoints_element + k] =
          -tau[eN * nQuadraturePoints_element + k] *
          pdeResidual[eN * nQuadraturePoints_element + k];
      for (j = 0; j < nDOF_trial_element; j++)
        dsubgridError[eN * nQuadraturePoints_element * nDOF_trial_element +
                      k * nDOF_trial_element + j] =
            -tau[eN * nQuadraturePoints_element + k] *
            dpdeResidual[eN * nQuadraturePoints_element * nDOF_trial_element +
                         k * nDOF_trial_element + j];
    }
}


x86 assembly (clang, -O0):
<calculateSubgridError_tauRes>:
push   %rbp
mov    %rsp,%rbp
mov    0x18(%rbp),%rax
mov    0x10(%rbp),%rax
mov    %edi,-0x4(%rbp)
mov    %esi,-0x8(%rbp)
mov    %edx,-0xc(%rbp)
mov    %rcx,-0x18(%rbp)
mov    %r8,-0x20(%rbp)
mov    %r9,-0x28(%rbp)
movl   $0x0,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    125a <calculateSubgridError_tauRes+0x15a>
movl   $0x0,-0x30(%rbp)
mov    -0x30(%rbp),%eax
cmp    -0x8(%rbp),%eax
jge    1247 <calculateSubgridError_tauRes+0x147>
mov    -0x18(%rbp),%rax
mov    -0x2c(%rbp),%ecx
imul   -0x8(%rbp),%ecx
add    -0x30(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mov    -0x20(%rbp),%rax
mov    -0x2c(%rbp),%ecx
imul   -0x8(%rbp),%ecx
add    -0x30(%rbp),%ecx
movslq %ecx,%rcx
mulsd  (%rax,%rcx,8),%xmm0
mov    0x10(%rbp),%rax
mov    -0x2c(%rbp),%ecx
imul   -0x8(%rbp),%ecx
add    -0x30(%rbp),%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
movl   $0x0,-0x34(%rbp)
mov    -0x34(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    1234 <calculateSubgridError_tauRes+0x134>
mov    -0x18(%rbp),%rax
mov    -0x2c(%rbp),%ecx
imul   -0x8(%rbp),%ecx
add    -0x30(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mov    -0x28(%rbp),%rax
mov    -0x2c(%rbp),%ecx
imul   -0x8(%rbp),%ecx
imul   -0xc(%rbp),%ecx
mov    -0x30(%rbp),%edx
imul   -0xc(%rbp),%edx
add    %edx,%ecx
add    -0x34(%rbp),%ecx
movslq %ecx,%rcx
mulsd  (%rax,%rcx,8),%xmm0
mov    0x18(%rbp),%rax
mov    -0x2c(%rbp),%ecx
imul   -0x8(%rbp),%ecx
imul   -0xc(%rbp),%ecx
mov    -0x30(%rbp),%edx
imul   -0xc(%rbp),%edx
add    %edx,%ecx
add    -0x34(%rbp),%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    11a7 <calculateSubgridError_tauRes+0xa7>
jmp    1239 <calculateSubgridError_tauRes+0x139>
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
jmp    113b <calculateSubgridError_tauRes+0x3b>
jmp    124c <calculateSubgridError_tauRes+0x14c>
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    1128 <calculateSubgridError_tauRes+0x28>
pop    %rbp
ret

x86 assembly (clang, -O2):
<calculateSubgridError_tauRes>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x30,%rsp
mov    %r9,-0x60(%rsp)
mov    %esi,-0x7c(%rsp)
test   %edi,%edi
jle    14b2 <calculateSubgridError_tauRes+0x3b2>
mov    0x70(%rsp),%rbp
mov    0x68(%rsp),%r14
mov    -0x7c(%rsp),%eax
movslq %eax,%r9
mov    %edi,%esi
mov    %rsi,-0x38(%rsp)
mov    %edx,%ebx
mov    %edx,%esi
imul   %eax,%esi
mov    %rsi,-0x58(%rsp)
mov    -0x60(%rsp),%rsi
lea    (%rsi,%rbx,8),%rax
mov    %ebx,%edi
and    $0xfffffffc,%edi
mov    %rdi,-0x8(%rsp)
add    $0xfffffffffffffffc,%rdi
mov    %rdi,-0x10(%rsp)
shr    $0x2,%rdi
add    $0x1,%rdi
cmp    %rbp,%rax
seta   %r10b
lea    0x0(,%rbx,8),%rax
add    %rbp,%rax
mov    %rax,0x10(%rsp)
cmp    %rsi,%rax
seta   %al
mov    %r9,-0x40(%rsp)
mov    %r9d,%r12d
and    %r10b,%al
mov    %al,-0x7d(%rsp)
mov    %rdi,%rax
mov    %rdi,-0x18(%rsp)
and    $0xfffffffffffffffe,%rdi
mov    %rdi,-0x20(%rsp)
mov    %rbx,%rax
neg    %rax
mov    %rax,0x28(%rsp)
lea    0x30(%rbp),%rax
mov    %rax,-0x28(%rsp)
movapd 0xe40(%rip),%xmm0
movapd 0xe38(%rip),%xmm1
lea    0x30(%rsi),%rax
mov    %rax,-0x30(%rsp)
lea    0x8(%rsi),%rax
mov    %rax,0x8(%rsp)
lea    0x8(%rbp),%rax
mov    %rax,(%rsp)
xor    %eax,%eax
mov    %rax,-0x70(%rsp)
xor    %eax,%eax
mov    %rax,-0x78(%rsp)
mov    %r12,-0x68(%rsp)
jmp    122a <calculateSubgridError_tauRes+0x12a>
nopw   0x0(%rax,%rax,1)
mov    -0x78(%rsp),%rsi
add    $0x1,%rsi
mov    -0x70(%rsp),%rax
add    -0x58(%rsp),%eax
mov    %rax,-0x70(%rsp)
mov    %rsi,%rax
mov    %rsi,-0x78(%rsp)
cmp    -0x38(%rsp),%rsi
je     14b2 <calculateSubgridError_tauRes+0x3b2>
cmpl   $0x0,-0x7c(%rsp)
jle    1200 <calculateSubgridError_tauRes+0x100>
mov    -0x78(%rsp),%r15
mov    %r15,%rsi
imul   -0x58(%rsp),%rsi
mov    %rsi,0x20(%rsp)
imul   -0x40(%rsp),%r15
lea    0x1(%r15),%rax
mov    %rax,0x18(%rsp)
mov    -0x70(%rsp),%rax
mov    %eax,%r11d
xor    %r10d,%r10d
mov    %r15,-0x50(%rsp)
jmp    127c <calculateSubgridError_tauRes+0x17c>
cs nopw 0x0(%rax,%rax,1)
nop
add    $0x1,%r10
add    %edx,%r11d
cmp    %r12,%r10
je     1200 <calculateSubgridError_tauRes+0x100>
lea    (%r15,%r10,1),%r13
movsd  (%rcx,%r13,8),%xmm2
xorpd  %xmm0,%xmm2
mulsd  (%r8,%r13,8),%xmm2
movsd  %xmm2,(%r14,%r13,8)
test   %edx,%edx
jle    1270 <calculateSubgridError_tauRes+0x170>
movslq %r11d,%rax
mov    %rax,-0x48(%rsp)
mov    %r13d,%eax
imul   %edx,%eax
movslq %eax,%r9
cmp    $0x4,%edx
jae    12c0 <calculateSubgridError_tauRes+0x1c0>
xor    %eax,%eax
jmp    1400 <calculateSubgridError_tauRes+0x300>
nopw   0x0(%rax,%rax,1)
mov    %r10d,%eax
imul   %ebx,%eax
add    0x20(%rsp),%eax
mov    0x18(%rsp),%rsi
add    %r10,%rsi
cltq
mov    %rbp,%r12
lea    0x0(%rbp,%rax,8),%rdi
mov    0x10(%rsp),%rbp
lea    0x0(,%rax,8),%r15
add    %rbp,%r15
lea    (%rcx,%r13,8),%rbp
lea    (%rcx,%rsi,8),%rsi
cmp    %rsi,%rdi
setb   %sil
cmp    %r15,%rbp
setb   %al
and    %sil,%al
or     -0x7d(%rsp),%al
je     131e <calculateSubgridError_tauRes+0x21e>
xor    %eax,%eax
mov    %r12,%rbp
mov    -0x68(%rsp),%r12
mov    -0x50(%rsp),%r15
jmp    1400 <calculateSubgridError_tauRes+0x300>
cmpq   $0x0,-0x10(%rsp)
je     1498 <calculateSubgridError_tauRes+0x398>
mov    -0x28(%rsp),%rax
mov    -0x48(%rsp),%rdi
lea    (%rax,%rdi,8),%rax
mov    -0x30(%rsp),%rsi
lea    (%rsi,%rdi,8),%r15
movsd  (%rcx,%r13,8),%xmm2
unpcklpd %xmm2,%xmm2
xorpd  %xmm0,%xmm2
mov    -0x20(%rsp),%rsi
xor    %edi,%edi
mov    %r12,%rbp
mov    -0x68(%rsp),%r12
xchg   %ax,%ax
movupd -0x30(%r15,%rdi,8),%xmm3
movupd -0x20(%r15,%rdi,8),%xmm4
mulpd  %xmm2,%xmm3
mulpd  %xmm2,%xmm4
movupd %xmm3,-0x30(%rax,%rdi,8)
movupd %xmm4,-0x20(%rax,%rdi,8)
movupd -0x10(%r15,%rdi,8),%xmm3
movupd (%r15,%rdi,8),%xmm4
mulpd  %xmm2,%xmm3
mulpd  %xmm2,%xmm4
movupd %xmm3,-0x10(%rax,%rdi,8)
movupd %xmm4,(%rax,%rdi,8)
add    $0x8,%rdi
add    $0xfffffffffffffffe,%rsi
jne    1360 <calculateSubgridError_tauRes+0x260>
testb  $0x1,-0x18(%rsp)
je     13e8 <calculateSubgridError_tauRes+0x2e8>
movsd  (%rcx,%r13,8),%xmm2
unpcklpd %xmm2,%xmm2
xorpd  %xmm0,%xmm2
add    %r9,%rdi
mov    -0x60(%rsp),%rax
movupd (%rax,%rdi,8),%xmm3
movupd 0x10(%rax,%rdi,8),%xmm4
mulpd  %xmm2,%xmm3
mulpd  %xmm2,%xmm4
movupd %xmm3,0x0(%rbp,%rdi,8)
movupd %xmm4,0x10(%rbp,%rdi,8)
mov    -0x8(%rsp),%rsi
mov    %rsi,%rax
cmp    %rbx,%rsi
mov    -0x50(%rsp),%r15
je     1270 <calculateSubgridError_tauRes+0x170>
xchg   %ax,%ax
mov    %rax,%rdi
test   $0x1,%bl
je     142e <calculateSubgridError_tauRes+0x32e>
movsd  (%rcx,%r13,8),%xmm2
xorpd  %xmm0,%xmm2
add    %rax,%r9
mov    -0x60(%rsp),%rsi
mulsd  (%rsi,%r9,8),%xmm2
movsd  %xmm2,0x0(%rbp,%r9,8)
mov    %rax,%rdi
or     $0x1,%rdi
not    %rax
cmp    0x28(%rsp),%rax
je     1270 <calculateSubgridError_tauRes+0x170>
mov    0x8(%rsp),%rax
mov    %rbp,%r9
mov    -0x48(%rsp),%rbp
lea    (%rax,%rbp,8),%rax
mov    (%rsp),%rsi
lea    (%rsi,%rbp,8),%rsi
mov    %r9,%rbp
nopl   0x0(%rax,%rax,1)
movsd  (%rcx,%r13,8),%xmm2
xorpd  %xmm1,%xmm2
mulsd  -0x8(%rax,%rdi,8),%xmm2
movsd  %xmm2,-0x8(%rsi,%rdi,8)
movsd  (%rcx,%r13,8),%xmm2
xorpd  %xmm1,%xmm2
mulsd  (%rax,%rdi,8),%xmm2
movsd  %xmm2,(%rsi,%rdi,8)
add    $0x2,%rdi
cmp    %rdi,%rbx
jne    1460 <calculateSubgridError_tauRes+0x360>
jmp    1270 <calculateSubgridError_tauRes+0x170>
xor    %edi,%edi
mov    %r12,%rbp
mov    -0x68(%rsp),%r12
testb  $0x1,-0x18(%rsp)
jne    13b3 <calculateSubgridError_tauRes+0x2b3>
jmp    13e8 <calculateSubgridError_tauRes+0x2e8>
add    $0x30,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

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