void d_scale_pmat(int row, int col, double alpha, int offset, double *pA,
                  int sda) {
  const int bs = 4;
  int ii, jj;
  int na = (bs - offset % bs) % bs;
  na = row < na ? row : na;
  ii = 0;
  if (na > 0) {
    for (; ii < na; ii++) {
      for (jj = 0; jj < col; jj++) {
        pA[jj * bs] *= alpha;
      }
      pA += 1;
    }
    pA += bs * (sda - 1);
  }
  for (; ii < row - 3; ii += 4) {
    for (jj = 0; jj < col; jj++) {
      pA[0 + jj * bs] *= alpha;
      pA[1 + jj * bs] *= alpha;
      pA[2 + jj * bs] *= alpha;
      pA[3 + jj * bs] *= alpha;
    }
    pA += bs * sda;
  }
  for (; ii < row; ii++) {
    for (jj = 0; jj < col; jj++) {
      pA[jj * bs] *= alpha;
    }
    pA += 1;
  }
}


x86 assembly (clang, -O0):
<d_scale_pmat>:
push   %rbp
mov    %rsp,%rbp
mov    %edi,-0x4(%rbp)
mov    %esi,-0x8(%rbp)
movsd  %xmm0,-0x10(%rbp)
mov    %edx,-0x14(%rbp)
mov    %rcx,-0x20(%rbp)
mov    %r8d,-0x24(%rbp)
movl   $0x4,-0x28(%rbp)
mov    -0x14(%rbp),%eax
mov    $0x4,%ecx
cltd
idiv   %ecx
mov    $0x4,%eax
sub    %edx,%eax
mov    $0x4,%ecx
cltd
idiv   %ecx
mov    %edx,-0x34(%rbp)
mov    -0x4(%rbp),%eax
cmp    -0x34(%rbp),%eax
jge    1155 <d_scale_pmat+0x55>
mov    -0x4(%rbp),%eax
mov    %eax,-0x38(%rbp)
jmp    115b <d_scale_pmat+0x5b>
mov    -0x34(%rbp),%eax
mov    %eax,-0x38(%rbp)
mov    -0x38(%rbp),%eax
mov    %eax,-0x34(%rbp)
movl   $0x0,-0x2c(%rbp)
cmpl   $0x0,-0x34(%rbp)
jle    11f5 <d_scale_pmat+0xf5>
jmp    1177 <d_scale_pmat+0x77>
mov    -0x2c(%rbp),%eax
cmp    -0x34(%rbp),%eax
jge    11da <d_scale_pmat+0xda>
movl   $0x0,-0x30(%rbp)
mov    -0x30(%rbp),%eax
cmp    -0x8(%rbp),%eax
jge    11c0 <d_scale_pmat+0xc0>
movsd  -0x10(%rbp),%xmm0
mov    -0x20(%rbp),%rax
mov    -0x30(%rbp),%ecx
shl    $0x2,%ecx
movslq %ecx,%rcx
mulsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
jmp    118a <d_scale_pmat+0x8a>
mov    -0x20(%rbp),%rax
add    $0x8,%rax
mov    %rax,-0x20(%rbp)
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    1177 <d_scale_pmat+0x77>
mov    -0x24(%rbp),%ecx
sub    $0x1,%ecx
shl    $0x2,%ecx
mov    -0x20(%rbp),%rax
movslq %ecx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    %rax,-0x20(%rbp)
jmp    11fa <d_scale_pmat+0xfa>
mov    -0x2c(%rbp),%eax
mov    -0x4(%rbp),%ecx
sub    $0x3,%ecx
cmp    %ecx,%eax
jge    12ce <d_scale_pmat+0x1ce>
movl   $0x0,-0x30(%rbp)
mov    -0x30(%rbp),%eax
cmp    -0x8(%rbp),%eax
jge    12a8 <d_scale_pmat+0x1a8>
movsd  -0x10(%rbp),%xmm0
mov    -0x20(%rbp),%rax
mov    -0x30(%rbp),%ecx
shl    $0x2,%ecx
add    $0x0,%ecx
movslq %ecx,%rcx
mulsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x10(%rbp),%xmm0
mov    -0x20(%rbp),%rax
mov    -0x30(%rbp),%ecx
shl    $0x2,%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
mulsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x10(%rbp),%xmm0
mov    -0x20(%rbp),%rax
mov    -0x30(%rbp),%ecx
shl    $0x2,%ecx
add    $0x2,%ecx
movslq %ecx,%rcx
mulsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x10(%rbp),%xmm0
mov    -0x20(%rbp),%rax
mov    -0x30(%rbp),%ecx
shl    $0x2,%ecx
add    $0x3,%ecx
movslq %ecx,%rcx
mulsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
jmp    1212 <d_scale_pmat+0x112>
mov    -0x24(%rbp),%ecx
shl    $0x2,%ecx
mov    -0x20(%rbp),%rax
movslq %ecx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    %rax,-0x20(%rbp)
mov    -0x2c(%rbp),%eax
add    $0x4,%eax
mov    %eax,-0x2c(%rbp)
jmp    11fa <d_scale_pmat+0xfa>
jmp    12d3 <d_scale_pmat+0x1d3>
mov    -0x2c(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    1336 <d_scale_pmat+0x236>
movl   $0x0,-0x30(%rbp)
mov    -0x30(%rbp),%eax
cmp    -0x8(%rbp),%eax
jge    131c <d_scale_pmat+0x21c>
movsd  -0x10(%rbp),%xmm0
mov    -0x20(%rbp),%rax
mov    -0x30(%rbp),%ecx
shl    $0x2,%ecx
movslq %ecx,%rcx
mulsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
jmp    12e6 <d_scale_pmat+0x1e6>
mov    -0x20(%rbp),%rax
add    $0x8,%rax
mov    %rax,-0x20(%rbp)
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    12d3 <d_scale_pmat+0x1d3>
pop    %rbp
ret

x86 assembly (clang, -O2):
<d_scale_pmat>:
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
neg    %edx
and    $0x3,%edx
cmp    %edi,%edx
cmovg  %edi,%edx
mov    %esi,%r11d
test   %edx,%edx
jle    120c <d_scale_pmat+0x10c>
lea    -0x1(%r11),%r9
mov    %r11d,%r10d
and    $0x3,%r10d
mov    %r11d,%r12d
and    $0xfffffffc,%r12d
lea    0x60(%rcx),%r14
mov    %r10,%r13
shl    $0x5,%r13
xor    %r15d,%r15d
jmp    1155 <d_scale_pmat+0x55>
add    $0x8,%rcx
add    $0x1,%r15d
add    $0x8,%r14
cmp    %edx,%r15d
je     11fc <d_scale_pmat+0xfc>
test   %esi,%esi
jle    1140 <d_scale_pmat+0x40>
cmp    $0x3,%r9
jae    1170 <d_scale_pmat+0x70>
xor    %eax,%eax
jmp    11c3 <d_scale_pmat+0xc3>
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
mov    %r14,%rbx
xor    %eax,%eax
cs nopw 0x0(%rax,%rax,1)
nop
movsd  -0x60(%rbx),%xmm1
mulsd  %xmm0,%xmm1
movsd  %xmm1,-0x60(%rbx)
movsd  -0x40(%rbx),%xmm1
mulsd  %xmm0,%xmm1
movsd  %xmm1,-0x40(%rbx)
movsd  -0x20(%rbx),%xmm1
mulsd  %xmm0,%xmm1
movsd  %xmm1,-0x20(%rbx)
movsd  (%rbx),%xmm1
mulsd  %xmm0,%xmm1
movsd  %xmm1,(%rbx)
add    $0x4,%rax
sub    $0xffffffffffffff80,%rbx
cmp    %rax,%r12
jne    1180 <d_scale_pmat+0x80>
test   %r10,%r10
je     1140 <d_scale_pmat+0x40>
shl    $0x5,%rax
add    %rcx,%rax
xor    %ebx,%ebx
cs nopw 0x0(%rax,%rax,1)
nop
movsd  (%rax,%rbx,1),%xmm1
mulsd  %xmm0,%xmm1
movsd  %xmm1,(%rax,%rbx,1)
add    $0x20,%rbx
cmp    %rbx,%r13
jne    11e0 <d_scale_pmat+0xe0>
jmp    1140 <d_scale_pmat+0x40>
lea    -0x4(,%r8,4),%eax
ff
cltq
lea    (%rcx,%rax,8),%rcx
jmp    120e <d_scale_pmat+0x10e>
xor    %edx,%edx
lea    -0x3(%rdi),%r9d
cmp    %r9d,%edx
jge    132d <d_scale_pmat+0x22d>
shl    $0x2,%r8d
movslq %r8d,%r8
mov    %r11d,%r10d
and    $0xfffffffe,%r10d
movapd %xmm0,%xmm1
unpcklpd %xmm0,%xmm1
shl    $0x3,%r8
jmp    124f <d_scale_pmat+0x14f>
nopw   0x0(%rax,%rax,1)
add    %r8,%rcx
add    $0x4,%edx
cmp    %r9d,%edx
jge    132d <d_scale_pmat+0x22d>
test   %esi,%esi
jle    1240 <d_scale_pmat+0x140>
cmp    $0x1,%esi
jne    1260 <d_scale_pmat+0x160>
xor    %eax,%eax
jmp    12e8 <d_scale_pmat+0x1e8>
nop
mov    %r10,%rax
mov    %rcx,%rbx
cs nopw 0x0(%rax,%rax,1)
movupd (%rbx),%xmm2
movupd 0x10(%rbx),%xmm3
movupd 0x20(%rbx),%xmm4
movupd 0x30(%rbx),%xmm5
movapd %xmm2,%xmm6
unpcklpd %xmm4,%xmm6
unpckhpd %xmm4,%xmm2
movapd %xmm3,%xmm4
unpcklpd %xmm5,%xmm4
unpckhpd %xmm5,%xmm3
mulpd  %xmm1,%xmm6
mulpd  %xmm1,%xmm2
mulpd  %xmm1,%xmm4
mulpd  %xmm1,%xmm3
movapd %xmm4,%xmm5
unpcklpd %xmm3,%xmm5
movapd %xmm6,%xmm7
unpcklpd %xmm2,%xmm7
movhlps %xmm4,%xmm3
movhlps %xmm6,%xmm2
movups %xmm2,0x20(%rbx)
movups %xmm3,0x30(%rbx)
movupd %xmm7,(%rbx)
movupd %xmm5,0x10(%rbx)
add    $0x40,%rbx
add    $0xfffffffffffffffe,%rax
jne    1270 <d_scale_pmat+0x170>
mov    %r10,%rax
cmp    %r11,%r10
je     1240 <d_scale_pmat+0x140>
mov    %r11,%rbx
sub    %rax,%rbx
shl    $0x5,%rax
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
movupd (%rcx,%rax,1),%xmm2
movupd 0x10(%rcx,%rax,1),%xmm3
mulpd  %xmm1,%xmm2
movupd %xmm2,(%rcx,%rax,1)
mulpd  %xmm1,%xmm3
movupd %xmm3,0x10(%rcx,%rax,1)
add    $0x20,%rax
add    $0xffffffffffffffff,%rbx
jne    1300 <d_scale_pmat+0x200>
jmp    1240 <d_scale_pmat+0x140>
cmp    %edi,%edx
jge    141c <d_scale_pmat+0x31c>
lea    -0x1(%r11),%r8
mov    %r11d,%r9d
and    $0x3,%r9d
and    $0xfffffffc,%r11d
lea    0x60(%rcx),%r10
mov    %r9,%r14
shl    $0x5,%r14
jmp    1373 <d_scale_pmat+0x273>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
add    $0x8,%rcx
add    $0x1,%edx
add    $0x8,%r10
cmp    %edi,%edx
je     141c <d_scale_pmat+0x31c>
test   %esi,%esi
jle    1360 <d_scale_pmat+0x260>
cmp    $0x3,%r8
jae    1390 <d_scale_pmat+0x290>
xor    %ebx,%ebx
jmp    13e3 <d_scale_pmat+0x2e3>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
mov    %r10,%rax
xor    %ebx,%ebx
cs nopw 0x0(%rax,%rax,1)
nop
movsd  -0x60(%rax),%xmm1
mulsd  %xmm0,%xmm1
movsd  %xmm1,-0x60(%rax)
movsd  -0x40(%rax),%xmm1
mulsd  %xmm0,%xmm1
movsd  %xmm1,-0x40(%rax)
movsd  -0x20(%rax),%xmm1
mulsd  %xmm0,%xmm1
movsd  %xmm1,-0x20(%rax)
movsd  (%rax),%xmm1
mulsd  %xmm0,%xmm1
movsd  %xmm1,(%rax)
add    $0x4,%rbx
sub    $0xffffffffffffff80,%rax
cmp    %rbx,%r11
jne    13a0 <d_scale_pmat+0x2a0>
test   %r9,%r9
je     1360 <d_scale_pmat+0x260>
shl    $0x5,%rbx
add    %rcx,%rbx
xor    %eax,%eax
cs nopw 0x0(%rax,%rax,1)
nop
movsd  (%rbx,%rax,1),%xmm1
mulsd  %xmm0,%xmm1
movsd  %xmm1,(%rbx,%rax,1)
add    $0x20,%rax
cmp    %rax,%r14
jne    1400 <d_scale_pmat+0x300>
jmp    1360 <d_scale_pmat+0x260>
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
ret

LLVM IR (-O0):
define void @d_scale_pmat(i32 %0, i32 %1, double %2, i32 %3, double* %4, i32 %5) {
%7 = alloca i32
%8 = alloca i32
%9 = alloca double
%10 = alloca i32
%11 = alloca double*
%12 = alloca i32
%13 = alloca i32
%14 = alloca i32
%15 = alloca i32
%16 = alloca i32
store i32 %0, i32* %7
store i32 %1, i32* %8
store double %2, double* %9
store i32 %3, i32* %10
store double* %4, double** %11
store i32 %5, i32* %12
store i32 4, i32* %13
%17 = load i32, i32* %10
%18 = srem i32 %17, 4
%19 = sub nsw i32 4, %18
%20 = srem i32 %19, 4
store i32 %20, i32* %16
%21 = load i32, i32* %7
%22 = load i32, i32* %16
%23 = icmp slt i32 %21, %22
br i1 %23, label %24, label %26
24: ; preds = %6
%25 = load i32, i32* %7
br label %28
26: ; preds = %6
%27 = load i32, i32* %16
br label %28
28: ; preds = %26, %24
%29 = phi i32 [ %25, %24 ], [ %27, %26 ]
store i32 %29, i32* %16
store i32 0, i32* %14
%30 = load i32, i32* %16
%31 = icmp sgt i32 %30, 0
br i1 %31, label %32, label %67
32: ; preds = %28
br label %33
33: ; preds = %57, %32
%34 = load i32, i32* %14
%35 = load i32, i32* %16
%36 = icmp slt i32 %34, %35
br i1 %36, label %37, label %60
37: ; preds = %33
store i32 0, i32* %15
br label %38
38: ; preds = %51, %37
%39 = load i32, i32* %15
%40 = load i32, i32* %8
%41 = icmp slt i32 %39, %40
br i1 %41, label %42, label %54
42: ; preds = %38
%43 = load double, double* %9
%44 = load double*, double** %11
%45 = load i32, i32* %15
%46 = mul nsw i32 %45, 4
%47 = sext i32 %46 to i64
%48 = getelementptr inbounds double, double* %44, i64 %47
%49 = load double, double* %48
%50 = fmul double %49, %43
store double %50, double* %48
br label %51
51: ; preds = %42
%52 = load i32, i32* %15
%53 = add nsw i32 %52, 1
store i32 %53, i32* %15
br label %38
54: ; preds = %38
%55 = load double*, double** %11
%56 = getelementptr inbounds double, double* %55, i64 1
store double* %56, double** %11
br label %57
57: ; preds = %54
%58 = load i32, i32* %14
%59 = add nsw i32 %58, 1
store i32 %59, i32* %14
br label %33
60: ; preds = %33
%61 = load i32, i32* %12
%62 = sub nsw i32 %61, 1
%63 = mul nsw i32 4, %62
%64 = load double*, double** %11
%65 = sext i32 %63 to i64
%66 = getelementptr inbounds double, double* %64, i64 %65
store double* %66, double** %11
br label %67
67: ; preds = %60, %28
br label %68
68: ; preds = %124, %67
%69 = load i32, i32* %14
%70 = load i32, i32* %7
%71 = sub nsw i32 %70, 3
%72 = icmp slt i32 %69, %71
br i1 %72, label %73, label %127
73: ; preds = %68
store i32 0, i32* %15
br label %74
74: ; preds = %115, %73
%75 = load i32, i32* %15
%76 = load i32, i32* %8
%77 = icmp slt i32 %75, %76
br i1 %77, label %78, label %118
78: ; preds = %74
%79 = load double, double* %9
%80 = load double*, double** %11
%81 = load i32, i32* %15
%82 = mul nsw i32 %81, 4
%83 = add nsw i32 0, %82
%84 = sext i32 %83 to i64
%85 = getelementptr inbounds double, double* %80, i64 %84
%86 = load double, double* %85
%87 = fmul double %86, %79
store double %87, double* %85
%88 = load double, double* %9
%89 = load double*, double** %11
%90 = load i32, i32* %15
%91 = mul nsw i32 %90, 4
%92 = add nsw i32 1, %91
%93 = sext i32 %92 to i64
%94 = getelementptr inbounds double, double* %89, i64 %93
%95 = load double, double* %94
%96 = fmul double %95, %88
store double %96, double* %94
%97 = load double, double* %9
%98 = load double*, double** %11
%99 = load i32, i32* %15
%100 = mul nsw i32 %99, 4
%101 = add nsw i32 2, %100
%102 = sext i32 %101 to i64
%103 = getelementptr inbounds double, double* %98, i64 %102
%104 = load double, double* %103
%105 = fmul double %104, %97
store double %105, double* %103
%106 = load double, double* %9
%107 = load double*, double** %11
%108 = load i32, i32* %15
%109 = mul nsw i32 %108, 4
%110 = add nsw i32 3, %109
%111 = sext i32 %110 to i64
%112 = getelementptr inbounds double, double* %107, i64 %111
%113 = load double, double* %112
%114 = fmul double %113, %106
store double %114, double* %112
br label %115
115: ; preds = %78
%116 = load i32, i32* %15
%117 = add nsw i32 %116, 1
store i32 %117, i32* %15
br label %74
118: ; preds = %74
%119 = load i32, i32* %12
%120 = mul nsw i32 4, %119
%121 = load double*, double** %11
%122 = sext i32 %120 to i64
%123 = getelementptr inbounds double, double* %121, i64 %122
store double* %123, double** %11
br label %124
124: ; preds = %118
%125 = load i32, i32* %14
%126 = add nsw i32 %125, 4
store i32 %126, i32* %14
br label %68
127: ; preds = %68
br label %128
128: ; preds = %152, %127
%129 = load i32, i32* %14
%130 = load i32, i32* %7
%131 = icmp slt i32 %129, %130
br i1 %131, label %132, label %155
132: ; preds = %128
store i32 0, i32* %15
br label %133
133: ; preds = %146, %132
%134 = load i32, i32* %15
%135 = load i32, i32* %8
%136 = icmp slt i32 %134, %135
br i1 %136, label %137, label %149
137: ; preds = %133
%138 = load double, double* %9
%139 = load double*, double** %11
%140 = load i32, i32* %15
%141 = mul nsw i32 %140, 4
%142 = sext i32 %141 to i64
%143 = getelementptr inbounds double, double* %139, i64 %142
%144 = load double, double* %143
%145 = fmul double %144, %138
store double %145, double* %143
br label %146
146: ; preds = %137
%147 = load i32, i32* %15
%148 = add nsw i32 %147, 1
store i32 %148, i32* %15
br label %133
149: ; preds = %133
%150 = load double*, double** %11
%151 = getelementptr inbounds double, double* %150, i64 1
store double* %151, double** %11
br label %152
152: ; preds = %149
%153 = load i32, i32* %14
%154 = add nsw i32 %153, 1
store i32 %154, i32* %14
br label %128
155: ; preds = %128
ret void
}

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