void feat_gettraj(int N, float *feats, int D, int d, float *traj) {
  for (int n = 0; n < N; n++)
    traj[n] = feats[n * D + d];
}


x86 assembly (clang, -O0):
<feat_gettraj>:
push   %rbp
mov    %rsp,%rbp
mov    %edi,-0x4(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %edx,-0x14(%rbp)
mov    %ecx,-0x18(%rbp)
mov    %r8,-0x20(%rbp)
movl   $0x0,-0x24(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    1159 <feat_gettraj+0x59>
mov    -0x10(%rbp),%rax
mov    -0x24(%rbp),%ecx
imul   -0x14(%rbp),%ecx
add    -0x18(%rbp),%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x20(%rbp),%rax
movslq -0x24(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    111c <feat_gettraj+0x1c>
pop    %rbp
ret

x86 assembly (clang, -O2):
<feat_gettraj>:
push   %r15
push   %r14
push   %rbx
test   %edi,%edi
jle    12bc <feat_gettraj+0x1bc>
movslq %edx,%r9
movslq %ecx,%r10
mov    %edi,%r14d
xor    %ecx,%ecx
cmp    $0x8,%edi
jb     11f1 <feat_gettraj+0xf1>
cmp    $0x1,%edx
jne    11f1 <feat_gettraj+0xf1>
lea    (%r10,%r14,1),%rax
lea    (%rsi,%rax,4),%rax
cmp    %r8,%rax
jbe    114b <feat_gettraj+0x4b>
lea    (%r8,%r14,4),%rax
lea    (%rsi,%r10,4),%rcx
cmp    %rax,%rcx
jae    114b <feat_gettraj+0x4b>
xor    %ecx,%ecx
jmp    11f1 <feat_gettraj+0xf1>
mov    %r14d,%ecx
and    $0xfffffff8,%ecx
lea    -0x8(%rcx),%rax
mov    %rax,%r11
shr    $0x3,%r11
add    $0x1,%r11
test   %rax,%rax
je     12c2 <feat_gettraj+0x1c2>
mov    %r11,%rdx
and    $0xfffffffffffffffe,%rdx
mov    %r9,%rbx
shl    $0x5,%rbx
lea    (%rsi,%r10,4),%rax
add    $0x10,%rax
mov    %r9,%r15
shl    $0x6,%r15
xor    %edi,%edi
nopl   0x0(%rax,%rax,1)
movups -0x10(%rax),%xmm0
movups (%rax),%xmm1
movups %xmm0,(%r8,%rdi,4)
movups %xmm1,0x10(%r8,%rdi,4)
movups -0x10(%rax,%rbx,1),%xmm0
movups (%rax,%rbx,1),%xmm1
movups %xmm0,0x20(%r8,%rdi,4)
movups %xmm1,0x30(%r8,%rdi,4)
add    $0x10,%rdi
add    %r15,%rax
add    $0xfffffffffffffffe,%rdx
jne    1190 <feat_gettraj+0x90>
test   $0x1,%r11b
je     11e8 <feat_gettraj+0xe8>
mov    %rdi,%rax
imul   %r9,%rax
add    %r10,%rax
movups (%rsi,%rax,4),%xmm0
movups 0x10(%rsi,%rax,4),%xmm1
movups %xmm0,(%r8,%rdi,4)
movups %xmm1,0x10(%r8,%rdi,4)
cmp    %r14,%rcx
je     12bc <feat_gettraj+0x1bc>
mov    %rcx,%rdx
not    %rdx
add    %r14,%rdx
mov    %r14,%rax
and    $0x3,%rax
je     1237 <feat_gettraj+0x137>
mov    %rcx,%rdi
imul   %r9,%rdi
add    %r10,%rdi
lea    (%rsi,%rdi,4),%rdi
lea    0x0(,%r9,4),%rbx
nopl   0x0(%rax)
movss  (%rdi),%xmm0
movss  %xmm0,(%r8,%rcx,4)
add    $0x1,%rcx
add    %rbx,%rdi
add    $0xffffffffffffffff,%rax
jne    1220 <feat_gettraj+0x120>
cmp    $0x3,%rdx
jb     12bc <feat_gettraj+0x1bc>
lea    0x3(%rcx),%rax
imul   %r9,%rax
lea    0x0(,%r10,4),%r15
lea    (%r15,%rax,4),%r11
mov    %r9,%rdi
shl    $0x4,%rdi
lea    0x2(%rcx),%rax
imul   %r9,%rax
lea    (%r15,%rax,4),%rbx
mov    %rcx,%rax
imul   %r9,%rax
add    %r10,%rax
lea    0x1(%rcx),%rdx
imul   %r9,%rdx
lea    (%r15,%rdx,4),%rdx
nopw   0x0(%rax,%rax,1)
movss  (%rsi,%rax,4),%xmm0
movss  %xmm0,(%r8,%rcx,4)
movss  (%rsi,%rdx,1),%xmm0
movss  %xmm0,0x4(%r8,%rcx,4)
movss  (%rsi,%rbx,1),%xmm0
movss  %xmm0,0x8(%r8,%rcx,4)
movss  (%rsi,%r11,1),%xmm0
movss  %xmm0,0xc(%r8,%rcx,4)
add    $0x4,%rcx
add    %rdi,%rsi
cmp    %rcx,%r14
jne    1280 <feat_gettraj+0x180>
pop    %rbx
pop    %r14
pop    %r15
ret
xor    %edi,%edi
test   $0x1,%r11b
jne    11ca <feat_gettraj+0xca>
jmp    11e8 <feat_gettraj+0xe8>

LLVM IR (-O0):
define void @feat_gettraj(i32 %0, float* %1, i32 %2, i32 %3, float* %4) {
%6 = alloca i32
%7 = alloca float*
%8 = alloca i32
%9 = alloca i32
%10 = alloca float*
%11 = alloca i32
store i32 %0, i32* %6
store float* %1, float** %7
store i32 %2, i32* %8
store i32 %3, i32* %9
store float* %4, float** %10
store i32 0, i32* %11
br label %12
12: ; preds = %30, %5
%13 = load i32, i32* %11
%14 = load i32, i32* %6
%15 = icmp slt i32 %13, %14
br i1 %15, label %16, label %33
16: ; preds = %12
%17 = load float*, float** %7
%18 = load i32, i32* %11
%19 = load i32, i32* %8
%20 = mul nsw i32 %18, %19
%21 = load i32, i32* %9
%22 = add nsw i32 %20, %21
%23 = sext i32 %22 to i64
%24 = getelementptr inbounds float, float* %17, i64 %23
%25 = load float, float* %24
%26 = load float*, float** %10
%27 = load i32, i32* %11
%28 = sext i32 %27 to i64
%29 = getelementptr inbounds float, float* %26, i64 %28
store float %25, float* %29
br label %30
30: ; preds = %16
%31 = load i32, i32* %11
%32 = add nsw i32 %31, 1
store i32 %32, i32* %11
br label %12
33: ; preds = %12
ret void
}

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