void init_alpha(float *b_d, float *pi_d, int nstates, float *alpha_d,
                float *ones_n_d, int obs_t) {
  int i = 0;
  for (i = 0; i < nstates; ++i) {
    alpha_d[i] = pi_d[i] * b_d[(obs_t * nstates) + i];
    ones_n_d[i] = 1.0f;
  }
}


x86 assembly (clang, -O0):
<init_alpha>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %edx,-0x14(%rbp)
mov    %rcx,-0x20(%rbp)
mov    %r8,-0x28(%rbp)
mov    %r9d,-0x2c(%rbp)
movl   $0x0,-0x30(%rbp)
movl   $0x0,-0x30(%rbp)
mov    -0x30(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    1188 <init_alpha+0x88>
mov    -0x10(%rbp),%rax
movslq -0x30(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x2c(%rbp),%ecx
imul   -0x14(%rbp),%ecx
add    -0x30(%rbp),%ecx
movslq %ecx,%rcx
mulss  (%rax,%rcx,4),%xmm0
mov    -0x20(%rbp),%rax
movslq -0x30(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x28(%rbp),%rax
movslq -0x30(%rbp),%rcx
movss  0xe8b(%rip),%xmm0
movss  %xmm0,(%rax,%rcx,4)
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
jmp    1129 <init_alpha+0x29>
pop    %rbp
ret

x86 assembly (clang, -O2):
<init_alpha>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
test   %edx,%edx
jle    134b <init_alpha+0x24b>
imul   %edx,%r9d
movslq %r9d,%r9
mov    %edx,%r10d
cmp    $0x8,%edx
jae    1128 <init_alpha+0x28>
xor    %edx,%edx
jmp    12d5 <init_alpha+0x1d5>
lea    (%rcx,%r10,4),%rdx
lea    (%r8,%r10,4),%rbp
lea    (%rsi,%r10,4),%rbx
cmp    %rcx,%rbp
seta   -0x11(%rsp)
lea    (%rdi,%r9,4),%r15
cmp    %r8,%rdx
seta   %r14b
mov    %r9,-0x10(%rsp)
lea    (%r9,%r10,1),%rax
cmp    %rcx,%rbx
seta   %r12b
mov    %rdi,-0x8(%rsp)
lea    (%rdi,%rax,4),%rax
cmp    %rsi,%rdx
seta   %dil
cmp    %rcx,%rax
seta   %r13b
cmp    %rdx,%r15
setb   -0x12(%rsp)
cmp    %r8,%rbx
seta   %r11b
cmp    %rsi,%rbp
seta   %bl
cmp    %r8,%rax
seta   %r9b
cmp    %rbp,%r15
setb   %bpl
xor    %edx,%edx
test   %r14b,-0x11(%rsp)
jne    12cb <init_alpha+0x1cb>
and    %dil,%r12b
jne    12cb <init_alpha+0x1cb>
and    -0x12(%rsp),%r13b
jne    12cb <init_alpha+0x1cb>
and    %bl,%r11b
mov    -0x8(%rsp),%rdi
jne    12d0 <init_alpha+0x1d0>
and    %bpl,%r9b
mov    -0x10(%rsp),%r9
jne    12d5 <init_alpha+0x1d5>
mov    %r10d,%edx
and    $0xfffffff8,%edx
lea    -0x8(%rdx),%rax
mov    %rax,%r11
shr    $0x3,%r11
add    $0x1,%r11
test   %rax,%rax
je     12c1 <init_alpha+0x1c1>
mov    %r11,%rbx
and    $0xfffffffffffffffe,%rbx
lea    (%rdi,%r9,4),%r14
add    $0x30,%r14
xor    %eax,%eax
movaps 0xdfc(%rip),%xmm0
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
movups (%rsi,%rax,4),%xmm1
movups 0x10(%rsi,%rax,4),%xmm2
movups -0x30(%r14,%rax,4),%xmm3
mulps  %xmm1,%xmm3
movups -0x20(%r14,%rax,4),%xmm1
mulps  %xmm2,%xmm1
movups %xmm3,(%rcx,%rax,4)
movups %xmm1,0x10(%rcx,%rax,4)
movups %xmm0,(%r8,%rax,4)
movups %xmm0,0x10(%r8,%rax,4)
movups 0x20(%rsi,%rax,4),%xmm1
movups 0x30(%rsi,%rax,4),%xmm2
movups -0x10(%r14,%rax,4),%xmm3
mulps  %xmm1,%xmm3
movups (%r14,%rax,4),%xmm1
mulps  %xmm2,%xmm1
movups %xmm3,0x20(%rcx,%rax,4)
movups %xmm1,0x30(%rcx,%rax,4)
movups %xmm0,0x20(%r8,%rax,4)
movups %xmm0,0x30(%r8,%rax,4)
add    $0x10,%rax
add    $0xfffffffffffffffe,%rbx
jne    1210 <init_alpha+0x110>
test   $0x1,%r11b
je     12b7 <init_alpha+0x1b7>
movups (%rsi,%rax,4),%xmm0
movups 0x10(%rsi,%rax,4),%xmm1
lea    (%rax,%r9,1),%rbp
movups (%rdi,%rbp,4),%xmm2
mulps  %xmm0,%xmm2
movups 0x10(%rdi,%rbp,4),%xmm0
mulps  %xmm1,%xmm0
movups %xmm2,(%rcx,%rax,4)
movups %xmm0,0x10(%rcx,%rax,4)
movaps 0xd54(%rip),%xmm0
movups %xmm0,(%r8,%rax,4)
movups %xmm0,0x10(%r8,%rax,4)
cmp    %r10,%rdx
jne    12d5 <init_alpha+0x1d5>
jmp    134b <init_alpha+0x24b>
xor    %eax,%eax
test   $0x1,%r11b
jne    1280 <init_alpha+0x180>
jmp    12b7 <init_alpha+0x1b7>
mov    -0x8(%rsp),%rdi
mov    -0x10(%rsp),%r9
mov    %rdx,%rax
not    %rax
test   $0x1,%r10b
je     1300 <init_alpha+0x200>
movss  (%rsi,%rdx,4),%xmm0
lea    (%rdx,%r9,1),%rbp
mulss  (%rdi,%rbp,4),%xmm0
movss  %xmm0,(%rcx,%rdx,4)
movl   $0x3f800000,(%r8,%rdx,4)
3f
or     $0x1,%rdx
add    %r10,%rax
je     134b <init_alpha+0x24b>
lea    (%rdi,%r9,4),%rax
add    $0x4,%rax
nopl   (%rax)
movss  (%rsi,%rdx,4),%xmm0
mulss  -0x4(%rax,%rdx,4),%xmm0
movss  %xmm0,(%rcx,%rdx,4)
movl   $0x3f800000,(%r8,%rdx,4)
3f
movss  0x4(%rsi,%rdx,4),%xmm0
mulss  (%rax,%rdx,4),%xmm0
movss  %xmm0,0x4(%rcx,%rdx,4)
movl   $0x3f800000,0x4(%r8,%rdx,4)
80 3f
add    $0x2,%rdx
cmp    %rdx,%r10
jne    1310 <init_alpha+0x210>
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

LLVM IR (-O0):
define void @init_alpha(float* %0, float* %1, i32 %2, float* %3, float* %4, i32 %5) {
%7 = alloca float*
%8 = alloca float*
%9 = alloca i32
%10 = alloca float*
%11 = alloca float*
%12 = alloca i32
%13 = alloca i32
store float* %0, float** %7
store float* %1, float** %8
store i32 %2, i32* %9
store float* %3, float** %10
store float* %4, float** %11
store i32 %5, i32* %12
store i32 0, i32* %13
store i32 0, i32* %13
br label %14
14: ; preds = %42, %6
%15 = load i32, i32* %13
%16 = load i32, i32* %9
%17 = icmp slt i32 %15, %16
br i1 %17, label %18, label %45
18: ; preds = %14
%19 = load float*, float** %8
%20 = load i32, i32* %13
%21 = sext i32 %20 to i64
%22 = getelementptr inbounds float, float* %19, i64 %21
%23 = load float, float* %22
%24 = load float*, float** %7
%25 = load i32, i32* %12
%26 = load i32, i32* %9
%27 = mul nsw i32 %25, %26
%28 = load i32, i32* %13
%29 = add nsw i32 %27, %28
%30 = sext i32 %29 to i64
%31 = getelementptr inbounds float, float* %24, i64 %30
%32 = load float, float* %31
%33 = fmul float %23, %32
%34 = load float*, float** %10
%35 = load i32, i32* %13
%36 = sext i32 %35 to i64
%37 = getelementptr inbounds float, float* %34, i64 %36
store float %33, float* %37
%38 = load float*, float** %11
%39 = load i32, i32* %13
%40 = sext i32 %39 to i64
%41 = getelementptr inbounds float, float* %38, i64 %40
store float 1.000000e+00, float* %41
br label %42
42: ; preds = %18
%43 = load i32, i32* %13
%44 = add nsw i32 %43, 1
store i32 %44, i32* %13
br label %14
45: ; preds = %14
ret void
}

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