void getPM(double *PMat, double *P, int N, int lS, int tS, int t, int l,
           int a) {
  for (int i = 0; i < N; i++) {
    for (int j = 0; j < N; j++) {
      PMat[a * tS * lS * N * N + t * lS * N * N + l * N * N + i * N + j] =
          P[i * N + j];
    }
  }
}


x86 assembly (clang, -O0):
<getPM>:
push   %rbp
mov    %rsp,%rbp
mov    0x18(%rbp),%eax
mov    0x10(%rbp),%eax
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %edx,-0x14(%rbp)
mov    %ecx,-0x18(%rbp)
mov    %r8d,-0x1c(%rbp)
mov    %r9d,-0x20(%rbp)
movl   $0x0,-0x24(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    11c6 <getPM+0xc6>
movl   $0x0,-0x28(%rbp)
mov    -0x28(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    11b3 <getPM+0xb3>
mov    -0x10(%rbp),%rax
mov    -0x24(%rbp),%ecx
imul   -0x14(%rbp),%ecx
add    -0x28(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x8(%rbp),%rax
mov    0x18(%rbp),%ecx
imul   -0x1c(%rbp),%ecx
imul   -0x18(%rbp),%ecx
imul   -0x14(%rbp),%ecx
imul   -0x14(%rbp),%ecx
mov    -0x20(%rbp),%edx
imul   -0x18(%rbp),%edx
imul   -0x14(%rbp),%edx
imul   -0x14(%rbp),%edx
add    %edx,%ecx
mov    0x10(%rbp),%edx
imul   -0x14(%rbp),%edx
imul   -0x14(%rbp),%edx
add    %edx,%ecx
mov    -0x24(%rbp),%edx
imul   -0x14(%rbp),%edx
add    %edx,%ecx
add    -0x28(%rbp),%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x28(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x28(%rbp)
jmp    113a <getPM+0x3a>
jmp    11b8 <getPM+0xb8>
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    1127 <getPM+0x27>
pop    %rbp
ret

x86 assembly (clang, -O2):
<getPM>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
test   %edx,%edx
jle    1368 <getPM+0x268>
mov    %rsi,%r11
imul   0x40(%rsp),%r8d
add    %r9d,%r8d
imul   %ecx,%r8d
add    0x38(%rsp),%r8d
mov    %edx,%r10d
imul   %r10d,%r10d
imul   %r8d,%r10d
mov    %edx,%esi
lea    -0x1(%rsi),%rcx
lea    (%rdi,%rsi,8),%rax
mov    %rax,-0x18(%rsp)
mov    %esi,%eax
and    $0xfffffffc,%eax
mov    %rax,-0x28(%rsp)
add    $0xfffffffffffffffc,%rax
mov    %rax,-0x30(%rsp)
mov    %rax,%rbp
shr    $0x2,%rbp
add    $0x1,%rbp
mov    %rbp,%rax
mov    %rbp,-0x38(%rsp)
and    $0xfffffffffffffffe,%rbp
mov    %rbp,-0x40(%rsp)
mov    %rsi,%r8
neg    %r8
lea    0x4(%r10),%r15d
lea    0x30(%r11),%rax
lea    0x0(,%rsi,8),%rbp
mov    %rbp,-0x10(%rsp)
xor    %r14d,%r14d
mov    %r10,%r12
mov    %r11,%r13
mov    %r11,-0x20(%rsp)
mov    %edx,-0x44(%rsp)
mov    %rcx,-0x8(%rsp)
mov    %r10,-0x50(%rsp)
jmp    11d3 <getPM+0xd3>
nopl   0x0(%rax)
add    $0x1,%r14
add    %edx,%r15d
mov    -0x10(%rsp),%rcx
add    %rcx,%rax
add    %rcx,%r13
add    %rsi,%r12
cmp    %rsi,%r14
mov    -0x8(%rsp),%rcx
je     1368 <getPM+0x268>
mov    %r14,%r9
imul   %rsi,%r9
cmp    $0x4,%edx
jae    1260 <getPM+0x160>
xor    %ebp,%ebp
mov    %rbp,%rbx
test   $0x1,%sil
je     1210 <getPM+0x110>
mov    %r9,%rcx
add    %rbp,%rcx
movsd  (%r11,%rcx,8),%xmm0
lea    (%r9,%rbp,1),%ecx
add    %r10d,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rdi,%rcx,8)
mov    %rbp,%rbx
or     $0x1,%rbx
not    %rbp
cmp    %r8,%rbp
je     11b0 <getPM+0xb0>
lea    0x1(%rbx),%rbp
add    %r12d,%ebx
nop
movsd  -0x8(%r13,%rbp,8),%xmm0
movslq %ebx,%rbx
movsd  %xmm0,(%rdi,%rbx,8)
movsd  0x0(%r13,%rbp,8),%xmm0
lea    (%r12,%rbp,1),%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rdi,%rcx,8)
lea    (%r8,%rbp,1),%rcx
add    $0x2,%rcx
add    $0x2,%rbp
add    $0x2,%ebx
cmp    $0x1,%rcx
jne    1220 <getPM+0x120>
jmp    11b0 <getPM+0xb0>
nopl   0x0(%rax)
mov    %r14d,%ebp
imul   %edx,%ebp
add    %r10d,%ebp
lea    (%rcx,%rbp,1),%ebx
cmp    %ebp,%ebx
jl     11e3 <getPM+0xe3>
mov    %rcx,%rbp
shr    $0x20,%rbp
mov    $0x0,%ebp
jne    11e5 <getPM+0xe5>
lea    (%r10,%r9,1),%ebp
lea    (%r11,%r9,8),%r10
movslq %ebp,%rbp
mov    -0x18(%rsp),%rcx
lea    (%rcx,%rbp,8),%rbx
cmp    %rbx,%r10
jae    12bc <getPM+0x1bc>
lea    (%r9,%rsi,1),%rbx
lea    (%r11,%rbx,8),%rbx
lea    (%rdi,%rbp,8),%rbp
cmp    %rbx,%rbp
jae    12bc <getPM+0x1bc>
xor    %ebp,%ebp
mov    -0x50(%rsp),%r10
jmp    11e5 <getPM+0xe5>
cmpq   $0x0,-0x30(%rsp)
je     1312 <getPM+0x212>
mov    -0x40(%rsp),%rbx
mov    %r15d,%ebp
xor    %r10d,%r10d
nop
movups -0x30(%rax,%r10,8),%xmm0
movups -0x20(%rax,%r10,8),%xmm1
lea    -0x4(%rbp),%r11d
movslq %r11d,%rcx
movups %xmm0,(%rdi,%rcx,8)
movups %xmm1,0x10(%rdi,%rcx,8)
movups -0x10(%rax,%r10,8),%xmm0
movups (%rax,%r10,8),%xmm1
movslq %ebp,%rbp
movups %xmm0,(%rdi,%rbp,8)
movups %xmm1,0x10(%rdi,%rbp,8)
add    $0x8,%r10
add    $0x8,%ebp
add    $0xfffffffffffffffe,%rbx
jne    12d0 <getPM+0x1d0>
jmp    1315 <getPM+0x215>
xor    %r10d,%r10d
testb  $0x1,-0x38(%rsp)
mov    -0x20(%rsp),%r11
je     1349 <getPM+0x249>
mov    %r10,%rcx
add    %r9,%rcx
movups (%r11,%rcx,8),%xmm0
movups 0x10(%r11,%rcx,8),%xmm1
lea    (%r10,%r9,1),%ecx
mov    -0x50(%rsp),%rdx
add    %edx,%ecx
movslq %ecx,%rcx
movups %xmm0,(%rdi,%rcx,8)
movups %xmm1,0x10(%rdi,%rcx,8)
mov    -0x28(%rsp),%rcx
mov    %rcx,%rbp
cmp    %rsi,%rcx
mov    -0x44(%rsp),%edx
mov    -0x50(%rsp),%r10
je     11b0 <getPM+0xb0>
jmp    11e5 <getPM+0xe5>
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

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