void weight_matrix_mean(float **weight_matrix, int motif_length,
                        float *mean_weight_matrix) {
  int i;
  int j;
  float current_mean;
  for (i = 1; i <= motif_length; i++) {
    current_mean = 0;
    for (j = 1; j <= 4; j++)
      current_mean += weight_matrix[j][i];
    mean_weight_matrix[i] = current_mean;
  }
}


x86 assembly (clang, -O0):
<weight_matrix_mean>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
mov    %rdx,-0x18(%rbp)
movl   $0x1,-0x1c(%rbp)
mov    -0x1c(%rbp),%eax
cmp    -0xc(%rbp),%eax
jg     1188 <weight_matrix_mean+0x88>
xorps  %xmm0,%xmm0
movss  %xmm0,-0x24(%rbp)
movl   $0x1,-0x20(%rbp)
cmpl   $0x4,-0x20(%rbp)
jg     1168 <weight_matrix_mean+0x68>
mov    -0x8(%rbp),%rax
movslq -0x20(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x1c(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
addss  -0x24(%rbp),%xmm0
movss  %xmm0,-0x24(%rbp)
mov    -0x20(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x20(%rbp)
jmp    1131 <weight_matrix_mean+0x31>
movss  -0x24(%rbp),%xmm0
mov    -0x18(%rbp),%rax
movslq -0x1c(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x1c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1c(%rbp)
jmp    1116 <weight_matrix_mean+0x16>
pop    %rbp
ret

x86 assembly (clang, -O2):
<weight_matrix_mean>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
test   %esi,%esi
jle    12bc <weight_matrix_mean+0x1bc>
mov    %esi,%ebp
lea    0x1(%rbp),%r8
mov    0x8(%rdi),%r9
mov    0x10(%rdi),%r10
mov    0x18(%rdi),%rcx
mov    0x20(%rdi),%rdi
mov    $0x1,%eax
cmp    $0x8,%esi
jb     1230 <weight_matrix_mean+0x130>
mov    %rbp,-0x8(%rsp)
lea    0x4(%rdx),%rbp
lea    (%rdx,%r8,4),%rsi
lea    0x4(%r9),%r14
lea    (%r9,%r8,4),%rbx
lea    0x4(%r10),%r15
lea    (%r10,%r8,4),%r12
cmp    %rbx,%rbp
setb   -0x9(%rsp)
cmp    %rsi,%r14
lea    0x4(%rcx),%r11
setb   %r14b
cmp    %r12,%rbp
setb   %r12b
cmp    %rsi,%r15
lea    (%rcx,%r8,4),%rbx
setb   %r15b
cmp    %rbx,%rbp
setb   %r13b
cmp    %rsi,%r11
lea    (%rdi,%r8,4),%rbx
setb   %r11b
cmp    %rbx,%rbp
lea    0x4(%rdi),%rbx
setb   %bpl
cmp    %rsi,%rbx
setb   %sil
test   %r14b,-0x9(%rsp)
jne    1230 <weight_matrix_mean+0x130>
and    %r15b,%r12b
jne    1230 <weight_matrix_mean+0x130>
and    %r11b,%r13b
jne    1230 <weight_matrix_mean+0x130>
and    %sil,%bpl
jne    1230 <weight_matrix_mean+0x130>
mov    -0x8(%rsp),%rbp
mov    %ebp,%esi
and    $0xfffffff8,%esi
lea    0x1(%rsi),%rax
xor    %ebx,%ebx
xorps  %xmm0,%xmm0
nopl   (%rax)
movups 0x4(%r9,%rbx,4),%xmm1
movups 0x14(%r9,%rbx,4),%xmm2
addps  %xmm0,%xmm1
addps  %xmm0,%xmm2
movups 0x4(%r10,%rbx,4),%xmm3
addps  %xmm1,%xmm3
movups 0x14(%r10,%rbx,4),%xmm1
addps  %xmm2,%xmm1
movups 0x4(%rcx,%rbx,4),%xmm2
addps  %xmm3,%xmm2
movups 0x14(%rcx,%rbx,4),%xmm3
addps  %xmm1,%xmm3
movups 0x4(%rdi,%rbx,4),%xmm1
addps  %xmm2,%xmm1
movups 0x14(%rdi,%rbx,4),%xmm2
addps  %xmm3,%xmm2
movups %xmm1,0x4(%rdx,%rbx,4)
movups %xmm2,0x14(%rdx,%rbx,4)
add    $0x8,%rbx
cmp    %rbx,%rsi
jne    11d0 <weight_matrix_mean+0xd0>
cmp    %rsi,%rbp
je     12bc <weight_matrix_mean+0x1bc>
mov    %rax,%rsi
not    %rsi
test   $0x1,%r8b
jne    125e <weight_matrix_mean+0x15e>
xorps  %xmm0,%xmm0
addss  (%r9,%rax,4),%xmm0
addss  (%r10,%rax,4),%xmm0
addss  (%rcx,%rax,4),%xmm0
addss  (%rdi,%rax,4),%xmm0
movss  %xmm0,(%rdx,%rax,4)
add    $0x1,%rax
add    %r8,%rsi
je     12bc <weight_matrix_mean+0x1bc>
xorps  %xmm0,%xmm0
cs nopw 0x0(%rax,%rax,1)
movss  (%r9,%rax,4),%xmm1
addss  %xmm0,%xmm1
addss  (%r10,%rax,4),%xmm1
addss  (%rcx,%rax,4),%xmm1
addss  (%rdi,%rax,4),%xmm1
movss  %xmm1,(%rdx,%rax,4)
movss  0x4(%r9,%rax,4),%xmm1
addss  %xmm0,%xmm1
addss  0x4(%r10,%rax,4),%xmm1
addss  0x4(%rcx,%rax,4),%xmm1
addss  0x4(%rdi,%rax,4),%xmm1
movss  %xmm1,0x4(%rdx,%rax,4)
add    $0x2,%rax
cmp    %rax,%r8
jne    1270 <weight_matrix_mean+0x170>
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

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