int solcon_extrinsic_product(int a, int b, int n1, int n2, double *s1,
                             double *s2, int pn, double *ps) {
  int i, inds1;
  ps[0] = 0.0;
  ps[1] = 0.0;
  for (i = 2; i < 2 * (n1 - a) + 2; i++)
    ps[i] = s1[i];
  inds1 = 2 * (n1 - a);
  for (i = 2; i < 2 * (n2 - b) + 2; i++)
    ps[inds1 + i] = s2[i];
  i = inds1 + i;
  while (i < 2 * pn + 5)
    ps[i++] = 0.0;
  return 0;
}


x86 assembly (clang, -O0):
<solcon_extrinsic_product>:
push   %rbp
mov    %rsp,%rbp
mov    0x18(%rbp),%rax
mov    0x10(%rbp),%eax
mov    %edi,-0x4(%rbp)
mov    %esi,-0x8(%rbp)
mov    %edx,-0xc(%rbp)
mov    %ecx,-0x10(%rbp)
mov    %r8,-0x18(%rbp)
mov    %r9,-0x20(%rbp)
mov    0x18(%rbp),%rax
xorps  %xmm0,%xmm0
movsd  %xmm0,(%rax)
mov    0x18(%rbp),%rax
xorps  %xmm0,%xmm0
movsd  %xmm0,0x8(%rax)
movl   $0x2,-0x24(%rbp)
mov    -0x24(%rbp),%eax
mov    -0xc(%rbp),%ecx
sub    -0x4(%rbp),%ecx
shl    $0x1,%ecx
add    $0x2,%ecx
cmp    %ecx,%eax
jge    117c <solcon_extrinsic_product+0x7c>
mov    -0x18(%rbp),%rax
movslq -0x24(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    0x18(%rbp),%rax
movslq -0x24(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    113d <solcon_extrinsic_product+0x3d>
mov    -0xc(%rbp),%eax
sub    -0x4(%rbp),%eax
shl    $0x1,%eax
mov    %eax,-0x28(%rbp)
movl   $0x2,-0x24(%rbp)
mov    -0x24(%rbp),%eax
mov    -0x10(%rbp),%ecx
sub    -0x8(%rbp),%ecx
shl    $0x1,%ecx
add    $0x2,%ecx
cmp    %ecx,%eax
jge    11d3 <solcon_extrinsic_product+0xd3>
mov    -0x20(%rbp),%rax
movslq -0x24(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    0x18(%rbp),%rax
mov    -0x28(%rbp),%ecx
add    -0x24(%rbp),%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    118f <solcon_extrinsic_product+0x8f>
mov    -0x28(%rbp),%eax
add    -0x24(%rbp),%eax
mov    %eax,-0x24(%rbp)
mov    -0x24(%rbp),%eax
mov    0x10(%rbp),%ecx
shl    $0x1,%ecx
add    $0x5,%ecx
cmp    %ecx,%eax
jge    120f <solcon_extrinsic_product+0x10f>
mov    0x18(%rbp),%rax
mov    -0x24(%rbp),%ecx
mov    %ecx,%edx
add    $0x1,%edx
mov    %edx,-0x24(%rbp)
movslq %ecx,%rcx
xorps  %xmm0,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
jmp    11dc <solcon_extrinsic_product+0xdc>
xor    %eax,%eax
pop    %rbp
ret

x86 assembly (clang, -O2):
<solcon_extrinsic_product>:
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
mov    0x38(%rsp),%r11
mov    %edx,%eax
sub    %edi,%eax
lea    (%rax,%rax,1),%r10d
xorps  %xmm0,%xmm0
movups %xmm0,(%r11)
test   %eax,%eax
jle    12a7 <solcon_extrinsic_product+0x197>
mov    %r10d,%eax
or     $0x1,%eax
movslq %eax,%r15
cmp    $0x3,%r15
mov    $0x2,%ebx
mov    $0x2,%r14d
cmovge %r15,%r14
add    $0xffffffffffffffff,%r14
cmp    $0x4,%r14
jb     1285 <solcon_extrinsic_product+0x175>
lea    0x10(%r11),%r12
cmp    $0x3,%r15
mov    $0x2,%ebx
mov    $0x2,%r13d
cmovge %r15,%r13
lea    (%r8,%r13,8),%rax
add    $0x8,%rax
cmp    %rax,%r12
jae    1198 <solcon_extrinsic_product+0x88>
lea    (%r11,%r13,8),%r12
add    $0x8,%r12
lea    0x10(%r8),%rax
cmp    %r12,%rax
jb     1285 <solcon_extrinsic_product+0x175>
mov    %r14,%rbx
and    $0xfffffffffffffffc,%rbx
lea    -0x4(%rbx),%rax
mov    %rax,%r13
shr    $0x2,%r13
add    $0x1,%r13
mov    %r13d,%r12d
and    $0x3,%r12d
cmp    $0xc,%rax
jae    11c2 <solcon_extrinsic_product+0xb2>
xor    %eax,%eax
jmp    1240 <solcon_extrinsic_product+0x130>
and    $0xfffffffffffffffc,%r13
xor    %eax,%eax
nopl   0x0(%rax,%rax,1)
movups 0x10(%r8,%rax,8),%xmm0
movups 0x20(%r8,%rax,8),%xmm1
movups %xmm0,0x10(%r11,%rax,8)
movups %xmm1,0x20(%r11,%rax,8)
movups 0x30(%r8,%rax,8),%xmm0
movups 0x40(%r8,%rax,8),%xmm1
movups %xmm0,0x30(%r11,%rax,8)
movups %xmm1,0x40(%r11,%rax,8)
movups 0x50(%r8,%rax,8),%xmm0
movups 0x60(%r8,%rax,8),%xmm1
movups %xmm0,0x50(%r11,%rax,8)
movups %xmm1,0x60(%r11,%rax,8)
movups 0x70(%r8,%rax,8),%xmm0
movups 0x80(%r8,%rax,8),%xmm1
movups %xmm0,0x70(%r11,%rax,8)
movups %xmm1,0x80(%r11,%rax,8)
add    $0x10,%rax
add    $0xfffffffffffffffc,%r13
jne    11d0 <solcon_extrinsic_product+0xc0>
test   %r12,%r12
je     127c <solcon_extrinsic_product+0x16c>
cs nopw 0x0(%rax,%rax,1)
nop
lea    0x0(,%rax,8),%r13
or     $0x10,%r13
movups (%r8,%r13,1),%xmm0
movups 0x10(%r8,%r13,1),%xmm1
movups %xmm0,(%r11,%r13,1)
movups %xmm1,0x10(%r11,%r13,1)
add    $0x4,%rax
add    $0xffffffffffffffff,%r12
jne    1250 <solcon_extrinsic_product+0x140>
cmp    %rbx,%r14
je     12a7 <solcon_extrinsic_product+0x197>
or     $0x2,%rbx
add    $0xffffffffffffffff,%rbx
nopl   0x0(%rax)
movsd  0x8(%r8,%rbx,8),%xmm0
movsd  %xmm0,0x8(%r11,%rbx,8)
add    $0x1,%rbx
cmp    %r15,%rbx
jl     1290 <solcon_extrinsic_product+0x180>
mov    0x30(%rsp),%r8d
sub    %esi,%ecx
jle    13e9 <solcon_extrinsic_product+0x2d9>
movslq %r10d,%r15
lea    (%rcx,%rcx,1),%eax
add    $0x1,%eax
cltq
cmp    $0x3,%rax
mov    $0x2,%ebx
mov    $0x2,%r13d
cmovge %rax,%r13
add    $0xffffffffffffffff,%r13
cmp    $0x4,%r13
jb     13bc <solcon_extrinsic_product+0x2ac>
lea    (%r11,%r15,8),%r14
add    $0x10,%r14
cmp    $0x3,%rax
mov    $0x2,%ebx
mov    $0x2,%ecx
cmovge %rax,%rcx
lea    (%r9,%rcx,8),%rsi
add    $0x8,%rsi
cmp    %rsi,%r14
jae    131f <solcon_extrinsic_product+0x20f>
add    %r15,%rcx
lea    (%r11,%rcx,8),%rcx
add    $0x8,%rcx
lea    0x10(%r9),%rsi
cmp    %rcx,%rsi
jb     13bc <solcon_extrinsic_product+0x2ac>
mov    %r13,%r14
and    $0xfffffffffffffffc,%r14
lea    -0x4(%r14),%rcx
mov    %rcx,%r12
shr    $0x2,%r12
add    $0x1,%r12
test   %rcx,%rcx
je     1434 <solcon_extrinsic_product+0x324>
movabs $0x7ffffffffffffffc,%rbx
ff ff 7f
or     $0x2,%rbx
and    %r12,%rbx
lea    (%r11,%r15,8),%rcx
mov    $0x2,%esi
nopl   0x0(%rax,%rax,1)
movups (%r9,%rsi,8),%xmm0
movups 0x10(%r9,%rsi,8),%xmm1
movups %xmm0,(%rcx,%rsi,8)
movups %xmm1,0x10(%rcx,%rsi,8)
movups 0x20(%r9,%rsi,8),%xmm0
movups 0x30(%r9,%rsi,8),%xmm1
movups %xmm0,0x20(%rcx,%rsi,8)
movups %xmm1,0x30(%rcx,%rsi,8)
add    $0x8,%rsi
add    $0xfffffffffffffffe,%rbx
jne    1360 <solcon_extrinsic_product+0x250>
test   $0x1,%r12b
je     13b3 <solcon_extrinsic_product+0x2a3>
movups (%r9,%rsi,8),%xmm0
movups 0x10(%r9,%rsi,8),%xmm1
add    %r15,%rsi
movups %xmm0,(%r11,%rsi,8)
movups %xmm1,0x10(%r11,%rsi,8)
lea    0x2(%r14),%rbx
cmp    %r14,%r13
je     13ee <solcon_extrinsic_product+0x2de>
lea    (%r11,%r15,8),%rcx
mov    %rbx,%rsi
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
movsd  (%r9,%rsi,8),%xmm0
movsd  %xmm0,(%rcx,%rsi,8)
lea    0x1(%rsi),%rbx
cmp    %rax,%rsi
mov    %rbx,%rsi
jl     13d0 <solcon_extrinsic_product+0x2c0>
jmp    13ee <solcon_extrinsic_product+0x2de>
mov    $0x2,%ebx
add    %ebx,%r10d
lea    (%r8,%r8,1),%eax
add    $0x5,%eax
cmp    %eax,%r10d
jge    1428 <solcon_extrinsic_product+0x318>
movslq %r10d,%rax
lea    (%r11,%rax,8),%rax
add    %edi,%r8d
add    %r8d,%r8d
sub    %ebx,%r8d
add    %edx,%edx
sub    %edx,%r8d
add    $0x4,%r8d
lea    0x8(,%r8,8),%rdx
mov    %rax,%rdi
xor    %esi,%esi
call   1030 <memset@plt>
xor    %eax,%eax
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
ret
mov    $0x2,%esi
test   $0x1,%r12b
jne    139a <solcon_extrinsic_product+0x28a>
jmp    13b3 <solcon_extrinsic_product+0x2a3>

LLVM IR (-O0):
define i32 @solcon_extrinsic_product(i32 %0, i32 %1, i32 %2, i32 %3, double* %4, double* %5, i32 %6, double* %7) {
%9 = alloca i32
%10 = alloca i32
%11 = alloca i32
%12 = alloca i32
%13 = alloca double*
%14 = alloca double*
%15 = alloca i32
%16 = alloca double*
%17 = alloca i32
%18 = alloca i32
store i32 %0, i32* %9
store i32 %1, i32* %10
store i32 %2, i32* %11
store i32 %3, i32* %12
store double* %4, double** %13
store double* %5, double** %14
store i32 %6, i32* %15
store double* %7, double** %16
%19 = load double*, double** %16
%20 = getelementptr inbounds double, double* %19, i64 0
store double 0.000000e+00, double* %20
%21 = load double*, double** %16
%22 = getelementptr inbounds double, double* %21, i64 1
store double 0.000000e+00, double* %22
store i32 2, i32* %17
br label %23
23: ; preds = %41, %8
%24 = load i32, i32* %17
%25 = load i32, i32* %11
%26 = load i32, i32* %9
%27 = sub nsw i32 %25, %26
%28 = mul nsw i32 2, %27
%29 = add nsw i32 %28, 2
%30 = icmp slt i32 %24, %29
br i1 %30, label %31, label %44
31: ; preds = %23
%32 = load double*, double** %13
%33 = load i32, i32* %17
%34 = sext i32 %33 to i64
%35 = getelementptr inbounds double, double* %32, i64 %34
%36 = load double, double* %35
%37 = load double*, double** %16
%38 = load i32, i32* %17
%39 = sext i32 %38 to i64
%40 = getelementptr inbounds double, double* %37, i64 %39
store double %36, double* %40
br label %41
41: ; preds = %31
%42 = load i32, i32* %17
%43 = add nsw i32 %42, 1
store i32 %43, i32* %17
br label %23
44: ; preds = %23
%45 = load i32, i32* %11
%46 = load i32, i32* %9
%47 = sub nsw i32 %45, %46
%48 = mul nsw i32 2, %47
store i32 %48, i32* %18
store i32 2, i32* %17
br label %49
49: ; preds = %69, %44
%50 = load i32, i32* %17
%51 = load i32, i32* %12
%52 = load i32, i32* %10
%53 = sub nsw i32 %51, %52
%54 = mul nsw i32 2, %53
%55 = add nsw i32 %54, 2
%56 = icmp slt i32 %50, %55
br i1 %56, label %57, label %72
57: ; preds = %49
%58 = load double*, double** %14
%59 = load i32, i32* %17
%60 = sext i32 %59 to i64
%61 = getelementptr inbounds double, double* %58, i64 %60
%62 = load double, double* %61
%63 = load double*, double** %16
%64 = load i32, i32* %18
%65 = load i32, i32* %17
%66 = add nsw i32 %64, %65
%67 = sext i32 %66 to i64
%68 = getelementptr inbounds double, double* %63, i64 %67
store double %62, double* %68
br label %69
69: ; preds = %57
%70 = load i32, i32* %17
%71 = add nsw i32 %70, 1
store i32 %71, i32* %17
br label %49
72: ; preds = %49
%73 = load i32, i32* %18
%74 = load i32, i32* %17
%75 = add nsw i32 %73, %74
store i32 %75, i32* %17
br label %76
76: ; preds = %82, %72
%77 = load i32, i32* %17
%78 = load i32, i32* %15
%79 = mul nsw i32 2, %78
%80 = add nsw i32 %79, 5
%81 = icmp slt i32 %77, %80
br i1 %81, label %82, label %88
82: ; preds = %76
%83 = load double*, double** %16
%84 = load i32, i32* %17
%85 = add nsw i32 %84, 1
store i32 %85, i32* %17
%86 = sext i32 %84 to i64
%87 = getelementptr inbounds double, double* %83, i64 %86
store double 0.000000e+00, double* %87
br label %76
88: ; preds = %76
ret i32 0
}

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